SEARCH · Search NASA

Results for “3D Simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 739 records · Page 41

All sources NASA NTRS DOE OSTI

NASA NTRS1994

Core magnetic field enhancement in single X line, multiple X line and patchy reconnection

Magnetic flux transfer events often show a significant increase of the magnetic field strength at the center of the events. Similar magnetic field observations have been reported for structures in or near the plasma sheet of the magnetotail at about 20 R(sub E). We have carried out two-dimensional (2D) and three-dimensional (3D) and simulations of single X line reconnection (SXR), multiple X line reconnection (MXR), and patchy reconnection to determine and compare the amplification of the magnetic field in the center of the developing flux tubes. The various processes are achieved by appropriate choices of 2D or 3D resistivity models. The simulations show that the increase in magnetic field strength depends on both the property of the initial configuration and the particular reconnection geometry. For the chosen initial conditions the MXR process leads to a larger increase in the core magnetic field than the patchy reconnection and SXR caused by larger magnetic tensions in the MXR process. The 3D processes always lead to a larger amplification than the corresponding 2D processes. In the 3D cases, force imbalance in the y direction will accelerate plasma out of the flux tube. This process reduces the thermal pressure and leads to a further compression of the flux rope, which yields an additional increase in the interior magnetic field strength.

Ma, Z. W.↗

DOE OSTI2017

Nonlinear simulations of beam-driven Compressional Alfv´en Eigenmodes in NSTX

We present results for the 3D nonlinear simulations of neutral-beam-driven compressional Alfv´en eigenmodes (CAEs) in the National Spherical Torus Experiment (NSTX). Hybrid MHD-particle simulations for the H-mode NSTX discharge (shot 141398) using the HYM code show unstable CAE modes for a range of toroidal mode numbers, n = 4 - 9, and frequencies below the ion cyclotron frequency. It is found that the essential feature of CAEs is their coupling to kinetic Alfv´en wave (KAW) that occurs on the high-field side at the Alfv´en resonance location. We frequently observe high-frequency Alfv´en eigenmodes in beam-heated NSTX plasmas, and have been linked to flattening of the electron temperature profiles at high beam power. Coupling between CAE and KAW suggests an energy channeling mechanism to explain these observations, in which beam driven CAEs dissipate their energy at the resonance location, therefore significantly modifying the energy deposition profile. Nonlinear simulations demonstrate that CAEs can channel the energy of the beam ions from the injection region near the magnetic axis to the location of the resonant mode conversion at the edge of the beam density profile. Furthermore, a set of nonlinear simulations show that the CAE instability saturates due to nonlinear particle trapping, and a large fraction of beam energy can be transferred to several unstable CAEs of relatively large amplitudes and absorbed at the resonant location. Absorption rate shows a strong scaling with the beam power.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2017

On the Preconditioning of a Newton-Krylov Solver for a High-Order reconstructed Discontinuous Galerkin Discretization of All-Speed Compressible Flow with Phase Change for Application in Laser-Based Additive Manufacturing

This dissertation focuses on the development of a fully-implicit, high-order compressible ow solver with phase change. The work is motivated by laser-induced phase change applications, particularly by the need to develop large-scale multi-physics simulations of the selective laser melting (SLM) process in metal additive manufacturing (3D printing). Simulations of the SLM process require precise tracking of multi-material solid-liquid-gas interfaces, due to laser-induced melting/ solidi cation and evaporation/condensation of metal powder in an ambient gas. These rapid density variations and phase change processes tightly couple the governing equations, requiring a fully compressible framework to robustly capture the rapid density variations of the ambient gas and the melting/evaporation of the metal powder. For non-isothermal phase change, the velocity is gradually suppressed through the mushy region by a variable viscosity and Darcy source term model. The governing equations are discretized up to 4th-order accuracy with our reconstructed Discontinuous Galerkin spatial discretization scheme and up to 5th-order accuracy with L-stable fully implicit time discretization schemes (BDF2 and ESDIRK3-5). The resulting set of non-linear equations is solved using a robust Newton-Krylov method, with the Jacobian-free version of the GMRES solver for linear iterations. Due to the sti nes associated with the acoustic waves and thermal and viscous/material strength e ects, preconditioning the GMRES solver is essential. A robust and scalable approximate block factorization preconditioner was developed, which utilizes the velocity-pressure (vP) and velocity-temperature (vT) Schur complement systems. This multigrid block reduction preconditioning technique converges for high CFL/Fourier numbers and exhibits excellent parallel and algorithmic scalability on classic benchmark problems in uid dynamics (lid-driven cavity ow and natural convection heat transfer) as well as for laser-induced phase change problems in 2D and 3D.

42 ENGINEERING↗

DOE OSTI2014

Measurement of the internal stress and electric field in a resonating piezoelectric transformer for high-voltage applications using the electro-optic and photoelastic effects

The high output voltages from piezoelectric transformers are currently being used to accelerate charged particle beams for x-ray and neutron production. Traditional methods of characterizing piezoelectric transformers (PTs) using electrical probes can decrease the voltage transformation ratio of the device due to the introduction of load impedances on the order of hundreds of kiloohms to hundreds of megaohms. Consequently, an optical diagnostic was developed that used the photoelastic and electro-optic effects present in piezoelectric materials that are transparent to a given optical wavelength to determine the internal stress and electric field. The combined effects of the piezoelectric, photoelastic, and electro-optic effects result in a time-dependent change the refractive indices of the material and produce an artificially induced, time-dependent birefringence in the piezoelectric material. This induced time-dependent birefringence results in a change in the relative phase difference between the ordinary and extraordinary wave components of a helium-neon laser beam. The change in phase difference between the wave components was measured using a set of linear polarizers. The measured change in phase difference was used to calculate the stress and electric field based on the nonlinear optical properties, the piezoelectric constitutive equations, and the boundary conditions of the PT. Maximum stresses of approximately 10 MPa and electric fields of as high as 6 kV/cm were measured with the optical diagnostic. Measured results were compared to results from both a simple one-dimensional (1D) model of the piezoelectric transformer and a three-dimensional (3D) finite element model. Measured stresses and electric fields along the length of an operating length-extensional PT for two different electrical loads were within at least 50 % of 3D finite element simulated results. Additionally, the 3D finite element results were more accurate than the results from the 1D model for a wider range of electrical load impedances under test.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2025

PyCMG-based Simulation of Volumetric Concrete Microstructure

Concrete is a complex, heterogeneous material with a microstructure composed of aggregates, cement paste, and pores spanning multiple length scales. Understanding this microstructure is critical for advancing the performance, durability, and modeling of concrete-based systems. While experimental imaging such as X-ray computed tomography (XCT) provides valuable insights, generating large datasets with detailed ground truth annotations is both costly and labor-intensive due to challenges in segmenting similar phases, such as aggregates and cement paste, that often share similar attenuation properties. To address this, we developed a pipeline to simulate realistic 3D concrete microstructures using the open-source Python package PyCMG. This simulation effort focuses on generating high-fidelity, annotated microstructures that can serve as training or benchmarking datasets for image analysis, segmentation algorithms, and machine learning models, particularly in scenarios where experimental data is scarce.

Ziabari, Amir [Oak Ridge National Laboratory; ORNL↗

NASA NTRS1993

Direct simulations of premixed turbulent flames with nonunity Lewis numbers

A principal effect of turbulence on premixed flames in the flamelet regime is to wrinkle the flame fronts. For nonunity Lewis numbers, Le is not equal to 1, the local flame structure is altered in curved regions. This effect is examined using direct numerical simulations of 3D isotropic turbulence with constant density, single-step Arrhenius kinetics chemistry. Simulations of Lewis numbers 0.8, 1.0, and 1.2 are compared. At the local level, curvature effects dominated changes to the flame structure while strain effects were insignificant. A strong Lewis-number-dependent correlation was found between surface curvature and the local flame speed. The correlation was positive for Le less than 1 and negative for Le greater than 1. At the global level, strain-related effects were more significant than curvature effects. The turbulent flame speed changed significantly with Lewis number, increasing as Le decreased. This was found to be due to strain effects that have a nonzero mean over the flame surface, rather than to curvature effects that have a nearly zero mean. The mean product temperature was also found to vary with Lewis number, being higher for Le greater than 1 and lower for Le less than 1.

Rutland, C. J.↗

NASA NTRS2007

Simulation Studies of Forest Structure using 3D Lidar and Radar Models

The use of lidars and radars to measure forest structure attributes such as height and biomass are being considered for future Earth Observation missions. Large footprint lidar makes a direct measurement of the heights of scatterers in the illuminated footprint and can yield information about the vertical profile of the canopy. Synthetic Aperture Radar (SAR) is known to sense the canopy volume, especially at longer wavelengths and is useful for estimating biomass. Interferometric SAR (InSAR) has been shown to yield forest canopy height information. For example, the height of scattering phase retrieved from InSAR data is considered to be correlated with the three height and the spatial structure of the forest stand. There is much interest in exploiting these technologies separately and together to get important information for carbon cycle and ecosystem science. More detailed information of the electromagnetic radiation interactions within forest canopies is needed. And backscattering models can be of much utility here. As part of a NASA funded project to explore data fusion, a three-dimensional (3D) coherent radar backscattering model and a 3D lidar backscatter models were used to investigate the use of large footprint lidar, SAR and InSAR for characterizing realistic forest scenes. For this paper, we use stem maps and other forest measurements to develop a realistic spatial structure of a spruce-hemlock forest canopy found in Maine, USA. The radar and lidar models used measurements of the 3D forest scene as input and simulated the coherent radar backscattering signature and 1064nm energy backscatter, respectively. The relationships of backscatter derived forest structure were compared with field measurements. In addition, we also had detailed airborne lidar (Laser Imaging Vegetation Sensor, LVIS) data available over the stem map sites that was used to study the accuracies of tree height derived from modeled SAR backscatter and the scattering phase center retrieved from the simulated InSAR data will be compared with the height indices, or other structure parameters derived from the lidar data. These results will address the possible synergies between lidar and radar in data in terms of forest structural information.

Sun, Guoqing↗

DOE OSTI2016

Determination of Experimental Fuel Rod Parameters using 3D Modelling of PCMI with MPS Defect

An in-reactor experiment is being designed in order to validate the pellet-cladding mechanical interaction (PCMI) behavior of the BISON fuel performance code. The experimental parameters for the test rod being placed in the Halden Research Reactor are being determined using BISON simulations. The 3D model includes a missing pellet surface (MPS) defect to generate large local cladding deformations, which should be measureable after typical burnup times. The BISON fuel performance code is being developed at Idaho National Laboratory (INL) and is built on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. BISON supports both 2D and 3D finite elements and solves the fully coupled equations for solid mechanics, heat conduction and species diffusion. A number of fuel performance effects are included using models for swelling, densification, creep, relocation and fission gas production & release. In addition, the mechanical and thermal contact between the fuel and cladding is explicitly modelled using a master-slave based contact algorithm. In order to accurately predict PCMI effects, the BISON code includes the relevant physics involved and provides a scalable and robust solution procedure. The depth of the proposed MPS defect is being varied in the BISON model to establish an optimum value for the experiment. The experiment will be interrupted approximately every 6 months to measure cladding radial deformation and provide data to validate BISON. The complete rodlet (~20 discrete pellets) is being simulated using a 180° half symmetry 3D model with MPS defects at two axial locations. In addition, annular pellets will be used at the top and bottom of the pellet stack to allow thermocouples within the rod to measure the fuel centerline temperature. Simulation results will be presented to illustrate the expected PCMI behavior and support the chosen experimental design parameters.

BISON↗

DOE OSTI2016

Analysis of Experimental Fuel Rod Parameters using 3D Modelling of PCMI with MPS Defect

BISON↗

DOE OSTI2025

Capability in Theory, Modeling, and Validation for a Range of Innovative Fusion Concepts using High-Fidelity Moment-Kinetic Models

A computational modeling capability is created and available to the fusion community to understand and design lower-cost and innovative fusion concepts. The approach uses high- fidelity kinetic, moment-kinetic, and moment models and includes sophisticated plasma- boundary interactions. A majority of fusion-relevant simulations are performed with magnetohydrodynamic models and hybrid particle-in-cell codes, with limited-fidelity electron and kinetic physics. However, in fusion configurations like Z-pinches, field-reversed- configurations, plasma jet magneto-inertial fusion, spinning mirrors, and others, kinetic effects (both electron and ions) are critical to understand the physics and design scaling into the highly kinetic regime of a burning fusion plasma. Furthermore, as present fusion machines move towards a burning plasma regime, liquid-metal blankets are needed to handle first-wall heat- flux, reduce erosion, and eventually for energy conversion and fuel breeding. The work performed under this ARPA-E BETHE Capability Team advances the state-of-the-art in modeling and understanding plasma dynamics in fusion devices and its coupling with liquid-metal dynamics. These are critical areas of research for fusion energy to become realizable. To address these complex problems, we have leveraged and extended computational capabilities through the code, Gkeyll (developed jointly with Princeton Plasma Physics Laboratory and academic partners), for kinetic and moment modeling of fusion plasmas. The Concept Teams supported by this Capability Team include the Wisconsin High-field Axisymmetric Mirror (WHAM), Centrifugal Mirror Experiment (CFME), Plasma-Jet Magneto- Inertial Fusion (PJMIF), and solid and liquid wall plasma-material interaction studies relevant to a number of fusion concepts including Zap Energy’s Z-pinch. This software is open-source and available to the fusion community as a high-fidelity tool for the design of lower-cost fusion experiments. 3D gyrokinetic simulations of WHAM are now possible for long enough time scales to understand the evolution of interchange instabilities. 3D multi-fluid simulations of CMFE at higher Mach numbers are now possible for detailed design iterations with the goal of stability. The state-of-the-art in understanding shock formation and shock mitigation regimes in merging liners for PJMIF have been furthered by our kinetic simulations. Our novel models and frameworks studying plasma-material interaction by incorporating wall emission for various solid wall materials of relevance to pulsed and steady fusion concepts have advanced the state-of-the-art in our understanding of particle fluxes, heat fluxes, and other quantities at cathodes and anodes. The results from this work may explain discrepancies between experimental and theoretical predictions of achieved current densities in pulsed concepts such as Z-pinches. Another significant contribution of this Capability Team is the development and deployment of a novel experimental platform, LEX (Liquid Electrode eXperiment), at Virginia Tech to understand liquid metal free-surface response to electromagnetic pulses. The novel experiments along with model validation quantified the effect of different materials and sizes of liquid metal droplets on the radiative power balance of fusion plasmas for pulsed concepts. Furthermore, these experiments provided mitigation strategies for violent liquid metal response for high current pulses as would be expected in fusion regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2016

NON-LTE INVERSIONS OF THE Mg ii h and k AND UV TRIPLET LINES

The Mg ii h and k lines are powerful diagnostics for studying the solar chromosphere. They have become particularly popular with the launch of the Interface Region Imaging Spectrograph ( IRIS ) satellite, and a number of studies that include these lines have lead to great progress in understanding chromospheric heating, in many cases thanks to the support from 3D MHD simulations. In this study, we utilize another approach to analyze observations: non-LTE inversions of the Mg ii h and k and UV triplet lines including the effects of partial redistribution. Our inversion code attempts to construct a model atmosphere that is compatible with the observed spectra. We have assessed the capabilities and limitations of the inversions using the FALC atmosphere and a snapshot from a 3D radiation-MHD simulation. We find that Mg ii h and k allow reconstructing a model atmosphere from the middle photosphere to the transition region. We have also explored the capabilities of a multi-line/multi-atom setup, including the Mg ii h and k, the Ca ii 854.2 nm, and the Fe i 630.25 lines to recover the full stratification of physical parameters, including the magnetic field vector, from the photosphere to the chromosphere. Finally, we present the first inversions of observed IRIS spectra from quiet-Sun, plage, and sunspot, with very promising results.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1991

Three-dimensional direct particle simulation on the Connection Machine

This paper presents the algorithms necessary for an efficient data parallel implementation of a 3D particle simulation. In particular, a general master/slave algorithm and a fast sorting algorithm are described and the use of these algorithms in a particle simulation is outlined. A particle simulation using these algorithms has been implemented on a 32768 processor Connection Machine that is capable of simulating over 30 million particles at an average rate of 2.4-microsec/particle/step. Results are presented from the simulation of flow over an Aeroassisted Flight Experiment geometry at 100 km altitude.

Dagum, Leonardo↗

DOE OSTI2013

Geometrical investigation of the kinetic evolution of the magnetic field in a periodic flux rope

Flux ropes are bundles of magnetic field wrapped around an axis. Many laboratory, space, and astrophysics processes can be represented using this idealized concept. Here, a massively parallel 3D kinetic simulation of a periodic flux rope undergoing the kink instability is studied. The focus is on the topology of the magnetic field and its geometric structures. The analysis considers various techniques such as Poincaré maps and the quasi-separatrix layer (QSL). These are used to highlight regions with expansion or compression and changes in the connectivity of magnetic field lines and consequently to outline regions where heating and current may be generated due to magnetic reconnection. The present study is, to our knowledge, the first QSL analysis of a fully kinetic 3D particle in cell simulation and focuses the existing QSL method of analysis to periodic systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

An embedded 3D fracture modeling approach for simulating fracture-dominated fluid flow and heat transfer in geothermal reservoirs

In this paper, we describe an efficient modeling approach, named embedded discrete fracture method (EDFM), for incorporating arbitrary 3D, discrete fractures, such as hydraulic fractures or faults, into modeling fracture-dominated fluid flow and heat transfer in fractured geothermal reservoirs. This technique allows 3D discrete fractures to be discretized independently from surrounding rock volume and inserted explicitly into a primary fracture/matrix grid, generated without including 3D discrete fractures in prior. An effective computational algorithm is developed to discretize these 3D discrete fractures and construct local connections between 3D fractures and fracture/matrix grid blocks representing the surrounding rock volume. The constructed gridding information on 3D fractures is then added to the primary grid. This embedded fracture modeling approach can be directly implemented into a developed geothermal reservoir simulator via the integral finite difference (IFD) method or with TOUGH2 technology. This embedded fracture modeling approach is very promising and computationally efficient to handle realistic 3D discrete fractures with complicated geometries, connections, and spatial distributions. Compared with other fracture modeling approaches, it avoids cumbersome 3D unstructured, local refining procedures, and increases computational efficiency by simplifying Jacobian matrix size and sparsity, while maintaining enough accuracy. Several numeral simulations are presented to demonstrate the utility and robustness of the proposed technique. Our numerical experiments show that this approach captures all the key patterns about fluid flow and heat transfer dominated by fractures in these cases. Thus, this approach is readily available to the simulation of fractured geothermal reservoirs with both artificial and natural fractures.

15 GEOTHERMAL ENERGY↗

DOE OSTI2025

Magnetically Driven Neutron-rich Ejecta Unleashed: Global 3D Neutrino–General Relativistic Magnetohydrodynamic Simulations of Collapsars Probe the Conditions for r -process Nucleosynthesis

Abstract Collapsars—rapidly rotating stellar cores that form black holes—can power gamma-ray bursts and are proposed to be key contributors to the production of heavy elements in the Universe via the rapid neutron capture process ( r -process). Previous neutrino-transport collapsar simulations have been unable to unbind neutron-rich material from the disk. However, these simulations have not included sufficiently strong magnetic fields and the black hole (BH), both of which are essential for launching mass outflows. We present ν h-amr , a novel neutrino-transport general relativistic magnetohydrodynamic ( ν GRMHD) code, which we use to perform the first 3D global ν GRMHD collapsar simulations. We find a self-consistent formation of a weakly magnetized dense accretion disk, which has sufficient time to neutronize. Eventually, substantial magnetic flux accumulates near the BH, becomes dynamically important, leads to a magnetically arrested disk (MAD), and unbinds some of the neutron-rich material. However, the strong flux also hinders accretion, lowers density, and increases neutrino-cooling timescale, which prevents further disk neutronization. Typical collapsar progenitors with mass accretion rates, M ̇ ∼ 0.1 − 1 M ⊙ s - 1 , do not produce significant neutron-rich ( Y e < 0.25) ejecta. However, we find that MADs at higher mass accretion rates, M ̇ ≳ few M ⊙ s - 1 (e.g., for more centrally concentrated progenitors), can unbind M ej ≲ M ⊙ of neutron-rich ejecta. The outflows inflate a shocked cocoon that mixes with the infalling neutron-poor stellar gas and raises the final outflow Y e ; however, the final r -process yield may be determined earlier at the point of neutron capture freeze-out. Future work will explore under what conditions more typical collapsar engines become r -process factories.

Issa, Danat (ORCID:0009000524787631)↗

DOE OSTI2016

Global MHD simulations of accretion disks in cataclysmic variables. I. The importance of spiral shocks

We present results from the first global 3D MHD simulations of accretion disks in cataclysmic variable (CV) systems in order to investigate the relative importance of angular momentum transport via turbulence driven by the magnetorotational instability (MRI) compared with that driven by spiral shock waves. Remarkably, we find that even with vigorous MRI turbulence, spiral shocks are an important component of the overall angular momentum budget, at least when temperatures in the disk are high (so that Mach numbers are low). In order to understand the excitation, propagation, and damping of spiral density waves in our simulations more carefully, we perform a series of 2D global hydrodynamical simulations with various equation of states, both with and without mass inflow via the Lagrangian point (L1). Compared with previous similar studies, we find the following new results. (1) The linear wave dispersion relation fits the pitch angles of spiral density waves very well. (2) We demonstrate explicitly that mass accretion is driven by the deposition of negative angular momentum carried by the waves when they dissipate in shocks. (3) Using Reynolds stress scaled by gas pressure to represent the effective angular momentum transport rate α{sub eff} is not accurate when mass accretion is driven by non-axisymmetric shocks. (4) Using the mass accretion rate measured in our simulations to directly measure α defined in standard thin-disk theory, we find 0.02≲α{sub eff}≲0.05 for CV disks, consistent with observed values in quiescent states of dwarf novae. In this regime, the disk may be too cool and neutral for the MRI to operate and spiral shocks are a possible accretion mechanism. However, we caution that our simulations use unrealistically low Mach numbers in this regime and, therefore, future models with more realistic thermodynamics and non-ideal MHD are warranted.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2019

Cluster Analysis of IRIS Spectroscopic Line Profiles and SDO/AIA EUV Emission in Observations and RMHD Simulations of the Solar Atmosphere

Spatially-resolved observations from the IRIS and SDO/AIA satellites, especially when coupled with realistic 3D RMHD simulations, are a powerful tool for analysis of processes in the solar chromosphere, transition region, and corona. However, the complexity of the data makes understanding the observations and modeling results difficult. In this work, we apply unsupervised clustering algorithms for analysis of observational and synthetic chromospheric Mg II h&k 2796Å&2803Å and transition region C II 1334Å&1335Å line profiles observed by IRIS, and extreme ultraviolet (EUV) emission observed by SDO/AIA, for various types of problems. The synthetic line profiles are computed for simulations of the quiescent solar atmosphere (using the StellarBox and RH1.5 codes). The K-Means clustering algorithm is applied, and the selection of an optimal number of clusters is supported by the average silhouette width technique. We discuss applications of the line profile clustering method to 1) visualization of computational and observational spectroscopic imaging data; 2) understanding of evolutionary trends and behavior patterns of quiet Sun emission and during solar flares; and 3) recognition of heating events and shock waves.

Sadykov, Viacheslav↗

NASA NTRS2022

Interrogating Solar Flare Loop Models With IRIS Observations 1: Overview of the Models, and Mass Flows

Solar flares are transient yet dramatic events in the atmosphere of the Sun, during which a vast amount of magnetic energy is liberated. This energy is subsequently transported through the solar atmosphere or into the heliosphere, and together with coronal mass ejections flares comprise a fundamental component of space weather. Thus, understanding the physical processes at play in flares is vital. That understanding often requires the use of forward modelling in order to predict the hydrodynamic and radiative response of the solar atmosphere. Those predictions must then be critiqued by observations to show us where our models are missing ingredients. While flares are of course 3D phenomenon, simulating the flaring atmosphere including an accurate chromosphere with the required spatial scales in 3D is largely beyond current computational capabilities, and certainly performing parameter studies of energy transport mechanisms is not yet tractable in 3D. Therefore, field-aligned 1D loop models that can resolve the relevant scales have a crucial role to play in advancing our knowledge of flares. In recent years, driven in part by the spectacular observations from the Interface Region Imaging Spectrograph (IRIS), flare loop models have revealed many interesting features of flares. For this review I highlight some important results that illustrate the utility of attacking the problem of solar flares with a combination of high quality observations, and state-of-the-art flare loop models, demonstrating: 1) how models help to interpret flare observations from IRIS, 2) how those observations show us where we are missing physics from our models, and 3) how the ever increasing quality of solar observations drives model improvements. Here in Paper one of this two part review I provide an overview of modern flare loop models, and of electron-beam driven mass flows during solar flares.

solar flares↗