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 919 records · Page 51

Gyrodampers for large space structures

The problem of controlling the vibrations of a large space structures by the use of actively augmented damping devices distributed throughout the structure is addressed. The gyrodamper which consists of a set of single gimbal control moment gyros which are actively controlled to extract the structural vibratory energy through the local rotational deformations of the structure, is described and analyzed. Various linear and nonlinear dynamic simulations of gyrodamped beams are shown, including results on self-induced vibrations due to sensor noise and rotor imbalance. The complete nonlinear dynamic equations are included. The problem of designing and sizing a system of gyrodampers for a given structure, or extrapolating results for one gyrodamped structure to another is solved in terms of scaling laws. Novel scaling laws for gyro systems are derived, based upon fundamental physical principles, and various examples are given.

Aubrun, J. N.↗

The adiabatic motion of charged dust grains in rotating magnetospheres

Adiabatic equations of motion are derived for the micrometer-sized dust grains detected in the Jovian and Saturn magnetospheres by the Pioneer 10 and 11 spacecraft. The adiabatic theory of charged particle motion is extended to the case of variable grain charge. Attention is focused on the innermost and outermost limits to the grain orbit evolution, with all orbits tending to become circular with time. The parameters such as the center equation of motion, the drift velocity, and the parallel equation of motion are obtained for grains in a rotating magnetosphere. Consideration is given to the effects of periodic grain charge-discharge, which are affected by the ambient plasma properties and the grain plasma velocity. The charge-discharge process at the gyrofrequency is determined to eliminate the invariance of the magnetic moment and cause the grain to exhibit radial movement. The magnetic moment increases or decreases as a function of the gyrophase of the charge variation.

Northrop, T. G.↗

On the application of transonic similarity rules to wings of finite span

The transonic aerodynamic characteristics of wings of finite span are discussed from the point of view of a unified small perturbation theory for subsonic, transonic, and supersonic flows about thin wings. This approach avoids certain ambiguities which appear if one studies transonic flows by means of equations derived under the more restrictive assumption that the local velocities are everywhere close to sonic velocity. The relation between the two methods of analysis of transonic flow is examined, the similarity rules and known solutions of transonic flow theory are reviewed, and the asymptotic behavior of the lift, drag, and pitching-moment characteristics of wings of large and small aspect ratio is discussed. It is shown that certain methods of data presentation are advantageous for the effective display of these characteristics.

Spreiter, John R↗

Analysis of the effects of wing interference on the tail contributions to the rolling derivatives

An analysis of the effects of wing interference on the tail contributions to the rolling stability derivatives of complete airplane configurations is made by calculating the angularity of the air stream at the vertical tail due to rolling and determining the resulting forces and moments. Some of the important factors which affect the resultant angularity on the vertical tail are wing aspect ratio and sweepback, vertical-tail span, and considerations associated with angle of attack and airplane geometry. Some calculated sidewash results for a limited range of plan forms and vertical-tail sizes are presented. Equations taking into account the sidewash results are given for determining the tail contributions to the rolling derivatives. Comparisons of estimated and experimental results indicate that a consideration of wing interference effects improves the estimated values of the tail contributions to the rolling derivatives and that fair agreement with available experimental data is obtained.

Michael, William H , Jr↗

A multi-ion generalized transport model of the polar wind

The higher-order generalizations of the equations of standard hydrodynamics, known collectively as generalized transport theories, have been used since the early 1980s to describe the terrestrial polar wind. Inherent in the structure of generalized transport theories is the ability to describe not only interparticle collisions but also certain non-Maxwellian processes, such as heat flow and viscous stress, that are characteristic of any plasma flow that is not collision dominated. Because the polar wind exhibits a transition from collision-dominated to collisionless flow, generalized transport theories possess advantages for polar wind modeling not shared by either collision-dominated models (such as standard hydrodynamics) or collisionless models (such as those based on solving the collisionless Boltzmann equation). In general, previous polar wind models have used generalized transport equations to describe electrons and only one species of ion (H(+)). If other ion species were included in the models at all, it was in a simplified or semiempirical manner. The model described in this paper is the first polar wind model that uses a generalized transport theory (bi-Maxwellian-based 16-moment theory) to describe all of the species, both major and minor, in the polar wind plasma. In the model, electrons and three ion species (H(+), He(+), O(+)) are assumed to be major and several ion species are assumed to be minor (NO(+), Fe(+), O(++)). For all species, a complete 16-moment transport formulation is used, so that profiles of density, drift velocity, parallel and perpendicular temperatures, and the field-aligned parallel and perpendicular energy flows are obtained. In the results presented here, emphasis is placed on describing those constituents of the polar wind that have received little attention in past studies. In particular, characteristic solutions are presented for supersonic H(+) outflow and for both supersonic and subsonic outflows of the major ion He(+). Solutions are also presented for various minor ions, both atomic and molecular and both singly and multiply charged.

Demars, H. G.↗

Application of AWE Along with a Combined FEM/MoM Technique to Compute RCS of a Cavity-Backed Aperture in an Infinite Ground Plane Over a Frequency Range

A hybrid Finite Element Method (FEM)/Method of Moments (MoM) technique in conjunction with the Asymptotic Waveform Evaluation (AWE) technique is applied to obtain radar cross section (RCS) of a cavity-backed aperture in an infinite ground plane over a frequency range. The hybrid FEM/MoM technique when applied to the cavity-backed aperture results in an integro-differential equation with electric field as the unknown variable, the electric field obtained from the solution of the integro-differential equation is expanded in Taylor series. The coefficients of the Taylor series are obtained using the frequency derivatives of the integro-differential equation formed by the hybrid FEM/MoM technique. The series is then matched via the Pade approximation to a rational polynomial, which can be used to extrapolate the electric field over a frequency range. The RCS of the cavity-backed aperture is calculated using the electric field at different frequencies. Numerical results for a rectangular cavity, a circular cavity, and a material filled cavity are presented over a frequency range. Good agreement between AWE and the exact solution over the frequency range is obtained.

Reddy, C.J.↗

Conservative DG method for the micro-macro decomposition of the Vlasov–Poisson–Lenard–Bernstein model

The micro-macro (mM) decomposition approach is considered for the numerical solution of the Vlasov– Poisson–Lenard–Bernstein (VPLB) system, which is relevant for plasma physics applications. In the mM approach, the kinetic distribution function is decomposed as $f$ = $\mathscr{E}$[$ρ_f$] + g , where $\mathscr{E}$ is a local equilibrium distribution, depending on the macroscopic moments $ρ_f$ = $∫_{\mathbb{R}}$ $efdv$ = $\langle$$ef$$\rangle$$_{\mathbb{R}}$, where e = (1, $v$, $\frac{1}{2}$ $v^2$) T , and $g$, the microscopic distribution, is defined such that $\langle$$eg$$\rangle$$_{\mathbb{R}}$ = 0. We aim to design numerical methods for the mM decomposition of the VPLB system, which consists of coupled equations for $ρ_f$ and $g$. To this end, we use the discontinuous Galerkin (DG) method for phase-space discretization, and implicit-explicit (IMEX) time integration, where the phase-space advection terms are integrated explicitly and the collision operator is integrated implicitly. We give special consideration to ensure that the resulting mM method maintains the $\langle$$eg$$\rangle$$_{\mathbb{R}}$ = 0 constraint, which may be necessary for obtaining (i) satisfactory results in the collision dominated regime with coarse velocity resolution, and (ii) unambiguous conservation properties. The constraint-preserving property is achieved through a consistent discretization of the equations governing the micro and macro components. Here, we present numerical results that demonstrate the performance of the mM method. The mM method is also compared against a corresponding DG-IMEX method solving directly for $f$.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Opening of the magnetic field lines in a fast rotating magnetosphere, with an application to Jupiter

A simple model of a magnetosphere rotating around a rotating Jupiter-like object with a spin-aligned dipolar moment mu is considered. A low-energy plasma is released by inner sources located beyond the corotation radius and is diffused outward through closed lines. The lines form a closed thin equatorial disk, maintaining a quasi-static balance between the centrifugal force and the magnetic tension. Solutions are obtained for the analytical equations characterizing the plasma in the closed region; the structure of the magnetic field outside the disk; and the self-consistent value of the critical radius r sub zero. The analytical solutions are compared to the observed properties of the Jovian magnetosphere, and the results are discussed in detail.

Aly, J. J.↗

Navier-Stokes analysis of a delta wing in static and dynamic roll

The three-dimensional, Reynolds-averaged, Navier-Stokes equations are used to numerically simulate non-steady high-incidence vortical flow about a 65 degree sweep delta wing under static roll and forced periodic roll motions. These computations have been previously reported in the literature, where emphasis was placed on validating the computational results with nonsteady experimental surface pressures, forces and moments. The emphasis of this research is to revisit these earlier computations and use nonsteady numerical flow visualization to analyze the nonlinear surface flow and vortex breakdown dynamics. Most notable are the periodic formation of surface-flow separation lines, and the periodic formation of vortex breakdown near the delta wing trailing edge.

Chaderjian, N.↗

Chimera Grid Tools

Chimera Grid Tools (CGT) is a software package for performing computational fluid dynamics (CFD) analysis utilizing the Chimera-overset-grid method. For modeling flows with viscosity about geometrically complex bodies in relative motion, the Chimera-overset-grid method is among the most computationally cost-effective methods for obtaining accurate aerodynamic results. CGT contains a large collection of tools for generating overset grids, preparing inputs for computer programs that solve equations of flow on the grids, and post-processing of flow-solution data. The tools in CGT include grid editing tools, surface-grid-generation tools, volume-grid-generation tools, utility scripts, configuration scripts, and tools for post-processing (including generation of animated images of flows and calculating forces and moments exerted on affected bodies). One of the tools, denoted OVERGRID, is a graphical user interface (GUI) that serves to visualize the grids and flow solutions and provides central access to many other tools. The GUI facilitates the generation of grids for a new flow-field configuration. Scripts that follow the grid generation process can then be constructed to mostly automate grid generation for similar configurations. CGT is designed for use in conjunction with a computer-aided-design program that provides the geometry description of the bodies, and a flow-solver program.

Chan, William M.↗

Neutrino transport in general relativistic neutron star merger simulations

Abstract Numerical simulations of neutron star–neutron star and neutron star–black hole binaries play an important role in our ability to model gravitational-wave and electromagnetic signals powered by these systems. These simulations have to take into account a wide range of physical processes including general relativity, magnetohydrodynamics, and neutrino radiation transport. The latter is particularly important in order to understand the properties of the matter ejected by many mergers, the optical/infrared signals powered by nuclear reactions in the ejecta, and the contribution of that ejecta to astrophysical nucleosynthesis. However, accurate evolutions of the neutrino transport equations that include all relevant physical processes remain beyond our current reach. In this review, I will discuss the current state of neutrino modeling in general relativistic simulations of neutron star mergers and of their post-merger remnants. I will focus on the three main types of algorithms used in simulations so far: leakage, moments, and Monte-Carlo scheme. I will review the advantages and limitations of each scheme, as well as the various neutrino–matter interactions that should be included in simulations. We will see that the quality of the treatment of neutrinos in merger simulations has greatly increased over the last decade, but also that many potentially important interactions remain difficult to take into account in simulations (pair annihilation, oscillations, inelastic scattering).

Foucart, Francois (ORCID:0000000346174738)↗

Relative advantages of thin-layer Navier-Stokes and interactive boundary-layer procedures

Numerical procedures for solving the thin-shear-layer Navier-Stokes equations and for the interaction of solutions to inviscid and boundary-layer equations are described and evaluated. To allow appraisal of the numerical and fluid dynamic abilities of the two schemes, they have been applied to one airfoil as a function of angle of attack at two slightly different Reynolds numbers. The NACA 0012 airfoil has been chosen because it allows comparison with measured lift, drag, and moment and with surface-pressure distributions. Calculations have been performed with algebraic eddy-viscosity formulations, and they include consideration of transition. The results are presented in a form that allows easy appraisal of the accuracy of both procedures and of the relative costs. The interactive procedure is computationally efficient but restrictive relative to the thin-layer Navier-Stokes procedure. The latter procedure does a better job of predicting drag than does the former. In both procedures, the location of transition is crucial for accurate or detailed computations, particularly at high angles of attack. When the upstream influence of pressure field through the shear layer is important, the thin-layer Navier-Stokes procedure has an edge over the interactive procedure.

Mehta, U.↗

A comparison of interactive boundary-layer and thin-layer Navier-Stokes procedures

Numerical procedures for solving the thin-shear-layer Navier-Stokes equations and for the interaction of solutions to inviscid and boundary-layer equations are described and evaluated. To allow appraisal of the numerical and fluid dynamic abilities of the two schemes, they have been applied to one airfoil as a function of angle of attack at two slightly different Reynolds numbers. The NACA 0012 airfoil has been chosen because it allows comparison with measured lift, drag, and moment and with surface-pressure distributions. Calculations have been performed with algebraic eddy-viscosity formulations, and they include consideration of transition. The results are presented in a form that allows easy appraisal of the accuracy of both procedures and of the relative costs. The interactive procedure is computationally efficient but restrictive relative to the thin-layer Navier-Stokes procedure. The latter procedure does a better job of predicting drag than does the former. In both procedures, the location of transition is crucial for accurate or detailed computations, particularly at high angles of attack. When the upstream influence of pressure field through the shear layer is important, the thin-layer Navier-Stokes procedure has an edge over the interactive procedure.

Mehta, U.↗

Fast Fermi and gradient drift acceleration of electrons at nearly perpendicular collisionless shocks

Electrons can be reflected and effectively energized at quasi-perpendicular shocks. At a curved shock this process is most prominent where the upstream magnetic field and the shock surface are tangent. A theoretical explanation of the underlying physical mechanism in terms of fast Fermi acceleration has been proposed in the literature. The theory uses properties of the de Hoffmann-Teller frame and assumes conservation of magnetic moment in a static, simplified shock profile. Here, the discussion is extended in order to clarify certain aspects of the reflection process and to pinpoint the physical mechanisms that are operative and dominant from the viewpoint of the normal incidence frame. By reducing the analysis to the essential physical content and solving the pertinent energy equation, the equivalence of fast Fermi and gradient drift acceleration is shown.

Krauss-Varban, D.↗

A mapping closure for turbulent scalar mixing using a time-evolving reference field

A general mapping-closure approach for modeling scalar mixing in homogeneous turbulence is developed. This approach is different from the previous methods in that the reference field also evolves according to the same equations as the physical scalar field. The use of a time-evolving Gaussian reference field results in a model that is similar to the mapping closure model of Pope (1991), which is based on the methodology of Chen et al. (1989). Both models yield identical relationships between the scalar variance and higher-order moments, which are in good agreement with heat conduction simulation data and can be consistent with any type of epsilon(phi) evolution. The present methodology can be extended to any reference field whose behavior is known. The possibility of a beta-pdf reference field is explored. The shortcomings of the mapping closure methods are discussed, and the limit at which the mapping becomes invalid is identified.

Girimaji, Sharath S.↗

A kinetic study of solar wind electrons in the transition region from collision dominated to collisionless flow

We have studied the evolution of the velocity distribution function of a test population of electrons in the solar corona and inner solar wind region, using a recently developed kinetic model. The model solves the time dependent, linear transport equation, with a Fokker-Planck collision operator to describe Coulomb collisions between the 'test population' and a thermal background of charged particles, using a finite differencing scheme. The model provides information on how non-Maxwellian features develop in the distribution function in the transition region from collision dominated to collisionless flow. By taking moments of the distribution the evolution of higher order moments, such as the heat flow, can be studied.

Lie-Svendsen, O.↗

Anomalous Hall effect in conical helimagnetic crystals

Spin-spiral texture can substantially change charge transport properties in helimagnets. Here we find the anomalous Hall effect (AHE) exhibiting the dramatic behavior with respect to chemical potential μ in conical magnetic structures. The direct conductivity demonstrates kinks, and the anomalous Hall current exhibits minima and maxima changing the sign. We analytically derive the expression for energy bands and eigenstates in the most general case. Because of the conical potential, the energy bands are split into two nonparabolic bands where the lower band can have one- or two-minima shapes in the k z direction ($\hat{z}$ is a direction of the spiral axis). We prove that the origin of the anomalous Hall effect is not topological and is due to the interplay between the asymmetry of energy bands in the x and z directions and spin restrictions in the phase space due to the conical potential. We also investigate the dependence of transport properties on cone half- angle θ, and find that the effects are most pronounced at θ=π/2 (a helical state). Electric current is calculated using the Boltzmann equation where the relaxation is caused by electron-acoustic phonon interaction. The transition probability is found to be a 2×2 matrix with nonvanishing off-diagonal elements indicating the strong interband transitions. The origin of interband transitions is because of the nature of the conical potential where conduction electron spins interact with localized magnetic moments. To verify the proposed theory, we calculate the temperature dependence of resistivity for MnSi crystals and find the discontinuity at the phase transition between conical and paramagnetic phases. The calculations are in the excellent agreement with the experimental data. In addition, we predict the discontinuity behavior for the anomalous Hall resistivity at the phase transition where the resistivity exhibits the abrupt change at T=T C , (a) to zero if the relativistic effects for the conduction electrons are small or (b) to a nonzero value if Rashba/Dresselhaus effects are taken into account.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Limit analysis for combined edge and pressure loading on a cylindrical shell.

Equations describing the stress field and velocity field occurring in a circular cylindrical shell at plastic collapse are derived corresponding to stress states lying on each face of a yield surface for a uniform shell of material obeying the Tresca yield condition. They are then applied to the case of a shell under combined axisymmetric loadings (moment, shear force, and axial force) at one end and uniform internal or external pressure on the lateral surface. For a sufficiently long shell, complete solutions are obtained for a fixed far end, and for a certain range of values of axial force and pressure, they are obtained for a free far end. All the solutions are represented by either closed form or by quadratures. It is shown that in many cases the radial velocity field is proportional to the shear force.

Ho, H. S.↗