Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure model”

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 1,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Newly Available Single Field-of-view Sounder Atmospheric Product (SiFSAP) and Its Derivative Product

The Single Field-of-view Sounder Atmospheric product (SiFSAP) has been developed and delivered to NASA GES DISC. The SiFSAP Algorithm Theoretical Basis Documents (ATBD) and users manuals are ready to be released to public. This novel data product supplements other operational products such as AIRS version 7 and the Community Long-term Infrared Microwave Combined Atmospheric Product System (CLIMCAPS) from two main perspectives: 1) improving the spatial resolution of sounder Level-2 data to extend its usage in weather and dynamics focus area; 2) establishing radiance closure between Level-2 data and directly measured Level-1 radiances to facilitate the climate trend analysis. SiFSAP has 3-times higher spatial resolution and 9-times denser data products comparing to current NASA and NOAA operational IR sounder products. SiFSAP is derived using the optimal estimation method based physical retrieval algorithm. The Principal Component-based Radiative Transfer Model (PCRTM) which includes the cloud scattering simulation is used for the forward model so that the solution can fit the spectral radiances under all-sky conditions for individual single field-of-view (SFOV) measurements. A general introduction of the SiFSAP algorithm and corresponding validation work will be presented. Also introduced here is the Climate Fingerprinting Sounder product (ClimFiSP) that is the derivative product of SiFSAP and will be released in the near future. The ClimFiSP algorithm uses pre-constructed fingerprinting relationship to achieve a low-latency Level-3 data production and facilitate the fusion of data of different sounders.

Wan Wu↗
DOE OSTI2025

Direct statistical simulation of the Lorenz96 system in model reduction approaches

Direct statistical simulation (DSS) of nonlinear dynamical systems bypasses the traditional route of accumulating statistics by lengthy direct numerical simulations by solving the equations that govern the statistics themselves. DSS suffers, however, from the curse of dimensionality as the statistics (such as correlations) generally have higher dimensions than the underlying dynamical variables. Here we investigate two approaches to reduce the dimensionality of DSS, illustrating each method with numerical experiments with the Lorenz96 dynamical system. The forms of DSS chosen here involve approximate closures at second and third order in the equal-time cumulants. We demonstrate significant reduction in computational effort that can be achieved without sacrificing the accuracy of DSS. The methods developed here can be applied to turbulent fluid and magnetohydrodynamical systems. Published by the American Physical Society 2025

Li, Kuan↗
DOE OSTI2023

Predicting future well performance for environmental remediation design using deep learning

Here in this study, we developed a deep learning (DL) framework with a multi-channel three-dimensional convolutional neural network (MC3D-CNN) to predict well performance and thereby assist future environmental remediation design. Such prediction of extraction well performance at designated locations is critical for configuring pump-and-treat (P&T) well network design and operation, setting reasonable target closure dates for overall remedying, and estimating remedy costs. The framework is developed with operational and monitoring data routinely collected during P&T remedy operations, including well extraction and injection rates as well as in situ contaminant concentrations. Traditionally, the collected data were rarely used for purposes other than assessing past well performance and the accuracy of the conceptual site model. However, recent advances in data-driven computational approaches enable better use of the large datasets to inform future well performance, enhance site characterization, and improve remediation planning. In this study, we established a DL framework to integrate transient three-dimensional contaminant plumes and multiple aquifer properties (e.g., hydraulic conductivity and hydrostratigraphic maps) to identify characteristic patterns controlling and representing extraction well mass recovery, aiming at providing future mass recovery estimates for existing wells and candidate wells at any proposed locations. We evaluated our framework by using a realistic synthetic dataset generated from a well-calibrated flow and transport model used in the 200 West Area of the U.S. Department of Energy’s Hanford Site in southeastern Washington state. The multi-channel feature in our framework allows integration of various types and temporal densities of training datasets for DL model development. Overall, we found that the trained DL model achieved an accuracy of over 90% in ranking extraction well performance in validation datasets, and over 80% in predicting high-performance-ranking well locations. This data-informed approach provides a flexible tool to support adaptive site management, streamline decision-making, and potentially reduce remediation time and costs. Our DL framework can be used as a filtering tool to improve the current P&T network optimization design by reducing the number of candidate well locations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Parallel expansion of a fuel pellet plasmoid

The problem of the assimilation of a cryogenic fuel pellet injected into a hot plasma is considered. Due to the transparency to ambient particles of the plasmoid, the localised region of high-density plasma created by ionisation of the ablated pellet material, electrons reach a ‘quasiequilibrium’ (QE) state which is characterised by a steady-state on the fastest collisional time scale. The simplified electron kinetic equation of the QE state is solved. Taking a velocity moment of the higher-order electron kinetic equation, which is valid on the expansion time scale, permits a fluid closure, yielding an evolution equation for the macroscopic parameters describing the QE distribution function. In contrast to the Braginskii equations, the closure does not require that electrons have a short mean free path compared with the size of density perturbations, and permits an anisotropic and highly non-Maxwellian distribution function. As the QE distribution function accounts for both trapped and passing electrons, the self-consistent electric potential that causes the expansion can be properly described, in contrast to earlier models of pellet plasmoid expansion with an unbounded potential. The plasmoid expansion is simulated using both a Vlasov model and a cold-fluid model for the ions. During the expansion plasmoid ions and electrons obtain nearly equal amounts of energy; as hot ambient electrons provide this energy in the form of collisional heating of plasmoid electrons, the expansion of a pellet plasmoid is expected to be a potent mechanism for the transfer of energy from electrons to ions on a time scale shorter than that of ion–electron thermalisation.

Physics↗
NASA NTRS1986

Catalytic processes for space station waste conversion

Catalytic techniques for processing waste products onboard space vehicles were evaluated. The goal of the study was the conversion of waste to carbon, wash water, oxygen and nitrogen. However, the ultimate goal is conversion to plant nutrients and other materials useful in closure of an ecological life support system for extended planetary missions. The resulting process studied involves hydrolysis at 250 C and 600 psia to break down and compact cellulose material, distillation at 100 C to remove water, coking at 450 C and atmospheric pressure, and catalytic oxidation at 450 to 600 C and atmospheric pressure. Tests were conducted with a model waste to characterize the hydrolysis and coking processes. An oxidizer reactor was sized based on automotive catalytic conversion experience. Products obtained from the hydrolysis and coking steps included a solid residue, gases, water condensate streams, and a volatile coker oil. Based on the data obtained, sufficient component sizing was performed to make a preliminary comparison of the catalytic technique with oxidation for processing waste for a six-man spacecraft. Wet oxidation seems to be the preferred technique from the standpoint of both component simplicity and power consumption.

Schoonover, M. W.↗
NASA NTRS1994

The whistler heat flux instability: Threshold conditions in the solar wind

Solar wind electrons are observed often to consist of two components: a core and a halo. The anisotropics and relative average speeds of these components correspond to a heat flux that has the potential to excite several different electromagetic instabilities; wave-particle scattering by the resulting enhanced fluctuations can limit this heat flux. This manuscript describes theoretical studies using the linear Vlasco dispersion equation for drifting bi-Maxwellian component distributions in a homogeneous plasma to examine the threshold of the whistler heat flux instability. Expressions for this threshold are obtained from two different parametric baselines: a local model that yields scalings as functions of local dimensionless plasma paramaters, and a global model based on average electron properties observed during the in-eliptic phase of the Ulysses mission. The latter model yields an expression for the heat flux at threshold of the whistler instability as a function of helisopheric radius that scales in the same way as the average heat flux observed form Ulysses and that provides an approximate upper bound for that same quantity. This theoretical scaling is combined with the observational results to yield a semi-empirical closure relation for the average electron heat flux in the solar wind between 1 and 5 AU.

Gary, S. Peter↗
NASA NTRS2011

Constraint Embedding Technique for Multibody System Dynamics

Multibody dynamics play a critical role in simulation testbeds for space missions. There has been a considerable interest in the development of efficient computational algorithms for solving the dynamics of multibody systems. Mass matrix factorization and inversion techniques and the O(N) class of forward dynamics algorithms developed using a spatial operator algebra stand out as important breakthrough on this front. Techniques such as these provide the efficient algorithms and methods for the application and implementation of such multibody dynamics models. However, these methods are limited only to tree-topology multibody systems. Closed-chain topology systems require different techniques that are not as efficient or as broad as those for tree-topology systems. The closed-chain forward dynamics approach consists of treating the closed-chain topology as a tree-topology system subject to additional closure constraints. The resulting forward dynamics solution consists of: (a) ignoring the closure constraints and using the O(N) algorithm to solve for the free unconstrained accelerations for the system; (b) using the tree-topology solution to compute a correction force to enforce the closure constraints; and (c) correcting the unconstrained accelerations with correction accelerations resulting from the correction forces. This constraint-embedding technique shows how to use direct embedding to eliminate local closure-loops in the system and effectively convert the system back to a tree-topology system. At this point, standard tree-topology techniques can be brought to bear on the problem. The approach uses a spatial operator algebra approach to formulating the equations of motion. The operators are block-partitioned around the local body subgroups to convert them into aggregate bodies. Mass matrix operator factorization and inversion techniques are applied to the reformulated tree-topology system. Thus in essence, the new technique allows conversion of a system with closure-constraints into an equivalent tree-topology system, and thus allows one to take advantage of the host of techniques available to the latter class of systems. This technology is highly suitable for the class of multibody systems where the closure-constraints are local, i.e., where they are confined to small groupings of bodies within the system. Important examples of such local closure-constraints are constraints associated with four-bar linkages, geared motors, differential suspensions, etc. One can eliminate these closure-constraints and convert the system into a tree-topology system by embedding the constraints directly into the system dynamics and effectively replacing the body groupings with virtual aggregate bodies. Once eliminated, one can apply the well-known results and algorithms for tree-topology systems to solve the dynamics of such closed-chain system.

Woo, Simon S.↗
NASA NTRS1994

Separation of crack extension modes in composite delamination problems

This work concerns fracture mechanics modeling of composite delamination problems. In order to predict delamination resistance, an applied stress intensity factor, K, or energy release rate, G, must be compared to a mode-dependent critical value of K or G from experiment. In the interfacial fracture analysis of most applications and some tests, the mode of crack extension is not uniquely defined. It is instead a function of distance from the crack tip due to the oscillating singularity existing at the tip. In this work, a consistent method is presented of extracting crack extension modes in such cases. In particular, use of the virtual crack closure technique (VCCT) to extract modes of crack extension is studied for cases of a crack along the interface between two in-plane orthotropic materials. Modes of crack extension extracted from oscillatory analyses using VCCT are a function of the virtual crack extension length, delta. Most existing efforts to obtain delta-independent modes of crack extension involve changing the analysis in order to eliminate its oscillatory nature. One such method involves changing one or more properties of the layers to make the oscillatory exponent parameter, epsilon, equal zero. Standardized application of this method would require consistent criteria for identifying which properties can be altered without changing the physical aspects of the problem. Another method involves inserting a thin homogeneous layer (typically referred to as a resin interlayer) along the interface and placing the crack within it. The drawbacks of this method are that it requires increased modeling effort and introduces the thickness of the interlayer as an additional length parameter. The approach presented here does not attempt to alter the interfacial fracture analysis to eliminate its oscillatory behavior. Instead, the argument is made that the oscillatory behavior is non-physical and that if its effects were separated from VCCT quantities, then consistent, delta-independent modes of crack extension could be defined. Knowledge of the near-tip fields in a planar orthotropic material interfacial fracture analysis is used to determine the explicit delta dependence of VCCT parameters. Once this delta dependence is determined, energy release rates are defined with this delta dependence factored out. This modified VCCT method is applied to results from two finite element test cases. It is shown that, as predicted, delta-independent modes of crack extension result. The modified VCCT approach shows potential as a consistent method of extracting crack extension modes. It uses the same information from a finite element analysis (i.e., nodal forces and displacements) as the traditional VCCT method does. The A-independent modes extracted using the modified VCCT approach can also be used as guides to test the convergence of finite element solutions.

Beuth, Jack L.↗
DOE OSTI2025

Predictive Assessment of the Chemical Composition of Coal Ash in Reserve at U.S. Disposal Sites

In the United States, more than 2 Gt of coal combustion residuals (i.e., coal ash) are stored in hundreds of disposal units. Recent federal regulations mandate the closure or retrofitting of most coal ash impoundments, presenting significant challenges for waste management. These regulatory pressures also present opportunities to reuse coal ash. However, the quality and quantity of discarded coal ash across the U.S. are not well known, even though this information is crucial for spurring its reuse for conventional and new material applications. This study describes a predictive model for the major element composition of coal ash in reserve at disposal sites of major U.S. coal-fired power plants. This model was constructed from coal purchase records of 705 power stations from 1973 to 2022 and was trained on coal ash composition data, showing that coal ash elemental composition is strongly associated with the source of feedstock coal. The model showed regional shifts in the major element contents of ash produced by power plants in the last 50 years, particularly for calcium and iron (expressed as %CaO and %Fe2O3), as power stations changed their source of coal over this time frame. Our approach enables an estimation of chemical composition for ash stored in waste impoundments at individual power stations. Such information can help to delineate the regional market resource potential of supplementary cements for concrete and other material innovations that would utilize coal ash harvested from disposal sites across the U.S.

01 COAL, LIGNITE, AND PEAT↗
NASA NTRS1980

Effect of propellant deformation on ignition and combustion processes in solid propellant cracks

A comprehensive theoretical model was formulated to study the development of convective burning in a solid propellant crack which continually deforms due to burning and pressure loading. In the theoretical model, the effect of interrelated structural deformation and combustion processes was taken into account by considering (1) transient, one dimensional mass, momentum, and energy conservation equations in the gas phase; (2) a transient, one dimensional heat conduction equation in the solid phase; and (3) quasi-static deformation of the two dimensional, linear viscoelastic propellant crack caused by pressure loading. Partial closures may generate substantial local pressure peaks along the crack, implying a strong coupling between chamber pressurization, crack combustion, and propellant deformation, especially when the cracks are narrow and the chamber pressurization rates high. The maximum pressure in the crack cavity is generally higher than that in the chamber. The initial flame-spreading process is not affected by propellant deformation.

Kumar, M.↗
NASA NTRS1981

Glass shell manufacturing in space

Highly-uniform, hollow glass spheres (shells), which are used for inertial confinement fusion targets, were formed from metal-organic gel powder feedstock in a vertical furnace. As a result of the rapid pyrolysis caused by the furnace, the gel is transformed to a shell in five distinct stages: (a) surface closure of the porous gel; (b) generation of a closed-cell foam structure in the gel; (c) spheridization of the gel and further expansion of the foam; (d) coalescence of the closed-cell foam to a single-void shell; and (e) fining of the glass shell. The heat transfer from the furnace to the falling gel particle was modeled to determine the effective heating rate of the gel. The model predicts the temperature history for a particle as a function of mass, dimensions, specific heat, and absorptance as well as furnace temperature profile and thermal conductivity of the furnace gas. A model was developed that predicts the gravity-induced degradation of shell concentricity in falling molten shells as a function of shell characteristics and time.

Downs, R. L.↗
NASA NTRS1984

Experimental investigation of shock-cell noise reduction for dual-stream nozzles in simulated flight comprehensive data report. Volume 1: Test nozzles and acoustic data

Parameters which contribute to supersonic jet shock noise were investigated for the purpose of determining means to reduce such noise generation to acceptable levels. Six dual-stream test nozzles with varying flow passage and plug closure designs were evaluated under simulated flight conditions in an anechoic chamber. All nozzles had combined convergent-divergent or convergent flow passages. Acoustic behavior as a function of nozzle flow passage geometry was measured. The acoustic data consist primarily of 1/3 octave band sound pressure levels and overall sound pressure levels. Detailed schematics and geometric characteristics of the six scale model nozzle configurations and acoustic test point definitions are presented. Tabulation of aerodynamic test conditions and a computer listing of the measured acoustic data are displayed.

Yamamoto, K.↗
DOE OSTI2020

Formation of phenanthrene via H–assisted isomerization of 2–ethynylbiphenyl produced in the reaction of phenyl with phenylacetylene

Model chemistry G3(MP2,CC)//B3LYP/6–311G(d,p) calculations of the potential energy surface for the reaction of phenyl radical (C 6 H 5 ) with phenylacetylene (C 8 H 6 ) have been carried out and combined with Rice–Ramsperger–Kassel–Marcus/Master Equation calculations of temperature– and pressure–dependent rate constants. The results showed that the reaction can serve as a viable source for the formation of phenanthrene via an indirect route involving a primary reaction of phenyl addition to the ortho carbon in the ring of phenylacetylene and H elimination producing 2–ethynylbiphenyl followed by secondary H–assisted isomerization of 2–ethynylbiphenyl to phenanthrene. In the secondary reaction, the H atom adds to the α carbon of the ethynyl side chain, then a six–member ring closure takes place followed by aromatization via an H loss. The channel of H addition to the side chain of 2–ethynylbiphenyl appears to be much faster than H addition to the ortho carbon in the ethynyl–substituted ring leading back to the initial C 6 H 5 + C 8 H 6 reactants. Rate constants for the primary C 6 H 5 + C 8 H 6 ⇌2–ethynylbiphenyl (p1) + H and secondary p1 + H⇌phenanthrene (p2) + H reactions have been computed in the temperature range of 500–2500 K at pressures of 30 Torr, 1, 10, and 100 atm and fitted to modified Arrhenius expressions. Here, the suggested kinetic scheme and rate constants are proposed as a prototype for the modeling of the growth of polycyclic aromatic hydrocarbons via the phenyl addition–dehydrocyclization (PAC) mechanism involving an addition of a PAH radical to an ethynyl–substituted PAH molecule.

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

Effects of negative triangularity shaping on energetic particle driven Alfvén eigenmodes in DIII-D

Shape variations from positive to negative triangularity may lead to improved performance regimes in a tokamak with a reduction in microturbulence as well as improved thermal confinement. The motivation of this investigation is to model and compare the neutral beam driven Alfvén eigenmode (AE) activity in two DIII-D discharges with positive and negative triangularity shaping of plasma. The simulations are performed using the linear version of the FAR3d code, which solves the reduced MHD equations for thermal plasma with addition of moment equations for the energetic ion density and parallel velocity with appropriate Landau closure relations. Overall, our numerical results indicate that for similar physical parameters, the unstable AEs observed in the negative triangularity case have lower growth rates as compared to the positive triangularity regime. Our findings may be useful to analyze the influence of the reverse-D like geometry on the AE instabilities in DIII-D and may lead to better configurations for minimizing fast ion losses in a tokamak device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Transonic Navier-Stokes flow computations over wing-fuselage geometries

The steady-state transonic viscous flows over two wing-fuselage configurations are solved numerically. The 3D compressible thin-layer Navier-Stokes equations are solved by means of multistage timestepping similar to the Runge-Kutta method. The timestepping technique is applied with a multigrid scheme to make the transonic Navier-Stokes computations using the Baldwin-Lomax turbulence for closure and converging the solution on a grid of 1.52 million grid points with a C-O topology. Both generic and fighter wing-fuselage geometries are solved, and the results are compared to experimental studies. Solution convergence is found to be fast, and the fine grid generates accurate results for both configurations by means of the Baldwin-Lomax turbulence model. Results from the Johnson-King model calculations produce results that are more consistent with experimental results for the fighter configuration.

Agrawal, S.↗
NASA NTRS1992

Large eddy simulations and direct numerical simulations of high speed turbulent reacting flows

The basic objective of this research is to extend the capabilities of Large Eddy Simulations (LES) and Direct Numerical Simulations (DNS) for the computational analyses of high speed reacting flows. In the efforts related to LES, we were primarily involved with assessing the performance of the various modern methods based on the Probability Density Function (PDF) methods for providing closures for treating the subgrid fluctuation correlations of scalar quantities in reacting turbulent flows. In the work on DNS, we concentrated on understanding some of the relevant physics of compressible reacting flows by means of statistical analysis of the data generated by DNS of such flows. In the research conducted in the second year of this program, our efforts focused on the modeling of homogeneous compressible turbulent flows by PDF methods, and on DNS of non-equilibrium reacting high speed mixing layers. Some preliminary work is also in progress on PDF modeling of shear flows, and also on LES of such flows.

Givi, P.↗
DOE OSTI2022

Effect of the Grell-Freitas Deep Convection Scheme in Quasi-Uniform and Variable-Resolution Aquaplanet CAM Simulations

The Grell-Freitas deep convection scheme (GF) has been implemented in the Community Atmosphere Model (CAM) with the nonhydrostatic Model for Prediction Across Scales (MPAS) dynamical core. Aquaplanet simulations have been performed with a global quasi-uniform resolution mesh with ~120 km grid spacing and a variable resolution mesh employing a circular refined region with ~30 km grid spacing centering at the equator and ~120 km elsewhere. GF produces overall patterns of the general circulations comparable to those of other aquaplanet experiments, including those produced using the default Zhang-McFarlane convection scheme (ZM) in CAM Version 5.4. GF alleviates the overestimated occurrence of weak intensity precipitation characteristic of ZM due to GF’s different closure methods and deep-convection trigger conditions. Compared to ZM, the GF characteristics of the equatorial Kelvin waves are closer to observations. GF also better simulates precipitation at variable resolution, with no unrealistic wave distortion by the refinement. Overall, GF shows less sensitivity of convective precipitation and heating/moistening tendencies to changes in horizontal resolution compared to ZM. However, vertically-overlapped total cloud fraction simulated with GF in the tropics is more sensitive to the resolution than ZM due to the larger resolution-sensitivity of high-level cloud fraction in the tropics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Separating divertor closure effects on divertor detachment and pedestal shape in DIII-D

Comparison between an open divertor and a more-closed divertor in DIII-D demonstrates detachment at an up to 40% lower pedestal density (n e,ped ) in the closed divertor, due to a combination of decreased fueling of the pedestal and increased dissipation in the scrape off layer (SOL) in the closed divertor, both resulting from increased neutral trapping in the divertor. Predicting whether the relationship between divertor closure and detachment will hold for an opaque SOL, in which the contribution of ionizing neutrals to fueling the pedestal is lessened, requires separating out the different mechanisms contributing to the density difference at detachment. A series of experiments on DIII-D characterizes matched discharges using various divertor configurations to isolate the effects of divertor closure. These experiments show detachment at up to 25% lower n e,sep in the closed divertor than the open, supported by simulations showing increased neutral trapping, and hence increased dissipation, in the closed divertor. A difference in n e,ped /n e,sep is also seen: for matched n e,sep , the closed divertor has up to 20% lower ne,ped, consistent with modeling showing a smaller ionization fraction inside the separatrix in this case. Understanding how these pieces fit together will help in develop of predictive models of pedestal density and detached divertors compatible with a high-performance core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗