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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

A BBGKY framework for fluid turbulence

A framework is presented for a systematic kinetic theory of turbulence originating from the Liouville equation for the Fourier coefficients of fluid variables. The real and imaginary parts of these Fourier coefficients play the role that particle coordinates (positions and momenta) play in the BBGKY theory. The basic relations of the problem are the incompressible Navier-Stokes equations in two dimensions with zero viscosity, with the probability distributions of Fourier coefficients rather than moments being the basic variables of the theory. A kinetic equation is derived and shown to possess a number of requirements that any reasonable kinetic equation must have: conservation laws, positive-definite spectral densities, and an H-theorem. The major lack in the theory is any reliable information on the relaxation predicted by the complicated linear operator H. Closure of the hierarchy is achieved by the hypothesis that the five-coefficient correlation function is negligible. Problems associated with inclusion of viscosity and external driving forces are discussed.

Montgomery, D.↗
NASA NTRS2012

A New and General Formulation of the Parametric HFGMC Micromechanical Method for Three-Dimensional Multi-Phase Composites

The recent two-dimensional (2-D) parametric formulation of the high fidelity generalized method of cells (HFGMC) reported by the authors is generalized for the micromechanical analysis of three-dimensional (3-D) multiphase composites with periodic microstructure. Arbitrary hexahedral subcell geometry is developed to discretize a triply periodic repeating unit-cell (RUC). Linear parametric-geometric mapping is employed to transform the arbitrary hexahedral subcell shapes from the physical space to an auxiliary orthogonal shape, where a complete quadratic displacement expansion is performed. Previously in the 2-D case, additional three equations are needed in the form of average moments of equilibrium as a result of the inclusion of the bilinear terms. However, the present 3-D parametric HFGMC formulation eliminates the need for such additional equations. This is achieved by expressing the coefficients of the full quadratic polynomial expansion of the subcell in terms of the side or face average-displacement vectors. The 2-D parametric and orthogonal HFGMC are special cases of the present 3-D formulation. The continuity of displacements and tractions, as well as the equilibrium equations, are imposed in the average (integral) sense as in the original HFGMC formulation. Each of the six sides (faces) of a subcell has an independent average displacement micro-variable vector which forms an energy-conjugate pair with the transformed average-traction vector. This allows generating symmetric stiffness matrices along with internal resisting vectors for the subcells which enhances the computational efficiency. The established new parametric 3-D HFGMC equations are formulated and solution implementations are addressed. Several applications for triply periodic 3-D composites are presented to demonstrate the general capability and varsity of the present parametric HFGMC method for refined micromechanical analysis by generating the spatial distributions of local stress fields. These applications include triply periodic composites with inclusions in the form of a cavity, spherical inclusion, ellipsoidal inclusion, discontinuous aligned short fiber. A 3-D repeating unit-cell for foam material composite is simulated.

Haj-Ali, Rami↗
NASA NTRS1998

Application of Model Based Parameter Estimation for RCS Frequency Response Calculations Using Method of Moments

An implementation of the Model Based Parameter Estimation (MBPE) technique is presented for obtaining the frequency response of the Radar Cross Section (RCS) of arbitrarily shaped, three-dimensional perfect electric conductor (PEC) bodies. An Electric Field Integral Equation (EFTE) is solved using the Method of Moments (MoM) to compute the RCS. The electric current is expanded in a rational function and the coefficients of the rational function are obtained using the frequency derivatives of the EFIE. Using the rational function, the electric current on the PEC body is obtained over a frequency band. Using the electric current at different frequencies, RCS of the PEC body is obtained over a wide frequency band. Numerical results for a square plate, a cube, and a sphere are presented over a bandwidth. Good agreement between MBPE and the exact solution over the bandwidth is observed.

Reddy, C. J.↗
NASA NTRS1998

The Real-Time Wall Interference Correction System of the NASA Ames 12-Foot Pressure Wind Tunnel

An improved version of the Wall Signature Method was developed to compute wall interference effects in three-dimensional subsonic wind tunnel testing of aircraft models in real-time. The method may be applied to a full-span or a semispan model. A simplified singularity representation of the aircraft model is used. Fuselage, support system, propulsion simulator, and separation wake volume blockage effects are represented by point sources and sinks. Lifting effects are represented by semi-infinite line doublets. The singularity representation of the test article is combined with the measurement of wind tunnel test reference conditions, wall pressure, lift force, thrust force, pitching moment, rolling moment, and pre-computed solutions of the subsonic potential equation to determine first order wall interference corrections. Second order wall interference corrections for pitching and rolling moment coefficient are also determined. A new procedure is presented that estimates a rolling moment coefficient correction for wings with non-symmetric lift distribution. Experimental data obtained during the calibration of the Ames Bipod model support system and during tests of two semispan models mounted on an image plane in the NASA Ames 12 ft. Pressure Wind Tunnel are used to demonstrate the application of the wall interference correction method.

Ulbrich, Norbert↗
DOE OSTI2020

Physical, numerical, and computational challenges of modeling neutrino transport in core-collapse supernovae

The proposal that core collapse supernovae are neutrino driven is still the subject of active investigation more than 50 years after the seminal paper by Colgate and White. The modern version of this paradigm, which we owe to Wilson, proposes that the supernova shock wave is powered by neutrino heating, mediated by the absorption of electron-flavor neutrinos and antineutrinos emanating from the proto-neutron star surface, or neutrinosphere. Neutrino weak interactions with the stellar core fluid, the theory of which is still evolving, are flavor and energy dependent. The associated neutrino mean free paths extend over many orders of magnitude and are never always small relative to the stellar core radius. Thus, neutrinos are never always fluid like. Instead, a kinetic description of them in terms of distribution functions that determine the number density of neutrinos in the six-dimensional phase space of position, direction, and energy, for both neutrinos and antineutrinos of each flavor, or in terms of angular moments of these neutrino distributions that instead provide neutrino number densities in the four-dimensional phase-space subspace of position and energy, is needed. In turn, the computational challenge is twofold: (i) to map the kinetic equations governing the evolution of these distributions or moments onto discrete representations that are stable, accurate, and, perhaps most important, respect physical laws such as conservation of lepton number and energy and the Fermi–Dirac nature of neutrinos and (ii) to develop efficient, supercomputer-architecture-aware solution methods for the resultant nonlinear algebraic equations. In this review, we present the current state of the art in attempts to meet this challenge.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2004

Radiative Transfer Modeling of a Large Pool Fire by Discrete Ordinates, Discrete Transfer, Ray Tracing, Monte Carlo and Moment Methods

Five computational methods for solution of the radiative transfer equation in an absorbing-emitting and non-scattering gray medium were compared on a 2 m JP-8 pool fire. The temperature and absorption coefficient fields were taken from a synthetic fire due to the lack of a complete set of experimental data for fires of this size. These quantities were generated by a code that has been shown to agree well with the limited quantity of relevant data in the literature. Reference solutions to the governing equation were determined using the Monte Carlo method and a ray tracing scheme with high angular resolution. Solutions using the discrete transfer method, the discrete ordinate method (DOM) with both S(sub 4) and LC(sub 11) quadratures, and moment model using the M(sub 1) closure were compared to the reference solutions in both isotropic and anisotropic regions of the computational domain. DOM LC(sub 11) is shown to be the more accurate than the commonly used S(sub 4) quadrature technique, especially in anisotropic regions of the fire domain. This represents the first study where the M(sub 1) method was applied to a combustion problem occurring in a complex three-dimensional geometry. The M(sub 1) results agree well with other solution techniques, which is encouraging for future applications to similar problems since it is computationally the least expensive solution technique. Moreover, M(sub 1) results are comparable to DOM S(sub 4).

Jensen, K. A.↗
NASA NTRS1930

Flow and Force Equations for a Body Revolving in a Fluid

Part I gives a general method for finding the steady-flow velocity relative to a body in plane curvilinear motion, whence the pressure is found by Bernoulli's energy principle. Integration of the pressure supplies basic formulas for the zonal forces and moments on the revolving body. Part II, applying this steady-flow method, finds the velocity and pressure at all points of the flow inside and outside an ellipsoid and some of its limiting forms, and graphs those quantities for the latter forms. Part III finds the pressure, and thence the zonal force and moment, on hulls in plane curvilinear flight. Part IV derives general equations for the resultant fluid forces and moments on trisymmetrical bodies moving through a perfect fluid, and in some cases compares the moment values with those found for bodies moving in air. Part V furnishes ready formulas for potential coefficients and inertia coefficients for an ellipsoid and its limiting forms. Thence are derived tables giving numerical values of those coefficients for a comprehensive range of shapes.

Zahm, A F↗
NASA NTRS1981

Vibration-rotation intensities of SiO

Dipole moment matrix elements have been computed for a large number of transitions of astrophysical interest for the more abundant isotopes of SiO. The wave functions utilized were obtained from a direct solution of the Schroedinger equation with an accurate RKR potential. The dipole moment function, in the form of a Pade approximant, was chosen to reproduce the experimental measurements near equilibrium, to have the proper united and separated atom limits, and to have the correct long-range asymptotic dependence on internuclear separation. Because of the large number of transitions involved, and to facilitate applications, the squares of the dipole moment matrix elements were fitted by a least-squares procedure to polynomials in v and J. In addition, Einstein A coefficients are given for observed maser transitions and for selected vibrational bands. These latter are compared with previous calculations, and it is concluded that for the higher Delta v transitions, the present results represent a significant improvement.

Tipping, R. H.↗
NASA NTRS1993

A comparative study of the time-dependent standard 8-, 13- and 16-moment transport formulations of the polar wind

The ionosphere, composed of O(+), H(+), and electrons is modeled with four different transport formulations. The equations corresponding to the standard set, the 8-, 13-, and 16-moment approximations are presented, and the collision terms are expressed. Using a time-dependent technique, the ionosphere is studied between altitudes of 200 and 8600 km. The production of electrons and O(+) ions is described by a neutral atmosphere simplified photoionization scheme, while the energy deposition is supported by a downward electron heat flow of -0.005 erg/sq cm per s imposed at the topside boundary. When the models reach a steady state equilibrium, the electron solutions show differences due to the introduction of temperature anisotropies and heat flows between the components parallel and perpendicular to the magnetic field. As a corollary, the ions show structures depending on the level of approximation. A depletion of a factor of 10 is then applied to the ion densities above a certain altitude, and the development of the perturbation is followed for 1000 s for all the models.

Blelly, P. L.↗
NASA NTRS2005

Application of an Extended Parabolic Equation to the Calculation of the Mean Field and the Transverse and Longitudinal Mutual Coherence Functions Within Atmospheric Turbulence

Solutions are derived for the generalized mutual coherence function (MCF), i.e., the second order moment, of a random wave field propagating through a random medium within the context of the extended parabolic equation. Here, "generalized" connotes the consideration of both the transverse as well as the longitudinal second order moments (with respect to the direction of propagation). Such solutions will afford a comparison between the results of the parabolic equation within the pararaxial approximation and those of the wide-angle extended theory. To this end, a statistical operator method is developed which gives a general equation for an arbitrary spatial statistical moment of the wave field. The generality of the operator method allows one to obtain an expression for the second order field moment in the direction longitudinal to the direction of propagation. Analytical solutions to these equations are derived for the Kolmogorov and Tatarskii spectra of atmospheric permittivity fluctuations within the Markov approximation.

Manning, Robert M.↗
NASA NTRS1988

A closure model of diffusion transport of the Reynolds-stress equations and its application to a turbulent step flow

A study is made for the development and computations of the separating and reattaching shear flows. The highlight of the study is an attempt to predict the third-moments of turbulent velocity which is responsible for the diffusion transport of the Reynolds stresses. The present computations show that the third-moments obtained by employing the low-Reynolds number model of transport equations improve the prediction of the third-moments.

Amano, R. S.↗
NASA NTRS1993

Large eddy simulation and direct numerical simulation of high speed turbulent reacting flows

The objective of this research is to make use of Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) for the computational analyses of high speed reacting flows. Our efforts in the first phase of this research conducted within the past three years have been directed in several issues pertaining to intricate physics of turbulent reacting flows. In our previous 5 semi-annual reports submitted to NASA LaRC, as well as several technical papers in archival journals, the results of our investigations have been fully described. In this progress report which is different in format as compared to our previous documents, we focus only on the issue of LES. The reason for doing so is that LES is the primary issue of interest to our Technical Monitor and that our other findings were needed to support the activities conducted under this prime issue. The outcomes of our related investigations, nevertheless, are included in the appendices accompanying this report. The relevance of the materials in these appendices are, therefore, discussed only briefly within the body of the report. Here, results are presented of a priori and a posterior analyses for validity assessments of assumed Probability Density Function (PDF) methods as potential subgrid scale (SGS) closures for LES of turbulent reacting flows. Simple non-premixed reacting systems involving an isothermal reaction of the type A + B yields Products under both chemical equilibrium and non-equilibrium conditions are considered. A priori analyses are conducted of a homogeneous box flow, and a spatially developing planar mixing layer to investigate the performance of the Pearson Family of PDF's as SGS models. A posteriori analyses are conducted of the mixing layer using a hybrid one-equation Smagorinsky/PDF SGS closure. The Smagorinsky closure augmented by the solution of the subgrid turbulent kinetic energy (TKE) equation is employed to account for hydrodynamic fluctuations, and the PDF is employed for modeling the effects of scalar fluctuations. The implementation of the model requires the knowledge of the local values of the first two SGS moments. These are provided by additional modeled transport equations. In both a priori and a posteriori analyses, the predicted results are appraised by comparison with subgrid averaged results generated by DNS. Based on these results, the paths to be followed in future investigations are identified.

Adumitroaie, V.↗
NASA NTRS1974

Reaction formulation for radiation and scattering from plates, corner reflectors and dielectric-coated cylinders

The reaction concept is employed to formulate an integral equation for radiation and scattering from plates, corner reflectors, and dielectric-coated conducting cylinders. The surface-current density on the conducting surface is expanded with subsectional bases. The dielectric layer is modeled with polarization currents radiating in free space. Maxwell's equation and the boundary conditions are employed to express the polarization-current distribution in terms of the surface-current density on the conducting surface. By enforcing reaction tests with an array of electric test sources, the moment method is employed to reduce the integral equation to a matrix equation. Inversion of the matrix equation yields the current distribution, and the scattered field is then obtained by integrating the current distribution. The theory, computer program and numerical results are presented for radiation and scattering from plates, corner reflectors, and dielectric-coated conducting cylinders.

Wang, N. N.↗
NASA NTRS1994

Modeling of aircraft unsteady aerodynamic characteristics. Part 1: Postulated models

A short theoretical study of aircraft aerodynamic model equations with unsteady effects is presented. The aerodynamic forces and moments are expressed in terms of indicial functions or internal state variables. The first representation leads to aircraft integro-differential equations of motion; the second preserves the state-space form of the model equations. The formulations of unsteady aerodynamics is applied in two examples. The first example deals with a one-degree-of-freedom harmonic motion about one of the aircraft body axes. In the second example, the equations for longitudinal short-period motion are developed. In these examples, only linear aerodynamic terms are considered. The indicial functions are postulated as simple exponentials and the internal state variables are governed by linear, time-invariant, first-order differential equations. It is shown that both approaches to the modeling of unsteady aerodynamics lead to identical models.

Klein, Vladislav↗
NASA NTRS1956

Theoretical span load distributions and rolling moments for sideslipping wings of arbitrary plan form in incompressible flow

A method of computing span loads and the resulting rolling moments for sideslipping wings of arbitrary plan form in incompressible flow is derived. The method requires that the span load at zero sideslip be known for the wing under consideration. Because this information is available for a variety of wings, this requirement should not seriously restrict the application of the present method. The basic method derived herein requires a mechanical differentiation and integration to obtain the rolling moment for the general wing in sideslip. For wings having straight leading and trailing edges over each semispan, the rolling moment due to sideslip is given by a simple equation in terms of plan-form parameters and the lateral center of pressure of the lift due to angle of attack. Calculated span loads and rolling-moment parameters are compared with experimental values. The comparison indicates good agreement between calculations and available experimental data.

Queijo, M J↗
NASA NTRS1974

Approach to the origin of turbulence on the basis of two-point kinetic theory

Equations for the fluctuation correlation in an incompressible shear flow are derived on the basis of kinetic theory, utilizing the two-point distribution function which obeys the BBGKY hierarchy equation truncated with the hypothesis of 'ternary' molecular chaos. The step from the molecular to the hydrodynamic description is accomplished by a moment expansion which is a two-point version of the thirteen-moment method, and which leads to a series of correlation equations, viz., the two-point counterparts of the continuity equation, the Navier-Stokes equation, etc. For almost parallel shearing flows the two-point equation is separable and reduces to two Orr-Sommerfeld equations with different physical implications.

Tsuge, S.↗
NASA NTRS2005

Diffraction of Harmonic Flexural Waves in a Cracked Elastic Plate Carrying Electrical Current

The scattering effect of harmonic flexural waves at a through crack in an elastic plate carrying electrical current is investigated. In this context, the Kirchhoffean bending plate theory is extended as to include magnetoelastic interactions. An incident wave giving rise to bending moments symmetric about the longitudinal z-axis of the crack is applied. Fourier transform technique reduces the problem to dual integral equations, which are then cast to a system of two singular integral equations. Efficient numerical computation is implemented to get the bending moment intensity factor for arbitrary frequency of the incident wave and of arbitrary electrical current intensity. The asymptotic behaviour of the bending moment intensity factor is analysed and parametric studies are conducted.

Ambur, Damodar R.↗
NASA NTRS1981

Computation of the modes and polar moment of inertial of the blades of an HAWT

The coupled differential equations of motion of the blades of a horizontal axis wind turbine are solved numerically, permitting the optimization of the design at relatively low cost. The equation of motion is transformed into a set of first order equations and solved with fourth order Runge-Kutta integrators. This technique is applied to a twisted, tapered blade of variable cross section and stiffness including discontinuities. The first six natural frequencies and mode shapes are obtained. The polar moment of inertia of the blades is obtained as a function of frequency and rotational speed.

Beaulieu, G.↗