Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation framework”

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.

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

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Image-Driven Hybrid Structural Analysis Based on Continuum Point Cloud Method with Boundary Capturing Technique

Conventional approaches for the structural health monitoring of infrastructures often rely on physical sensors or targets attached to structural members, which require considerable preparation, maintenance, and operational effort, including continuous on-site adjustments. This paper presents an image-driven hybrid structural analysis technique that combines digital image processing (DIP) and regression analysis with a continuum point cloud method (CPCM) built on a particle-based strong formulation. Polynomial regressions capture the boundary shape change due to the structural loading and precisely identify the edge and corner coordinates of the deformed structure. The captured edge profiles are transformed into essential boundary conditions. This allows the construction of a strongly formulated boundary value problem (BVP), classified as the Dirichlet problem. Capturing boundary conditions from the digital image is novel, although a similar approach was applied to the point cloud data. It was shown that the CPCM is more efficient in this hybrid simulation framework than the weak-form-based numerical schemes. Unlike the finite element method (FEM), it can avoid aligning boundary nodes with regression points. A three-point bending test of a rubber beam was simulated to validate the developed technique. The simulation results were benchmarked against numerical results by ANSYS and various relevant numerical schemes. The technique can effectively solve the Dirichlet-type BVP, yielding accurate deformation, stress, and strain values across the entire problem domain when employing a linear strain model and increasing the number of CPCM nodes. In addition, comparative analysis with conventional displacement tracking techniques verifies the developed technique’s robustness. The proposed technique effectively circumvents the inherent limitations of traditional monitoring methods resulting from the reliance on physical gauges or target markers so that a robust and non-contact solution for remote structural health monitoring in real-scale infrastructures can be provided, even in unfavorable experimental environments.

Chemistry↗
DOE OSTI2024

Electron energization in reconnection: Eulerian vs Lagrangian perspectives

Particle energization due to magnetic reconnection is an important unsolved problem for myriad space and astrophysical plasmas. Electron energization in magnetic reconnection has traditionally been examined from a particle, or Lagrangian, perspective using particle-in-cell (PIC) simulations. Guiding-center analyses of ensembles of PIC particles have suggested that Fermi (curvature drift) acceleration and direct acceleration via the reconnection electric field are the primary electron energization mechanisms. However, both PIC guiding-center ensemble analyses and spacecraft observations are performed in an Eulerian perspective. For this work, we employ the continuum Vlasov–Maxwell solver within the Gkeyll simulation framework to reexamine electron energization from a kinetic continuum, Eulerian, perspective. We separately examine the contribution of each drift energization component to determine the dominant electron energization mechanisms in a moderate guide-field Gkeyll reconnection simulation. In the Eulerian perspective, we find that the diamagnetic and agyrotropic drifts are the primary electron energization mechanisms away from the reconnection x-point, where direct acceleration dominates. We compare the Eulerian (Vlasov Gkeyll) results with the wisdom gained from Lagrangian (PIC) analyses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Advancing coupled fire-atmosphere research with HIGRAD/FIRETEC (w20_firetec) (CY 2021 Institutional Computing Annual Progress Report) [Slides]

LANL’s HIGRAD and FIRETEC models were used in a multi-scale simulation framework to investigate impacts of small-scale fuel heterogeneities on large-scale Pyro-Cb formation. Forest structural diversity is found to increase canopy fire spread, and thus risk of tree mortality in FIRETEC simulations. Structural diversity is associated with an increase in ladder fuels, which facilitate fire spread into the canopy. Low-intensity fire burning in marginal conditions is highly sensitive to small perturbation in atmospheric turbulence. These results highlight predictive limitations of individual fire simulations, and emphasize the need for ensemble simulations, especially for low-intensity, prescribed fire.

36 MATERIALS SCIENCE↗
NASA NTRS2005

Java-based Graphical User Interface for MAVERIC-II

A computer program entitled "Marshall Aerospace Vehicle Representation in C II, (MAVERIC-II)" is a vehicle flight simulation program written primarily in the C programming language. It is written by James W. McCarter at NASA/Marshall Space Flight Center. The goal of the MAVERIC-II development effort is to provide a simulation tool that facilitates the rapid development of high-fidelity flight simulations for launch, orbital, and reentry vehicles of any user-defined configuration for all phases of flight. MAVERIC-II has been found invaluable in performing flight simulations for various Space Transportation Systems. The flexibility provided by MAVERIC-II has allowed several different launch vehicles, including the Saturn V, a Space Launch Initiative Two-Stage-to-Orbit concept and a Shuttle-derived launch vehicle, to be simulated during ascent and portions of on-orbit flight in an extremely efficient manner. It was found that MAVERIC-II provided the high fidelity vehicle and flight environment models as well as the program modularity to allow efficient integration, modification and testing of advanced guidance and control algorithms. In addition to serving as an analysis tool for techno logy development, many researchers have found MAVERIC-II to be an efficient, powerful analysis tool that evaluates guidance, navigation, and control designs, vehicle robustness, and requirements. MAVERIC-II is currently designed to execute in a UNIX environment. The input to the program is composed of three segments: 1) the vehicle models such as propulsion, aerodynamics, and guidance, navigation, and control 2) the environment models such as atmosphere and gravity, and 3) a simulation framework which is responsible for executing the vehicle and environment models and propagating the vehicle s states forward in time and handling user input/output. MAVERIC users prepare data files for the above models and run the simulation program. They can see the output on screen and/or store in files and examine the output data later. Users can also view the output stored in output files by calling a plotting program such as gnuplot. A typical scenario of the use of MAVERIC consists of three-steps; editing existing input data files, running MAVERIC, and plotting output results.

Seo, Suk Jai↗
NASA NTRSDate not supplied

Advanced Materials for the Lunar Surface: Multiscale Computational Design of Refractory Alloys and Carbides

Emerging operational environments, such as the lunar surface, present novel challenges for NASA and drive the need for advanced materials in applications like fission surface power systems. To address these demands, computational materials science is rapidly evolving to augment or replace costly and hazardous empirical testing. Although materials selection at NASA remains predominantly experimentally driven, advanced simulation methodologies are being steadily integrated into the engineering lifecycle. This work details the application of multiscale simulation techniques—including first-principles calculations, CALPHAD, dislocation dynamics, and molecular dynamics—at NASA's Ames Research Center to evaluate advanced materials for extreme environments. First, we present contributions to the Space Nuclear Propulsion Project. Be-cause propellant channel coatings in nuclear thermal rockets must withstand high-pressure, high-temperature hydro-gen, optimizing these materials is critical. First-principles calculations were employed to establish a rigorous quantitative and qualitative understanding of the behavior of the refractory carbides ZrC, NbC, and their mixtures in high-enthalpy hydrogen environments. This necessitated the generation of high-fidelity thermodynamic models for both stoichiometric and carbon-depleted carbides, both with and without the presence of hydrogen. Furthermore, we highlight efforts under the Refractory Alloy Additive Manufacturing Build Optimization (RAAMBO) project, where existing and novel alloy compositions were assessed for additive manufacturing printability and subsequent performance in applications such as heat pipes and rocket nozzle extensions. This was accomplished through a comprehensive multiscale simulation framework that bridged the gap from the nanometer to the millimeter scale. Across both initiatives, rigorous validation against empirical data was prioritized. By systematically employing a verified and validated computational frame-work, we demonstrate how simulation effectively supports multidisciplinary engineering efforts, builds project-wide confidence, and drives critical materials development.

computational materials↗
DOE OSTI2022

Simulation study of fast ion losses associated with the rotating n = 1 resonant magnetic perturbations in KSTAR

In this work we describe a simulation framework for testing the fast ion loss mechanism associated with the experimentally observed beam ion losses when an externally applied toroidally rotating perturbed magnetic field is used to control edge localized modes in the KSTAR tokamak. The simulations reproduce the key qualitative features of neutral beam injection (NBI) ion detection by a fast ion loss detector. The NBI ion losses in the simulation mainly occur for passing particles due to orbit stochastization, which is caused by orbit resonance with the 3D field perturbations. The relative toroidal angle of NBI ion deposition to the orbit island determines the radial path of the particles between confinement or loss. The fast ion loss quantity depends on the toroidal angle rotation of the 3D field with respect to the beam deposition position. The averaged transport of toroidal canonical angular momentum reveals that linear and nonlinear resonance of the NBI ions within the 3D field are the main factors determining fast ion transport and loss.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Machine Protection Studies of Fast Instability-Driven Beam Losses in the EIC Electron Storage Ring

Fast transverse beam instabilities represent one of the most challenging protection scenarios for the Electron-Ion Collider (EIC) electron storage ring (ESR). Coherent betatron oscillations driven by collective effects can grow rapidly and lead to substantial beam losses within only a few tens of turns, posing a threat to collimators, superconducting (SC) magnets, cryogenic systems, and detector components. A dedicated simulation framework based on distributed transverse excitations was developed to evaluate instability-driven beam losses in the ESR. A conservative instability rise time on the order of 10 turns was used to define a bounding protection scenario. Simulations show that the original two-collimator configuration is insufficient to localize losses under all operating conditions, whereas an optimized four-collimator configuration reduces loss leakage into protected regions by more than two orders of magnitude. Nevertheless, the most severe instability scenarios produce catastrophic beam impacts, with nearly half of the stored beam being intercepted by the collimators within a single turn, exceeding the preliminary estimated survivability limits of the baseline collimator design. The results demonstrate that collimation alone is insufficient to guarantee safe ESR operation and must be complemented by a fast machine-protection system. Dedicated turn-by-turn beam-position monitors and fast beam-loss monitors provide sufficient advance warning to detect instability growth and initiate a beam abort before equipment-damaging loss levels are reached. These findings establish quantitative requirements for the ESR machine-protection architecture and support the implementation of a multilayer protection strategy for both accelerator and detector systems.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Regional Subsidence Analysis Through a Multi-Scale Modeling Framework Based on Breakage Mechanics

Although poroelastic models are often used to explain the delay between subsurface fluid depletion and ground subsidence, inelastic compaction involving permanent changes of rock microstructure may exacerbate hydro-mechanical coupling, thus influencing the interpretation of measurements and long-term forecasts. Here, a multi-scale modeling approach is discussed, which accounts for the inherent connection between rock microstructure, hydraulic conductivity, and pore compaction. A constitutive model built within the framework of breakage mechanics is proposed to link the hydraulic conductivity of granular rocks with inelastic deformations and changes in grading caused by injection-depletion cycles at stress levels far from yielding. The proposed model has been incorporated into large-scale simulation frameworks, thus enabling the spatiotemporal mapping of regional subsidence through a hybrid, semi-analytical approach. Numerical results based on this strategy show that the model allows isolating near-field and far-field effects into the computation of land subsidence and can generate forecasts for different modeling scenarios (e.g., elastic and inelastic compaction, constant permeability, and concurrent change of compressibility and permeability). In particular, examples of simulations for the case of the Groningen gas field are discussed, showing the model capabilities to use both field measurements and laboratory tests for the generation of reasonable subsidence maps, without expensive computational costs. Results indicate that ignoring coupled inelastic effects has major consequences on the predicted timescale of subsidence. Specifically, while all the model scenarios produced similar long-term ground settlements, those ignoring breakage-dependent permeability changes result in a variation of the temporal window of residual subsidence of the order of several decades.

58 GEOSCIENCES↗
DOE OSTI2024

COSMIC DAWN: Distributed Analysis of Wireless at Nextscale

Distributed Analysis of Wireless at Nextscale (DAWN) is a novel simulation framework for large-scale design-space exploration (DSE) of unmodified software-defined radio (SDR) applications interacting in a scalable, high-fidelity, virtual physics environment. The software-defined nature of the coupled software-physics simulation leverages hardware emulation to permit in-depth examination and modification of not only the electromagnetic environment, including each signal in flight, but also the precise state of system software and components. DAWN supports modular, customizable physics environments allowing realistic propagation effects so that computationally efficient empirical models, reduced order/surrogate models, or large-scale, high-fidelity, site-specific simulations can be used as a propagation medium based on scenario requirements. This paper introduces DAWN’s design and initial implementation, detailing key architectural components, including the Physics Realization Engine (PhyRE), Runtime Infrastructure for Simulation Environments (RISE), and the design space exploration (DSE) suite. It concludes with demonstrations using unmodified 4G/LTE software available from srsRAN on computing resources ranging from a small cluster to ORNL’s Frontier Exascale system.

Wise, Mike [ORNL] (ORCID:0000000266120641)↗
NASA NTRSDate not supplied

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A Liever↗
NASA NTRSDate not supplied

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A. Liever↗
NASA NTRS2006

AEROELASTIC SIMULATION TOOL FOR INFLATABLE BALLUTE AEROCAPTURE

A multidisciplinary analysis tool is under development for predicting the impact of aeroelastic effects on the functionality of inflatable ballute aeroassist vehicles in both the continuum and rarefied flow regimes. High-fidelity modules for continuum and rarefied aerodynamics, structural dynamics, heat transfer, and computational grid deformation are coupled in an integrated multi-physics, multi-disciplinary computing environment. This flexible and extensible approach allows the integration of state-of-the-art, stand-alone NASA and industry leading continuum and rarefied flow solvers and structural analysis codes into a computing environment in which the modules can run concurrently with synchronized data transfer. Coupled fluid-structure continuum flow demonstrations were conducted on a clamped ballute configuration. The feasibility of implementing a DSMC flow solver in the simulation framework was demonstrated, and loosely coupled rarefied flow aeroelastic demonstrations were performed. A NASA and industry technology survey identified CFD, DSMC and structural analysis codes capable of modeling non-linear shape and material response of thin-film inflated aeroshells. The simulation technology will find direct and immediate applications with NASA and industry in ongoing aerocapture technology development programs.

Liever, P. A.↗
NASA NTRSDate not supplied

Generic Urban Air Mobility Simulation

This research presents a simulation framework for autonomous research for a UAM vehicle using the NASA Revolutionary Vertical Lift Technology Lift+Cruise concept vehicle. Our research results were produced using the open-source, six degree of freedom, rigid-body, nonlinear generic urban air mobility (GUAM) simulation. The intent of this paper is to demonstrate the GUAM simulation and a series of Challenge Problems that our researchers have posed to the broader autonomous vehicle research community. Our team has developed the GUAM simulation for the express purpose of providing a high-fidelity transition vehicle dynamics model to foster collaboration and algorithm performance comparison across research teams. In this paper, we demonstrate some of the autonomous flight research challenges and some of our current approaches to tackling basic autonomous flight tasks (e.g., trajectory following, stationary and moving obstacle avoidance). Additionally, we propose some flight metrics to assess autonomous algorithm performance while accomplishing these basic autonomous tasks.

autonomous flight↗
DOE OSTI2025

Numerical Simulation Studies of Ultrasonic De-Icing for Heating, Ventilation, Air Conditioning, and Refrigeration Structures

Ice accumulation on heating, ventilation, air conditioning, and refrigeration (HVACR) structures presents significant operational challenges. These challenges include reduced efficiency, increased energy consumption, and potential damage to equipment. Traditional de-icing methods, such as chemical treatments, mechanical scraping, or heating-based techniques, are often labor-intensive, costly, and environmentally harmful. Here, this study uniquely investigates ultrasonic de-icing as an energy-efficient alternative for HVACR applications, focusing on the specific structural geometries found in these systems. A comprehensive numerical simulation framework was developed using finite element analysis to explore ultrasonic wave propagation across four distinct HVACR structures. Key parameters such as ultrasonic frequency, power levels, and the number and placement of actuators were examined for their impact on ice detachment efficiency. Results from simulations on a plate structure reveal that ultrasonic excitation can propagate effectively across large areas (at least 150 × 150 mm), enhancing the de-icing coverage. Lower frequency (e.g., 30 to 45 kHz) excitation results in greater displacement, improving de-icing performance, while increased actuator numbers with the same total power input also enhance effectiveness. Two actuators seem sufficient for the de-icing of a 300 × 300 mm plate. For tube-and-fin structures, specific high-power ultrasonic frequencies selectively excite the fin plates, demonstrating efficient ice removal when actuated on the tube. However, optimal performance requires careful design of actuator placement and vibration modes to accommodate the irregular shapes of these structures.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Computational studies of impurity migration during induction stirring of molten uranium

Understanding and controlling impurity behavior in actinide metal casting processes are foundational for efficient part production yet remain major challenges for researchers and industry. To help provide insight regarding impurity distribution during actinide metal casting, we have developed computational fluid dynamics (CFD) models for a laboratory-scale system using commercial and open-source code. Multiple simulation frameworks allow for improved confidence in the resulting outputs, while taking advantage of the maturity and convenience offered by commercial platforms and simultaneously maintaining the transparency and flexibility often provided by open-source software. Here, we describe multiple experiment-informed models designed to simulate a specific laboratory system in which uranium melt, containing a known starting concentration of carbon impurity, is electromagnetically stirred in an induction furnace. Here, the goal of the simulations is to predict the motion of uranium carbide microparticles in the velocity field of the melt. Prior to simulating the uranium-carbon system, numerical models were validated using a previously published nonradioactive experimental system. Effects of the size and shape of impurity particles in the models were investigated and agree with experimental findings. Simulation of smaller particles (< 50 µm) shows more homogenous distribution throughout the stirred melt. With increased particle size (100 µm), the body forces, which include the buoyancy force, dominate over the drag force, causing larger particles to move toward the crucible walls and upward in the system.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Roadmap on methods and software for electronic structure based simulations in chemistry and materials

This Roadmap article provides a succinct, comprehensive overview of the state of electronic structure methods and software for molecular and materials simulations. Seventeen distinct sections collect insights by 51 leading scientists in the field. Each contribution addresses the status of a particular area, as well as current challenges and anticipated future advances, with a particular eye towards software related aspects and providing key references for further reading. Foundational sections cover density functional theory and its implementation in real-world simulation frameworks, Green's function based many-body perturbation theory, wave-function based and stochastic electronic structure approaches, relativistic effects and semiempirical electronic structure theory approaches. Subsequent sections cover nuclear quantum effects, real-time propagation of the electronic structure, challenges for computational spectroscopy simulations, and exploration of complex potential energy surfaces. The final sections summarize practical aspects, including computational workflows for complex simulation tasks, the impact of current and future high-performance computing architectures, software engineering practices, education and training to maintain and broaden the community, as well as the status of and needs for electronic structure based modeling from the vantage point of industry environments. Overall, the field of electronic structure software and method development continues to unlock immense opportunities for future scientific discovery, based on the growing ability of computations to reveal complex phenomena, processes and properties that are determined by the make-up of matter at the atomic scale, with high precision.

36 MATERIALS SCIENCE↗
DOE OSTI2023

SILCC – VII. Gas kinematics and multiphase outflows of the simulated ISM at high gas surface densities

ABSTRACT We present magnetohydrodynamic (MHD) simulations of the star-forming multiphase interstellar medium (ISM) in stratified galactic patches with gas surface densities Σgas = 10, 30, 50, and 100 $\mathrm{M_\odot \, pc^{-2}}$. The silcc project simulation framework accounts for non-equilibrium thermal and chemical processes in the warm and cold ISM. The sink-based star formation and feedback model includes stellar winds, hydrogen-ionizing UV radiation, core-collapse supernovae, and cosmic ray (CR) injection and diffusion. The simulations follow the observed relation between Σgas and the star formation rate surface density ΣSFR. CRs qualitatively change the outflow phase structure. Without CRs, the outflows transition from a two-phase (warm and hot at 1 kpc) to a single-phase (hot at 2 kpc) structure. With CRs, the outflow always has three phases (cold, warm, and hot), dominated in mass by the warm phase. The impact of CRs on mass loading decreases for higher Σgas and the mass loading factors of the CR-supported outflows are of order unity independent of ΣSFR. Similar to observations, vertical velocity dispersions of the warm ionized medium (WIM) and the cold neutral medium (CNM) correlate with the star formation rate as $\sigma _\mathrm{z} \propto \Sigma _\mathrm{SFR}^a$, with a ∼ 0.20. In the absence of stellar feedback, we find no correlation. The velocity dispersion of the WIM is a factor ∼2.2 higher than that of the CNM, in agreement with local observations. For $\Sigma _\mathrm{SFR} \gtrsim 1.5 \times 10^{-2}\, \mathrm{M}_\odot \, \mathrm{yr}^{-1}\, \mathrm{kpc}^{-2}$ the WIM motions become supersonic.

Astronomy & Astrophysics↗