Engineering Papers⌕ Search

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.

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

At least 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Analysis of the Harrier forebody/inlet design using computational techniques

Under the support of this Cooperative Agreement, computations of transonic flow past the complex forebody/inlet configuration of the AV-8B Harrier II have been performed. The actual aircraft configuration was measured and its surface and surrounding domain were defined using computational structured grids. The thin-layer Navier-Stokes equations were used to model the flow along with the Chimera embedded multi-grid technique. A fully conservative, alternating direction implicit (ADI), approximately-factored, partially flux-split algorithm was employed to perform the computation. An existing code was altered to conform with the needs of the study, and some special engine face boundary conditions were developed. The algorithm incorporated the Chimera technique and an algebraic turbulence model in order to deal with the embedded multi-grids and viscous governing equations. Comparison with experimental data has yielded good agreement for the simplifications incorporated into the analysis. The aim of the present research was to provide a methodology for the numerical solution of complex, combined external/internal flows. This is the first time-dependent Navier-Stokes solution for a geometry in which the fuselage and inlet share a wall. The results indicate the methodology used here is a viable tool for transonic aircraft modeling.

Chow, Chuen-Yen↗
NASA NTRS1993

Modeling of Non-Spherical Droplet Dynamics

A two-dimensional time-dependent computer code based on the modified Arbitrary Lagrangian Eulerian (ALE) technique, has been developed to simulate non-spherical droplet dynamics and evaporation under convective flows at real rocket combustion chamber conditions. The equations of mass, momentum, energy and species are simultaneously solved for both liquid and gas phases with an accurate dynamic interface tracking. The jump boundary conditions across the deforming droplet surface are obtained by applying the integral forms of conservation of mass, momentum, and energy. At each time step, the interface geometry and flow properties at the droplet surface are implicitly solved by satisfying the interface boundary conditions. A Lagrangian technique was developed to track the arbitrarily moving interface between the liquid droplet and the external gas. An elliptic grid generator is adopted to dynamically reconstruct grids both inside and outside the droplet surface. This code has been used to study droplet oscillation, droplet deformation/breakup, nonspherical droplet evaporation in both low and high pressure convective flows. This presentation briefly describes the numerical algorithm for modeling of the nonspherical droplet dynamics and demonstrates the representative simulation results of nonspherical droplet evaporation at low and high pressure convective flows. Potential applications of this code to rocket combustor design and performance predictions are discussed.

Deng, Zheng-Tao↗
NASA NTRS1989

Numerical simulation of flow through a two-strut scramjet inlet

A three-dimensional, Reynolds-averaged Navier-Stokes code has been used to numerically analyze flow through a two-strut, supersonic combustion ramjet (scramjet) inlet configuration. It solves the governing equations in full conservation form using either a fully explicit or explicit-implicit method. An algebraic, two-layer eddy-viscosity model is used for turbulent flow calculations. The analysis allows inclusion of end effects that are caused by the aft placement of the cowl on the underside of the inlet. A special grid has been developed to accommodate the struts embedded in the inlet flowfield. Detailed numerical results are presented here for the two-strut configuration, and a comparison is made with the available experimental results.

Kumar, Ajay↗
DOE OSTI2024

Predicting turbulent wake flow of marine hydrokinetic turbine arrays in large-scale waterways via physics-enhanced convolutional neural networks

We present a physics-enhanced convolutional neural network (PECNN) algorithm for reconstructing the mean flow and turbulence statistics in the wake of marine hydrokinetic (MHK) turbine arrays installed in large-scale meandering rivers. The algorithm embeds the mass and momentum conservation equations into the loss function of the PECNN algorithm to improve the physical realism of the reconstructed flow fields. The PECNN is trained using large eddy simulation (LES) results of the wake flow of a single row of turbines in a virtual meandering river. Subsequently, the trained PECNN is applied to predict the wake flow of MHK turbines with arrangements and positionings different than those considered during the training process. The PECNN predictions are validated using the results of separately performed LES. The results show that the PECNN algorithm can accurately predict the wake flow of MHK turbine farms at a small fraction of the cost of LES. The PECNN can improve the accuracy by around 1% and reduce the physical constraint indices by around 50% compared to the CNN without physical constraints. This work underscores the potential of PECNN to develop reduced-order models for control co-design and optimization of MHK turbine arrays in natural riverine environments.

Mechanics↗
NASA NTRS1979

Numerical modeling of two-dimensional confined flows

A numerical model of two-dimensional confined flows is presented. The flow in the duct is partitioned into finite streams. The difference equations are then obtained by applying conservation principles directly to the individual streams. A listing of a computer code based on this approach in FORTRAN 4 language is presented. The code computes two dimensional compressible turbulent flows in ducts when the duct area along the flow is specified and the pressure gradient is unknown.

Greywall, M. S.↗
NASA NTRS1993

A coupled implicit method for chemical non-equilibrium flows at all speeds

The present time-accurate coupled-solution procedure addresses the chemical nonequilibrium Navier-Stokes equations over a wide Mach-number range uses, in conjunction with the strong conservation form of the governing equations, five unknown primitive variables. The numerical tests undertaken address steady convergent-divergent nozzle flows with air dissociation/recombination, dump combustor flows with n-pentane/air chemistry, and unsteady nonreacting cavity flows.

Shuen, Jian-Shun↗
DOE OSTI2020

Understanding the relative importance of vertical and horizontal flow in ice-wedge polygons

Abstract. Ice-wedge polygons are common Arctic landforms. The future of these landforms in a warming climate depends on the bidirectional feedback between the rate of ice-wedge degradation and changes in hydrological characteristics. This work aims to better understand the relative roles of vertical and horizontal water fluxes in the subsurface of polygonal landscapes, providing new insights and data to test and calibrate hydrological models. Field-scale investigations were conducted at an intensively instrumented location on the Barrow Environmental Observatory (BEO) near Utqiaġvik, AK, USA. Using a conservative tracer, we examined controls of microtopography and the frost table on subsurface flow and transport within a low-centered and a high-centered polygon. Bromide tracer was applied at both polygons in July 2015 and transport was monitored through two thaw seasons. Sampler arrays placed in polygon centers, rims, and troughs were used to monitor tracer concentrations. In both polygons, the tracer first infiltrated vertically until encountering the frost table and was then transported horizontally. Horizontal flow occurred in more locations and at higher velocities in the low-centered polygon than in the high-centered polygon. Preferential flow, influenced by frost table topography, was significant between polygon centers and troughs. Estimates of horizontal hydraulic conductivity were within the range of previous estimates of vertical conductivity, highlighting the importance of horizontal flow in these systems. This work forms a basis for understanding complexity of flow in polygonal landscapes.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1987

Numerical analysis of flow through scramjet engine inlets

A set of computer programs has been developed to analyze flow through supersonic combustion ramjet (scramjet) inlets. These programs solve either the two- or three-dimensional Euler/Navier-Stokes equations in full conservation form by MacCormack's explicit or explicit-implicit method. An algebraic two-layer eddy viscosity model is used for turbulent flow calculations. The programs are operational on Control Data CYBER-200 series vector-processing computer system and have been optimized to take maximum advantage of the vector processing capability of the system. Since their development, the programs have been extensively verified and used to analyze a number of very complex inlet configurations. In this paper, results are presented from two-dimensional, quasi-three-dimensional, and three-dimensional analyses of the inlet flow field to illustrate the use of the programs.

Kumar, Ajay↗
DOE OSTI2023

Development, calibration, and validation of a novel gray-box energy model for residential split air conditioners

Energy models for vapor compression refrigeration in residential air conditioners have been developed through white-box, gray-box, and black-box methods in decades. However, existing white-box and gray-box models require complicated equations with detailed geometries while black-box models require substantial experimental data. Further, this paper aims to develop and validate a simple gray-box steady-state energy model without the need for detailed geometries, which can accurately predict the cooling capacity and electrical power input based on outdoor and indoor air conditions, and supply airflow rates. First three state variables, including the evaporation and condensation temperatures, and refrigerant mass flowrate, are applied to develop the energy model and are solved by three physical equations, including the energy conservations at the evaporator and condenser, and the refrigerant volume-mass flow correlation of the compressor. Secondly, seven performance property functions related to three state variables and three physical equations are identified and calibrated by simple temperature and power measurements. Finally, field experiments are conducted on a residential air conditioner to calibrate these performance property functions and validate the developed model. The validated results reveal the model can accurately predict the cooling capacity and electrical power input, with the normalized root mean square errors of 2.3% and 0.87% respectively.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

A review of spatial Markov models for predicting pre-asymptotic and anomalous transport in porous and fractured media

Heterogeneity across a broad range of scales in geologic porous media often manifests in observations of non-Fickian or anomalous transport. While traditional anomalous transport models can successfully make predictions in certain geological systems, increasing evidence suggests that assumptions relating to independent and identically distributed increments constrain where and when they can be reliably applied. A relatively novel model, the Spatial Markov model (SMM), relaxes the assumption of independence. The SMM belongs to the family of correlated continuous time random walks and has shown promise across a wide range of transport problems relevant to natural porous media. It has been successfully used to model conservative as well as more recently reactive transport in highly complex flows ranging from pore scales to much larger scales of interest in geology and subsurface hydrology. In this review paper we summarize its original development and provide a comprehensive review of its advances and applications as well as lay out a vision for its future development.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Probing the high-density nuclear symmetry energy with the Ξ – / Ξ 0 ratio in heavy-ion collisions at s N N ≈ 3 GeV

Recent beam energy scan experiments at the BNL Relativistic Heavy Ion Collider by the STAR Collaboration found that hadronic interactions dominate the collective flow and the proton cumulant ratios are driven by baryon number conservation in a region of high baryon density in √ s NN = 3 GeV Au + Au reactions, indicating that the dense medium formed in such collisions is likely hadronic matter. Within an updated a relativistic transport model with momentum dependent isoscalar and isovector single-nucleon mean-field potentials corresponding to different symmetry energies at suprasaturation densities, the n/p, π – /π + , $K$ $^{0}_{s}$/K + , Σ – /Σ + , and Ξ – /Ξ 0 ratios are studied for central Au + Au collisions at √ s NN = 3 GeV, where the maximum central density reaches about (3.6–4.0) ρ 0 . The doubly strange Ξ – /Ξ 0 ratio is found to have the strongest sensitivity to the variation of high-density nuclear symmetry energy. Furthermore, the Ξ – /Ξ 0 ratio in relativistic heavy-ion reactions at √ s NN~3 GeV may help probe sensitively the poorly known symmetry energy of dense neutron-rich matter critically important for understanding various properties of neutron stars.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Extracting the Quantum Hall Conductance from a Single Bulk Wave Function

In this work, we propose a new formula that extracts the quantum Hall conductance from a single (2+1)⁢ D gapped wave function. The formula applies to general many-body systems that conserve particle number, and is based on the concept of modular flow, i.e., unitary dynamics generated from the entanglement structure of the wave function. The formula is shown to satisfy all formal properties of the Hall conductance: it is odd under time reversal and reflection, even under charge conjugation, universal and topologically rigid in the thermodynamic limit. Further evidence for relating the formula to the Hall conductance is obtained from conformal field theory arguments. Finally, we numerically check the formula by applying it to a noninteracting Chern band where excellent agreement is obtained.

2-dimensional systems↗
DOE OSTI2022

Co-flow injection for serial crystallography at X-ray free-electron lasers

Serial femtosecond crystallography (SFX) is a powerful technique that exploits X-ray free-electron lasers to determine the structure of macromolecules at room temperature. Despite the impressive exposition of structural details with this novel crystallographic approach, the methods currently available to introduce crystals into the path of the X-ray beam sometimes exhibit serious drawbacks. Samples requiring liquid injection of crystal slurries consume large quantities of crystals (at times up to a gram of protein per data set), may not be compatible with vacuum configurations on beamlines or provide a high background due to additional sheathing liquids present during the injection. Proposed and characterized here is the use of an immiscible inert oil phase to supplement the flow of sample in a hybrid microfluidic 3D-printed co-flow device. Co-flow generation is reported with sample and oil phases flowing in parallel, resulting in stable injection conditions for two different resin materials experimentally. A numerical model is presented that adequately predicts these flow-rate conditions. The co-flow generating devices reduce crystal clogging effects, have the potential to conserve protein crystal samples up to 95% and will allow degradation-free light-induced time-resolved SFX.

3D printing↗
DOE OSTI2021

Computing water flow through complex landscapes – Part 3: Fill–Spill–Merge: flow routing in depression hierarchies

Abstract. Depressions – inwardly draining regions – are common to many landscapes. When there is sufficient moisture, depressions take the form of lakes and wetlands; otherwise, they may be dry. Hydrological flow models used in geomorphology, hydrology, planetary science, soil and water conservation, and other fields often eliminate depressions through filling or breaching; however, this can produce unrealistic results. Models that retain depressions, on the other hand, are often undesirably expensive to run. In previous work we began to address this by developing a depression hierarchy data structure to capture the full topographic complexity of depressions in a region. Here, we extend this work by presenting the Fill–Spill–Merge algorithm that utilizes our depression hierarchy data structure to rapidly process and distribute runoff. Runoff fills depressions, which then overflow and spill into their neighbors. If both a depression and its neighbor fill, they merge. We provide a detailed explanation of the algorithm and results from two sample study areas. In these case studies, the algorithm runs 90–2600 times faster (with a reduction in compute time of 2000–63 000 times) than the commonly used Jacobi iteration and produces a more accurate output. Complete, well-commented, open-source code with 97 % test coverage is available on GitHub and Zenodo.

58 GEOSCIENCES↗
NASA NTRS1975

Numerical solutions for inviscid supersonic corner flows

Analytical solutions for inviscid supersonic corner flows are virtually nonexistent due to the complexity of the interference geometry. In view of this, numerical solutions for compressive-compressive and expansive-compressive corner flows are obtained. The governing equations are written in strong conservation-law form and are solved iteratively in nonorthogonal conical coordinates by use of a second-order, shock-capturing, finite-difference technique. The computed wave structure and surface pressure distributions are compared with high Reynolds number (Re greater than 2,500,000 ft) experimental data. The results clearly show that the wave structure in the corner is dominated by the inviscid field.

Shankar, V. S. V.↗
NASA NTRS1979

A new flow model for highly separated airfoil flows at low speeds

An analytical model for separated airfoil flows is presented which is based on experimentally observed physical phenomena. These include a free stagnation point aft of the airfoil and a standing vortex in the separated region. A computer program is described which iteratively matches the outer potential flow, the airfoil turbulent boundary layer, the separated jet entrainment, mass conservation in the separated bubble, and the rear stagnation pressure. Separation location and pressure are not specified a priori. Results are presented for surface pressure coefficient and compared with experiment for three angles of attack for a GA(W)-1, 17% thick airfoil.

Zumwalt, G. W.↗
NASA NTRS1983

A new stream function formulation for the Euler equations

A new stream function formulation is developed for the solution of Euler's equations in the transonic flow region. The stream function and the density are the dependent variables in this method, while the governing equations for adiabatic flow are the momentum equations which are solved in the strong conservation law form. The application of this method does not require a knowledge of the vorticity. The algorithm is combined with the automatic grid solver (GRAPE) of Steger and Sorenson (1979) in order to study arbitrary geometries. Results of the application of this method are presented for the NACA 0012 airfoil at various Mach numbers and angles of attack, and cylinders. In addition, detailed comparisons are made with other solutions of the Euler equations.

Atkins, H. L.↗
NASA NTRS1994

A time-dependent gyro-kinetic model of thermal ion upflows in the high-latitude F region

Ample evidence supports the significance of the high-latitude ionospheric contribution to magnetospheric plasma. Assuming flux conservation along a flux tube, the upward field-aligned ion flows observed in the magnetosphere require high-latitude ionospheric field-aligned ion upflows of the order of 10(exp 8) to 10(exp 9)/sq cm/s. Since radar and satellite observations of high-latitude F region flows at times exceed this flux requirement by an order of magnitude, the thermal ionospheric upflows are not simply the ionospheric response to a magnetospheric flux requirement. Several ionospheric ion upflow mechanisms have been proposed, but simulations based on fluid theory do not reproduce all the observed features of ionospheric ion upflows. Certain asymmetries in the statistical morphology of high-latitude F region ion upflows suggest that the ion upflows may be generated by ion-neutral frictional heating. We developed a single-component (O(+)), time-dependent gyro-kinetic model of the high-latitude F region response to frictional heating in which the neutral exobase is a discontinuous boundary between fully collisional and collisionless plasmas. The concept of a discontinuous neutreal exobase and the assumption of a constant and uniform polarization electric field reduce the ion velocity distribution function, from which we can compute the ion density, parallel velocity, parallel and perpendicular temperature, and parallel flux. Using our model, we simulated the response of a convecting flux tube between 500 km and 2500 km to various frictional heating inputs; the results were both qualitatively and quantitatively different from fluid model results, which may indicate an inadequacy of the fluid theory approach. The gyro-kinetic frictional heating model responses to the various simulations were qualitatively similar: (1) initial perturbations of all the modeled parameters propagated rapidly up the flux tube, (2) transient values of the ion parallel velocity, temperature, and flux exceeded 3 km/s, 2 x 10(exp 4) K, and 10(exp 9)/sq cm/s, respectively, (3) a second transient regime developed wherein the parallel temperature drops to very low values (a few hundred Kelvins), and (4) well after heating ceased, large parallel temperatures and large downward parallel velocities and fluxes developed as the flux tube slowly returned to diffusive equilibrium. The ion velocity distributions during the simulation are often non-Maxwellian and are sometimes composed of two distinct ion populations.

Loranc, M.↗