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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Magnetic fields of the outer planets

It is difficult to imagine a group of planetary dynamos more diverse than those visited by the Pioneer and Voyager spacecraft. The magnetic field of Jupiter is large in magnitude and has a dipole axis within 10 deg of its rotation axis, comfortably consistent with the paleomagnetic history of the geodynamo. Saturn's remarkable (zonal harmonic) magnetic field has an axis of symmetry that is indistinguishable from its rotation axis (mush less than 1 deg angular separation); it is also highly antisymmetric with respect to the equator plane. According to one hypothesis, the spin symmetry may arise from the differential rotation of an electrically conducting and stably stratified layer above the dynamo. The magnetic fields of Uranus and Neptune are very much alike, and equally unlike those of the other known magnetized planets. These two planets are characterized by a large dipole tilts (59 deg and 47 deg, respectively) and quadrupole moments (Schmidt-normalized quadrupole/dipole ratio approximately equal 1.0). These properties may be characteristic of dynamo generation in the relatively poorly conducting 'ice' interiors of Uranus and Neptune. Characteristics of these planetary magnetic fields are illustrated using contour maps of the field on the planet's surface and discussed in the context of planetary interiors and dynamo generation.

Connerney, J. E. P.↗
NASA NTRS1995

The second-moment climatology of the GATE rain rate data

The first part of this paper presents the description of the GARP (Global Atmospheric Research Program) Atlantic Tropical Experiment (GATE) 1 rain-rate data and its two-dimensional spectral and correlation characteristics, which has made it possible to accomplish the following: to show the concentration of a significant power along the frequency axis in the spatiotemporal spectra; to detect a diurnal cycle (which has a range of variation of about 3.4-5.4 mm/n) as one of the sources of bias in the rain statistics of satellite data; to study the distinction between the north-south and east-west transport of spatial rain-rate field and character of its anisotropy; to evaluate the scales of the distinction between second-moment estimates associated with ground and satellite samples; and to determine the appropriate spatial and temporal scales of simple linear stochastic models fitted to averaged rain-rate fields. The second part of this paper is devoted to an analysis of the diffusion of the rain rate by establishing a relationship between the parameters of the multivariate autoregressive model and the coefficients of a diffusion equation. This analysis led to the use of rain data to estimate the rain advection velocity as well as other coefficients of the diffusion equation of the corresponding field. The results obtained can be used for comparison with corresponding estimates of other sources of data (satellite, Tropical Oceans Global Atmosphere Coupled Ocean - Atmosphere Response Experiment (TOGA, COARE) or simulated by physical models), for generating multiple samples of any size, for solving the inverse problems of some of the hydrodynamic equations, and in some other areas of rain data analysis and modeling.

Polyak, Ilya↗
DOE OSTI2022

Simple new methods for deducing lifetimes in recoil distance Doppler-shift measurements

In this work, new approaches for lifetime determination using data from recoil distance Doppler-shift experiments are presented based on the fundamental properties of the functions describing the time evolution of the population of excited nuclear states. To some extent, one of them represents a contraction of the well-known Differential decay-curve method (DDCM) by using the most reliable data point [the maximum of the $n_i(t)$ function describing the population of level $\textit{i}$ in time] and a purely numerical procedure avoiding any fitting of decay curves. The combination with the standard DDCM analysis is promising for improving the reliability and the precision of the results for the lifetimes obtained. The novel part of the approach consists of using a chain of equations at the consecutive maxima of the ni(t) functions, which allow us to precisely determine the ratio of the lifetimes of two consecutive levels and, in the case where one of these lifetimes is known, to determine the unknown one. In addition, a simple integral derivation of the lifetime is presented involving the peak areas measured at different distances, and an application of the first moments (expectation values and centroids in time) of the $n_i(t)$ functions for determining lifetimes is also demonstrated to be useful.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

A tale of two polarization paradoxes: The diamagnetic polarization paradox

An accurate calculation of the total polarization charge density in a plasma is essential for a self-consistent determination of the electric field. Yet, the “diamagnetic polarization paradox” refers to the fact that there is a paradoxical factor of 1/2 difference between the pressure-driven “diamagnetic polarization” density calculated using real space drift theory vs action-angle space guiding center and gyrokinetic theory that has not been explained before. Here, we show that both results can be made consistent with one another. Half of the diamagnetic polarization is due to the transformation from the guiding center density to the real space density. The other half is due to the fact that, within the drift kinetic ordering assumptions, the guiding center density should be expressed as the gyroaverage of the density in the limit of vanishing Larmor radius. A comprehensive review of polarization is presented to complete the derivation, and then we derive results that are required for the polarization calculations to agree. Expressions for the diamagnetic polarization density are given that are accurate to first order in amplitude and all orders in gyroradius within the gyrokinetic theory for a constant magnetic field. Applications to Maxwell–Boltzmann particle distribution functions (PDFs), including anisotropic temperature, are presented. Local invariants, like total energy and toroidal momentum, do not generate net polarization effects; the electric and thermodynamic polarizations must precisely cancel. In contrast, anisotropic dependence on the magnetic moment generates a net polarization proportional to the temperature anisotropy. Finally, when sources are present, the equilibrium PDF is approximately the ratio of two orbit averages.

adiabatic theory↗
NASA NTRS1978

Electromagnetic fields radiated from a lightning return stroke - Application of an exact solution to Maxwell's equations

A solution is presented for the electromagnetic fields radiated by an arbitrarily oriented current filament over a conducting ground plane in the case where the current propagates along the filament at the speed of light, and this solution is interpreted in terms of radiation from lightning return strokes. The solution is exact in the fullest sense; no mathematical approximations are made, and the governing differential equations and boundary conditions are satisfied. The solution has the additional attribute of being specified in closed form in terms of elementary functions. This solution is discussed from the point of view of deducing lightning current wave forms from measurements of the electromagnetic fields and understanding the effects of channel tortuosity on the radiated fields. In addition, it is compared with two approximate solutions, the traditional moment approximation and the Fraunhofer approximation, and a set of criteria describing their applicability are presented and interpreted.

D M Le Vine↗
NASA NTRS1991

Properties of planetary fluids at high pressures and temperatures

Observational data obtained by the Voyager space probes to the giant planets Jupiter, Saturn, Uranus, and Neptune have provided valuable information, which is used to refine the picture of the nature of the interiors of these planets. Major results from the Voyager missions include observations of substantial magnetic fields and improved models of internal density distributions. The goal is to obtain equations of state and electrical conductivity data for planetary gases (H2 and He) and the ices (H2O, CH4, and NH3, and their mixtures), which are considered to be the major constituents of the giant planets. These data are needed to test theoretical data bases used to construct models of the chemical composition of planetary interiors, models which are consistent with observables such as mass, diameter, gravitational moments, rotation rate, and magnetic field. The 100 GPa (1 Mbar) pressures and several 1000 K temperatures in the giant planets can be achieved in the lab by the shock compression of liquid specimens. Results are briefly examined.

Nellis, W. J.↗
NASA NTRS1992

Extending a transonic small disturbance code to treat swept vertical surfaces

A flexible-swept vertical surface capability has been developed and implemented within the CAP-TSD transonic small disturbance (TSD) code. The new capability required a modification to the TSD equation and a grid transformation for swept vertical surfaces. Modifications to the vertical surface boundary conditions allow it to be treated as a flexible surface. The new capability extends the range of problems which the code can treat. In order to assess the accuracy of the modifications, calculations were performed for a rectangular T-tail configuration and an AGARD T-tail configuration. Unsteady forces and moments are presented for the rectangular T-tail oscillating in yaw for a range of reduced frequencies. Comparisons are presented with linear theory and experiment. Steady and unsteady surface pressures are presented for the AGARD T-tail along with generalized aerodynamic forces. Comparisons are made with linear theory. The comparisons demonstrate the accuracy of the vertical surface modifications.

Gibbons, Michael D.↗
NASA NTRS2003

Multibody Parachute Flight Simulations for Planetary Entry Trajectories Using "Equilibrium Points"

A method has been developed to reduce numerical stiffness and computer CPU requirements of high fidelity multibody flight simulations involving parachutes for planetary entry trajectories. Typical parachute entry configurations consist of entry bodies suspended from a parachute, connected by flexible lines. To accurately calculate line forces and moments, the simulations need to keep track of the point where the flexible lines meet (confluence point). In previous multibody parachute flight simulations, the confluence point has been modeled as a point mass. Using a point mass for the confluence point tends to make the simulation numerically stiff, because its mass is typically much less that than the main rigid body masses. One solution for stiff differential equations is to use a very small integration time step. However, this results in large computer CPU requirements. In the method described in the paper, the need for using a mass as the confluence point has been eliminated. Instead, the confluence point is modeled using an "equilibrium point". This point is calculated at every integration step as the point at which sum of all line forces is zero (static equilibrium). The use of this "equilibrium point" has the advantage of both reducing the numerical stiffness of the simulations, and eliminating the dynamical equations associated with vibration of a lumped mass on a high-tension string.

Raiszadeh, Ben↗
NASA NTRS2008

Numerical Determination of Critical Conditions for Thermal Ignition

The determination of ignition or thermal explosion in an oxidizing porous body of material, as described by a dimensionless reaction-diffusion equation of the form .tu = .2u + .e-1/u over the bounded region O, is critically reexamined from a modern perspective using numerical methodologies. First, the classic stationary model is revisited to establish the proper reference frame for the steady-state solution space, and it is demonstrated how the resulting nonlinear two-point boundary value problem can be reexpressed as an initial value problem for a system of first-order differential equations, which may be readily solved using standard algorithms. Then, the numerical procedure is implemented and thoroughly validated against previous computational results based on sophisticated path-following techniques. Next, the transient nonstationary model is attacked, and the full nonlinear form of the reaction-diffusion equation, including a generalized convective boundary condition, is discretized and expressed as a system of linear algebraic equations. The numerical methodology is implemented as a computer algorithm, and validation computations are carried out as a prelude to a broad-ranging evaluation of the assembly problem and identification of the watershed critical initial temperature conditions for thermal ignition. This numerical methodology is then used as the basis for studying the relationship between the shape of the critical initial temperature distribution and the corresponding spatial moments of its energy content integral and an attempt to forge a fundamental conjecture governing this relation. Finally, the effects of dynamic boundary conditions on the classic storage problem are investigated and the groundwork is laid for the development of an approximate solution methodology based on adaptation of the standard stationary model.

Luo, W.↗
NASA NTRS2010

Analysis of Sting Balance Calibration Data Using Optimized Regression Models

Calibration data of a wind tunnel sting balance was processed using a candidate math model search algorithm that recommends an optimized regression model for the data analysis. During the calibration the normal force and the moment at the balance moment center were selected as independent calibration variables. The sting balance itself had two moment gages. Therefore, after analyzing the connection between calibration loads and gage outputs, it was decided to choose the difference and the sum of the gage outputs as the two responses that best describe the behavior of the balance. The math model search algorithm was applied to these two responses. An optimized regression model was obtained for each response. Classical strain gage balance load transformations and the equations of the deflection of a cantilever beam under load are used to show that the search algorithm s two optimized regression models are supported by a theoretical analysis of the relationship between the applied calibration loads and the measured gage outputs. The analysis of the sting balance calibration data set is a rare example of a situation when terms of a regression model of a balance can directly be derived from first principles of physics. In addition, it is interesting to note that the search algorithm recommended the correct regression model term combinations using only a set of statistical quality metrics that were applied to the experimental data during the algorithm s term selection process.

Ulbrich, N.↗
DOE OSTI2020

thornado-transport: Anderson- and GPU-accelerated nonlinear solvers for neutrino-matter coupling

Algorithms for neutrino-matter coupling in core-collapse supernovae (CCSNe) are investigated in the context of a spectral two-moment model, which is discretized in space with the discontinuous Galerkin method, integrated in time with implicit-explicit (IMEX) methods, and implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado). The model considers electron neutrinos and antineutrinos and tabulated opacities from Bruenn (1985), which includes neutrino-electron scattering and pair processes. The nonlinear system arising from implicit time discretization of the equations governing neutrino-matter coupling is iterated to convergence using Anderson-accelerated fixed-point methods, which avoid formation of Jacobians and inversion of dense linear systems. Numerical experiments show that, for a given tolerance, a nested iteration scheme which aims to reduce opacity evaluations can lower the computational cost. Our initial port to GPUs, using both OpenMP and OpenACC, shows an overall speedup of up to ~ 100× when compared to results using a single CPU core. These results indicate that the algorithms implemented in thornado are well-suited to GPU acceleration.

Laiu, Paul↗
DOE OSTI2023

Physics-based adaptivity of a spectral method for the Vlasov–Poisson equations based on the asymmetrically-weighted Hermite expansion in velocity space

We propose a spectral method for the 1D-1V Vlasov–Poisson system where the discretization in velocity space is based on asymmetrically-weighted Hermite functions, dynamically adapted via a scaling α and shifting u of the velocity variable. Specifically, at each time instant an adaptivity criterion selects new values of α and u based on the numerical solution of the discrete Vlasov–Poisson system obtained at that time step. Once the new values of the Hermite parameters α and u are fixed, the Hermite expansion is updated and the discrete system is further evolved for the next time step. The procedure is applied iteratively over the desired temporal interval. The key aspects of the adaptive algorithm are: the map between approximation spaces associated with different values of the Hermite parameters that preserves total mass, momentum and energy; and the adaptivity criterion to update α and u based on physics considerations relating the Hermite parameters to the average velocity and temperature of each plasma species. For the discretization of the spatial coordinate, we rely on Fourier functions and use the implicit midpoint rule for time stepping. The resulting numerical method possesses intrinsically the property of fluid-kinetic coupling, where the low-order terms of the expansion are akin to the fluid moments of a macroscopic description of the plasma, while kinetic physics is retained by adding more spectral terms. Moreover, the scheme features conservation of total mass, momentum and energy associated in the discrete, for periodic boundary conditions. A set of numerical experiments confirms that the adaptive method outperforms the non-adaptive one in terms of accuracy and stability of the numerical solution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Strong ultrafast demagnetization due to the intraband transitions

Abstract Demagnetization in ferromagnetic transition metals driven by a femtosecond laser pulse is a fundamental problem in solid state physics, and its understanding is essential to the development of spintronic devices.Ab initiocalculation of time-dependent magnetic moment in the velocity gauge so far has not been successful in reproducing the large amount of demagnetization observed in experiments. In this work, we propose a method to incorporate intraband transitions within the velocity gauge through a convective derivative in the crystal momentum space. Our results for transition-element bulk crystals (bcc Fe, hcp Co and fcc Ni) based on the time-dependent quantum Liouville equation show a dramatic enhancement in the amount of demagnetization after the inclusion of an intraband term, in agreement with experiments. We also find that the effect of intraband transitions on each ferromagnetic material is distinctly different because of their band structure and spin property differences. Our finding has a far-reaching impact on understanding of ultrafast demagnetization.

Physics↗
NASA NTRS1974

Application of a parameter identification technique to a hingeless helicopter rotor

A mathematical model of a gyro-controlled, three-bladed hingeless helicopter rotor was developed and parameters of the model were estimated using a parameter identification technique. The flapping and feathering degrees of freedom of the blades were modeled. The equations of the model contain time-varying, periodic coefficients due to the forward speed of the rotor. A digital simulation of the analytical model was compared with wind-tunnel measurements to establish the validity of the model. Comparisons of steady-state and transient solutions of the analytical model with the tunnel measurements gave reasonably good matching of gyro angle but less satisfactory matching of hub moment measurements. Further improvements were obtained by use of a parameter identification technique to adjust as many as 10 parameters of the analytical model. The sensitivity of the blade response to small changes in the parameters was also calculated.

Kanning, G.↗
NASA NTRS1978

Theory of flux anisotropies in a guiding center plasma

The one particle distribution function f on the scale of the bounce motion of particles in a magnetic field B is considered. The Vlasov equation is expanded through O(epsilon) in the adiabatic parameter which is the ratio of particle gyroradius to scale length of the magnetic field. Because f is directly proportional to particle flux differential in kinetic energy and solid angle, f is in principle measurable in space experiments, and the analysis is tailored to be explicitly applicable to space problems. To O(1), f is gyrotropic; its first velocity moment is (if non-vanishing) parallel to B, and hence macroscopic parallel flow is included in this term. The O(epsilon) contribution is non-gyrotropic and macroscopic flow parallel to B plus additional parallel flow results from these terms. The degree of non-gyrotropy and the amount of cross-field macroscopic flow depend on the perpendicular component of the electric field, on curvature and shear in the magnetic field, and on the spatial gradient, pitch angle derivative, and speed derivative of the lowest order distribution function.

Birmingham, T. J.↗
NASA NTRS1952

Matrix Method of Determining the Longitudinal-Stability Coefficients and Frequency Response of an Aircraft from Transient Flight Data

A matrix method is presented for determining the longitudinal-stability coefficients and frequency response of an aircraft from arbitrary maneuvers. The method is devised so that it can be applied to time-history measurements of combinations of such simple quantities as angle of attack, pitching velocity, load factor, elevator angle, and hinge moment to obtain the over-all coefficients. Although the method has been devised primarily for the evaluation of stability coefficients which are of primary interest in most aircraft loads and stability studies, it can be used also, with a simple additional computation, to determine the frequency-response characteristics. The entire procedure can be applied or extended to other problems which can be expressed by linear differential equations.

Donegan, James J↗
NASA NTRS2006

Direct Numerical Simulation of Transitional Multicomponent-Species Gaseous and Multicomponent-Liquid Drop-Laden Mixing

A model of multicomponent-liquid (MC-liquid) drop evaporation in a three-dimensional mixing layer is here exercised at larger Reynolds numbers than in a previous study, and transitional states are obtained. The gas phase is followed in an Eulerian frame and the multitude of drops is described in a Lagrangian frame. Complete coupling between phases is included with source terms in the gas conservation equations accounting for the drop/flow interaction in terms of drop drag, drop heating and species evaporation. The liquid composition, initially specified as a single-Gamma (SG) probability distribution function (PDF) depending on the molar mass is allowed to evolve into a linear combination of two SGPDFs, called the double-Gamma PDF (DGPDF). The compositions of liquid and vapor emanating from the drops are calculated through four moments of the DGPDFs, which are drop-specific and location-specific, respectively. The mixing layer is initially excited to promote the double pairing of its four initial spanwise vortices into an ultimate vortex in which small scales proliferate. Simulations are performed for four liquids of different compositions and the effect of the initial mass loading and initial free-stream gas temperature are explored. For reference, Simulations are also performed for gaseous multicomponent mixing layers for which the effect of Reynolds number is investigated. The results encompass examination of the global layer characteristics, flow visualizations and homogeneous-plane statistics at transition. Comparisons are performed with previous pre-transitional MC-liquid simulations and with transitional single-component (SC) liquid studies. It is found that MCC flows at transition, the classical energy cascade is of similar strength, but that the smallest scales contain orders of magnitude less energy than SC flows, which is confirmed by the larger viscous dissipation in the former case. Contrasting to pre-transitional MC flows, the vorticity and drop organization depend on the initial gas temperature, this being due to the drop/turbulence coupling. The vapor-composition mean molar mass and standard deviation distributions strongly correlate with the initial liquid-composition PDF; such a correlation only exists for the magnitude of the mean but not for that of the standard deviation. Unlike in pre-transitional situations, regions of large composition standard deviation no longer necessarily coincide with regions of large mean molar mass. The kinetic energy, rotational and composition characteristics, and dissipation are liquid specific and the variation among liquids is amplified with increasing free-stream gas temperature. Eulerian and Lagrangian statistics of gas-phase quantities show that the different. Observation framework may affect the perception of the flow characteristics. The gas composition, of which the first four moments are calculated, is shown to be close to, but distinct from a SGPDF. The PDF of the scalar dissipation rate is calculated for drop-laden layers and is shown to depart more significantly from the typically assumed Gaussian in gaseous flows than experimentally measured gaseous scalar dissipation rates, this being attributed to the increased heterogeneity due to drop/flow interactions.

tubulent multicompenent-species↗