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 829 records · Page 46

A new consistent spatial differencing scheme for the transonic full-potential equation

A new spatial differencing scheme for the transonic full-potential equation in conservative form has been developed. This scheme guarantees zero truncation error on any curvilinear mesh for freestream flows in either two- or three-space dimensions. Solutions obtained with this new differencing scheme, away from freestream regions, exhibit greatly improved accuracy, especially for nonsmooth or singular meshes.

Flores, J.↗

Aerothermal loads analysis for high speed flow over a quilted surface configuration

Attention is given to hypersonic laminar flow over a quilted surface configuration that simulates an array of Space Shuttle Thermal Protection System panels bowed in a spherical shape as a result of thermal gradient through the panel thickness. Pressure and heating loads to the surface are determined. The flow field over the configuration was mathematically modeled by means of time-dependent, three-dimensional conservation of mass, momentum, and energy equations. A boundary mapping technique was then used to obtain a rectangular, parallel piped computational domain, and an explicit MacCormack (1972) explicit time-split predictor corrector finite difference algorithm was used to obtain steady state solutions. Total integrated heating loads vary linearly with bowed height when this value does not exceed the local boundary layer thickness.

Olsen, G. C.↗

Aerothermal loads analysis for high speed flow over a quilted surface configuration

Attention is given to hypersonic laminar flow over a quilted surface configuration that simulates an array of Space Shuttle Thermal Protection System panels bowed in a spherical shape as a result of thermal gradients through the panel thickness. Pressure and heating loads to the surface are determined. The flow field over the configuration was mathematically modeled by means of time-dependent, three-dimensional conservation of mass, momentum, and energy equations. A boundary mapping technique was then used to obtain a rectangular, parallelepiped computational domain, and an explicit MacCormack (1972) explicit time-split predictor-corrector finite difference algorithm was used to obtain steady state solutions. Total integrated heating loads vary linearly with bowed height when this value does not exceed the local boundary layer thickness.

Olsen, G. C.↗

Conservative multizonal interface algorithm for the 3-D Navier-Stokes equations

A conservative zonal interface algorithm using features of both structured and unstructured mesh CFD technology is presented. The flow solver within each of the zones is based on structured mesh CFD technology. The interface algorithm was implemented into two three-dimensional Navier-Stokes finite volume codes and was found to yield good results.

Klopfer, G. H.↗

Simulation of cryogenic turbopump annular seals

San Andres employed the NBS software package MIPROPS to account for density's dependence on pressure in the simulation of liquid annular seals. His example on a LH2 seal showed a significant change in the mass coefficient compared to a constant density model. San Andres, Yang, and Childs extended this analysis by including the pressure and temperature dependence of density, specific heat, viscosity, volumetric expansion, and thermal conductivity in a coupled solution of the energy, momentum, and continuity equations. Their example showed very significant changes in stiffness and inertia for a high speed (38,000 rpm), large L/D ratio (0.5) LOX seal, as compared to their constant temperature results. The current research rederived the San Andres-Yang-Childs (SYC) analysis and extended it to include not only the Moody friction model of SYC but also the Hir's friction model. The derivation begins with obtaining the local differential equations of continuity, momentum, and energy conservation in the seal. These equations are averaged across the film thickness to obtain the resulting 'bulk flow' differential equations. Shear stress and convective heat loss through the stator (seal) and rotor are related to the Moody and Hir's friction factor model. The Holman analogy is employed to relate heat conduction in or out of the fluid film's boundary layer to the friction induced shear stress.

Palazzolo, Alan B.↗

Power Controller

The device called the Power Factor Controller (PFC) offers exceptional energy conservation potential by virtue of its ability to sense shifts in the relationship between voltage and current flow, and to match them with the motor's need. Originating from the solar heating/cooling program, the PFC senses a light load, it cuts the voltage level to the minimum needed which in turn reduces current flow and heat loss. Laboratory tests showed that the PFC could reduce power used by six to eight percent under normal motor loads, and as much as 65 percent when the motor was idling. Over 150 companies have been granted NASA licenses for commercial use of this technology. One system that utilizes this technology is the Vectrol Energy System, (VES) produced by Vectrol, Inc. a subsidiary of Westinghouse. The VES is being used at Woodward & Lothrop, on their escalators. Energy use is regulated according to how many people are on the escalator at any time. It is estimated that the energy savings are between 30 to 40 percent.

Source record↗

The Application of Nonstandard Analysis to the Study of Inviscid Shock Wave Jump Conditions

The use of conservation laws in nonconservative form for deriving shock jump conditions by Schwartz distribution theory leads to ambiguous products of generalized functions. Nonstandard analysis is used to define a class of Heaviside functions where the jump from zero to one occurs on an infinitesimal interval. These Heaviside functions differ by their microstructure near x = 0, i.e., by the nature of the rise within the infinitesimal interval it is shown that the conservation laws in nonconservative form can relate the different Heaviside functions used to define jumps in different flow parameters. There are no mathematical or logical ambiguities in the derivation of the jump conditions. An important result is that the microstructure of the Heaviside function of the jump in entropy has a positive peak greater than one within the infinitesimal interval where the jump occurs. This phenomena is known from more sophisticated studies of the structure of shock waves using viscous fluid assumption. However, the present analysis is simpler and more direct.

Farassat, F.↗

Hypersonic Viscous Flow Over Large Roughness Elements

Viscous flow over discrete or distributed surface roughness has great implications for hypersonic flight due to aerothermodynamic considerations related to laminar-turbulent transition. Current prediction capability is greatly hampered by the limited knowledge base for such flows. To help fill that gap, numerical computations are used to investigate the intricate flow physics involved. An unstructured mesh, compressible Navier-Stokes code based on the space-time conservation element, solution element (CESE) method is used to perform time-accurate Navier-Stokes calculations for two roughness shapes investigated in wind tunnel experiments at NASA Langley Research Center. It was found through 2D parametric study that at subcritical Reynolds numbers of the boundary layers, absolute instability resulting in vortex shedding downstream, is likely to weaken at supersonic free-stream conditions. On the other hand, convective instability may be the dominant mechanism for supersonic boundary layers. Three-dimensional calculations for a rectangular or cylindrical roughness element at post-shock Mach numbers of 4.1 and 6.5 also confirm that no self-sustained vortex generation is present.

Chang, Chau-Lyan↗

Passive Rocket Diffuser Theory: A Re-Examination of Minimum Second Throat Size

Second-throat diffusers serve to isolate rocket engines from the effects of ambient back pressure during testing without using active control systems. Among the most critical design parameters is the relative area of the diffuser throat to that of the nozzle throat. A smaller second throat is generally desirable because it decreases the stagnation-to-ambient pressure ratio the diffuser requires for nominal operation. There is a limit, however. Below a certain size, the second throat can cause pressure buildup within the diffuser and prevent it from reaching the start condition that protects the nozzle from side-load damage. This paper presents a method for improved estimation of the minimum second throat area which enables diffuser start. The new 3-zone model uses traditional quasi-one-dimensional compressible flow theory to approximate the structure of two distinct diffuser flow fields observed in Computational Fluid Dynamics (CFD) simulations and combines them to provide a less-conservative estimate of the second throat size limit. It is unique among second throat sizing methods in that it accounts for all major conical nozzle and second throat diffuser design parameters within its limits of application. The performance of the 3-zone method is compared to the historical normal shock and force balance methods, and verified against a large number of CFD simulations at specific heat ratios of 1.4 and 1.25. Validation is left as future work, and the model is currently intended to function only as a first-order design tool.

Jones, Daniel R.↗

Implementation of Manifold-Based Combustion Models in a Highly Scalable Low Mach Number Reacting Flow Solver: Preprint

Manifold-based representations of the thermochemistry are often employed in conjunction with large eddy eimulation (LES) to lower the cost of combustion simulations. This work describes steps taken to implement this modeling approach in PeleLM, a scalable and performance-portable low Mach number flow solver. Most significantly, this includes adapting the projection method used by PeleLM to satisfy the mass conservation constraint for use with manifold-based models. The implementation is designed to be general across manifold-based models, including both those that employ traditional tabulation and those that employ neural networks. An initial demonstration for simple test cases is presented and will be used for performance assessment.

high-performance computing↗

A Computational Icing Effects Study for a Three-Dimensional Wing: Comparison with Experimental Data and Investigation of Spanwise Variation

Steady-state solutions of the Reynolds-averaged Navier-Stokes (RANS) equations were computed using the Colbalt flow solver for a constant-section, rectangular wing based on an extruded two-dimensional glaze ice shape. The one equation Spalart-Allmaras turbulence model was used. The results were compared with data obtained from a recent wind tunnel test. Computed results indicate that the steady RANS solutions do not accurately capture the recirculating region downstream of the ice accretion, even after a mesh refinement. The resulting predicted reattachment is farther downstream than indicated by the experimental data. Additionally, the solutions computed on a relatively coarse baseline mesh had detailed flow characteristics that were different from those computed on the refined mesh or the experimental data. Steady RANS solutions were also computed to investigate the effects of spanwise variation in the ice shape. The spanwise variation was obtained via a bleeding function that merged the ice shape with the clean wing using a sinusoidal spanwise variation. For these configurations, the results predicted for the extruded shape provided conservative estimates for the performance degradation of the wing. Additionally, the spanwise variation in the ice shape and the resulting differences in the flow fields did not significantly change the location of the primary reattachment.

Thompson, D.↗

Semi-coupled flow and thermal analysis of the field joint during rapid pressurization of the redesigned solid rocket motor

A transient, semi-coupled, multi-dimensional thermal and flow analysis methodology was developed to predict the thermal/gas dynamic conditions in the field joint region of the Redesigned Solid Rocket Motor (RSRM). Transient temperature response, pressure history, and combustion gas flow rate (within the field joint region), were of principle interest, in the course of this study. The thermal environment in the field joint was modeled using SINDA, a finite difference based thermal network analyzer. The combustion gas flow boundary condition was generated using the FLAP code; this code performs a transient, lumped-parameter, control volume analysis to solve the mass, momentum, and energy conservation equations. The FLAP computer code was modified to account for erosion of the NBR insulation material, following ignition. An independent grid sensitivity study was conducted to determine an appropriate grid distribution near the wall. The predicted results, obtained using an optimum grid distribution and computer generated flow boundary condition, were compared with subscale test data.

Ghaffarian, Benny↗

Preconditioning for Numerical Simulation of Low Mach Number Three-Dimensional Viscous Turbomachinery Flows

A preconditioning scheme has been implemented into a three-dimensional viscous computational fluid dynamics code for turbomachine blade rows. The preconditioning allows the code, originally developed for simulating compressible flow fields, to be applied to nearly-incompressible, low Mach number flows. A brief description is given of the compressible Navier-Stokes equations for a rotating coordinate system, along with the preconditioning method employed. Details about the conservative formulation of artificial dissipation are provided, and different artificial dissipation schemes are discussed and compared. The preconditioned code was applied to a well-documented case involving the NASA large low-speed centrifugal compressor for which detailed experimental data are available for comparison. Performance and flow field data are compared for the near-design operating point of the compressor, with generally good agreement between computation and experiment. Further, significant differences between computational results for the different numerical implementations, revealing different levels of solution accuracy, are discussed.

Tweedt, Daniel L.↗

An investigation of shock wave turbulent boundary layer interaction with bleed through slanted slots

Flow-field characteristics are simulated numerically in an oblique shock wave/turbulent boundary layer interactions with six different bleed slot configurations. The strong conservation-law form of the 2D compressible Navier-Stokes equations and the k-epsilon equations are solved throughout the interaction region and inside the bleed slot. The computed results are presented for a normal and 20-deg slanted bleed slots at three different locations, upstream, across and downstream of the impingement point of an oblique shock of sufficient strength to cause boundary layer separation without bleed. The detailed flow characteristics in the interaction zone and inside the bleed slot are compared for the different bleed slot configurations. The resulting surface pressure and shear stress distributions as well as the boundary layer characteristics downstream of the interaction region are also presented for the mix bleed configurations at different bleed mass flows up to choking.

Hamed, A.↗

Physiological water model development

The water of the human body can be categorized as existing in two main compartments: intracellular water and extracellular water. The intracellular water consists of all the water within the cells and constitutes over half of the total body water. Since red blood cells are surrounded by plasma, and all other cells are surrounded by interstitial fluid, the intracellular compartment has been subdivided to represent these two cell types. The extracellular water, which includes all of the fluid outside of the cells, can be further subdivided into compartments which represent the interstitial fluid, circulating blood plasma, lymph, and transcellular water. The interstitial fluid surrounds cells outside of the vascular system whereas plasma is contained within the blood vessels. Avascular tissues such as dense connective tissue and cartilage contain interstitial water which slowly equilibrates with tracers used to determine extracellular fluid volume. For this reason, additional compartments are sometimes used to represent these avascular tissues. The average size of each compartment, in terms of percent body weight, has been determined for adult males and females. These compartments and the forces which cause flow between them are presented. The kidneys, a main compartment, receive about 25 percent of the cardiac output and filters out a fluid similar to plasma. The composition of this filtered fluid changes as it flows through the kidney tubules since compounds are continually being secreted and reabsorbed. Through this mechanism, the kidneys eliminate wastes while conserving body water, electrolytes, and metabolites. Since sodium accounts for over 90 percent of the cations in the extracellular fluid, and the number of cations is balanced by the number of anions, considering the renal handling sodium and water only should sufficiently describe the relationship between the plasma compartment and kidneys. A kidney function model is presented which has been adapted from a previous model of normal renal function in man. To test the validity of the proposed kidney model, results predicted by the model will be compared to actual data involving injected or ingested fluids and subsequent urine flow rates. Comparison of the model simulation to actual data following the ingestion of 1 liter of water is shown. The model simulation is also shown with actual data following the intravenous infusion of hypertonic saline.

Doty, Susan↗

Advanced air transport concepts

The concepts of laminar flow control, very large all-wing aircraft, an aerial relay transportation system and alternative fuels, which would enable large improvements in fuel conservation in air transportation in the 1990's are discussed. Laminar boundary layer control through suction would greatly reduce skin friction and has been reported to reduce fuel consumption by up to 29%. Distributed load aircraft, in which all fuel and payload are carried in the wing and the fuselage is absent, permit the use of lighter construction materials and the elimination of fuselage and tail drag. Spanloader aircraft with laminar flow control could be used in an aerial relay transportation system which would employ a network of continuously flying liners supplied with fuel, cargo and crews by smaller feeder aircraft. Liquid hydrogen and methane fuels derived from coal are shown to be more weight efficient and less costly than coal-derived synthetic jet fuels.

Molloy, J. K.↗

The CE/SE Method: a CFD Framework for the Challenges of the New Millennium

The space-time conservation element and solution element (CE/SE) method, which was originated and is continuously being developed at NASA Glenn Research Center, is a high-resolution, genuinely multidimensional and unstructured-mesh compatible numerical method for solving conservation laws. Since its inception in 1991, the CE/SE method has been used to obtain highly accurate numerical solutions for 1D, 2D and 3D flow problems involving shocks, contact discontinuities, acoustic waves, vortices, shock/acoustic waves/vortices interactions, shock/boundary layers interactions and chemical reactions. Without the aid of preconditioning or other special techniques, it has been applied to both steady and unsteady flows with speeds ranging from Mach number = 0.00288 to 10. In addition, the method has unique features that allow for (i) the use of very simple non-reflecting boundary conditions, and (ii) a unified wall boundary treatment for viscous and inviscid flows. The CE/SE method was developed with the conviction that, with a solid foundation in physics, a robust, coherent and accurate numerical framework can be built without involving overly complex mathematics. As a result, the method was constructed using a set of design principles that facilitate simplicity, robustness and accuracy. The most important among them are: (i) enforcing both local and global flux conservation in space and time, with flux evaluation at an interface being an integral part of the solution procedure and requiring no interpolation or extrapolation; (ii) unifying space and time and treating them as a single entity; and (iii) requiring that a numerical scheme be built from a nondissipative core scheme such that the numerical dissipation can be effectively controlled and, as a result, will not overwhelm the physical dissipation. Part I of the workshop will be devoted to a discussion of these principles along with a description of how the ID, 2D and 3D CE/SE schemes are constructed. In Part II, various applications of the CE/SE method, particularly those involving chemical reactions and acoustics, will be presented. The workshop will be concluded with a sketch of the future research directions.

Chang, Sin-Chung↗

Flight-Management Algorithm for Fuel-Conservative Descents

Federal Aviation Administration has developed an automated time-based metering form of air traffic control for arrivals into terminal area called local flow management/profile descent (LFM/PD). LFM/PD saves fuel by matching airplane arrival flow to airport acceptance rate through time-control computations and by allowing pilot to descend at his discretion from cruise altitude to metering fix in an idle-thrust, clean configuration (landing gear up, flaps zero, speed brakes retracted).

Knox, C. E.↗