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 1,351 records · Page 75

Numerical studies of the formation and destruction of vortices in a motored four-stroke piston-cylinder configuration

A finite-difference procedure which solves the conservation equations of mass, momentum, and energy is used to investigate the effects of the compression ratio, engine speed, bore-to-stroke ratio, and air intake flow angle on the turbulent flow field within an axisymmetric piston-cylinder configuration. It is shown that in a four-stroke piston-cylinder configuration, the intake stroke is characterized by the formation of a piston vortex. The piston vortex is stretched during the intake stroke, and the head vortex has an almost constant diameter. For a 0-deg air intake flow angle, both vortices disappear by the end of the compression stroke; for an air intake flow angle of 45 deg, the flow field within the cylinder shows three elongated vortices which persist into the compression stroke and then break up and merge. It is also shown that larger bore-to-stroke ratios give rise to lower turbulent levels than smaller bore-to-stroke ratios and that the turbulent intensity is almost independent of the rpm.

Schock, H. J.↗

Fresh attack on laminar flow

Laminar flow control is a technology with great potential for aircraft drag reduction. Stabilization of laminar boundary layers became known as natural laminar flow (NLF) and research led to the development of NLF airfoils. Research was also conducted on stabilization by suction, referred to as laminar flow control (LFC). Experiments demonstrated that extensive laminar flow could be achieved in flight. However, there remained doubts regarding the practicality of producing, with the technology then available, wing surfaces sufficiently smooth and wavefree to meet laminar-flow criteria and maintaining the wing surface quality in normal service. In 1976, the Aircraft Energy Efficiency (ACEE) program was begun by NASA to develop fuel-conservative technology for commercial transports. The progress of the ACEE program is discussed. Attention is given to LFC wing structures, and LFC leading-edge systems.

Wagner, R. D.↗

Drag characteristics of unsteady, perturbed boundary layer flows

A series of time-dependent numerical computations have been performed for the flow of an unsteady boundary layer over an embedded cavity, which is aligned normal to the flow direction. The unsteady flow is the result of the imposition of a Stuart vortex, at the inflow boundary of the computational domain, onto an otherwise unperturbed two-dimensional laminar boundary layer flow. This produces a strongly vortical, but spatially monochromatic, motion at inflow which is allowed to evolve downstream through the time-dependent numerical solution of appropriate conservation equations. As the imposed vortex evolves over the solid boundary, a weaker induced vortical motion is produced and together these perturbing motions interact with the fluid motion in the embedded cavity. An analysis of the pertinent dynamic variables as well as the time-dependent drag characteristics for this type of non-planar geometry is performed; with particular emphasis on the relative contribution of both pressure drag and frictional drag at various times in the evolution process.

Gatski, T. B.↗

Numerical modeling of multidimensional flow in seals and bearings used in rotating machinery

The rotordynamic behavior of turbomachinery is critically dependent on fluid dynamic rotor forces developed by various types of seals and bearings. The occurrence of self-excited vibrations often depends on the rotor speed and load. Misalignment and rotor wobbling motion associated with differential clearance were often attributed to stability problems. In general, the rotative character of the flowfield is a complex three dimensional system with secondary flow patterns that significantly alter the average fluid circumferential velocity. A multidimensional, nonorthogonal, body-fitted-grid fluid flow model is presented that describes the fluid dynamic forces and the secondary flow pattern development in seals and bearings. Several numerical experiments were carried out to demonstrate the characteristics of this complex flowfield. Analyses were performed by solving a conservation form of the three dimensional Navier-Stokes equations transformed to those for a rotating observer and using the general-purpose computer code PHOENICS with the assumptions that the rotor orbit is circular and that static eccentricity is zero. These assumptions have enabled a precise steady-state analysis to be used. Fluid injection from ports near the seal or bearing center increased fluid-film direct dynamic stiffness and, in some cases, significantly increased quadrature dynamic stiffness. Injection angle and velocity could be used for active rotordynamic control; for example, injection, when compared with no injection, increased direct dynamic stiffness, which is an important factor for hydrostatic bearings.

Hendricks, R. C.↗

Numerical and analytical study of fluid dynamic forces in seals and bearings

Study results for rotor/bearing/seal system models are compared with both a numerical model based on a transformed, conservative version of the three-dimensional Navier-Stokes equation, and an analytical model based on 'lumped' fluid parameters. Rotor-destabilizing factors are related to flow field rotation; it is shown that these factors can be reduced through a decrease in the fluid average circumferential velocity. A transformation of the Navier-Stokes equations to those for a rotating observer and a numerical code with a nonorthogonal body-fitted grid are used to carry out numerical experiments demonstrating the character of the complex flow field.

Tam, L. T.↗

The effect of surface radiation on flame spread in a quiescent, microgravity environment

In most theories of laminar flame spread over solid combustibles, radiation is neglected. However, when gas motion is completely absent, which may occur in the quiescent, microgravity environment of a spacecraft, the importance of radiation compared to convection is enhanced. Here, a theoretical model is developed of flame spread over a thin solid fuel into an opposing flow of oxidizer, the opposing flow being present in the quiescent environment in flame-fixed coordinates, including the effects of radiative heat transfer from the fuel surface. Numerical solutions to the conservation equations in the gas and the solid phase that describe the spreading flame are presented for a variety of ambient conditions and surface emittances. When surface radiation is significant, the solid surface acts as a heat sink, causing the flame to cool, shrink in size, and spread more slowly. Better agreement with experiment for net heat transfer to the surface is obtained when surface radiation is included. For low ambient oxygen levels, radiative effects lead to flame extinction. A radiation/conduction parameter is identified that adequately describes the importance of surface radiation.

Bhattacharjee, S.↗

Computational Aeroacoustics by the Space-Time CE/SE Method

The Space-Time Conservation Element and Solution Element Method, or CE/SE Method for short, is a newly developed numerical method for conservation laws. Despite its second order accuracy, it possesses low dispersion errors and low dissipation. The method is robust enough to cover a wide spectrum of compressible flows: from weak linear acoustic waves to strong discontinuous waves (shocks). An outstanding feature of the CE/SE scheme is its novel, simple but effective non-reflecting boundary condition (NRBC), which is particularly valuable for CAA (computational aeroacoustics). In this seminar, the 1-D and 2-D unstructured version of the CE/SE schemes are first briefly described. Secondly, some discussions on the NRBC are given. Then, various examples for linear, nonlinear aeroacoustics are presented.

Loh, Ching Y.↗

Numerical Modeling of Saturated Boiling in a Heated Tube

This paper describes a mathematical formulation and numerical solution of boiling in a heated tube. The mathematical formulation involves a discretization of the tube into a flow network consisting of fluid nodes and branches and a thermal network consisting of solid nodes and conductors. In the fluid network, the mass, momentum and energy conservation equations are solved and in the thermal network, the energy conservation equation of solids is solved. A pressure-based, finite-volume formulation has been used to solve the equations in the fluid network. The system of equations is solved by a hybrid numerical scheme which solves the mass and momentum conservation equations by a simultaneous Newton-Raphson method and the energy conservation equation by a successive substitution method. The fluid network and thermal network are coupled through heat transfer between the solid and fluid nodes which is computed by Chen's correlation of saturated boiling heat transfer. The computer model is developed using the Generalized Fluid System Simulation Program and the numerical predictions are compared with test data.

Majumdar, Alok↗

Entropy Stable Method for the Euler Equations Revisited: Central Differencing via Entropy Splitting and SBP

The two decades old high order central differencing via entropy splitting and summation-by-parts (SBP) difference boundary closure of Ols- son & Oliger (1994), Gerritsen & Olsson (1996), and Yee et al. (2000) is revisited. The entropy splitting is a form of skew-symmetric splitting of the nonlinear Euler flux derivatives. Central differencing applied to the entropy splitting form of the Euler flux derivatives together with SBP difference operators will, hereafter, be referred to as entropy split schemes. This study is prompted by the recent growing interest in numerical methods for which a discrete entropy conservation law holds, a discrete global entropy conservation can be proved and/or the numerical method possesses a stable entropy in the framework of SBP difference operators and L2-energy norm estimate. The objective of this paper is to recast the entropy split scheme as the re- cent definition of an entropy stable method for central differencing with SBP operators for both periodic and non-periodic boundary conditions for non- linear Euler equations. Standard high order spatial central differencing as well as high order central spatial DRP (dispersion relation preserving) spatial differencing is part of the entropy stable methodology framework. Long time integration of 2D and 3D test cases is included to show the comparison of this efficient entropy stable method with the Tadmor-type of entropy conservative methods. Studies also include the comparison among the three skew-symmetric splittings on their nonlinear stability and accuracy performance without added numerical dissipations for smooth flows. These are, namely, entropy splitting, Ducros et al. splitting and the Kennedy & Grub- ber splitting.

Sjogreen, Bjorn↗

Entropy Stable Method for the Euler Equations Revisited: Central Differencing via Entropy Splitting and SBP

The two decades old high order central differencing via entropy splitting and summation-by-parts (SBP) difference boundary closure of Olsson & Oliger, Gerritsen & Olsson, and Yee et al. (15, 7, 37) is revisited. The objective of this paper is to prove for the first time that the entropy split scheme is an entropy stable method for central differencing with SBP operators for both periodic and non-periodic boundary conditions for nonlinear Euler equations. Standard high order spatial central differencing as well as high order central spatial DRP (dispersion relation preserving) spatial differencing is part of the entropy stable methodology framework. The proof is to replace the spatial derivatives by summation-by-parts (SBP) difference operators in the entropy split form of the equations using the physical entropy of the Euler equations. The numerical boundary closure follows directly from the SBP operator. No additional numerical boundary procedure is required. In contrast, Tadmor-type entropy conserving schemes (31) using mathematical entropies and more recently in (35], do not naturally come with a numerical boundary closure and a generalized SBP operator has to be developed (18). Long time integration of 2D and 3D test cases is included to show the comparison of this efficient entropy stable method with the Tadmor-type of entropy conservative methods. Studies also include the comparison among the three skew-symmetric splittings on their nonlinear stability and accuracy performance without added numerical dissipations for smooth flows. These are, namely, entropy splitting, Ducros et al. splitting and the Kennedy & Grubber splitting.

Sjogreen, Bjorn↗

Effect of a Diffuser on Conditioning Flow Field Fluctuations at the Exit of a Rotating Detonation Combustor

Pressure gain combustion (PGC) seeks to convert fuel’s chemical energy simultaneously into thermal energy and mechanical energy, thereby reducing the entropy production in the process. Recent research has shown that the rotating detonation combustion or combustor (RDC) can provide excellent specific thrust, specific impulse, and large pressure gain through rapid energy release by continuous detonation in the circumferential direction. In this study, we present experimental results from RDC operated on methane and oxygen-enriched air to represent reactants used in gas turbines for land-based power generation. A thorough understanding of the velocity flow field exiting the combustor is important to develop future rotating detonation engine (RDE) designs, as well as to successfully implement the flow conditioning devices to power generating turbine hardware located downstream. The RDC is operated at high pressures in two ways through this study: 1) placing a back-pressure plate (BPP) downstream of the annular combustor, and 2) placing a diffuser downstream of the annular combustor. Time-resolved particle image velocimetry (TR-PIV) is applied at a frame rate of 50 kHz to quantify the RDC exit flow. A qualitative analysis is presented with vector data obtain from experimental testing. Results show that for the BPP configuration the highly periodic exit flow tilts in the azimuthal direction, indicating a significant circumferential velocity component. For the diffuser configuration the flow is dominantly axial, with an insignificant circumferential velocity component. The flow field undergoes cyclic fluctuations with the BPP, but not so with the diffuser. Results also show the challenges associated with achieving uniform seed distribution in this highly transient and rapidly changing flow field.

03 NATURAL GAS↗

Implicit solution of three-dimensional internal turbulent flows

The scalar form of the approximate factorization method was used to develop a new code for the solution of three-dimensional internal laminar and turbulent compressible flows. The Navier-Stokes equations in their Reynolds-averaged form are iterated in time until a steady solution is reached. Evidence is given to the implicit and explicit artificial damping schemes that proved to be particularly efficient in speeding up convergence and enhancing the algorithm robustness. A conservative treatment of these terms at domain boundaries is proposed in order to avoid undesired mass and/or momentum artificial fluxes. Turbulence effects are accounted for by the zero-equation Baldwin-Lomax turbulence model and the q-omega two-equation model. For the first, an investigation on the model behavior in case of multiple boundaries is performed. The flow in a developing S-duct is then solved in the laminar regime at Reynolds number (Re) 790 and in the turbulent regime at Re=40,000 using the Baldwin-Lomax model . The Stanitz elbow is then solved using an inviscid version of the same code at M(sub inlet)=0.4. Grid dependence and convergence rate are investigated showing that for this solver the implicit damping scheme may play a critical role for convergence characteristics. The same flow at Re=2.5x10(exp 6) is solved with the Baldwin-Lomax and the q-omega models. Both approaches showed satisfactory agreement with experiments, although the q-omega model is slightly more accurate.

Michelassi, V.↗

Gulf stream velocity structure through combined inversion of hydrographic and acoustic Doppler data

Near-surface velocities from an acoustic Doppler instrument are used in conjunction with CTD/O2 data to produce estimates of the absolute flow field off Cape Hatteras. The data set consists of two transects across the Gulf Stream made by the R/V Endeavor cruise EN88 in August 1982. An inverse procedure is applied which makes use of both the acoustic Doppler data and property conservation constraints. Velocity sections at approximately 73 deg. W and 71 deg. W are presented with formal errors of 1-2 cm/s. The net Gulf Stream transports are estimated to be 116 + or - 2 Sv across the south leg and 161 + or - 4 Sv across the north. A Deep Western Boundary Current transport of 4 + or - 1 Sv is also estimated. While these values do not necessarily represent the mean, they are accurate estimates of the synoptic flow field in the region.

Pierce, S. D.↗

Numerical Study of Base Pressure Characteristic Curve for a Four-Engine Clustered Nozzle Configuration

Excessive base heating has been a problem for many launch vehicles. For certain designs such as the direct dump of turbine exhaust in the nozzle section and at the nozzle lip of the Space Transportation Systems Engine (STME), the potential burning of the turbine exhaust in the base region has caused tremendous concern. Two conventional approaches have been considered for predicting the base environment: (1) empirical approach, and (2) experimental approach. The empirical approach uses a combination of data correlations and semi-theoretical calculations. It works best for linear problems, simple physics and geometry. However, it is highly suspicious when complex geometry and flow physics are involved, especially when the subject is out of historical database. The experimental approach is often used to establish database for engineering analysis. However, it is qualitative at best for base flow problems. Other criticisms include the inability to simulate forebody boundary layer correctly, the interference effect from tunnel walls, and the inability to scale all pertinent parameters. Furthermore, there is a contention that the information extrapolated from subscale tests with combustion is not conservative. One potential alternative to the conventional methods is computational fluid dynamics (CFD), which has none of the above restrictions and is becoming more feasible due to maturing algorithms and advancing computer technology. It provides more details of the flowfield and is only limited by computer resources. However, it has its share of criticisms as a predictive tool for base environment. One major concern is that CFD has not been extensively tested for base flow problems. It is therefore imperative that CFD be assessed and benchmarked satisfactorily for base flows. In this study, the turbulent base flowfield of a experimental investigation for a four-engine clustered nozzle is numerically benchmarked using a pressure based CFD method. Since the cold air was the medium, accurate prediction of the base pressure distributions at high altitudes is the primary goal. Other factors which may influence the numerical results such as the effects of grid density, turbulence model, differencing scheme, and boundary conditions are also being addressed. Preliminary results of the computed base pressure agreed reasonably well with that of the measurement. Basic base flow features such as the reverse jet, wall jet, recompression shock, and static pressure field in plane of impingement have been captured.

Wang, Ten-See↗

Longitudinal Laminar Flow Between Cylinders Arranged in Regular Array

The increasing complexity of heat transfer and process situations which involve fluid flow has demanded the frequent use of flow passages of unusual geometrical configuration. The present investigation is concerned with one such novel configuration, namely the longitudinal flow between solid cylindrical rods which are arranged in regular array. A schematic diagram of the situation under study. The rods may be located either in triangular or square array. The flow will be taken to be laminar and fully developed. The aim of this analysis is to determine the pressure drop, shear stress, and velocity-distribution characteristics of the system. The starting point of this study is the basic law of momentum conservation. The resulting differential equation has been solved in an approximate, but almost exact, manner by the use of truncated trigonometric series. Results are obtained over a wide range of porosity values for both the triangular and square arrays. Heat transfer has not been considered. The configuration under investigation has potential application in compact heat exchangers for nuclear reactors and other situations. Further the results should also be of interest in the theory of flow through unconsolidated porous beds (ia, 9a). The only related analytical work known to the authors is that of Emersleben (S), who considered only the square array. His rather involved solution, based on complex zeta functions, appears to be valid only at high porosities. Experiments covering a porosity range of 0.093 to 0.984 have been made by Sullivan (4) using parallel-oriented fibers, most of the tests being for fibers in random array. These previous investigations will be compared with the present theory in a later section.

Sparrow, E. M.↗

Viscous shock-layer flowfield analysis by an explicit-implicit method

The present work extends the recently reported implicit analogue of MacCormack's earlier widely-used explicit method to external axisymmetric laminar flows with strong entropy gradients. The details of the 'numerics' of the implicit part are provided in a body-oriented coordinate system with a moving outer (shock) boundary during the transient part of the solutions. The limiting values of the Courant number are obtained when the shock boundary is treated explicitly. The solution algorithm outlined includes the treatment of the source term associated with the equations in weak conservation form. From the results obtained for two sample problems, it becomes clear that accuracy of predictions is, indeed, very good at higher values of the Courant number. There is a significant saving in overall computing time, depending on the Courant number used and the flow Reynolds number. These properties combined with the simplicity of programming the implicit analogue may appeal to researchers for using it in the analysis of three-dimensional flow problems.

Gupta, R. N.↗

Development of a Regional Groundwater Flow Model for the U.S. Department of Energy Oak Ridge Reservation: Numerical Model using MODFLOW-USG - 17368

A regional scale groundwater flow model (GFM) has been developed for the U.S. Department of Energy (DOE) Oak Ridge Reservation (ORR). Development of the model implements a key recommendation of the ORR Groundwater Strategy. The model will be used as the single calibrated flow model for the ORR and as the framework for future, smaller scale, modeling efforts to support cleanup actions and decisions under the Comprehensive Environmental Response, Compensation, and Liability Act of 1980 (CERCLA). The Oak Ridge site is located in a geologically complex region and encompasses three large government facilities, including the Y-12 National Security Complex (Y-12), Oak Ridge National Laboratory (ORNL), and the East Tennessee Technology Park (ETTP). An interagency approach for addressing legacy groundwater contamination from past operations at these facilities has been developed, resulting in an ORR Groundwater Strategy that was agreed to by DOE, the U.S. Environmental Protection Agency (EPA), and Tennessee Department of Environment and Conservation (TDEC). The ORR Groundwater Strategy provides a comprehensive approach to addressing and prioritizing groundwater issues across the DOE reservation. One component of the strategy is to 'develop and maintain an ORR-wide regional groundwater flow model to ensure a single, regional, calibrated model to support groundwater characterization, decision-making, and remediation.' As part of the cleanup plans for the Oak Ridge site under CERCLA and the Federal Facility Agreement (FFA), there are six final watershed-scale groundwater Records of Decision (RODs) currently planned. In preparation for future projects, the regional GFM has been constructed to serve as the single, calibrated regional flow model to be used as the hydrologic base for the groundwater plume-specific modeling developed for the RODs. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The Space-Time Conservation Element and Solution Element Method-A New High-Resolution and Genuinely Multidimensional Paradigm for Solving Conservation Laws: Numerical Simulation of Shock Waves and Contact Discontinuities - 2

Without resorting to special treatment for each individual test case, the 1D and 2D CE/SE shock-capturing schemes described previously (in Part I) are used to simulate flows involving phenomena such as shock waves, contact discontinuities, expansion waves and their interactions. Five 1D and six 2D problems are considered to examine the capability and robustness of these schemes. Despite their simple logical structures and low computational cost (for the 2D CE/SE shock-capturing scheme, the CPU time is about 2 micro-secs per mesh point per marching step on a Cray C90 machine), the numerical results, when compared with experimental data, exact solutions or numerical solutions by other methods, indicate that these schemes can accurately resolve shock and contact discontinuities consistently.

Wang, Xiao-Yen↗