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 487 records · Page 27

New Third-Order Moments for the PBL

Turbulent convection is inherently non-local and a primary condition for a successful treatment of the PBL is a reliable model of non-locality. In the dynamic equations governing the convective flux, turbulent kinetic energy, etc, non-locality enters through the third-order moments, TOMs. Since the simplest form, the so-called down gradient approximation (DGA , severely underestimates the TOMs (by up to an order of magnitude), a more physical model is needed. In 1994, an analytical model was presented which was derived directly from the dynamical equations for the TOMs. It considerably improved the DGA but was a bit cumbersome to use. Here, we present a new analytic expression for the TOMs which is considerably simpler than the 1994 expression and which at the same time yields a much better fit to the LES data.

Canuto, V. M.↗

A Second Moment Method for k -Eigenvalue Acceleration with Continuous Diffusion and Discontinuous Transport Discretizations

The second moment method is a linear acceleration technique that couples the transport equation to a diffusion equation with transport-dependent additive closures. The resulting low-order diffusion equation can be discretized independent of the transport discretization, unlike diffusion synthetic acceleration, and is symmetric positive definite, unlike quasidiffusion. While this method has been shown to be comparable to quasidiffusion in iterative performance for fixed source and time-dependent problems, it is largely unexplored as an eigenvalue problem acceleration scheme due to the belief that the resulting inhomogeneous source makes the problem ill posed. Recently, a preliminary feasibility study was performed on the second moment method for eigenvalue problems. The results suggested comparable performance to quasidiffusion and more robust performance than diffusion synthetic acceleration. This work extends the initial study to more realistic reactor problems using state-of-the-art discretization techniques. Finally, the results in this paper show that the second moment method is more computationally efficient than its alternatives on complex reactor problems with unstructured meshes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modal analysis of turborotors using planar modes - Theory

The generalized dynamic equations of motion have been obtained by the direct stiffness method for multimass flexible rotor bearing systems including the effects of gyroscopic moments, disk skew, and rotor acceleration. A set of undamped critical speed mode shapes calculated from the average horizontal and vertical bearing stiffness is used to transform the equations of motion into a set of coupled modal equations of motion. The modal equations are coupled by the generalized bearing coefficients and the gyroscopic moments. An analysis using only undamped critical speeds or decoupled modal analysis assuming proportional damping may lead to erroneous results. This paper presents a rapid method of calculating rotor resonance speeds with their corresponding amplification factors, stability and unbalance response of turborotors. Examples of the application of this modal approach are presented and results are compared to those of other methods such as matrix transfer analysis.

Gunter, E. J.↗

New Third-Order Moments for the CBL

Turbulent convection is an inherently non-local phenomenon and a primary condition for a successful treatment of the CBL (convective boundary layer) is a reliable model of non-locality. In the dynamic equations governing the convective flux, the turbulent kinetic energy, etc., non-locality is represented by the third-order moments, TOMs. Since the simplest form, the so-called down gradient approximation (DGA), severely underestimates the TOMs (up to an order of magnitude), a more physical model is needed. In 1994, an analytical model was presented which was derived directly from the dynamical equations for the TOMs. It considerably improved the DGA but was a bit cumbersome to use and, more importantly, it was based on the quasi-normal (QN) approximation for the fourth-order moments. Here, we present a new analytic expression for the TOMs which is structurally simpler than the 1994 expression and which avoids the QN approximation. The resulting fit to the LES data is superior to that of the 1994 model.

Canuto, V. M.↗

MOCMC: Method of Characteristics Moment Closure, a Numerical Method for Covariant Radiation Magnetohydrodynamics

In this work, we present a conservative numerical method for radiation magnetohydrodynamics with frequency-dependent full transport in stationary spacetimes. This method is stable and accurate for both large and small optical depths and radiation pressures. The radiation stress–energy tensor is evolved in flux-conservative form, and closed with a swarm of samples that each transport a multigroup representation of the invariant specific intensity along a null geodesic. In each zone, the enclosed samples are used to efficiently construct a Delaunay triangulation of the unit sphere in the comoving frame, which in turn is used to calculate the Eddington tensor, average source terms, and adaptively refine the sample swarm. Furthermore, radiation four-forces are evaluated in the moment sector in a semi-implicit fashion. The radiative transfer equation is solved in invariant form deterministically for each sample. Since each sample carries a discrete representation of the full spectrum, the cost of evaluating the transport operator is independent of the number of frequency groups, representing a significant reduction of algorithmic complexity for transport in frequency-dependent problems. The major approximation we make in this work is performing scattering in an angle-averaged way. Local adaptivity in samples also makes this scheme more amenable to nonuniform meshes than a traditional Monte Carlo method. We describe the method and present results on a suite of test problems. We find that Method of Characteristics Moment Closure converges at least as ~N -1 , rather than the canonical Monte Carlo N -1/2 , where N is the number of samples per zone.

79 ASTRONOMY AND ASTROPHYSICS↗

Research study on stabilization and control modern sampled-data control theory

The methods of continuous and discrete describing function analysis were applied to predicting the existence of self-sustained oscillations in the single-axis model of the large space telescope system with nonlinear control moment gyroscope friction characteristics. It is shown that the stability equations may be solved by a numerical-iterative technique using the describing function analysis, instead of the usual graphical methods. The numerical method is found to be effective in leading to a convergent solution rapidly, with an appropriate guess of the initial condition.

Kuo, B. C.↗

Analysis of shuttle oscillation in the Mach Number = 1.7 to Mach Number = 1.0 range

The first five flights of the Space Transportation System (STS) have exhibited unexpected lateral oscillations in the Mach Number = 1.7 to Mach Number = 1.0 region of the descent trajectory. These oscillations can be split into two parts: a predominantly rolling oscillation referred to as the quarter hertz motion and a long-priod, large-amplitude oscillation in sideslips, roll rate, and yaw rate. These motions are analyzed in this paper in order to determine their source and the mechanism driving them. They will first be examined to determine how much of the motion can be explained using the set of equations that contain linear aerodynamic parameters. Then differences between measured moments and those calculated using linear aerodynamics which are considered to be due to unmodeled moments will be exained. Possible error sources such as not inluding nonlinear aerodynamics in the vehicle model, measurement errors, and the effect of pressure gradients will be considered. Based on these analyses, conclusions as to the nature and cause of the Shuttle lateral oscillations will be drawn.

Suit, W. T.↗

Pulsars as cosmic ray particle accelerators: Proton orbits

Proton orbits are calculated in the electromagnetic vacuum field of a magnetic point dipole rotating with its angular velocity omega perpendicular to its dipole moment mu by numerical integration of the Lorentz-Dirac equation. Trajectories are shown and discussed for various initial conditions. A critical surface is shown separating initial positions of protons which finally hit the pulsar in the polar region from those which finally recede to infinity.

Thielheim, K. O.↗

Analog simulation of stochastic satellite response

A numerical study was performed on the effects of stochastic geometry and random environmental torques on the pointing accuracy of spinning and three-axis stabilized satellites. The Euler equations for the motions of satellites yielded the stochastic principal moments and a Fokker-Planck analog simulation was employed to model the response vector of a rigid satellite. White noise inputs were obtained from a random number generator. Attention was given to both fast and slow spin satellites. Responses grew with time in all cases. The growth was initially exponential and eventually reached stable values. In comparison with other theoretical techniques, the Fokker-Planck model detailed the most responses and showed that spin-stabilization is preferable for satellites with low inertial noise levels.

Huang, T. C.↗

A Comparative Study of Some Dynamic Stall Models

Three semi-empirical aerodynamic stall models are compared with respect to their lift and moment hysteresis loop prediction, limit cycle behavior, easy implementation, and feasibility in developing the parameters required for stall flutter prediction of advanced turbines. For the comparison of aeroelastic response prediction including stall, a typical section model and a plate structural model are considered. The response analysis includes both plunging and pitching motions of the blades. In model A, a correction to the angle of attack is applied when the angle of attack exceeds the static stall angle. In model B, a synthesis procedure is used for angles of attack above static stall angles and the time history effects are accounted through the Wagner function. In both models the life and moment coefficients for angle of attack below stall are obtained from tabular data for a given Mach number and angle of attack. In model C, referred to an the ONERA model, the life and moment coefficients are given in the form of two differential equations, one for angles below stall, and the other for angles above stall. The parameters of those equations are nonlinear functions of the angle of attack.

Reddy, T. S. R.↗

Maxwell's second- and third-order equations of transfer for non-Maxwellian gases

Condensed algebraic forms for Maxwell's second- and third-order equations of transfer are developed for the case of molecules described by either elastic hard spheres, inverse-power potentials, or by Bird's variable hard-sphere model. These hardly reduced, yet exact, equations provide a new point of origin, when using the moment method, in seeking approximate solutions in the kinetic theory of gases for molecular models that are physically more realistic than that provided by the Maxwell model. An important by-product of the analysis when using these second- and third-order relations is that a clear mathematical connection develops between Bird's variable hard-sphere model and that for the inverse-power potential.

Baganoff, D.↗

Approximate relations and charts for low-speed stability derivatives of swept wings

Contains derivations, based on a simplified theory, of approximate relations for low-speed stability derivatives of swept wings. Method accounts for the effects and, in most cases, taper ratio. Charts, based on the derived relations, are presented for the stability derivatives of untapered swept wings. Calculated values of the derivatives are compared with experimental results.

THEORIES - WINGS SWEPTBACK & SWEPTFORWARD↗

Simulation of Gas Mixtures in RF Discharges for Nitride Deposition

Tungsten and titanium nitride films have long been grown by chemical vapor deposition (CVD) methods. However, there has been recent interest in low temperature growth using plasma enhanced CVD. For the present work, we focus on the radio frequency (BE) discharge characteristics of gas mixtures used in nitride deposition (for example, WF6 and ammonia). Because the radial variations for a standard 200 mm, parallel plate reactor are limited to a small zone near the edges of the electrodes, a 1-D (one-dimensional) analysis is considered. This model consists of a self-consistent, 3-D (three-dimensional) moment fluid simulation that solves the continuity, momentum, and energy equations for neutral and charged species. The results in terms of plasma structure, radical concentrations, and local deposition rate will be presented. We will also compare the 1-D results with those obtained from a 2-D hybrid plasma equipment model (HPEM) developed at the University of Illinois.

Hwang, Helen H.↗

Application Program Interface for the Orion Aerodynamics Database

The Application Programming Interface (API) for the Crew Exploration Vehicle (CEV) Aerodynamic Database has been developed to provide the developers of software an easily implemented, fully self-contained method of accessing the CEV Aerodynamic Database for use in their analysis and simulation tools. The API is programmed in C and provides a series of functions to interact with the database, such as initialization, selecting various options, and calculating the aerodynamic data. No special functions (file read/write, table lookup) are required on the host system other than those included with a standard ANSI C installation. It reads one or more files of aero data tables. Previous releases of aerodynamic databases for space vehicles have only included data tables and a document of the algorithm and equations to combine them for the total aerodynamic forces and moments. This process required each software tool to have a unique implementation of the database code. Errors or omissions in the documentation, or errors in the implementation, led to a lengthy and burdensome process of having to debug each instance of the code. Additionally, input file formats differ for each space vehicle simulation tool, requiring the aero database tables to be reformatted to meet the tool s input file structure requirements. Finally, the capabilities for built-in table lookup routines vary for each simulation tool. Implementation of a new database may require an update to and verification of the table lookup routines. This may be required if the number of dimensions of a data table exceeds the capability of the simulation tools built-in lookup routines. A single software solution was created to provide an aerodynamics software model that could be integrated into other simulation and analysis tools. The highly complex Orion aerodynamics model can then be quickly included in a wide variety of tools. The API code is written in ANSI C for ease of portability to a wide variety of systems. The input data files are in standard formatted ASCII, also for improved portability. The API contains its own implementation of multidimensional table reading and lookup routines. The same aerodynamics input file can be used without modification on all implementations. The turnaround time from aerodynamics model release to a working implementation is significantly reduced

Robinson, Philip E.↗

Investigation of Multiphase Flow in a Packed Bed Reactor Under Microgravity Conditions

In this paper we study the two-phase flow phenomena in a packed bed reactor using an integrated experimental and numerical method. The cylindrical bed is filled with uniformly sized spheres. In the experiment water and air are injected into the bed simultaneously. The pressure distribution along the bed will be measured. The numerical simulation is based on a two-phase flow solver which solves the Navier-Stokes equations on Cartesian grids. A novel coupled level set and moment of fluid method is used to construct the interface. A sequential method is used to position spheres in the cylinder. Preliminary experimental results showed that the tested flow rates resulted in pulse flow. The numerical simulation revealed that air bubbles could merge into larger bubbles and also could break up into smaller bubbles to pass through the pores in the bed. Preliminary results showed that flow passed through regions where the porosity is high. Comparison between the experimental and numerical results in terms of pressure distributions at different flow injection rates will be conducted. Comparison of flow phenomena under terrestrial gravity and microgravity will be made.

computational fluid dynamics↗

A Mercury Lander Mission Concept Study for the Next Decadal Survey

Mariner 10 provided our first closeup reconnaissance of Mercury during its three flybys in 1974 and 1975. MESSENGER’s 2011–2015 orbital investigation enabled numerous discoveries, several of which led to substantial or complete changes in our fundamental understanding of the planet. Among these were the unanticipated, widespread presence of volatile elements (e.g., Na, K, S); a surface with extremely low Fe abundance whose darkening agent is likely C; a previously unknown landform—hollows— that may form by volatile sublimation from within rocks exposed to the harsh conditions on the surface; a history of expansive effusive and explosive volcanism; substantial radial contraction of the planet from interior cooling; offset of the dipole moment of the internal magnetic field northward from the geographic equator by ~20% of the planet’s radius; crustal magnetization, attributed at least in part to an ancient field; unexpected seasonal variability and relationships among exospheric species and processes; and the presence in permanently shadowed polar terrain of water ice and other volatile materials, likely to include complex organic compounds. Mercury’s highly chemically reduced and unexpectedly volatile-rich composition is unique among the terrestrial planets and was not predicted by earlier hypotheses for the planet’s origin. As an end-member of terrestrial planet formation, Mercury holds unique clues about the original distribution of elements in the earliest stages of the Solar System and how planets (and exoplanets) form and evolve in close proximity to their host stars. The BepiColombo mission promises to expand our knowledge of this planet and to shed light on some of the mysteries revealed by the MESSENGER mission. However, several fundamental science questions raised by MESSENGER’s pioneering exploration of Mercury can only be answered with in situ measurements from the planet’s surface.

Ernst, C. M.↗

Theoretical Evaluation of the Pressures, Forces, and Moments at Hypersonic Speeds Acting on Arbitrary Bodies of Revolution Undergoing Separate and Combined Angle-of-attack and Pitching Motions

Equations based on Newtonian impact theory have been derived and a computational procedure developed with the aid of several design-type charts which enable the determination of the aerodynamic forces and moments acting on arbitrary bodies of revolution undergoing either separate or combined angle-of-attack and pitching motions. Bodies with axially increasing and decreasing cross-sectional area distributions are considered; nose shapes may be sharp, blunt, or flat faced. The analysis considers variations in angle of attack from -90 degrees to 90 degrees and allows for both positive and negative pitching rates of arbitrary magnitude. The results are also directly applicable to bodies in either separate or combined sideslip and yawing maneuvers.

Margolis, Kenneth↗

Contact dynamics math model

The Space Station Mechanism Test Bed consists of a hydraulically driven, computer controlled six degree of freedom (DOF) motion system with which docking, berthing, and other mechanisms can be evaluated. Measured contact forces and moments are provided to the simulation host computer to enable representation of orbital contact dynamics. This report describes the development of a generalized math model which represents the relative motion between two rigid orbiting vehicles. The model allows motion in six DOF for each body, with no vehicle size limitation. The rotational and translational equations of motion are derived. The method used to transform the forces and moments from the sensor location to the vehicles' centers of mass is also explained. Two math models of docking mechanisms, a simple translational spring and the Remote Manipulator System end effector, are presented along with simulation results. The translational spring model is used in an attempt to verify the simulation with compensated hardware in the loop results.

Glaese, John R.↗