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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Prediction of an axisymmetric combusting flow

A numerical model for turbulent, recirculating combusting flow is developed and applied to a research combustor. The model is based on the time-averaged Navier-Stokes equations with k-epsilon turbulence closure and the compositional fluctuations at each point are given probabilistically in terms of the mixture fraction. The probability density function is derived from transport equations for its first two moments along with an assumption regarding its shape. The resulting equations are solved using a standard line relaxation algorithm. It is found that the predictions of the model are in good agreement with data from an axisymmetric, bluff-body stabilized research combustor, while the major discrepancies are similar to those found in isothermal flow comparisons. The peculiar features of this flow which contribute to the errors are examined. The agreement between the theory and data deteriorates as the central jet velocity is increased which indicates an enhanced role for unsteady effects.

Correa, S. M.↗
NASA NTRS2000

Theory of the lattice Boltzmann Method: Dispersion, Dissipation, Isotropy, Galilean Invariance, and Stability

The generalized hydrodynamics (the wave vector dependence of the transport coefficients) of a generalized lattice Boltzmann equation (LBE) is studied in detail. The generalized lattice Boltzmann equation is constructed in moment space rather than in discrete velocity space. The generalized hydrodynamics of the model is obtained by solving the dispersion equation of the linearized LBE either analytically by using perturbation technique or numerically. The proposed LBE model has a maximum number of adjustable parameters for the given set of discrete velocities. Generalized hydrodynamics characterizes dispersion, dissipation (hyper-viscosities), anisotropy, and lack of Galilean invariance of the model, and can be applied to select the values of the adjustable parameters which optimize the properties of the model. The proposed generalized hydrodynamic analysis also provides some insights into stability and proper initial conditions for LBE simulations. The stability properties of some 2D LBE models are analyzed and compared with each other in the parameter space of the mean streaming velocity and the viscous relaxation time. The procedure described in this work can be applied to analyze other LBE models. As examples, LBE models with various interpolation schemes are analyzed. Numerical results on shear flow with an initially discontinuous velocity profile (shock) with or without a constant streaming velocity are shown to demonstrate the dispersion effects in the LBE model; the results compare favorably with our theoretical analysis. We also show that whereas linear analysis of the LBE evolution operator is equivalent to Chapman-Enskog analysis in the long wave-length limit (wave vector k = 0), it can also provide results for large values of k. Such results are important for the stability and other hydrodynamic properties of the LBE method and cannot be obtained through Chapman-Enskog analysis.

Lallemand, Pierre↗
NASA NTRS1975

Structural Properties and Estimation of Delay Systems

Two areas in the theory of delay systems were studied: structural properties and their applications to feedback control, and optimal linear and nonlinear estimation. The concepts of controllability, stabilizability, observability, and detectability were investigated. The property of pointwise degeneracy of linear time-invariant delay systems is considered. Necessary and sufficient conditions for three dimensional linear systems to be made pointwise degenerate by delay feedback were obtained, while sufficient conditions for this to be possible are given for higher dimensional linear systems. These results were applied to obtain solvability conditions for the minimum time output zeroing control problem by delay feedback. A representation theorem is given for conditional moment functionals of general nonlinear stochastic delay systems, and stochastic differential equations are derived for conditional moment functionals satisfying certain smoothness properties.

Kwong, R. H. S.↗
NASA NTRS1989

Solutions to bi-Maxwellian transport equations for the polar wind

In this study, polar wind solutions are obtained for a broad range of O(+) density, H(+) drift velocity, electron temperature and H(+) temperature boundary conditions. The bi-Maxwellian-based 16-moment set of transport equations is used, since this set is expected to be superior to Maxwellian-based equations in describing large temperature anisotropies and heat flows. The present solutions corroborate earlier results when similar boundary conditions are used. Also, for previously unexplored combinations of boundary conditions, the present solutions are often qualitatively different from any obtained before.

Demars, H. G.↗
NASA NTRS1978

Extensions of guiding center motion to higher order

In a static magnetic field, some well-known guiding-center equations maintain their form when extended to next order in gyroradius. In these cases, it is only necessary to include the next order term in the magnetic moment series. The differential equation for guiding-center motion which describes both the parallel and perpendicular velocities correctly through first order in gyroradius is given. The question of how to define the guiding center position through second order arises and is discussed, and second order drifts are derived for one usual definition. The toroidal canonical angular momentum, P-phi, of the guiding center in an axisymmetric field is shown to be conserved using the guiding center velocity correct through first order. When second-order motion is included, P-phi is no longer a constant. The above extensions of guiding-center theory help to resolve the different tokamak orbits obtained either by using the guiding-center equations of motion or by using conservation of P-phi.

Northrop, T. G.↗
NASA NTRS1988

Preshock region acceleration of implanted cometary H(+) and O(+)

A self-consistent, three-fluid model of plasma transport and implanted ion acceleration in the unshocked solar wind is presented. The solar wind plasma is depleted by charge exchange with the expanding cometary exosphere, while implanted protons and heavy ions are produced by photoionization and charge transfer and lost by charge exchange. A generalized transport equation describing convection, adiabatic and diffusive velocity change, and the appropriate production terms is used to describe the evolution of the two cometary ion components, while the moments of the Boltzmann equation are used to calculate the solar wind density and pressure. The flow velocity is obtained self-consistently by combining the conservation equations of the three ion species. The results imply that second-order Fermi acceleration can explain the implanted spectra observed in the unshocked solar wind. Comparison of measured and calculated distribution indicates that spatial diffusion of implanted ions probably plays an important role in forming the energetic particle environment in the shock vicinity.

Gombosi, Tamas I.↗
NASA NTRS1951

Equations and charts for the rapid estimation of hinge-moment and effectiveness parameters for trailing-edge controls having leading and trailing edges swept ahead of the Mach lines

Existing conical-flow solutions have been used to calculate the hinge-moments and effectiveness parameters of trailing-edge controls having leading and trailing edges swept ahead of the Mach lines and having streamwise root and tip chords. Equations and detailed charts are presented for the rapid estimation of these parameters. Also included is an approximate method by which these parameters may be corrected for airfoil-section thickness.

Goin, Kennith L↗
NASA NTRSDate not supplied

Unsteady Aerodynamic Modeling of Atmospheric Entry Vehicles in Subsonic and Incompressible flow: A Frequency Response Approach

The determination of time-varying lift force and pitch moment generated by a purely pitching Earth-entry capsule is investigated. Experiments were conducted in the 12-foot Low-Speed Tunnel at NASA Langley Research Center, testing a range of oscillation frequencies at a zero-mean angle of attack with a pitching amplitude of 10 degrees. In light of these measurements, a closed-form set of analytically derived equations for lift and moment was used to develop a semi-empirical formulation, incorporating empirically determined values from the experimental runs. The equations are grounded in potential flow theory, Theodorsen's classical theory of unsteady aerodynamics, and the Joukowski theorem of conformal mapping. The unsteady aerodynamics generated by the oscillating body are then modeled by constructing frequency response functions, with quasi-steady forces and moments serving as inputs and unsteady forces and moments as outputs. The experimentally determined gain and phase variations characterize the unsteady nature of the flow and the system's response and flow time-lag to input flow parameters for a blunt-body entry vehicle. The final semi-empirical model is validated with a set of parameters beyond the initial test matrix.

Entry Systems Modeling↗
NASA NTRSDate not supplied

Unsteady Aerodynamic Modeling of Atmospheric Entry Vehicles in Subsonic and Incompressible Flow: A Frequency Response Approach

The determination of time-varying lift force and pitch moment generated by a purely pitching Earth-entry capsule is investigated. Experiments were conducted in the 12-foot Low-Speed Tunnel at NASA Langley Research Center, testing a range of oscillation frequencies at a zero-mean angle of attack with a pitching amplitude of 10 degrees. In light of these measurements, a closed-form set of analytically derived equations for lift and moment was used to develop a semi-empirical formulation, incorporating empirically determined values from the experimental runs. The equations are grounded in potential flow theory, Theodorsen's classical theory of unsteady aerodynamics, and the Joukowski theorem of conformal mapping. The unsteady aerodynamics generated by the oscillating body are then modeled by constructing frequency response functions, with quasi-steady forces and moments serving as inputs and unsteady forces and moments as outputs. The experimentally determined gain and phase variations characterize the unsteady nature of the flow and the system's response and flow time-lag to input flow parameters for a blunt-body entry vehicle. The final semi-empirical model is validated with a set of parameters beyond the initial test matrix.

Entry Systems Modeling↗
DOE OSTI2021

Formulation of 8-moment plasma transport with application to the Nernst effect

8-moment plasma models using two different distribution functions are used to study the Nernst effect and heat transport in dense plasma. These models are presented in hyperbolic form in contrast to traditional parabolic systems derived from perturbing the distribution function, as in Braginskii [Rev. Plasma Phys. 1, 205 (1965)]. The hyperbolic moment formulation can be solved implicitly in time with straightforward and fast local solvers. The numerical implementation of 8-moment models with the relaxation method in the PERSEUS code is also presented. To test 8-moment PERSEUS compared to Braginskii’s transport equations, a verification test for the Nernst thermomagnetic wave by Velikovich et al. [Phys. Plasmas 26, 112702 (2019)] is performed that confirms the presence of the same physics, but with slight differences in the transport coefficients, which are tabulated in the limits of high and low magnetization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1988

Creep rupture analysis of a beam resting on high temperature foundation

A simplified uniaxial strain controlled creep damage law is deduced with the use of experimental observation from a more complex strain dependent law. This creep damage law correlates the creep damage, which is interpreted as the density variation in the material, directly with the accumulated creep strain. Based on the deduced uniaxial strain controlled creep damage law, a continuum mechanical creep rupture analysis is carried out for a beam resting on a high temperature elastic (Winkler) foundation. The analysis includes the determination of the nondimensional time for initial rupture, the propagation of the rupture front with the associated thinning of the beam, and the influence of creep damage on the deflection of the beam. Creep damage starts accumulating in the beam as soon as the load is applied, and a creep rupture front develops at and propagates from the point at which the creep damage first reaches its critical value. By introducing a series of fundamental assumptions within the framework of technical Euler-Bernoulli type beam theory, a governing set of integro-differential equations is derived in terms of the nondimensional bending moment and the deflection. These governing equations are subjected to a set of interface conditions at the propagating rupture front. A numerical technique is developed to solve the governing equations together with the interface equations, and the computed results are presented and discussed in detail.

Gu, Randy J.↗
NASA NTRS1995

A new approach to the formulation of scalar flux closure

This report shows that if a stochastic differential equation (Langevin equation) for velocity fluctuation vector is known, it is possible to derive the equations for scalar flux transport. Durbin and Speziale (1994) showed that the second moment of this stochastic differential equation gives an equation for the evolution of Reynolds stress tensor. Similarly, the stochastic equation will give an equation for scalar flux. Therefore, a coupling between these two is present. The basis for the present work is that there should be Langevin equations that can produce acceptable models for both the Reynolds stress tensor and the scalar flux vector. Having found this basic Langevin equation, the amount of work needed to model the second order closure problems is reduced; using the well developed models for Reynolds stress equations, it will be possible to derive corresponding models for scalar flux equation.

Shabany, Y.↗
NASA NTRS2005

Deflection-Based Aircraft Structural Loads Estimation with Comparison to Flight

Traditional techniques in structural load measurement entail the correlation of a known load with strain-gage output from the individual components of a structure or machine. The use of strain gages has proved successful and is considered the standard approach for load measurement. However, remotely measuring aerodynamic loads using deflection measurement systems to determine aeroelastic deformation as a substitute to strain gages may yield lower testing costs while improving aircraft performance through reduced instrumentation weight. With a reliable strain and structural deformation measurement system this technique was examined. The objective of this study was to explore the utility of a deflection-based load estimation, using the active aeroelastic wing F/A-18 aircraft. Calibration data from ground tests performed on the aircraft were used to derive left wing-root and wing-fold bending-moment and torque load equations based on strain gages, however, for this study, point deflections were used to derive deflection-based load equations. Comparisons between the strain-gage and deflection-based methods are presented. Flight data from the phase-1 active aeroelastic wing flight program were used to validate the deflection-based load estimation method. Flight validation revealed a strong bending-moment correlation and slightly weaker torque correlation. Development of current techniques, and future studies are discussed.

Lizotte, Andrew M.↗
NASA NTRS2005

Deflection-Based Structural Loads Estimation From the Active Aeroelastic Wing F/A-18 Aircraft

Traditional techniques in structural load measurement entail the correlation of a known load with strain-gage output from the individual components of a structure or machine. The use of strain gages has proved successful and is considered the standard approach for load measurement. However, remotely measuring aerodynamic loads using deflection measurement systems to determine aeroelastic deformation as a substitute to strain gages may yield lower testing costs while improving aircraft performance through reduced instrumentation weight. This technique was examined using a reliable strain and structural deformation measurement system. The objective of this study was to explore the utility of a deflection-based load estimation, using the active aeroelastic wing F/A-18 aircraft. Calibration data from ground tests performed on the aircraft were used to derive left wing-root and wing-fold bending-moment and torque load equations based on strain gages, however, for this study, point deflections were used to derive deflection-based load equations. Comparisons between the strain-gage and deflection-based methods are presented. Flight data from the phase-1 active aeroelastic wing flight program were used to validate the deflection-based load estimation method. Flight validation revealed a strong bending-moment correlation and slightly weaker torque correlation. Development of current techniques, and future studies are discussed.

Lizotte, Andrew M.↗
NASA NTRS1982

Particle trajectory computer program for icing analysis of axisymmetric bodies

General aviation aircraft and helicopters exposed to an icing environment can accumulate ice resulting in a sharp increase in drag and reduction of maximum lift causing hazardous flight conditions. NASA Lewis Research Center (LeRC) is conducting a program to examine, with the aid of high-speed computer facilities, how the trajectories of particles contribute to the ice accumulation on airfoils and engine inlets. This study, as part of the NASA/LeRC research program, develops a computer program for the calculation of icing particle trajectories and impingement limits relative to axisymmetric bodies in the leeward-windward symmetry plane. The methodology employed in the current particle trajectory calculation is to integrate the governing equations of particle motion in a flow field computed by the Douglas axisymmetric potential flow program. The three-degrees-of-freedom (horizontal, vertical, and pitch) motion of the particle is considered. The particle is assumed to be acted upon by aerodynamic lift and drag forces, gravitational forces, and for nonspherical particles, aerodynamic moments. The particle momentum equation is integrated to determine the particle trajectory. Derivation of the governing equations and the method of their solution are described in Section 2.0. General features, as well as input/output instructions for the particle trajectory computer program, are described in Section 3.0. The details of the computer program are described in Section 4.0. Examples of the calculation of particle trajectories demonstrating application of the trajectory program to given axisymmetric inlet test cases are presented in Section 5.0. For the examples presented, the particles are treated as spherical water droplets. In Section 6.0, limitations of the program relative to excessive computer time and recommendations in this regard are discussed.

Frost, Walter↗
NASA NTRS1986

Group-kinetic theory of turbulence

The two phases are governed by two coupled systems of Navier-Stokes equations. The couplings are nonlinear. These equations describe the microdynamical state of turbulence, and are transformed into a master equation. By scaling, a kinetic hierarchy is generated in the form of groups, representing the spectral evolution, the diffusivity and the relaxation. The loss of memory in formulating the relaxation yields the closure. The network of sub-distributions that participates in the relaxation is simulated by a self-consistent porous medium, so that the average effect on the diffusivity is to make it approach equilibrium. The kinetic equation of turbulence is derived. The method of moments reverts it to the continuum. The equation of spectral evolution is obtained and the transport properties are calculated. In inertia turbulence, the Kolmogoroff law for weak coupling and the spectrum for the strong coupling are found. As the fluid analog, the nonlinear Schrodinger equation has a driving force in the form of emission of solitons by velocity fluctuations, and is used to describe the microdynamical state of turbulence. In order for the emission together with the modulation to participate in the transport processes, the non-homogeneous Schrodinger equation is transformed into a homogeneous master equation. By group-scaling, the master equation is decomposed into a system of transport equations, replacing the Bogoliubov system of equations of many-particle distributions. It is in the relaxation that the memory is lost when the ensemble of higher-order distributions is simulated by an effective porous medium. The closure is thus found. The kinetic equation is derived and transformed into the equation of spectral flow.

Tchen, C. M.↗
NASA NTRS1949

Response of a Rotating Propeller to Aerodynamic Excitation

The flexural vibration of a rotating propeller blade with clamped shank is analyzed with the object of presenting, in matrix form, equations for the elastic bending moments in forced vibration resulting from aerodynamic forces applied at a fixed multiple of rotational speed. Matrix equations are also derived which define the critical speeds end mode shapes for any excitation order and the relation between critical speed and blade angle. Reference is given to standard works on the numerical solution of matrix equations of the forms derived. The use of a segmented blade as an approximation to a continuous blade provides a simple means for obtaining the matrix solution from the integral equation of equilibrium, so that, in the numerical application of the method presented, the several matrix arrays of the basic physical characteristics of the propeller blade are of simple form, end their simplicity is preserved until, with the solution in sight, numerical manipulations well-known in matrix algebra yield the desired critical speeds and mode shapes frame which the vibration at any operating condition may be synthesized. A close correspondence between the familiar Stodola method and the matrix method is pointed out, indicating that any features of novelty are characteristic not of the analytical procedure but only of the abbreviation, condensation, and efficient organization of the numerical procedure made possible by the use of classical matrix theory.

Arnoldi, Walter E.↗