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 685 records · Page 38

Narrowing of the Upwelling Branch of the Brewer-Dobson Circulation and Hadley Cell in Chemistry-Climate Model Simulations of the 21st Century

Changes in the width of the upwelling branch of the Brewer-Dobson circulation and Hadley cell in the 21st Century are investigated using simulations from a coupled chemistry-climate model. In these model simulations the tropical upwelling region narrows in the troposphere and lower stratosphere. The narrowing of the Brewer-Dobson circulation is caused by an equatorward shift of Rossby wave critical latitudes and Eliassen-Palm flux convergence in the subtropical lower stratosphere. In the troposphere, the model projects an expansion of the Hadley cell's poleward boundary, but a narrowing of the Hadley cell's rising branch. Model results suggest that eddy forcing may also play a part in the narrowing of the rising branch of the Hadley cell.

Li, Feng↗

Numerical Investigation of Combustor-Turbine Interactions with Cooling Air Flows Included in High-Pressure Turbine

In this paper, we would like to include the cooling airflows inside the turbine and perform the fully-coupled combustor–turbine simulation at a realistic sea-level takeoff (SLTO) condition. The new results are compared with the previous study,in which we performed the fully-coupled combustor–turbine simulation without the cooling airflows and the purge flows inside the HPT. This paper also serves to provide a better understanding about the effect of the unsteadiness in the turbine inflow condition on the turbine performance in the presence of the cooling airflows.

LES↗

Clouds in GEOS-5

The GEOS-5 atmospheric model is being developed as a weather-and-climate capable model. It must perform well in assimilation mode as well as in weather and climate simulations and forecasts and in coupled chemistry-climate simulations. In developing GEOS-5, attention has focused on the representation of moist processes. The moist physics package uses a single phase prognostic condensate and a prognostic cloud fraction. Two separate cloud types are distinguished by their source: "anvil" cloud originates in detraining convection, and large-scale cloud originates in a PDF-based condensation calculation. Ice and liquid phases for each cloud type are considered. Once created, condensate and fraction from the anvil and statistical cloud types experience the same loss processes: evaporation of condensate and fraction, auto-conversion of liquid or mixed phase condensate, sedimentation of frozen condensate, and accretion of condensate by falling precipitation. The convective parameterization scheme is the Relaxed Arakawa-Schubert, or RAS, scheme. Satellite data are used to evaluate the performance of the moist physics packages and help in their tuning. In addition, analysis of and comparisons to cloud-resolving models such as the Goddard Cumulus Ensemble model are used to help improve the PDFs used in the moist physics. The presentation will show some of our evaluations including precipitation diagnostics.

Bacmeister, Julio↗

Assessing Fan Flutter Stability in the Presence of Inlet Distortion Using One-way and Two-way Coupled Methods

Concerns regarding noise, propulsive efficiency, and fuel burn are inspiring aircraft designs wherein the propulsive turbomachines are partially (or fully) embedded within the airframe; such designs present serious concerns with regard to aerodynamic and aeromechanic performance of the compression system in response to inlet distortion. Previously, a preliminary design of a forward-swept high-speed fan exhibited flutter concerns in cleaninlet flows, and the present author then studied this fan further in the presence of off-design distorted in-flows. Continuing this research, a three-dimensional, unsteady, Navier-Stokes computational fluid dynamics code is again applied to analyze and corroborate fan performance with clean inlet flow and now with a simplified, sinusoidal distortion of total pressure at the aerodynamic interface plane. This code, already validated in its application to assess aerodynamic damping of vibrating blades at various flow conditions using a one-way coupled energy-exchange approach, is modified to include a two-way coupled time-marching aeroelastic simulation capability. The two coupling methods are compared in their evaluation of flutter stability in the presence of distorted in-flows.

aeroelasticity↗

Application of OpenFOAM to Plume Impingement in Space Environments

After 30 years of continuous human presence in low-earth orbit, NASA is returning to the moon and eventually will go to Mars. Travelling beyond low earth orbit requires NASA to learn how humans can live in Deep Space environments – beyond the protection of Earth’s magnetosphere and at distances from Earth that prevent a quick return in case of trouble. To this end, NASA is constructing the Lunar Gateway, an ISS-like space station to be put in orbit around the moon to act as a home base for Lunar exploration for NASA astronauts. The Gateway Lunar outpost will be built incrementally, via modules which will arrive at separate times and dock to the existing structure. The incremental addition of Gateway modules, and the docking of visiting vehicles, is achieved via a sequence of firings from the approaching body’s onboard reaction control system (RCS) thrusters to achieve the required approach trajectory. The typical hypergolic chemical RCS thrusters work by firing hot gases to produce adverse thrust and the needed change in velocity to safely finish the docking process. The exhaust gas from the RCS thrusters form plumes that expand into the vacuum of space and can impinge onto the outer surfaces of the Lunar Gateway, causing unwanted forces and moments, heat loads, sediment deposition, and in extreme cases, even surface erosion - all mechanisms that can damage the Lunar Gateway and must be minimized. Both permanent and visiting modules will have this RCS thruster exhaust impingement problem. This research aims to establish existing OpenFOAM solvers as a methodology for improving simulation techniques of rocket exhaust plume impingement in space environments. The flow structure of a plume in a space environment is complex; a plume that originates from a hypergolic chemical RCS thruster and expands into a vacuum will experience several regimes of rarefication. This range includes the continuum flow in the rocket nozzle through the fully rarefied free molecular flow further from the nozzle. The flow physics is different at these two extremes, and as such, the simulation approach for plumes is generally divided into a traditional computational fluid dynamics (CFD) simulation in and near the nozzle which is coupled to a subsequent direct simulation Monte Carlo (DSMC) simulation. At this time, the scope of this research is developing, verifying, and validating a method using existing solvers in the OpenFOAM framework for performing coupled CFD/DSMC calculations to determine the extent of plume impingement loading on generic space structures. This presentation will detail code-to-code comparisons between the hyStrath dsmcFoam+ solver, developed using OpenFOAM and available as open-source, and NASA’s in-house DSMC Analysis Code (DAC). Comparisons to several open-source publication findings using DAC [3,4] are presented, and advantages of using an OpenFOAM based solver are also discussed. The presentation concludes with a discussion of future work, and a plan for coupling the dsmcFoam+ solver with CFD simulations of chemical rocket engines for unified coupled plume simulation.

DSMC↗

Forecasting Dynamic Line Rating with Spatial Variation Considerations

Dynamic line rating (DLR) is a technology that allows the ampacity of an electrical conductor to be calculated using real-time or forecasted weather conditions. Historically, the ampacity of a conductor has been determined using a static line rating method which assumes conservative weather assumptions. Therefore, not only can DLR give a more accurate measurement of the true ampacity of a conductor, but it can also increase its ampacity during weather conditions with greater thermal mitigations. The two primary cooling factors in the ampacity calculations are wind speed and direction. In complex terrain, wind speed and direction can have large variations over short distances. Therefore, accurately identifying the limiting span of a transmission line requires high spatial resolution of the wind along its path. One solution is to install dense weather stations along their path, though this can become costly over long distances. Therefore, researchers have investigated the use of Computation Fluid Dynamic (CFD) simulations to accurately compute the wind field along the path of a transmission line and use these results to identify the limiting section of the conductor. This work presents a case study that evaluates the coupling of CFD simulations and forecasted weather simulations using the High-Resolution Rapid Refresh (HRRR) model points over a 2-year span within a region in south eastern Idaho. The primary goal of the work is the evaluation of the number of HRRR model points used, i.e., weather stations, along the path of the line and the accuracy of the resulting ampacity. This was done using 4, 10, 17, 26, and 35 HRRR model points along two transmission line paths. The results indicate that as the number of model points are increased, the DLR ampacity of the lines decrease, yet converge as more points are added and demonstrate little change with additional HRRR points. It is expected that these results can help transmission line operators identify the number of weather stations that must be installed when coupled with CFD simulations and DLR ampacity to ensure accurate ratings and safe operations.

17 WIND ENERGY↗

Forecasting Dynamic Line Rating with Spatial Variation Considerations

Dynamic line rating (DLR) is a technology that allows the ampacity of an electrical conductor to be calculated using real-time or forecasted weather conditions. Historically, the ampacity of a conductor has been determined using a static line rating method which assumes conservative weather assumptions. Therefore, not only can DLR give a more accurate measurement of the true ampacity of a conductor, but it can also increase its ampacity during weather conditions with greater thermal mitigations. The two primary cooling factors in the ampacity calculations are wind speed and direction. In complex terrain, wind speed and direction can have large variations over short distances. Therefore, accurately identifying the limiting span of a transmission line requires high spatial resolution of the wind along its path. One solution is to install dense weather stations along their path, though this can become costly over long distances. Therefore, researchers have investigated the use of Computation Fluid Dynamic (CFD) simulations to accurately compute the wind field along the path of a transmission line and use these results to identify the limiting section of the conductor. This work presents a case study that evaluates the coupling of CFD simulations and forecasted weather simulations using the High-Resolution Rapid Refresh (HRRR) model points over a 2-year span within a region in south eastern Idaho. The primary goal of the work is the evaluation of the number of HRRR model points used, i.e., weather stations, along the path of the line and the accuracy of the resulting DLR ampacity. This was done using 4, 10, 17, 26, and 35 HRRR model points along two transmission line paths. The results indicate that as the number of model points are increased, the DLR ampacity of the lines decrease, yet converge as more points are added and demonstrate little change with additional HRRR points. It is expected that these results can help transmission line operators identify the number of weather stations that must be installed when coupled with CFD simulations and DLR ampacity to ensure accurate ratings and safe operations.

17 WIND ENERGY↗

Chirality transfer and chiral turbulence in gauge theories

Chirality transfer between fermions and gauge fields plays a crucial role for understanding the dynamics of anomalous transport phenomena such as the Chiral Magnetic Effect. In this proceeding we present a first principles study of these processes based on classical-statistical real-time lattice simulations of strongly coupled QED (e 2 N f = 64). Furthermore, our simulations demonstrate that a chirality imbalance in the fermion sector triggers chiral plasma instabilities in the gaugefield sector, which ultimately lead to the generation of long range helical magnetic fields via an self-similar turbulent cascade of the magnetic helicity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Spin-lattice dynamics of surface vs core magnetization in Fe nanoparticles

Magnetization of clusters is often simulated using atomistic spin dynamics for a fixed lattice. Coupled spinlattice dynamics simulations of the magnetization of nanoparticles have up to now neglected the change in the size of the atomic magnetic moments near surfaces. We show that the introduction of variable magnetic moments leads to a better description of experimental data for the magnetization of small Fe nanoparticles. To this end we divide atoms into a surface-near shell and a core with bulk properties. Here, it is demonstrated that both the magnitude of the shell magnetic moment and the exchange interactions need to be modified to obtain a fair representation of the experimental data. This allows for a reasonable description of the average magnetic moment versus cluster size, and also the cluster magnetization versus temperature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Optimization of simulated high-field side lower hybrid current drive coupling using machine learning predictions of scrape-off layer density

Lower hybrid current drive (LHCD) is a potential source of non-inductive off-axis current drive (CD) for tokamaks. Although LHCD has been successfully deployed on a number of tokamaks, it is highly sensitive to the scrape-off layer (SOL) conditions local to the LHCD launcher. Large gaps between the launcher and plasma core, SOL turbulence, or edge density perturbations due to edge-localized modes can hamper CD or cause large reflected power. These coupling issues in part motivated the installation of an LHCD launcher on the high-field side (HFS) of DIII-D. On the HFS, the SOL is less turbulent and more controllable compared to the low-field side. This quiescence may result in more predictable edge conditions and thus a more predictable CD. Here, in this work, HFS SOL reflectometry measurements are predicted from global plasma parameters using machine learning models. The SOL predictions coupled with the full-wave simulation of the LHCD launcher allow for the prediction of reflected power, directivity, and arcing risk before the discharge. Launcher performance is then optimized using multi-objective Bayesian optimization, finding the shot parameters that result in an optimal SOL density that maximizes CD while minimizing the risk of arcing. The predictions and optimizations of LHCD performance are then accelerated using a surrogate model of the full-wave LHCD simulation.

Bayesian optimization↗

Simulations of Radiation Hydrodynamics in Strongly Coupled Plasmas (Final Report)

Radiation propagates throughout the Universe and is often the only information we have from distant astronomical objects. Additionally, radiation plays a key role in systems ranging from the interior of stars to the establishment of shocks at the edges of forming galaxies to neutron star mergers. Here radiative processes, absorption, emission, and scattering, are fundamental to transport of momentum and energy, as well as, the overall evolution of these and other astrophysical systems. In some cases, radiation-dominant systems are also sufficiently dense and ionized, such that, the ions are strongly coupled (ie electrostatic energy >> thermal energy or Γ >> 1), which will affect radiative processes and the transport of radiation. During this project, we developed a computational ray-tracing algorithm to study radiation transport through a simulated neutron star capsule implosion under conditions relevant to Omega-60 laser experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

Ice Sheet and Ice Shelf Freshwater Forcing of the NASA GISS ModelE GCM

Recent mass loss from ice sheets and ice shelves is now persistent and prolonged enough that it impacts downstream oceanographic conditions. To demonstrate this, we use an ensemble of coupled GISS-E2.1-G (a CMIP6 model) simulations forced with historical estimates of anomalous freshwater, in addition to other climate forcings, from 1990 through 2019. In this ensemble there are detectable differences in zonal-mean sea surface temperatures (SST) and sea ice in the Southern Ocean, and in regional sea level around Antarctica and in the western North Atlantic. These impacts mostly improve the model's representation of historical changes, including reversing the forced trends in Southern Ocean surface temperature and Antarctic sea ice. The changes in SST may have implications for estimates of the SST pattern effect on climate sensitivity and for cloud feedbacks. We conclude that the changes are sufficiently large that model groups should strive to include more accurate estimates of these drivers in all-forcing historical simulations in future coupled model intercomparisons.

Ice sheet↗

Analysis of Inlet-Compressor Acoustic Interactions Using Coupled CFD Codes

A problem that arises in the numerical simulation of supersonic inlets is the lack of a suitable boundary condition at the engine face. In this paper, a coupled approach, in which the inlet computation is coupled dynamically to a turbomachinery computation, is proposed as a means to overcome this problem. The specific application chosen for validation of this approach is the collapsing bump experiment performed at the University of Cincinnati. The computed results are found to be in reasonable agreement with experimental results. The coupled simulation results could also be used to aid development of a simplified boundary condition.

Suresh, A.↗

Simulation of runaway electron production with CQL3D coupled to NIMROD

Abstract A coupling between two distinctly different codes—one magnetohydrodynamic (MHD) and another kinetic—is achieved and applied for simulation of runaway electron (RE) production. The 3D initial value MHD code NIMROD simulates a DIII-D pure neon shattered pellet injection plasma quench including the propagation and ablation of the fragments, ionization and recombination of the impurities, and the radiated and transported energies. The field data from NIMROD is then used by the bounce-averaged Fokker–Planck Collisional QuasiLinear 3D (CQL3D) kinetic code to simulate the production of REs and their radial transport. The coupling procedure involves mapping of data between different grids and adjustment of the NIMROD toroidal electric field when REs appear. It is shown that without the radial transport, a large RE current is generated, up to 30% of the pre-pellet ohmic current. However, when the radial transport is included in CQL3D, the RE current is reduced to undetectable level, consistent with experiment. Various forms of the radial diffusion are surveyed to determine conditions when the fast electrons would not have time to be accelerated to relativistic energies before they are lost to chamber wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Coupled radiation transport and plasma physics for ICF simulations [Slides]

Describing high-temperature plasmas requires inclusion of several nonlinear coupled physics, including ion kinetics - for each isotope species present in the plasma, electron kinetics, thermal radiative transfer (TRT), electromagnetic description, and nuclear reactions (fusion burn). Simulating plasmas in ICF hohlraums aims to solve kinetic ions (particles), fluid electrons, energetic electrons, electromagnetics, kinetic radiation (MG DP solver), and laser-plasma interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sensitivity Analysis of Multidisciplinary Rotorcraft Simulations

A multidisciplinary sensitivity analysis of rotorcraft simulations involving tightly coupled high-fidelity computational fluid dynamics and comprehensive analysis solvers is presented and evaluated. An unstructured sensitivity-enabled Navier-Stokes solver, FUN3D, and a nonlinear flexible multibody dynamics solver, DYMORE, are coupled to predict the aerodynamic loads and structural responses of helicopter rotor blades. A discretely-consistent adjoint-based sensitivity analysis available in FUN3D provides sensitivities arising from unsteady turbulent flows and unstructured dynamic overset meshes, while a complex-variable approach is used to compute DYMORE structural sensitivities with respect to aerodynamic loads. The multidisciplinary sensitivity analysis is conducted through integrating the sensitivity components from each discipline of the coupled system. Numerical results verify accuracy of the FUN3D/DYMORE system by conducting simulations for a benchmark rotorcraft test model and comparing solutions with established analyses and experimental data. Complex-variable implementation of sensitivity analysis of DYMORE and the coupled FUN3D/DYMORE system is verified by comparing with real-valued analysis and sensitivities. Correctness of adjoint formulations for FUN3D/DYMORE interfaces is verified by comparing adjoint-based and complex-variable sensitivities. Finally, sensitivities of the lift and drag functions obtained by complex-variable FUN3D/DYMORE simulations are compared with sensitivities computed by the multidisciplinary sensitivity analysis, which couples adjoint-based flow and grid sensitivities of FUN3D and FUN3D/DYMORE interfaces with complex-variable sensitivities of DYMORE structural responses.

Wang, Li↗

A coupled fluid-mechanical workflow to simulate the directed energy deposition additive manufacturing process

Simulation of additive manufacturing processes can provide essential insight into material behavior, residual stress, and ultimately, the performance of additively manufactured parts. In this work, we describe a new simulation based workflow utilizing both solid mechanics and fluid mechanics based formulations within the finite element software package SIERRA (Sierra Solid Mechanics Team in Sierra/Solid Mechanics 4.52 User’s Guide SAND2019-2715. Technical report, Sandia National Laboratories, 2011) to enable integrated simulations of directed energy deposition (DED) additive manufacturing processes. In this methodology, a high-fidelity fluid mechanics based model of additive manufacturing is employed as the first step in a simulation workflow. This fluid model uses a level set field to track the location of the boundary between the solid material and background gas and precisely predicts temperatures and material deposition shapes from additive manufacturing process parameters. Next, the resulting deposition shape and temperature field from the fluid model are then mapped into a solid mechanics formulation to provide a more accurate surface topology for radiation and convection boundary conditions and a prescribed temperature field. Solid mechanics simulations are then conducted to predict the evolution of material stresses and microstructure within a part. By combining thermal history and deposition shape from fluid mechanics with residual stress and material property evolutions from solid mechanics, additional fidelity and precision are incorporated into additive manufacturing process simulations providing new insight into complex DED builds.

36 MATERIALS SCIENCE↗

Analysis of Coupled Seals, Secondary and Powerstream Flow Fields in Aircraft and Aerospace Turbomachines

Higher power, high efficiency gas turbine engines require optimization of the seals and secondary flow systems as well as their impact on the powerstream. This work focuses on two aspects: 1. To apply the present day CFD tools (SCISEAL) to different real-life secondary flow applications from different original equipment manufacturers (OEM s) to provide feedback data and 2. Develop a computational methodology for coupled time-accurate simulation of the powerstream and secondary flow with emphasis on the interaction between the disk-cavity and rim seals flows with the powerstream (SCISEAL-MS-TURBO). One OEM simulation was of the Allison Engine Company T-56 turbine drum cavities including conjugate heat transfer with good agreement with data and provided design feedback information. Another was the GE aspirating seal where the 3-D CFD simulations played a major role in analysis and modification of that seal configuration. The second major objective, development of a coupled flow simulation capability was achieved by using two codes MS-TURBO for the powerstream and SCISEAL for the secondary flows with an interface coupling algorithm. The coupled code was tested against data from three differed configurations: 1. bladeless-rotor-stator-cavity turbine test rig, 2. UTRC high pressure turbine test rig, and, 3. the NASA Low-Speed-Air Compressor rig (LSAC) with results and limitations discussed herein.

Athavale, M. M.↗