Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation decomposition”

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 415 records · Page 23

Hicks Pries, Sulman, et al 2018. In situ incubation of 13C-labeled litter in three soil pits at Blodgett Forest, CA from 2013 to 2016

Even though over half of the world's soil organic carbon (SOC) is stored in subsoils (>20 cm deep), and the old ages of subsoil OC indicate its cycling differs from surface SOC, there are few studies examining in situ decomposition processes in deep soils. The purpose of this dataset is to elucidate these processes. In this study we added 13C-labeled fine roots to 15, 55, and 95 cm depths of a well-characterized coniferous forest Alfisol and monitored the amount of root-derived C remaining over 6, 12, and 30 months. We recovered the root-derived C in microbial phospholipid fatty acids (PLFAs) after 6 months and in coarse (>2 mm) particulate, fine (<2 mm) particulate, and dense, mineral-associated pools after 6, 12, and 30 months. Overall, root decomposition in the first 6 months was similar among all depths but significantly diverged at 30 months with faster decomposition at 15 cm than at 95 cm. There were more fungal and Gram negative-associated PLFAs at 15 cm than at 95 cm, and 13C analysis revealed those microbial groups preferred the added root carbon to native SOC. Mineral-associations were not the cause of slower decomposition at depth because similar amounts of applied root C was recovered in the dense fraction at all depths. The largest difference among depths was in the amount of root C recovered in the coarse particulate fraction, which was greater at 95 cm (50%) than at 15 cm (15%). Slower decomposition of the particulate pool at depth likely contributed to the increase in C:N ratios and depletion of δ13C values below 60 cm depth in our soil profiles. Simulations of these soils using the CORPSE model, which incorporates microbial priming effects and mineral stabilization of SOC, reproduced patterns of particulate and mineral-associated SOC over both time and depth and suggested that a lack of priming by root exudates at depth could account for the slower decomposition rate of particulate root material. Decomposition of deep particulate SOC may increase if root exudation or dissolved OC transport to depth increases.This dataset includes (1) characterization of the soil pits prior to litter addition, (2) characterization of how the litter changed after 6, 12, and 30 months, and (3) the CORPSE model structure, parameterization, and output.

54 ENVIRONMENTAL SCIENCES↗

A dual potential formulation of the Navier-Stokes equations

A dual potential formulation for numerically solving the Navier-Stokes equations is developed and presented. The velocity field is decomposed using a scalar and vector potential. Vorticity and dilatation are used as the dependent variables in the momentum equations. Test cases in two dimensions verify the capability to solve flows using approximations from potential flow to full Navier-Stokes simulations. A three-dimensional incompressible flow formulation is also described. An interesting feature of this approach to solving the Navier-Stokes equations is the decomposition of the velocity field into a rotational part (vector potential) and an irrotational part (scalar potential). The Helmholtz decomposition theorem allows this splitting of the velocity field. This approach has had only limited use since it increases the number of dependent variables in the solution. However, it has often been used for incompressible flows where the solution scheme is known to be fast and accurate. This research extends the usage of this method to fully compressible Navier-Stokes simulations by using the dilatation variable along with vorticity. A time-accurate, iterative algorithm is used for the uncoupled solution of the governing equations. Several levels of flow approximation are available within the framework of this method. Potential flow, Euler and full Navier-Stokes solutions are possible using the dual potential formulation. Solution efficiency can be enhanced in a straightforward way. For some flows, the vorticity and/or dilatation may be negligible in certain regions (e.g., far from a viscous boundary in an external flow). It is possible to drop the calculation of these variables then and optimize the solution speed. Also, efficient Poisson solvers are available for the potentials. The relative merits of non-primitive variables versus primitive variables for solution of the Navier-Stokes equations are also discussed.

Gegg, S. G.↗

Efficient Modal Basis Selection Criteria for Reduced-Order Nonlinear Simulation

A modal basis selection technique for a reduced-order nonlinear numerical simulation with application to two-dimensional structures is presented as a two-step procedure. A system identification analysis is first performed using proper orthogonal decomposition. Using these results, a set of load-invariant bases consisting of the normal modes is next selected. Two criteria for making the basis selection are offered; one using the modal assurance criterion and the other using the modal expansion theorem. The quality of the subsequent reduced-order analyses are examined through comparison with computationally intensive finite element nonlinear simulations in physical degrees-of-freedom. A clamped flat isotropic plate under a random acoustic loading is considered to demonstrate the procedure. It is found that the subject procedure enables formation of an accurate and computationally efficient reduced-order system applicable to a broad range of loading conditions.

Przekop, Adam↗

On the Existence of Fast Modes in Compressible Magnetohydrodynamic Turbulence

Abstract We study the existence and properties of fast magnetosonic modes in 3D compressible MHD turbulence by carrying out a number of simulations with compressible and incompressible driving conditions. We use two approaches to determine the presence of fast modes: mode decomposition based on spatial variations only and spatio-temporal 4D fast Fourier transform (4D FFT) analysis of all fluctuations. The latter method enables us to quantify fluctuations that satisfy the dispersion relation of fast modes with finite frequency. Overall, we find that the fraction of fast modes identified via the spatio-temporal 4D FFT approach in total fluctuation power is either tiny with nearly incompressible driving or ∼2% with highly compressible driving. We discuss the implications of our results for understanding the compressible fluctuations in space and astrophysical plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗

Stoichiometrically coupled carbon and nitrogen cycling in the MIcrobial-MIneral Carbon Stabilization model version 1.0 (MIMICS-CN v1.0)

Explicit consideration of microbial physiology in soil biogeochemical models that represent coupled carbon–nitrogen dynamics presents opportunities to deepen understanding of ecosystem responses to environmental change. The MIcrobial-MIneral Carbon Stabilization (MIMICS) model explicitly represents microbial physiology and physicochemical stabilization of soil carbon (C) on regional and global scales. Here we present a new version of MIMICS with coupled C and nitrogen (N) cycling through litter, microbial, and soil organic matter (SOM) pools. The model was parameterized and validated against C and N data from the Long-Term Inter-site Decomposition Experiment Team (LIDET; six litter types, 10 years of observations, and 13 sites across North America). The model simulates C and N losses from litterbags in the LIDET study with reasonable accuracy (C: R 2 =0.63; N: R 2 =0.29), which is comparable with simulations from the DAYCENT model that implicitly represents microbial activity (C: R 2 =0.67; N: R 2 =0.30). Subsequently, we evaluated equilibrium values of stocks (total soil C and N, microbial biomass C and N, inorganic N) and microbial process rates (soil heterotrophic respiration, N mineralization) simulated by MIMICS-CN across the 13 simulated LIDET sites against published observations from other continent-wide datasets. We found that MIMICS-CN produces equilibrium values in line with measured values, showing that the model generates plausible estimates of ecosystem soil biogeochemical dynamics across continental-scale gradients. MIMICS-CN provides a platform for coupling C and N projections in a microbially explicit model, but experiments still need to identify the physiological and stoichiometric characteristics of soil microbes, especially under environmental change scenarios.

58 GEOSCIENCES↗

MPAS-Seaice (v1.0.0): sea-ice dynamics on unstructured Voronoi meshes

Abstract. We present MPAS-Seaice, a sea-ice model which uses the Model for Prediction Across Scales (MPAS) framework and spherical centroidal Voronoi tessellation (SCVT) unstructured meshes. As well as SCVT meshes, MPAS-Seaice can run on the traditional quadrilateral grids used by sea-ice models such as CICE. The MPAS-Seaice velocity solver uses the elastic–viscous–plastic (EVP) rheology and the variational discretization of the internal stress divergence operator used by CICE, but adapted for the polygonal cells of MPAS meshes, or alternatively an integral (“finite-volume”) formulation of the stress divergence operator. An incremental remapping advection scheme is used for mass and tracer transport. We validate these formulations with idealized test cases, both planar and on the sphere. The variational scheme displays lower errors than the finite-volume formulation for the strain rate operator but higher errors for the stress divergence operator. The variational stress divergence operator displays increased errors around the pentagonal cells of a quasi-uniform mesh, which is ameliorated with an alternate formulation for the operator. MPAS-Seaice shares the sophisticated column physics and biogeochemistry of CICE and when used with quadrilateral meshes can reproduce the results of CICE. We have used global simulations with realistic forcing to validate MPAS-Seaice against similar simulations with CICE and against observations. We find very similar results compared to CICE, with differences explained by minor differences in implementation such as with interpolation between the primary and dual meshes at coastlines. We have assessed the computational performance of the model, which, because it is unstructured, runs with 70 % of the throughput of CICE for a comparison quadrilateral simulation. The SCVT meshes used by MPAS-Seaice allow removal of equatorial model cells and flexibility in domain decomposition, improving model performance. MPAS-Seaice is the current sea-ice component of the Energy Exascale Earth System Model (E3SM).

58 GEOSCIENCES↗

Queueing Network Models for Parallel Processing of Task Systems: an Operational Approach

Computer performance modeling of possibly complex computations running on highly concurrent systems is considered. Earlier works in this area either dealt with a very simple program structure or resulted in methods with exponential complexity. An efficient procedure is developed to compute the performance measures for series-parallel-reducible task systems using queueing network models. The procedure is based on the concept of hierarchical decomposition and a new operational approach. Numerical results for three test cases are presented and compared to those of simulations.

Mak, Victor W. K.↗

Unsteady aerodynamic simulation of multiple bodies in relative motion

A prototype method for time-accurate simulation of multiple aerodynamic bodies in relative motion is presented. The method is general and features unsteady chimera domain decomposition techniques and an implicit approximately factored finite-difference procedure to solve the time-dependent thin-layer Navier-Stokes equations. The method is applied to a set of two- and three-dimensional test problems to establish spatial and temporal accuracy, quantify computational efficiency, and begin to test overall code robustness.

Meakin, Robert L.↗

Unsteady aerodynamic simulation of multiple bodies in relative motion: A prototype method

A prototype method for time-accurate simulation of multiple aerodynamic bodies in relative motion is presented. The method is general and features unsteady chimera domain decomposition techniques and an implicit approximately factored finite-difference procedure to solve the time-dependent thin-layer Navier-Stokes equations. The method is applied to a set of two- and three- dimensional test problems to establish spatial and temporal accuracy, quantify computational efficiency, and begin to test overall code robustness.

Meakin, Robert L.↗

Turbulence modeling in supersonic combusting flows

To support the National Aerospace Plane project, the RPLUS3D CFD code has been developed at NASA Lewis. The code has the ability to solve three-dimensional flowfields with finite rate combustion of hydrogen and air. The combustion processes of the hydrogen-air system are simulated by an 18-reaction path, 8-species chemical kinetic mechanism. The code uses a Lower-Upper (LU) decomposition numerical algorithm as its basis, making it a very efficient and robust code. Except for the Jacobian matrix for the implicit chemistry source terms, there is no inversion of a matrix even though it uses a fully implicit numerical algorithm. A k-epsilon (two equation) turbulence model is incorporated into the RPLUS3D code.

Chitsomboon, Tawit↗

Implementation of a kappa-epsilon turbulence model to RPLUS3D code

The RPLUS3D code has been developed at the NASA Lewis Research Center to support the National Aerospace Plane (NASP) project. The code has the ability to solve three dimensional flowfields with finite rate combustion of hydrogen and air. The combustion process of the hydrogen-air system are simulated by an 18 reaction path, 8 species chemical kinetic mechanism. The code uses a Lower-Upper (LU) decomposition numerical algorithm as its basis, making it a very efficient and robust code. Except for the Jacobian matrix for the implicit chemistry source terms, there is no inversion of a matrix even though a fully implicit numerical algorithm is used. A k-epsilon turbulence model has recently been incorporated into the code. Initial validations have been conducted for a flow over a flat plate. Results of the validation studies are shown. Some difficulties in implementing the k-epsilon equations to the code are also discussed.

Chitsomboon, Tawit↗

Progress in the Simulation of Steady and Time-Dependent Flows with 3D Parallel Unstructured Cartesian Methods

The proposed paper will present recent extensions in the development of an efficient Euler solver for adaptively-refined Cartesian meshes with embedded boundaries. The paper will focus on extensions of the basic method to include solution adaptation, time-dependent flow simulation, and arbitrary rigid domain motion. The parallel multilevel method makes use of on-the-fly parallel domain decomposition to achieve extremely good scalability on large numbers of processors, and is coupled with an automatic coarse mesh generation algorithm for efficient processing by a multigrid smoother. Numerical results are presented demonstrating parallel speed-ups of up to 435 on 512 processors. Solution-based adaptation may be keyed off truncation error estimates using tau-extrapolation or a variety of feature detection based refinement parameters. The multigrid method is extended to for time-dependent flows through the use of a dual-time approach. The extension to rigid domain motion uses an Arbitrary Lagrangian-Eulerlarian (ALE) formulation, and results will be presented for a variety of two- and three-dimensional example problems with both simple and complex geometry.

Aftosmis, M. J.↗

Limits on the UV Photodecomposition of Carbonates on Mars

The effect of UV (ultraviolet) light on the stability of calcium carbonate in a simulated martian atmosphere was experimentally investigated. Sample cells containing C-13 labeled calcite were irradiated with a Xe arc lamp in 10 mbar of simulated martian atmosphere and a quadrupole mass spectrometer was used to monitor the headspace gases for the production of (13)CO2. We found no experimental evidence of the UV photodecomposition of calcium carbonate in a simulated martian atmosphere. Extrapolating the lower limit of detection of our experimental system to an upper limit of carbonate decomposition on Mars yields a quantum efficiency of 3.5 x 10(exp -8) molecules/photon over the wavelength interval of 190-390 nm and a maximum UV photodecomposition rate of 1.2 x 10(exp -13) kg m(exp -2) s(exp -1) from a calcite surface. The maximum loss of bulk calcite due to this process would be 2.5 nm yr(exp -1). However, calcite is expected to be thermodynamically stable on the surface of Mars and potential UV photodecomposition reaction mechanisms indicate that while calcium carbonate may decompose under vacuum, it would be stable in a CO2 atmosphere. Given the expected stability of carbonate on Mars and our inability to detect carbonate decomposition, we conclude that it is unlikely that the apparent absence of carbonate on the martian surface is due to UV photo decomposition of calcite in the current environment.

Quinn, Richard↗

On the Existence of Fast Modes in Compressible Magnetohydrodynamic Turbulence

We study the existence and properties of fast magnetosonic modes in 3D compressible MHD turbulence by carrying out a number of simulations with compressible and incompressible driving conditions. We use two approaches to determine the presence of fast modes: mode decomposition based on spatial variations only and spatio-temporal 4D fast Fourier transform (4D FFT) analysis of all fluctuations. The latter method enables us to quantify fluctuations that satisfy the dispersion relation of fast modes with finite frequency. Overall, we find that the fraction of fast modes identified via the spatio-temporal 4D FFT approach in total fluctuation power is either tiny with nearly incompressible driving or ∼2% with highly compressible driving. We discuss the implications of our results for understanding the compressible fluctuations in space and astrophysical plasmas.

Zhaoming Gan↗

38th JANNAF Combustion Subcommittee Meeting

This volume, the first of two volumes, is a collection of 55 unclassified/unlimited-distribution papers which were presented at the Joint Army-Navy-NASA-Air Force (JANNAF) 38th Combustion Subcommittee (CS), 26 th Airbreathing Propulsion Subcommittee (APS), 20th Propulsion Systems Hazards Subcommittee (PSHS), and 21 Modeling and Simulation Subcommittee. The meeting was held 8-12 April 2002 at the Bayside Inn at The Sandestin Golf & Beach Resort and Eglin Air Force Base, Destin, Florida. Topics cover five major technology areas including: 1) Combustion - Propellant Combustion, Ingredient Kinetics, Metal Combustion, Decomposition Processes and Material Characterization, Rocket Motor Combustion, and Liquid & Hybrid Combustion; 2) Liquid Rocket Engines - Low Cost Hydrocarbon Liquid Rocket Engines, Liquid Propulsion Turbines, Liquid Propulsion Pumps, and Staged Combustion Injector Technology; 3) Modeling & Simulation - Development of Multi- Disciplinary RBCC Modeling, Gun Modeling, and Computational Modeling for Liquid Propellant Combustion; 4) Guns Gun Propelling Charge Design, and ETC Gun Propulsion; and 5) Airbreathing - Scramjet an Ramjet- S&T Program Overviews.

Fry, Ronald S.↗

The Surface Albedo of Sea Ice in CMIP6 and the Implications for the Surface Albedo Feedback

The Arctic has experienced rapid sea ice loss and a substantial decline in surface albedo, significantly impacting its radiation budget. Climate models from the Coupled Model Intercomparison Project (CMIP6) reproduce these changes. However, inconsistencies remain among models regarding the magnitude, spatial distribution, and seasonal patterns of Arctic surface albedo evolution. This study investigates these discrepancies by comparing model outputs with observation from the Clouds and the Earth's Radiant Energy System (CERES). We develop a decomposition method to assess the contributions of sea ice albedo, sea ice concentration, and sea ice extent to Arctic surface albedo. Over land, differences in snow cover account for the substantial inter-model spread in surface albedo, while over the ocean, sea ice albedo, concentration, and extent all contribute. Comparisons between CMIP6 and the Atmospheric Model Intercomparison Project (AMIP) simulations, which use identically prescribed sea ice concentrations, reveals considerable inter-model spread in Arctic Ocean surface albedo due to differences in sea ice albedo. Applying the decomposition method to projections shows that models predicting larger decreases in sea ice concentration and extent, especially in the Central Arctic, exhibit lower surface albedo and stronger sea ice albedo feedback. Beyond 2050, Arctic Ocean surface albedo decline is mainly influenced by sea ice extent indicating that the retreat of the ice edge is the most important process to constrain the surface albedo feedback. This study provides insights into factors contributing to the spread and changes in Arctic surface albedo and the associated sea ice albedo feedback.

Patrick C Taylor↗

Stochastic quantum Krylov protocol with double-factorized Hamiltonians

Here we propose a class of randomized quantum Krylov diagonalization (rQKD) algorithms capable of solving the eigenstate estimation problem with modest quantum resource requirements. Compared to previous real-time evolution quantum Krylov subspace methods, our approach expresses the time evolution operator e –i$\widehat{H}$$\tau$ as a linear combination of unitaries and subsequently uses a stochastic sampling procedure to reduce circuit depth requirements. While our methodology applies to any Hamiltonian with fast-forwardable subcomponents, we focus on its application to the explicitly double-factorized electronic-structure Hamiltonian. To demonstrate the potential of the proposed rQKD algorithm on near-term quantum devices, we provide numerical benchmarks for a variety of molecular systems with circuit-based state-vector simulators including the effects of sampling noise, achieving ground-state energy errors of less than 1 kcal mol -1 with circuit depths orders of magnitude shallower than those required for low-rank deterministic Trotter-Suzuki decompositions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Aerodynamic Data Fusion Toward the Digital Twin Paradigm

This paper considers the fusion of two aerodynamic data sets originating from differing types of physical or computer experiments. This paper specifically addresses the fusion of 1) noisy and in-complete fields from wind-tunnel measurements and 2) deterministic but biased fields from numerical simulations. These two data sources are fused in order to estimate the true field that best matches measured quantities that serve as the ground truth. For example, two sources of pressure fields about an aircraft are fused based on measured forces and moments from a wind-tunnel experiment. A fundamental challenge in this problem is that the true field is unknown and cannot be estimated with 100% certainty. A Bayesian framework is employed to infer the true fields conditioned on measured quantities of interest; essentially a statistical correction to the data is performed. The fused data may then be used to construct more accurate surrogate models suitable for early stages of aerospace design. An extension of the proper orthogonal decomposition with constraints is also introduced to solve the same problem. In this work, both methods are demonstrated on fusing the pressure distributions for flow past the RAE2822 airfoil and the Common Research Model wing at transonic conditions. Comparison of both methods reveals that the Bayesian method is more robust when data are scarce and capable of also accounting for uncertainties in the data. Furthermore, given adequate data, the proper-orthogonal-decomposition-based and Bayesian approaches lead to surprisingly similar results.

42 ENGINEERING↗