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.

SourceResearch connectionFrom yearThrough year

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

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

NASA NTRS1995

Numerical modeling of spray combustion with an advanced VOF method

This paper summarizes the technical development and validation of a multiphase computational fluid dynamics (CFD) numerical method using the volume-of-fluid (VOF) model and a Lagrangian tracking model which can be employed to analyze general multiphase flow problems with free surface mechanism. The gas-liquid interface mass, momentum and energy conservation relationships are modeled by continuum surface mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed flow regimes. The objectives of the present study are to develop and verify the fractional volume-of-fluid cell partitioning approach into a predictor-corrector algorithm and to demonstrate the effectiveness of the present approach by simulating benchmark problems including laminar impinging jets, shear coaxial jet atomization and shear coaxial spray combustion flows.

Chen, Yen-Sen↗

NASA NTRS2003

Robust and Simple Non-Reflecting Boundary Conditions for the Euler Equations: A New Approach Based on the Space-Time CE/SE Method

This paper reports on a significant advance in the area of non-reflecting boundary conditions (NRBCs) for unsteady flow computations. As a part of the development of the space-time conservation element and solution element (CE/SE) method, sets of NRBCs for 1D Euler problems are developed without using any characteristics-based techniques. These conditions are much simpler than those commonly reported in the literature, yet so robust that they are applicable to subsonic, transonic and supersonic flows even in the presence of discontinuities. In addition, the straightforward multidimensional extensions of the present 1D NRBCs have been shown numerically to be equally simple and robust. The paper details the theoretical underpinning of these NRBCs, and explains their unique robustness and accuracy in terms of the conservation of space-time fluxes. Some numerical results for an extended Sod's shock-tube problem, illustrating the effectiveness of the present NRBCs are included, together with an associated simple Fortran computer program. As a preliminary to the present development, a review of the basic CE/SE schemes is also included.

Chang, Sin-Chung↗

NASA NTRS1991

Grid generation about complex three-dimensional aircraft configurations

The problem of obtaining three dimensional grids with sufficient resolution to resolve all the flow or other physical features of interest is addressed. The generation of a computational grid involves a series of compromises to resolve several conflicting requirements. On one hand, one would like the grid to be fine enough and not too skewed to reduce the numerical errors and to adequately resolve the pertinent physical features of the flow field about the aircraft. On the other hand, the capabilities of present or even future supercomputers are finite and the number of mesh points must be limited to a reasonable number: one which is usually much less than desired for numerical accuracy. One technique to overcome this limitation is the 'zonal' grid approach. In this method, the overall field is subdivided into smaller zones or blocks in each of which an independent grid is generated with enough grid density to resolve the flow features in that zone. The zonal boundaries or interfaces require special boundary conditions such that the conservation properties of the governing equations are observed. Much work was done in 3-D zonal approaches with nonconservative zonal interfaces. A 3-D zonal conservative interfacing method that is efficient and easy to implement was developed during the past year. During the course of the work, it became apparent that it would be much more feasible to do the conservative interfacing with cell-centered finite volume codes instead of the originally planned finite difference codes. Accordingly, the CNS code was converted to finite volume form. This new version of the code is named CNSFV. The original multi-zonal interfacing capability of the CNS code was enhanced by generalizing the procedure to allow for completely arbitrarily shaped zones with no mesh continuity between the zones. While this zoning capability works well for most flow situations, it is, however, still nonconservative. The conservative interface algorithm was also implemented but was not completely validated.

Klopfer, Goetz H.↗

NASA NTRS1988

Circumferential flow analysis at the aft field joint of the Space Shuttle solid rocket motor

Flow analyses have been performed to determine the nature of the three-dimensional flow field in the vicinity of the aft-most field joint of the Space Shuttle Redesigned Solid Rocket Motor (RSRM). Specific objectives included the quantification of the circumferential pressure and velocity gradients at the joint location which might be caused by the non-uniform erosion of the rubber inhibitor which protrudes from the wall into the flow field. Three-dimensional Navier-Stokes equations have been solved in conjunction with the conservation equation for the turbulence energy and the dissipation rate. The numerical predictions have been compared with the measurements from a 7.5 percent scale cold flow model of the redesigned solid rocket motor.

Majumdar, Alok K.↗

NASA NTRS1979

Axisymmetric internal flows with shocks

Internal flow problems with supersonic entrance conditions and subsonic exit conditions are studied. The suitability of applying the transonic-small disturbance theory to internal flows is examined and found to be very limited. The full inviscid equations in nonconservative form are solved by relaxation. Jump conditions which conserve mass and normal momentum are applied explicitly at shock waves. This method is suitable for solving flows which include supersonic-subsonic shock waves nearly normal to the streamlines.

Rizk, M. H.↗

NASA NTRS1992

Low-to-high altitude predictions of three-dimensional ablative reentry flowfields

A new coupled-chemistry PNS scheme has been developed to study 3D hypersonic reentry flows over sphere cones at small angles of attack. A new quasi-steady ablation model for Teflon was used to predict ablative reentry flow fields under laminar, transitional, and turbulent conditions. This solution scheme uses a new two-step solution scheme for solving the species conservation equations under low-altitude conditions. As an example, the scheme is used to study hypersonic flow over a 6-deg sphere-cone configuration at 50 and 125 kft conditions under zero and nonzero angle of attack conditions.

Bhutta, Bilal A.↗

NASA NTRS1988

Extensions of the CSCM methodology for nonequilibrium reacting gas flows

The paper describes the extension of the CSCM flux difference splitting to produce conservative unconditionally stable implicit upwind methods for the extended equation set. The CSCM-S single data level relaxation algorithm is argued to be particularly advantageous for efficiently solving the set of equations. Attention is given to issues of alternative state variable representations and to complementary thermodynamic relations that parameterize and need to be effectively coupled with the solution procedure for the partial differential equations. Within the framework of iteration, consideration is given to more efficient approximate solution procedures than the full block coupled LU decomposition. Results of early numerical experiments are given to illustrate the developments.

Lombard, C. K.↗

NASA NTRS1987

Evaluation of a research circulation control airfoil using Navier-Stokes methods

The compressible Reynolds time averaged Navier-Stokes equations were used to obtain solutions for flows about a two dimensional circulation control airfoil. The governing equations were written in conservation form for a body-fitted coordinate system and solved using an Alternating Direction Implicit (ADI) procedure. A modified algebraic eddy viscosity model was used to define the turbulent characteristics of the flow, including the wall jet flow over the Coanda surface at the trailing edge. Numerical results are compared to experimental data obtained for a research circulation control airfoil geometry. Excellent agreement with the experimental results was obtained.

Shrewsbury, George D.↗

NASA NTRSDate not supplied

MEDLI2: MISP Inferred Aerothermal Environment and Flow Transition Assessment

The Mars Entry, Descent, and Landing Instrumentation 2 (MEDLI2) sensor suite on the Mars2020 mission contained multiple sensors on the aeroshell to measure the aerothermal environment during entry into the Martian atmosphere. These sensors performed superbly and successfully returned forebody and aftbody heating measurements. Analysis of MEDLI2 data indicated flow transitioning from a laminar to turbulent state on the heatshield. No evidence of flow transition was observed on the backshell. Two methods were used to estimate flow transition times on the heatshield: (1) temperature gradient of near-surface thermocouple data and (2) heat flux gradient from an inverse reconstruction approach using thermocouple data and material response modeling. Both methods produced similar transition times with an estimated accuracy of ±1 s. To assess various transition criteria, transition parameters were evaluated at each sensor location using flow field solutions from computational fluid dynamics (CFD) simulations. The idea was to use conservative values inferred from MEDLI2 data as transition criteria for other Mars missions. To test this hypothesis, MEDLI data from the Mars Science Laboratory (MSL) mission was used to compare predicted vs. actual flow transition times. The comparisons suggest smooth wall transition criteria are not well-suited in modeling the rapid progression of a turbulent transition front. Transition criteria containing a roughness element parameter agreed better with the flight data. In summary, critical transition values derived from MEDLI2 data may be used as a starting point in constructing a flow transition model for future Mars missions.

Chun Y Tang↗

NASA NTRS1994

Prediction of Undsteady Flows in Turbomachinery Using the Linearized Euler Equations on Deforming Grids

A linearized Euler solver for calculating unsteady flows in turbomachinery blade rows due to both incident gusts and blade motion is presented. The model accounts for blade loading, blade geometry, shock motion, and wake motion. Assuming that the unsteadiness in the flow is small relative to the nonlinear mean solution, the unsteady Euler equations can be linearized about the mean flow. This yields a set of linear variable coefficient equations that describe the small amplitude harmonic motion of the fluid. These linear equations are then discretized on a computational grid and solved using standard numerical techniques. For transonic flows, however, one must use a linear discretization which is a conservative linearization of the non-linear discretized Euler equations to ensure that shock impulse loads are accurately captured. Other important features of this analysis include a continuously deforming grid which eliminates extrapolation errors and hence, increases accuracy, and a new numerically exact, nonreflecting far-field boundary condition treatment based on an eigenanalysis of the discretized equations. Computational results are presented which demonstrate the computational accuracy and efficiency of the method and demonstrate the effectiveness of the deforming grid, far-field nonreflecting boundary conditions, and shock capturing techniques. A comparison of the present unsteady flow predictions to other numerical, semi-analytical, and experimental methods shows excellent agreement. In addition, the linearized Euler method presented requires one or two orders-of-magnitude less computational time than traditional time marching techniques making the present method a viable design tool for aeroelastic analyses.

Clark, William S.↗

NASA NTRS1988

Temporal and spatial adaptive algorithm for reacting flows

A numerical integration scheme for quasi-one-dimensional unsteady flows with finite-rate chemistry is developed and demonstrated. The governing Euler and species-conservation equations are derived; the integration method and the spatial and temporal grid embedding techniques are explained; and results for sample problems involving stream-tube flow with one dissociating gas, shock-tube flow with one dissociating gas, and diverging-channel flow with multiple reactions are presented in extensive graphs and briefly characterized. Accuracy comparable to that of globally fine grid solutions is obtained with significant savings in CPU time.

Pervaiz, Mehtab M.↗

NASA NTRS2018

Tetrahedral-Mesh Simulations of Shock-Turbulence Interaction

Despite decades of development of unstructured mesh methods, direct numerical simulations (DNS) of turbulent flows are still predominantly performed on structured or unstructured hexahedral meshes with high-order finite-difference methods, weighted essentially nonoscillatory (WENO) schemes, or hybrid schemes formed by their combinations. Tetrahedral meshes offer easy mesh generation and adaptation around complex geometries and the potential of an orientation-free grid that would benefit the isotropic nature of small-scale dissipation, as well as the solution accuracy of intermediate scales. To advance the state of the art of unstructured-mesh simulation capabilities for shock/turbulence interaction, DNS using pure tetrahedral meshes are carried out with the space-time conservation element, solution element (CESE) method in this research. By its design, the CESE method is constructed based on a non-dissipative scheme and is a genuinely multidimensional numerical framework that is free from the use of an approximate Riemann-solver. The numerical framework also provides the ability to add numerical dissipation (the nondissipative scheme acts as the reference state like that of the reversible state in thermodynamics) when needed (with justification from mathematics/physics). The above-mentioned features along with the CESE method's consistent shock-capturing approach and strong enforcement of flux conservation in spacetime offers a novel method to accurately simulate turbulent flows and their interaction with shocks using tetrahedral meshes. Two canonical problems, namely, isotropic turbulence interaction with a normal shock and a Mach 2.9 turbulent boundary layer flow over a 24deg compression corner are investigated in this study. Computational results show reasonably good agreement with experimental data and results from structured-mesh, high-order simulations available in the literature. Successful validation of these canonical problems demonstrated here paves the way for future high-fidelity supersonic flow simulations involving complex-geometries.

Venkatachari, Balaji Shankar↗

NASA NTRS1989

Dynamic grid deformation using Navier-displacement equation for deforming wings

For dynamic and aeroelastic applications of maneuvering wings, the solid boundaries undergo rigid-body motion and aeroelastic deformation. For rigid-body motion, the conservative fluid dynamics equation, in terms of the Eulerian description, is written in terms of a moving frame of reference, and the problem is solved on a time-independent body-conformed grid. For both rigid-body motion and aeroelastic deformation, the Navier-displacement equation, in terms of the Lagrangian coordinates, is modified for fluid-flow problems. It is used along with the Eulerian description of the conservative fluid dynamics equations to account for the grid deformation.

Kandil, Osama A.↗

NASA NTRS1986

Analysis of Scramjet Inlets

NASCRIN analyzes two-dimensional flow fields in supersonic combustion ramjet (scramjet) inlets. Solves two-dimensional Euler or Navier-Strokes equations in conservative form by unsplit, explicit, two-step finitedifference method. More recent explicit/implicit, two-step scheme also incorporated for viscous flow analysis. Algebraic, two-layer eddy-viscosity model used for turbulent flow calculations.

Kumar, A.↗

NASA NTRS1993

Simplified, inverse, ejector design tool

A simple lumped parameter based inverse design tool has been developed which provides flow path geometry and entrainment estimates subject to operational, acoustic, and design constraints. These constraints are manifested through specification of primary mass flow rate or ejector thrust, fully-mixed exit velocity, and static pressure matching. Fundamentally, integral forms of the conservation equations coupled with the specified design constraints are combined to yield an easily invertible linear system in terms of the flow path cross-sectional areas. Entrainment is computed by back substitution. Initial comparison with experimental and analogous one-dimensional methods show good agreement. Thus, this simple inverse design code provides an analytically based, preliminary design tool with direct application to High Speed Civil Transport (HSCT) design studies.

Dechant, Lawrence J.↗

NASA NTRS1974

A comparison of numerical solutions of the advective equation

Second- and third-order finite-difference methods recently applied to problems in high-speed fluid flow are applied to the model advection equation cast in conservative form. The differencing methods considered use forward time differencing with both centered and preferential space differences employed in the predictor-corrector sequences. The free parameter required for stability in the third-order method is adjusted to cause the solution to be either minimum dispersive or minimum dissipative in nature. Results indicate that the third-order method using minimum dispersion is the most accurate method tested. Computer time requirements are approximately twice those needed for second-order techniques.

Anderson, D.↗

NASA NTRS1976

Hybrid finite element methods

The purpose of this paper is to show how Lagrange multipliers can be used with finite elements to achieve a number of desirable properties in the underlying approximation. For elliptic boundary value problems, variational principles can be developed in which all boundary conditions are natural. In fluid flow problems, one can endow the approximations with physically essential conservation laws.

Fix, G. M.↗

NASA NTRS1976

Diffraction of a shock wave by a compression corner. I - Regular reflection

The unsteady, two-dimensional flowfield resulting from the interaction of a moving planar shock wave with a compression corner is determined using a second-order, discontinuity-fitting, finite-difference approach. The time-dependent Euler equations are transformed to normalize the distance between the body and peripheral shock and to include the existing self-similar property of the flow. The resulting set of partial differential equations in conservation-law form is then solved in a time-dependent fashion using MacCormack's scheme. The vortical singularity, which lies on the body surface, and the single reflected shock are both treated as discontinuities in the numerical procedure. The results of the numerical simulation compare quite favorably with existing experimental interferograms and yield better flowfield resolution than previous first-order, shock-capturing, numerical solutions.

Kutler, P.↗