Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Moment 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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Comparative Study of Advanced Turbulence Models for Turbomachinery

A computational study has been undertaken to study the performance of advanced phenomenological turbulence models coded in a modular form to describe incompressible turbulent flow behavior in two dimensional/axisymmetric and three dimensional complex geometry. The models include a variety of two equation models (single and multi-scale k-epsilon models with different near wall treatments) and second moment algebraic and full Reynolds stress closure models. These models were systematically assessed to evaluate their performance in complex flows with rotation, curvature and separation. The models are coded as self contained modules that can be interfaced with a number of flow solvers. These modules are stand alone satellite programs that come with their own formulation, finite-volume discretization scheme, solver and boundary condition implementation. They will take as input (from any generic Navier-Stokes solver) the velocity field, grid (structured H-type grid) and computational domain specification (boundary conditions), and will deliver, depending on the model used, turbulent viscosity, or the components of the Reynolds stress tensor. There are separate 2D/axisymmetric and/or 3D decks for each module considered. The modules are tested using Rocketdyn's proprietary code REACT. The code utilizes an efficient solution procedure to solve Navier-Stokes equations in a non-orthogonal body-fitted coordinate system. The differential equations are discretized over a finite-volume grid using a non-staggered variable arrangement and an efficient solution procedure based on the SIMPLE algorithm for the velocity-pressure coupling is used. The modules developed have been interfaced and tested using finite-volume, pressure-correction CFD solvers which are widely used in the CFD community. Other solvers can also be used to test these modules since they are independently structured with their own discretization scheme and solver methodology. Many of these modules have been independently tested by Professor C.P. Chen and his group at the University of Alabama at Huntsville (UAH) by interfacing them with own flow solver (MAST).

Hadid, Ali H.↗
NASA NTRS1999

Some Continuum Aspects of Data Assimilation

A long-sought goal in data assimilation is to build more fully the time dimension into the data assimilation process. Truly four-dimensional data assimilation requires evolving second-moment information, namely the estimation error covariance, along with the state estimate itself. The continuum evolution equations for the estimation error covariance constitute a system of partial differential equations in six space dimensions, forced by a model error covariance and by the observations themselves. The high dimensionality of this system poses a difficult computational problem. Numerous methods have been proposed for approximate, discrete solution of this system of equations. Four-dimensional variational schemes solve these equations indirectly, while schemes based on Kalman filtering solve the equations more directly, usually on a small subspace of the full six-dimensional space. In both cases, most proposed solution methods are not derived from the continuum covariance evolution equations themselves. Instead, most methods simply inherit the discretization used for evolving the state estimate, either in the form of a discrete tangent linear model or, in the case of ensemble schemes, in the form of a fully nonlinear discrete model. In this lecture we show that solutions of the continuum covariance evolution equations possess simple properties that are not readily satisfied upon discretization of the equations. we give numerical examples illustrating that "inherited" discretizations sometimes fail to produce meaningfully accurate solutions. Finally, we suggest alternative discretization methods that may yield more faithful approximate solutions of the continuum problem.

Cohn, Stephen E.↗
NASA NTRS2001

Computation of Flow Over a Drag Prediction Workshop Wing/Body Transport Configuration Using CFL3D

A Drag Prediction Workshop was held in conjunction with the 19th AIAA Applied Aerodynamics Conference in June 2001. The purpose of the workshop was to assess the prediction of drag by computational methods for a wing/body configuration (DLR-F4) representative of subsonic transport aircraft. This report details computed results submitted to this workshop using the Reynolds-averaged Navier-Stokes code CFL3D. Two supplied grids were used: a point-matched 1-to-1 multi-block grid, and an overset multi-block grid. The 1-to-1 grid, generally of much poorer quality and with less streamwise resolution than the overset grid, is found to be too coarse to adequately resolve the surface pressures. However, the global forces and moments are nonetheless similar to those computed using the overset grid. The effect of three different turbulence models is assessed using the 1-to-1 grid. Surface pressures are very similar overall, and the drag variation due to turbulence model is 18 drag counts. Most of this drag variation is in the friction component, and is attributed in part to insufficient grid resolution of the 1-to-1 grid. The misnomer of 'fully turbulent' computations is discussed; comparisons are made using different transition locations and their effects on the global forces and moments are quantified. Finally, the effect of two different versions of a widely used one-equation turbulence model is explored.

Rumsey, Christopher L.↗
NASA NTRS2002

Elliptic Relaxation of a Tensor Representation for the Redistribution Terms in a Reynolds Stress Turbulence Model

A formulation to include the effects of wall proximity in a second-moment closure model that utilizes a tensor representation for the redistribution terms in the Reynolds stress equations is presented. The wall-proximity effects are modeled through an elliptic relaxation process of the tensor expansion coefficients that properly accounts for both correlation length and time scales as the wall is approached. Direct numerical simulation data and Reynolds stress solutions using a full differential approach are compared for the case of fully developed channel flow.

Carlson, J. R.↗
NASA NTRS2016

Development of Constraint Force Equation Methodology for Application to Multi-Body Dynamics Including Launch Vehicle Stage Seperation

The objective of this report is to develop and implement a physics based method for analysis and simulation of multi-body dynamics including launch vehicle stage separation. The constraint force equation (CFE) methodology discussed in this report provides such a framework for modeling constraint forces and moments acting at joints when the vehicles are still connected. Several stand-alone test cases involving various types of joints were developed to validate the CFE methodology. The results were compared with ADAMS(Registered Trademark) and Autolev, two different industry standard benchmark codes for multi-body dynamic analysis and simulations. However, these two codes are not designed for aerospace flight trajectory simulations. After this validation exercise, the CFE algorithm was implemented in Program to Optimize Simulated Trajectories II (POST2) to provide a capability to simulate end-to-end trajectories of launch vehicles including stage separation. The POST2/CFE methodology was applied to the STS-1 Space Shuttle solid rocket booster (SRB) separation and Hyper-X Research Vehicle (HXRV) separation from the Pegasus booster as a further test and validation for its application to launch vehicle stage separation problems. Finally, to demonstrate end-to-end simulation capability, POST2/CFE was applied to the ascent, orbit insertion, and booster return of a reusable two-stage-to-orbit (TSTO) vehicle concept. With these validation exercises, POST2/CFE software can be used for performing conceptual level end-to-end simulations, including launch vehicle stage separation, for problems similar to those discussed in this report.

Pamadi, Bandu N.↗
NASA NTRS1992

Turbulence modelling in CFD: Present status, future prospects

Information is given in viewgraph form for turbulence modeling in computational fluid dynamics (CFD). The Eddy Viscosity Models (EVM), Algebraic Second Moment Closures (ASM), and Differential Second-Moment Closures (DSM) are considered. It is concluded that EVM's, ASM's, and DSM's will remain in use, though with a steady decline in importance of EVM's and ASM's in favor of DSM's. Improved versions of low-Re two-equation EVM's should lead to more reliable predictions of separated flows than are achievable at present. Further refinement of sub-models in second moment closures can be expected throughout this decade. There will be increasing attention given to interfacing SMC with higher order approaches such as LES, and an increased use of two-time-scale schemes providing distinct time scales for large and fairly small eddies.

Launder, Brian E.↗
DOE OSTI2025

Impact of Magnetic-field-driven Anisotropies on the Equation of State Probed in Neutron Star Mergers

Binary neutron star mergers can produce extreme magnetic fields, some of which can lead to strong magnetar-like remnants. While strong magnetic fields have been shown to affect the dynamics of outflows and angular momentum transport in the remnant, they can also crucially alter the properties of nuclear matter probed in the merger. In this work, we provide a first assessment of the latter, determining the strength of the pressure anisotropy caused by Landau-level quantization and the anomalous magnetic moment. To this end, we perform the first numerical relativity simulation with a magnetic polarization tensor and a magnetic-field-dependent equation of state using a new algorithm we present here, which also incorporates a mean-field dynamo model to control the magnetic field strength present in the merger remnant. Our results show that—in the most optimistic case—corrections to the anisotropy can be in excess of 10% and are potentially largest in the outer layers of the remnant. This work paves the way for a systematic investigation of these effects.

General relativity↗
NASA NTRS1991

Navier-Stokes predictions of pitch damping for axisymmetric shell using steady coning motion

Previous theoretical investigations have proposed that the side force and moment acting on a body of revolution in steady coning motion could be related to the pitch-damping force and moment. In the current research effort, this approach is applied to produce predictions of the pitch damping for axisymmetric shell. The flow fields about these projectiles undergoing steady coning motion are successfully computed using a parabolized Navier-Stokes computational approach which makes use of a rotating coordinate frame. The governing equations are modified to include the centrifugal and Coriolis force terms due to the rotating coordinate frame. From the computed flow field, the side moments due to coning motion, spinning motion, and combined spinning and coning motion are used to determine the pitch-damping coefficients. Computations are performed for two generic shell configurations, a secant-ogive-cylinder and a secant-ogive-cylinder-boattail.

Weinacht, Paul↗
NASA NTRS2001

Pluto and Triton: Interactions Between Volatiles and Dynamics

Volatiles moving across the surfaces of Pluto and Triton can give rise to interesting dynamical consequences. Conversely, measurement of dynamical states can help constrain the movement of volatiles and interior structure of both bodies. Polar wander may theoretically occur on both Triton and Pluto. Triton's obliquity is low, so that the equatorial regions receive more insolation than the poles. Hence there is a tendency for nitrogen ice to sublime at the equator and condense at the poles, creating polar caps. If the nitrogen supply is large enough, then these caps could move in approximately 10(exp 5) years the global equivalent of 200 m of ice to the poles. At this point the equatorial moment of inertia becomes larger than the moment of inertia measured about the rotation axis, so that Triton overbalances and becomes dynamically unstable. The satellite then undergoes polar wander, restoring stability when the new equator contains the excess matter. Hence the pole may be continually wandering. Neptune raises a permanent tidal bulge on Triton, so that the satellite's surface is elongated like a football, with the long axis pointing at Neptune. This is expected to be the axis about which the pole wanders. Volatile migration would resurface the satellite to some depth and wandering would disturb leading side/trailing side crater statistics. Additional information is contained in the original extended abstract.

Rubincam, D. P.↗
NASA NTRS2003

Quadrature Moments Method for the Simulation of Turbulent Reactive Flows

A sub-filter model for reactive flows, namely the DQMOM model, was formulated for Large Eddy Simulation (LES) using the filtered mass density function. Transport equations required to determine the location and size of the delta-peaks were then formulated for a 2-peak decomposition of the FDF. The DQMOM scheme was implemented in an existing structured-grid LES solver. Simulations of scalar shear layer using an experimental configuration showed that the first and second moments of both reactive and inert scalars are in good agreement with a conventional Lagrangian scheme that evolves the same FDF. Comparisons with LES simulations performed using laminar chemistry assumption for the reactive scalar show that the new method provides vast improvements at minimal computational cost. Currently, the DQMOM model is being implemented for use with the progress variable/mixture fraction model of Pierce. Comparisons with experimental results and LES simulations using a single-environment for the progress-variable are planned. Future studies will aim at understanding the effect of increase in environments on predictions.

Raman, Venkatramanan↗
NASA NTRS1976

Computation and sensitivity considerations of the Ricatti equation in spacecraft CMG steering laws

The steady-state matrix Ricatti equation is studied with respect to sensitivity and computational burden for a linear regulator controller in spacecraft. The class of systems studied use control moment gyros (CMGs) as the torque driver, assuming that the CMGs are free of gimbal stops. Different solution techniques are evaluated relative to computation time and accuracy to determine which solution method is most acceptable with respect to the prescribed requirements. It is shown that it is not feasible to eliminate off-diagonal Ricatti equation terms in the relevent updata equation due to their high sensitivity, that Potter's algebraic technique requires the least computation time of the methods studied, and that the simple Euler method yields the best results when accuracy and computation time are equally weighted. If an optimal linear regulator formulation is employed, a full Ricatti equation solution of n(n + 1)2 elements of the symmetric positive definite gain matrix which is the solution to the matrix Ricatti equation must be computed at each update cycle.

Colburn, B. K.↗
NASA NTRS1994

Neutron star moments of inertia

An approximation for the moment of inertia of a neutron star in terms of only its mass and radius is presented, and insight into it is obtained by examining the behavior of the relativistic structural equations. The approximation is accurate to approximately 10% for a variety of nuclear equations of state, for all except very low mass stars. It is combined with information about the neutron-star crust to obtain a simple expression (again in terms only of mass and radius) for the fractional moment of inertia of the crust.

Ravenhall, D. G.↗
DOE OSTI2020

Generalized Kompaneets formalism for inelastic neutrino-nucleon scattering in supernova simulations

Based on the Kompaneets approximation, we develop a robust methodology to calculate spectral redistribution via inelastic neutrino-nucleon scattering in the context of core-collapse supernova simulations. The resulting equations conserve lepton number to machine precision and scale linearly, not quadratically, with number of energy groups. The formalism also provides an elegant means to derive the rate of energy transfer to matter which, as it must, automatically goes to zero when the neutrino radiation field is in thermal equilibrium. Furthermore, we derive the next higher order in $ϵ/mc^2$ correction to the neutrino Kompaneets equation. Unlike other Kompaneets schema, ours also generalizes to the case of anisotropic angular distributions, while retaining the conservative form that is a hallmark of the classical Kompaneets equation. In this work, our formalism enables immediate incorporation into supernova codes that follow the spectral angular moments of the neutrino radiation fields.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1981

Scaling from Jupiter to pulsar and the acceleration of cosmic ray particles by pulsars

The total energy output per second from the 73 known pulsars is estimated. By following the principle of similitude (Rayleigh, 1915), an expression is derived for the energy generated by a rotating magnetized celestial body, and the dipole moment of the body is related linearly with its angular momentum. The scaling constant of the two equations combined is determined, and then checked by comparing the deduced masses of the 73 pulsars with the theoretical limits for neutron stars. By using the scaling constant and deduced pulsar masses, the rates of acceleration of high-energy particles by pulsars, thence the energy output from pulsars, can be estimated. Results show that the Crab pulsar outshines all the others combined.

Fan, Y. C.↗
NASA NTRS1983

Hypersonic flows over biconics using a variable-effective-gamma, Parabolized-Navier-Stokes code

Hypersonic flows over straight and bent biconics are calculated for a range of freestream conditions in which the gas behind the shock is treated as either perfect or real. The Parabolized-Navier-Stokes (PNS) equations form the basis of the approximation scheme. Good comparisons with experimental data for pressures, forces and moments, heat transfer, and oil-flow patterns serve to validate the perfect-gas version of the code. Circumferential velocity vector plots further aid in the interpretation of leeside oil-flow patterns. A variable-effective-gamma (VEG) option is implemented for the real-gas calculations. Gamma, now defined as the ratio of enthalpy to internal energy, is determined from a locally valid linear relation in enthalpy and pressure at every mesh point which in turn is calculated from a benchmark equilibrium code. The VEG option is easily incorporated into a host PNS code because it uses the underlying perfect-gas structure of the code. Comparisons to experimental data for heat transfer and shock shape in high enthalpy air have been obtained using the VEG option.

Gnoffo, P. A.↗
NASA NTRS1985

Flight-Stability Extremes Calculated From the Uncertainty Ellipse

Tangents and normals to ellipse derived for situations with two correlated aerodynamic force or moment variables having normally (Gaussian) distributed uncertainty. Although transformations and vector calculations tedious when done by hand, equations in closed form and give rapid results when performed by computers or programable desk calculators. Probability extremes of ratios and magnitude obtained directly, without repeated table lookup or other complication.

Dominik, D. F.↗
NASA NTRS1988

Scaling laws for homogeneous turbulent shear flows in a rotating frame

The scaling properties of plane homogeneous turbulent shear flows in a rotating frame are examined mathematically by a direct analysis of the Navier-Stokes equations. It is proved that two such shear flows are dynamically similar if and only if their initial dimensionless energy spectrum E star (k star, 0), initial dimensionless shear rate SK sub 0/epsilon sub 0, initial Reynolds number K squared sub 0/nu epsilon sub 0, and the ration of the rotation rate to the shear rate omega/S are identical. Consequently, if universal equilibrium states exist, at high Reynolds numbers, they will only depend on the single parameter omega/S. The commonly assumed dependence of such equilibrium states on omega/S through the Richardson number Ri=-2(omega/S)(1-2 omega/S) is proven to be inconsistent with the full Navier-Stokes equations and to constitute no more than a weak approximation. To be more specific, Richardson number similarity is shown to only rigorously apply to certain low-order truncations of the Navier-Stokes equations (i.e., to certain second-order closure models) wherein closure is achieved at the second-moment level by assuming that the higher-order moments are a small perturbation of their isotropic states. The physical dependence of rotating turbulent shear flows on omega/S is discussed in detail along with the implications for turbulence modeling.

Speziale, Charles G.↗
NASA NTRS1989

Scaling laws for homogeneous turbulent shear flows in a rotating frame

The scaling properties of plane homogeneous turbulent shear flows in a rotating frame are examined mathematically by a direct analysis of the Navier-Stokes equations. It is proved that two such shear flows are dynamically similar if and only if their initial dimensionless energy spectrum E star (k star, 0), initial dimensionless shear rate SK sub 0/epsilon sub 0, initial Reynolds number K squared sub 0/nu epsilon sub 0, and the ration of the rotation rate to the shear rate omega/S are identical. Consequently, if universal equilibrium states exist, at high Reynolds numbers, they will only depend on the single parameter omega/S. The commonly assumed dependence of such equilibrium states on omega/S through the Richardson number Ri = -2(omega/S)(1-2 omega/S) is proven to be inconsistent with the full Navier-Stokes equations and to constitute no more than a weak approximation. To be more specific, Richardson number similarity is shown to only rigorously apply to certain low-order truncations of the Navier-Stokes equations (i.e., to certain second-order closure models) wherein closure is achieved at the second-moment level by assuming that the higher-order moments are a small perturbation of their isotropic states. The physical dependence of rotating turbulent shear flows on omega/S is discussed in detail along with the implications for turbulence modeling.

Speziale, Charles G.↗