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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Runge-Kutta methods combined with compact difference schemes for the unsteady Euler equations

Recent development using compact difference schemes to solve the Navier-Stokes equations show spectral-like accuracy. A study was made of the numerical characteristics of various combinations of the Runge-Kutta (RK) methods and compact difference schemes to calculate the unsteady Euler equations. The accuracy of finite difference schemes is assessed based on the evaluations of dissipative error. The objectives are reducing the numerical damping and, at the same time, preserving numerical stability. While this approach has tremendous success solving steady flows, numerical characteristics of unsteady calculations remain largely unclear. For unsteady flows, in addition to the dissipative errors, phase velocity and harmonic content of the numerical results are of concern. As a result of the discretization procedure, the simulated unsteady flow motions actually propagate in a dispersive numerical medium. Consequently, the dispersion characteristics of the numerical schemes which relate the phase velocity and wave number may greatly impact the numerical accuracy. The aim is to assess the numerical accuracy of the simulated results. To this end, the Fourier analysis is to provide the dispersive correlations of various numerical schemes. First, a detailed investigation of the existing RK methods is carried out. A generalized form of an N-step RK method is derived. With this generalized form, the criteria are derived for the three and four-step RK methods to be third and fourth-order time accurate for the non-linear equations, e.g., flow equations. These criteria are then applied to commonly used RK methods such as Jameson's 3-step and 4-step schemes and Wray's algorithm to identify the accuracy of the methods. For the spatial discretization, compact difference schemes are presented. The schemes are formulated in the operator-type to render themselves suitable for the Fourier analyses. The performance of the numerical methods is shown by numerical examples. These examples are detailed. described. The third case is a two-dimensional simulation of a Lamb vortex in an uniform flow. This calculation provides a realistic assessment of various finite difference schemes in terms of the conservation of the vortex strength and the harmonic content after travelling a substantial distance. The numerical implementation of Giles' non-refelctive equations coupled with the characteristic equations as the boundary condition is discussed in detail. Finally, the single vortex calculation is extended to simulate vortex pairing. For the distance between two vortices less than a threshold value, numerical results show crisp resolution of the vortex merging.

Yu, Sheng-Tao↗
NASA NTRS1995

An MHD simulation study of the poloidal mode field line resonance in the Earth's dipole magnetosphere

The poloidal mode field line resonance in the Earth's dipole magnetic field is investigated using cold plasma ideal MHD simulations in dipole geometry. In order to excite the poloidal mode resonance, we use either an initial or a continuous velocity perturbation to drive the system. The perturbation is localized at magnetic shell L = 7 with plasma flow in the radial direction (electric field component in the azimuthal direction). It is found that with the initial perturbation alone, no polodial mode resonance can be obtained and the initially localized perturbation spreads out across all magnetic L shells. With the continuous perturbation, oscillating near the poloidal resonance frequency, a global-scale poloidal cavity mode can be obtained. For the first time, a localized guided poloidal mode resonance is obtained when a radial component of electric field is added to the initial perturbation such that the curl of the electric field is everywhere perpendicular to the background dipole magnetic field. During the localized poloidal resonance, plasma vortices parallel/antiparallel to the background dipole magnetic field B(sub 0). This circular flow, elongated radially, results in twisting of magnetic field flux tubes, which, in turn, leads to the slowdown of the circular plasma flow and reversal of the plasma vortices. The energy associated with the localized poloidal resonance is conserved as it shifts back and forth between the oscillating plasma vortices and the alternately twisted magnetic flux tubes. In the simulations the eigenfunctions associated with the localized poloidal resonance are grid-scale singular functions. This result indicates that ideal MHD is inadequate to describe the underlying problem and nonideal MHD effects are needed for mode broadening.

Ding, D. Q.↗
NASA NTRS1998

A Second Law Based Unstructured Finite Volume Procedure for Generalized Flow Simulation

An unstructured finite volume procedure has been developed for steady and transient thermo-fluid dynamic analysis of fluid systems and components. The procedure is applicable for a flow network consisting of pipes and various fittings where flow is assumed to be one dimensional. It can also be used to simulate flow in a component by modeling a multi-dimensional flow using the same numerical scheme. The flow domain is discretized into a number of interconnected control volumes located arbitrarily in space. The conservation equations for each control volume account for the transport of mass, momentum and entropy from the neighboring control volumes. In addition, they also include the sources of each conserved variable and time dependent terms. The source term of entropy equation contains entropy generation due to heat transfer and fluid friction. Thermodynamic properties are computed from the equation of state of a real fluid. The system of equations is solved by a hybrid numerical method which is a combination of simultaneous Newton-Raphson and successive substitution schemes. The paper also describes the application and verification of the procedure by comparing its predictions with the analytical and numerical solution of several benchmark problems.

Majumdar, Alok↗
NASA NTRS1998

Theory of the Shape of the Heliospheric Termination Shock: Axial Symmetry

The heliospheric termination shock must exhibit asymmetry in its shape, due in part to internal latitude variations in the solar wind, and in part to the special directions defined by the external interstellar flow and/or the galactic magnetic field. This asymmetry shows up naturally in numerical simulations of the interaction between the heliosphere and local interstellar medium. However, to date only one analytical treatment of the theory has appeared, which discussed the modification of the shock shape due to solar wind latitude variations in the presence of spherically symmetric outer boundary conditions. In the present report, we discuss an extension and generalization of this theory to an axially symmetric gasdynamic system, in which departure from spherical symmetry may be due either to internal solar wind variations or to the directional properties of the external interstellar medium. It is shown that for steady flow the post-shock region is characterized by an infinite set of quantities that are conserved along streamlines; among these invariants are the stagnation pressure and a quantity closely related to vorticity. Moreover, for a given latitude profile of the (supersonic) solar wind, the geometry of the termination shock uniquely determines the valued of these invariants at the points where the streamlines emerge from the shock. A self-consistent solution of the problem thus requires a match between the shock geometry and distant boundary conditions. The theory is applied to give estimates of termination-shock asymmetry for various internal and external conditions.

Barnes, Aaron↗
NASA NTRS2002

The Dynamics of Miscible Interfaces: Simulations

The goal of this experimental/computational investigation (joint with Prof Maxworthy at USC) has been to study the dynamics of miscible interfaces, both from a scientific and a practical point of view, and to prepare a related experiment to be flown on the International Space Station. In order to address these effects, we have focused experimental and computational investigations on miscible displacements in cylindrical capillary tubes, as well as in Hele-Shaw cells. Regarding the flow in a capillary tube, the question was addressed as to whether Korteweg stresses and/or divergence effects can potentially account for discrepancies observed between conventional Stokes flow simulations and experiments for miscible flows in capillary tubes. An estimate of the vorticity and streamfunction fields induced by the Kortewegs stresses was derived, which shows these stresses to result in the formation of a vortex ring structure near the tip of the concentration front. Through this mechanism the propagation velocity of the concentration front is reduced, in agreement with the experimental observations. Divergence effects, on the other hand, were seen to be very small, and they have a negligible influence on the tip velocity. As a result, it can be concluded that they are not responsible for the discrepancies between experiments and conventional Stokes simulations. A further part of our investigation focussed on the development of high-accuracy three-dimensional spectral element simulation techniques for miscible flows in capillary tubes, including the effects of variable density and viscosity. Towards this end, the conservation equations are treated in cylindrical coordinates.

Meiburg, Eckart↗
NASA NTRS2011

Performance of the Extravehicular Mobility Unit (EMU) Airlock Coolant Loop Remediation (A/L CLR) Hardware - Final

An EMU water processing kit (Airlock Coolant Loop Recovery -- A/L CLR) was developed as a corrective action to Extravehicular Mobility Unit (EMU) coolant flow disruptions experienced on the International Space Station (ISS) in May of 2004 and thereafter. A conservative duty cycle and set of use parameters for A/L CLR use and component life were initially developed and implemented based on prior analysis results and analytical modeling. Several initiatives were undertaken to optimize the duty cycle and use parameters of the hardware. Examination of post-flight samples and EMU Coolant Loop hardware provided invaluable information on the performance of the A/L CLR and has allowed for an optimization of the process. The intent of this paper is to detail the evolution of the A/L CLR hardware, efforts to optimize the duty cycle and use parameters, and the final recommendations for implementation in the post-Shuttle retirement era.

Steele, John W.↗
NASA NTRS1990

Flowfield computation of entry vehicles

The equations governing the multidimensional flow of a reacting mixture of thermally perfect gasses were derived. The modeling procedures for the various terms of the conservation laws are discussed. A numerical algorithm, based on the finite-volume approach, to solve these conservation equations was developed. The advantages and disadvantages of the present numerical scheme are discussed from the point of view of accuracy, computer time, and memory requirements. A simple one-dimensional model problem was solved to prove the feasibility and accuracy of the algorithm. A computer code implementing the above algorithm was developed and is presently being applied to simple geometries and conditions. Once the code is completely debugged and validated, it will be used to compute the complete unsteady flow field around the Aeroassist Flight Experiment (AFE) body.

Prabhu, Dinesh K.↗
NASA NTRS1987

A new PNS code for three-dimensional chemically reacting flows

A new parabolized Navier-Stokes (PNS) code has been developed to compute the hypersonic, viscous, chemically reacting flow fields around three-dimensional bodies. The flow medium is assumed to be a multicomponent mixture of thermally perfect but calorically imperfect gases. The new PNS code solves the gasdynamic and species conservation equations in a coupled manner using a noniterative, implicit, approximately-factored, finite-difference algorithm. The space-marching method is made well-posed by special treatment of the streamwise pressure gradient term. The code has been used to compute hypersonic laminar flow of chemically reacting air over cones at angles of attack. The results of the computations are compared with the results of reacting boundary-layer computations and show excellent agreement.

Prabhu, D. K.↗
NASA NTRS1986

An analysis of finite-difference and finite-volume formulations of conservation laws

Finite-difference and finite-volume formulations are analyzed in order to clear up the confusion concerning their application to the numerical solution of conservation laws. A new coordinate-free formulation of systems of conservation laws is developed, which clearly distinguishes the role of physical vectors from that of algebraic vectors which characterize the system. The analysis considers general types of equations--potential, Euler, and Navier-Stokes. Three-dimensional unsteady flows with time-varying grids are described using a single, consistent nomeclature for both formulations. Grid motion due to a non-inertial reference frame as well as flow adaptation is covered. In comparing the two formulations, it is found useful to distinguish between differences in numerical methods and differences in grid definition. The former plays a role for non-Cartesian grids, and results in only cosmetic differences in the manner in which geometric terms are handled. The differences in grid definition for the two formulations is found to be more important, since it affects the manner in which boundary conditions, zonal procedures, and grid singularities are handled at computational boundaries. The proper interpretation of strong and weak conservation-law forms for quasi-one-dimensional and axisymmetric flows is brought out.

Vinokur, Marcel↗
NASA NTRS1989

An analysis of finite-difference and finite-volume formulations of conservation laws

Finite-difference and finite-volume formulations are analyzed in order to clear up the confusion concerning their application to the numerical solution of conservation laws. A new coordinate-free formulation of systems of conservation laws is developed, which clearly distinguishes the role of physical vectors from that of algebraic vectors which characterize the system. The analysis considers general types of equations: potential, Euler, and Navier-Stokes. Three-dimensional unsteady flows with time-varying grids are described using a single, consistent nomenclature for both formulations. Grid motion due to a non-inertial reference frame as well as flow adaptation is covered. In comparing the two formulations, it is found useful to distinguish between differences in numerical methods and differences in grid definition. The former plays a role for non-Cartesian grids, and results in only cosmetic differences in the manner in which geometric terms are handled. The differences in grid definition for the two formulations is found to be more important, since it affects the manner in which boundary conditions, zonal procedures, and grid singularities are handled at computational boundaries. The proper interpretation of strong and weak conservation-law forms for quasi-one-dimensional and axisymmetric flows is brought out.

Vinokur, Marcel↗
NASA NTRS1973

Time-split finite-volume method for three-dimensional blunt-body flow

An efficient numerical method for calculating plane, axisymmetric, and fully three-dimensional blunt-body flow is presented. It is a second-order-accurate, time-dependent finite-volume procedure that solves the Euler equations in integral conservation-law form. These equations are written with respect to a Cartesian coordinate system in which an embedded mesh adjusts in time to the motion of the bow shock that is automatically captured as part of the weak solution. With such an adjusting mesh, oscillations in flow properties near the shock are shown to be virtually eliminated. The scheme uses a time-splitting concept that accelerates the convergence appreciably. Comparisons are made between computed and experimental results.

Rizzi, A. W.↗
NASA NTRS1988

Numerical computation of space shuttle orbiter flow field

A new parabolized Navier-Stokes (PNS) code has been developed to compute the hypersonic, viscous chemically reacting flow fields around 3-D bodies. The flow medium is assumed to be a multicomponent mixture of thermally perfect but calorically imperfect gases. The new PNS code solves the gas dynamic and species conservation equations in a coupled manner using a noniterative, implicit, approximately factored, finite difference algorithm. The space-marching method is made well-posed by special treatment of the streamwise pressure gradient term. The code has been used to compute hypersonic laminar flow of chemically reacting air over cones at angle of attack. The results of the computations are compared with the results of reacting boundary-layer computations and show excellent agreement.

Tannehill, John C.↗
NASA NTRS1987

Computation of transonic potential flow on helicopter rotor blades

Two computer codes, the full-potential three-dimensional (FP3D) code and the full-potential rotor (FPR) code have recently been developed. Both of these codes solve the three-dimensional conservative formulation of the full potential equation. The FPR code was developed at the U.S. Army Aeroflightdynamics Directorate (AFDD) while the FP3D code was a joint development by ONERA and AFDD. Both of these codes were used to predict the nonlifting, unsteady flow over a rotor operating at high advance ratio and tip speed. Three different rotor tip planform shapes were studied: a rectangular tip, a 30 deg aft swept tip and a 30 deg forward swept tip. Results of these computations are compared to results obtained using an earlier small-disturbances code. Also, the lifting flow over a rectangular tip operating at a slightly different condition was computed. These results are also compared with the small disturbances computations and with experimental results.

Costes, M.↗
NASA NTRS1972

Numerical computation of multishocked, three-dimensional supersonic flow fields with real gas effects.

A computational procedure is presented which is capable of determining the supersonic flow field surrounding three-dimensional wing-body configurations such as a delta-wing space shuttle. The governing equations in conservation-law form are solved by a finite difference method using a second-order noncentered algorithm between the body and the outermost shock wave, which is treated as a sharp discontinuity. Secondary shocks which form between these boundaries are captured automatically, and the intersection of these shocks with the bow shock posed no difficulty. Resulting flow fields about typical blunt nose shuttle-like configurations at angle of attack are presented. The differences between perfect and real gas effects for high Mach number flows are shown.

Kutler, P.↗
NASA NTRS1993

A multiblock/multigrid Euler method to simulate 2D and 3D compressible flow

An efficient and flexible multiblock/multigrid Euler solver is described which uses a cell-centered finite-volume method based on a multigrid implementation of the alternating direction implicit algorithm. The fully conservative treatment of the inter-block boundary allows the passage of discontinuities across block boundaries with minimum distortion of the solution. Information is exchanged between blocks by using surface arrays which contain all the data needed to update the inter-block boundary conditions. The proposed multiblock/multigrid approach is considered to be feasible for solving flow problems involving complex geometries.

Wang, Lixia↗
NASA NTRS1999

Semiempirical Two-Dimensional Magnetohydrodynamic Model of the Solar Corona and Interplanetary Medium

We have developed a two-dimensional semiempirical MHD model of the solar corona and solar wind. The model uses empirically derived electron density profiles from white-light coronagraph data measured during the Skylub period and an empirically derived model of the magnetic field which is fitted to observed streamer topologies, which also come from the white-light coronagraph data The electron density model comes from that developed by Guhathakurta and coworkers. The electron density model is extended into interplanetary space by using electron densities derived from the Ulysses plasma instrument. The model also requires an estimate of the solar wind velocity as a function of heliographic latitude and radial component of the magnetic field at 1 AU, both of which can be provided by the Ulysses spacecraft. The model makes estimates as a function of radial distance and latitude of various fluid parameters of the plasma such as flow velocity V, effective temperature T(sub eff), and effective heat flux q(sub eff), which are derived from the equations of conservation of mass, momentum, and energy, respectively. The term effective indicates that wave contributions could be present. The model naturally provides the spiral pattern of the magnetic field far from the Sun and an estimate of the large-scale surface magnetic field at the Sun, which we estimate to be approx. 12 - 15 G. The magnetic field model shows that the large-scale surface magnetic field is dominated by an octupole term. The model is a steady state calculation which makes the assumption of azimuthal symmetry and solves the various conservation equations in the rotating frame of the Sun. The conservation equations are integrated along the magnetic field direction in the rotating frame of the Sun, thus providing a nearly self-consistent calculation of the fluid parameters. The model makes a minimum number of assumptions about the physics of the solar corona and solar wind and should provide a very accurate empirical description of the solar corona and solar wind Once estimates of mass density rho, flow velocity V, effective temperature T(sub eff), effective heat flux q(sub eff), and magnetic field B are computed from the model and waves are assumed unimportant, all other plasma parameters such as Mach number, Alfven speed, gyrofrequency, etc. can be derived as a function of radial distance and latitude from the Sun. The model can be used as a planning tool for such missions as Slar Probe and provide an empirical framework for theoretical models of the solar corona and solar wind The model will be used to construct a semiempirical MHD description of the steady state solar corona and solar wind using the SOHO Large Angle Spectrometric Coronagraph (LASCO) polarized brightness white-light coronagraph data, SOHO Extreme Ultraviolet Imaging Telescope data, and Ulysses plasma data.

Sittler, Edward C., Jr.↗
NASA NTRS1993

A time-accurate finite volume method valid at all flow velocities

A finite volume method to solve the Navier-Stokes equations at all flow velocities (e.g., incompressible, subsonic, transonic, supersonic and hypersonic flows) is presented. The numerical method is based on a finite volume method that incorporates a pressure-staggered mesh and an incremental pressure equation for the conservation of mass. Comparison of three generally accepted time-advancing schemes, i.e., Simplified Marker-and-Cell (SMAC), Pressure-Implicit-Splitting of Operators (PISO), and Iterative-Time-Advancing (ITA) scheme, are made by solving a lid-driven polar cavity flow and self-sustained oscillatory flows over circular and square cylinders. Calculated results show that the ITA is the most stable numerically and yields the most accurate results. The SMAC is the most efficient computationally and is as stable as the ITA. It is shown that the PISO is the most weakly convergent and it exhibits an undesirable strong dependence on the time-step size. The degenerated numerical results obtained using the PISO are attributed to its second corrector step that cause the numerical results to deviate further from a divergence free velocity field. The accurate numerical results obtained using the ITA is attributed to its capability to resolve the nonlinearity of the Navier-Stokes equations. The present numerical method that incorporates the ITA is used to solve an unsteady transitional flow over an oscillating airfoil and a chemically reacting flow of hydrogen in a vitiated supersonic airstream. The turbulence fields in these flow cases are described using multiple-time-scale turbulence equations. For the unsteady transitional over an oscillating airfoil, the fluid flow is described using ensemble-averaged Navier-Stokes equations defined on the Lagrangian-Eulerian coordinates. It is shown that the numerical method successfully predicts the large dynamic stall vortex (DSV) and the trailing edge vortex (TEV) that are periodically generated by the oscillating airfoil. The calculated streaklines are in very good comparison with the experimentally obtained smoke picture. The calculated turbulent viscosity contours show that the transition from laminar to turbulent state and the relaminarization occur widely in space as well as in time. The ensemble-averaged velocity profiles are also in good agreement with the measured data and the good comparison indicates that the numerical method as well as the multipletime-scale turbulence equations successfully predict the unsteady transitional turbulence field. The chemical reactions for the hydrogen in the vitiated supersonic airstream are described using 9 chemical species and 48 reaction-steps. Consider that a fast chemistry can not be used to describe the fine details (such as the instability) of chemically reacting flows while a reduced chemical kinetics can not be used confidently due to the uncertainty contained in the reaction mechanisms. However, the use of a detailed finite rate chemistry may make it difficult to obtain a fully converged solution due to the coupling between the large number of flow, turbulence, and chemical equations. The numerical results obtained in the present study are in good agreement with the measured data. The good comparison is attributed to the numerical method that can yield strongly converged results for the reacting flow and to the use of the multiple-time-scale turbulence equations that can accurately describe the mixing of the fuel and the oxidant.

Kim, S.-W.↗
NASA NTRS1977

Altitude engine test of a turbofan exhaust gas mixer to conserve fuel

A comparison of the specific fuel consumption was made with and without an internal mixer installed in a low bypass ratio, confluent flow turbofan engine. Tests were conducted at several Mach numbers and altitudes for core to fan stream total temperature ratios of 2.0 and 2.5 and mixing lengths of L/D = 0.95 and 1.74. For these test conditions, the specific fuel consumption improvement varied from 2.5 to 4.0 percent.

Cullom, R. R.↗