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 415 records · Page 23

Strength evaluation of socket joints

This report documents the development of a set of equations that can be used to provide a relatively simple solution for identifying the strength of socket joints and for most cases avoid the need of more lengthy analyses. The analytical approach was verified by comparison of the contact load distributions to results obtained from a finite element analysis. The contacting surfaces for the specific joint in this analysis are in the shape of frustrums of a cone and are representative of the tapered surfaces in the socket-type joints used to join segments of model support systems for wind tunnels. The results are in the form of equations that can be used to determine the contact loads and stresses in the joint from the given geometry and externally applied loads. Equations were determined to define the bending moments and stresses along the length of the joints based on strength and materials principles. The results have also been programmed for a personal computer and a copy of the program is included.

Rash, Larry C.↗

Active load control during rolling maneuvers

A rolling maneuver load alleviation (RMLA) system has been demonstrated on the active flexible wing (AFW) wind tunnel model in the Langley Transonic Dynamics Tunnel (TDT). The objective was to develop a systematic approach for designing active control laws to alleviate wing loads during rolling maneuvers. Two RMLA control laws were developed that utilized outboard control-surface pairs (leading and trailing edge) to counteract the loads and that used inboard trailing-edge control-surface pairs to maintain roll performance. Rolling maneuver load tests were performed in the TDT at several dynamic pressures that included two below and one 11 percent above open-loop flutter dynamic pressure. The RMLA system was operated simultaneously with an active flutter suppression system above open-loop flutter dynamic pressure. At all dynamic pressures for which baseline results were obtained, torsion-moment loads were reduced for both RMLA control laws. Results for bending-moment load reductions were mixed; however, design equations developed in this study provided conservative estimates of load reduction in all cases.

Woods-Vedeler, Jessica A.↗

Scattering from perfectly conducting and resistive strips on a grounded dielectric slab

The scattering properties of perfectly conducting and resistive strips are predicted for strips which are located on a dielectric slab backed by a perfectly conducting ground plane. The spectral domain Green's function is used to relate the currents and fields on the strip, and the resulting integral equation is solved using the method of moments. Both TE and TM strips are examined using piecewise linear and pulse subdomain basis functions, respectively, to model the current on the strip. Calculated results are compared with results measured at the NASA Langley Research Center.

PULSE SUBDOMAIN BASIS FUNCTION↗

Strain-Gage Loads Calibration Parametric Study

This paper documents a parametric study of various aircraft wing-load test features that affect the quality of the resultant derived shear, bending-moment, and torque strain-gage load equations. The effect of the following on derived strain-gage equation accuracy are compared: single-point loading compared with distributed loading, variation in applied test load magnitude, number of applied load cases, and wing-box-only compared with control-surface loading. The subject of this study is an extensive wing-load calibration test of the Active Aeroelastic Wing F/A-18 airplane. Selected subsets of the available test data were used to derive load equations using the linear regression method. Results show the benefit of distributed loading and the diminishing-return benefits of test load magnitudes and number of load cases. The use of independent check cases as a quality metric for the derived load equations is shown to overcome blind extrapolating beyond the load data used to derive the load equations.

Lokos, William A.↗

Excitation of a Parallel Plate Waveguide by an Array of Rectangular Waveguides

This work addresses the problem of excitation of a parallel plate waveguide by an array of rectangular waveguides that arises in applications such as the continuous transverse stub (CTS) antenna and dual-polarized parabolic cylindrical reflector antennas excited by a scanning line source. In order to design the junction region between the parallel plate waveguide and the linear array of rectangular waveguides, waveguide sizes have to be chosen so that the input match is adequate for the range of scan angles for both polarizations. Electromagnetic wave scattered by the junction of a parallel plate waveguide by an array of rectangular waveguides is analyzed by formulating coupled integral equations for the aperture electric field at the junction. The integral equations are solved by the method of moments. In order to make the computational process efficient and accurate, the method of weighted averaging was used to evaluate rapidly oscillating integrals encountered in the moment matrix. In addition, the real axis spectral integral is evaluated in a deformed contour for speed and accuracy. The MoM results for a large finite array have been validated by comparing its reflection coefficients with corresponding results for an infinite array generated by the commercial finite element code, HFSS. Once the aperture electric field is determined by MoM, the input reflection coefficients at each waveguide port, and coupling for each polarization over the range of useful scan angles, are easily obtained. Results for the input impedance and coupling characteristics for both the vertical and horizontal polarizations are presented over a range of scan angles. It is shown that the scan range is limited to about 35 for both polarizations and therefore the optimum waveguide is a square of size equal to about 0.62 free space wavelength.

Rengarajan, Sembiam↗

Parallelization of a Six Degree of Freedom Entry Vehicle Trajectory Simulation Using OpenMP and OpenACC

The art and science of writing parallelized software, using methods such as Open Multi-Processing (OpenMP) and Open Accelerators (OpenACC), is dominated by computer scientists. Engineers and non-computer scientists looking to apply these techniques to their project applications face a steep learning curve, especially when looking to adapt their original single threaded software to run multi-threaded on graphics processing units (GPUs). There are significant changes in mindset that must occur; such as how to manage memory, the organization of instructions, and the use of if statements (also known as branching). The purpose of this work is twofold: 1) to demonstrate the applicability of parallelized coding methodologies, OpenMP and OpenACC, to tasks outside of the typical large scale matrix mathematics; and 2) to discuss, from an engineer’s perspective, the lessons learned from parallelizing software using these computer science techniques. This work applies OpenMP, on both multi-core central processing units (CPUs) and Intel® Xeon Phi™ 7210, and OpenACC on GPUs. These parallelization techniques are used to tackle the simulation of thousands of entry vehicle trajectories through the integration of six degree of freedom (DoF) equations of motion (EoM). The forces and moments acting on the entry vehicle, and used by the EoM, are estimated using multiple models of varying levels of complexity. Several benchmark comparisons are made on the execution of six DoF trajectory simulation: single thread Intel® Xeon® E5-2670 CPU, multi-thread CPU using OpenMP, multi-thread Xeon Phi™ 7210 using OpenMP, and multi-thread NVIDIA® Tesla® K40 GPU using OpenACC. These benchmarks are run on the Pleiades Supercomputer Cluster at the National Aeronautics and Space Administration (NASA) Ames Research Center (ARC), and a Xeon Phi™ 7210 node at NASA Langley Research Center (LaRC).

Green, Justin S.↗

Lectures on statistical mechanics

Presented here is a transcription of the lecture notes from Professor Allan N. Kaufman’s graduate statistical mechanics course Physics 212A and 212B at the University of California Berkeley from the 1972–1973 academic year. 212A addressed equilibrium statistical mechanics with topics: fundamentals (micro-canonical and sub-canonical ensembles, adiabatic law and action conservation, fluctuations, pressure, and virial theorem), classical fluids and other systems (equation of state, deviations from ideality, virial coefficients and van der Waals potential, canonical ensemble and partition function, quasistatic evolution, grand-canonical ensemble and partition function, chemical potential, simple model of a phase transition, quantum virial expansion, numerical simulation of equations of state, and phase transition), chemical equilibrium (systems with multiple species and chemical reactions, law of mass action, Saha equation, chemical equilibrium including ionization and excited states), and long-range interactions (including Coulomb, dipole, and gravitational interactions, Debye–Hückel theory, and shielding). 212B addressed nonequilibrium statistical mechanics with topics: fundamentals (definitions: realizations, moments, characteristic function, and discrete variables), Brownian motion (Langevin equation, fluctuation–dissipation theorem, spatial diffusion, Boltzmann’s H-theorem), Liouville and Klimontovich equations, Landau equation (derivation, elaboration, and H-theorem, and irreversibility), Markov processes and Fokker–Planck equation (derivations of the Fokker–Planck equation and a master equation), linear response and transport theory (linear Boltzmann equation, linear response theory of Kubo and Mori, relation of entropy production to electrical conductivity, transport relations and coefficients, normal mode solutions of the transport equations, sketch of a generalized Langevin equation method for transport theory), and an introduction to nonequilibrium quantum statistical mechanics.

plasma dynamics↗

Nonlinear aerodynamic modeling of flap oscillations in transonic flow - A numerical validation

The regime of validity of a nonlinear aerodynamic force and moment formulation, based on concepts from nonlinear functional analysis and applicable to a transonic airfoil with a deflecting flap, is investigated. A time-dependent finite difference technique is used to evaluate the aerodynamic data of the formulation in terms of specified, characteristic motions. Flap-motion histories are generated from the flap inertial equations of motion, with aerodynamic reactions specified by the moment formulation. The motion histories depicting the cases of decaying and growing flap oscillations are compared with histories generated through simultaneous, coupled solution of the fluid-dynamic equations and flap inertial equations of motion. The range of applicability of the formulation is discussed.

Chyu, W. J.↗

Hidden Variables and Magnetic Relaxation

Nonlinear generalization of Schroedinger equation describing hidden variables effect solved for spin-1/2 magnetic moment relaxing in constant magnetic field

MAGNETIC RELAXATION↗

Probability Density Functions of Observed Rainfall in Montana

The question of whether a rain rate probability density function (PDF) can vary uniformly between precipitation events is examined. Image analysis on large samples of radar echoes is possible because of advances in technology. The data provided by such an analysis easily allow development of radar reflectivity factors (and by extension rain rate) distribution. Finding a PDF becomes a matter of finding a function that describes the curve approximating the resulting distributions. Ideally, one PDF would exist for all cases; or many PDF's that have the same functional form with only systematic variations in parameters (such as size or shape) exist. Satisfying either of theses cases will, validate the theoretical basis of the Area Time Integral (ATI). Using the method of moments and Elderton's curve selection criteria, the Pearson Type 1 equation was identified as a potential fit for 89 percent of the observed distributions. Further analysis indicates that the Type 1 curve does approximate the shape of the distributions but quantitatively does not produce a great fit. Using the method of moments and Elderton's curve selection criteria, the Pearson Type 1 equation was identified as a potential fit for 89% of the observed distributions. Further analysis indicates that the Type 1 curve does approximate the shape of the distributions but quantitatively does not produce a great fit.

Larsen, Scott D.↗

Some Remarks on GMRES for Transport Theory

We review some work on the application of GMRES to the solution of the discrete ordinates transport equation in one-dimension. We note that GMRES can be applied directly to the angular flux vector, or it can be applied to only a vector of flux moments as needed to compute the scattering operator of the transport equation. In the former case we illustrate both the delights and defects of ILU right-preconditioners for problems with anisotropic scatter and for problems with upscatter. When working with flux moments we note that GMRES can be used as an accelerator for any existing transport code whose solver is based on a stationary fixed-point iteration, including transport sweeps and DSA transport sweeps. We also provide some numerical illustrations of this idea. We finally show how space can be traded for speed by taking multiple transport sweeps per GMRES iteration. Key Words: transport equation, GMRES, Krylov subspace

Patton, Bruce W.↗

Scalar flux transport models for self-similar turbulent mixing

A common approach to closing turbulent species flux in multicomponent Reynolds-averaged Navier-Stokes models is to use the standard gradient diffusion approximation. While such an approach has been shown to work well when applied to many canonical turbulent mixing configurations, a gradient diffusion approach is fundamentally limited in its ability to capture complex phenomena such as countergradient transport. For this reason, complicated mixing applications may benefit by treating the turbulent diffusivity with a model transport equation in a manner analogous to second-moment momentum closure in Reynolds-stress transport models. Here, the present work explores the development and application of two different scalar flux transport (SFT) models. Self-similarity constraints are derived for these models, and they are evaluated against gradient-diffusion-based models in several one- and two-dimensional problems of turbulent mixing. It is found that the new SFT models out-perform gradient diffusion models in problems involving rapid acceleration reversal and in problems involving anisotropic transport of materials. In addition, it is found that even a hybrid-SFT approach, in which an SFT equation is utilized along with a gradient diffusion closure, provides some measure of improvement over models that transport the mass flux rather than the scalar flux.

Reynolds-averaged Navier Stokes↗

Pitch angle diffusion in the Jovian magnetodisc

Using a linear approximation to the field variation and a mathematical technique for distorting the orbit integration into the complex gyrophase plane, a formula is derived for the nonadiabatic change in magnetic moment. This change varies as the cosine of the particle gyrophase at the equator, so that if a particle traverses the equator several times with different values of this gyrophase, it suffers randomly correlated values for the nonadiabatic change in magnetic moment. The algebraic formula for this nonadiabatic change is found to be in excellent agreement with previous numerical computations. The phase space density, averaged over several bounce periods, is found to satisfy a diffusion equation in pitch angle. The diffusion coefficient is dependent on particle energy, pitch angle, and the field line along which the particles are moving.

Birmingham, T. J.↗

Second Order Closures for the Radiative Transfer Equation: Some Are Unstable

The largest existing simulations of cosmic reionization model radiative transfer with moment methods that require a closure relation. The two most commonly used closure relations are M1 and OTVET; both close the moment hierarchy at the first moment. We explore the properties of a higher, second-order closure. We show that direct generalizations of M1 and OTVET to one higher order are physically unstable - i.e., the closure equations themselves result in unstable solutions, not just their numerical implementation. In fact, a generalization of OTVET to any order higher than the first one is unstable. We are also able to show that any local (i.e., depending only on the local moments of the radiation field, like M1) second-order closure that depends only on the radiation intensity and radiation flux, but does not explicitly depend on the radiation pressure, is physically unstable. This result restricts the choice of possible second-order closure relations.

Gnedin, Nickolay Y. [Fermilab; Chicago U., Astron.↗

Stochastic Transport Theory for Investigating the Three-Dimensional Canopy Structure from Space Measurements

Radiation reflected from vegetation canopies exhibits high spatial variation. Satellite-borne sensors measure the mean intensities emanating from heterogeneous vegetated pixels. The theory of radiative transfer in stochastic media provides the most logical linkage between satellite observations and the three-dimensional canopy structure through a closed system of simple equations which contains the mean intensity and higher statistical moments directly as its unknowns. Although this theory has been a highly active research field in recent years, its potential for satellite remote sensing of vegetated surfaces has not been fully realized because of the lack of models of a canopy pair-correlation function that the stochastic radiative transfer equations require. The pair correlation function is defined as the probability of finding simultaneously phytoelements at two points. This paper presents analytical and Monte Carlo generated pair correlation functions. Theoretical and numerical analyses show that the spatial correlation between phytoelements is primarily responsible for the effects of the three-dimensional canopy structure on canopy reflective and absorptive properties. The pair correlation function, therefore, is the most natural and physically meaningful measure of the canopy structure over a wide range of scales. The stochastic radiative transfer equations naturally admit this measure and thus provide a powerful means to investigate the three-dimensional canopy structure from space. Canopy reflectances predicted by the stochastic equations are assessed by comparisons with the PARABOLA measurements from coniferous and broadleaf forest stands in the BOREAS Southern Study Areas. The pair correlation functions are derived from data on tree structural parameters collected during field campaigns conducted at these sites. The simulated canopy reflectances compare well with the PARABOLA data.

Huang, Dong↗

Inertial and anisotropic pressure effects on cross-field electron transport in low-temperature magnetized plasmas

In this paper, a one-dimensional (1D) particle-in-cell Monte Carlo collision (PIC-MCC) model is developed to investigate the effects of anisotropic pressure and inertial terms due to non-Maxwellian velocity distribution functions on cross-field electron transport. The conservation of momentum is evaluated by taking the moments of the first-principles gas-kinetic equation. A steady-state discharge is obtained without any low-frequency ionization oscillations by considering an anomalous electron scattering profile. Further, the results obtained from the 1D PIC-MCC model are compared with fluid models, including the quasi-neutral drift-diffusion (DD), non-neutral DD, and full fluid moment models. The discharge current obtained from the PIC-MCC model is in good agreement with the fluid models. The cross-field electron transport due to the inertial terms, i.e. the gradient of axial and azimuthal drift, is evaluated. Moreover, PIC-MCC simulation results show non-zero, anisotropic, off-diagonal pressure tensor terms due to asymmetric non-Maxwellian electron velocity distribution function, potentially contributing to cross-field electron transport.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

General equations for the stress analysis of rings

In this report it is shown that the shear, axial force, and moment at one point in a simple ring subjected to any loading condition can be given by three independent equations involving certain integrals that must be evaluated regardless of the method of analysis used. It is also shown how symmetry of the ring alone or of the ring and the loading about 1 or 2 axes makes it possible to simplify the three equations and greatly reduces the number of integrals that must be evaluated. Application of the general equations presented in this report to practical problems in the stress analysis of rings makes it possible to shorten, simplify, and systematize the calculations for both simple and braced rings. Three illustrative problems are included to demonstrate the application of the general equations to a simple ring with different loadings.

Lundquist, Eugene E↗