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 685 records · Page 38

Relativistic two-moment neutrino transport

We implement a general relativistic neutrino transport scheme appropriate to problems in stellar collapse and neutron star formation. In this scheme the spectral energy and momentum balance equations are solved as a coupled set. We perform transport calculations on two frozen stellar background models neglecting the feedback of the neutrinos on the matter. We compare with schemes that ignore either relativity or momentum balance or both. The work performed by the neutrino stress on the matter is taken into account and has an important effect. In particular, if there is a hydrodynamical shock in the system, this work contribution boosts the energy deposition behind the shock. This may enhance the chances for a successful neutrino driven supernova explosion.

Cernohorsky, J.↗

Transient scattering from two-dimensional dielectric cylinders of arbitrary shape

In this work, the problem of transient scattering by arbitrary shaped two-dimensional dielectric cylinders is solved using the marching-on-in-time (MOT) technique. The dielectric problem is approached via the surface equivalence principle. A pair of coupled integral equations are derived by enforcing the continuity of the electric and magnetic fields which are solved by using the method of moments. Numerical results are presented for two cross sections, viz. a circle and a square, and compared with inverse discrete Fourier transform (IDFT) techniques. In each case, good agreement is obtained with the IDFT solution.

Vechinski, Douglas A.↗

Electrokinetic transport of heterogeneous particles in suspension

The Smoluchowski equation for electrophoresis predicts that the electrophoretic velocity of a particle is proportional to its zeta potential but not its size, shape, or orientation. Furthermore, the equation predicts that the rotation rate is identically zero. The Smoluchowski equation fails for heterogeneous particles (i.e., those with nonuniform zeta potentials). Recent theories and experiments show that particles with a dipole moment of zeta potential rotate into alignment with an externally applied electric field. For doublets (particles composed of two spheres) the rotation rate depends on (1) whether the spheres are rigidly rocked or freely rotating, and (2) the gap distance between the spheres. The relative configuration of two coagulated spheres is determined by the colloidal forces of the system. The goal of our research is to use measurements of electrophoretic rotation to determine the gap between two spheres of a colloidal doublet and also to determine whether or not the doublet is rigid.

Velegol, D.↗

Spectral Density of Laser Beam Scintillation in Wind Turbulence: Theory - Part 1

The temporal spectral density of the log-amplitude scintillation of a laser beam wave due to a spatially dependent vector-valued crosswind (deterministic as well as random) is evaluated. The path weighting functions for normalized spectral moments are derived, and offer a potential new technique for estimating the wind velocity profile. The Tatarskii-Klyatskin stochastic propagation equation for the Markov turbulence model is used with the solution approximated by the Rytov method. The Taylor 'frozen-in' hypothesis is assumed for the dependence of the refractive index on the wind velocity, and the Kolmogorov spectral density is used for the refractive index field.

Balakrishnan, A. V.↗

Neutron Star Equation of State in Light of GW190814

The observation of gravitational waves from an asymmetric binary opens the possibility for heavy neutron stars, but these pose challenges to models of the neutron star equation of state. Here, we construct heavy neutron stars by introducing nontrivial structure in the speed of sound sourced by deconfined QCD matter, which cannot be well recovered by spectral representations. Their moment of inertia, Love number, and quadrupole moment are very small, so a tenfold increase in sensitivity may be needed to test this possibility with gravitational waves, which is feasible with third generation detectors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Equilibrium expectations for non-Gaussian fluctuations near a QCD critical point

With the highly anticipated results from the Beam Energy Scan II program at RHIC being recently revealed, an understanding of particle-number fluctuations and their significance as a potential signature of a possible QCD critical point is crucial. Early works that embarked on this endeavor sought to estimate the fluctuations due to the presence of a critical point assuming they stay in equilibrium. From these results came the proposal to focus efforts on higher, non-Gaussian, moments of the event-by-event distributions, in particular of the number of protons. These non-Gaussian moments are especially sensitive to critical fluctuations, as their magnitudes are proportional to high powers of the critical correlation length. As the equation of state provides key input for hydrodynamical simulations of heavy-ion collisions, we estimate equilibrium fluctuations from the BEST equation of state (EoS) that includes critical features from the 3D Ising Model. In particular, the proton factorial cumulants and their dependence on non-universal mapping parameters is investigated within the BEST EoS. Furthermore, the correlation length, as a central quantity for the assessment of fluctuations in the vicinity of a critical point, is also calculated in a consistent manner with the scaling equation of state. An understanding of the equilibrium estimates of proton factorial cumulants will be useful for further comparison to estimates of out-of-equilibrium fluctuations in order to determine the magnitude of the observable fluctuations to be expected in heavyion collision experiments, in which the time spent near a critical point is short.

Karthein, Jamie M. [Massachusetts Institute of Tec↗

Deterministic errors in the Magellan orbit due to attitude control thruster activity

The pitch and yaw attitude control thrusters of the Magellan spacecraft generate uncoupled moments about the spacecraft center of mass, perturbing the spacecraft's orbit. A strategy to model these perturbing forces in the spacecraft equations of motion is presented. This strategy can model the thruster forces occurring in the nominal and contingency modes of attitude control operation. In the nominal mode, the thrusters fire only to balance the daily unloading of the momentum wheels. The contingency case will occur if a pitch or yaw axis momentum wheel fails, and the thrusters could fire up to 1200 times around the orbit to replace the failed wheel's momentum contribution. Deterministic errors are computed and shown to be under specified error requirements.

Engelhardt, Douglas B.↗

An Euler analysis of leading-edge vortex flows on a forebody-strake at supersonic speeds

The flowfield, surface pressure, and integrated forces and moments for a fighter-type forebody with sharp leading-edge strakes have been obtained by numerical solution of the Euler equations. The method is found to correctly predict the leading-edge vortices and embedded shocks which arise at higher angles of attack. Results are presented for the effects of crossflow grid density, artificial viscosity, angle of attack, streamwise station, and camber on the flow characteristics.

Rose, O. J.↗

Blade mistuning coupled with shaft flexibility effects in rotor aeroelasticity

The effect of bladed-disk polar dissymmetry, resulting from variations in mass from one blade to another, on aeroelastic stability boundaries for a fan stage is presented. In addition to both in-plane and out-of-plane deformations of the bladed-disk, bending of the supporting shaft in two planes is considered, and the resulting Coriolis forces and gyroscopic moments are included in the analysis. A quasi-steady aerodynamics approach is combined with the Lagrangian method to develop the governing equations of motion for the flexible bladed-disk-shaft assembly. Calculations are performed for an actual fan stage.

Khader, Naim↗

Jet-boundary corrections for reflection-plane models in rectangular wind tunnel

A detailed method for determining the jet-boundary corrections for reflection-plane models in rectangular wind tunnels is presented. The method includes the determination of the tunnel span local distribution and the derivation of equations for the corrections to the angle of attack, the lift and drag coefficients, and the pitching-, rolling-, yawing-, and hinge-moment coefficients. The principle effects of aerodynamic induction and of the boundary-induced curvature of the streamlines have been considered. An example is included to illustrate the method. Numerical values of the more important corrections for reflection-plane models in 7 by 10-foot closed wind tunnels are presented.

Swanson, Robert S↗

Turbulence modeling for separated flow

Two projects are described in this report. The first involves assessing turbulence models in separated flow. The second addresses the anomalous behavior of certain turbulence models in stagnation point flow. The primary motivation for developing turbulent transport models is to provide tools for computing non-equilibrium, or complex, turbulent flows. Simple flows can be analyzed using data correlations or algebraic eddy viscosities, but in more complicated flows such as a massively separated boundary layer, a more elaborate level of modeling is required. It is widely believed that at least a two-equation transport model is required in such cases. The transport equations determine the evolution of suitable velocity and time-scales of the turbulence. The present study included assessment of second-moment closures in several separated flows, including sharp edge separation; smooth wall, pressure driven separation; and unsteady vortex shedding. Flows with mean swirl are of interest for their role in enhancing mixing both by turbulent and mean motion. The swirl can have a stabilizing effect on the turbulence. An axi-symmetric extension to the INS-2D computer program was written adding the capability of computing swirling flow. High swirl can produce vortex breakdown on the centerline of the jet and it occurs in various combustors.

Durbin, Paul A.↗

Mass Anomalies on Ganymede

Radio Doppler data from two Ganymede encounters (G1 and G2) on the first two orbits in the Galileo mission have been analyzed previously for gravity information . For a satellite in hydrostatic equilibrium, its gravitational field can be modeled adequately by a truncated spherical harmonic series of degree two. However, a fourth degree field is required in order to fit the second Galileo flyby (G2). This need for a higher degree field strongly suggests that Ganymede s gravitational field is perturbed by a gravity anomaly near the G2 closest approach point (79.29 latitude, 123.68 west longitude). In fact, a plot of the Doppler residuals , after removal of the best-fit model for the zero degree term (GM) and the second degree moments (J2 and C22), suggests that if an anomaly exists, it is located downtrack of the closest approach point, closer to the equator.

Schubert, G.↗

The Cauchy Problem for Boltzmann Bi-linear Systems: The Mixing of Monatomic and Polyatomic Gases

Abstract From a unified vision of vector valued solutions in weighted Banach spaces, this paper establishes the existence and uniqueness for space homogeneous Boltzmann bi-linear systems with conservative collisional forms arising in complex gas dynamical structures. This broader vision is directly applied to dilute multi-component gas mixtures composed of both monatomic and polyatomic gases. Such models can be viewed as extensions of scalar Boltzmann binary elastic flows, as much as monatomic gas mixtures with disparate masses and single polyatomic gases, providing a unified approach for vector valued solutions in weighted Banach spaces. Novel aspects of this work include developing the extension of a general ODE theory in vector valued weighted Banach spaces, precise lower bounds for the collision frequency in terms of the weighted Banach norm, energy identities, angular or compact manifold averaging lemmas which provide coerciveness resulting into global in time stability, a new combinatorics estimate for p -binomial forms producing sharper estimates for the k -moments of bi-linear collisional forms. These techniques enable the Cauchy problem improvement that resolves the model with initial data corresponding to strictly positive and bounded initial vector valued mass and total energy, in addition to only a $$2^+$$ 2 + moment determined by the hard potential rates discrepancy, a result comparable in generality to the classical Cauchy theory of the scalar homogeneous Boltzmann equation.

Physics↗

Thermally driven mass flows in the convection zone of the sun

A formulation of the fluid dynamics of convective regions is developed which leads to an analytical description of the solar rotation, the Evershed flow, and the supergranulation. The starting point of the present formulation is the mixing length picture of convective equilibrium, but the earlier point mass model for convective molecules is replaced here by a model with both inertia and intrinsic moment of inertia. This extension introduces three rotational degrees of freedom into the dynamics of individual convective molecules, which enter into the dynamical equations for a mixing length fluid in the form of a separate vector field which we term the spin field. It is shown that for convective molecules having a spherically symmetric mass distribution, the spin field is proportional to the local vorticity.

Dijkhuis, G. C.↗

Magnetization and magnetic susceptibility of DyH3

The magnetization and differential magnetic susceptibility of powdered samples of DyH3 have been measured at 4.2 K in applied magnetic fields ranging to 9 Teslas. The differential magnetic susceptibility has also been studied in zero applied field as a function of temperature. The magnetization data are described by an equation of the form M = aB/(1 + bB + cB. The ratio a/b is a measure of the saturation magnetization and gives an effective moment of 5.12 Bohr magnetons per ion. The zero field susceptibility exhibits a maximum at T = 3.45 K, and an inflection point near 2.85 K. The susceptibility at 4.2 K has a 1/B squared dependence on the applied magnetic field for B approximately greater than 0.3 Teslas.

Flood, D. J.↗

Hypercube matrix computation task

The Hypercube Matrix Computation (Year 1986-1987) task investigated the applicability of a parallel computing architecture to the solution of large scale electromagnetic scattering problems. Two existing electromagnetic scattering codes were selected for conversion to the Mark III Hypercube concurrent computing environment. They were selected so that the underlying numerical algorithms utilized would be different thereby providing a more thorough evaluation of the appropriateness of the parallel environment for these types of problems. The first code was a frequency domain method of moments solution, NEC-2, developed at Lawrence Livermore National Laboratory. The second code was a time domain finite difference solution of Maxwell's equations to solve for the scattered fields. Once the codes were implemented on the hypercube and verified to obtain correct solutions by comparing the results with those from sequential runs, several measures were used to evaluate the performance of the two codes. First, a comparison was provided of the problem size possible on the hypercube with 128 megabytes of memory for a 32-node configuration with that available in a typical sequential user environment of 4 to 8 megabytes. Then, the performance of the codes was anlyzed for the computational speedup attained by the parallel architecture.

Calalo, R.↗

A new adaptive control approach for aerospace vehicles with parameter uncertainties

A new stochastic adaptive control structure is developed for the problem of combined parameter estimation and control of aerospace vehicles with changing parameters. Parameter uncertainties are modeled as first-order Gauss-Markov processes, and are introduced to the system dynamics through a small parameter. It is assumed that an accurate inertial measurement unit gives perfect measurements of the state variables. Since the stochastic system is assumed to be Gauss-Markov, the density function of the parameters given these measurements is conditionally Gaussian. Based on this conditionally Gaussian density, the problem of minimizing a quadratic cost over an infinite time horizon can be set up within the framework of stochastic optimal control theory. The optimal feedback control law is derived from a straightforward expansion of the Hamilton-Jacobi-Bellman equation, based on the LQG solution. The resulting nonlinear controller is applied to the pitch axis control of a space platform with uncertain moments of inertia and is shown to produce marked improvement over a fixed controller.

Hahn, Yungsun↗

The Effect of Systematic Error in Forced Oscillation Testing

One of the fundamental problems in flight dynamics is the formulation of aerodynamic forces and moments acting on an aircraft in arbitrary motion. Classically, conventional stability derivatives are used for the representation of aerodynamic loads in the aircraft equations of motion. However, for modern aircraft with highly nonlinear and unsteady aerodynamic characteristics undergoing maneuvers at high angle of attack and/or angular rates the conventional stability derivative model is no longer valid. Attempts to formulate aerodynamic model equations with unsteady terms are based on several different wind tunnel techniques: for example, captive, wind tunnel single degree-of-freedom, and wind tunnel free-flying techniques. One of the most common techniques is forced oscillation testing. However, the forced oscillation testing method does not address the systematic and systematic correlation errors from the test apparatus that cause inconsistencies in the measured oscillatory stability derivatives. The primary objective of this study is to identify the possible sources and magnitude of systematic error in representative dynamic test apparatuses. Sensitivities of the longitudinal stability derivatives to systematic errors are computed, using a high fidelity simulation of a forced oscillation test rig, and assessed using both Design of Experiments and Monte Carlo methods.

Williams, Brianne Y.↗