Search NASA⌕ Search

SEARCH · Search NASA

Results for “Conservative variables”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Mixing of Multiple Jets With a Confined Subsonic Crossflow

Results from a recently completed enhanced mixing program are summarized in the two technical papers. These studies were parts of a High Speed Research (HSR)-supported joint Government/industry/university program that involved, in addition to the NASA Lewis Research Center, researchers at United Technologies Research Center, Allison Engine Company, CFD Research Corporation, and the University of California, Irvine. The studies investigated the mixing of jets injected normal to a confined subsonic mainsteam in both rectangular and cylindrical ducts. Experimental and computational studies were performed in both nonreacting and reacting flows. The orifice geometries and flow conditions were selected as typical of the complex three-dimensional flows in the combustion chambers in low-emission gas turbine engines. The principal conclusion from both the experiments and modeling was that the momentum-flux ratio J and orifice spacing S/H were the most significant flow and geometry variables, respectively. Conserved scalar distributions were similar-independent of reaction, orifice diameter H/d, and shape-when the orifice spacing and the square root of the momentum-flux ratio were inversely proportional. Jet penetration was critical, and penetration decreased as either momentum-flux ratio or orifice spacing decreased. We found that planar averages must be considered in context with the distributions. The mass-flow ratios and the orifices investigated were often very large. The jet-to-mainstream mass-flow ratio was varied from significantly less than 1 to greater than 1. The orifice-area to mainstream-cross-sectional-area was varied from approx. 0 to 0.5, and the axial planes of interest were often just downstream of the orifice trailing edge. Three-dimensional flow was a key part of efficient mixing and was observed for all configurations. As an example of the results, the accompanying figure shows the effects of different rates of mass addition on the opposite walls of a rectangular duct.

Holdeman, James D.↗
NASA NTRS2006

Accuracy Quantification of the Loci-CHEM Code for Chamber Wall Heat Transfer in a GO2/GH2 Single Element Injector Model Problem

A robust rocket engine combustor design and development process must include tools which can accurately predict the multi-dimensional thermal environments imposed on solid surfaces by the hot combustion products. Currently, empirical methods used in the design process are typically one dimensional and do not adequately account for the heat flux rise rate in the near-injector region of the chamber. Computational Fluid Dynamics holds promise to meet the design tool requirement, but requires accuracy quantification, or validation, before it can be confidently applied in the design process. This effort presents the beginning of such a validation process for the Loci-CHEM CFD code. The model problem examined here is a gaseous oxygen (GO2)/gaseous hydrogen (GH2) shear coaxial single element injector operating at a chamber pressure of 5.42 MPa. The GO2/GH2 propellant combination in this geometry represents one the simplest rocket model problems and is thus foundational to subsequent validation efforts for more complex injectors. Multiple steady state solutions have been produced with Loci-CHEM employing different hybrid grids and two-equation turbulence models. Iterative convergence for each solution is demonstrated via mass conservation, flow variable monitoring at discrete flow field locations as a function of solution iteration and overall residual performance. A baseline hybrid was used and then locally refined to demonstrate grid convergence. Solutions were obtained with three variations of the k-omega turbulence model.

West, Jeff↗
NASA NTRS2006

Accuracy Quantification of the Loci-CHEM Code for Chamber Wall Heat Fluxes in a G02/GH2 Single Element Injector Model Problem

A robust rocket engine combustor design and development process must include tools which can accurately predict the multi-dimensional thermal environments imposed on solid surfaces by the hot combustion products. Currently, empirical methods used in the design process are typically one dimensional and do not adequately account for the heat flux rise rate in the near-injector region of the chamber. Computational Fluid Dynamics holds promise to meet the design tool requirement, but requires accuracy quantification, or validation, before it can be confidently applied in the design process. This effort presents the beginning of such a validation process for the Loci- CHEM CPD code. The model problem examined here is a gaseous oxygen (GO2)/gaseous hydrogen (GH2) shear coaxial single element injector operating at a chamber pressure of 5.42 MPa. The GO2/GH2 propellant combination in this geometry represents one the simplest rocket model problems and is thus foundational to subsequent validation efforts for more complex injectors. Multiple steady state solutions have been produced with Loci-CHEM employing different hybrid grids and two-equation turbulence models. Iterative convergence for each solution is demonstrated via mass conservation, flow variable monitoring at discrete flow field locations as a function of solution iteration and overall residual performance. A baseline hybrid grid was used and then locally refined to demonstrate grid convergence. Solutions were also obtained with three variations of the k-omega turbulence model.

West, Jeff↗
NASA NTRS2012

Calculating Soil Wetness, Evapotranspiration and Carbon Cycle Processes Over Large Grid Areas Using a New Scaling Technique

Soil wetness typically shows great spatial variability over the length scales of general circulation model (GCM) grid areas (approx 100 km ), and the functions relating evapotranspiration and photosynthetic rate to local-scale (approx 1 m) soil wetness are highly non-linear. Soil respiration is also highly dependent on very small-scale variations in soil wetness. We therefore expect significant inaccuracies whenever we insert a single grid area-average soil wetness value into a function to calculate any of these rates for the grid area. For the particular case of evapotranspiration., this method - use of a grid-averaged soil wetness value - can also provoke severe oscillations in the evapotranspiration rate and soil wetness under some conditions. A method is presented whereby the probability distribution timction(pdf) for soil wetness within a grid area is represented by binning. and numerical integration of the binned pdf is performed to provide a spatially-integrated wetness stress term for the whole grid area, which then permits calculation of grid area fluxes in a single operation. The method is very accurate when 10 or more bins are used, can deal realistically with spatially variable precipitation, conserves moisture exactly and allows for precise modification of the soil wetness pdf after every time step. The method could also be applied to other ecological problems where small-scale processes must be area-integrated, or upscaled, to estimate fluxes over large areas, for example in treatments of the terrestrial carbon budget or trace gas generation.

Sellers, Piers↗
NASA NTRS2026

Spacecraft Disposal Rosetta Stone: Parametric Tool for Orbital Lifetime, Disposal, and Cost Assessment

This Technical Memorandum documents a simplified, parametric method for evaluating spacecraft orbital lifetime, disposal compliance, and disposal-related cost impacts during early mission formulation and preliminary design. The method captures the dominant drivers of orbital decay—effective ballistic coefficient, operating altitude, and solar-cycle variability—using conservative bounding assumptions. Solar maximum conditions are used to bound achievable mission lifetime, while solar minimum conditions are used to bound disposal timelines and compliance with orbital debris requirements. A single tabulated dataset provides orbital lifetime under both solar-cycle extremes together with representative disposal ΔV required to ensure compliant disposal, enabling rapid assessment of disposal feasibility, cost sensitivity, and system-level impacts prior to committing to higher-fidelity analyses.

Orbital debris mitigation↗
NASA NTRS1993

Variable-Volume Flushing (V-VF) device for water conservation in toilets

Thirty five percent of residential indoor water used is flushed down the toilet. Five out of six flushes are for liquid waste only, which requires only a fraction of the water needed for solid waste. Designers of current low-flush toilets (3.5-gal. flush) and ultra-low-flush toilets (1.5-gal. flush) did not consider the vastly reduced amount of water needed to flush liquid waste versus solid waste. Consequently, these toilets are less practical than desired and can be improved upon for water conservation. This paper describes a variable-volume flushing (V-VF) device that is more reliable than the currently used flushing devices (it will not leak), is simple, more economical, and more water conserving (allowing one to choose the amount of water to use for flushing solid and liquid waste).

Jasper, Louis J., Jr.↗
DOE OSTI2021

The minor pilin PilV provides a conserved adhesion site throughout the antigenically variable meningococcal type IV pilus

Type IV pili (T4P) are among the most widespread adhesive factors in prokaryotes. In pathogenic Neisseria, the major pilin, which provides the structural framework for the filamentous T4P, undergoes antigenic variation allowing the bacteria to evade the humoral immune response without impacting host-cell adhesion. Here, we show that a minor pilin, PilV, is distributed throughout the pilus and contributes to Neisseria meningitidis adhesion and that antibodies to PilV block meningococcal adhesion in vivo. Our results provide a mechanism whereby N. meningitidis varies its immunodominant major pilin to escape antibody recognition while maintaining conserved sites throughout the pilus for host receptor binding. They further suggest a strategy to prevent or block deadly N. meningitidis infections by targeting this minor pilin.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1981

Numerical approximation of boundary conditions with applications to inviscid equations of gas dynamics

A comprehensive overview of the state of the art of well-posedness and stability analysis of difference approximations for initial boundary value problems of the hyperbolic type is presented. The applicability of recent theoretical development to practical calculations for nonlinear gas dynamics is examined. The one dimensional inviscid gas dynamics equations in conservation law form are selected for numerical experiments. The class of implicit schemes developed from linear multistep methods in ordinary differential equations is chosen and the use of linear extrapolation as an explicit or implicit boundary scheme is emphasized. Specification of boundary data in the primitive variables and computation in terms of the conservative variables in the interior is discussed. Some numerical examples for the quasi-one-dimensional nozzle are given.

Yee, H. C.↗
DOE OSTI2023

Predicting variable gene content in Escherichia coli using conserved genes

Having the ability to predict the protein-encoding gene content of an incomplete genome or metagenome-assembled genome is important for a variety of bioinformatic tasks. In this study, as a proof of concept, we built machine learning classifiers for predicting variable gene content in Escherichia coli genomes using only the nucleotide k-mers from a set of 100 conserved genes as features. Protein families were used to define orthologs, and a single classifier was built for predicting the presence or absence of each protein family occurring in 10%–90% of all E. coli genomes. The resulting set of 3,259 extreme gradient boosting classifiers had a per-genome average macro F1 score of 0.944 [0.943–0.945, 95% CI]. We show that the F1 scores are stable across multi-locus sequence types and that the trend can be recapitulated by sampling a smaller number of core genes or diverse input genomes. Surprisingly, the presence or absence of poorly annotated proteins, including “hypothetical proteins” was accurately predicted (F1 = 0.902 [0.898–0.906, 95% CI]). Models for proteins with horizontal gene transfer-related functions had slightly lower F1 scores but were still accurate (F1s = 0.895, 0.872, 0.824, and 0.841 for transposon, phage, plasmid, and antimicrobial resistance-related functions, respectively). Finally, using a holdout set of 419 diverse E. coli genomes that were isolated from freshwater environmental sources, we observed an average per-genome F1 score of 0.880 [0.876–0.883, 95% CI], demonstrating the extensibility of the models. Overall, this study provides a framework for predicting variable gene content using a limited amount of input sequence data.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2003

A Discrete Constraint for Entropy Conservation and Sound Waves in Cloud-Resolving Modeling

Ideal cloud-resolving models contain little-accumulative errors. When their domain is so large that synoptic large-scale circulations are accommodated, they can be used for the simulation of the interaction between convective clouds and the large-scale circulations. This paper sets up a framework for the models, using moist entropy as a prognostic variable and employing conservative numerical schemes. The models possess no accumulative errors of thermodynamic variables when they comply with a discrete constraint on entropy conservation and sound waves. Alternatively speaking, the discrete constraint is related to the correct representation of the large-scale convergence and advection of moist entropy. Since air density is involved in entropy conservation and sound waves, the challenge is how to compute sound waves efficiently under the constraint. To address the challenge, a compensation method is introduced on the basis of a reference isothermal atmosphere whose governing equations are solved analytically. Stability analysis and numerical experiments show that the method allows the models to integrate efficiently with a large time step.

Zeng, Xi-Ping↗
DOE OSTI2021

Machine Learning for Conservative-to-Primitive in Relativistic Hydrodynamics

The numerical solution of relativistic hydrodynamics equations in conservative form requires root-finding algorithms that invert the conservative-to-primitive variables map. These algorithms employ the equation of state of the fluid and can be computationally demanding for applications involving sophisticated microphysics models, such as those required to calculate accurate gravitational wave signals in numerical relativity simulations of binary neutron stars. This work explores the use of machine learning methods to speed up the recovery of primitives in relativistic hydrodynamics. Artificial neural networks are trained to replace either the interpolations of a tabulated equation of state or directly the conservative-to-primitive map. The application of these neural networks to simple benchmark problems shows that both approaches improve over traditional root finders with tabular equation-of-state and multi-dimensional interpolations. In particular, the neural networks for the conservative-to-primitive map accelerate the variable recovery by more than an order of magnitude over standard methods while maintaining accuracy. Neural networks are thus an interesting option to improve the speed and robustness of relativistic hydrodynamics algorithms.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2001

High Order Difference Method for Low Mach Number Aeroacoustics

A high order finite difference method with improved accuracy and stability properties for computational aeroacoustics (CAA) at low Mach numbers is proposed. The Euler equations are split into a conservative and a symmetric non- conservative portion to allow the derivation of a generalized energy estimate. Since the symmetrization is based on entropy variables, that splitting of the flux derivatives is referred to as entropy splitting. Its discretization by high order central differences was found to need less numerical dissipation than conventional conservative schemes. Owing to the large disparity of acoustic and stagnation quantities in low Mach number aeroacoustics, the split Euler equations are formulated in perturbation form. The unknowns are the small changes of the conservative variables with respect to their large stagnation values. All nonlinearities and the conservation form of the conservative portion of the split flux derivatives can be retained, while cancellation errors are avoided with its discretization opposed to the conventional conservative form. The finite difference method is third-order accurate at the boundary and the conventional central sixth-order accurate stencil in the interior. The difference operator satisfies the summation by parts property analogous to the integration by parts in the continuous energy estimate. Thus, strict stability of the difference method follows automatically. Spurious high frequency oscillations are suppressed by a characteristic-based filter similar to but without limiter. The time derivative is approximated by a 4-stage low-storage second-order explicit Runge-Kutta method. The method has been applied to simulate vortex sound at low Mach numbers. We consider the Kirchhoff vortex, which is an elliptical patch of constant vorticity rotating with constant angular frequency in irrotational flow. The acoustic pressure generated by the Kirchhoff vortex is governed by the 2D Helmholtz equation, which can be solved analytically using separation of variables.

Mueller, B.↗
NASA NTRS2008

Comparison of TID Effects in Space-Like Variable Dose Rates and Constant Dose Rates

The degradation of the LM193 dual voltage comparator has been studied at different TID dose rate profiles, including several different constant dose rates and a variable dose rate that simulates the behavior of a solar flare. A comparison of results following constant dose rate vs. variable dose rates is made to explore how well the constant dose rates used for typical part testing predict the performance during a simulated space-like mission. Testing at a constant dose rate equal to the lowest dose rate seen during the simulated flare provides an extremely conservative estimate of the overall amount of degradation. A constant dose rate equal to the average dose rate is also more conservative than the variable rate. It appears that, for this part, weighting the dose rates by the amount of total dose received at each rate (rather than the amount of time at each dose rate) results in an average rate that produces an amount of degradation that is a reasonable approximation to that received by the variable rate.

Solar Flare↗
NASA NTRS1983

Flow in coronal loops with a mass source

This research studies the flow of plasma inside a coronal loop in which an injection of plasma through the lateral surface is permitted. The flow is assumed steady and polytropic. The problem covers two cases: (a) upflow at one footpoint, downflow at the other; (b) downflow at both footpoints. The first case can be shown to be quite similar to that of a mass-conserving flow with variable cross section; the second, instead, is characterized by solutions with a different type of topology; its main new feature is the obvious fact that all the solutions pass through a single point going from negative to positive velocities. In this second case the density ratio between footpoints and top can be much smaller than in a mass conserving flow. This can explain some properties of observed loops.

Noci, G.↗
DOE OSTI2025

Geometric adiabatic angle in anisotropic oscillators

We discuss a classical anisotropic oscillator and the Foucault pendulum as examples illustrating non-conservation of action variables in integrable classical mechanical systems with adiabatically slow evolution. We also emphasize the importance of the mass parameter of a harmonic oscillator, alongside its frequency, in explicitly time-dependent situations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A moving discontinuous Galerkin finite element method with interface condition enforcement for compressible flows

A variation of moving discontinuous Galerkin finite element method with interface condition enforcement (MDG-ICE) is developed for solving the compressible Euler equations. The MDG-ICE method, originating from the work of Corrigan et al. [1], [2], [3], [4], is based on the space-time DG formulation, where both flow field and grid geometry are considered as independent variables and the conservation laws are enforced both on discrete elements and element interfaces. The element conservation laws are solved in the standard discontinuous solution space to determine conservative quantities, while the interface conservation is enforced using a variational formulation in a continuous space to determine discrete grid geometry. The resulting over-determined system of nonlinear equations arising from the MDG-ICE formulation can then be solved in a least-squares sense, leading to an unconstrained nonlinear least-squares problem that is regularized and solved by Levenberg-Marquardt method. A number of numerical experiments for both 1D unsteady and 2D steady state compressible flow problems are conducted to assess the accuracy and robustness of the MDG-ICE method. Numerical results obtained indicate that the MDG-ICE method is able to implicitly detect and track all types of discontinuities via interface conservation enforcement and satisfy the conservation law on both elements and interfaces via grid movement and grid management, demonstrating that an exponential rate of convergence for Sod and Lax-Harden shock tube problems can be achieved and highly accurate solutions without overheating to both double-rarefaction wave and Noh problems can be obtained.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Neoclassical transport in strong gradient regions of large aspect ratio tokamaks

We present a new neoclassical transport model for large aspect ratio tokamaks where the gradient scale lengths are of the size of the ion poloidal gyroradius. Previous work on neoclassical transport across transport barriers assumed large density and potential gradients but a small temperature gradient, or neglected the gradient of the mean parallel flow. Using large aspect ratio and low collisionality expansions, we relax these restrictive assumptions. We define a new set of variables based on conserved quantities, which simplifies the drift kinetic equation whilst keeping strong gradients, and derive equations describing the transport of particles, parallel momentum and energy by ions in the banana regime. The poloidally varying parts of density and electric potential are included. Studying contributions from both passing and trapped particles, we show that the resulting transport is dominated by trapped particles. We find that a non-zero neoclassical particle flux requires parallel momentum input which could be provided through interaction with turbulence or impurities. We derive upper and lower bounds for the energy flux across a transport barrier in both temperature and density and present example profiles and fluxes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Disentangling Alzheimer’s disease neurodegeneration from typical brain ageing using machine learning

Abstract Neuroimaging biomarkers that distinguish between changes due to typical brain ageing and Alzheimer’s disease are valuable for determining how much each contributes to cognitive decline. Supervised machine learning models can derive multivariate patterns of brain change related to the two processes, including the Spatial Patterns of Atrophy for Recognition of Alzheimer’s Disease (SPARE-AD) and of Brain Aging (SPARE-BA) scores investigated herein. However, the substantial overlap between brain regions affected in the two processes confounds measuring them independently. We present a methodology, and associated results, towards disentangling the two. T1-weighted MRI scans of 4054 participants (48–95 years) with Alzheimer’s disease, mild cognitive impairment (MCI), or cognitively normal (CN) diagnoses from the Imaging-based coordinate SysTem for AGIng and NeurodeGenerative diseases (iSTAGING) consortium were analysed. Multiple sets of SPARE scores were investigated, in order to probe imaging signatures of certain clinically or molecularly defined sub-cohorts. First, a subset of clinical Alzheimer’s disease patients (n = 718) and age- and sex-matched CN adults (n = 718) were selected based purely on clinical diagnoses to train SPARE-BA1 (regression of age using CN individuals) and SPARE-AD1 (classification of CN versus Alzheimer’s disease) models. Second, analogous groups were selected based on clinical and molecular markers to train SPARE-BA2 and SPARE-AD2 models: amyloid-positive Alzheimer’s disease continuum group (n = 718; consisting of amyloid-positive Alzheimer’s disease, amyloid-positive MCI, amyloid- and tau-positive CN individuals) and amyloid-negative CN group (n = 718). Finally, the combined group of the Alzheimer’s disease continuum and amyloid-negative CN individuals was used to train SPARE-BA3 model, with the intention to estimate brain age regardless of Alzheimer’s disease-related brain changes. The disentangled SPARE models, SPARE-AD2 and SPARE-BA3, derived brain patterns that were more specific to the two types of brain changes. The correlation between the SPARE-BA Gap (SPARE-BA minus chronological age) and SPARE-AD was significantly reduced after the decoupling (r = 0.56–0.06). The correlation of disentangled SPARE-AD was non-inferior to amyloid- and tau-related measurements and to the number of APOE ε4 alleles but was lower to Alzheimer’s disease-related psychometric test scores, suggesting the contribution of advanced brain ageing to the latter. The disentangled SPARE-BA was consistently less correlated with Alzheimer’s disease-related clinical, molecular and genetic variables. By employing conservative molecular diagnoses and introducing Alzheimer’s disease continuum cases to the SPARE-BA model training, we achieved more dissociable neuroanatomical biomarkers of typical brain ageing and Alzheimer’s disease.

Hwang, Gyujoon↗