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 523 records · Page 29

A Theoretical Study of the Angular Motions of Spinning Bodies in Space

Euler's dynamic equations were linearized and solved analytically. Analytical expressions which relate angular motions to spin-rate and inertia distributions were obtained and found to be in good agreement with machine solutions of the nonlinear equations for the case of a rectangular-pulse pitching-moment disturbance. Consideration was given to the effects produced by having artificial damping in the system.

Suddath, Jerrold H.↗

Strain Gage Loads Calibration Testing of the Active Aeroelastic Wing F/A-18 Aircraft

This report describes strain-gage calibration loading through the application of known loads of the Active Aeroelastic Wing F/A-18 airplane. The primary goal of this test is to produce a database suitable for deriving load equations for left and right wing root and fold shear; bending moment; torque; and all eight wing control-surface hinge moments. A secondary goal is to produce a database of wing deflections measured by string potentiometers and the onboard flight deflection measurement system. Another goal is to produce strain-gage data through both the laboratory data acquisition system and the onboard aircraft data system as a check of the aircraft system. Thirty-two hydraulic jacks have applied loads through whiffletrees to 104 tension-compression load pads bonded to the lower wing surfaces. The load pads covered approximately 60 percent of the lower wing surface. A series of 72 load cases has been performed, including single-point, double-point, and distributed load cases. Applied loads have reached 70 percent of the flight limit load. Maximum wingtip deflection has reached nearly 16 in.

Lokos, William A.↗

Application of CFE/POST2 for Simulation of Launch Vehicle Stage Separation

The constraint force equation (CFE) methodology provides a framework for modeling constraint forces and moments acting at joints that connect multiple vehicles. With implementation in Program to Optimize Simulated Trajectories II (POST 2), the CFE provides a capability to simulate end-to-end trajectories of launch vehicles, including stage separation. In this paper, the CFE/POST2 methodology is applied to the Shuttle-SRB separation problem as a test and validation case. The CFE/POST2 results are compared with STS-1 flight test data.

Pamadi, Bandu N.↗

A Stochastic Fractional Dynamics Model of Space-time Variability of Rain

Rainfall varies in space and time in a highly irregular manner and is described naturally in terms of a stochastic process. A characteristic feature of rainfall statistics is that they depend strongly on the space-time scales over which rain data are averaged. A spectral model of precipitation has been developed based on a stochastic differential equation of fractional order for the point rain rate, that allows a concise description of the second moment statistics of rain at any prescribed space-time averaging scale. The model is thus capable of providing a unified description of the statistics of both radar and rain gauge data. The underlying dynamical equation can be expressed in terms of space-time derivatives of fractional orders that are adjusted together with other model parameters to fit the data. The form of the resulting spectrum gives the model adequate flexibility to capture the subtle interplay between the spatial and temporal scales of variability of rain but strongly constrains the predicted statistical behavior as a function of the averaging length and times scales. We test the model with radar and gauge data collected contemporaneously at the NASA TRMM ground validation sites located near Melbourne, Florida and in Kwajalein Atoll, Marshall Islands in the tropical Pacific. We estimate the parameters by tuning them to the second moment statistics of radar data. The model predictions are then found to fit the second moment statistics of the gauge data reasonably well without any further adjustment.

TRMM↗

Birkeland currents in an anisotropic, magnetostatic plasma

The paper derives an expression for the parallel current density for a plasma characterized by negligible bulk flow (magnetostatic velocity and a two-component (anisotropic) pressure tensor by expanding the equilibrium Vlasov equation for each species in the adiabatic parameter until such point as a nonvanishing moment j-parallel = integral d-cubed vv-parallel f is identified. The result is a nonlocal one: it relates j-parallel at one point s along a field line to j-parallel at another (reference) point s0 plus an integral function of the pressure and magnetic field between them. The equation derived by Vasyliunas (1970) is generalized and extended to a plasma in which the pressure tensor is comprised of two separate elements, P-perpendicular and P-parallel. This equation follows when P-perpendicular is set equal to P-parallel and s and s0 are taken to be at the ionosphere and the equator. These results are compared to others in the literature.

Birmingham, Thomas J.↗

Tidal-friction theory of the earth-moon system

Serious errors contained in Jeffreys' (1952, 1959, 1970, 1976) discussion of tidal friction in the earth-moon system are identified and their consequences are discussed. A direct solution of the dynamical tidal equations for the couple from the earth acting upon the moon and the couple from the earth acting upon the sun, which were left unsolved by Jeffreys, is found to be incompatible with observations and the predictions of linear or quadratic friction theory, due to his failure to take into account the possible change of the moment of inertia of the earth with time in the derivation of the dynamical equations. Consideration of this factor leads to the conclusion that the earth must be contracting at a rate of 14.7 x 10 to the -11th/year, which can be accounted for only by the Ramsey theory, in which the terrestrial core is considered as a phase change rather than a change in chemical composition. Implications of this value for the rates of changes in day length and lunar distance are also indicated.

Lyttleton, R. A.↗

Macroscopic Models of Radiative Transfer as Applied to Computation of the Radiation Field in the Solar Atmosphere

In this paper a new three dimensional half-moment model for radiative transfer is presented for a gray medium. It describes the evolution of the zeroth and first directional half moments of the radiative intensity. The closure is provided, similarly to Dubroca and Klar, by the maximum entropy concept. This work generalizes that model to three dimensions. The model presented here (the derivation being done in Ripoll and Wray, called the M(sup 1/2)(sub 1) model, is a hyperbolic system consisting of a total of eight equations in three dimensions, four equations for each direction. Each half model has the classical form of a macroscopic moment model in which the pressure tensor is constructed from the well-known Eddington tensor with a particular Eddington factor. Moreover, different source and border terms occur. The latter introduce couplings between the macroscopic and microscopic quantities and between the + and - streams, through the intensity in the plane perpendicular to the flux. The main theoretical application of the half moment model, treated in this paper, is its reduction to a full moment model, called M(sup +)(sub 1), for the particular but important case of a hot, opaque source radiating in a cold transparent (or semi-transparent) medium for very specific applications, such as stellar interiors or atmospheres, or combustion problems. The structure of the paper is as follows. In section 2, the model M(sup 1/2)(sub 1) is presented. In section 3, for the particular case of a hot, opaque source radiating into a cold medium, the half moment model is reduced to the M(sup +)(sub 1) model. In section 4, we first solve a simple and academic problem to validate the models, followed by a simplified solar atmosphere.

Ripoll, J.-F.↗

Dynamic Stability Analysis of Blunt Body Entry Vehicles Using Time-Lagged Aftbody Pitching Moments

This analysis defines an analytic model for the pitching motion of blunt bodies during atmospheric entry. The proposed model is independent of the pitch damping sum coefficient present in the standard formulation of the equations of motion describing pitch oscillations of a decelerating blunt body, instead using the principle of a time-lagged aftbody moment as the forcing function for oscillation divergence. Four parameters, all with intuitive physical relevance, are introduced to fully define the aftbody moment and the associated time delay. It is shown that the dynamic oscillation responses typical to blunt bodies can be produced using hysteresis of the aftbody moment in place of the pitch damping coefficient. The approach used in this investigation is shown to be useful in understanding the governing physical mechanisms for blunt body dynamic stability and in guiding vehicle and mission design requirements. A validation case study using simulated ballistic range test data is conducted. From this, parameter identification is carried out through the use of a least squares optimizing routine. Results show good agreement with the limited existing literature for the parameters identified, suggesting that the model proposed could be validated by an experimental ballistic range test series. The trajectories produced by the identified parameters were found to match closely those from the MER ballistic range tests for a wide array of initial conditions and can be identified with a reasonable number of ballistic range shots and computational effort.

Supersonic↗

A Method of Moments Wide Band Adaptive Rational Interpolation Method for High-Quality Factor Resonant Cavities

A new adaptive rational interpolation method is proposed to obtain the wideband frequency response of a resonant cavity simulated with the method of moments (MoM). This interpolation method uses both the Loewner matrix to construct a rational expression for the solution vector of MoM’s matrix system and an error estimator generated by the solution vectors and their derivatives. This error estimator is implemented in the adaptive procedure to gain a minimum set of frequencies and solution vectors required in the interpolation. The resulting set of frequencies and solution vectors is applied to interpolate other system variables, such as shielding effectiveness and input impedance. Here numerical results of a slotted cylindrical cavity supporting high-quality factor resonances are presented, showing that the new rational interpolation method is accurate and efficient in interpolating the complicated resonant response of the solution vector functions.

42 ENGINEERING↗

ISS Double-Gimbaled CMG Subsystem Simulation Using the Agile Development Method

This paper presents an evolutionary approach in simulating a cluster of 4 Control Moment Gyros (CMG) on the International Space Station (ISS) using a common sense approach (the agile development method) for concurrent mathematical modeling and simulation of the CMG subsystem. This simulation is part of Training systems for the 21st Century simulator which will provide training for crew members, instructors, and flight controllers. The basic idea of how the CMGs on the space station are used for its non-propulsive attitude control is briefly explained to set up the context for simulating a CMG subsystem. Next different reference frames and the detailed equations of motion (EOM) for multiple double-gimbal variable-speed control moment gyroscopes (DGVs) are presented. Fixing some of the terms in the EOM becomes the special case EOM for ISS's double-gimbaled fixed speed CMGs. CMG simulation development using the agile development method is presented in which customer's requirements and solutions evolve through iterative analysis, design, coding, unit testing and acceptance testing. At the end of the iteration a set of features implemented in that iteration are demonstrated to the flight controllers thus creating a short feedback loop and helping in creating adaptive development cycles. The unified modeling language (UML) tool is used in illustrating the user stories, class designs and sequence diagrams. This incremental development approach of mathematical modeling and simulating the CMG subsystem involved the development team and the customer early on, thus improving the quality of the working CMG system in each iteration and helping the team to accurately predict the cost, schedule and delivery of the software.

Inampudi, Ravi↗

Analysis of aerodynamic coefficients using gradient data: Spanwise turbulence effects on airplane response

The influence of spanwise turbulence on airplane dynamic behavior is determined. Calculations are based on data collected from test flights with the NASA B-57 research aircraft. The approach is to first compute aerodynamic forces and moments due to a spanwise distribution of angle of attack and airspeed. Secondly, these quantities are incorporated into the equations of motion. Simulation of flights done with the effects of spanwise turbulence included are compared to simulations without any spanwise turbulence. The findings of the study are that the moments developed by turbulence along the span are significant and that more realistic flight simulation can be achieved by including the spanwise turbulence terms.

Ringnes, E. A.↗

A moment-conserving discontinuous Galerkin representation of the relativistic Maxwellian distribution

Kinetic simulations of relativistic gases and plasmas are critical for understanding diverse astrophysical and terrestrial systems, but the accurate construction of the relativistic Maxwellian, the Maxwell–Jüttner distribution, on a discrete simulation grid is challenging. Difficulties arise from the finite velocity bounds of the domain, which may not capture the entire distribution function, as well as errors introduced by projecting the function onto a discrete grid. Here, we present a novel scheme for iteratively correcting the moments of the projected distribution applicable to all grid-based discretizations of the relativistic kinetic equation. In addition, we describe how to compute the needed nonlinear quantities, such as Lorentz boost factors, in a discontinuous Galerkin scheme through a combination of numerical quadrature and weak operations. The resulting method accurately captures the distribution function and ensures that the moments match the desired values to machine precision.

astrophysical plasmas↗

Neutron scattering off Weyl semimetals

We present how to detect type-I Weyl nodes in a material by inelastic neutron scattering. Such an experiment first of all allows one to determine the dispersion of the Weyl fermions. We extend the reasoning to produce a quantitative test of the Weyl equation, taking into account realistic anisotropic properties. These anisotropies are mostly contained in the form of the emergent magnetic moment of the excitations, which determines how they couple to the neutrons. Although there are many material parameters, we find several quantitative predictions that are universal and demonstrate that the excitations are described by solutions to the Weyl equation. The anisotropic coupling between electrons and neutrons implies that even fully unpolarized neutrons can reveal the spin-momentum locking of the Weyl fermions because the neutrons will couple to some components of the Weyl fermion pseudospin more strongly. On the other hand, in an experiment with polarized neutrons, the scattered neutron beam remains fully polarized in a direction that varies as a function of momentum transfer (within the range of validity of the Weyl equation). This allows measurement of the chirality of Weyl fermions for inversion-symmetric nodes. Furthermore, in this work we estimate that the scattering rate may be large enough for such experiments to be practical; in particular, the magnetic moment may be larger than the ordinary Bohr magneton, compensating for a small density of states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A B-B-G-K-Y framework for fluid turbulence

A kinetic theory for fluid turbulence is developed from the Liouville equation and the associated BBGKY hierarchy. Real and imaginary parts of Fourier coefficients of fluid variables play the roles of particles. Closure is achieved by the assumption of negligible five-coefficient correlation functions and probability distributions of Fourier coefficients are the basic variables of the theory. An additional approximation leads to a closed-moment description similar to the so-called eddy-damped Markovian approximation. A kinetic equation is derived for which conservation laws and an H-theorem can be rigorously established, the H-theorem implying relaxation of the absolute equilibrium of Kraichnan. The equation can be cast in the Fokker-Planck form, and relaxation times estimated from its friction and diffusion coefficients. An undetermined parameter in the theory is the free decay time for triplet correlations. Some attention is given to the inclusion of viscous damping and external driving forces.

Montgomery, D.↗

Electromagnetic scattering by coated convex surfaces and wedges simulated by approximate boundary conditions

Asymptotic/high-frequency solutions are developed for analyzing the non-specular scattering mechanisms associated with coated convex surfaces and edges simulated by approximate boundary conditions. In particular, the standard impedance boundary conditions (SIBC's) and the second order generalized impedance boundary conditions (GIBC's) are employed for a characterization of the edge diffraction, creeping wave, and surface diffracted wave contributions. To study the creeping wave and surface diffracted wave mechanisms, rigorous UTD (uniform geometrical theory of diffraction) diffraction coefficients are developed for a convex coated cylinder simulated with SIBC's and GIBC's. The ray solutions obtained remain valid in the transition region and reduce uniformly to those in the deep lit and shadow regions. A uniform asymptotic solution is also presented for observations in the close vicinity of the cylinder. The diffraction coefficient for a convex cylinder are obtained via a generalization of the corresponding ones of the circular cylinder. To validate the asymptotic/high-frequency solution, integral equations are derived for both E and H-polarization and solved numerically using the method of moments. Results are presented for a single and three layered coated convex cylinder. Some insights are also provided on the accuracy of the employed GIBC's versus SIBC's for application to curved surfaces. To characterize the scattering by impedance wedges illuminated at skew incidence, diffraction coefficients are derived from an approximate solution of the governing functional difference equations. This solution exactly recovers the known ones for an impedance half plane or an arbitrary wedge at normal incidence, and to validate it for other wedge angles, a moment method code was used. Finally, to test the usefulness of the approximate skew incidence impedance wedge diffraction coefficient for three dimensional structures, equivalent currents are derived in the context of the physical theory of diffraction (PTD) for a finite length impedance wedge of arbitrary internal angles. These are incorporated in a standard general purpose PTD code and results are presented for a number of different impedance structures.

Syed, H. H.↗

TAILSIM Users Guide

The TAILSIM program uses a 4th order Runge-Kutta method to integrate the standard aircraft equations-of-motion (EOM). The EOM determine three translational and three rotational accelerations about the aircraft's body axis reference system. The forces and moments that drive the EOM are determined from aerodynamic coefficients, dynamic derivatives, and control inputs. Values for these terms are determined from linear interpolation of tables that are a function of parameters such as angle-of-attack and surface deflections. Buildup equations combine these terms and dimensionalize them to generate the driving total forces and moments. Features that make TAILSIM applicable to studies of tailplane stall include modeling of the reversible control System, modeling of the pilot performing a load factor and/or airspeed command task, and modeling of vertical gusts. The reversible control system dynamics can be described as two hinged masses connected by a spring. resulting in a fifth order system. The pilot model is a standard form of lead-lag with a time delay applied to an integrated pitch rate and/or airspeed error feedback. The time delay is implemented by a Pade approximation, while the commanded pitch rate is determined by a commanded load factor. Vertical gust inputs include a single 1-cosine gust and a continuous NASA Dryden gust model. These dynamic models. coupled with the use of a nonlinear database, allow the tailplane stall characteristics, elevator response, and resulting aircraft response, to be modeled. A useful output capability of the TAILSIM program is the ability to display multiple post-run plot pages to allow a quick assessment of the time history response. There are 16 plot pages currently available to the user. Each plot page displays 9 parameters. Each parameter can also be displayed individually. on a one plot-per-page format. For a more refined display of the results the program can also create files of tabulated data. which can then be used by other plotting programs. The TAILSIM program was written straightforwardly assuming the user would want to change the database tables, the buildup equations, the output parameters. and the pilot model parameters. A separate database file and input file are automatically read in by the program. The use of an include file to set up all common blocks facilitates easy changing of parameter names and array sizes.

Hiltner, Dale W.↗

A theory of nonlocal mixing-length convection. I - The moment formalism

A flexible and potentially powerful theory of convection, based on the mixing length picture, is developed to make unbiased self-consistent predictions about overshooting and other complicated phenomena in convection. The basic formalism is set up, and the method's power is demonstrated by showing that a simplified version of the theory reproduces all the standard results of local convection. The second-order equations of the theory are considered in the limit of a steady state and vanishing third moments, and it is shown that they reproduce all the standard results of local mixing-length convection. There is a particular value of the superadiabatic gradient, below which the only possible steady state of a fluid is nonconvecting. Above this critical value, a fluid is convectively unstable. Two distinct regimes of convection, which are identified as efficient and inefficient convection, are determined.

Grossman, Scott A.↗

Multivariate space - time analysis of PRE-STORM precipitation

This paper presents the methodologies and results of the multivariate modeling and two-dimensional spectral and correlation analysis of PRE-STORM rainfall gauge data. Estimated parameters of the models for the specific spatial averages clearly indicate the eastward and southeastward wave propagation of rainfall fluctuations. A relationship between the coefficients of the diffusion equation and the parameters of the stochastic model of rainfall fluctuations is derived that leads directly to the exclusive use of rainfall data to estimate advection speed (about 12 m/s) as well as other coefficients of the diffusion equation of the corresponding fields. The statistical methodology developed here can be used for confirmation of physical models by comparison of the corresponding second-moment statistics of the observed and simulated data, for generating multiple samples of any size, for solving the inverse problem of the hydrodynamic equations, and for application in some other areas of meteorological and climatological data analysis and modeling.

Polyak, Ilya↗