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 307 records · Page 17

Finite element analysis of inviscid subsonic boattail flow

A finite element code for analysis of inviscid subsonic flows over arbitrary nonlifting planar or axisymmetric bodies is described. The code solves a novel primitive variable formulation of the coupled irrotationality and compressible continuity equations. Results for flow over a cylinder, a sphere, and a NACA 0012 airfoil verify the code. Computed subcritical flows over an axisymmetric boattailed afterbody compare well with finite difference results and experimental data. Interative coupling with an integral turbulent boundary layer code shows strong viscous effects on the inviscid flow. Improvements in code efficiency and extensions to transonic flows are discussed.

Chima, R. V.↗

Finite element analysis of inviscid subsonic boattail flow

A finite element code for analysis of inviscid subsonic flows over arbitrary nonlifting planar or axisymmetric bodies is described. The code solves a novel primitive variable formulation of the coupled irrotationality and compressible continuity equations. Results for flows over a cylinder, a sphere, and a NACA 0012 airfoil verify the code. Computed subcritical flows over an axisymmetric boattailed afterbody compare well with finite difference results and experimental data. Iterative coupling with an integral turbulent boundary layer code shows strong viscous effects on the inviscid flow. Improvements in code efficiency and extensions to transonic flows are discussed.

Chima, R. V.↗

A study of non-unique solutions of the two-dimensional boundary layer equations at laminar separation and reattachment points

Nonunique laminar boundary layer equation solutions in direct problem calculations are identified for the case of accelerating flow. As a separation or reattachment point is approached, the multiple solutions approach each other and become identical. The computer code used to generate these results was developed for the solution of compressible, laminar or turbulent boundary layer, and free wake problems, in either direct or inverse mode. Similarity solutions in either a primitive variable or a stream function form are possible, and the resulting equations are solved by means of a modified Keller's Box scheme in which the energy equation and turbulence modeling equations are solved simultaneously with the continuity and momentum equations. Examples illustrating the nature of the solutions at the separation and reattachment points are presented.

Drela, M.↗

Probabilistic lifetime strength of aerospace materials via computational simulation

The results of a second year effort of a research program are presented. The research included development of methodology that provides probabilistic lifetime strength of aerospace materials via computational simulation. A probabilistic phenomenological constitutive relationship, in the form of a randomized multifactor interaction equation, is postulated for strength degradation of structural components of aerospace propulsion systems subjected to a number of effects of primitive variables. These primitive variables often originate in the environment and may include stress from loading, temperature, chemical, or radiation attack. This multifactor interaction constitutive equation is included in the computer program, PROMISS. Also included in the research is the development of methodology to calibrate the constitutive equation using actual experimental materials data together with the multiple linear regression of that data.

Boyce, Lola↗

Computational simulation of probabilistic lifetime strength for aerospace materials subjected to high temperature, mechanical fatigue, creep, and thermal fatigue

The results of a fourth year effort of a research program conducted for NASA-LeRC by The University of Texas at San Antonio (UTSA) are presented. The research included on-going development of methodology that provides probabilistic lifetime strength of aerospace materials via computational simulation. A probabilistic material strength degradation model, in the form of a randomized multifactor interaction equation, is postulated for strength degradation of structural components of aerospace propulsion systems subjected to a number of effects or primitive variables. These primitive variables may include high temperature, fatigue, or creep. In most cases, strength is reduced as a result of the action of a variable. This multifactor interaction strength degradation equation was randomized and is included in the computer program, PROMISC. Also included in the research is the development of methodology to calibrate the above-described constitutive equation using actual experimental materials data together with regression analysis of that data, thereby predicting values for the empirical material constants for each effect or primitive variable. This regression methodology is included in the computer program, PROMISC. Actual experimental materials data were obtained from industry and the open literature for materials typically for applications in aerospace propulsion system components. Material data for Inconel 718 was analyzed using the developed methodology.

Boyce, Lola↗

Computational simulation of probabilistic lifetime strength for aerospace materials subjected to high temperature, mechanical fatigue, creep and thermal fatigue

This report presents the results of a fourth year effort of a research program, conducted for NASA-LeRC by the University of Texas at San Antonio (UTSA). The research included on-going development of methodology that provides probabilistic lifetime strength of aerospace materials via computational simulation. A probabilistic material strength degradation model, in the form of a randomized multifactor interaction equation, is postulated for strength degradation of structural components of aerospace propulsion systems subject to a number of effects or primitive variables. These primitive variables may include high temperature, fatigue or creep. In most cases, strength is reduced as a result of the action of a variable. This multifactor interaction strength degradation equation has been randomized and is included in the computer program, PROMISS. Also included in the research is the development of methodology to calibrate the above-described constitutive equation using actual experimental materials data together with regression analysis of that data, thereby predicting values for the empirical material constants for each effect or primitive variable. This regression methodology is included in the computer program, PROMISC. Actual experimental materials data were obtained from industry and the open literature for materials typically for applications in aerospace propulsion system components. Material data for Inconel 718 has been analyzed using the developed methodology.

Boyce, Lola↗

Application of an Affine Nonlinear Galerkin Reduced-order Model to Compressible Fluid Flows

Galerkin reduced-order models (ROMs) often struggle to accurately capture multiscale fluid physics in challenging flow regimes such as flows experiencing compressibility effects. This in part stems from the global nature of both the basis construction problem and the spectral formulation itself. In this work, a multi-basis ROM is developed in an affine space based on proper orthogonal decomposition (POD) by projecting the full Navier-Stokes equations expressed in terms of the specific volume, velocity, and pressure primitive variables. The model is applied to high-fidelity numerical simulation datasets obtained for a canonical compressible flow configuration: the flow over a backward facing step at different subsonic Mach numbers. It is observed that application of an eigenvalue reassignment (ER) stabilization method is required to avoid early divergence of the ROM predictions for this configuration in the three Mach numbers tested. The sensitivity of the POD-ROM results to the choice of parameters in the stabilization algorithm is discussed.

reduced-order model↗

Application of an Affine Nonlinear Galerkin Reduced-order Model to Compressible Fluid Flows

Galerkin reduced-order models (ROMs) often struggle to accurately capture multiscale fluid physics in challenging flow regimes such as flows experiencing compressibility effects. This in part stems from the global nature of both the basis construction problem and the spectral formulation itself. In this work, a multi-basis ROM is developed in an affine space based on proper orthogonal decomposition (POD) by projecting the full Navier-Stokes equations expressed in terms of the specific volume, velocity, and pressure primitive variables. The model is applied to high-fidelity numerical simulation datasets obtained for a canonical compressible flow configuration: the flow over a backward facing step at different subsonic Mach numbers. It is observed that application of an eigenvalue reassignment (ER) stabilization method is required to avoid early divergence of the ROM predictions for this configuration in the three Mach numbers tested. The sensitivity of the POD-ROM results to the choice of parameters in the stabilization algorithm is discussed.

reduced-order model↗

A fractional step solution method for the unsteady incompressible Navier-Stokes equations in generalized coordinate systems

The time-dependent, three-dimensional incompressible Navier-Stokes equations are presently solved in generalized coordinate systems by means of a fractional-step method whose primitive variable formulation uses as dependent variables, in place of the Cartesian components of the velocity: (1) pressure (defined at the center of the computational cell), and (2) volume fluxes across the faces of the cells. The momentum equations are solved by means of an approximate factorization method. A novel 'ZEBRA' scheme incorporating four-color ordering efficiently solves the Poisson equation. Illustrative two- and three-dimensional laminar flow test cases are computed and evaluated relative to extant numerical and experimental results, and good agreement is obtained.

Rosenfeld, Moshe↗

Implicit finite-difference procedures for the computation of vortex wakes

Implicit finite-difference procedures for the primitive form of the incompressible Navier-Stokes and the compressible Euler equations are used to compute vortex wake flows. The partial differential equations in strong conservation-law form are transformed to cluster grid points in regions with large changes in vorticity. In addition to clustering, fourth-order accurate, spatial difference operators are used to help resolve the flow-field gradients. The use of implicit time-differencing permits large time steps to be taken since temporal variations are typically small. Computational efficiency is achieved by approximate factorization. Both two-dimensional and preliminary three-dimensional calculations are described and qualitatively compared with existing experimental data.

Steger, J. L.↗

Toward higher order tests of the gravitational interaction

Analyses and interpretations of experiments which test post-Newtonian gravity are usually done under the assumption that gravity is a metric field phenomenon - a manifestation of space-time geometry. This, however, is unnecessary and one can start at a more primitive level - that there simply exists a phenomenological, gravitational, many-body equation of motion which must be determined by a package of observations. In fact, over the last couple decades, a diverse collection of solar system interbody tracking observations, supplemented by data from the binary pulsar system PSR 1913 + 16, has completely mapped out the first post-Newtonian order. After the fact, using empirically determined equations of motion, along with some observed properties of nongravitational clocks and rulers and conservation laws for energy, momentum and angular momentum, a post-Newtonian Lagrangian can be constructed, a geometrical space-time metric field conceptual interpretation can be developed, Lorentz invariance of the equations of motion can be shown, and the equations of motion are found to agree with the predictions of Einstein's gravitational theory, General Relativity, within experimental accuracy.

Nordtvedt, Ken↗

Numerical solution of the Navier-Stokes equations for high Reynolds number incompressible turbulent flow

The full Navier-Stokes equations for incompressible turbulent flow must be solved to accurately represent all flow phenomena which occur in a high Reynolds number incompressible flow. A two layer algebraic eddy viscosity turbulence model is used to represent the Reynolds stress in the primitive variable formulation. The development of the boundary-fitted coordinate systems makes the numerical solution of these equations feasible for arbitrarily shaped bodies. The nondimensional time averaged Navier-Stokes equations, including the turbulence mode, are represented by finite difference approximations in the transformed plane. The resulting coupled system of nonlinear algebraic equations is solved using a point successive over relaxation iteration. The test case considered was a NACA 64A010 airfoil section at an angle of attack of two degrees and a Reynolds number of 2,000,000.

Thompson, D. S.↗

Random field estimation approach to robot dynamics

The difference equations of Kalman filtering and smoothing recursively factor and invert the covariance of the output of a linear state-space system driven by a white-noise process. Here it is shown that similar recursive techniques factor and invert the inertia matrix of a multibody robot system. The random field models are based on the assumption that all of the inertial (D'Alembert) forces in the system are represented by a spatially distributed white-noise model. They are easier to describe than the models based on classical mechanics, which typically require extensive derivation and manipulation of equations of motion for complex mechanical systems. With the spatially random models, more primitive locally specified computations result in a global collective system behavior equivalent to that obtained with deterministic models. The primary goal of applying random field estimation is to provide a concise analytical foundation for solving robot control and motion planning problems.

Rodriguez, Guillermo↗

The MOOSE fluid properties module

The fluid properties module in MOOSE~\cite{lindsay2022moose} serves a variety of fluid simulation applications based on MOOSE, including the MOOSE Navier Stokes module~\cite{moose_ns}, Pronghorn~\cite{pgh}, the MOOSE thermal hydraulics module, SAM~\cite{sam}, RELAP-7~\cite{relap7}, Sockeye~\cite{sockeye}, Pronghorn-subchannel~\cite{subchannel} and the MOOSE porous flow module~\cite{porous}. These applications are used to solve coarse mesh multi-dimensional thermal-hydraulics~\cite{pgh}, 1D systems analysis~\cite{sam,relap7} in nuclear reactor analysis, heat pipe modeling~\cite{sockeye} and porous flow simulations~\cite{porous} for underground gas storage and water seepage. The use of consistent fluid properties across fluid flow applications facilitates coupled flow simulations~\cite{anl_sam_pgh,osti_1889653}. Each application has historically driven the implementation of several fluid properties, which were later extended to be compatible with other applications. The unique diversity of applications of the module, due to its presence in MOOSE, has driven its expansion to new fluids, such as advanced nuclear reactor coolants and, more recently, arbitrary functions or tables-based property definitions, as detailed in section~\ref{content}, as well as numerous thermophysical properties and variable sets, as detailed in subsection~\ref{sec:prop}. The need for different discretizations of flow equations based on the compressibility of the fluid has motivated support for both primitive (pressure- and temperature-based) and conservative (internal energy- and specific volume-based) flow variables in the module. Thermodynamic relations are used to convert between these two formulations, as needed. The module serves a dual purpose of providing fluid properties for direct use in flow simulations and facilitating the implementation of user-specific fluid properties. Contributions of new properties for existing fluids or new fluids are strongly encouraged.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

An incompressible Navier-Stokes flow solver in three-dimensional curvilinear coordinate systems using primitive variables

An implicit, finite-difference computer code has been developed to solve the incompressible Navier-Stokes equations in a three-dimensional, curvilinear coordinate system. The pressure-field solution is based on the pseudo compressibility approach in which the time derivative pressure term is introduced into the mass conservation equation to form a set of hyperbolic equations. The solution procedure employs an implicit, approximate factorization scheme. The Reynolds stresses, that are uncoupled from the implicit scheme, are lagged by one time-step to facilitate implementing various levels of the turbulence model. Test problems for external and internal flows are computed, and the results are compared with existing experimental data. The application of this technique for general three-dimensional problems is then demonstrated.

Kwak, D.↗

The rotation of the Uranian system

The rotation of Uranus is examined for clues as to the origin of the Solar System. Both theories based on the formation of planets through the accretion of small planetesimals, and theories based on the formation of giant gaseous protoplanets through a gravitational instability in the primitive solar nebula allow for qualitative explanations of the large tilt of Uranus's equator to the orbital plane, and the fact that its satellites lie in the equatorial plane. Models of the planetary interior show that the mass ratio of ice-forming materials to rock in Uranus's interior must be more than about three if the rotation period is about 16 h. Such a large ratio seems to exclude those accretional theories that require most of the nebular gas to be heated to relatively high temperatures before being accreted into the planet.

Podolak, M.↗

Mathematical aspects of finite element methods for incompressible viscous flows

Mathematical aspects of finite element methods are surveyed for incompressible viscous flows, concentrating on the steady primitive variable formulation. The discretization of a weak formulation of the Navier-Stokes equations are addressed, then the stability condition is considered, the satisfaction of which insures the stability of the approximation. Specific choices of finite element spaces for the velocity and pressure are then discussed. Finally, the connection between different weak formulations and a variety of boundary conditions is explored.

Gunzburger, M. D.↗