SEARCH · Engineering Papers

Results for “flow 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,279 records · Page 71

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Validation of EDQNM for subgrid and supergrid models

From preliminary calculations it was concluded that the two-point Eddy Damping Quasi-Normal Markovian (EDQNM) closure accurately describes the behavior of second order moments. This closure can be applied as subgrid and supergrid models for Large Eddy Simulations at higher Reynolds numbers. In the case of homogeneous anisotropic turbulence, when the nonlinear terms are introduced the calculation becomes quite onerous but is still considerably less expensive than the calculation of a DS. The major merit of two-point closure models is that they can be easily applied to flows at Reynolds numbers that are unreachable by a DS. Work is in progress to derive expressions for the nonlinear terms that give good global conservation properties.

Orlandi, P.↗

NASA NTRS1990

Application of a streamwise upwind algorithm for unsteady transonic computations over oscillating wings

A new streamwise upwind algorithm has been derived to compute unsteady flows with a moving grid system and applied to compute flows over oscillating wings at transonic Mach numbers. Comparisons have been made between results obtained from this upwind algorithm, using both temporally nonconservative- and conservative-implicit methods, with the results obtained from a central-difference method, and also with experimental data. The results show (1) the efficiency and practicality of the temporally nonconservative implicit solver and (2) the robustness and accuracy of the upwind method for unsteady computations compared to the central-difference method.

Obayashi, Shigeru↗

NASA NTRS1994

Lox droplet vaporization in a supercritical forced convective environment

A systematic investigation has been conducted to study the effects of ambient flow conditions (i.e. pressure and velocity) on supercritical droplet gasification in a forced-convective environment. The model is based on the time-dependent conservation equations in axisymmetric coordinates, and accommodates thermodynamic nonidealities and transport anomalies. In addition, an efficient scheme for evaluating thermophysical properties over the entire range of fluid thermodynamic states is established. The analysis allows a thorough examination of droplet behavior during its entire lifetime, including transient gasification, dynamic deformation, and shattering. A parametric study of droplet vaporization rate in terms of ambient pressure and Reynolds number is also conducted.

Hsiao, Chia-Chun↗

NASA NTRS1995

Useful Scaling Parameters for the Pulse Tube

A set of eight non-dimensional scaling parameters for use in evaluating the performance of Pulse Tube Refrigerators is presented. The parameters result after scaling the mass, momentum and energy conservation equations for an axisymmetric, two-dimensional system. The physical interpretation of the parameters are described, and their usefulness is outlined for the enthalpy flow tube (open tube of the pulse tube). The scaling parameters allow the experimentalist to characterize three types of transport: enthalpy flow, mass streaming and heat transfer between the gas and the tube. Also reported are the results from a flow visualization experiment in which steady mass streaming in compressible oscillating flow is observed.

Lee, J. M.↗

NASA NTRSDate not supplied

Applying Machine Learning Tools for Runway Configuration Decision Support

Determining optimal runway configurations at airports, a responsibility assigned to air traffic controllers, is a challenging task. The decision-making process is intricate and involves consideration of many factors such as prevailing wind condition, convective weather, visibility, cloud ceilings, departure and arrival demand, traffic flow, equipment status, and other airport constraints. In a previous work, we developed a Runway Configuration Assistance tool using an offline reinforcement learning method called conservative Q-learning. In this paper, we evaluate and validate our Runway Configuration Assistance tool as a decision support for air traffic controllers. We validated our tool using three airports with differing levels of complexity: Charlotte Douglas International Airport, Denver International Airport, and Dallas Fort Worth International Airport. We quantified the performance of the Runway Configuration Assistance tool based on (1) agreement with historical air traffic controller decisions and (2) violation of decisions that would be obvious to subject-matter experts. Our tool showed promising results in both performance metrics for the three airports, despite the complexities in the runway configuration decision-making process. We also discuss challenges in using machine learning in general to aid air traffic management and identify deployment considerations for the Runway Configuration Assistance tool.

Runway Configuration Management↗

NASA NTRSDate not supplied

Machine Learning Tools for Runway Configuration Decision Support

Runway Configuration Management↗

NASA NTRS1997

Predictive Modeling in Plasma Reactor and Process Design

Research continues toward the improvement and increased understanding of high-density plasma tools. Such reactor systems are lauded for their independent control of ion flux and energy enabling high etch rates with low ion damage and for their improved ion velocity anisotropy resulting from thin collisionless sheaths and low neutral pressures. Still, with the transition to 300 mm processing, achieving etch uniformity and high etch rates concurrently may be a formidable task for such large diameter wafers for which computational modeling can play an important role in successful reactor and process design. The inductively coupled plasma (ICP) reactor is the focus of the present investigation. The present work attempts to understand the fundamental physical phenomena of such systems through computational modeling. Simulations will be presented using both computational fluid dynamics (CFD) techniques and the direct simulation Monte Carlo (DSMC) method for argon and chlorine discharges. ICP reactors generally operate at pressures on the order of 1 to 10 mTorr. At such low pressures, rarefaction can be significant to the degree that the constitutive relations used in typical CFD techniques become invalid and a particle simulation must be employed. This work will assess the extent to which CFD can be applied and evaluate the degree to which accuracy is lost in prediction of the phenomenon of interest; i.e., etch rate. If the CFD approach is found reasonably accurate and bench-marked with DSMC and experimental results, it has the potential to serve as a design tool due to the rapid time relative to DSMC. The continuum CFD simulation solves the governing equations for plasma flow using a finite difference technique with an implicit Gauss-Seidel Line Relaxation method for time marching toward a converged solution. The equation set consists of mass conservation for each species, separate energy equations for the electrons and heavy species, and momentum equations for the gas. The sheath is modeled by imposing the Bohm velocity to the ions near the walls. The DSMC method simulates each constituent of the gas as a separate species which would be analogous in CFD to employing separate species mass, momentum, and energy equations. All particles including electrons are moved and allowed to collide with one another with the stipulation that the electrons remain tied to the ions consistent with the concept of ambipolar diffusion. The velocities of the electrons are allowed to be modified during collisions and are not confined to a Maxwellian distribution. These benefits come at a price in terms of computational time and memory. The DSMC and CFD are made as consistent as possible by using similar chemistry and power deposition models. Although the comparison of CFD and DSMC is interesting, the main goal of this work is the increased understanding of high-density plasma flowfields that can then direct improvements in both techniques. This work is unique in the level of the physical models employed in both the DSMC and CFD for high-density plasma reactor applications. For example, the electrons are simulated in the present DSMC work which has not been done before for low temperature plasma processing problems. In the CFD approach, for the first time, the charged particle transport (discharge physics) has been self-consistently coupled to the gas flow and heat transfer.

Hash, D. B.↗

NASA NTRS1976

Computation of the inviscid supersonic flow over an external axial corner

A second-order finite-difference procedure is used to evaluate the inviscid supersonic flowfield surrounding an external axial corner composed of swept planar compression surfaces and representing the inlets on existing high-speed aircraft. The governing partial differential equations in conservation-law form are hyperbolic with respect to the axial coordinate and are solved iteratively by means of MacCormack's algorithm. The procedure treats both the peripheral shock wave and vortical singularities as discontinuities. Numerical results are presented for two parametric studies regarding the effects on the flowfield of varying the free-stream Mach number and the leading edge sweep of the horizontal wedge. Results of parametric Mach number study agree with the Mach number independence principle in that as the Mach number increases, such characteristics as shock shape, cross-flow sonic line location, and vortical singularity position approach an asymptote.

Kutler, P.↗

NASA NTRS1983

Analysis of aerothermal loads on spherical dome protuberances

Hypersonic flow over spherical dome protuberances was investigated to determine increased pressure and heating loads to the surface. The configuration was mathematically modeled in a time-dependant three-dimensional analysis of the conservation of mass, momentum (Navier-Stokes), and energy equations. A boundary mapping techique was used to obtain a rectangular parallelepiped computational domain, and a MacCormack explicit time-split predictor-corrector finite difference algorithm was used to obtain solutions. Results show local pressures and heating rates for domes one-half, one, and two boundary layer thicknesses high were increased by factors on the order of 1.4, 2, and 6, respectively. However, because lee-side pressure and thermal loads were reduced the two lower height domes did not experience any net increase in total loads. The total loads on the higher dome were increased by twenty-five percent. Flow over the lower dome was everywhere attached while flow over the intermediate dome had small windward and leeside separations. The higher dome had an unsteady windward separation region and a large leeside separation region. Trailing vortices form on all domes with intensity increasing with dome height. Discussions of applying the results to a thermally bowed thermal protection system are presented.

Olsen, G. C.↗

DOE OSTI2018

Ambient and building condition effects modeling of air-cooled natural circulation systems

Throughout the testing program at Natural convection Shutdown heat removal Test Facility (NSTF), strong influences of ambient conditions were observed in the experimental data when baseline tests were repeated under identical test procedures. Significant analysis efforts were conducted to gain understanding of these influences, and quantify the system response of the test facility. In this paper, an empirical model was developed based on theoretical considerations and using experimental data to correlate zero-power system flow rates with ambient meteorological conditions, including wind and temperatures. Coefficients in the model were determined based on best fitting of the experimental data set. The predictive capability of the empirical model was demonstrated by applying it to new sets of experimental data. Finally, after implementing the empirical model in both system-level thermal fluids and computational fluid dynamics modeling of the facility, simulations from both approaches resulted in accurate predictions in system natural circulation flow rates.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2016

Upwind methods for the Baer–Nunziato equations and higher-order reconstruction using artificial viscosity

This article is devoted to the numerical discretisation of the hyperbolic two-phase flow model of Baer and Nunziato. A special attention is paid on the discretisation of intercell flux functions in the framework of Finite Volume and Discontinuous Galerkin approaches, where care has to be taken to efficiently approximate the non-conservative products inherent to the model equations. Various upwind approximate Riemann solvers have been tested on a bench of discontinuous test cases. New discretisation schemes are proposed in a Discontinuous Galerkin framework following the criterion of Abgrall and the path-conservative formalism. A stabilisation technique based on artificial viscosity is applied to the high-order Discontinuous Galerkin method and compared against classical TVD-MUSCL Finite Volume flux reconstruction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1992

Estimation of excitation energy of diatomic molecules in expanding nonequilibrium flows

The energy contained in the highly excited vibrational and rotational states in a diatomic gas in a thermochemical nonequilibrium state during expansion is estimated. The estimation is made on the assumption that the populations of the vibrational and rotational states, when normalized by their respective equilibrium values, are describable by simple functions containing no more than four arbitrary parameters. A cubic polynomial, a logarithmic-cubic polynomial, and a bimodal step function are used for this purpose. The four parameters are determined by imposing conditions known at the ground state and the dissociation limit and the mass conservation law. The energy in excess of that accounted for by assuming a Boltzmann distribution of these states, defined here as excess excitation energy, is calculated for N2, O2, NO, CO, OH, and H2. A calculation made for a typical nozzle flow shows that the excess energy may reach 6 percent of the total enthalpy of the flow, and that the flow velocity may decrease by as much as 4 percent due to the nonequilibrium excitation phenomenon.

Park, Chul↗

DOE OSTI2017

KIVA-hpFE: Predictive turublent reactive and multiphase flow in engines : Science Supporting Mission of the Laboratory

Research and development of KIVA-hpFE for turbulent reactive and multiphase flow particularly as related to engine modeling program has relevance to National energy security and climate change. Climate change is a source problem, and energy national security is consumption of petroleum products problem. Accurately predicting engine processes leads to, lower greenhouse gas (GHG) emission, where engines in the transportation sector currently account for 26% of the U.S. GHG emissions. Less dependence on petroleum products leads to greater energy security. By Environmental Protection Agency standards, some vehicles are now reaching 42 to the 50 mpg mark. These are conventional gasoline engines. Continued investment and research into new technical innovations, the potential exists to save more than 4 million barrels of oil per day or approximately $200 to $400 million per day. This would be a significant decrease in emission and use of petroleum and a very large economic stimulus too! It is estimated with further advancements in combustion, the current emissions can be reduced up to 40%. Enabling better understanding of fuel injection and fuel-air mixing, thermodynamic combustion losses, and combustion/emission formation processes enhances our ability to help solve both problems. To provide adequate capability for accurately simulating these processes, minimize time and labor for development of engine technology, are the goals of our KIVA development program.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2018

A Statistical Approach to Safety Margin Evaluation of Irradiation Experiments in the ATR

The current safety basis of the INL Advanced Test Reactor requires that irradiation experiments meet minimum safety margins to critical heat flux and flow instability during anticipated reactor fault conditions. Compliance to these requirements may be demonstrated by calculating safety margins using safety factors that account for uncertainties, or using a statistical approach to quantify the uncertainty in the safety margins. Safety margins obtained from the first approach are conservative since the individual events leading to the cumulative safety factor are assumed to occur simultaneously at their maximum uncertainty. The statistical approach is less conservative in the case where the uncertainties are not too large, and may reduce the conservatism resulting from an overabundance of safety factors. This report describes the procedure for safety margin evaluation using a statistical approach and identifies the important thermal-hydraulic variables to consider in an uncertainty analysis. In particular, safety margins to critical heat flux and flow instability are evaluated for three recent high power fueled experiments (MP-1, FSP-1 and AFIP-6) using RELAP5 coupled to RAVEN to perform statistical analysis and uncertainty quantification. The results demonstrate adequate safety margins to critical heat flux and flow instability during an anticipated loss of commercial power to the reactor.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2014

A projection gradient method for computing ground state of spin-2 Bose–Einstein condensates

In this paper, a projection gradient method is presented for computing ground state of spin-2 Bose–Einstein condensates (BEC). We first propose the general projection gradient method for solving energy functional minimization problem under multiple constraints, in which the energy functional takes real functions as independent variables. We next extend the method to solve a similar problem, where the energy functional now takes complex functions as independent variables. We finally employ the method into finding the ground state of spin-2 BEC. The key of our method is: by constructing continuous gradient flows (CGFs), the ground state of spin-2 BEC can be computed as the steady state solution of such CGFs. We discretized the CGFs by a conservative finite difference method along with a proper way to deal with the nonlinear terms. We show that the numerical discretization is normalization and magnetization conservative and energy diminishing. Numerical results of the ground state and their energy of spin-2 BEC are reported to demonstrate the effectiveness of the numerical method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Invariant regimes of Spencer scaling law for magnetic compression of rotating FRC plasma

Abstract The scaling laws for the magnetic compression of a toroidally rotating field reversed configuration (FRC) have been investigated in this work. The magnetohydrodynamics (MHD) simulations of the magnetic compression on rotating FRCs employing the NIMROD code (Sovinec et al 2004 J. Comput. Phys. 195 355), are compared with the Spencer’s one-dimensional (1D) theory (Spencer et al 1983 Phys. Fluids 26 1564) for a wide range of initial flow speeds and profiles. The toroidal flow can influence the scalings directly through the alteration of the compressional work as also evidenced in the 1D adiabatic model, and indirectly by reshaping the initial equilibrium. However, in comparison to the static initial FRC equilibrium cases, the pressure and the radius scalings remain invariant for the magnetic compression ratio B w 2 / B w 1 up to 6 in presence of the initial equilibrium flow, suggesting a broader applicable regime of the Spencer scaling law for FRC magnetic compression. The invariant scaling has been proven a natural consequence of the conservation of angular momentum of both fluid and magnetic field during the dynamic compression process.

Ma, Yiming↗

DOE OSTI2020

Thermal Image Processing for Feature Extraction from Encapsulated Phase Change Materials

Encapsulated inorganic particles with high melting points (>300 °C) are desired as high-temperature Phase Change Materials (PCMs) for next-generation Latent Heat Thermal Energy Storage (LHTES) systems. One of the many challenges during the development of PCMs is to achieve a high throughput that in turn depends on accurately modeling the relation between process parameters and geometric & thermal properties of the PCMs particle. During the production of the PCMs, a high-speed infrared camera is used to acquire images of the encapsulated material under controlled illumination conditions. This research article focuses on the development of image processing techniques for both geometric and thermal feature extraction during the development of the PCMs. A user-friendly GUI has been designed in MATLAB and preliminary experimental results have demonstrated that the method is fast, accurate and reliable for a high throughput production. The extracted features will be used to develop Machine Learning (ML) models to predict the geometric and thermal properties of the PCM based on the process parameter settings. The ML model will accelerate the search for the optimized process settings to boost the throughput of the production.

25 ENERGY STORAGE↗

DOE OSTI2023

An Analytic Formula for Entraining CAPE in Midlatitude Storm Environments

Abstract This article introduces an analytic formula for entraining convective available potential energy (ECAPE) with an entrainment rate that is determined directly from an environmental sounding, rather than prescribed by the formula user. Entrainment is connected to the background environment using an eddy diffusivity approximation for lateral mixing, updraft geometry assumptions, and mass continuity. These approximations result in a direct correspondence between the storm-relative flow and the updraft radius and an inverse scaling between the updraft radius squared and entrainment rate. The aforementioned concepts, combined with the assumption of adiabatic conservation of moist static energy, yield an explicit analytic equation for ECAPE that depends entirely on state variables in an atmospheric profile and a few constant parameters with values that are established in past literature. Using a simplified Bernoulli-like equation, the ECAPE formula is modified to account for updraft enhancement via kinetic energy extracted from the cloud’s background environment. CAPE and ECAPE can be viewed as predictors of the maximum vertical velocity w max in an updraft. Hence, these formulas are evaluated using w max from past numerical modeling studies. Both of the new formulas improve predictions of w max substantially over commonly used diagnostic parameters, including undiluted CAPE and ECAPE with a constant prescribed entrainment rate. The formula that incorporates environmental kinetic energy contribution to the updraft correctly predicts instances of exceedance of by w max , and provides a conceptual explanation for why such exceedance is rare among past simulations. These formulas are potentially useful in nowcasting and forecasting thunderstorms and as thunderstorm proxies in climate change studies. Significance Statement Substantial mixing occurs between the upward-moving air currents in thunderstorms (updrafts) and the surrounding comparatively dry environmental air, through a process called entrainment. Entrainment controls thunderstorm intensity via its diluting effect on the buoyancy of air within updrafts. A challenge to representing entrainment in forecasting and predictions of the intensity of updrafts in future climates is to determine how much entrainment will occur in a given thunderstorm environment without a computationally expensive high-resolution simulation. To address this gap, this article derives a new formula that computes entrainment from the properties of a single environmental profile. This formula is shown to predict updraft vertical velocity more accurately than past diagnostics, and can be used in forecasting and climate prediction to improve predictions of thunderstorm behavior and impacts.

54 ENVIRONMENTAL SCIENCES↗