Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Primitive equations”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

A defect stream function, law of the wall/wake method for compressible turbulent boundary layers

The application of the defect stream function to the solution of the two-dimensional, compressible boundary layer is examined. A law of the wall/law of the wake formulation for the inner part of the boundary layer is presented which greatly simplifies the computational task near the wall and eliminates the need for an eddy viscosity model in this region. The eddy viscosity model in the outer region is arbitrary. The modified Crocco temperature-velocity relationship is used as a simplification of the differential energy equation. Formulations for both equilibrium and nonequilibrium boundary layers are presented including a constrained zero-order form which significantly reduces the computational workload while retaining the significant physics of the flow. A formulation for primitive variables is also presented. Results are given for the constrained zero-order and second-order equilibrium formulations and are compared with experimental data. A compressible wake function valid near the wall has been developed from the present results.

Barnwell, Richard W.↗
NASA NTRS1994

Three-Dimensional Navier-Stokes Method with Two-Equation Turbulence Models for Efficient Numerical Simulation of Hypersonic Flows

A new computational efficient 3-D compressible Reynolds-averaged implicit Navier-Stokes method with advanced two equation turbulence models for high speed flows is presented. All convective terms are modeled using an entropy satisfying higher-order Total Variation Diminishing (TVD) scheme based on implicit upwind flux-difference split approximations and arithmetic averaging procedure of primitive variables. This method combines the best features of data management and computational efficiency of space marching procedures with the generality and stability of time dependent Navier-Stokes procedures to solve flows with mixed supersonic and subsonic zones, including streamwise separated flows. Its robust stability derives from a combination of conservative implicit upwind flux-difference splitting with Roe's property U to provide accurate shock capturing capability that non-conservative schemes do not guarantee, alternating symmetric Gauss-Seidel 'method of planes' relaxation procedure coupled with a three-dimensional two-factor diagonal-dominant approximate factorization scheme, TVD flux limiters of higher-order flux differences satisfying realizability, and well-posed characteristic-based implicit boundary-point a'pproximations consistent with the local characteristics domain of dependence. The efficiency of the method is highly increased with Newton Raphson acceleration which allows convergence in essentially one forward sweep for supersonic flows. The method is verified by comparing with experiment and other Navier-Stokes methods. Here, results of adiabatic and cooled flat plate flows, compression corner flow, and 3-D hypersonic shock-wave/turbulent boundary layer interaction flows are presented. The robust 3-D method achieves a better computational efficiency of at least one order of magnitude over the CNS Navier-Stokes code. It provides cost-effective aerodynamic predictions in agreement with experiment, and the capability of predicting complex flow structures in complex geometries with good accuracy.

Bardina, J. E.↗
NASA NTRS1990

Finite-volume application of high order ENO schemes to multi-dimensional boundary-value problems

The finite volume approach in developing multi-dimensional, high-order accurate essentially non-oscillatory (ENO) schemes is considered. In particular, a two dimensional extension is proposed for the Euler equation of gas dynamics. This requires a spatial reconstruction operator that attains formal high order of accuracy in two dimensions by taking account of cross gradients. Given a set of cell averages in two spatial variables, polynomial interpolation of a two dimensional primitive function is employed in order to extract high-order pointwise values on cell interfaces. These points are appropriately chosen so that correspondingly high-order flux integrals are obtained through each interface by quadrature, at each point having calculated a flux contribution in an upwind fashion. The solution-in-the-small of Riemann's initial value problem (IVP) that is required for this pointwise flux computation is achieved using Roe's approximate Riemann solver. Issues to be considered in this two dimensional extension include the implementation of boundary conditions and application to general curvilinear coordinates. Results of numerical experiments are presented for qualitative and quantitative examination. These results contain the first successful application of ENO schemes to boundary value problems with solid walls.

Casper, Jay↗
NASA NTRS1990

A defect stream function formulation for compressible turbulent boundary layers

Progress to date on the development of a method for turbulent, wall-bounded flow which uses the defect stream function formulation in the outer layer and an analytic law of the wall and wake formulation in the inner region is reviewed. This two-formulation approach avoids the need to computationally resolve the high-gradient inner layer. One of the most appealing recent developments is the transformation of the compressible governing equation for the defect stream function into a linear, second-order differential equation which has analytic solutions for many problems of practical interest. Numerical and analytic results for incompressible and compressible flows are shown to be in excellent agreement with experimental results. In this paper the two-formulation approach is applied to primitive-variable computations. Excellent comparisons with experiment are presented for two compressible flat plate flows.

Barnwell, Richard W.↗
NASA NTRS1998

A Semantic Analysis Method for Scientific and Engineering Code

This paper develops a procedure to statically analyze aspects of the meaning or semantics of scientific and engineering code. The analysis involves adding semantic declarations to a user's code and parsing this semantic knowledge with the original code using multiple expert parsers. These semantic parsers are designed to recognize formulae in different disciplines including physical and mathematical formulae and geometrical position in a numerical scheme. In practice, a user would submit code with semantic declarations of primitive variables to the analysis procedure, and its semantic parsers would automatically recognize and document some static, semantic concepts and locate some program semantic errors. A prototype implementation of this analysis procedure is demonstrated. Further, the relationship between the fundamental algebraic manipulations of equations and the parsing of expressions is explained. This ability to locate some semantic errors and document semantic concepts in scientific and engineering code should reduce the time, risk, and effort of developing and using these codes.

Stewart, Mark E. M.↗
NASA NTRS1995

Probabilistic simulation of long term behavior in polymer matrix composites

A methodology to compute cumulative probability distribution functions (CDF) of fatigue life for different ratios, r of applied stress to the laminate strength based on first ply failure criteria has been developed and demonstrated. Degradation effects due to long term environmental exposure and mechanical cyclic loads are considered in the simulation process. A unified time-stress dependent multi-factor interaction equation model developed at NASA Lewis Research Center has been used to account for the degradation/aging of material properties due to cyclic loads. Fast probability integration method is used to perform probabilistic simulation of uncertainties. Sensitivity of fatigue life reliability to uncertainties in the primitive random variables are computed and their significance in the reliability based design for maximum life is discussed. The results show that the graphite/epoxy (0/+45/90) deg laminate with ply thickness 0.125 in. has 500,000 cycles life for applied stress to laminate strength ratio of 0.6 and a reliability of 0.999. Also, the fatigue life reliability has been found to be most sensitive to the ply thickness and matrix tensile strength. Tighter quality controls must therefore be enforced on ply thickness and matrix strength in order to achieve high reliability of the structure.

Shah, A. R.↗
NASA NTRS1993

Segmented multigrid domain decomposition solutions for three dimensional viscous recirculating flows

A segmented multigrid domain decomposition strategy is combined with a pressure-based form of flux-vector discretization for 3D incompressible and compressible viscous flow applications. A pressure-based form of flux-vector splitting is applied to the Navier-Stokes (NS) equations, which are represented by an implicit lowest-order reduced NS system and a purely diffusive higher-order deferred corrector. A trapezoidal or boxlike form of discretization insures that all mass conservation properties are satisfied at interfacial and outflow boundaries, even for this primitive-variable nonstaggered grid computation. Improvements in gridding strategy are presented by allowing for disjoint subdomains that provide optimal resolution of disparate flow features.

Srinivasan, Kumar↗
NASA NTRS1993

Probabilistic simulation of the human factor in structural reliability

A formal approach is described in an attempt to computationally simulate the probable ranges of uncertainties of the human factor in structural probabilistic assessments. A multi-factor interaction equation (MFIE) model has been adopted for this purpose. Human factors such as marital status, professional status, home life, job satisfaction, work load and health, are considered to demonstrate the concept. Parametric studies in conjunction with judgment are used to select reasonable values for the participating factors (primitive variables). Suitability of the MFIE in the subsequently probabilistic sensitivity studies are performed to assess the validity of the whole approach. Results obtained show that the uncertainties for no error range from five to thirty percent for the most optimistic case.

Chamis, C. C.↗
NASA NTRS1990

On computations of the integrated space shuttle flowfield using overset grids

Numerical simulations using the thin-layer Navier-Stokes equations and chimera (overset) grid approach were carried out for flows around the integrated space shuttle vehicle over a range of Mach numbers. Body-conforming grids were used for all the component grids. Testcases include a three-component overset grid - the external tank (ET), the solid rocket booster (SRB) and the orbiter (ORB), and a five-component overset grid - the ET, SRB, ORB, forward and aft attach hardware, configurations. The results were compared with the wind tunnel and flight data. In addition, a Poisson solution procedure (a special case of the vorticity-velocity formulation) using primitive variables was developed to solve three-dimensional, irrotational, inviscid flows for single as well as overset grids. The solutions were validated by comparisons with other analytical or numerical solution, and/or experimental results for various geometries. The Poisson solution was also used as an initial guess for the thin-layer Navier-Stokes solution procedure to improve the efficiency of the numerical flow simulations. It was found that this approach resulted in roughly a 30 percent CPU time savings as compared with the procedure solving the thin-layer Navier-Stokes equations from a uniform free stream flowfield.

Chiu, I-T.↗
NASA NTRS1987

The role of turbulent convection in the primitive solar nebula. I - Theory. II - Results

A model of convective turbulence which takes radiative dissipation, rotation, and convective motion anisotropy into account, on the basis of a closure for the nonlinear interactions that employs the growth rates of hydrodynamic instabilities, is used to obtain a theoretical framework for modeling the primordial solar nebula. It is assumed that convection is the sole source of turbulence causing the solar nebula to evolve. Vertical structure equations in the thin disk approximation are developed and a detailed comparison with the previous solar nebula convective models of such workers as Lin et al. (1981, 1982) is undertaken. The present values for the turbulent efficiency are much lower and more sensitive to opacity and surface density, resulting in low turbulent speeds, a more massive disk, a lower accretion rate 'best value', and a longer characteristic dispersal time for the disk. It is concluded that convection may not be the dominant source of turbulence needed to evolve young solar/stellar nebulae.

Cabot, W.↗
NASA NTRS2010

Physics of Life: A Model for Non-Newtonian Properties of Living Systems

This innovation proposes the reconciliation of the evolution of life with the second law of thermodynamics via the introduction of the First Principle for modeling behavior of living systems. The structure of the model is quantum-inspired: it acquires the topology of the Madelung equation in which the quantum potential is replaced with the information potential. As a result, the model captures the most fundamental property of life: the progressive evolution; i.e. the ability to evolve from disorder to order without any external interference. The mathematical structure of the model can be obtained from the Newtonian equations of motion (representing the motor dynamics) coupled with the corresponding Liouville equation (representing the mental dynamics) via information forces. All these specific non-Newtonian properties equip the model with the levels of complexity that matches the complexity of life, and that makes the model applicable for description of behaviors of ecological, social, and economical systems. Rather than addressing the six aspects of life (organization, metabolism, growth, adaptation, response to stimuli, and reproduction), this work focuses only on biosignature ; i.e. the mechanical invariants of life, and in particular, the geometry and kinematics of behavior of living things. Living things obey the First Principles of Newtonian mechanics. One main objective of this model is to extend the First Principles of classical physics to include phenomenological behavior on living systems; to develop a new mathematical formalism within the framework of classical dynamics that would allow one to capture the specific properties of natural or artificial living systems such as formation of the collective mind based upon abstract images of the selves and non-selves; exploitation of this collective mind for communications and predictions of future expected characteristics of evolution; and for making decisions and implementing the corresponding corrections if the expected scenario is different from the originally planned one. This approach postulates that even a primitive living species possesses additional, non-Newtonian properties that are not included in the laws of Newtonian or statistical mechanics. These properties follow from a privileged ability of living systems to possess a self-image (a concept introduced in psychology) and to interact with it. The proposed mathematical system is based on the coupling of the classical dynamical system representing the motor dynamics with the corresponding Liouville equation describing the evolution of initial uncertainties in terms of the probability density and representing the mental dynamics. The coupling is implemented by the information-based supervising forces that can be associated with self-awareness. These forces fundamentally change the pattern of the probability evolution, and therefore, lead to a major departure of the behavior of living systems from the patterns of both Newtonian and statistical mechanics. This innovation is meant to capture the signature of life based only on observable behavior, not on any biochemistry. This will not prevent the use of this model for developing artificial living systems, as well as for studying some general properties of behavior of natural, living systems.

Zak, Michail↗
NASA NTRS1992

Adaptive multigrid domain decomposition solutions for viscous interacting flows

Several viscous incompressible flows with strong pressure interaction and/or axial flow reversal are considered with an adaptive multigrid domain decomposition procedure. Specific examples include the triple deck structure surrounding the trailing edge of a flat plate, the flow recirculation in a trough geometry, and the flow in a rearward facing step channel. For the latter case, there are multiple recirculation zones, of different character, for laminar and turbulent flow conditions. A pressure-based form of flux-vector splitting is applied to the Navier-Stokes equations, which are represented by an implicit lowest-order reduced Navier-Stokes (RNS) system and a purely diffusive, higher-order, deferred-corrector. A trapezoidal or box-like form of discretization insures that all mass conservation properties are satisfied at interfacial and outflow boundaries, even for this primitive-variable, non-staggered grid computation.

Rubin, Stanley G.↗
NASA NTRS1993

Two-dimensional CFD modeling of wave rotor flow dynamics

A two-dimensional Navier-Stokes solver developed for detailed study of wave rotor flow dynamics is described. The CFD model is helping characterize important loss mechanisms within the wave rotor. The wave rotor stationary ports and the moving rotor passages are resolved on multiple computational grid blocks. The finite-volume form of the thin-layer Navier-Stokes equations with laminar viscosity are integrated in time using a four-stage Runge-Kutta scheme. The Roe approximate Riemann solution scheme or the computationally less expensive Advection Upstream Splitting Method (AUSM) flux-splitting scheme are used to effect upwind-differencing of the inviscid flux terms, using cell interface primitive variables set by MUSCL-type interpolation. The diffusion terms are central-differenced. The solver is validated using a steady shock/laminar boundary layer interaction problem and an unsteady, inviscid wave rotor passage gradual opening problem. A model inlet port/passage charging problem is simulated and key features of the unsteady wave rotor flow field are identified. Lastly, the medium pressure inlet port and high pressure outlet port portion of the NASA Lewis Research Center experimental divider cycle is simulated and computed results are compared with experimental measurements. The model accurately predicts the wave timing within the rotor passage and the distribution of flow variables in the stationary inlet port region.

Welch, Gerard E.↗
NASA NTRS1994

Two-dimensional CFD modeling of wave rotor flow dynamics

A two-dimensional Navier-Stokes solver developed for detailed study of wave rotor flow dynamics is described. The CFD model is helping characterize important loss mechanisms within the wave rotor. The wave rotor stationary ports and the moving rotor passages are resolved on multiple computational grid blocks. The finite-volume form of the thin-layer Navier-Stokes equations with laminar viscosity are integrated in time using a four-stage Runge-Kutta scheme. Roe's approximate Riemann solution scheme or the computationally less expensive advection upstream splitting method (AUSM) flux-splitting scheme is used to effect upwind-differencing of the inviscid flux terms, using cell interface primitive variables set by MUSCL-type interpolation. The diffusion terms are central-differenced. The solver is validated using a steady shock/laminar boundary layer interaction problem and an unsteady, inviscid wave rotor passage gradual opening problem. A model inlet port/passage charging problem is simulated and key features of the unsteady wave rotor flow field are identified. Lastly, the medium pressure inlet port and high pressure outlet port portion of the NASA Lewis Research Center experimental divider cycle is simulated and computed results are compared with experimental measurements. The model accurately predicts the wave timing within the rotor passages and the distribution of flow variables in the stationary inlet port region.

Welch, Gerard E.↗
NASA NTRS1997

Reduced Navier Stokes Relaxation Procedures for Internal Flows

In spite of significant advancement in the field of high speed computing, flow calculations involving complex geometries and/or flow behavior still require large amounts of CPU time and memory. In order to predict such flows without sacrificing grid convergence and accuracy, adaptive gridding techniques that provide optimal resolution are highly desirable. The present work combines multigrid techniques and domain decomposition concepts to provide local, solution adaptive, grid refinement. Several viscous compressible and incompressible, two and three-dimensional, flows with strong inviscid interaction and/or axial flow reversal, are considered with a segmented multigrid domain decomposition (SMGDD) procedure for which uniform meshes result in each domain. A pressure-based form of flux-vector splitting is applied to the Navier-Stokes equations, which are represented by an implicit lowest-order reduced Navier-Stokes (RNS) system and a purely diffusive, higher-order, deferred-corrector. A trapezoidal or box-like form of discretization insures that all mass conservation properties are satisfied at interfacial and outflow boundaries, even for this primitive-variable non-staggered grid computation. The SMGDD technique presented herein has previously been applied for incompressible two dimensional flows. The present work offers improvement in the gridding strategy, by allowing for disjoint subdomains that provide optimal resolution of disparate flow features. It also extends the SMGDD technique to three dimensional compressible flows. Laminar and turbulent flow in a backward facing step channel is considered; although the procedure is applicable to more severe geometries. The standard K-epsilon model is applied for turbulence closure. For Re greater than 400, differences between two-dimensional theory and experiment are resolved through a three dimensional simulation, which confirms the experimentally observed three dimensionality of the recirculation patterns on the upper and lower surfaces.

Rubin, Stanley G.↗
NASA NTRS1982

Flow aerodynamics modeling of an MHD swirl combustor - Calculations and experimental verification

This paper describes a computer code for calculating the flow dynamics of constant density flow in the second stage trumpet shaped nozzle section of a two stage MHD swirl combustor for application to a disk generator. The primitive pressure-velocity variable, finite difference computer code has been developed to allow the computation of inert nonreacting turbulent swirling flows in an axisymmetric MHD model swirl combustor. The method and program involve a staggered grid system for axial and radial velocities, and a line relaxation technique for efficient solution of the equations. Turbulence simulation is by way of a two-equation Kappa-epsilon model. The code produces as output the flowfield map of the nondimensional stream function, axial, and swirl velocity. Good agreement was obtained between the theoretical predictions and the qualitative experimental results. The best seed injector location for uniform seed distribution at combustor exit is with injector located centrally on the combustor axis at entrance to the second stage combustor.

Gupta, A. K.↗
NASA NTRS1989

Time-dependent viscous incompressible Navier-Stokes equations - The finite difference Galerkin formulation and streamfunction algorithms

Numerical techniques are developed to solve the Navier-Stokes equations for unsteady incompressible flow. The extension of the finite-difference Galerkin (FDG) method of Stephens et al. (1984) to the continuous-time case in two or three space dimensions is explained, and the numerical implementation of the method is discussed with particular attention to the staggered-MAC-grid primitive-variable discretization, the application of discrete mass balance to avoid problems inherent in FDG schemes, the direct interpretation of the FDG expansion variables as a discrete streamfunction, and a mass-balance approach to two-dimensional problems with throughflow or obstacles. Numerical results are presented graphically for the evolution of asymptotic steady flow in a driven cavity at Reynolds number 400, 1000, or 3200; good agreement with published experimental data is demonstrated, with accurate predictions of secondary-vortex formation from wall bubble recirculations at Reynolds number 1000.

Goodrich, John W.↗
NASA NTRS1981

Flow aerodynamics modeling of an MHD swirl combustor - Calculations and experimental verification

The paper describes a computer code for calculating the flow dynamics of a constant-density flow in the second-stage trumpet shaped nozzle section of a two-stage MHD swirl combustor for application to a disk generator. The primitive pressure-velocity variable, finite-difference computer code has been developed for the computation of inert nonreacting turbulent swirling flows in an axisymmetric MHD model swirl combustor. The method and program involve a staggered grid system for axial and radial velocities, and a line relaxation technique for the efficient solution of the equations. The code produces as output the flow field map of the nondimensional stream function, axial and swirl velocity. It was found that the best location for seed injection to obtain a uniform distribution at the combustor exit is in the central location for seed injected at the entrance to the second stage combustor.

Gupta, A. K.↗