Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “current conservation”

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 343 records · Page 19

Time-Accurate Local Time Stepping and High-Order Time CESE Methods for Multi-Dimensional Flows Using Unstructured Meshes

With the wide availability of affordable multiple-core parallel supercomputers, next generation numerical simulations of flow physics are being focused on unsteady computations for problems involving multiple time scales and multiple physics. These simulations require higher solution accuracy than most algorithms and computational fluid dynamics codes currently available. This paper focuses on the developmental effort for high-fidelity multi-dimensional, unstructured-mesh flow solvers using the space-time conservation element, solution element (CESE) framework. Two approaches have been investigated in this research in order to provide high-accuracy, cross-cutting numerical simulations for a variety of flow regimes: 1) time-accurate local time stepping and 2) highorder CESE method. The first approach utilizes consistent numerical formulations in the space-time flux integration to preserve temporal conservation across the cells with different marching time steps. Such approach relieves the stringent time step constraint associated with the smallest time step in the computational domain while preserving temporal accuracy for all the cells. For flows involving multiple scales, both numerical accuracy and efficiency can be significantly enhanced. The second approach extends the current CESE solver to higher-order accuracy. Unlike other existing explicit high-order methods for unstructured meshes, the CESE framework maintains a CFL condition of one for arbitrarily high-order formulations while retaining the same compact stencil as its second-order counterpart. For large-scale unsteady computations, this feature substantially enhances numerical efficiency. Numerical formulations and validations using benchmark problems are discussed in this paper along with realistic examples.

Chang, Chau-Lyan↗

Southern Bhutan Ecological Forecasting Tech Paper: Modeling Asian Elephant (Elephas maximus) Habitat Suitability along the Southern Bhutan Border with NASA Earth Observations

Asian elephants (Elephas maximus) are a flagship species essential for the functioning of forest ecosystems, and they also have cultural significance in Bhutan. Elephants receive the highest legal protection as listed under Schedule I of the Bhutan Forests and Nature Conservation Act, 1995. Yet, they face threats of extinction due to poaching for ivory as well as the loss and fragmentation of their habitat. Due to the recent clearing of forests and the growing populations in these areas, there has been an increase in incidents of human-elephant conflict. These conflicts have been detrimental to farmer’s annual harvests and livelihoods and have led to retaliatory killing and injury of elephants in southern Bhutan. The DEVELOP Southern Bhutan Ecological Forecasting team partnered with the Bhutan Foundation and Bhutan Tiger Center to help address this problem. The team integrated NASA Earth observations, including Landsat 5 Thematic Mapper (TM), Landsat 8 Operational Land Imager (OLI), Terra Moderate Resolution Imaging Spectroradiometer (MODIS), and the Shuttle Radar Topography Mission (SRTM) to acquire information on land cover change and elephant habitat suitability along the southern border of Bhutan. The team utilized Esri ArcGIS Pro and Software for Assisted Habitat Modeling (SAHM) for data analysis, modeling, and visualization. The team used elephant occurrence data and environmental variables to model current habitat suitability for migrating elephant populations. This analysis provided partners with maps to inform decisions about the placement and conservation of elephant corridors and helped build their capacity to use satellite data for future studies and project planning.

Palchen Wangchuk↗

Southern Bhutan Ecological Forecasting Project Summary - Modeling Asian Elephant (Elephas Maximus) Habitat Suitability along the Southern Bhutan Border with NASA Earth Observations

Asian elephants (Elephas maximus) are a flagship species essential for the functioning of forest ecosystems, and they also have cultural significance in Bhutan. Elephants receive the highest legal protection as listed under Schedule I of the Bhutan Forests and Nature Conservation Act, 1995. Yet, they face threats of extinction due to poaching for ivory as well as the loss and fragmentation of their habitat. Due to the recent clearing of forests and the growing populations in these areas, there has been an increase in incidents of human-elephant conflict. These conflicts have been detrimental to farmer’s annual harvests and livelihoods and have led to retaliatory killing and injury of elephants in southern Bhutan. The DEVELOP Southern Bhutan Ecological Forecasting team partnered with the Bhutan Foundation and Bhutan Tiger Center to help address this problem. The team integrated NASA Earth observations, including Landsat 5 Thematic Mapper (TM), Landsat 8 Operational Land Imager (OLI), Terra Moderate Resolution Imaging Spectroradiometer (MODIS), and the Shuttle Radar Topography Mission (SRTM) to acquire information on land cover change and elephant habitat suitability along the southern border of Bhutan. The team utilized Esri ArcGIS Pro and Software for Assisted Habitat Modeling (SAHM) for data analysis, modeling, and visualization. The team used elephant occurrence data and environmental variables to model current habitat suitability for migrating elephant populations. This analysis provided partners with maps to inform decisions about the placement and conservation of elephant corridors and helped build their capacity to use satellite data for future studies and project planning.

DEVELOP Project Summary↗

An Evaluation of the Impact of Glycolate and Glycolate Mitigation on the Defense Waste Processing Facility Recycle Diversion Project Flowsheet

The Savannah River Site (SRS) Defense Waste Processing Facility (DWPF) processes radioactive High Level Waste (HLW) sludge solids from the Concentration, Storage, and Transfer Facilities (CSTF); which includes the SRS Tank Farm and Evaporator facilities, and a concentrated Cs-137 laden stream and a Monosodium Titanate (MST) and sludge solids stream from the SRS Salt Waste Processing Facility (SWPF). The waste is chemically adjusted with acids and reductant (currently with 50 wt.% nitric acid and ~90 wt.% formic acid, but eventually formic acid will be substituted with ~70 wt.% glycolic acid), and frit is added so that a durable, borosilicate glass waste form can be produced when the material is vitrified in the melter. As a result of the evaporation of water during both the melter feed preparation and the melter feed vitrification steps in DWPF, a recycle waste stream is generated, neutralized, and sent back to the CSTF. The recycle waste is a dilute aqueous stream originating from the collection of condensate liquids containing some minor sludge, MST, and frit solids and other waste components resulting from melter feed entrainment during foamover events and transfer of volatile species into the condensate. The recycle stream volume is significant and is expected to approach 3 million gallons per year once SWPF reaches full operation, requiring the use of multiple large CSTF tanks for storage. The recycle waste is currently collected in the SRS Tank Farm and periodically evaporated in the 242-16H (2H) Evaporator to conserve storage space.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

User Manual - HydraGNN v5.0: Distributed Implementation of Multi-Tasking Graph Neural Networks

This document serves as the user manual for HydraGNN v5.0, a scalable graph neural network (GNN) architecture for simultaneous prediction of multiple target properties using multi-task learning (MTL). This version of HydraGNN has been developed primarily to support the development, training, and deployment of predictive graph-based deep learning (DL) models for atomistic materials modeling. HydraGNN is templated over 13 message-passing policies, including invariant models (GIN, PNA, PNAPlus, GAT, MFC, CGCNN, SAGE, SchNet, DimeNet) and equivariant models (EGNN, PNAEq, PAINN, MACE), and supports distributed training via distributed data parallelism (DDP), DeepSpeed, and Fully Sharded Data Parallelism (FSDP) on leadership-class supercomputers. Although HydraGNN can be applied to problems beyond atomistic materials modeling, its current use is confined to homogeneous graphs. Additional capabilities include machine-learned interatomic potentials with energy-conserving forces, General, Powerful, and Scalable Graph Transformer (GraphGPS) global attention, periodic boundary conditions, hyperparameter optimization, mixed-precision training, and uncertainty quantification.

97 MATHEMATICS AND COMPUTING↗

National and international capital shortages slowing, perhaps preventing, American economic growth and leading to a decline in the domestic standard of living

It is shown that many of the basic industries that the U.S. has relied upon in the past for economic growth and development are now so obsolete, so old, and so technologically inferior to that of foreign competitors that the U.S. is losing its international competitive position. The most conservative estimate suggests that it will require $325 billion between now and 1982 merely to meet existing and currently anticipated pollution requirements and that it would take an additional $197 billion to replace outmoded existing facilities.

Source record↗

Influence of land-surface evapotranspiration on the earth's climate

Land-surface evapotranspiration is shown to strongly influence global fields of rainfall, temperature and motion by calculations using a numerical model of the atmosphere, confirming the general belief in the importance of evapotranspiration-producing surface vegetation for the earth's climate. The current version of the Goddard Laboratory atmospheric general circulation model is used in the present experiment, in which conservation equations for mass, momentum, moisture and energy are expressed in finite-difference form for a spherical grid to calculate (1) surface pressure field evolution, and (2) the wind, temperature, and water vapor fields at nine levels between the surface and a 20 km height.

Shukla, J.↗

The Decompression Sickness and Venous Gas Emboli Consequences of Air Breaks During 100% Oxygen Prebreathe

Not enough is known about the increased risk of hypobaric decompression sickness (DCS) and production of venous (VGE) and arterial (AGE) gas emboli following an air break in an otherwise normal 100% resting oxygen (O2) prebreathe (PB), and certainly a break in PB when exercise is used to accelerate nitrogen (N2) elimination from the tissues. Current Aeromedical Flight Rules at the Johnson Space Center about additional PB payback times are untested, possibly too conservative, and therefore not optimized for operational use.

Conkin, Johnny↗

The Decompression Sickness and Venous Gas Emboli Consequences of Air Breaks During 100% Oxygen Prebreathe

Not enough is known about the increased risk of hypobaric decompression sickness (DCS) and production of venous (VGE) and arterial (AGE) gas emboli following an air break in an otherwise normal 100% resting oxygen (O2) prebreathe (PB), and certainly a break in PB when exercise is used to accelerate nitrogen (N2) elimination from the tissues. Current Aeromedical Flight Rules at the Johnson Space Center about additional PB payback times are untested, possibly too conservative, and therefore not optimized for operational use. A 10 min air break at 90 min into a 120 min PB that includes initial dual-cycle ergometry for 10 min will show a measurable increase in the risk of DCS and VGE after ascent to 4.3 psia compared to a 10 min break at 15 min into the PB, or when there is no break in PB. Data collection with humans begins in 2005, but here we first evaluate the hypothesis using three models of tissue N2 kinetics: Model I is a simple single half-time compartment exponential model, Model II is a three compartment half-time exponential model, and Model III is a variable half-time compartment model where the percentage of maximum O2 consumption for the subject during dual-cycle ergometry exercise defines the half-time compartment. Model I with large rate constants to simulate an exercise effect always showed a late break in PB had the greatest consequence. Model II showed an early break had the greatest consequence. Model III showed there was no difference between early or late break in exercise PB. Only one of these outcomes will be observed when humans are tested. Our results will favor one of these models, and so advance our understanding of tissue N2 kinetics, and of altitude DCS after an air break in PB.

Conkin, J.↗

Global Aspects of Charged Particle Motion in Axially Symmetric Multipole Magnetic Fields

The motion of a single charged particle in the space outside of a compact region of steady currents is investigated. The charged particle is assumed to produce negligible electromagnetic radiation, so that its energy is conserved. The source of the magnetic field is represented as a point multipole. After a general description, attention is focused on magnetic fields with axial symmetry. Lagrangian dynamical theory is utilized to identify constants of the motion as well as the equations of motion themselves. The qualitative method of Stonner is used to examine charged particle motion in axisymmetric multipole fields of all orders. Although the equations of motion generally have no analytical solutions and must be integrated numerically to produce a specific orbit, a topological examination of dynamics is possible, and can be used, d la Stonner, to completely describe the global aspects of the motion of a single charged particle in a space with an axisymmetric multipole magnetic field.

Shebalin, John V.↗

The Impact of Back-Sputtered Carbon on the Accelerator Grid Wear Rates of the NEXT and NSTAR Ion Thrusters

A study was conducted to quantify the impact of back-sputtered carbon on the downstream accelerator grid erosion rates of the NEXT (NASA's Evolutionary Xenon Thruster) Long Duration Test (LDT1). A similar analysis that was conducted for the NSTAR (NASA's Solar Electric Propulsion Technology Applications Readiness Program) Life Demonstration Test (LDT2) was used as a foundation for the analysis developed herein. A new carbon surface coverage model was developed that accounted for multiple carbon adlayers before complete surface coverage is achieved. The resulting model requires knowledge of more model inputs, so they were conservatively estimated using the results of past thin film sputtering studies and particle reflection predictions. In addition, accelerator current densities across the grid were rigorously determined using an ion optics code to determine accelerator current distributions and an algorithm to determine beam current densities along a grid using downstream measurements. The improved analysis was applied to the NSTAR test results for evaluation. The improved analysis demonstrated that the impact of back-sputtered carbon on pit and groove wear rate for the NSTAR LDT2 was negligible throughout most of eroded grid radius. The improved analysis also predicted the accelerator current density for transition from net erosion to net deposition considerably more accurately than the original analysis. The improved analysis was used to estimate the impact of back-sputtered carbon on the accelerator grid pit and groove wear rate of the NEXT Long Duration Test (LDT1). Unlike the NSTAR analysis, the NEXT analysis was more challenging because the thruster was operated for extended durations at various operating conditions and was unavailable for measurements because the test is ongoing. As a result, the NEXT LDT1 estimates presented herein are considered preliminary until the results of future posttest analyses are incorporated. The worst-case impact of carbon back-sputtering was determined to be the full power operating condition, but the maximum impact of back-sputtered carbon was only a four percent reduction in wear rate. As a result, back-sputtered carbon is estimated to have an insignificant impact on the first failure mode of the NEXT LDT at all operating conditions.

Soulas, George C.↗

The Impact of Back-Sputtered Carbon on the Accelerator Grid Wear Rates of the NEXT and NSTAR Ion Thrusters

A study was conducted to quantify the impact of back-sputtered carbon on the downstream accelerator grid erosion rates of the NASA's Evolutionary Xenon Thruster (NEXT) Long Duration Test (LDT1). A similar analysis that was conducted for the NASA's Solar Electric Propulsion Technology Applications Readiness Program (NSTAR) Life Demonstration Test (LDT2) was used as a foundation for the analysis developed herein. A new carbon surface coverage model was developed that accounted for multiple carbon adlayers before complete surface coverage is achieved. The resulting model requires knowledge of more model inputs, so they were conservatively estimated using the results of past thin film sputtering studies and particle reflection predictions. In addition, accelerator current densities across the grid were rigorously determined using an ion optics code to determine accelerator current distributions and an algorithm to determine beam current densities along a grid using downstream measurements. The improved analysis was applied to the NSTAR test results for evaluation. The improved analysis demonstrated that the impact of back-sputtered carbon on pit and groove wear rate for the NSTAR LDT2 was negligible throughout most of eroded grid radius. The improved analysis also predicted the accelerator current density for transition from net erosion to net deposition considerably more accurately than the original analysis. The improved analysis was used to estimate the impact of back-sputtered carbon on the accelerator grid pit and groove wear rate of the NEXT Long Duration Test (LDT1). Unlike the NSTAR analysis, the NEXT analysis was more challenging because the thruster was operated for extended durations at various operating conditions and was unavailable for measurements because the test is ongoing. As a result, the NEXT LDT1 estimates presented herein are considered preliminary until the results of future post-test analyses are incorporated. The worst-case impact of carbon back-sputtering was determined to be the full power operating condition, but the maximum impact of back-sputtered carbon was only a 4 percent reduction in wear rate. As a result, back-sputtered carbon is estimated to have an insignificant impact on the first failure mode of the NEXT LDT1 at all operating conditions.

Soulas, George C.↗

Modeling Two-Phase Flow and Vapor Cycles Using the Generalized Fluid System Simulation Program

This work presents three new applications for the general purpose fluid network solver code GFSSP developed at NASA's Marshall Space Flight Center: (1) cooling tower, (2) vapor-compression refrigeration system, and (3) vapor-expansion power generation system. These systems are widely used across engineering disciplines in a variety of energy systems, and these models expand the capabilities and the use of GFSSP to include fluids and features that are not part of its present set of provided examples. GFSSP provides pressure, temperature, and species concentrations at designated locations, or nodes, within a fluid network based on a finite volume formulation of thermodynamics and conservation laws. This paper describes the theoretical basis for the construction of the models, their implementation in the current GFSSP modeling system, and a brief evaluation of the usefulness of the model results, as well as their applicability toward a broader spectrum of analytical problems in both university teaching and engineering research.

Smith, Amanda D.↗

Global benefits of non-continuous flooding to reduce greenhouse gases and irrigation water use without rice yield penalty

Non-continuous flooding is an effective practice for reducing greenhouse gas emissions (GHGs) and irrigation water use (IRR) in rice fields. However, advancing global implementation is hampered by the lack of comprehensive understanding of GHGs and IRR reduction benefits without compromising rice yield. Here, we present the largest observational data set for such effects as of yet. By using Random Forest regression models based on 636 field trials at 105 globally georeferenced sites, we identified the key drivers of effects of non-continuous flooding practices and mapped maximum GHGs or IRR reduction benefits under optimal non-continuous flooding strategies. The results show that variation in effects of non-continuous flooding practices are primarily explained by the UnFlooded days Ratio (UFR, that is the ratio of the number of days without standing water in the field to total days of the growing period). Non-continuous flooding practices could be feasible to be adopted in 76% of global rice harvested areas. This would reduce the global warming potential (GWP) of CH4 and N2O combined from rice production by 47% or the total GWP by 7% and alleviate irrigation water use by 25%, while maintaining yield levels. The identified UFR targets far exceed currently observed levels particularly in South and Southeast Asia, suggesting large opportunities for climate mitigation and water use conservation, associated with the rigorous implementation of non-continuous flooding practices in global rice cultivation.

climate change mitigation↗

Development of a Performance Portable Non-Equilibrium Plasma Fluid Solver on Adaptive Grids

This presentation will describe the numerical techniques, programming paradigms, verification, and performance of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures. Our plasma fluid model solves the conservation equations for self-consistent electrostatic Poisson, electron and heavy species transport, and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive mesh management library, AMReX (Zhang et al., JOSS, 4 (37) 1370, 2019), and can be built and run on widely available vendor specific GPU architectures (NVIDIA/AMD/Intel). We utilize a non-subcycled second order semi-implicit time-stepping method where all adaptive mesh refinement (AMR) levels are advanced with the same time step. The composite multi-level multigrid solver from within AMReX is used for each of the governing equations that are cast into a Helmholtz equation form. We have also developed a python based chemical mechanism parser framework that uses a similar format as CANTERA (Goodwin et al., Zenodo, 2018) yaml files as input. Our custom parser reads the yaml file and provides C++ files with transport and production rate functions that can be executed on both host (CPU) and device (GPU). We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on low-pressure capacitive and high-pressure streamer discharges. Our initial performance studies indicate 10X speed-up using 20 NVIDIA GPUs versus 200 CPUs for an atmospheric streamer discharge problem solved on a 512 x 1024 x 512 grid.

graphics processing units↗

Modeling the Effect of Damage on the Electrical Resistivity of Melt-Infiltrated SiC/SiC Composites

Electrical resistivity (ER) measurements are a possible health monitoring technique for CMC components in future aerospace applications. In order to use ER measurements to detect and identify damage, it is necessary to understand how each specific damage state will affect the ER response. In this paper, finite element models are developed and applied to quantify the effect of specific damage states on the electrical resistivity response in a melt-infiltrated silicon carbide (SiC) fiber-reinforced silicon carbide composite. The electrical resistivity of several damage states are calculated by simulating the electric current flow through the damaged micro-structure. This is achieved by performing the numerical solution of the steady-state conservation of charge density equation. Numerical results reveal that cracking of the inter-tow matrix has the most profound effect on the composite electrical resistivity. Also, fiber/matrix debonding at matrix cracks in the 0° tows (tows aligned with the loading direction) may cause a significant increase in the electrical resistivity, but only if the fiber coating resistivity is 1000 Ω-cm or less. Cracks in the 90° tows and the crack opening displacement have very little effect on the composite electrical resistivity.

finite element analysis↗

Conserving approximations for strongly correlated electron systems - Bethe-Salpeter equation and dynamics for the two-dimensional Hubbard model

A semianalytical approach is described for strongly correlated electronic systems which satisfies microscopic conservation laws, treats strong frequency and momentum dependences, and provides information on both static and dynamic properties. This approach may be used to treat large systems and temperatures lower than those currently accessible to finite-temperature quantum Monte Carlo techniques. Examples of such systems include heavy-electron compounds, organic Bechegaard salts, bis-(ethylenedithiolo)-TTF superconductors, and the oxide superconductors. The technique is based on the derivation and self-consistent solution of infinite-order conserving approximations. The technique is used to derive a low-temperature phase diagram and dynamic correlation functions for the two-dimensional Hubbard lattice model.

Bickers, N. E.↗