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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

Energy and economic trade offs for advanced technology subsonic aircraft

Changes in future aircraft technology which conserve energy are studied, along with the effect of these changes on economic performance. Among the new technologies considered are laminar-flow control, composite materials with and without laminar-flow control, and advanced airfoils. Aircraft design features studied include high-aspect-ratio wings, thickness ratio, and range. Engine technology is held constant at the JT9D level. It is concluded that wing aspect ratios of future aircraft are likely to significantly increase as a result of new technology and the push of higher fuel prices. Whereas current airplanes have been designed for AR = 7, supercritical technology and much higher fuel prices will drive aspect ratio to the AR = 9-10 range. Composite materials may raise aspect ratio to about 11-12 and practical laminar flow-control systems may further increase aspect ratio to 14 or more. Advanced technology provides significant reductions in aircraft take-off gross weight, energy consumption, and direct operating cost.

Maddalon, D. V.↗
NASA NTRS1977

A parametric investigation of an existing supersonic relative tip speed propeller noise model

A high tip speed turboprop is being considered as a future energy conservative airplane. The high tip speed of the propeller combined with the cruise speed of the airplane may result in supersonic relative flow on the propeller tips. These supersonic blade sections could generate noise that is a cabin environment problem. An existing supersonic propeller noise model was parametrically investigated to identify and evaluate the noise reduction variables. Both independent and interdependent parameter variations (constant propeller thrust) were performed. The noise reductions indicated by the independent investigation varied from sizable in the case of reducing Mach number to minimal for adjusting the thickness and loading distributions. The noise reduction possibilities of decreasing relative Mach number were further investigated during the interdependent variations. The interdependent investigation indicated that significant noise reductions could be achieved by increasing the propeller diameter and/or increasing the number of propeller blades while maintaining a constant propeller thrust.

Dittmar, J. H.↗
NASA NTRS1979

Calculated global conditions for the unsteady Euler's flow equations

This paper considers the numerical solution of a fluid dynamical system with time-dependent oscillating boundaries. The problem detailed, a compressible and inviscid flow in a one-dimensional pipe with vibrating ends, is examined to develop techniques needed to solve the more complex unsteady flows in turbomachine cascades. The problem is shown to be mathematically ill-posed but by complementing the governing equations with energy conservation, a unique solution is obtained. An artificial damping term is introduced into the governing equations to improve the performance of the time-marching integration.

Kuzo, D. M.↗
NASA NTRS1985

Measurement of acoustic shielding by a turbulent jet

The acoustic shielding properties of a turbulent jet have been investigated experimentally. The experimental arrangement consisted of an acoustic point source and a turbulent shielding jet. The source and jet parameters investigated include the source frequency, source spectrum, jet velocity, jet heating by simulation and the lateral and longitudinal source positions with respect to the shielding jet. It is found that the maximum sound attenuation provided by the shielding jet depends on the balance between refraction and diffraction. Over the frequency range investigated, the redistribution of sound by the shielding air jet is power conserving. Comparison between measurement and prediction based on an idealized cylindrical uniform jet model indicates that the spreading and decay of the flow field in real jets are important. Comparisons between the present data trends and those reported for jet-by-jet shielding suggest that the major effects observed in the latter are acoustical rather than aerodynamical.

Yu, J. C.↗
DOE OSTI2022

Preliminary MITR Design-Demonstration Element Thermal-Hydraulics and Structural Analyses

The United State High Performance Research Reactor (USHPRR) program aims to eliminate more than 200kg of High Enriched Uranium (HEU) from commerce annually by converting five U.S. high-performance research reactors and one associated critical assembly to Low Enriched Uranium (LEU) fuel using a high-density alloy of uranium-10 wt% molybdenum (U-10Mo). Massachusetts Institute of Technology Reactor (MITR) is one of five research reactors selected for this program. The objective of this report is to provide preliminary thermal-hydraulic and mechanical analyses of the hydrodynamic effects in the MITR Design Demonstration Element (DDE) under conservative approximations for the plate power distribution. This report provides details on the modeling approach and the simulation results obtained, including pressure, flow velocity, temperature, and oxide layer over the design demonstration experiment for the irradiation cycles in the Belgian Reactor (BR)-2.

42 ENGINEERING↗
DOE OSTI2023

Preliminary NBSR Design-Demonstration Element Thermal-Hydraulics and Structural Analyses

1.1 OBJECTIVE The United State High Performance Research Reactor (USHPRR) program aims to eliminate more than 200kg of High Enriched Uranium (HEU) from commerce annually by converting five U.S. high-performance research reactors and one associated critical assembly to Low Enriched Uranium (LEU) fuel using a high-density alloy of uranium-10 wt% molybdenum (U-10Mo). The National Bureau of Standards Reactor (NBSR) is one of five research reactors selected for this program. The objective of this report is to provide preliminary thermal-hydraulic and mechanical analyses of the hydrodynamic effects in the NBSR Design Demonstration Element (DDE) under conservative approximations for the plate power distribution. This report provides details on the modeling approach and the simulation results obtained, including pressure, flow velocity, temperature, and oxide layer over the design demonstration experiment for the irradiation cycles in the Belgian Reactor (BR)-2.

42 ENGINEERING↗
DOE OSTI2021

A physics-based model for frost buildup under turbulent flow using direct numerical simulations

We present a new model for frost buildup under turbulent (and laminar) flow using direct numerical simulations. The physical model consists of two layers, the air and the frost. The air layer is fully resolved and consists of solving for the velocity, temperature, and vapor mass fraction fields. The frost layer thickness is resolved using conservation of mass and energy. Both phases are dynamically coupled using the immersed boundary method. Three-dimensional simulations are conducted in an open-channel configuration. A number of challenges need to be overcome to make these simulations feasible. First, to enforce far-field conditions of zero gradient and prescribed mean temperature and humidity, a source term is added to the energy and transport equations in the flow solver. Second, the mean frost thickness is subtracted after each time step to ensure a constant mean flow thickness and level of turbulence in the numerical domain. Third, a slow-time acceleration approach, which accelerates the frost buildup by a predetermined factor, is employed to bridge the gap between the fast turbulent and slow frost buildup time scales. Finally, a frost densification scheme is used to overcome the difficulties of vertically varying frost properties. The model is validated by comparing the frost thickness and frost thickness buildup rate over a period of one hour from a cooled flat plate experiment. As a result, both quantities compare favorably with experiments.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1993

Energy transfer and constrained simulations in isotropic turbulence

The defining characteristic of turbulent flows is their ability to dissipate energy, even in the limit of zero viscosity. The Euler equations, if constrained in such a way that the velocity derivatives remain bounded, conserve energy. But when they arise as the limit of the Navier-Stokes (NS) equations, when the Reynolds number goes to infinity, there is persuasive empirical evidence that the gradients become singular as just the right function of Re for the dissipation to remain non-zero and to approach a well defined limit. It is generally believed that this limiting value of the dissipation is a property of the Euler equations themselves, independent of the particular dissipative mechanism involved, and that it can be normalized with the large scale properties of the turbulent flow (e.g. the kinetic energy per unit volume u'(exp 2)/2, and the integral scale L) without reference to the Reynolds number or to other dissipative quantities. This is usually taken to imply that the low wave number end of the energy spectrum, far from the dissipative range, is also independent of the particular mechanism chosen to dispose of the energy transfer. In the following sections, we present some numerical experiments on the effect of substituting different dissipation models into the truncated Euler equations. We will see that the effect is mainly felt in the 'near dissipation' range of the energy spectrum, but that this range can be quite wide in some cases, contaminating a substantial range of wave numbers. In the process, we will develop a 'practical' approximation to the subgrid energy transfer in isotropic turbulence, and we will gain insight into the structure of the nonlinear interactions among turbulent scales of comparable size, and into the nature of energy backscatter. Some considerations on future research directions are offered at the end.

Jimenez, Javier↗
DOE OSTI2021

Flow condensation heat transfer performance of natural and emerging synthetic refrigerants

There is great interest in predicting flow condensation heat transfer for lower global warming potential (GWP) fluids. This paper analyzes the efficacy of common flow condensation correlations developed for particular fluids in order to identify their suitability to predict heat transfer performance of low GWP fluids. Condensation heat transfer data were extracted from the literature, including 19 papers and 1,473 data points for natural refrigerants [i.e., ammonia (R717), CO 2 (R744), propane (R290), isobutane (R600a)] and 35 papers and 5,030 data points for synthetic refrigerants [i.e., R12, R1234yf, R1234ze(E), R1234ze(Z), R22, R32, R41, R123, R125, R134a, R142b, R152a, R161, R404A, R410A, R448A, R449A, R450A, R452B, R454C, R455A, R513A] encompassing tube diameters of 0.1–11.5 mm, mass fluxes of 55–1200 kg/m2s, and saturation temperatures of -25°C–65°C. Correlations analyzed included Akers et al. (1959), Cavallini et al. (2006, 2011), Kim and Mudawar (2013), Macdonald and Garimella (2016), Shah (1979, 2009, 2013, 2016) and Traviss et al. (1973) for smooth tubes and Chamra et al. (2005) and Kedzierski and Goncalves (1999) for enhanced tubes. Since most studies did not report wall temperature, correlations which relied on wall temperature directly or indirectly were excluded from the analysis. For synthetic refrigerants, mean average error (MAE) ranged from 6%–257%, and Cavallini et al. (2011) and Kim and Mudawar (2013) were the best predictors for emerging synthetic refrigerants. Furthermore, the Kim and Mudawar (2013) correlation was found to best predict the heat transfer performance for propane and R600a data, but most correlations did not accurately predict ammonia and CO 2 flow condensation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
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-dependent three-dimensional analysis of the conservation of mass, momentum (Navier-Stokes), and energy equations. A boundary mapping technique was used to obtain a rectangular parallelepiped computational domain, 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. 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. Discussion of applying the results to a thermally bowed thermal protection system are presented.

Olsen, G. C.↗
NASA NTRS2002

Numerical Investigation of 'Transonic Resonance' with a Convergent-Divergent Nozzle

At pressure ratios lower than the design value, convergent-divergent (C-D) nozzles often undergo a flow resonance accompanied by the emission of acoustic tones. The phenomenon, driven by the unsteady shock within the divergent section of the nozzle, has been studied experimentally by Zaman et al. In this paper, the space-time conservation element solution element (CE/SE) method is employed to numerically investigate the phenomenon. The computations are performed for a given nozzle geometry for several different pressure ratios. Sustained 'limit cycle' oscillations are encountered in all cases. The oscillation frequencies, their variation with pressure ratio including a 'stage jump', agree well with the experimental results. The unsteady flow data confirm that stage 1 of the resonance (fundamental) involves a one-quarter standing wave while stage 2 (third harmonic) involves a three-quarter standing wave within the divergent section of the nozzle. Details of the shock motion, and the flow and near acoustic field, are documented for one case each of stages 1 and 2.

Loh, Ching Y.↗
NASA NTRS1984

A solution procedure for three-dimensional incompressible Navier-Stokes equation and its application

An implicit, finite-difference procedure is presented for numerically solving viscous incompressible flows. For convenience of applying the present method to three-dimensional problems, primitive variables, namely the pressure and velocities, are used. One of the major difficulties in solving incompressible flows that use primitive variables is caused by the pressure field solution method which is used as a mapping procedure to obtain a divergence-free velocity field. The present method is designed to accelerate the pressure-field solution procedure. This is achieved by the method of pseudocompressibility in which the time derivative pressure term is introduced into the mass conservation equation. The pressure wave propagation and the spreading of the viscous effect is investigated using simple test problems. The present study clarifies physical and numerical characteristics of the pseudo-compressible approach in simulating incompressible flows. Computed results for external and internal flows are presented to verify the present procedure. The present algorithm has been shown to be very robust and accurate if the selection of the pseudo-compressibility parameter has been made according to the guidelines given.

Kwak, D.↗
DOE OSTI2022

Conservative Numerical Schemes with Optimal Dispersive Wave Relations: Part II. Numerical Evaluations

A new energy and enstrophy conserving scheme (EEC) for the shallow water equations is proposed and evaluated using a suite of test cases over the global spherical or bounded domain. The evaluation is organized around a set of pre-defined properties: accuracy of individual operators, accuracy of the whole scheme, conservation of key quantities, control of the divergence variable, representation of the energy and enstrophy spectra, and simulation of nonlinear dynamics. The results confirm that the scheme is between the first and second order accurate, and conserves the total energy and potential enstrophy up to the time truncation errors. Here, the scheme is capable of producing more physically realistic energy and enstrophy spectra, indicating that it can help prevent the unphysical energy cascade towards the finest resolvable scales. With an optimal representation of the dispersive wave relations, the scheme is able to keep the flow close to being non-divergent, and maintain the geostrophically balanced structures with large-scale geophysical flows over long-term simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Baseline predictions of elliptic flow and fluctuations for the RHIC Beam Energy Scan using response coefficients

Currently the RHIC Beam Energy Scan is exploring a new region of the Quantum Chromodynamic phase diagram at large baryon densities that approaches nuclear astrophysics regimes. This provides an opportunity to study relativistic hydrodynamics in a regime where the net conserved charges of baryon number, strangeness, and electric charge play a role, which will significantly change the theoretical approach to simulating the baryon-dense Quark-Gluon Plasma. In this work, we detail many of the important changes needed to adapt both initial conditions and the medium to baryon-rich matter. Then, we make baseline predictions for the elliptical flow and fluctuations based on extrapolating the physics at LHC and top RHIC energies to support future analyses of where and how the new baryon-dense physics causes these extrapolations to break down. First we compare eccentricities across beam energies, exploring their underlying assumptions; we find the the extrapolated initial state is predicted to be nearly identical to that at AuAu $\sqrt{s_{NN}}$ = 200 GeV. Then the final flow harmonic predictions are based on linear+cubic response. We discuss preliminary STAR results in order to determine the implications that they have for linear+cubic response coefficients at the lowest beam energy of AuAu $\sqrt{s_{NN}}$ = 7 GeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1999

Dynamics and Structure of Dusty Reacting Flows: Inert Particles in Strained, Laminar, Premixed Flames

A detailed numerical study was conducted on the dynamics and thermal response of inert, spherical particles in strained, laminar, premixed hydrogen/air flames. The modeling included the solution of the steady conservation equations for both the gas and particle phases along and around the stagnation streamline of an opposed-jet configuration, and the use of detailed descriptions of chemical kinetics and molecular transport, For the gas phase, the equations of mass, momentum, energy, and species are considered, while for the particle phase, the model is based on conservation equations of the particle momentum balance in the axial and radial direction, the particle number density, and the particle thermal energy equation. The particle momentum equation includes the forces as induced by drag, thermophoresis, and gravity. The particle thermal energy equation includes the convective/conductive heat exchange between the two phases, as well as radiation emission and absorption by the particle. A one-point continuation method is also included in the code that allows for the description of turning points, typical of ignition and extinction behavior. As expected, results showed that the particle velocity can be substantially different than the gas phase velocity, especially in the presence of large temperature gradients and large strain rates. Large particles were also found to cross the gas stagnation plane, stagnate, and eventually reverse as a result of the opposing gas phase velocity. It was also shown that the particle number density varies substantially throughout the flowfield, as a result of the straining of the flow and the thermal expansion. Finally, for increased values of the particle number density, substantial flame cooling to extinction states and modification of the gas phase fluid mechanics were observed. As also expected, the effect of gravity was shown to be important for low convective velocities and heavy particles. Under such conditions, simulations indicate that the magnitude and direction of the gravitational force can substantially affect the profiles of the particle velocity, number density, mass flux, and temperature.

Egolfopoulos, Fokion N.↗
NASA NTRS1999

Dynamics and Structure of Dusty Reacting Flows: Inert Particles in Strained, Laminar, Premixed Flames

A detailed numerical study was conducted on the dynamics and thermal response of inert spherical particles in strained, laminar, premixed hydrogen/air flames. The modeling included the solution of the steady conservation equations for both the gas and particle phases along and around the stagnation streamline of an opposed-jet configuration, and the use of detailed descriptions of chemical kinetics and molecular transport. For the gas phase, the equations of mass, momentum, energy, and species are considered, while for the particle phase, the model is based on conservation equations of the particle momentum balance in the axial and radial direction, the particle number density, and the particle thermal energy equation. The particle momentum equation includes the forces as induced by drag, thermophoresis, and gravity. The particle thermal energy equation includes the convective/conductive heat exchange between the two phases, as well as radiation emission and absorption by the particle. A one-point continuation method is also included in the code that allows for the description of turning points, typical of ignition and extinction behavior. As expected, results showed that the particle velocity can be substantially different than the gas phase velocity, especially in the presence of large temperature gradients and large strain rates. Large particles were also found to cross the gas stagnation plane, stagnate, and eventually reverse as a result of the opposing gas phase velocity. It was also shown that the particle number density varies substantially throughout the flowfield, as a result of the straining of the flow and the thermal expansion. Finally, for increased values of the particle number density, substantial flame cooling to extinction states and modification of the gas phase fluid mechanics were observed. As also expected, the effect of gravity was shown to be important for low convective velocities and heavy particles. Under such conditions, simulations indicate that the magnitude and direction of the gravitational force can substantially affect the profiles of the particle velocity, number density, mass flux, and temperature.

Egolfopoulos, Fokion N.↗
DOE OSTI2020

A GPU-based 2D shallow water quality model

In this paper, a 2D shallow water flow solver integrated with a water quality model is presented. The interaction between the main water quality constituents included is based on the Water Quality Analysis Simulation Program. Efficiency is achieved by computing with a combination of a Central Processing Unit (CPU) and a Graphics Processing Unit (GPU) device. This technique is intended to provide robust and accurate simulations with high computation speedups with respect to a single-core CPU in real events. The proposed numerical model is evaluated in cases that include the transport and reaction of water quality components over irregular bed topography and dry–wet fronts, verifying that the numerical solution in these situations conserves the required properties (C-property and positivity). The model can operate in any steady or unsteady form allowing an efficient assessment of the environmental impact of water flows. The field data from an unsteady river reach test case are used to show that the model is capable of predicting the measured temporal distribution of dissolved oxygen and water temperature, proving the robustness and computational efficiency of the model, even in the presence of noisy signals such as wind speed.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1986

Divergence study of a high-aspect ratio, forward-swept wing

An experimental wind-tunnel study to determine the divergence characteristics of a high-aspect ratio, forward-swept wing has been conducted in the NASA Langley Research Center (LaRC) Transonic Dynamics Tunnel (TDT). The rectangular wing used for this study had a panel aspect ratio of 9.16 (lambda = 0 deg.) and the sweep angle could be set at lambda = 0 deg., -15 deg., -30 deg., -45 deg., or -60 deg. A rectangular wing tip shape was tested at each of these sweep angles. In addition, a tip shape parallel to the freestream flow was tested for a wing sweep angle of lambda = -45 deg. The root of the wing was cantilever mounted to the wall of the wind tunnel. Divergence conditions were measured at M = 0.4 for each sweep angle and tip configuration tested. Subcritical response techniques were used to extrapolate to the divergence conditions during the wind-tunnel test. The primary objective of this test was to obtain data which could be used to verify for this configuration the divergence prediction capability of an aeroelastic analysis code. Subsonic lifting surface theory (kernel function) aerodynamics are utilized by this particular code. The analytical predictions of divergence were found to be significantly conservative at all forward sweep angles. At lambda = -45 deg., the analysis was 14 percent conservative. The effect of the two tip shapes on the divergence dynamic pressure was predicted accurately by the analysis. The divergence condition for the tip shape parallel to the flow occurred at a dynamic pressure 14 percent higher than the divergence condition with a rectangular tip shape.

Cole, S. R.↗