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 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

An immersed interface method for the 2D vorticity-velocity Navier-Stokes equations with multiple bodies

We present an immersed interface method for the vorticity-velocity form of the 2D Navier Stokes equations that directly addresses challenges posed by nonconvex immersed bodies, multiply connected domains, and the calculation of force distributions on immersed surfaces. The immersed interface method is re-interpreted as a polynomial extrapolation of flow quantities and boundary conditions into the immersed solid bodies, reducing computational cost and enabling simulations with nonconvex bodies that could not be discretized with previous immersed interface methods. In the flow, the vorticity transport equation is discretized using a conservative finite difference scheme and explicit Runge-Kutta time integration. The velocity reconstruction problem is transformed to a scalar Poisson equation that is discretized with conservative finite differences, and solved using an FFT-accelerated iterative algorithm. The use of conservative differencing throughout leads to exact enforcement of a discrete Kelvin's theorem, allowing for simulations with multiply connected domains and outflow boundaries that have challenged other immersed interface vortex methods. We also explore novel methods for recovering time-dependent pressure distributions on immersed bodies within a vorticity-based method and present a novel control volume formulation for recovering aerodynamic moments from only the vorticity and velocity fields. The method achieves second order spatial accuracy and third order temporal accuracy, and is validated on a variety of 2D flows in internal and free-space domains.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1979

Gust alleviation using direct turbulence measurements

The research reported upon in this paper describes an effective method of gust alleviation using direct measurements of atmospheric turbulence to drive the aircraft control surfaces in a way that attempts to directly counter or cancel those forces and moments produced on the aircraft by gusts. The method yields a feedforward or open loop control law, simple to mechanize and relatively insensitive to changes in flight condition. When applied directly, the resulting control law effectively gust-alleviates in the low frequency phugoid and short period range but has a tendency to amplify structural mode vehicle motions due to the phase lag of the actuators. A method of design based upon the use of a diagonal or Jordan form of the equations of motion enables the designer to avoid this problem of structural mode excitation.

Rynaski, E. G.↗
NASA NTRS1989

Connection forces in deformable multibody dynamics

In the dynamic formulation of holonomic and nonholonomic systems based on D'Alembert-Lagrange equation, the forces of constraints are maintained in the dynamic equations by introducing auxiliary variables, called Lagrange multipliers. This approach introduces a set of generalized reaction forces associated with the system generalized coordinates. Different sets of variables can be used as generalized coordinates and accordingly, the generalized reactions associated with these generalized coordinates may not be the actual reaction forces at the joints. In rigid body dynamics, the generalized reaction forces and the actual reaction forces at the joints represent equipollent systems of forces since they produce the same total forces and moments at and about any point on the rigid body. This is not, however, the case in deformable body analyses wherein the generalized reaction forces depend on the system generalized reference and elastic coordinates. In this paper, a method for determining the actual reaction forces at the joints from the generalized reaction forces in deformable multibody systems is presented.

Shabana, A. A.↗
NASA NTRS1996

Second moment closure analysis of the backstep flow database

A Second Moment Closure computation (SMC) is compared in detail with the Direct Numerical Simulation (DNS) data of Le and Moin for the backstep flow at Re = 5,000 in an attempt to understand why the intensity of the backflow and, consequently, the friction coefficient in the recirculation bubble are severely underestimated. The data show that this recirculation bubble is far from being laminar except in the very near wall layer. A novel 'differential a priori' procedure was used, in which the full transport equation for one isolated component of the Reynolds stress tensor was solved using DNS data as input. Conclusions are then different from what would have been deduced by comparing a full simulation to a DNS. One cause of discrepancy was traced back to insufficient transfer of energy to the normal stress by pressure strain, but was not cured. A significant finding, confirmed by the DNS data in the core region of a channel flow, is that the coefficient that controls destruction of dissipation, C epsilon(sub 2), should be decreased by a factor of 2 when production is vanishing. This is also the case in the recirculation bubble, and a new formulation has cured 25% of the backflow discrepancy.

Parneix, S.↗
NASA NTRS2009

Exploring Alternative Parameterizations for Snowfall with Validation from Satellite and Terrestrial Radars

Increases in computational resources have allowed operational forecast centers to pursue experimental, high resolution simulations that resolve the microphysical characteristics of clouds and precipitation. These experiments are motivated by a desire to improve the representation of weather and climate, but will also benefit current and future satellite campaigns, which often use forecast model output to guide the retrieval process. The combination of reliable cloud microphysics and radar reflectivity may constrain radiative transfer models used in satellite simulators during future missions, including EarthCARE and the NASA Global Precipitation Measurement. Aircraft, surface and radar data from the Canadian CloudSat/CALIPSO Validation Project are used to check the validity of size distribution and density characteristics for snowfall simulated by the NASA Goddard six-class, single moment bulk water microphysics scheme, currently available within the Weather Research and Forecast (WRF) Model. Widespread snowfall developed across the region on January 22, 2007, forced by the passing of a mid latitude cyclone, and was observed by the dual-polarimetric, C-band radar King City, Ontario, as well as the NASA 94 GHz CloudSat Cloud Profiling Radar. Combined, these data sets provide key metrics for validating model output: estimates of size distribution parameters fit to the inverse-exponential equations prescribed within the model, bulk density and crystal habit characteristics sampled by the aircraft, and representation of size characteristics as inferred by the radar reflectivity at C- and W-band. Specified constants for distribution intercept and density differ significantly from observations throughout much of the cloud depth. Alternate parameterizations are explored, using column-integrated values of vapor excess to avoid problems encountered with temperature-based parameterizations in an environment where inversions and isothermal layers are present. Simulation of CloudSat reflectivity is performed by adopting the discrete-dipole parameterizations and databases provided in literature, and demonstrate an improved capability in simulating radar reflectivity at W-band versus Mie scattering assumptions.

Molthan, Andrew L.↗
NASA NTRSDate not supplied

Developing Urban Air Mobility Vehicle Models to Support Air Traffic Management Concept Development

To support Urban Air Mobility (UAM) research efforts at NASA, the Airspace Target Generator (ATG) software used in the FutureFlight Central (FFC) air traffic control tower simulator is undergoing updates to support physics-based UAM vehicle models. A process was developed to integrate UAM vertical takeoff and landing (VTOL) aircraft into the fixed-wing ATG modeling environment without significant change to the underlying equations of motion and vehicle model database. The VTOL aircraft models were converted from a six degrees-of-freedom (6-DOF) representation into a four degrees-of-freedom (4-DOF) representation for integration within ATG. Three vehicle designs from the NASA Revolutionary Vertical-Lift Technologies (RVLT) project were selected: a lift-plus-cruise (LPC) aircraft model and quadrotor, electric-powered (QEP) 1-seater and 6-seater models. With the LPC model comprised of a nonlinear force and moment build-up, and the QEP models comprised of linearized stability derivatives, two separate processes were developed to convert the lift, drag, and propulsion characteristics of each model into the ATG model database. Key aircraft performance characteristics including climb, cruise, and descent performance were preserved during the conversion process. Because ATG simulates fixed-wing aircraft through ground taxi and takeoff to approach and landing, acceleration command algorithms were developed to model the vertical takeoff and vertical landing phase of UAM operations. A strategy was then developed to transition the aircraft model to- and from- the new control mode.

urban air mobility↗
DOE OSTI2023

A New Hybrid Mass-Flux/High-Order Turbulence Closure for Ocean Vertical Mixing

While various parameterizations of vertical turbulent fluxes at different levels of complexity have been proposed, each has its own limitations. For example, simple first-order closure schemes such as the K-Profile Parameterization (KPP) lack energetic constraints; two-equation models $k-ε$ like directly solve an equation for the turbulent kinetic energy but do not account for non-diffusive fluxes, and high-order closures that include the high-order transport terms are computationally expensive. To address these, we extend the Assumed-Distribution Higher-Order Closure (ADC) framework originally proposed for the atmospheric boundary layer and apply it to the ocean surface boundary layer. By assuming a probability distribution function relationship between the vertical velocity and tracers, all second-order and higher-order moments are exactly constructed and turbulence closure is achieved in the ADC scheme. In addition, this ADC parameterization has full energetic constraints and includes non-diffusive fluxes without the computational cost of a full higher-order closure scheme. We have tested the ADC scheme against a combination of large eddy simulation (LES), KPP, and $k-ε$ for surface buoyancy-driven convective mixing and found that the ADC scheme is robust with different vertical resolutions and compares well to the LES results.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1999

Recent Development in the Analytical Investigation of Shipboard Rotorcraft Engage and Disengage Operations

This paper presents an overview of recent developments in an effort to predict transient aeroelastic rotor response during shipboard engage and disengage sequences. The blade is modeled as an elastic beam undergoing in flap, lag, extension and torsion. The blade equations of motion are formulated using Hamilton's principle and they are spatially discretized using the finite element method. The discretized blade equations of motion are integrated for a specified rotor speed run-up or run-down profile. Blade element theory is used to calculate quasi-steady or unsteady aerodynamic loads in linear and nonlinear regimes. The analysis is capable of simulating both articulated, hingeless, and gimballed rotor systems. Validation of the rotor code is discussed, including correlation with droop stop impact tests and wind tunnel experiments. Predictions of safe engagement and disengagement envelopes, limited by excessive blade tip deflections or hub moments, are presented. Future directions of study are also discussed.

Smith, Edward C.↗
DOE OSTI2022

Magnetic field mapping of inaccessible regions using physics-informed neural networks

A difficult problem concerns the determination of magnetic field components within an experimentally inaccessible region when direct field measurements are not feasible. In this paper, we propose a new method of accessing magnetic field components using non-disruptive magnetic field measurements on a surface enclosing the experimental region. Magnetic field components in the experimental region are predicted by solving a set of partial differential equations (Ampere’s law and Gauss’ law for magnetism) numerically with the aid of physics-informed neural networks (PINNs). Prediction errors due to noisy magnetic field measurements and small number of magnetic field measurements are regularized by the physics information term in the loss function. We benchmark our model by comparing it with an older method. The new method we present will be of broad interest to experiments requiring precise determination of magnetic field components, such as searches for the neutron electric dipole moment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1997

Generalized Kinetic Description of Steady-State Collisionless Plasmas

We present a general solution to the collisionless Boltzmann (Vlasov) equation for a free-flowing plasma along a magnetic field line using Liouville's theorem, allowing for an arbitrary potential structure including non-monotonicities. The constraints of the existing collisionless kinetic transport models are explored, and the need for a more general approach to the problem of self- consistent potential energy calculations is described. Then a technique that handles an arbitrary potential energy distribution along the field line is presented and discussed. For precipitation of magnetospherically trapped hot plasma, this model yields moment calculations that vary by up to a factor of two for various potential energy structures with the same total potential drop. The differences are much greater for the high-latitude outflow scenario, giving order of magnitude variations depending on the shape of the potential energy distribution.

Khazanov, G. V.↗
NASA NTRS1995

Progress in Navier-Stokes Simulation of Delta Wing Maneuver Aerodynamics

This paper will summarize recent progress in the numerical simulation of high incidence vortical flow about a generic 65 degree sweep delta wing using the three dimensional, time-dependent, Reynolds averaged, Navier-Stokes (RANS) equations. Computations have been carried out at 15 and 30 degrees angle of attack under subsonic turbulent flow conditions, and compared with experimental data provided by Hanff, Jenkins, and their colleagues. This work has already been published elsewhere and widely disseminated. Computations carried out at 15 degrees angle of attack included static roll angles ranging up to 65 degrees, and a large-amplitude (40 degree), high rate (7 Hz), forced roll motion. There was very good agreement between computed and experimental forces and moments, and static surface pressures. There was a significant hysteresis of the dynamic rolling moment due to the high rate of roll motion. At this angle of attack, no vortex breakdown was observed in the computations or experiment. Computations were also carried out at 30 degrees angle of attack, where vortex breakdown was present in both the computations and experiment. There was overall good agreement in the computed and experimental forces and moments. The static rolling moment varied with roll angle in a highly nonlinear manner, and exhibited three stable trim points and two unstable trim points. This behavior was attributed to the presence of vortex breakdown. Two large-amplitude (30 degrees), high-rate (10 Hz) forced roll motions were computed. The dynamics of the vortex breakdown motion was dramatically visualized by tracking the time-dependent motion of particles released near the delta wing apex. This numerical visualization is analogous to experimental smoke flow techniques. In one of the dynamic cases the breakdown was found to move off the wing, convected downwind of the trailing edge, and later reformed near the trailing edge through an instability of the vortex core. A damped free-to-roll motion was also computed by releasing the wing from rest at 40 degrees of roll. The wing went to the same trim point as in the experiment.

Chaderjian, Neal M.↗
NASA NTRS1983

Estimation of parameters involved in high angle-of-attack aerodynamic theory using spin flight test data

The difficulty in applying parameter estimation techniques to spinning airplanes is due in part to the unwieldy number of possible combinations of terms in the equations of motion, when the model structure is unknown. The combination of high angle of attack and high rotation rate results in aerodynamic functions which are quite complex. For wing dominated configurations it is advantageous to use aerodynamic theory to generate the model structure. In this way, the number of unknown parameters is reduced and the model accuracy may be increased. Under conditions for which the theory is inadequate, however, model accuracy may be reduced. Strip theory, for example, is incapable of predicting autorotative rolling moments indicated by wind tunnel tests at angles of attack exceeding 40 degrees. An improved aerodynamic theory would be necessary to successfully apply the technique advanced for such regions.

Taylor, L. W., Jr.↗
NASA NTRS1993

Navier-Stokes prediction of a delta wing in roll with vortex breakdown

The three-dimensional, Reynolds-averaged, Navier-Stokes (RANS) equations are used to numerically simulate vortical flow about a 65 degree sweep delta wing. Subsonic turbulent flow computations are presented for this delta wing at 30 degrees angle of attack and static roll angles up to 42 degrees. This work is part of an on going effort to validate the RANS approach for predicting high-incidence vortical flows, with the eventual application to wing rock. The flow is unsteady and includes spiral-type vortex breakdown. The breakdown positions, mean surface pressures, rolling moments, normal forces, and streamwise center-of-pressure locations compare reasonably well with experiment. In some cases, the primary vortex suction peaks are significantly underpredicted due to grid coarseness. Nevertheless, the computations are able to predict the same nonlinear variation of rolling moment with roll angle that appeared in the experiment. This nonlinearity includes regions of local static roll instability, which is attributed to vortex breakdown.

Chaderjian, Neal M.↗
NASA NTRS1946

Stress Analysis of Columns and Beam Columns by the Photoelastic Method

Principles of similarity and other factors in the design of models for photoelastic testing are discussed. Some approximate theoretical equations, useful in the analysis of results obtained from photoelastic tests are derived. Examples of the use of photoelastic techniques and the analysis of results as applied to uniform and tapered beam columns, circular rings, and statically indeterminate frames, are given. It is concluded that this method is an effective tool for the analysis of structures in which column action is present, particularly in tapered beam columns, and in statically indeterminate structures in which the distribution of loads in the structures is influenced by bending moments due to axial loads in one or more members.

BEAM-COLUMNS - STRESSES↗
DOE OSTI2021

Light-wave control of correlated materials using quantum magnetism during time-periodic modulation of coherent transport

Light–wave quantum electronics utilizes the oscillating carrier wave to control electronic properties with intense laser pulses. Without direct light–spin interactions, however, magnetic properties can only be indirectly affected by the light electric field, mostly at later times. A grand challenge is how to establish a universal principle for quantum control of charge and spin fluctuations, which can allow for faster-than-THz clock rates. Using quantum kinetic equations for the density matrix describing non–equilibrium states of Hubbard quasiparticles, here we show that time–periodic modulation of electronic hopping during few cycles of carrier–wave oscillations can dynamically steer an antiferromagnetic insulating state into a metalic state with transient magnetization. While nonlinearities associated with quasi-stationary Floquet states have been achieved before, magneto–electronics based on quasiparticle acceleration by time–periodic multi–cycle fields and quantum femtosecond/attosecond magnetism via strongly–coupled charge–spin quantum excitations represents an alternative way of controlling magnetic moments in sync with quantum transport.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1979

General Potential Theory of Arbitrary Wing Sections

The problem of determining the two dimensional potential flow around wing sections of any shape is examined. The problem is condensed into the compact form of an integral equation capable of yielding numerical solutions by a direct process. An attempt is made to analyze and coordinate the results of earlier studies relating to properties of wing sections. The existing approximate theory of thin wing sections and the Joukowski theory with its numerous generalizations are reduced to special cases of the general theory of arbitrary sections, permitting a clearer perspective of the entire field. The method which permits the determination of the velocity at any point of an arbitrary section and the associated lift and moments is described. The method is also discussed in terms for developing new shapes of preassigned aerodynamical properties.

Theodorsen, T.↗
NASA NTRS1944

Ovalization of Tubes Under Bending and Compression

An empirical equation has been developed that gives the approximate amount of ovalization for tubes under bending loads. Tests were made on tubes in the d/t range from 6 to 14, the latter d/t ratio being in the normal landing gear range. Within the range of the series of tests conducted, the increase in ovalization due to a compression load in combination with a bending load was very small. The bending load, being the principal factor in producing the ovalization, is a rather complex function of the bending moment, d/t ratio, cantilever length, and distance between opposite bearing faces. (author)

Demer, L J↗
NASA NTRS1996

User's Manual for FOMOCO Utilities-Force and Moment Computation Tools for Overset Grids

In the numerical computations of flows around complex configurations, accurate calculations of force and moment coefficients for aerodynamic surfaces are required. When overset grid methods are used, the surfaces on which force and moment coefficients are sought typically consist of a collection of overlapping surface grids. Direct integration of flow quantities on the overlapping grids would result in the overlapped regions being counted more than once. The FOMOCO Utilities is a software package for computing flow coefficients (force, moment, and mass flow rate) on a collection of overset surfaces with accurate accounting of the overlapped zones. FOMOCO Utilities can be used in stand-alone mode or in conjunction with the Chimera overset grid compressible Navier-Stokes flow solver OVERFLOW. The software package consists of two modules corresponding to a two-step procedure: (1) hybrid surface grid generation (MIXSUR module), and (2) flow quantities integration (OVERINT module). Instructions on how to use this software package are described in this user's manual. Equations used in the flow coefficients calculation are given in Appendix A.

Chan, William M.↗