Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “First Order Second Moment”

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 55 records · Page 3

Carrying the mass flux terms exactly in the first and second moment equations of compressible turbulence

In compressible turbulence models, it is assumed that the Favre-mean velocities are suitable approximations to the Reynolds-mean velocities in order to close unknown terms. This neglects, in the mean momentum and energy equations, the contribution to the stress and work terms by the mean of the fluctuating Favre velocity, a quantity proportional to the turbulent mass flux. As the stress and work terms do not introduce any new unknown correlations requiring closure in either k-epsilon or Reynolds stress closures and because the exact form of the terms can, with little additional work, be carried there is no need to make any modeling assumptions. In the Reynolds stress equations the viscous terms appear naturally in Reynolds variables while the problem is posed in Favre variables. In the process of splitting the viscous terms into the viscous transport terms, carried in Favre variables, and the dissipation terms, carried in Reynolds variables, important contributions from the mass flux appear. The accurate accounting of these terms is important for any consistent near wall modeling and the retention of the mass flux terms is important in complex compressible turbulent flows.

Ristorcelli, J. R.↗

Optimal compensator structure for linear time-invariant plant with inaccessible states

The problem is considered of designing an optimal linear time-invariant dynamic compensator for the regulation of an n-th order linear time-invariant plant with m independent outputs. The initial plant state is characterized by its first and second moments, and the cost is usual quadratic infinite-time penalty on the state and control, averaged over the initial plant and compensator states. The compensator is based on a minimal-order Luenberger observer and consequently has fixed dimension n-m. Necessary and sufficient conditions are derived for optimality of the compensator gains. The optimal compensator is shown to be unique if the plant has a particular canonical form, and, in general, for any arbitrary plant, the class of all optimal compensators is precisely determined.

Blanvillain, P. J. P.↗

Coupling Between Waveguide-Fed Slot Arrays

Coupling between two waveguide-fed planar slot arrays has been investigated using full-wave analysis. The analysis employs the method-of-moments solution to the pertinent coupled integral equations for the aperture electric field of all slots. In order to compute coupling between two arrays, the input port of the first array is excited with a TE(sub 10) mode wave while the second one is match-terminated. After solving the moment method matrix equations, the aperture fields of all slots are obtained and thereby the TE(sub 10) mode wave received at the input port of the second array is determined. Coupling between two arrays is the ratio of the wave amplitude arriving in the second array port to the incident wave amplitude at the first array port. The coupling mechanism has been studied as a function of spacing between arrays in different directions, e.g. the electric field plane, the magnetic field plane, and the diagonal plane. Computed coupling values are presented for different array geometries. This work is novel since it provides a good understanding of coupling between waveguide-fed slot arrays as a function of spacing and orientation for different aperture distributions and array architectures. This serves as a useful tool for antenna design engineers and system engineers.

Rengarajan, Sembiam↗

A hybrid Monte Carlo-deterministic second moment method with efficient variance reduction

In this work, we present a hybrid method that combines Monte Carlo with deterministic finite element methods to solve a linear Boltzmann transport equation. Our hybrid method runs orders of magnitude faster than Monte Carlo, without sacrificing accuracy, for a proxy problem from radiative transfer that contains both optically-thick and optically-thin material. We believe that this is the first demonstration of a hybrid Second Moment Method in more than one spatial dimension, the first to consider more than one material, and the first to use variance reduction. Our variance reduction approach arises from an asymptotic analysis in which we show that the magnitude of the scattering source grows without bound. We transform the problem to compute the deviation of the radiation intensity from isotropy. The magnitude of the source in the transformed problem is bounded, and the quality of the hybrid method solution is dramatically improved by a substantial reduction in the variance.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modelling of Collision Induced Absorption Spectra Of H2-H2 Pairs for the Planetary Atmospheres Structure: The Second Overtone Band

The main objective of the proposal was to model the collision induced, second overtone band of gaseous hydrogen at low temperatures. The aim of this work is to assist planetary scientists in their investigation of planetary atmospheres, mainly those of Uranus and Neptune. The recently completed extended database of collision induced dipole moments of hydrogen pairs allowed us, for the first time, to obtain dipole moment matrix elements responsible for the roto-vibrational collision induced absorption spectra of H2-H2 in the second overtone band. Despite our numerous attempts to publish those data, the enormous volume of the database did not allow us to do this. Instead, we deposited the data on a www site. The final part of this work has been partially supported by NASA, Division for Planetary Atmospheres. In order to use our new data for modelling purpose, we first needed to test how well we can reproduce the existing experimental data from theory, when using our new input data. Two papers resulted from this work. The obtained agreement between theoretical results and the measurements appeared to be within 10-30%. The obviously poorer agreement than observed for the first H2 overtone, the fundamental, and the rototranslational bands can be attributed to the fact that dipole moments responsible for the second overtone are much weaker, therefore susceptible to larger numerical uncertainties. At the same time, the intensity of the second overtone band is much weaker and therefore it is much harder to be measured accurately in the laboratory. We need to point out that until now, no dependable model of the 2nd overtone band was available for modelling of the planetary atmospheres. The only one, often referred to in previous works on Uranian and Neptune's atmospheres, uses only one lineshape, with one (or two) parameter(s) deduced at the effective temperature of Uranus (by fitting the planetary observation). After that, the parameter(s) was(were) made temperature dependent according to some very simple relation. Summarizing, no reliable temperature-dependent model has been available yet. Our approach was a bit different from similar attempts done earlier, on account of the poorer agreement of theory with experiment. We needed to resort to some semi-empirical procedure. While we were in a favourable position to be able to rely on the physical input data, these, apparently, did not supply the most dependable predictions (simply because the results did not agree well enough with experimental data). On the other hand, the relative deviations between the theory and experiment were comparable at 77 and at 298 K. That fact indicated that theory is capable of predicting the temperature dependence of the absorption spectra well. We have thus chosen the "middle way". We have fitted the existing measurements with many 3- parameter lineshapes, in order to achieve the closest fit.

Borysow, Aleksandra↗

Anisotropic magnetism and electronic structure of trigonal EuAl 2 Ge 2 single crystals

Understanding the interplay between magnetic and electronic degrees of freedom is of profound recent interest in different Eu-based magnetic topological materials. Here, we studied the magnetic and electronic properties of the layered Zintl-phase compound EuAl 2 Ge 2 crystallizing in the trigonal CaAl 2 Si 2 -type structure. We report zero-field neutron diffraction, temperature T - and magnetic-field H-dependent magnetic susceptibility χ(T, H), isothermal magnetization M(T, H), heat capacity C p (T, H), and electrical resistivity ρ(T, H) measurements, together with T-dependent angle-resolved photoemission spectroscopy (ARPES) measurements complemented with first-principle calculations. EuAl 2 Ge 2 undergoes second-order A-type antiferromagnetic (AFM) ordering below T N = 27.5(5) K, with the Eu moments (Eu 2+ , S = 7/2) aligned ferromagnetically in the ab plane while these layers are stacked antiferromagnetically along the c axis. The critical fields at which all moments become parallel to the field are 37.5(5) and 52.5(5) kOe for H ∥ ab and H ∥ c, respectively. The H = 0 magnetic structure consists of trigonal AFM domains associated with ab-plane magnetic anisotropy and a field-induced reorientation of the Eu spins in the domains is also evident at T = 2 K below the critical field H c1 = 2.5(1) kOe. The ρ(T) measurements reveal metallic behavior transforming into a slight resistivity increase on cooling towards T N . A pronounced loss of spin-disorder scattering is observed below T N . The ARPES results show that EuAl 2 Ge 2 is metallic both above and below T N , and the Fermi surface is anisotropic with two hole pockets at the zone center and one small electron pocket at each M point. In the AFM phase, we directly observe folded bands in ARPES due to the doubling of the magnetic unit cell along the c axis with an enhancement of quasiparticle weight due to the complex change in the coupling between the magnetic moments and itinerant electrons on cooling below T N . The observed electronic structure is well reproduced by first-principle calculations, which also predict the presence of nontrivial electronic states near the Fermi level in the AFM phase with Z 2 topological numbers 1;(000)

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Sensitivity derivatives for advanced CFD algorithm and viscous modelling parameters via automatic differentiation

The computational technique of automatic differentiation (AD) is applied to a three-dimensional thin-layer Navier-Stokes multigrid flow solver to assess the feasibility and computational impact of obtaining exact sensitivity derivatives typical of those needed for sensitivity analyses. Calculations are performed for an ONERA M6 wing in transonic flow with both the Baldwin-Lomax and Johnson-King turbulence models. The wing lift, drag, and pitching moment coefficients are differentiated with respect to two different groups of input parameters. The first group consists of the second- and fourth-order damping coefficients of the computational algorithm, whereas the second group consists of two parameters in the viscous turbulent flow physics modelling. Results obtained via AD are compared, for both accuracy and computational efficiency with the results obtained with divided differences (DD). The AD results are accurate, extremely simple to obtain, and show significant computational advantage over those obtained by DD for some cases.

Green, Lawrence L.↗

Spectral Moments of Collision-Induced Absorption of CO2 Pairs: The Role of the Intermolecular Potential

In this paper we examine the role of the anisotropy of the intermolecular potential in the rototranslational collision-induced absorption of the CO2 pairs. Using newly developed formulas that include the effects of anisotropy of the potential to all orders, we calculate the two lowest spectral moments gamma(prime), and alpha(prime), for four different classes of C02 pair potentials and compare the results with the experimental values. We assumed only multipolar induction in the process of forming the induced dipole, with the second-order contributions included. Using a site-site LJ and a site-site semi-ab initio intermolecular potentials we were able to reproduce the experimental values of gamma(prime), and alpha(prime) moments over entire temperature range from 230 to 330 K. Also, the role of an electrostatic interaction between two C02 molecules and its impact on the spectral moments is thoroughly investigated. An isotropic core with a point quadrupole centered at each molecule is shown to be an inadequate representation of the C02-CO2 potential. Additionally, we show the results obtained with the first- and second-order perturbation theory to be more than twice too small.

Gruszka, Marcin↗

First-order antiferromagnetic transitions of SrMn 2 P 2 and CaMn 2 P 2 single crystals containing corrugated-honeycomb Mn sublattices

Significance With rare exceptions, an antiferromagnetic (AFM) transition in zero magnetic field is thermodynamically of second order where the thermal-average magnetic moments of the magnetic atoms (ordered moments) vary continuously on cooling through the AFM ordering temperature T N with no latent heat at the transition. Such materials include the AFM pnictides CaMn 2 As 2 , SrMn 2 As 2 , CaMn 2 Sb 2 , SrMn 2 Sb 2 , and CaMn 2 Bi 2 . Here we demonstrate that the closely related SrMn 2 P 2 and CaMn 2 P 2 insulators instead exhibit first-order AFM transitions at T N = 53 and 70 K, respectively, where the heat capacity exhibits a latent heat at T N . The mechanism causing these first-order transitions remains to be explained, but its understanding may lead to the development of novel magnetic materials of technological interest.

36 MATERIALS SCIENCE↗

Numerical solution of the Navier-Stokes equations with topography

A finite difference scheme for solving the equations of fluid motion in a generalized coordinate system has been constructed. The scheme conserves mass and all the first integral moments of the motion. The scheme also advectively 'almost conserves' second moments, in that the magnitude of implicit numerical smoothing is typically about an order smaller than explicit viscosity and diffusion. Calculations with the model support the theoretical conjecture that the difference scheme is stable whenever the analogous Cartesian scheme is stable. The scheme has been used to calculate dry atmospheric convection due to differential heating between top and bottom of mountainous terrain. The general small-scale characteristics of mountain up-slope winds have been simulated. In addition, the results have demonstrated the crucial role played by the eddy diffusivities and the environmental stability, in determining both the quantitative and the qualitative features of the circulation.

Gal-Chen, T.↗

Analyses of convective flows and spatial gradients in energetic ion observations

Based on the measured anisotropy of three-dimensional energetic ion distributions obtained from the Medium Energy Particle Instrument experiment of the ISEE 1 spacecraft, bulk flow, epsilon x B convection, and pressure gradients in magnetospheric plasmas are measured with two techniques. In the first technique, using a Compton-Getting transform into the appropriate rest frame in an iterative fashion and converging to the unique velocity that removes the first-order anisotropy from the flux angular distribution, the actual measured angular distribution, and spectrum at each angle, are employed. The second technique uses a first moment of the distribution function in which the shift in velocity space of the distribution function is obtained at two phase space density values. The convection contribution to the velocity shift is assumed to be perpendicular to the pressure gradient contribution if an MHD generator condition is assumed, rendering the two shifts separable, and a gradient scale and a lower limit to the convection electric field are derived in one calculation.

Mitchell, D. G.↗

High-Order Space-Time Methods for Conservation Laws

Current high-order methods such as discontinuous Galerkin and/or flux reconstruction can provide effective discretization for the spatial derivatives. Together with a time discretization, such methods result in either too small a time step size in the case of an explicit scheme or a very large system in the case of an implicit one. To tackle these problems, two new high-order space-time schemes for conservation laws are introduced: the first is explicit and the second, implicit. The explicit method here, also called the moment scheme, achieves a Courant-Friedrichs-Lewy (CFL) condition of 1 for the case of one-spatial dimension regardless of the degree of the polynomial approximation. (For standard explicit methods, if the spatial approximation is of degree p, then the time step sizes are typically proportional to 1/p(exp 2)). Fourier analyses for the one and two-dimensional cases are carried out. The property of super accuracy (or super convergence) is discussed. The implicit method is a simplified but optimal version of the discontinuous Galerkin scheme applied to time. It reduces to a collocation implicit Runge-Kutta (RK) method for ordinary differential equations (ODE) called Radau IIA. The explicit and implicit schemes are closely related since they employ the same intermediate time levels, and the former can serve as a key building block in an iterative procedure for the latter. A limiting technique for the piecewise linear scheme is also discussed. The technique can suppress oscillations near a discontinuity while preserving accuracy near extrema. Preliminary numerical results are shown

Huynh, H. T.↗

Dissipative magnetohydrodynamics for nonresistive relativistic plasmas: An implicit second-order flux-conservative formulation with stiff relaxation

Based on a 14-moment closure for nonresistive (general-) relativistic viscous plasmas, here we describe a new numerical scheme that is able to handle all first-order dissipative effects (heat conduction, bulk and shear viscosities), as well the anisotropies induced by the presence of magnetic fields. The latter is parametrized in terms of a thermal gyrofrequency or, equivalently, a thermal Larmor radius and allows to correctly capture the thermal Hall effect. By solving an extended Israel-Stewart-like system for the dissipative quantities that enforces algebraic constraints via stiff-relaxation, we are able to cast all first-order dissipative terms in flux-divergence form. This allows us to apply traditional high-resolution shock capturing methods to the equations, making the system suitable for the numerical study of highly turbulent flows. We present several numerical tests to assess the robustness of our numerical scheme in flat spacetime. The 14-moment closure can seamlessly interpolate between the highly collisional limit found in neutron star mergers, and the highly anisotropic limit of relativistic Braginskii magnetohydrodynamics appropriate for weakly collisional plasmas in black-hole accretion problems. We believe that this new formulation and numerical scheme will be useful for a broad class of relativistic magnetized flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Vacancy-tuned magnetism in LaMn x Sb 2

The layered AMP 2 (A = alkali-earth or rare-earth atom, M = transition metal, P = Sb, Bi) compounds are widely studied for their rich magnetism and electronic structure topology. We provide a detailed characterization of the magnetic and transport properties of LaMn x Sb 2 , an understudied member of the AMP 2 family. LaMn x Sb 2 forms with intrinsic Mn vacancies, and we demonstrate that by varying the starting ratio of La, Mn, and Sb, we can synthetically control the Mn occupancy and produce single crystals with x = 0.74 – 0.97. Magnetization and transport measurements indicate LaMn x Sb 2 has a rich temperature-composition (T–x) magnetic phase diagram with physical properties strongly influenced by the Mn occupancy. LaMn x Sb 2 orders antiferromagnetically at T 1 = 130–180 K, where T 1 increases with x. Below T 1 , the T–x phase diagram is complicated. At high x, there is a second transition T 2 that decreases in temperature as x is lowered, vanishing below x ≤ 0.85. A third, first-order, transition T 3 is detected at x ≈ 0.92, and the transition temperature increases as x is lowered, crossing above T 2 near x ≈ 0.9. On moving below x < 0.79, here we find the crystal structure changes from the P 4/nmm arrangement to an I$\bar{4}$2m structure with partially ordered Mn vacancies. The change in crystal structure results in the sudden appearance of two new low-temperature phases and a crossover between regimes of negative and positive magnetoresistance when x ≤ 0.78. Finally, we provide powder neutron diffraction for x = 0.93, and find that the high-x compositions first adopt a G-type antiferromagnetic structure with the Mn moments aligned within the ab plane, which is followed upon further cooling by a second transition to a different, noncollinear structure where the moments are rotated within the basal plane. Our results demonstrate that LaMn x Sb 2 is a highly tunable material with six unique magnetically ordered phases, depending on T and x.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A dissipation-induced superradiant transition in a strontium cavity-QED system

Driven-dissipative many-body systems are ubiquitous in nature and a fundamental resource for quantum technologies. However, they are also complex and hard to model because they cannot be described by the standard tools in equilibrium statistical mechanics. Probing nonequilibrium critical phenomena in pristine setups can illuminate fresh perspectives on these systems. Here, we use an ensemble of cold 88 Sr atoms coupled to a driven high-finesse cavity to study the cooperative resonance fluorescence (CRF) model, a classic driven-dissipative model describing coherently driven dipoles superradiantly emitting light. We observe its nonequilibrium phase diagram characterized by a second-order phase transition. Below a critical drive strength, the atoms quickly reach the so-called superradiant steady state featuring a macroscopic dipole moment; above the critical point, the atoms undergo persistent Rabi-like oscillations. At longer times, spontaneous emission transforms the second-order transition into a discontinuous first-order transition. Our observations pave the way for harnessing robust entangled states and exploring boundary time crystals in driven-dissipative systems.

Science & Technology - Other Topics↗

Measurements of multiplicity fluctuations of identified hadrons in inelastic proton–proton interactions at the CERN Super Proton Synchrotron

Measurements of multiplicity fluctuations of identified hadrons produced in inelastic p+p interactions at 31, 40, 80, and 158 $\text {Ge}\text {V}/c$ beam momentum are presented. Three different measures of multiplicity fluctuations are used: the scaled variance $\omega $ and strongly intensive measures $\Sigma $ and $\Delta $. These fluctuation measures involve second and first moments of joint multiplicity distributions. Data analysis is preformed using the Identity method which corrects for incomplete particle identification. Strongly intensive quantities are calculated in order to allow for a direct comparison to corresponding results on nucleus–nucleus collisions. The results for different hadron types are shown as a function of collision energy. A comparison with predictions of string-resonance Monte-Carlo models: Epos, Smash and Venus, is also presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Near surface analysis

This is a study to assist in the understanding of earth near surface structure. Higher order moments are used to detect the density distribution as well as to seek patterns found in geological structures. It is shown how higher order moments at points outside a mass structure are determined as well as how to recover the mass distribution from the higher order moments. It is interesting to note that the first moment at a point P outside the mass structure, V(sub O)(P), is the entire mass and the second moment, V(sub 1)(P), is the potential at P due to the mass structure. Usually only the mass and the potential function are used to determine the density distribution in a body. An infinite function sequence (V(sub n)(P))(sub n=0)(sup infinity) is required to uniquely determine the density distribution.

Johnson, Gordon G.↗

Dark Energy Survey Year 6 Results: Point-Spread Function Modeling

We present the point-spread function (PSF) modeling for weak lensing shear measurement using the full six years of the Dark Energy Survey (DES Y6) data. We review the PSF estimation procedure using the PIFF (PSFs In the Full FOV) software package and describe the key improvements made to PIFF and modeling diagnostics since the DES year three (Y3) analysis: (i) use of external Gaia and infrared photometry catalogs to ensure higher purity of the stellar sample used for model fitting, (ii) addition of color-dependent PSF modeling, the first for any weak lensing analysis, and (iii) inclusion of model diagnostics inspecting fourth-order moments, which can bias weak lensing measurements to a similar degree as second-order modeling errors. Through a comprehensive set of diagnostic tests, we demonstrate the improved accuracy of the Y6 models evident in significantly smaller systematic errors than those of the Y3 analysis, in which all g band data were excluded due to insufficiently accurate PSF models. For the Y6 weak lensing analysis, we include g band photometry data in addition to the riz bands, providing a fourth band for photometric redshift estimation. Looking forward to the next generation of wide-field surveys, we describe several ongoing improvements to PIFF, which will be the default PSF modeling software for weak lensing analyses for the Vera C. Rubin Observatory’s Legacy Survey of Space and Time.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗