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 505 records · Page 28

Sensitivity Analysis for Coupled Aero-structural Systems

A novel method has been developed for calculating gradients of aerodynamic force and moment coefficients for an aeroelastic aircraft model. This method uses the Global Sensitivity Equations (GSE) to account for the aero-structural coupling, and a reduced-order modal analysis approach to condense the coupling bandwidth between the aerodynamic and structural models. Parallel computing is applied to reduce the computational expense of the numerous high fidelity aerodynamic analyses needed for the coupled aero-structural system. Good agreement is obtained between aerodynamic force and moment gradients computed with the GSE/modal analysis approach and the same quantities computed using brute-force, computationally expensive, finite difference approximations. A comparison between the computational expense of the GSE/modal analysis method and a pure finite difference approach is presented. These results show that the GSE/modal analysis approach is the more computationally efficient technique if sensitivity analysis is to be performed for two or more aircraft design parameters.

Giunta, Anthony A.↗

A parallel-kinetic-perpendicular-moment model for magnetised plasmas

We describe a new model for the study of weakly collisional, magnetised plasmas derived from exploiting the separation of the dynamics parallel and perpendicular to the magnetic field. This unique system of equations retains the particle dynamics parallel to the magnetic field while approximating the perpendicular dynamics through a spectral expansion in the perpendicular degrees of freedom, analogous to moment-based fluid approaches. In so doing, a hybrid approach is obtained that is computationally efficient enough to allow for larger-scale modelling of plasma systems while eliminating a source of difficulty in deriving fluid equations applicable to magnetised plasmas. We connect this system of equations to historical asymptotic models and discuss advantages and disadvantages of this approach, including the extension of this parallel-kinetic-perpendicular moment beyond the typical region of validity of these more traditional asymptotic models. This paper forms the first of a multi-part series on this new model, covering the theory and derivation, alongside demonstration benchmarks of this approach that include shocks and magnetic reconnection.

astrophysical plasmas↗

Prediction of unsteady separated flows on oscillating airfoils

Techniques for calculating high Reynolds number flow around an airfoil undergoing dynamic stall are reviewed. Emphasis is placed on predicting the values of lift, drag, and pitching moments. Methods discussed include: the discrete potential vortex method; thin boundary layer method; strong interaction between inviscid and viscous flows; and solutions to the Navier-Stokes equations. Empirical methods for estimating unsteady airloads on oscillating airfoils are also described. These methods correlate force and moment data from wind tunnel tests to indicate the effects of various parameters, such as airfoil shape, Mach number, amplitude and frequency of sinosoidal oscillations, mean angle, and type of motion.

Mccroskey, W. J.↗

Strain Gage Load Calibration of the Wing Interface Fittings for the Adaptive Compliant Trailing Edge Flap Flight Test

The safety-of-flight parameters for the Adaptive Compliant Trailing Edge (ACTE) flap experiment require that flap-to-wing interface loads be sensed and monitored in real time to ensure that the structural load limits of the wing are not exceeded. This paper discusses the strain gage load calibration testing and load equation derivation methodology for the ACTE interface fittings. Both the left and right wing flap interfaces were monitored; each contained four uniquely designed and instrumented flap interface fittings. The interface hardware design and instrumentation layout are discussed. Twenty-one applied test load cases were developed using the predicted in-flight loads. Pre-test predictions of strain gage responses were produced using finite element method models of the interface fittings. Predicted and measured test strains are presented. A load testing rig and three hydraulic jacks were used to apply combinations of shear, bending, and axial loads to the interface fittings. Hardware deflections under load were measured using photogrammetry and transducers. Due to deflections in the interface fitting hardware and test rig, finite element model techniques were used to calculate the reaction loads throughout the applied load range, taking into account the elastically-deformed geometry. The primary load equations were selected based on multiple calibration metrics. An independent set of validation cases was used to validate each derived equation. The 2-sigma residual errors for the shear loads were less than eight percent of the full-scale calibration load; the 2-sigma residual errors for the bending moment loads were less than three percent of the full-scale calibration load. The derived load equations for shear, bending, and axial loads are presented, with the calculated errors for both the calibration cases and the independent validation load cases.

Miller, Eric J.↗

Higher-Order Statistical Moments of Predicted Sonic Boom Waveforms Through Turbulence

Sonic booms generated by supersonic aircraft are affected by turbulence in the atmospheric boundary layer through which they propagate. Turbulence effects lead to random variability of the sonic boom waveforms measured on the ground, complicating the prediction of such waveforms. As an initial effort to predict the waveform variability, the solution of the coherent or mean sonic boom waveform has previously been formulated by the author. The current paper extends the formulation to derive an expression of the second-order statistical moment necessary to calculate the variance of the spectral amplitudes. Since the derivation uses a full wave equation, results using the derived expression will be compared with those obtained using a parabolic approximation. In order to fully quantify the uncertainty of our predictions, the higher order statistical moments are also formulated and are shown to be approximately zero. Consequently, the probability density function (pdf) of the spectral amplitudes is predicted to be Gaussian. The formulation is further extended to determine the pdf of the loudness of sonic booms and quantify the uncertainties associated with the loudness prediction. Results from the formulation are compared with available flight test data.

sonic booms↗

Model for rotating and nonuniform planetary exospheres.

The model neutral exosphere with a uniformly rotating exobase is generalized by allowing variations in exobase density and temperature which characterize the thermosphere just below the base. The corresponding velocity distribution function, satisfying the collisionless Boltzmann equation, is constructed and used to form a general expression for the velocity moments. Resulting density profiles of rotating exospheres with nonuniform densities and temperatures on the exobase are compared with corresponding nonrotating exospheres. Density enhancements due to rotation are found to be greatest above regions of exobase density or temperature minima. Equatorial density enhancements of terrestrial hydrogen, resulting from rotation, are estimated to be 15%-17% at altitudes of 10-20 earth radii. Corresponding increases in terrestrial helium are 30%-50% on the equator at altitudes of 0.7 to 1 earth radii even when there is a polar density bulge in the barosphere.

Hartle, R. E.↗

Moment-based adaptive time integration for thermal radiation transport

Here, in this paper we develop a framework for moment-based adaptive time integration of deterministic multifrequency thermal radiation transpot (TRT). We generalize our recent semi-implicit-explicit (IMEX) integration framework for gray TRT to multifrequency TRT, and also introduce a semi-implicit variation that facilitates higher-order integration of TRT, where each stage is implicit in all components except opacities. To appeal to the broad literature on adaptivity with Runge–Kutta methods, we derive new embedded methods for four asymptotic preserving IMEX Runge–Kutta schemes we have found to be robust in our previous work on TRT and radiation hydrodynamics. We then use a moment-based high-order-low-order representation of the transport equations. Due to the high dimensionality, memory is always a concern in simulating TRT. We form error estimates and adaptivity in time purely based on temperature and radiation energy, for a trivial overhead in computational cost and memory usage compared with the base second order integrators. We then test the adaptivity in time on the tophat and Larsen problem, demonstrating the ability of the adaptive algorithm to naturally vary the timestep across 4–5 orders of magnitude, ranging from the dynamical timescales of the streaming regime to the thick diffusion limit.

97 MATHEMATICS AND COMPUTING↗

Bhatnagar-Gross-Krook subgrid model for neutrino quantum kinetics

Here, we present a new subgrid model for neutrino quantum kinetics, which is primarily designed to incorporate effects of collective neutrino oscillations into neutrino-radiation-hydrodynamic simulations for core-collapse supernovae and mergers of compact objects. We approximate the neutrino oscillation term in a quantum kinetic equation using the Bhatnagar-Gross-Krook (BGK) relaxation-time prescription, and the transport equation is directly applicable for classical neutrino transport schemes. The BGK model is motivated by recent theoretical indications that nonlinear phases of collective neutrino oscillations settle into quasisteady structures. We explicitly provide basic equations of the BGK subgrid model for both multiangle and moment-based neutrino transport to facilitate the implementation of the subgrid model in the existing neutrino transport schemes. We also show the capability of our BGK subgrid model by comparing it to fully quantum kinetic simulations for fast neutrino-flavor conversion. We find that the overall properties can be well reproduced in the subgrid model; the error of angular-averaged survival probability of neutrinos is within ~20%. By identifying the source of error, we also discuss perspectives to improve the accuracy of the subgrid model.

79 ASTRONOMY AND ASTROPHYSICS↗

Effects of Different Closure Choices in Core-collapse Supernova Simulations

The two-moment method is widely used to approximate the full neutrino transport equation in core-collapse supernova (CCSN) simulations, and different closures lead to subtle differences in the simulation results. In this paper, we compare the effects of closure choices on various physical quantities in 1D and 2D time-dependent CCSN simulations with our multigroup radiation hydrodynamics code FORNAX. We find that choices of the third-order closure relations influence the time-dependent simulations only slightly. Choices of the second-order closure relation have larger consequences than choices of the third-order closure, but these are still small compared to the remaining variations due to ambiguities in some physical inputs such as the nuclear equation of state. We also find that deviations in Eddington factors are not monotonically related to deviations in physical quantities, which means that simply comparing the Eddington factors does not inform one concerning which closure is better.

79 ASTRONOMY AND ASTROPHYSICS↗

Investigation of mass transfer between two parallel walls at different temperatures by a moment method

One-dimensional flow between two fixed parallel walls composed of the same substance but at different temperatures and spaced a distance 1 apart is considered. The hot plate is the evaporating surface (source) and the cold plate is the condensing surface (sink). The vapor between the two plates is assumed to be a monatomic gas consisting of Maxwell molecules. Lee's moment method is used to obtain a set of six nonlinear equations. Both the nonlinear equations and a linearized approximation to them are solved.

Sloat, T. N.↗

A simulation of saltation dynamics with application to terrestrial and Martian model surfaces.

The feasibility of predicting the saltation threshold from a dynamic model of the saltation process is demonstrated. The model is based on mechanical equations that describe the trajectories of particles (sand grains) near the ground-atmosphere interface. The equations of particle motion are formulated in terms of forces and moments. The forces employed are based on theoretical and experimental results obtained for continuum flow. Surface structuring analytics and simulated fluid turbulence further modify the particle behavior. The model is used to predict the threshold velocities on Mars.

Andres, R. M.↗

Continuous and discrete describing function analysis of the LST system

A describing function of the control moment gyros (CMG) frictional nonlinearity is derived using the analytic torque equation. Computer simulation of the simplified Large Space Telescope (LST) system with the analytic torque expression is discussed along with the transfer functions of the sampled-data LST system, and the discrete describing function of the GMC frictionality.

Kuo, B. C.↗

Analytical formulation of selected activities of the remote manipulator system

Existing analysis of Orbiter-RMS-Payload kinematics were surveyed, including equations dealing with the two body kinematics in the presence of a massless RMS and compares analytical explicit solutions with numerical solutions. For the following operational phases of the RMS numerical demonstration, problems are provided: (1) payload capture; (2) payload stowage and removal from cargo bay; and (3) payload deployment. The equation of motion provided accounted for RMS control forces and torque moments and could be extended to RMS flexibility and control loop simulation without increasing the degrees of freedom of the two body system.

K J Zimmerman↗

Effects of bearing offset and flexibility on the mesh force distribution of spiral bevel gears

For straight or spiral bevel gears of nominal spherical involute design, the resultant total force vector transmitted by the gear mesh, in the absence of friction, lies in the plane of tooth contact. This force vector can be characterized by three scalar components, two orthogonal force components lying in the plane of contact and the resultant moment taken about the nominal center of the zone of contact. Equations for these three generalized force components are derived. The equations are expressed in terms of tooth pair/gear body stiffnesses, bearing/bearing support flexibility influence coefficients, the shaft input torque, deviations of the tooth running surfaces from perfect spherical involute surfaces, and bearing centerline offsets from the positions occupied by the base cone axes of the perfect involute bevel gear counterparts to the actual gears under consideration. Inertial forces arising from transverse and axial vibrations of the gear bodies are assumed to be negligible in comparison with the bearing support reaction forces.

Mark, W. D.↗

Backscatter RCS for TE and TM excitations of dielectric-filled cavity-backed apertures in two-dimensional bodies

Transverse electric (TE) and transverse magnetic (TM) scattering from dielectric-filled, cavity-backed apertures in two-dimensional bodies are treated using the method of moments technique to solve a set of combined-field integral equations for the equivalent induced electric and magnetic currents on the exterior of the scattering body and on the associated aperture. Results are presented for the backscatter radar cross section (RCS) versus the electrical size of the scatterer for two different dielectric-filled cavity-backed geometries. The first geometry is a circular cylinder of infinite length which has an infinite length slot aperture along one side. The cavity inside the cylinder is dielectric filled and is also of circular cross section. The two cylinders (external and internal) are of different radii and their respective longitudinal axes are parallel but not collocated. The second is a square cylinder of infinite length which has an infinite length slot aperture along one side. The cavity inside the square cylinder is dielectric-filled and is also of square cross section.

Goggans, Paul M.↗

Navier-Stokes prediction of large-amplitude forced and free-to-roll delta-wing oscillations

The three-dimensional Reynolds-averaged, Navier-Stokes (RANS) equations are used to numerically simulate nonsteady vortical flow about a 65 degree sweep delta wing at 30 degrees angle of attack. Two large-amplitude, high-rate, forced-roll motions and a damped free-to-roll motion are presented. The free-to-roll motion is computed by coupling the time-dependent RANS equations to the flight dynamic equation of motion. The computed results are compared with experimental forces, moments, and roll-angle time histories. The overall agreement is good. Vortex breakdown is present in each case, which causes significant time lags in the vortex breakdown motions relative to the body motions. This behavior strongly influences the dynamic forces and moments.

Chaderjian, Neal M.↗

Theoretical Elastic Stress Distributions Arising from Discontinuities and Edge Loads in Several Shell-Type Structures

The deformation and complete stress distribution are determined for each of the following edge loaded thin shells of revolution: (1) a right circular cylinder, (2) a frustum of a right circular cone, and (3) a portion of a sphere. The locations of the maximum circumferential and meridional stresses on both the inner and outer surfaces are also found. The basic equations for the above were selected from the published literature on the subject and expanded to produce to resultant-stress equations in closed from where practicable to do so. Equations are also developed for the discontinuity shear force and bending moment at each of the following junction: (1) axial change of thickness in a circular cylinder, (2) axial change of thickness in a cone, (3) change of thickness in a portion of a sphere, (4) a cylinder and a cone, (5) a cylinder and a portion of a sphere(6) a cylinder and a flat head, and (7) a cone and a portion of a sphere.

Johns, Robert H.↗