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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Level 2 Milestone: Develop and Incorporate Novel Code Verification of Fundamental Equations in Gemma and Set Up Appropriate Tests

For computational physics simulations, code verification plays a major role in establishing the credibility of the results by assessing the correctness of the implementation of the underlying numerical methods. In computational electromagnetics, surface integral equations, such as the method-of-moments implementations of the electric-, magnetic-, and combinedfield integral equations, are frequently used to solve Maxwell’s equations on the surfaces of electromagnetic scatterers. These electromagnetic surface integral equations yield many code-verification challenges due to the various sources of numerical error and their possible interactions. In this report, we provide approaches to separately measure the numerical errors arising from these different error sources. We demonstrate the effectiveness of these approaches in Gemma.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2012

OVERSMART Reporting Tool for Flow Computations Over Large Grid Systems

Structured grid solvers such as NASA's OVERFLOW compressible Navier-Stokes flow solver can generate large data files that contain convergence histories for flow equation residuals, turbulence model equation residuals, component forces and moments, and component relative motion dynamics variables. Most of today's large-scale problems can extend to hundreds of grids, and over 100 million grid points. However, due to the lack of efficient tools, only a small fraction of information contained in these files is analyzed. OVERSMART (OVERFLOW Solution Monitoring And Reporting Tool) provides a comprehensive report of solution convergence of flow computations over large, complex grid systems. It produces a one-page executive summary of the behavior of flow equation residuals, turbulence model equation residuals, and component forces and moments. Under the automatic option, a matrix of commonly viewed plots such as residual histograms, composite residuals, sub-iteration bar graphs, and component forces and moments is automatically generated. Specific plots required by the user can also be prescribed via a command file or a graphical user interface. Output is directed to the user s computer screen and/or to an html file for archival purposes. The current implementation has been targeted for the OVERFLOW flow solver, which is used to obtain a flow solution on structured overset grids. The OVERSMART framework allows easy extension to other flow solvers.

Kao, David L.↗
DOE OSTI2024

Orbital angular momentum at small $x$ revisited

We revisit the problem of the small Bjorken-x asymptotics of the quark and gluon orbital angular momentum (OAM) distributions in the proton utilizing the revised small-x helicity evolution derived recently in [1]. We relate the quark and gluon OAM distributions at small x to the polarized dipole amplitudes and their (first) impact-parameter moments. To obtain the OAM distributions, we derive novel small-x evolution equations for the impact-parameter moments of the polarized dipole amplitudes in the double-logarithmic approximation (summing powers of α s ln 2 (1/x) with α s the strong coupling constant). We solve these evolution equations numerically and extract the leading large-N c , small-x asymptotics of the quark and gluon OAM distributions, which we determine to be ${L}_{q+\overline{q}}\left(x,{Q}^2\right)\sim {L}_G\left(x,{Q}^2\right)\sim \Delta \Sigma \left(x,{Q}^2\right)\sim \Delta G\left(x,{Q}^2\right)\sim {\left(\frac{1}{2}\right)}^{3.66\sqrt{\frac{\alpha_s{N}_c}{2_π}}}$in agreement with [2] within the precision of our numerical evaluation. (Here N c is the number of quark colors.) We also investigate the ratios of the quark and gluon OAM distributions to their helicity distribution counterparts in the small-x region.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Gravothermal evolution of dark matter halos with differential elastic scattering

Here, we study gravothermal evolution of dark matter halos in the presence of differential self-scattering that has strong velocity and angular dependencies. We design controlled N-body simulations to model Rutherford and Møller scatterings in the halo, and follow its evolution in both core-expansion and -collapse phases. The simulations show the commonly-used transfer cross section underestimates the effects of dark matter self-interactions, but the viscosity cross section provides an accurate approximation for modeling angular-dependent dark matter scattering. We investigate thermodynamic properties of the halo, and find that the three moments of the Boltzmann equation under the fluid approximation are satisfied. We further propose a constant effective cross section, which integrates over the halo's characteristic velocity dispersion with weighting kernels motivated by kinetic theory of heat conduction. The effective cross section provides a good approximation to differential self-scattering for most of the halo evolution. It indicates that we can map astrophysical constraints on a constant self-interacting cross section to an SIDM model with velocity- and angular-dependent scatterings.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Simulating dust grain-radiation coupling on a moving mesh

ABSTRACT We present a model for the interaction between dust and radiation fields in the radiation hydrodynamic code arepo-rt, which solves the moment-based radiative transfer equations on an unstructured moving mesh. Dust is directly treated using live simulation particles, each of which represent a population of grains that are coupled to hydrodynamic motion through a drag force. We introduce methods to calculate radiation pressure on and photon absorption by dust grains. By including a direct treatment of dust, we are able to calculate dust opacities and update radiation fields self-consistently based on the local dust distribution. This hybrid scheme coupling dust particles to an unstructured mesh for radiation is validated using several test problems with known analytic solutions, including dust driven via spherically symmetric flux from a constant luminosity source and photon absorption from radiation incident on a thin layer of dust. Our methods are compatible with the multifrequency scheme in arepo-rt, which treats UV, optical photons as single scattered and IR photons as multi scattered. At IR wavelengths, we model heating of and thermal emission from dust. Dust and gas are not assumed to be in local thermodynamic equilibrium but transfer energy through collisional exchange. We estimate dust temperatures by balancing these dust-radiation and dust-gas energy exchange rates. This framework for coupling dust and radiation can be applied in future radiation hydrodynamic simulations of galaxy formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Implementation of Monte Carlo Transport in the General Relativistic SpEC Code

Neutrino transport and neutrino-matter interactions are known to play an important role in the evolution of neutron star mergers and of their post-merger remnants. Neutrinos cool remnants, drive post-merger winds, and deposit energy in the low-density polar regions where relativistic jets may eventually form. Neutrinos also modify the composition of the ejected material, impacting the outcome of nucleosynthesis in merger outflows and the properties of the optical/infrared transients that they power (kilonovae). So far, merger simulations have largely relied on approximate treatments of the neutrinos (leakage, moments) that simplify the equations of radiation transport in a way that makes simulations more affordable but also introduces unquantifiable errors in the results. To improve on these methods, we recently published a first simulation of neutron star mergers using a low-cost Monte Carlo algorithm for neutrino radiation transport. Our transport code limits costs in optically thick regions by placing a hard ceiling on the value of the absorption opacity of the fluid, yet all approximations made within the code are designed to vanish in the limit of infinite numerical resolution. We provide here an in-depth description of this algorithm, of its implementation in the SpEC merger code, and of the expected impact of our approximations in optically thick regions. We argue that the last is a subdominant source of error at the accuracy reached by current simulations and for the interactions currently included in our code. We also provide tests of the most important features of this code.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1974

Analysis of the dynamics of a nutating body

The equations for the displacement, velocity, and acceleration of a point in a nutating body are developed. These are used to derive equations for the inertial moment developed by a nutating body of arbitrary shape. Calculations made for a previously designed nutating plate transmission indicate that that device is severely speed limited because of the very high magnitude inertial moment.

Anderson, W. J.↗
NASA NTRS1973

An unsteady wake model for a hingeless rotor

A simple nonsteady wake model derived from the unsteady moment of the momentum equation for zero advance ratio is correlated with cyclic pitch frequency response tests conducted with a small hingeless rotor model. Two and three or more bladed rotor analyses are presented.

Crews, S. T.↗
NASA NTRS1974

Estimation of material properties of a nonlinearly elastic bar

A method of determining certain characteristic flexural rigidities and elastic properties of nonlinearly elastic materials is presented. An estimation method utilizing perturbation methods and a least squares fitting technique is used to solve the nonlinear differential equation derived from the moment curvature relation, subject to boundary values representing deflections of the bar at discrete points. Deflection data from numerical simulations of a nonlinearly elastic, prismatic bar are used to demonstrate the estimation method. Numerical experiments relating the accuracy of the identification to the number and accuracy of the boundary values are presented. Conclusions based on the numerical experiments are included.

Baker, B. E.↗
NASA NTRS1975

Measurement of the inertial properties of the Helios F-1 spacecraft

A gravity pendulum method of measuring lateral moments of inertia of large structures with an error of less than 1% is outlined. The method is based on the fact that in a physical pendulum with a knife-edge support the distance from the axis of rotation to the system center of gravity determines the minimal period of oscillation and is equal to the system centroidal radius of gyration. The method is applied to results of a test procedure in which the Helios F-1 spacecraft was placed in a roll fixture with crossed flexure pivots as elastic constraints and system oscillation measurements were made with each of a set of added moment-of-inertia increments. Equations of motion are derived with allowance for the effect of the finite pivot radius and an error analysis is carried out to find the criterion for maximum accuracy in determining the square of the centroidal radius of gyration. The test procedure allows all measurements to be made with the specimen in upright position.

Gayman, W. H.↗
NASA NTRS1978

On the accelerations of the moon and sun, the constant of gravitation, and the origin of mountains

On the assumption of a constant value for the gravitational constant G, the lunar and solar tidal couples and the rate of change of angular velocity of the earth are recalculated on the basis of the improved values for the apparent accelerations of the moon and sun arrived at by Muller and Stephenson (1977). These new values are then used to derive general dynamical equations with a changing moment of inertia and changing G connecting the apparent accelerations to the lunar and solar couples. The effect of a changing G is considered, and a value of (dG/dt)/G = -3 times 10 to the -11th power per yr can be reconciled only weakly with the improved accelerations, especially for the case of linear friction. For quadratic friction the situation is more favorable. The new accelerations, however, are best reconciled with a nonchanging G, and certainly cannot be reconciled with a value of (dG/dt)/G = -6 times 10 to the -11th power per yr.

Lyttleton, R. A.↗
NASA NTRS1981

Electron temperature anisotropy in the polar wind

The steady state flow of a fully ionized H(+)-O(+)-electron plasma along geomagnetic field lines in the high-latitude topside ionosphere is investigated theoretically, with emphasis on the electron temperature anisotropy and heat flow in the polar wind. The 13-moment system of transport equations developed by Schunk (1975, 1977), which contains a continuity, momentum, internal energy, stress tensor and heat flow equation for each species, is employed to represent the electrons, with a simplified set of transport equations used for the ions, and the transport equations are solved at altitudes from 1500 to 12,000 km for both subsonic and supersonic H(+) outflows. For subsonic H(+) outflows, the electron gas is found to remain collision dominated to high altitudes. For supersonic H(+) outflows, the electron gas is also found to be collision dominated at altitudes below 2500 km, with an anisotropy in electron temperature distribution developing such that the temperature of the perpendicular electrons is greater than that of the parallel electrons at higher altitudes. In addition, the magnitude of downward electron heat flux at 1500 km is shown to have a dramatic effect on the individual parallel and perpendicular electron temperatures.

Schunk, R. W.↗
NASA NTRS1982

Proton temperature anisotropy in the polar wind

The steady state flow of a fully ionized H(+)-O(+) electron plasma along geomagnetic field lines in the high-latitude topside ionosphere is studied. The theoretical formulation is based on a 13-moment system of transport equations, and allows for different species temperatures parallel and perpendicular to the geomagnetic field and nonclassical heat flows. For subsonic and supersonic flows, an appreciable H(+) temperature anisotropy occurs at all altitudes above 1500 km, and tends to be regulated at high altitudes. The direction of the temperature anisotropy is related to the direction of the H(+) heat flow to some extent; for supersonic flow an upward flow of heat from the lower ionosphere is required, while for subsonic flow solutions can be obtained with a downward H(+) heat flow. For subsonic flow, H(+)-H(+) collisions have an important effect on the H(+) stress and heat flow balance at all altitudes between 1500 and 12,000 km.

Schunk, R. W.↗
NASA NTRS1986

Thermohydrodynamic analysis for laminar lubricating films

A Galerkin-type analysis to include thermal effects in laminar lubricating films was performed. The lubricant properties were assumed constant except for a temperature-dependent Newtonian viscosity. The cross-film temperature profile is established by collocation at the film boundaries and two interior Lobatto points. The interior temperatures are determined by requiring that the zeroth and first moment of the energy equation be satisfied across the film. The fluidity is forced to conform to a third--degree polynomial appropriate to the Lobatto-point temperatures. Preliminary indications are that the use of just two such sampling points enables satisfactory prediction of bearing performance even in the presence of substantial viscosity variation.

Elrod, H. G.↗
NASA NTRS1991

Krylov model reduction algorithm for undamped structural dynamics systems

Krylov vectors furnish an efficient basis for eigenvalue analysis and model reduction of structural dynamics systems. The reduced-order model obtained by the present Krylov model-reduction algorithm for an undamped structural-dynamics system is found to match low-frequency moments. The transformed system equation in Krylov coordinates reflects the structure of a tandem system.

Craig, Roy R., Jr.↗
NASA NTRS1992

Anisotropic ion heating and parallel O(+) acceleration in regions of rapid E x B convection

A numerical solution to the 20-moment set of transport equations has been found in order to study subauroral ionospheric outflows during periods of enhanced perpendicular ion drifts. The numerical model solves the time-dependent O(+) density, momentum, and both the parallel and perpendicular energy and heat flow equations in the 200-6000 km altitude range. Assuming perpendicular drifts of 3 km/s relative to the neutral atmosphere, we have found that anisotropic heating of O(+) (a result of ion-neutral collisions) leads to a temperature anisotropy, with perpendicular temperatures exceeding 8000 K and parallel temperatures greater than 5000 K (near 200 km altitude). Above approximately 2000 km, transport processes dominate the effects of collisions and wavelike oscillations in O(+) velocity, temperature and heat flux were noted.

Korosmezey, A.↗