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 955 records · Page 53

Numerical solutions of the linearized Euler equations for unsteady vortical flows around lifting airfoils

A linearized unsteady aerodynamic analysis is presented for unsteady, subsonic vortical flows around lifting airfoils. The analysis fully accounts for the distortion effects of the nonuniform mean flow on the imposed vortical disturbances. A frequency domain numerical scheme which implements this linearized approach is described, and numerical results are presented for a large variety of flow configurations. The results demonstrate the effects of airfoil thickness, angle of attack, camber, and Mach number on the unsteady lift and moment of airfoils subjected to periodic vortical gusts. The results show that mean flow distortion can have a very strong effect on the airfoil unsteady response, and that the effect depends strongly upon the reduced frequency, Mach number, and gust wave numbers.

Scott, James R.↗

Numerical solutions of the linearized Euler equations for unsteady vortical flows around lifting airfoils

A linearized unsteady aerodynamic analysis is presented for unsteady, subsonic vortical flows around lifting airfoils. The analysis fully accounts for the distortion effects of the nonuniform mean flow on the imposed vortical disturbances. A frequency domain numerical scheme which implements this linearized approach is described, and numerical results are presented for a large variety of flow configurations. The results demonstrate the effects of airfoil thickness, angle of attack, camber, and Mach number on the unsteady lift and moment of airfoils subjected to periodic vortical gusts. The results show that mean flow distortion can have a very strong effect on the airfoil unsteady response, and that the effect depends strongly upon the reduced frequency, Mach number, and gust wave numbers.

Scott, James R.↗

Absolute intensity measurement of the 4-0 vibration-rotation band of carbon monoxide

The absolute intensity of the 4-0 vibration band of CO is measured in spectra obtained using a 25-m base-path multiple-traversal absorption cell and a 5-m scanning spectrometer. The intensities of individual vibration-rotation lines in this band are determined from measurements of their equivalent widths, and absolute values for the rotationless transition moment and the vibration-rotation interaction factor are derived from the measured line strengths. The experimentally obtained vibration-rotation function is compared with a theoretical curve; agreement between theory and experiment is found to be good for the P-branch but poor for the R-branch. It is noted that numerical solutions to the radial Schroedinger equation lead to vibration-rotation function values that are in good agreement with the experiment.

Chackerian, C., Jr.↗

Uncertainty Quantification of Geophysical Inversion Using Stochastic Partial Differential Equations (LDRD #218329)

This report summarizes work completed under the Laboratory Directed Research and Development (LDRD) project "Uncertainty Quantification of Geophysical Inversion Using Stochastic Differential Equations." Geophysical inversions often require computationally expensive algorithms to find even one solution, let alone propagating uncertainties through to the solution domain. The primary purpose of this project was to find more computationally efficient means to approximate solution uncertainty in geophysical inversions. We found multiple computationally efficient methods of propagating Earth model uncertainty into uncertainties in solutions of full waveform seismic moment tensor inversions. However, the optimum method of approximating the uncertainty in these seismic source solutions was to use the Karhunen-Love theorem with data misfit residuals. This method was orders of magnitude more computationally efficient than traditional Monte Carlo methods and yielded estimates of uncertainty that closely approximated those of Monte Carlo. We will summarize the various methods we evaluated for estimating uncertainty in seismic source inversions as well as work toward this goal in the realm of 3-D seismic tomographic inversion uncertainty.

58 GEOSCIENCES↗

Attitude control of a spinning Skylab.

Active attitude control of a spinning Skylab is analyzed to determine a simple control law that will provide a satisfactory response, considering the dynamics of attached flexible appendages. A simplified model of the complex Skylab vehicle is selected to make it analytically tractable. The vehicle is modeled as a single rigid-core body with two attached flexible massless booms having masses on their tips. The equations of motion describing the attitude dynamics of the model are presented as a linear matrix-differential equation. The states of the vehicle are small perturbations about its steady-state spin. An analysis is performed to determine the domain of stability. Next, attitude dynamics are analyzed; both frequency domain (parameter plane) and time domain (an optimal linear quadratic loss program) techniques are compared. An analysis of the nonlinear effect of control torque saturation of Skylab's control moment gyroscopes is discussed. The results of the analysis compare favorably with a large-scale digital simulation of the Skylab.

Seltzer, S. M.↗

Electron Crescent Distributions as a Manifestation of Diamagnetic Drift in an Electron-Scale Current Sheet: Magnetospheric Multiscale Observations Using New 7.5 ms Fast Plasma Investigation Moments

We report Magnetospheric Multiscale observations of electron pressure gradient electric fields near a magnetic reconnection diffusion region using a new technique for extracting 7.5 ms electron moments from the Fast Plasma Investigation. We find that the deviation of the perpendicular electron bulk velocity from E × B drift in the interval where the out-of-plane current density is increasing can be explained by the diamagnetic drift. In the interval where the out-of-plane current is transitioning to in-plane current, the electron momentum equation is not satisfied at 7.5 ms resolution.

fast plasma↗

An Improved Theoretical Aerodynamic Derivatives Computer Program for Sounding Rockets

The paper outlines a Theoretical Aerodynamic Derivatives (TAD) computer program for computing the aerodynamics of sounding rockets. TAD outputs include normal force, pitching moment and rolling moment coefficient derivatives as well as center-of-pressure locations as a function of the flight Mach number. TAD is applicable to slender finned axisymmetric vehicles at small angles of attack in subsonic and supersonic flows. TAD improvement efforts include extending Mach number regions of applicability, improving accuracy, and replacement of some numerical integration algorithms with closed-form integrations. Key equations used in TAD are summarized and typical TAD outputs are illustrated for a second-stage Tomahawk configuration.

Barrowman, J. S.↗

Efficient modeling of interconnects and capacitive discontinuities in high-speed digital circuits

Modeling of interconnects and associated discontinuities with the recent advances high-speed digital circuits has gained a considerable interest over the last decade although the theoretical bases for analyzing these structures were well-established as early as the 1960s. Ongoing research at the present time is focused on devising methods which can be applied to more general geometries than the ones considered in earlier days and, at the same time, improving the computational efficiency and accuracy of these methods. In this thesis, numerically efficient methods to compute the transmission line parameters of a multiconductor system and the equivalent capacitances of various strip discontinuities are presented based on the quasi-static approximation. The presented techniques are applicable to conductors embedded in an arbitrary number of dielectric layers with two possible locations of ground planes at the top and bottom of the dielectric layers. The cross-sections of conductors can be arbitrary as long as they can be described with polygons. An integral equation approach in conjunction with the collocation method is used in the presented methods. A closed-form Green's function is derived based on weighted real images thus avoiding nested infinite summations in the exact Green's function; therefore, this closed-form Green's function is numerically more efficient than the exact Green's function. All elements associated with the moment matrix are computed using the closed-form formulas. Various numerical examples are considered to verify the presented methods, and a comparison of the computed results with other published results showed good agreement.

Oh, K. S.↗

Lunar Rotation and the Lunar Interior

Variations in rotation and orientation of the Moon are sensitive to solid-body tidal dissipation, dissipation due to relative motion at the fluid-core/ solid-mantle boundary, and tidal Love number k2. There is weaker sensitivity to flattening of the core-mantle boundary (CMB) and fluid core moment of inertia. Accurate Lunar Laser Ranging (LLR) measurements of the distance from observatories on the Earth to four retroreflector arrays on the Moon are sensitive to lunar rotation and orientation variations and tidal displacements. Past solutions using the LLR data have given results for dissipation due to solid-body tides and fluid core plus Love number. Past detection of CMB flattening has been marginal but is improving, while direct detection of the core moment has not yet been achieved. Three decades of Lunar Laser Ranging (LLR) data are analyzed using a weighted least-squares approach. The lunar solution parameters include dissipation at the fluid-core/solid-mantle boundary, tidal dissipation, dissipation-related coefficients for rotation and orientation terms, potential Love number k2, a correction to the constant term in the tilt of the equator to the ecliptic which is meant to approximate the influence of core-mantle boundary flattening, and displacement Love numbers h2 and l2. Several solutions, with different combinations of solution parameters and constraints, are considered.

Williams, J. G.↗

Formation of the low-mass solar nebula

We study an accretional stage of the formation and early evolution of the solar nebula with relatively small angular momentum. We investigate the evolution of the disk and its vertical structure, particularly the shock front between disk and infalling material. Calculations start at a moment when a low-mass star-like core surrounded by small embryo disk have been formed at the center of the presolar nebula and the bulk of mass remained in the envelope. The forming solar nebula is approximated as a thin viscous disk surrounded by accreting envelope. The distribution of temperature in the infalling envelope is determined by solving spherically symmetric equations of radiative transfer. As the energy source, we take into account all energy released within the centrifugal radius of the infalling matter. Other aspects of this study are discussed.

Ruzmaikina, T. V.↗

High-Fidelity Aeropropulsive Assessment of Distributed Electric Propulsion and Boundary Layer Ingestion for the SUSAN Electrofan

This paper presents a high-fidelity aeropropulsive assessment of the SUSAN Electrofan’s distributed electric propulsion (DEP) and boundary layer ingestion (BLI) aft fuselage turbofan technologies. This is achieved through the application of computational tools for aerodynamic and aeropropulsive design and analysis based on the Reynolds-averaged Navier-Stokes (RANS) equations. Representative aircraft models are first developed in 3D for the SUSAN Electrofan concept, a non-DEP variant, and a year 2020 technology level reference aircraft based on the Boeing 737-8. These include airframes developed through the application of RANS-based aerodynamic shape optimization, with the objective of drag minimization at cruise while subject to constant lift and zero pitching moment constraints. Wing propulsors are also included and modeled as free-flying ducted fans. These are developed through RANS-based aeropropulsive shape optimization, with the objective of shaft power minimization while subject to constant net thrust and maximum fan face Mach number constraints. An iterative power balance method is then used to trim each of the decoupled aircraft configurations, providing estimates of relative aerodynamic performance, shaft power, and range efficiency. Results also include power on and off analyses, which indicate reduced thrust-induced drag for the SUSAN Electrofan’s DEP and BLI aft fuselage turbofan propulsors when compared to those of the Boeing 737-8-like reference aircraft.

CAS↗

The planetary magnetic field and magnetosphere of Jupiter - Pioneer 10

Data obtained by the Pioneer 10 vector helium magnetometer are presented along with models of the intrinsic magnetic field of Jupiter and its magnetosphere. Data acquired between 2.84 and 6.0 Jupiter radii, where the intensity of the planetary field ranged between 1900 and 18,400 gamma, were used to develop a six-parameter eccentric dipole model of the field. The dipole so derived has a moment of 4.0 G (R sub J) cubed and a tilt angle with respect to Jupiter's rotation axis of 11 deg. A model of the Jovian magnetosphere is presented in which the essential feature is an eastward current sheet that forms an annulus with Jupiter at the center. At large distances from the planet the current sheet is nearly parallel to Jupiter's equator but, in general, does not lie in it. The current sheet is warped, so that it is above the equator on one side and below it on the other. The current sheet rotates with the planet, more or less like a rigid body, this behavior causes an apparent up and down motion and periodic crossings of the current sheet by Pioneer.

Smith, E. J.↗

A lifting surface theory in rotational flow

The partial differential equation for small disturbance steady rotational flow in three dimensions is solved through an integral equation approach. The solution is obtained by using the method of weighted residuals. Specific applications are directed to wings in nonuniform subsonic parallel streams with velocity varying in vertical and spanwise directions and to airfoils in nonuniform freestream. Comparison with limited known results indicates that the present method is reasonably accurate. Numerical results for the lifting pressure of airfoil, lift, induced drag, and pitching moments of airfoil, lift, induced drag, and pitching moments of elliptic, rectangular, and delta wings in a jet, wake, or monotonic sheared stream are presented. It is shown that, in addition to the effect of local dynamic pressures, a positive velocity gradient tends to enhance the lift.

Shiau, M. J.↗

Effect of canard deflection on close-coupled canard-wing-body aerodynamics

The thin-layer Navier-Stokes equations are solved for the flow about a canard-wing-body configuration at transonic Mach numbers of 0.85 and 0.90, angles of attack from -4 to 10 degrees and canard deflection angles from -10 to +10 degrees. Effects of canard deflection on aerodynamic performance, including canard-wing vortex interaction, are investigated. Comparisons with experimental measurements of surface pressures, lift, drag and pitching moments are made to verify the accuracy of the computations. The results of the study show that the deflected canard downwash not only influences the formation of the wing leading-edge vortex, but can cause the formation of an unfavorable vortex on the wing lower surface as well.

Tu, Eugene L.↗

High order interpolation of magnetic fields with vector potential reconstruction for particle simulations

We propose a method for interpolating divergence-free continuous magnetic fields via vector potential reconstruction using Hermite interpolation, which ensures high-order continuity for applications requiring adaptive, high-order ordinary differential equation (ODE) integrators, such as the Dormand-Prince method. The method provides C(m) continuity and achieves high-order accuracy, making it particularly suited for particle trajectory integration and Poincaré section analysis under optimal integration order and timestep adjustments. Through numerical experiments, we demonstrate that the Hermite interpolation method preserves volume and continuity, which are critical for conserving toroidal canonical momentum and magnetic moment in guiding center simulations, especially over long-term trajectory integration. Furthermore, we analyze the impact of insufficient derivative continuity on Runge-Kutta schemes and show how it degrades accuracy at low error tolerances, introducing discontinuity-induced truncation errors. Lastly, we demonstrate performant Poincaré section analysis in two relevant settings of field data collocated from finite element meshes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The motion and stability of spin stabilized spacecraft with hinged appendages

The dynamics and stability of a spin stabilized spacecraft with a hinged appendage system are treated analytically and numerically. The hinged system consists of a central hub with masses attached to massless booms of fixed length whose orientation relative to the main part can change. The equations of motion for the hinged system, with viscous damping at both hinge points, are linearized about the nominal equilibrium position where the booms are orthogonal to the nominal spin axis. Analytic stability criteria indicate that hinge damping must be present and that for limiting cases, where the spin axis is an axis of symmetry, certain conditions relating the hinge point offset coordinates to the moment of inertia ratio and end masses must be satisfied. The hinge damping is always required for the nominal deployment of hinge members.

Sellappan, R.↗

Orbiting chains and rings

Two very different, highly flexible, space structures proposed during the last decade which involve analytical solution of certain partial differential equations are discussed. The first structure is an array, or hanging chain of aluminum beads which would serve as a convenient communicator if, under the influence of the Earth's gravity gradient, it assumes a local vertical orientation as it circles the Earth. Two passive schemes have been proposed for damping the rigid modes: twist the wire at the ends to provide non-zero moment of inertia about the vertical, thereby inducing relative motion of the two tips during rigid pitch or roll, and thus exercising a damper; and introduce weak lossy springs between the end sections and the main section; thereby providing linear coupling between the springs and all the in-plane (pitch) modes. The second structure is a complete ring of satellites cabled together at synchronous altitude. By a slight increase in altitude, the cable is in tension. This configuration is, however, unstable, and an active feedback control scheme is required to stabilize it. A possible scheme involves local cable length adjustment based on measurement of local altitude and shape variations and their rates.

Breakwell, J. V.↗

Turbulent transport in premixed flames

Simulations of planar, premixed turbulent flames with heat release were used to study turbulent transport. Reynolds stress and Reynolds flux budgets were obtained and used to guide the investigation of important physical effects. Essentially all pressure terms in the transport equations were found to be significant. In the Reynolds flux equations, these terms are the major source of counter-gradient transport. Viscous and molecular terms were also found to be significant, with both dilatational and solenoidal terms contributing to the Reynolds stress dissipation. The BML theory of premixed turbulent combustion was critically examined in detail. The BML bimodal pdf was found to agree well with the DNS data. All BML decompositions, through the third moments, show very good agreement with the DNS results. Several BML models for conditional terms were checked using the DNS data and were found to require more extensive development.

Rutland, C. J.↗