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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Assimilation of AERONET and MODIS AOT Observations Using Variational and Ensemble Data Assimilation Methods and Its Impact on Aerosol Forecasting Skill

Data assimilation of Aerosol Robotic Network (AERONET) and Moderate Resolution Imaging Spectroradiometer (MODIS) aerosol optical thickness (AOT) for aerosol forecasting was tested within the Navy Aerosol Analysis Prediction System (NAAPS) framework, using variational and ensemble data assimilation methods. Navy aerosol forecasting currently makes use of a deterministic NAAPS simulation coupled to Navy Variational Data Assimilation System for aerosol optical depth, a two-dimensional variational data assimilation system, for MODIS AOT assimilation. An ensemble version of NAAPS (ENAAPS) coupled to an ensemble adjustment Kalman filter (EAKF) from the Data Assimilation Research Testbed was recently developed, allowing for a range of data assimilation and forecasting experiments to be run with deterministic NAAPS and ENAAPS. The main findings are that the EAKF, with its flow-dependent error covariances, makes better use of sparse observations such as AERONET AOT. Assimilating individual AERONET observations in the two-dimensional variational system can increase the analysis errors when observations are located in high AOT gradient regions. By including AERONET with MODIS AOT assimilation, the magnitudes of peak aerosol events (AOT> 1) were better captured with improved temporal variability, especially in India and Asia where aerosol prediction is a challenge. Assimilating AERONET AOT with MODIS had little impact on the 24 h forecast skill compared to MODIS assimilation only, but differences were found downwind of AERONET sites. The 24 h forecast skill was approximately the same for forecasts initialized with analyses from AERONET AOT assimilation alone compared to MODIS assimilation, particularly in regions where the AERONET network is dense; including the United States and Europe, indicating that AERONET could serve as a backup observation network for over-land synoptic-scale aerosol events.

Rubin, Juli I.↗
NASA NTRSDate not supplied

RDE Nozzle Computational Design Methodology Development and Application

The rotating detonation engine (RDE) is currently under investigation as an approach to achieving pressure gain combustion for propulsion and power systems. The RDE typically consists of an annulus with one end open (or having a throat and/or nozzle) and the other end valved using non-mechanical, fluidic means. In this study, a total of 6 nozzle design variations were studied for the NPS combustor geometry, using the validated quasi-2D to 3-D coupled simulation approach.

CFD↗
NASA NTRSDate not supplied

Comparison between Hayabusa 2 Spectral Measurements and Simulations

This paper compares recent Hayabusa 2 spectral measurements with state-of-the-art shock-layer radiation simulations resulting from the LAURA/HARA code suite. These simulations include coupled ablation, which accounts for the injection of ablation products into the flowfield, and coupled radiation, which accounts for radiative energy loss in the flowfield. To enable the coupled ablation simulations, a best-estimate model is developed for Hayabusa’s carbon-phenolic ablator, based on the limited available published information. The comparison between the simulations and measurements focuses on two atomic nitrogen lines and two atomic oxygen lines, as well as the CN Violet band system. For the atomic lines, the measurements and simulations agree within 25% over most of the trajectory. This excellent agreement is unprecedented for observed radiation measurements of atomic lines, where previous Stardust and Hayabusa 1 comparisons were significantly worse. The improved agreement for these Hayabusa 2 comparisons is both the result of improved measurement quality and enhanced flowfield/radiation modeling. To provide a link between these observed radiation measurements and the radiative heating to a vehicle surface, the recently developed flowfield property binning approach is used to identify flowfield properties that provide the dominant emission contribution to the measured spectrum. These identified flowfield properties are then shown to match those for the surface radiative heating to the currently developed Mars Sample Return (MSR) Earth Entry System (EES) at a specific trajectory point and surface location. This matching of flowfield properties indicates that this EES radiative heating and Hayabusa 2 observed radiation simulations are equivalent radiation problems. Therefore, the excellent agreement between Hayabusa 2 measurements and simulations may be leveraged to inform the radiation heating margin for the MSR EES.

Christopher O Johnston↗
NASA NTRS2023

Convective Impact on the Global Lower Stratospheric Water Vapor Budget

Water vapor in the stratosphere is primarily controlled by temperatures in the tropical upper troposphere and lower stratosphere. However, the direct impact of deep convection on the global lower stratospheric water vapor budget is still an actively debated issue. Two complementary modeling approaches are used to investigate the convective impact in boreal winter and summer. Convective influence is diagnosed by tracing trajectories through convective cloud top altitude fields derived from global rainfall and brightness temperature data. Backward trajectory model simulations coupled with a detailed treatment of cloud microphysical processes indicate that convection moistens the global lower stratosphere by approximately 0.3 ppmv in boreal winter and summer 2010. The diurnal peak in convection is responsible for about half of the total convective moistening during winter and nearly all of the convective moistening during summer. Deep convective clouds overshooting the tropopause have relatively minor effect on global lower stratospheric water vapor. A forward trajectory model coupled with a simplified cloud module is used to estimate the relative magnitude of the interannual variability of the convective impact. Combining the results from the two models, we find that the convective impact on the global lower stratospheric water vapor during 2006-2016 is approximately 0.3 ppmv with year-to-year variations of up to 0.1 ppmv. An important mechanism of convective hydration of the lower stratosphere is via the detrainment of saturated air and ice into the tropical uppermost troposphere and the subsequent upward transport of some of these moist air parcels across relatively warm and subsaturated tropopause.

CONVECTION↗
NASA NTRSDate not supplied

Analysis of Mars 2020 Entry with Coupled Material Response and CFD

When computing the response of an ablating thermal protection system during entry, the aerothermal environment and material response are traditionally computed independently, with the introduction of a blowing correction term in the material response model to account for the outgassing of char and pyrolysis gases. This study presents an approach, where the pyrolysis blowing gases, calculated within the PATO material response code [1], are integrated into the DPLR hypersonic CFD code [2] via a blowing boundary condition. This methodology employs an iterative process, whereby the blown pyrolysis gas products from PATO are incorporated into DPLR, refining surface heating predictions. The NEQAIR code is used to compute radiative heating [3]. The method is applied in computing 3D material response of the Mars 2020 entry. When compared with uncoupled material response, the coupled simulations show a lower surface heat flux initially and a higher heat flux at peak heating as shown in Fig. 1. This is in agreement with previous work using a sphere case with environments from MSL [4].

Heat Transfer↗
NASA NTRSDate not supplied

Analysis of Mars 2020 Entry with Coupled Material Response and CFD

When computing the response of an ablating thermal protection system during entry, the aerothermal environment and material response are traditionally computed independently, with the introduction of a blowing correction term in the material response model to account for the outgassing of char and pyrolysis gases. This study presents an approach, where the pyrolysis blowing gases, calculated within the PATO material response code, are integrated into the DPLR hypersonic CFD code via a blowing boundary condition. This methodology employs an iterative process, whereby the blown pyrolysis gas products from PATO are incorporated into DPLR, refining surface heating predictions. The NEQAIR code is used to compute radiative heating. The method is applied in computing 3D material response of the Mars 2020 entry. When compared with uncoupled material response, the coupled simulations show a lower surface heat flux initially and a higher heat flux at peak heating as shown in Fig. 1. This is in agreement with previous work using a sphere case with environments from MSL.

Heat Transfer↗
NASA NTRSDate not supplied

High-Fidelity Simulations of Lift+Cruise VTOL Urban Air Mobility Concept Aircraft in Hover

This paper presents progress of high-fidelity multidisciplinary simulations for the NASA lift+cruise vertical takeoff and landing (VTOL) urban air mobility concept aircraft in hover. The study focuses on the unsteady flow solutions featuring strong wake interactions between rotors and the airframe. The simulations couple a high-fidelity aerodynamic model with a comprehensive rotorcraft aeromechanics tool with rotor rotation-speed trim. The aerodynamic model is based on the Reynolds-averaged Navier-Stokes equations using the one-equation Spalart-Allmaras turbulence model with rotation correction. Unsteady aerodynamic flows are computed on a dynamic, deformable, unstructured, overset grid system. An integrated overset-grid assembler is used to construct the composite grid from 23 component grids and facilitate communications between individual component grids. Rotor performance and airframe forces and moments are computed and compared on two unstructured grids. The study demonstrates capabilities of high-fidelity multidisciplinary analysis tools to capture the strong unsteady flowfield around this multi-rotor aircraft in hover.

CFD↗
NASA NTRSDate not supplied

High-Fidelity Simulations of Lift+Cruise VTOL Urban Air Mobility Concept Aircraft in Hover

This presentation presents progress of high-fidelity multidisciplinary simulations for the NASA lift+cruise vertical takeoff and landing (VTOL) urban air mobility concept aircraft in hover. The study focuses on the unsteady flow solutions featuring strong wake interactions between rotors and the airframe. The simulations couple a high-fidelity aerodynamic model with a comprehensive rotorcraft aeromechanics tool with rotor rotation-speed trim. The aerodynamic model is based on the Reynolds-averaged Navier-Stokes equations using the one-equation Spalart-Allmaras turbulence model with rotation correction. Unsteady aerodynamic flows are computed on a dynamic, deformable, unstructured, overset grid system. An integrated overset-grid assembler is used to construct the composite grid from 23 component grids and facilitate communications between individual component grids. Rotor performance and airframe forces and moments are computed and compared on two unstructured grids. The study demonstrates capabilities of high-fidelity multidisciplinary analysis tools to capture the strong unsteady flowfield around this multi-rotor aircraft in hover.

CFD↗
DOE OSTI2020

Direct observation of 3D atomic packing in monatomic amorphous materials

Liquids and solids are two fundamental states of matter. Although the structure of crystalline solids has long been solved by crystallography, our understanding of the 3D atomic structure of liquids and amorphous materials remained speculative. In 1952, Frank hypothesized that icosahedral order is the prevalent atomic motif in monatomic liquids. Over the past six decades, there have been a great deal of experimental, computational and theoretical studies to understand the structure of liquids and amorphous materials. A polytetrahedral packing model was proposed to explain the 3D atomic structure of monatomic liquids and amorphous materials, in which icosahedral order is a key feature. The icosahedral order has also been found to play an important role in the structure of metallic glasses and quasicrystals. Here we advance atomic electron tomography to determine for the first time the 3D atomic positions in monatomic amorphous materials, including a Ta thin film and two Pd nanoparticles. Despite the different chemical composition and synthesis methods, we observe that pentagonal bipyramids are the most abundant atomic motifs in these amorphous materials. Contrary to the traditional understanding, the pentagonal bipyramids do not assemble full icosahedra, but are closely connected with each other to form networks that extend to medium-range scale. Our molecular dynamic simulations, coupled with experimental results, further reveal that the 3D atomic structure of monatomic liquids is similar to the experimental amorphous material. Moreover, we find that the pentagonal bipyramid networks, prevalent in monatomic liquids, rapidly grow in size and form a fraction of full icosahedra during the quench from a liquid to a metallic glass state, providing a possible new picture of structural evolution through the glass transition.

Yuan, Yakun↗
DOE OSTI2021

xSDK: Building an ecosystem of highly efficient math libraries for exascale

Current efforts to build increasingly powerful computer architectures are opening up new avenues for more complex and higher fidelity simulations coupled with data analytics and learning, leading to new scientific insights and deeper understanding. At one extreme, exascale computers will be much faster than previous computer generations (performing 10 18 operations per second—that is, 1,000 times faster than petascale). To achieve these performance improvements, computer architectures are becoming increasingly complex, with deep memory hierarchies, very high node and core counts, and heterogeneous features such as graphics processing units (GPUs). Such architectural changes impact the full breadth of computing scales, as heterogeneity pervades even current-generation laptops, workstations, and moderate-sized clusters. While emerging advanced architectures provide unprecedented opportunities, they also present significant challenges for developers of scientific applications, such as multiphysics and multiscale codes, who must adapt their software to handle disruptive changes in architectures and new programming models that have not yet stabilized. Developers must consider increasing concurrency while reducing communication and synchronization, and other complexities such as the potential for using mixed precision to leverage the compute power available in low-precision tensor cores. On one hand, developers must implement new scientific capabilities, which in turn increase code complexity. On the other hand, the codes must be ported to new architectures, requiring the inclusion of new programming models and the restructuring of code to achieve good performance. Addressing these issues is beyond the capability of any single person or team—leading to the need for collaboration among many teams, who encapsulate their expertise in reusable software and work together to create sustainable software ecosystems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Spatial coupling of gyrokinetic simulations, a generalized scheme based on first-principles

In this paper, we present a scheme that spatially couples two gyrokinetic codes using first-principles. Coupled equations are presented and a necessary and sufficient condition for ensuring accuracy is derived. This new scheme couples both the field and the particle distribution function. The coupling of the distribution function is only performed once every few time-steps, using a five-dimensional (5D) grid to communicate the distribution function between the two codes. This 5D grid interface enables the coupling of different types of codes and models, such as particle and continuum codes, or delta-f and total-f models. Transferring information from the 5D grid to the marker particle weights is achieved using a new resampling technique. Demonstration of the coupling scheme is shown using two XGC gyrokinetic simulations for both the core and edge. We also apply the coupling scheme to two continuum simulations for a one-dimensional advection–diffusion problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Multimode two-dimensional vibronic spectroscopy. II. Simulating and extracting vibronic coupling parameters from polarization-selective spectra

We report experimental demonstrations of polarization-selection two-dimensional Vibrational-Electronic (2D VE) and 2D Electronic-Vibrational (2D EV) spectroscopies aim to map the magnitudes and spatial orientations of coupled electronic and vibrational coordinates in complex systems. The realization of that goal depends on our ability to connect spectroscopic observables with molecular structural parameters. In this paper, we use a model Hamiltonian consisting of two anharmonically coupled vibrational modes in electronic ground and excited states with linear and bilinear vibronic coupling terms to simulate polarization-selective 2D EV and 2D VE spectra. We discuss the relationships between the linear vibronic coupling and two-dimensional Huang–Rhys parameters and between the bilinear vibronic coupling term and Duschinsky mixing. We develop a description of the vibronic transition dipoles and explore how the Hamiltonian parameters and non-Condon effects impact their amplitudes and orientations. Using simulated polarization-selective 2D EV and 2D VE spectra, we show how 2D peak positions, amplitudes, and anisotropy can be used to measure parameters of the vibronic Hamiltonian and non-Condon effects. This paper, along with the first in the series, provides the reader with a detailed description of reading, simulating, and analyzing multimode, polarization-selective 2D EV and 2D VE spectra with an emphasis on extracting vibronic coupling parameters from complex spectra.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

A 1D–0D–3D coupled model for simulating blood flow and transport processes in breast tissue

In this work, we present mixed dimensional models for simulating blood flow and transport processes in breast tissue and the vascular tree supplying it. These processes are considered, to start from the aortic inlet to the capillaries and tissue of the breast. Large variations in biophysical properties and flow conditions exist in this system necessitating the use of different flow models for different geometries and flow regimes. In total, we consider four different model types. First, a system of 1D nonlinear hyperbolic partial differential equations (PDEs) is considered to simulate blood flow in larger arteries with highly elastic vessel walls. Second, we assign 1D linearized hyperbolic PDEs to model the smaller arteries with stiffer vessel walls. The third model type consists of ODE systems (0D models). It is used to model the arterioles and peripheral circulation. Finally, homogenized 3D porous media models are considered to simulate flow and transport in capillaries and tissue within the breast volume. Sink terms are used to account for the influence of the venous and lymphatic systems. Combining the four model types, we obtain two different 1D–0D–3D coupled models for simulating blood flow and transport processes: The first model results in a fully coupled 1D–0D–3D model covering the complete path from the aorta to the breast combining a generic arterial network with a patient specific breast network and geometry. The second model is a reduced one based on the separation of the generic and patient specific parts. The information from a calibrated fully coupled model is used as inflow condition for the patient specific sub-model allowing a significant computational cost reduction. Several numerical experiments are conducted to calibrate the generic model parameters and to demonstrate realistic flow simulations compared to existing data on blood flow in the human breast and vascular system. Moreover, we use two different breast vasculature and tissue data sets to illustrate the robustness of our reduced sub-model approach.

multiscale blood flow models↗
NASA NTRS2021

Large-scale features of Last Interglacial climate: results from evaluating the lig127k simulations for the Coupled Model Intercomparison Project (CMIP6)–Paleoclimate Modeling Intercomparison Project (PMIP4)

The modeling of paleoclimate, using physically based tools, is increasingly seen as a strong out-of-sample test of the models that are used for the projection of future climate changes. New to the Coupled Model Intercomparison Project (CMIP6) is the Tier 1 Last Interglacial experiment for 127 000 years ago (lig127k), designed to address the climate responses to stronger orbital forcing than the midHolocene experiment, using the same state-of-the-art models as for the future and following a common experimental protocol. Here we present a first analysis of a multi-model ensemble of 17 climate models, all of which have completed the CMIP6 DECK (Diagnostic, Evaluation and Characterization of Klima) experiments. The equilibrium climate sensitivity (ECS) of these models varies from 1.8 to 5.6°C. The seasonal character of the insolation anomalies results in strong summer warming over the Northern Hemisphere continents in the lig127k ensemble as compared to the CMIP6 piControl and much-reduced minimum sea ice in the Arctic. The multi-model results indicate enhanced summer monsoonal precipitation in the Northern Hemisphere and reductions in the Southern Hemisphere. These responses are greater in the lig127k than the CMIP6 midHolocene simulations as expected from the larger insolation anomalies at 127 than 6 ka. New synthesis for surface temperature and precipitation, targeted for 127 ka, have been developed for comparison to the multi-model ensemble. The lig127k model ensemble and data reconstructions are in good agreement for summer temperature anomalies over Canada, Scandinavia, and the North Atlantic and for precipitation over the Northern Hemisphere continents. The model–data comparisons and mismatches point to further study of the sensitivity of the simulations to uncertainties in the boundary conditions and of the uncertainties and sparse coverage in current proxy reconstructions. The CMIP6–Paleoclimate Modeling Intercomparison Project (PMIP4) lig127k simulations, in combination with the proxy record, improve our confidence in future projections of monsoons, surface temperature, and Arctic sea ice, thus providing a key target for model evaluation and optimization.

Last Interglacial climate↗
NASA NTRS1990

Micro Blowing Simulations Using a Coupled Finite-Volume Lattice-Boltzman n L ES Approach

Three dimensional large-eddy simulations (LES) of single and multiple jet-in-cross-flow (JICF) are conducted using the 19-bit Lattice Boltzmann Equation (LBE) method coupled with a conventional finite-volume (FV) scheme. In this coupled LBE-FV approach, the LBE-LES is employed to simulate the flow inside the jet nozzles while the FV-LES is used to simulate the crossflow. The key application area is the use of this technique is to study the micro blowing technique (MBT) for drag control similar to the recent experiments at NASA/GRC. It is necessary to resolve the flow inside the micro-blowing and suction holes with high resolution without being restricted by the FV time-step restriction. The coupled LBE-FV-LES approach achieves this objectives in a computationally efficient manner. A single jet in crossflow case is used for validation purpose and the results are compared with experimental data and full LBE-LES simulation. Good agreement with data is obtained. Subsequently, MBT over a flat plate with porosity of 25% is simulated using 9 jets in a compressible cross flow at a Mach number of 0.4. It is shown that MBT suppresses the near-wall vortices and reduces the skin friction by up to 50 percent. This is in good agreement with experimental data.

Menon, S.↗
DOE OSTI2021

Consistent coupling algorithms for coupled core-edge simulations of plasma turbulence

Two schemes for coupling gyrokinetic simulations of microturbulence in tokamaks are proposed. The first scheme is based on an additive Schwarz domain decomposition. In this work, we show that, because the goal of turbulence is long-time averages of the dynamics rather than temporal accuracy, the iteration to self-consistency across domains, which is typically required by Schwarz schemes, can be avoided, thereby accelerating the computation. Second, we propose a coupling scheme that relies entirely on the addition of source terms, leaving the boundary conditions arbitrary. The practical motivations for such a scheme are discussed, and forms of the source terms that ensure consistency and stability are derived. The schemes are tested on a nonlinear, one-dimensional model problem, and the first scheme is further tested on the Hasegawa–Wakatani model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2006

Numerical Methodology for Coupled Time-Accurate Simulations of Primary and Secondary Flowpaths in Gas Turbines

Detailed information of the flow-fields in the secondary flowpaths and their interaction with the primary flows in gas turbine engines is necessary for successful designs with optimized secondary flow streams. Present work is focused on the development of a simulation methodology for coupled time-accurate solutions of the two flowpaths. The secondary flowstream is treated using SCISEAL, an unstructured adaptive Cartesian grid code developed for secondary flows and seals, while the mainpath flow is solved using TURBO, a density based code with capability of resolving rotor-stator interaction in multi-stage machines. An interface is being tested that links the two codes at the rim seal to allow data exchange between the two codes for parallel, coupled execution. A description of the coupling methodology and the current status of the interface development is presented. Representative steady-state solutions of the secondary flow in the UTRC HP Rig disc cavity are also presented.

Przekwas, A. J.↗