Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled 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.

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

At least 1,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Increased Asian Sulfate Aerosol Emissions Remarkably Enhance Sahel Summer Precipitation

Observational evidence shows that Sahel summer precipitation has experienced a considerable increase since the 1980s, coinciding with significant diverging trends of increased sulfate emissions in Asia and decreased emissions in Europe (dipole pattern of aerosols between Asia and Europe). The decrease in European sulfate aerosols has substantial effects on the Sahel summer precipitation increase, but the corresponding effect of increased Asian sulfate is unknown. Multi-model simulations in the Precipitation Driver and Response Model Intercomparison Project (PDRMIP) show, compared to decreased European aerosols, that increased Asian aerosols similarly enhance the Sahel summer precipitation but with different large-scale atmospheric circulation changes. Further analysis of the Sixth Coupled Model Intercomparison Project (CMIP6) simulations under historical attribution and various emission scenarios reinforces the results about the climate impacts of anthropogenic aerosols and suggests that in future scenarios with strong international cooperation and rapid climate mitigations (SSP2-45), the Sahel drought will be intensified likely due to the decline in Asian aerosol emissions. Our results suggest that Asian anthropogenic aerosols are likely a non-negligible driver of the recent recovery in Sahel precipitation amounts

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1974

Flight-path and airspeed control during landing approach for powered-lift aircraft

Manual control of flight path and airspeed during landing approach has been investigated for powered-lift transport aircraft. An analysis was conducted to identify the behavior of the aircraft which would be potentially significant to the pilot controlling flight path and airspeed during the approach. The response characteristics found to describe the aircraft behavior were (1) the initial flight-path response and flight-path overshoot for a step change in thrust, (2) the steady-state coupling of flight path and airspeed for a step change in thrust, and (3) the sensitivity of airspeed to changes in pitch attitude. The significance of these response characteristics was evaluated by pilots on a large-motion, ground-based simulator at Ames Research Center. Coupling between flight path and airspeed was considered by the pilot to be the dominant influence on handling qualities for the approach task. Results are compared with data obtained from flight tests of three existing powered-lift V/STOL aircraft.

Franklin, J. A.↗
NASA NTRSDate not supplied

A Methodology for Actuating RCS Jets in a Continuous, Time-Accurate CFD Simulation

Computational fluid dynamics simulations of Mars entry vehicle reaction control systems have been used to gain a better understanding of the wide range of flow phenomena encountered during various phases of flight. One aspect of this is simulating the system during controlled flight. Current practices model the system in a continuously actuated state, but this methodology is unable to accurately model the dynamic behavior of the system during EDL. One method to improve this is to create a coupled CFD-RBD-Control simulation. The first step to this is establishing a method of accurately actuating an RCS jet during a time accurate CFD simulation. This paper presents a method to achieve this actuation. There are three phases to the methodology: establishing the off state, creating a ramping function that will initialize the jet flow, and establishing the on state. A test case of a jet in a supersonic crossflow is implemented to establish the success of the methodology in being able to actuate a jet in the midst of an ongoing CFD simulation. It shows significant improvement in boundary condition responsiveness for linear and exponential ramping functions in comparison to the step function. A preliminary demonstration of the methodology is presented for the MSL vehicle which shows significant advantages for the implementation of a linear ramping function over a step function.

Alexandra Hickey↗
NASA NTRSDate not supplied

A Methodology for Actuating RCS Jets in a Continuous, Time-Accurate CFD Simulation

Computational fluid dynamics simulations of Mars entry vehicle reaction control systems have been used to gain a better understanding of the wide range of flow phenomena encountered during various phases of flight. One aspect of this is simulating the system during controlled flight. Current practices model the system in a continuously actuated state, but this methodology is unable to accurately model the dynamic behavior of the system during EDL. One method to improve this is to create a coupled CFD-RBD-Control simulation. The first step to this is establishing a method of accurately actuating an RCS jet during a time accurate CFD simulation. This paper presents a method to achieve this actuation. There are three phases to the methodology: establishing the off state, creating a ramping function that will initialize the jet flow, and establishing the on state. A test case of a jet in a supersonic crossflow is implemented to establish the success of the methodology in being able to actuate a jet in the midst of an ongoing CFD simulation. It shows significant improvement in boundary condition responsiveness for linear and exponential ramping functions in comparison to the step function. A preliminary demonstration of the methodology is presented for the MSL vehicle which shows significant advantages for the implementation of a linear ramping function over a step function.

Alexandra Hickey↗
DOE OSTI2021

Coupling of regional geophysics and local soil-structure models in the EQSIM fault-to-structure earthquake simulation framework

Accurate understanding and quantification of the risk to critical infrastructure posed by future large earthquakes continues to be a very challenging problem. Earthquake phenomena are quite complex and traditional approaches to predicting ground motions for future earthquake events have historically been empirically based whereby measured ground motion data from historical earthquakes are homogenized into a common data set and the ground motions for future postulated earthquakes are probabilistically derived based on the historical observations. This procedure has recognized significant limitations, principally due to the fact that earthquake ground motions tend to be dictated by the particular earthquake fault rupture and geologic conditions at a given site and are thus very site-specific. Historical earthquakes recorded at different locations are often only marginally representative. There has been strong and increasing interest in utilizing large-scale, physics-based regional simulations to advance the ability to accurately predict ground motions and associated infrastructure response. However, the computational requirements for simulations at frequencies of engineering interest have proven a major barrier to employing regional scale simulations. In a U.S. Department of Energy Exascale Computing Initiative project, the EQSIM application development is underway to create a framework for fault-to-structure simulations. This framework is being prepared to exploit emerging exascale platforms in order to overcome computational limitations. This article presents the essential methodology and computational workflow employed in EQSIM to couple regional-scale geophysics models with local soil-structure models to achieve a fully integrated, complete fault-to-structure simulation framework. Here, the computational workflow, accuracy and performance of the coupling methodology are illustrated through example fault-to-structure simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Toward full simulations for a liquid metal blanket: part 2. Computations of MHD flows with volumetric heating for a PbLi blanket prototype at Ha ~10 4 and Gr ~10 12

On the pathway toward full simulations for a liquid metal (LM) blanket, this part 2 extends a previous study of purely magnetohydrodynamic (MHD) flows in a DCLL blanket in reference Chen et al (2020 Nucl. Fusion 60 076003) to more general conditions when the MHD flow is coupled with heat transfer. The simulated prototypic blanket module includes all components of a real LM blanket system, such as supply ducts, inlet and outlet manifolds, multiple poloidal ducts and a U-turn zone. Volumetric heating generated by fusion neutrons is added to simulate thermal effects in the flowing lead–lithium (PbLi) breeder. The MHD flow equations and the energy equation are solved with a DNS-type finite-volume code ‘MHD-UCAS’ on a very fine mesh of 470 × 10 6 cells. The applied magnetic field is 5 T (Hartmann number Ha ~ 10 4 ), the PbLi velocity in the poloidal ducts is 10 cm s –1 (Reynolds number Re ~10 5 ), whereas the maximum volumetric heating is 30 MW m –3 (Grashof number Gr ~ 10 12 ). Four cases have been simulated, including forced- and mixed-convection flows, and either an electrically conducting or insulating blanket structure. Various comparisons are made between the four computed cases and also against the purely MHD flows computed earlier in reference Chen et al (2020 Nucl. Fusion 60 076003) with regards to the (1) MHD pressure drop, (2) flow balancing, (3) temperature field, (4) flows in particular blanket components, and (5) 3D and turbulent flow effects. The strongest buoyancy effects were found in the poloidal ducts. In the electrically non-conducting blanket, the buoyancy forces lead to significant modifications of the flow structure, such as formation of reverse flows, whereas their effect on the MHD pressure drop is relatively small. In the electrically conducting blanket case, the buoyancy effects on the flow and MHD pressure drop are almost negligible.

Physics↗
NASA NTRS2011

Influence of the Mt. Pinatubo Eruption on the Stratospheric Circulation

On June 15th, 1991 the eruption of Mt. Pinatubo in the Philippines injected about 20 Tg of sulfur dioxide in the stratosphere, which was transformed into sulfuric acid aerosol. Even though stratospheric winds climatologically tend to hinder the air mixing between the two hemispheres, observations have shown that a large part of the SO2 emitted by Mt. Pinatubo have been transported from the Northern to the Southern Hemisphere. We show how the absorption of radiation by sulfate aerosol is responsible for the spreading to the southern hemisphere through a middle stratospheric channel. We simulate the eruption of Mt. Pinatubo with the Goddard Earth Observing System (GEOS) version 5 general circulation model, coupled to the aerosol module GOCART and the stratospheric chemistry module StratChem. Our simulations are in good agreement with SAGE-II and AVHRR data. We perform two ensembles of simulations: the first ensemble consists of runs without coupling between aerosol and radiation. In these simulations the plume of aerosols is treated as a passive tracer and the atmosphere is unperturbed. In the second ensemble of simulations aerosols and radiation are coupled. We show that the set of runs with interactive aerosol produces a larger cross-equatorial transport of the Pinatubo cloud, in agreement with the observations. At first, the volcanic cloud is transported from the latitude of the eruption to both hemispheres through a lower stratospheric pathway. Additionally, in the interactive simulations the absorption of long wave radiation from the volcanic sulfate induces a lofting of the cloud to the middle atmosphere and, at the same time, a divergent motion from the center of the cloud. Such motion spreads the volcanic cloud across the equator and to the tropics, where the background circulation carry it to higher latitudes.

Aquila, V.↗
NASA NTRS2013

Stratospheric Ozone Response in Experiments G3 and G4 of the Geoengineering Model Intercomparison Project (GeoMIP)

Geoengineering with stratospheric sulfate aerosols has been proposed as a means of temporarily cooling the planet, alleviating some of the side effects of anthropogenic CO2 emissions. However, one of the known side effects of stratospheric injections of sulfate aerosols is a decrease in stratospheric ozone. Here we show results from two general circulation models and two coupled chemistry climate models that have simulated stratospheric sulfate aerosol geoengineering as part of the Geoengineering Model Intercomparison Project (GeoMIP). Changes in photolysis rates and upwelling of ozone-poor air in the tropics reduce stratospheric ozone, suppression of the NOx cycle increases stratospheric ozone, and an increase in available surfaces for heterogeneous chemistry modulates reductions in ozone. On average, the models show a factor 20-40 increase of the sulfate aerosol surface area density (SAD) at 50 hPa in the tropics with respect to unperturbed background conditions and a factor 3-10 increase at mid-high latitudes. The net effect for a tropical injection rate of 5 Tg SO2 per year is a decrease in globally averaged ozone by 1.1-2.1 DU in the years 2040-2050 for three models which include heterogeneous chemistry on the sulfate aerosol surfaces. GISS-E2-R, a fully coupled general circulation model, performed simulations with no heterogeneous chemistry and a smaller aerosol size; it showed a decrease in ozone by 9.7 DU. After the year 2050, suppression of the NOx cycle becomes more important than destruction of ozone by ClOx, causing an increase in total stratospheric ozone. Contribution of ozone changes in this experiment to radiative forcing is 0.23 W m-2 in GISS-E2-R and less than 0.1 W m-2 in the other three models. Polar ozone depletion, due to enhanced formation of both sulfate aerosol SAD and polar stratospheric clouds, results in an average 5 percent increase in calculated surface UV-B.

aerosol↗
DOE OSTI2023

A high-resolution large-eddy simulation framework for wildland fire predictions using TensorFlow

Background: Wildfires are becoming more severe, so we need improved tools to predict them over a wide range of conditions and scales. One approach towards this goal entails the use of coupled fire/atmosphere modelling tools. Although significant progress has been made in advancing their physical fidelity, existing tools have not taken full advantage of emerging programming paradigms and computing architectures to enable high-resolution wildfire simulations. Aims: The aim of this study was to present a new framework that enables landscape-scale wildfire simulations with physical representation of combustion at an affordable cost. Methods: We developed a coupled fire/atmosphere simulation framework using TensorFlow, which enables efficient and scalable computations on Tensor Processing Units. Key Results: Simulation results for a prescribed fire were compared with experimental data. Predicted fire behavior and statistical analysis for fire spread rate, scar area, and intermittency showed overall reasonable agreement. Scalability analysis was performed, showing close to linear scaling. Conclusions: While mesh refinement was shown to have less impact on global quantities, such as fire scar area and spread rate, it benefits predictions of intermittent fire behavior, buoyancy-driven dynamics, and small-scale turbulent motion. Implications: This new simulation framework is efficient in capturing both global quantities and unsteady dynamics of wildfires at high spatial resolutions.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2009

Self-Consistent Model of Magnetospheric Electric Field, Ring Current, Plasmasphere, and Electromagnetic Ion Cyclotron Waves: Initial Results

Further development of our self-consistent model of interacting ring current (RC) ions and electromagnetic ion cyclotron (EMIC) waves is presented. This model incorporates large scale magnetosphere-ionosphere coupling and treats self-consistently not only EMIC waves and RC ions, but also the magnetospheric electric field, RC, and plasmasphere. Initial simulations indicate that the region beyond geostationary orbit should be included in the simulation of the magnetosphere-ionosphere coupling. Additionally, a self-consistent description, based on first principles, of the ionospheric conductance is required. These initial simulations further show that in order to model the EMIC wave distribution and wave spectral properties accurately, the plasmasphere should also be simulated self-consistently, since its fine structure requires as much care as that of the RC. Finally, an effect of the finite time needed to reestablish a new potential pattern throughout the ionosphere and to communicate between the ionosphere and the equatorial magnetosphere cannot be ignored.

Gamayunov, K. V.↗
NASA NTRS2015

Electric Field Generation in Martian Dust Devils

Terrestrial dust devils are known to generate electric fields from the vertical separation of charged dust particles. The particles present within the dust devils on Mars may also be subject to similar charging processes and so likely contribute to electric field generation there as well. However, to date, no Marsin situ instrumentation has been deployed to measure electric field strength. In order to explore the electric environment of dust devils on Mars, the triboelectric dust charging physics from the MacroscopicTriboelectric Simulation (MTS) code has been coupled to the Mars Regional Atmospheric ModelingSystem (MRAMS). Using this model, we examine how macroscopic electric fields are generated within martian dust disturbances and attempt to quantify the time evolution of the electrodynamical system.Electric fields peak for several minutes within the dust devil simulations. The magnitude of the electric field is a strong function of the size of the particles present, the average charge on the particles and the number of particles lifted. Varying these parameters results in peak electric fields between tens of millivolts per meter and tens of kilovolts per meter.

Aeolian processes↗
DOE OSTI2024

Effects of hydrogen isotope species on ITG microturbulence in LHD

Abstract The linear and nonlinear effects of hydrogen isotope species on ion temperature gradient (ITG) instability in the Large Helical Device (LHD) stellarator are studied using radially global gyrokinetic simulation. We found that the coupling range of linear toroidal harmonics depends on the ion mass of the hydrogen isotope. The growth rates of ITG mode are almost the same for H, D, and T plasmas, indicating a gyro-Bohm scaling of ion-mass dependence. The nonlinear electrostatic simulations show that the zonal flow breaks the radially elongated eigenmode structures and reduces the size of the turbulence eddies, which suppresses the turbulence and the ion heat transport in the LHD. The turbulence amplitude without the zonal flow is almost the same for H, D, and T plasmas, while it decreases with increasing ion mass of the hydrogen isotope when the zonal flow is present. The reduction of the turbulent transport with larger ion mass is mostly due to the enhancement of zonal flows by larger ion mass. The ion heat conductivity deviates from the gyro-Bohm scaling for both cases with and without the zonal flow.

Physics↗
DOE OSTI2022

Quantum Markov chain Monte Carlo with digital dissipative dynamics on quantum computers

Modeling the dynamics of a quantum system connected to the environment is critical for advancing our understanding of complex quantum processes, as most quantum processes in nature are affected by an environment. Modeling a macroscopic environment on a quantum simulator may be achieved by coupling independent ancilla qubits that facilitate energy exchange in an appropriate manner with the system and mimic an environment. This approach requires a large, and possibly exponential number of ancillary degrees of freedom which is impractical. In contrast, we develop a digital quantum algorithm that simulates interaction with an environment using a small number of ancilla qubits. By combining periodic modulation of the ancilla energies, or spectral combing, with periodic reset operations, we are able to mimic interaction with a large environment and generate thermal states of interacting many-body systems. We evaluate the algorithm by simulating preparation of thermal states of the transverse Ising model. Our algorithm can also be viewed as a quantum Markov chain Monte Carlo process that allows sampling of the Gibbs distribution of a multivariate model. To illustrate this we evaluate the accuracy of sampling Gibbs distributions of simple probabilistic graphical models using the algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Deterministic modelling of the SPERT IV reactor transients using the multi-physics capability of WIMS

The ANSWERS WIMS reactor physics code is being developed for whole core multi-physics modelling. The established neutronics capability for lattice calculations has recently been extended to be suitable for whole core modelling of Small Modular Reactors (SMRs). A whole core transport, SP3 or diffusion flux solution is combined with fuel assembly resonance shielding. An integrated thermal hydraulic solver permits temperature and density variations to feedback to the neutronics calculation. This capability can be applied to both steady state and time dependent transient problems. Nuclear reactor design and safety case development requires assessment of a range of reactor transients, to inform both normal operation limitations and accident scenario analysis. This paper presents new methodology developed in WIMS to couple the core neutronics to the integrated core thermal hydraulics solver for the simulation of reactor transients in whole core models. To support the validation of the multi-physics capability of WIMS, this capability has been applied to the reactivity insertion transient experiments performed in the plate fuelled SPERT-IV reactor. This study employs WIMS using the whole core solver MERLIN, which calculates the time-dependent flux distribution and magnitude, coupled to ARTHUR, which solves for the time-dependent thermal hydraulics solution, and provides thermal feedback via the cross sections generated by GEOM, which performs resonance shielding calculations and generates cross section data for the plate geometry. This spatial kinetics model, with dynamic cross section generation, allows for variations in the temperature and neutron flux profile with time.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2009

WRF Simulations of the 20-22 January 2007 Snow Events over Eastern Canada: Comparison with In-Situ and Satellite Observations

One of the grand challenges of the Global Precipitation Measurement (GPM) mission is to improve cold season precipitation measurements in middle and high latitudes through the use of high-frequency passive microwave radiometry. For this, the Weather Research and Forecasting (WRF) model with the Goddard microphysics scheme is coupled with a satellite data simulation unit (WRF-SDSU) that has been developed to facilitate over-land snowfall retrieval algorithms by providing a virtual cloud library and microwave brightness temperature (Tb) measurements consistent with the GPM Microwave Imager (GMI). This study tested the Goddard cloud microphysics scheme in WRF for two snowstorm events, a lake effect and a synoptic event, that occurred between 20 and 22 January 2007 over the Canadian CloudSAT/CALIPSO Validation Project (C3VP) site in Ontario, Canada. The 24h-accumulated snowfall predicted by the WRF model with the Goddard microphysics was comparable to the observed accumulated snowfall by the ground-based radar for both events. The model correctly predicted the onset and ending of both snow events at the CARE site. WRF simulations capture the basic cloud properties as seen by the ground-based radar and satellite (i.e., CloudSAT, AMSU-B) observations as well as the observed cloud streak organization in the lake event. This latter result reveals that WRF was able to capture the cloud macro-structure reasonably well.

Shi, J. J.↗
DOE OSTI2021

Climate Change Projection in the Twenty-First Century Simulated by NIMS-KMA CMIP6 Model Based on New GHGs Concentration Pathways

Abstract The National Institute of Meteorological Sciences-Korea Meteorological Administration (NIMS-KMA) has participated in the Coupled Model Inter-comparison Project (CMIP) and provided long-term simulations using the coupled climate model. The NIMS-KMA produces new future projections using the ensemble mean of KMA Advanced Community Earth system model (K-ACE) and UK Earth System Model version1 (UKESM1) simulations to provide scientific information of future climate changes. In this study, we analyze four experiments those conducted following the new shared socioeconomic pathway (SSP) based scenarios to examine projected climate change in the twenty-first century. Present day (PD) simulations show high performance skill in both climate mean and variability, which provide a reliability of the climate models and reduces the uncertainty in response to future forcing. In future projections, global temperature increases from 1.92 °C to 5.20 °C relative to the PD level (1995–2014). Global mean precipitation increases from 5.1% to 10.1% and sea ice extent decreases from 19% to 62% in the Arctic and from 18% to 54% in the Antarctic. In addition, climate changes are accelerating toward the late twenty-first century. Our CMIP6 simulations are released to the public through the Earth System Grid Federation (ESGF) international data sharing portal and are used to support the establishment of the national adaptation plan for climate change in South Korea.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Simulation of Plasma Emission in Magnetized Plasmas

The recent Parker Solar Probe observations of type III radio bursts show that the effects of the finite background magnetic field can be an important factor in the interpretation of data. In the present paper, the effects of the background magnetic field on the plasma-emission process, which is believed to be the main emission mechanism for solar coronal and interplanetary type III radio bursts, are investigated by means of the particle-in-cell simulation method. The effects of the ambient magnetic field are systematically surveyed by varying the ratio of plasma frequency to electron gyrofrequency. The present study shows that for a sufficiently strong ambient magnetic field, the wave–particle interaction processes lead to a highly field-aligned longitudinal mode excitation and anisotropic electron velocity distribution function, accompanied by a significantly enhanced plasma emission at the second-harmonic plasma frequency. For such a case, the polarization of the harmonic emission is almost entirely in the sense of extraordinary mode. On the other hand, for moderate strengths of the ambient magnetic field, the interpretation of the simulation result is less clear. The underlying nonlinear-mode coupling processes indicate that to properly understand and interpret the simulation results requires sophisticated analyses involving interactions among magnetized plasma normal modes, including the two transverse modes of the magneto-active plasma, namely, the extraordinary and ordinary modes, as well as electron-cyclotron-whistler, plasma oscillation, and upper-hybrid modes. At present, a nonlinear theory suitable for quantitatively analyzing such complex-mode coupling processes in magnetized plasmas is incomplete, which calls for further theoretical research, but the present simulation results could provide a guide for future theoretical efforts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Evidence of nonlinear coupling in the edge harmonic oscillation sustaining quiescent high confinement in a tokamak plasma

In a tokamak plasma, we measure nonlinear coupling between harmonics comprising a saturated edge oscillation, as predicted by nonlinear magnetohydrodynamic simulations. The coherent structure formed by this coupling suppresses bursty edge-localized modes that cause energy and particle loss. We also measure nonlinear coupling of the edge harmonic oscillation to core-resonant magnetic tearing modes. When tearing modes are present, the edge oscillation harmonics decouple, and edge-localized modes return. In this article, we conclude that nonlinear interaction with tearing modes interrupts the formation of the edge harmonic oscillation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗