Beam Heating Research Articles

High Peak Power Mini-array Quantum Cascade Lasers Operating in Pulsed Mode

Abstract Broad area quantum cascade lasers (BA QCLs) have significant applications in many areas, but suffer from demanding pulse operating conditions and poor beam quality due to heat accumulation and generation of high order modes. A structure of mini-array is adopted to improve the heat dissipation capacity and beam quality of BA QCLs. The active region is etched to form a multi emitter and the channels are filled with InP: Fe, which acts as a lateral heat dissipation channel to improve the lateral heat dissipation efficiency. A device with λ~4.8 μm, a peak output power of 122 W at 1.2% duty cycle with a pulse of 1.5 μs is obtained in room temperature, with far-field single-lobed distribution. This result allows BA QCLs to obtain high peak power at wider pump pulse widths and higher duty cycle conditions, promotes the application of the mid-infrared laser operating in pulsed mode in the field of standoff photoacoustic chemical detection, space optical communication and so on.

The energy synergistic mechanism of coaxial heterogeneous-wavelength hybrid laser beam and its effect on welding molten pool dynamics for aluminum alloy

The coaxial heterogeneous-wavelength hybrid laser beam (HW-HLB) heat source featuring a “Gauss + Flat top” energy distribution, which integrates the 1080 nm and 915 nm laser beam, was employed in this paper to achieve excellent stability in the welding of aluminum alloy. The interaction mechanism between aluminum alloy and laser beams with different wavelengths was taken as a point cut to explore the synergistic enhancement mechanism of coaxial HW-HLB . A thermal-fluid coupling model for HW-HLB welding, considering the synergistic effect of dual lasers, was established according to the welding experiment results. The influence of fiber laser power (PF) and diode laser power (PD) on the molten pool flow behavior during welding process was investigated. The dynamic of the molten pool and the maintenance mechanism of the keyhole were elucidated through in-situ monitoring experiments and thermal-flow simulations. The results indicated that, in contrast to the single fiber laser beam (FLB), the power density of HW-HLB presents a significant increase, which results in the “avalanche-like” augmentation of the plasma plume, as well as deeper depth of the keyhole. The energy synergistic enhancement effect of coaxial HW-HLB with PF=2800 W and PD=2500 W performs stronger and thorough changes in the morphology and flow behavior of the molten pool. The introduction of 915 nm laser beam improves the multi-reflection behavior in the inner keyhole. The findings of our research provide a theoretical framework for improving the stability of the laser weld molten pool and the quality of welds in aluminum alloys through the implementation of HW-HLB.

Studies of beam ion confinement to enhance plasma performance on EAST

Abstract A key physics issue for achieving steady-state high-performance plasmas on EAST tokamak is to decrease beam-ion losses to improve plasma confinement during neutral beam injections (NBIs). To decrease the beam losses, previous counter-Ip NBI injections are upgraded to co-Ip injections. Analysis shows that due to the reversed direction of drift across the flux surfaces caused by the pitch angle, the beam prompt loss fraction decreases from about 49% to 3% after the upgrade. Moreover, because of the change of entire beam path, beam shine-through (ST) loss fraction for counter-Ip tangential and counter-Ip perpendicular injections is reversed to co-Ip tangential and co-Ip perpendicular injections, respectively. Due to change in the initial trapped-confined beam ion fraction caused by the peaked pitch profiles, the losses induced by toroidal ripple field are also reversed after the upgrade. To further improve the beam-ion confinement under the present NBI layout, the amplitudes of toroidal field are increased from 1.75 to 2.20 T. Result shows that, due to the smaller orbit width and peaked pitch angle profile, the beam prompt loss power is lower with higher toroidal field. Due to the synergy of higher initial trapped-confined beam ion fraction and narrower Goldston-White-Boozer (GWB) boundary, the loss induced by ripple diffusion is higher with higher toroidal field. The combined effect of beam ST loss, prompt loss, and ripple loss, contributes to the increase in beam ion density. The decrease in beam loss power enhances beam heating efficiency, especially the fraction of beam heating ions. Finally, comparison between simulation and measurement by 235U fission chamber (FC) indicates that the increase in neutron rate are mainly contributed by improvement of beam-ion confinement. This study can provide potential support for beam operation and high-Ti experiment on EAST tokamak.

Synthesis of high-entropy ceramics (Y0.2Yb0.2Lu0.2Eu0.2Er0.2)3Al5O12 by electron beam heating

High-entropy ceramic materials are promising for creating functional ceramics for various purposes. The unique properties of such materials make it possible to create new thermal barrier coatings based on them. Reducing costs, time and simplifying the technology of synthesis of high-entropy ceramics is the subject of numerous studies. In this work, the synthesis of high-entropy ceramics (Y0,2Yb0,2Lu0,2Eu0,2Er0,2)3Al5O12 in a powerful beam of fast electrons was carried out for the first time. A powder mixture of initial oxides, placed in the volume of a massive copper cell, was exposed in air to a short-term exposure to a beam of powerful fast electrons with an energy of 2 MeV and a beam current of 12 mA. The speed of movement of the cuvette with the powder mixture under the beam was 1 cm/s. The total irradiation time was 10 s. During the irradiation process, at least 96 % of the mass of the powder mixture melted, resulting in the formation of highly porous ceramic formations in the form of droplets. X-ray phase analysis showed that the melt droplets are high-entropy ceramics (Y0,2Yb0,2Lu0,2Eu0,2Er0,2)3Al5O12. The powder product remaining in the cuvette and not participating in the formation of droplets contains phases of the original oxides and intermediate phases of synthesized high-entropy ceramics. Electron microscopy with EDS showed a uniform distribution of elements on the surface and in the volume of the formed ceramic droplets.

Technologies to decontaminate aflatoxins in foods: a review

SummaryAflatoxins are toxic secondary metabolites produced by Aspergillus spp., found in staple food commodities. Aflatoxins are carcinogenic, teratogenic, and mutagenic and pose a serious threat to the health of humans. The identification and quantification of aflatoxins in foods is a major challenge to guarantee food safety. Therefore, developing feasible, sensitive, and robust methods for decontamination is paramount, with short processing time and negligible impact on quality. This review evaluates recent novel technologies for aflatoxins decontamination by physical methods (microwave heating, Gama and electron beam irradiation, pulse light and ultraviolet), chemical methods (ozone, natural plant extracts, and organic acids), and biological methods (atoxigenic Aspergillus strains, Trichoderma spp., and bacteria and yeast). The study highlights on the cutting‐edge technologies of smart packaging and artificial intelligence (AI). To achieve more efficiency and adaptability to different food matrices in aflatoxins decontamination, the study suggests integrating multiple strategies. The study also recommends integrating Partnership for Delivery (P4D) to share the responsibility to increase the chance for success and control aflatoxins in foods.

Kilohertz-linewidth and low-threshold single-frequency all-fiber laser utilizing self-developed Nd3+-doped fluoro-sulfo-phosphate fiber

Compared to Nd: YAG lasers, Nd3+-doped fiber lasers offer superior beam quality, compactness, and heat dissipation, especially in generating single-frequency lasers, which holds great promise for applications in optical atomic clocks, quantum computing, and high-precision bio-photonic imaging. In this study, theoretical simulations of the local environment and experimental analyses on the luminescent characteristics of what we believe to be a novel Nd3+-doped fluoro-sulfo-phosphate (FSP) laser glass were performed to mitigate the concentration and hydroxyl quenching effects. Based on that, a highly Nd3+-doped (4 mol%) FSP fiber with a large emission cross-section (3.24 × 10−20 cm2), wide bandwidth (33.7 nm), long lifetime (354 µs), and high gain coefficient (4.24 dB/cm) was designed. Utilizing this fiber, a 1065 nm SFFL with a low pump threshold of 18 mW, a narrow linewidth of 6.5 kHz, and a 0.9 µm compact all-fiber laser were demonstrated, highlighting the potential of Nd3+-doped FSP fiber in high-performance fiber lasers.

Microstructural BIB-SEM investigation of Upper Cretaceous Jordanian carbonate-rich oil shales bearing type II-S kerogen

AbstractIn this study, we use Broad Ion Beam polishing and Scanning Electron Microscopy (BIB-SEM) to characterize the microstructure of selected core samples of immature Upper Cretaceous carbonate-rich oil shales from Jordan and to link the observations to porosity and compositional and geochemical data. The aim of this study is to understand the distribution of pore space, primary organic matter, and organic sulfur on a sub-micron scale, particularly in carbonate- and silicate-dominated layers. The thermal maturity of these marine carbonate mudstone samples of pelagic origin was found to be influenced by the elevated sulfur contents in these Type II-S kerogen source rocks. This was confirmed through both organic geochemistry and BIB-SEM observations, which revealed high sulfur content. Porosity in the carbonate mudstone exists within foraminifera, and aggregates of microfossil fragments. Initially, these voids provided significant inter- and intra-particle porosity which were later filled by organic matter during diagenesis. This ‘mobile’ organic matter is interpreted as microscopic bitumen, which exists as a solid or highly viscous fluid at surface conditions. It is likely a residue of low-temperature (“early”) bitumen generation. By examining the samples before and after dichloromethane (DCM) extraction and subsequent BIB-SEM analyses, we observed that the specimens contained a significant amount of soluble organic matter (SOM), mostly present in the micropores associated with calcite. The microscopic solid bitumen is observed to remain stable even under various conditions, such as in vacuum oven conditions of 105 °C (24 h), or exposure to ultra-high vacuum, broad ion beam (heat > 70 °C) and an electron beam of 15 keV. This suggests that the solid bitumen acts as a solid at elevated temperatures and confining pressures (85 °C and 250 MPa), and its presence can lead to the buildup of significant fluid overpressures. Our observations indicate that the pores associated with calcite provide high storage capacity in the shales during the early stages of hydrocarbon generation. In contrast, it suggests that siliciclastic-rich samples are more prone to hydrofracturing as the (early) generated hydrocarbons cannot be expelled easily. These findings highlight the complex distribution and behavior of pore space, organic matter, and sulfur in shales, shedding light on their potential for hydrocarbon generation and storage. Graphical abstract

Investigation of the Composition of Samarium Monosulfide Films Obtained by Electron Beam Heating

Investigation of the Composition of Samarium Monosulfide Films Obtained by Electron Beam Heating

Designing ITER motional Stark effect line shift (MSE-LS) spectrometers.

As a part of ITER beam aided diagnostics, the design of Motional Stark Effect (MSE) diagnostic observing the emission from the Balmer-α line is underway. The physics of Stark splitting shows that the Stark manifold is polarization dependent, and the energy splitting results in a line shift proportional to the electric field. Due to the challenges of maintaining the calibration of the plasma facing mirrors in ITER, the conventional MSE polarimetry measurement technique is replaced with a spectral approach that is deemed more favorable in the ITER environment. The MSE line shift (LS) diagnostic is designed to quantify the Lorentz electric field magnitude by measuring the Stark manifold using visible spectroscopy. In the presence of large magnetic fields and high energy heating beams of 1MeV, the expected Stark splitting is much larger than in typical devices. The MSE-LS design has unique challenges requiring careful consideration and modeling of its viewing geometry and photon budget. The MSE-LS approach on ITER is promising but has stringent demands on the allowable errors for the statistical and systematic fitting uncertainties. In this study, a full system model and numerical simulations of data for each sightline are completed. For a range of optical transmission fractions, photon noise analysis is conducted to determine the statistical uncertainties. This provides guidance on the spectrometer throughput, dispersion at the detector, optics, and other design choices. A conceptual design of a high throughput spectrometer with a volume phase transmission grating is presented.

A calorimetric evaluation method for beam targets with IR imaging: Implementation for the negative ion source BATMAN Upgrade

A comprehensive calorimetric evaluation method is presented for the beam target installed in the testbed BATMAN Upgrade (BUG). BUG hosts the prototype RF-driven negative ion source for the ITER Neutral Beam Heating System. In nominal conditions BUG features an array of negative ion beamlets (5 horizontal × 14 vertical) with 20 mm spacing among them. Each beamlet is accelerated electrostatically up to a maximum of 45 kV and transports currents in the range from 0.01 to 0.05 A. BUG has two beam targets at different distances from the beam acceleration system. The furthest (2.08 m downstream) is a water-cooled diagnostic CW-calorimeter, which can measure the heat flux distribution of the beam with a resolution of 20 × 20 mm² over a surface of 600 (H) × 800 (V) mm² and water calorimetry. The closest (0.86 m downstream) beam target with sub-millimetric resolution consists of a (376 × 142 × 20 mm³) tile made of a composite anisotropic material with carbon fibers oriented in the thickness direction (1D-CFC). It can only be exposed to beams for limited time, therefore it is retractable. Infrared (IR) imaging of beam targets yields time-resolved high-resolution spatial temperature footprints. A simple algorithm is proposed to reconstruct the heat flux at the front of a beam target, from the IR footprints at the rear. The reconstructed heat flux profiles are lower and broader due to heat conduction effects in the beam target. A correction algorithm is proposed through simulation-based inference of 2D axisymmetrical FEM simulations that reduces the power and shape fit errors over an order of magnitude, to a few percent. The calorimetric evaluation method is applied to experimental shots and two advanced heat flux distribution analysis are shown. The most relevant error sources are discussed and quantified to assess systematic errors. Finally, the method is validated in BUG comparing total beam power estimations with the two calorimetric diagnostics from the CW-calorimeter, showing good agreement.

Enhancement of ionospheric heating effect by chemical release

The ionosphere can be artificially modified by employing ground-based high-power high-frequency electromagnetic waves to irradiate the ionosphere. This modification is achieved through the nonlinear interaction between the electromagnetic waves and the ionospheric plasma, leading to changes in the physical properties and structure of the ionosphere. The degree of artificial modification of the ionosphere is closely related to the heating energy density of high-frequency pump waves. Due to the high density of neutral constituents in the lower ionosphere and the high frequency of electron-neutral collisions, the energy of heating pump waves will be absorbed and attenuated during the penetration of the low ionosphere, seriously affecting the heating effect. This paper proposes a method to reduce the absorption of ionospheric heating pump waves by releasing electron attachment chemicals into low ionosphere to form a large-scale electron density hole. A model for mitigating pump waves absorption based on SF6 release is established, and the absorption at different frequencies is quantitatively calculated. The propagation characteristics of high-frequency signals in ionospheric holes are studied using a three-dimensional ray tracing method, and the results demonstrate that the chemical release method not only reduces the absorption attenuation of heating pump waves but also forms spherical electron density holes, which exhibit a focusing effect on the heating beam and enhance the heating effect. The results are of great significance for understanding the nonlinear interaction between electromagnetic wave and ionospheric plasma and improving the ionospheric heating efficiency.

Study of impurity C transport and plasma rotation in negative triangularity on the TCV tokamak

Carbon impurity transport is studied in the TCV tokamak using a charge exchange recombination diagnostic. TCVs flexible shaping capabilities were exploited to extend previous impurity transport studies to negative triangularity (δ < 0). A practical way of studying light impurity transport (like C, TCVs main impurity species due to graphite tiled walls) is to investigate the correlations between the impurity ion gradients that, in this study, highlighted significant differences between positive (PT) and negative δ (NT) plasma configurations. δ scans ( −0.6<δ<+0.6 ) were performed in limited configurations, but displayed little correlation between C temperature, rotation and density gradients for positive δ. This stiff response for δ > 0 changes for negative δ, where the evolution of ∇vtor was accompanied by variations of ∇nC over a range of negative δ, showing that transport, in NT, is affected by velocity gradients. Similar δ scans were performed with additional NBH (Neutral Beam Heating), with power steps ranging from 0.25 MW to 1.25 MW, highlighting increased momentum confinement in negative δ. Finally, the evolution of intrinsic plasma toroidal rotation across linear to saturated ohmic confinement regime (LOC/SOC) transitions was explored at δ < 0, expanding previous studies performed in TCV for δ> 0 (Bagnato et al 2023 Nucl. Fusion 63 056006). Toroidal rotation reversal was not observed for δ < 0, despite clear LOC/SOC transitions, confirming that these two phenomena occur concomitantly only in a restricted number of cases and under specific conditions.

Core density profile control by energetic ion anisotropy in LHD

Electron and impurity ion density profiles have been controlled by using tangential and perpendicular neutral beams for plasma heating in a stellarator/heliotron for the first time. Reduced anisotropy of stored energy for energetic ion En⊥/Enǁ has resulted in an inward electron and impurity transport, forming a core electron density peaking. Increased anisotropy leads to a flat or hollowed electron density profile with an impurity exhaust in a core region [Yoshinuma et al., Nucl. Fusion 49, 062002 (2009)]. A high confinement state of particles in LHD has yet to be achieved, except for a temporal state of an electron density peaking created by a pellet injection. As a pioneering and crucial research result, the operation of energetic ion anisotropy by neutral beams has newly demonstrated that the direction of the radial transport of bulk and impurity ions can be controlled. At the same time, the overall plasma performance rises in neutron flux and stored energy. On the other hand, the increase in the anisotropy flattens the density profile. This new finding holds promise for a control knob of nuclear fusion reactors to enhance fusion power output.

Development of pulsed plasma operation scenario and required conditions in JA DEMO

We have developed the pulsed plasma operation scenarios for JA DEMO, a design concept of the steady-state tokamak demonstration reactor, to clarify controls of the current profile and power required for the operation. We compare the scenarios when injecting electron cyclotron waves only and both neutral beam and electron cyclotron waves for external heating and current drive. We demonstrate current profile control that maintains the minimum value of the safety factor above one and avoids creating the local minima in the safety factor profile and power control by argon seeding that maintains the fusion power constant at the desired value and reduces the heat load on the divertor, performing long-time integrated modeling simulations. We clarify the conditions of the heating and current drive system and impurity injection system required for such control. The dependence of power control on argon anomalous transport coefficients is investigated. We have the prospect of maintaining the fusion power of 1 GW for more than two hours, i.e. obtaining the required plasma performance determined using a systems code.

Lattice destabilization in electron-beam-irradiated antimony in TEM

This study investigates the effect of electron-beam heating in transmission electron microscopy (TEM) on the structure of antimony (Sb) thin films. TEM of Sb irradiated with 300 keV electrons reveals a significant deviation (2–3°) from the theoretical interplanar angle between certain planes within the irradiated region. We relate this anomaly to electron-beam heating under electron-beam irradiation. Our FEA simulations reveal localized heating within the irradiated region caused by Auger electron excitation. These simulations predict temperatures as high as 949 K, accompanied by ultrafast heating and quenching rates. This extreme thermal cycling likely induces significant thermal stress in localized regions, possibly exceeding their elastic limits. This stress is likely related to the observed structural anomalies. This work highlights the importance of considering beam heating, particularly when interpreting TEM data from nanoscale materials with reduced thermal conductivity due to surface scattering effects.

Mass Flows in Expanding Coronal Loops

An expansion of the cross-sectional area directly impacts the mass flow along a coronal loop and significantly alters the radiative and hydrodynamic evolution of that loop as a result. Previous studies have found that an area expansion from the chromosphere to the corona significantly lengthens the cooling time of the corona and appears to suppress draining from the corona. In this work, we examine the fluid dynamics to understand how the mass flow rate, the energy balance, and the cooling and draining timescales are affected by a nonuniform area. We find that in loops with moderate or large expansion (cross-sectional area expansion factors of 2, 3, 10, 30, 100 from the photosphere to the apex), impulsive heating, for either direct thermal heating or electron beam heating, induces a steady flow into the corona, so that the coronal density continues to rise during the cooling phase, whereas a uniform loop drains during the cooling phase. The induced upflow carries energy into the corona, balancing the losses from thermal conduction, and continues until thermal conduction weakens enough so that it can no longer support the radiative losses of the transition region. As a result, the plasma cools primarily radiatively until the onset of catastrophic collapse. The speed and duration of the induced upflow both increase in proportion to the rate of area expansion. We argue that observations of blueshifted spectral lines, therefore, could place a constraint on a loop’s area expansion.

Multi-pellet injection into the NBI-heated phase of TJ-II plasmas

A pellet-induced enhanced confinement (PiEC) phase, with general characteristics similar to those reported for the stellarator W7-X, is observed after single pellet injection (>1019 H atoms) into the neutral beam injection heated phase of plasmas in the mid-sized heliac-type stellarator TJ-II. In addition to a step-like increase in density, plasma diamagnetic energy content rises significantly with respect to that of reference discharges, energy confinement time is similarly enhanced when compared to International Stellarator Scaling law predictions (Yamada et al 2005 Nucl. Fusion 45 1684) renormalized for TJ-II, and the triple product, n e · T i · τ E, exhibits a clear bifurcation towards an improved confinement branch when compared to the branch product predicted by the same law. In this work, multiple pellets are injected in series into NBI-heated plasmas in the TJ-II and post-injection plasma performance is reported and discussed. For instance, a charge-exchange recombination spectroscopy diagnostic reveals significantly increased core ion temperatures after pellet injection compared to temperatures achieved in comparable reference plasmas, this pointing to increased ion energy content and improved ion energy confinement during a PiEC phase. It is also found that enhanced performance is independent of whether co- or counter-NBI heating beam is employed. Finally, record stored diamagnetic energy content and plasma beta values are achieved when the largest available pellets are employed. The results indicate that pellet injections extend the operational regime well beyond limits previously achieved in TJ-II without pellets.

A three-probe method for accurate nanoscale thermal transport measurements

Measurements of transport phenomena are constantly plagued by contact resistance, prohibiting the sample's intrinsic electrical or thermal conductivity from being accurately determined. This predicament is particularly severe in thermal transport measurements due to the inability to meet similar impedance requirements for a four-probe method used in electrical resistance measurements. Here, we invent a three-probe measurement method that makes an accurate determination of thermal conductivity possible for nanomaterials. Incorporating electron beam heating provided by a scanning electron microscope (SEM) on a diffusive thermal conductor not only quantifies the thermal contact resistance, which may introduce an error of more than 270% to a sample's thermal conductivity, but also eliminates several device uncertainties that may contribute an additional 17% error in a measurement. The method also enables local temperature measurements, revealing nanoscale structural variations unfound by SEM. The high accuracy of the technique would make standardization of nanoscale thermal transport measurement possible.

Bayesian inference of electron density and ion temperature profiles from neutral beam and halo Balmer-α emission at Wendelstein 7-X

By employing Bayesian inference techniques, the full electron density profile from the plasma core to the edge of Wendelstein 7-X (W7-X) is inferred solely from neutral hydrogen beam and halo Balmer-α (Hα) emission data. The halo is a cloud of neutrals forming in the vicinity of the injected neutral beam due to multiple charge exchange reactions. W7-X is equipped with several neutral hydrogen beam heating sources and an Hα spectroscopy system that views these sources from different angles and penetration depths in the plasma. As the beam and halo emission form complex spectra for each spatial point that are non-linearly dependent on the plasma density profile and other parameters, a complete model from the neutral beam injection and halo formation through to the spectroscopic measurements is required. The model is used here to infer electron density profiles for a range of common W7-X plasma scenarios. The inferred profiles show good agreement with profiles determined by the Thomson scattering and interferometry diagnostics across a broad range of absolute densities without any changes to the input or fitting parameters. The time evolution of the density profile in a discharge with continuous core density peaking is successfully reconstructed, demonstrating sufficient spatial resolution to infer strongly shaped profiles. Furthermore, it is shown as a proof of concept that the model is also able to infer the main ion temperature profile using the same data set.

Exploring self-consistent 2.5D flare simulations with MPI-AMRVAC

Context. Multidimensional solar flare simulations have not yet included a detailed analysis of the lower atmospheric responses, such as downflowing chromospheric compressions and chromospheric evaporation processes. Aims. We present an analysis of multidimensional flare simulations, including an analysis of chromospheric upflows and downflows that provides important groundwork for comparing 1D and multidimensional models. Methods. We followed the evolution of a magnetohydrodynamic standard solar flare model that includes electron beams and in which localized anomalous resistivity initiates magnetic reconnection. We varied the background magnetic field strength to produce simulations that cover a large span of observationally reported solar flare strengths. Chromospheric energy fluxes and energy density maps were used to analyze the transport of energy from the corona to the lower atmosphere, and the resultant evolution of the flare. Quantities traced along 1D field lines allowed for detailed comparisons with 1D evaporation models. Results. The flares produced by varying the background coronal field strength between 20 G and 65 G have GOES classifications between B1.5 and M2.3. All produce a lobster claw reconnection outflow and a fast shock in the tail of this flow with a similar maximum Alfvén Mach number of ∼10. The impact of the reconnection outflow on the lower atmosphere and the heat conduction are the key agents driving the chromospheric evaporation and “downflowing chromospheric compressions”. The peak electron beam heating flux in the lower atmospheres varies between 1.4 × 109 and 4.7 × 1010 erg cm−2 s−1 across the simulations. The downflowing chromospheric compressions have kinetic energy signatures that reach the photosphere, but at subsonic speeds they would not generate sunquakes. The weakest flare generates a relatively dense flare loop system, despite having a negative net mass flux, through the top of the chromosphere, that is to say, more mass is supplied downward than is evaporated upward. The stronger flares all produce positive mass fluxes. Plasmoids form in the current sheets of the stronger flares due to tearing, and in all experiments the loop tops contain turbulent eddies that ring via a magnetic tuning fork process. Conclusions. The presented flares have chromospheric evaporation driven by thermal conduction and the impact and rebound of the reconnection outflow, in contrast to most 1D models where this process is driven by the beam electrons. Several multidimensional phenomena are critical in determining plasma behavior but are not generally considered in 1D flare simulations. They include loop-top turbulence, reconnection outflow jets, heat diffusion, compressive heating from the multidimensional expansion of the flux tubes due to changing pressures, and the interactions of upward and downward flows from the evaporation meeting the material squeezed downward from the loop tops.