Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density Laws”

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 91 records · Page 5

Preliminary scaling laws for plasma current, ion kinetic temperature, and plasma number density in the NASA Lewis bumpy torus plasma

Parametric variation of independent variables which may affect the characteristics of bumpy torus plasma have identified those which have a significant effect on the plasma current, ion kinetic temperature, and plasma number density, and those which do not. Empirical power law correlations of the plasma current, and the ion kinetic temperature and number density were obtained as functions of potential applied to the midplane electrode rings, the background neutral gas pressure, and the magnetic field strength. Additional parameters studied included the type of gas, the polarity of the midplane electrode rings, the mode of plasma operation, and the method of measuring the plasma number density. No significant departures from the scaling laws appear to occur at the highest ion kinetic temperatures or number densities obtained to date.

Roth, J. R.↗

Stress due to electric charge density distribution in a dielectric slab

The spatial distribution of electric field due to an imposed electric charge density profile in an infinite slab of dielectric material is derived analytically by integrating Gauss’s law. Various charge density distributions are considered, including exponential and power-law forms. Here, the Maxwell stress tensor is used to compute a notional static stress in the material due to the charge density and its electric field. Characteristics of the electric field and stress distributions are computed for example cases in polyethylene, showing that field magnitudes exceeding the dielectric strength would be required in order to achieve a stress exceeding the ultimate tensile strength.

42 ENGINEERING↗

Preliminary scaling laws for plasma current, ion kinetic temperature, and plasma number density in the NASA Lewis Bumpy Torus plasma

Parametric variation of independent variables which may affect the characteristics of the NASA Lewis Bumpy Torus plasma have identified those which have a significant effect on the plasma current, ion kinetic temperature, and plasma number density, and those which do not. Empirical power-law correlations of the plasma current, and the ion kinetic temperature and number density were obtained as functions of the potential applied to the midplane electrode rings, the background neutral gas pressure, and the magnetic field strength. Additional parameters studied include the type of gas, the polarity of the midplane electrode rings (and hence the direction of the radial electric field), the mode of plasma operation, and the method of measuring the plasma number density. No significant departures from the scaling laws appear to occur at the highest ion kinetic temperatures or number densities obtained to date.

Roth, J. R.↗

What Sets the Star Formation Rate of Molecular Clouds? The Density Distribution as a Fingerprint of Compression and Expansion Rates

We use a suite of 3D simulations of star-forming molecular clouds, with and without stellar feedback, magnetic fields, and driven turbulence, to study the compression and expansion rates of the gas as functions of density. We show that, around the mean density, supersonic turbulence promotes rough equilibrium between the amounts of compressing and expanding gas, consistent with continuous gas cycling between high- and low-density states. We find that the inclusion of protostellar jets produces rapidly expanding and compressing low-density gas. We find that the gas mass flux peaks at the transition between the lognormal and power-law forms of the density probability distribution function (PDF). This is consistent with the transition density tracking the post-shock density, which promotes an enhancement of mass at this density (i.e., shock compression and filament formation). At high densities, the gas dynamics are dominated by self-gravity: the compression rate in all of our runs matches the rate of the run with only gravity, suggesting that processes other than self-gravity have little effect at these densities. The net gas mass flux becomes constant at a density below the sink formation threshold, where it equals the star formation rate. The density at which the net gas mass flux equals the star formation rate is one order of magnitude lower than our sink threshold density, corresponds to the formation of the second power-law tail in the density PDF, and sets the overall star formation rates of these simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

The dynamics and stability of radiatively driven gas clouds. I - Plane-parallel slabs

A combination of numerical and analytical techniques has been used to investigate the dynamics and stability of optically thin plane-parallel radiatively driven slabs of gas confined by the thermal gas pressure of a high-temperature low-density medium. Scaling laws allow the individual model 'clouds' to be characterized by a single free parameter, chi, a normalized column density which measures the strength of the acceleration due to radiation pressure relative to that due to thermal gas pressure. It is found that these clouds are stable and coherently accelerated only when chi is small. In this regime a simple slab model is constructed which accurately reproduces the more complex gasdynamic results. The low-chi clouds are marginally able to reach the high velocities seen in the atmospheres of quasi-stellar objects, but only if their motion is subsonic with respect to the external confining medium. This implies either that the medium is extremely hot and tenuous or that it is moving outward with the clouds.

Haas, M. R.↗

Gapless spin liquid and pair density wave of the Hubbard model on three-leg triangular cylinders

We study the ground state properties of the Hubbard model on three-leg triangular cylinders using large-scale density-matrix renormalization group simulations. At half-filling, we identify an intermediate gapless spin liquid phase, which has one gapless spin mode and algebraic spin–spin correlations but exponential decay scalar chiral–chiral correlations, between a metallic phase at weak coupling and Mott insulating dimer phase at strong interaction. Upon light doping the gapless spin liquid, the system exhibits power-law charge-density-wave (CDW) correlations but short-range single-particle, spin–spin, and chiral–chiral correlations. Similar to CDW correlations, the superconducting correlations also decay in power-law but oscillate in sign as a function of distance, which is consistent with the striped pair-density wave. When further doping the gapless spin liquid phase or doping the dimer order phase, another phase takes over, which has similar CDW correlations but all other correlations decay exponentially.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Relativistically transparent magnetic filaments: scaling laws, initial results and prospects for strong-field QED studies

Abstract Relativistic transparency enables volumetric laser interaction with overdense plasmas and direct laser acceleration of electrons to relativistic velocities. The dense electron current generates a magnetic filament with field strength of the order of the laser amplitude (>10 5 T). The magnetic filament traps the electrons radially, enabling efficient acceleration and conversion of laser energy into MeV photons by electron oscillations in the filament. The use of microstructured targets stabilizes the hosing instabilities associated with relativistically transparent interactions, resulting in robust and repeatable production of this phenomenon. Analytical scaling laws are derived to describe the radiated photon spectrum and energy from the magnetic filament phenomenon in terms of the laser intensity, focal radius, pulse duration, and the plasma density. These scaling laws are compared to 3D particle-in-cell (PIC) simulations, demonstrating agreement over two regimes of focal radius. Preliminary experiments to study this phenomenon at moderate intensity ( a 0 ∼ 30) were performed on the Texas Petawatt Laser. Experimental signatures of the magnetic filament phenomenon are observed in the electron and photon spectra recorded in a subset of these experiments that is consistent with the experimental design, analytical scaling and 3D PIC simulations. Implications for future experimental campaigns are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Magnetohydrodynamic model of the Crab Nebula and its radiation

The interaction of the Crab pulsar with its surrounding nebula is modeled as a steady state, spherically symmetric relativistic wind which is confined by the outer layers of the supernova remnant. The pulsar wind is assumed to contain only positrons and electrons, and is terminated by a strong shock at a distance of 10 arcseconds from the pulsar, at the outer edge of the underluminous region. Down-stream of the shock, the flow decelerates and increases its pressure in order to match the boundary conditions imposed by outer nebula. Consistent flow solutions can be found if the magnetic luminosity of the pulsar wind is about 0.3 percent of its particle luminosity. The immediate postshock plasma is modeled as a power law distribution whose density and pressure moments satisfy the strong shock Rankine-Hugoniot relations. Using the downstream flow solution, the synchrotron continuum is calculated. The parameters of the pulsar wind and assumed power law spectral index are adjusted until the calculated total synchrotron power and the spectral distribution of the nebular continuum from optical to gamma-rays, are in reasonable accord with observations.

Coroniti, F. V.↗

The radial gradients and collisional properties of solar wind electrons

The plasma instrument on Mariner 10 carried out measurements of electron density and temperature in the interplanetary medium between heliocentric distances of 0.85 and 0.45 AU. Due to the stable coronal configuration and low solar activity during the period of observation, the radial variations of these quantities could be obtained. The power-law exponent of the core temperature was measured to be -0.3 + or - 0.04, and the halo temperature was found to be almost independent of heliocentric distance. The exponent of the power law for the density variation was 2.5 + or - 0.2 and the extrapolated value at 1 AU was consistent with measured values during the same period. Calculations of the core electron self-collision time, and the core-halo equipartition time were made as a function of radial distance. These measurements indicate a macroscale picture of a Coulomb-collisional core and a collisionless isothermal halo. Extrapolating back to the sun, core and halo temperatures become equal at a radial distance of approx. 2-15 radii.

Gilvie, K. W.↗

Thermodynamic properties of “near-perfect” gas

When applied thermodynamics requires simple representations of thermodynamic quantities, one often finds quantities such as pressure and the thermal energy density expressed as a monomial linear product of density and temperature power laws. This is a simple generalization of the product of density and temperature used for the perfect gas that admits a broader range of thermodynamic behavior into analytic fluid calculations or provides an analytic form that is readily fit to tabulated equation-of-state data for an arbitrary material. This paper reviews the thermodynamic properties of this generalized perfect-gas model, treating it as a class of “near-perfect-gas” models that are defined and unified here, for the first time, through a shared Helmholtz free energy. This generalization from perfect to near-perfect-gas models preserves some perfect-gas properties (e.g., constant adiabatic exponents) but not others (e.g., constant specific heats) and is constrained by thermodynamic consistency. It is important to be aware of the properties of this class of near-perfect-gas models when using them as surrogates for real nonperfect materials that may have properties that cannot be captured by near-perfect-gas models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Generation of coronal electric currents due to convective motions on the photosphere

Generation of electric currents in a magnetized plasma overlying a dense convective layer is studied, assuming that the magnetic field perturbation is small and satisfies the force-free equation. Currents are produced by rotational motions on the boundary in the case of a uniform equilibrium field. In a simple two-dimensional bipolar configuration, however, both irrotational and incompressible motions give rise to currents, and the current density has a peak at the magnetic neutral line. Scaling laws for the current density as well as for the stored magnetic energy are derived, and the possibility of heating the solar corona through the dissipation of coronal currents generated in this way is discussed.

Sakurai, T.↗

Orbital resonances and planetary formation sites

A cascaded resonance structure where planetesimal growth was accelerated at 2:1 interior and 1:2 exterior resonances, with an early-formed Jupiter producing runaway growth of planetary embryos, is hypothesized in a solar system formation model. The planetary embryos produce their own resonances, and these in turn lead to additional embryos in a process that successively propagates inwardly and outwardly to generate a resonant configuration of embryos. The early presence of Jupiter would in this way have imposed a harmonic structure on the accumulating planetesimal swarm. The positions of the planetary embryos can be moved into a degree of agreement with most of the present planetary positions which is comparable to that given by the Titius-Bode law, for the case of an accretion disk whose surface density obeys a power law of index -1.2.

Torbett, M.↗

Identification of crystal plasticity model parameters by multi-objective optimization integrating microstructural evolution and mechanical data

Crystal plasticity models evolve a polycrystalline yield surface using meso-scale descriptions of deformation mechanisms. The activation of deformation mechanisms is governed by crystallography and a set of model parameters, which are typically calibrated through the fitting of mechanical data such as stress–strain curves and elastic lattice strains. Microstructural data such as phase fractions and texture evolution are used for verifying crystal plasticity parameters. In this study, we use a multi-objective genetic algorithm to identify hardening parameters from flow stress curves with an option to incorporate texture into the optimization approach. Robust, generalized objective functions are developed and used to identify sets of parameters pertaining to dislocation density-based hardening laws in visco-plastic and elasto-plastic self-consistent (VPSC and EPSC) homogenization models. First, the parameters are identified for pure Nb directly from texture using an objective function based on generalized spherical harmonics. Since texture evolution is driven by the relative contribution of active slip systems, the parameters governing the evolution of slip resistance ratios can be recovered from fitting discrete textures at a series of strains. Next, a comprehensive set of load reversal data for dual phase (DP) 780 steel is used to fit a hardening law and a back-stress law in EPSC. Finally, parameters pertaining to a complex hardening law for the evolution of slip and twinning in pure α-Ti are identified. Remarkably, using texture as an objective in combination with stress–strain objectives constrains the model of Ti to fully reproduce not only stress–strain and texture evolution but also hierarchical twinning measurements as a function of initial grain size and texture. Furthermore, given an appropriate model fit to representative experimental texture evolution, underlying twin volume fractions contributing to texture evolution can be predicted.

42 ENGINEERING↗

Experimental verification of a crystal plasticity-based simulation framework for predicting microstructure and geometric shape changes: Application to bending and Taylor impact testing of Zr

This paper is concerned with experimental verification of a recently developed multi-scale simulation framework for plastic deformation of metallic materials from quasi-static to impact deformation conditions. The framework is a visco-plastic self-consistent (VPSC) polycrystalline model embedded in an implicit finite element method (FE-VPSC) to provide a microstructure-sensitive constitutive response at each material point. Each material point of the FEM model is a polycrystalline aggregate with crystallographic deformation mechanisms operating at the single crystal scale with their evolving activity based on a dislocation density-based hardening law and texture. Four beams and three cylinders machined in different orientations from a textured plate of high-purity zirconium are tested quasi-statically in 4-point bending and at speeds of 100 m/s, 170 m/s and 243 m/s during Taylor impact tests, respectively. The variation in dimensional changes resulting from different sample orientations in the plate with respect to loading directions is measured for each sample. Moreover, texture and twinning characterization is performed using electron backscattered diffraction (EBSD). The deformation processes and underlying evolution of microstructure are successfully simulated using the FE-VPSC framework. In doing so, the model parameters are optimized and validated across a broad range of strain rates and temperatures. Simulation results in terms of geometrical changes and microstructural evolution are compared with the experimental measurements. The model predicts anisotropic material flow resulting from the hard-to-deform crystallographic directions, the development of gradients in texture and twinning through the geometries, tension–compression asymmetry, as well as the extent of plasticity under impact.

42 ENGINEERING↗

Investigation of crossed-twin structure formation in magnesium and magnesium alloys

In this work, the effect of alloying addition on the propensity for twin-twin interactions to transform into crossed-twin structures in magnesium alloys is investigated. A full-field elasto-viscoplastic fast Fourier transform (EVP-FFT) framework combined with a discrete twin model and dislocation density-based hardening law for slip strengths is used to calculate the micromechanical fields in the crystalline matrix around the interacting twins. AZ31 and MgLi alloys are selected along with pure Mg to study the influence of plastic anisotropy in connection with alloying elements. These alloys were selected since their plastic anisotropy measure, which is defined as the ratio between the critical resolved shear stress for pyramidal $\langle$c +a$\rangle$ and basal $\langle$a$\rangle$ slip modes, spanned a wide range. To quantify the role of twin thicknesses, we probe a range of impinging twin thicknesses while fixing the recipient twin thickness. The analysis reveals that: (i) the local driving stress for crossed twin structure formation generated from the interaction of the two twins is lower in a low plastically anisotropic material, like a MgLi alloy, than a high plastically anisotropic material like pure Mg and (ii) the critical impinging twin thickness needed to form the crossed twin structure in pure Mg, AZ31 and MgLi alloys is ~0.5, ~0.75 and ~1.5 times the recipient twin thickness. We propose a relationship between the tendency for crossed twin structure formation and the experimentally observed higher ductility in MgLi alloys compared to pure Mg. Finally, one key implication of the findings is that crossed- twin structure formation can be hindered and the ductility of magnesium alloy thereby improved by properly choosing alloying elements that lower the slip strength for pyramidal $\langle$c + a$\rangle$ slip.

36 MATERIALS SCIENCE↗

Black hole to breakout: 3D GRMHD simulations of collapsar jets reveal a wide range of transients

We present a suite of the first 3D GRMHD collapsar simulations, which extend from the self-consistent jet launching by an accreting Kerr black hole (BH) to the breakout from the star. We identify three types of outflows, depending on the angular momentum, l, of the collapsing material and the magnetic field, B, on the BH horizon: (i) subrelativistic outflow (low l and high B), (ii) stationary accretion shock instability (SASI; high l and low B), (iii) relativistic jets (high l and high B). In the absence of jets, free-fall of the stellar envelope provides a good estimate for the BH accretion rate. Jets can substantially suppress the accretion rate, and their duration can be limited by the magnetization profile in the star. We find that progenitors with large (steep) inner density power-law indices (≳ 2), face extreme challenges as gamma-ray burst (GRB) progenitors due to excessive luminosity, global time evolution in the light curve throughout the burst and short breakout times, inconsistent with observations. Our results suggest that the wide variety of observed explosion appearances (supernova/supernova + GRB/low-luminosity GRBs) and the characteristics of the emitting relativistic outflows (luminosity and duration) can be naturally explained by the differences in the progenitor structure. Our simulations reveal several important jet features: (i) strong magnetic dissipation inside the star, resulting in weakly magnetized jets by breakout that may have significant photospheric emission and (ii) spontaneous emergence of tilted accretion disc-jet flows, even in the absence of any tilt in the progenitor.

79 ASTRONOMY AND ASTROPHYSICS↗

The impact of the $\text{WHIM}$ on the $\text{IGM}$ thermal state determined from the low- z Lyman $\alpha$ forest

At z ≲ 1, shock heating caused by large-scale velocity flows and possibly violent feedback from galaxy formation, converts a significant fraction of the cool gas (T ~ 10 4 K) in the intergalactic medium (IGM) into warm–hot phase (WHIM) with T > 10 5 K, resulting in a significant deviation from the previously tight power-law IGM temperature–density relationship, T = T 0 (P/$\bar{p}$) γ-1 ⁠. This study explores the impact of the WHIM on measurements of the low-z IGM thermal state, [T 0 , γ], based on the b–NH 1 distribution of the Ly α forest. Exploiting a machine learning-enabled simulation-based inference method trained on Nyx hydrodynamical simulations, we demonstrate that [T 0 , γ] can still be reliably measured from the b–NH 1 distribution at z = 0.1, notwithstanding the substantial WHIM in the IGM. To investigate the effects of different feedback, we apply this inference methodology to mock spectra derived from the IllustrisTNG and Illustris simulations at z = 0.1. The results suggest that the underlying [T 0 , γ] of both simulations can be recovered with biases as low as |Δlog(T 0 /K)| ≲ 0.05 dex, |Δγ| ≲ 0.1, smaller than the precision of a typical measurement. Given the large differences in the volume-weighted WHIM fractions between the three simulations (Illustris 38 percent, IllustrisTNG 10 percent, and Nyx 4 per cent), we conclude that the b–N H1 distribution is not sensitive to the WHIM under realistic conditions. Finally, we investigate the physical properties of the detectable Ly α absorbers, and discover that although their T and Δ distributions remain mostly unaffected by feedback, they are correlated with the photoionization rate used in the simulation.

79 ASTRONOMY AND ASTROPHYSICS↗

Numerical Analysis of Radiation Effects on Fiber Optic Sensors

Optical fiber sensors (OFS) are a potential candidate for monitoring physical parameters in nuclear environments. However, under an irradiation field the optical response of the OFS is modified via three primary mechanisms: (i) radiation-induced attenuation (RIA), (ii) radiation-induced emission (RIE), and (iii) radiation-induced compaction (RIC). For resonance-based sensors, RIC plays a significant role in modifying their performance characteristics. In this paper, we numerically investigate independently the effects of RIC and RIA on three types of OFS widely considered for radiation environments: fiber Bragg grating (FBG), long-period grating (LPG), and Fabry-Perot (F-P) sensors. In our RIC modeling, experimentally calculated refractive index (RI) changes due to low-dose radiation are extrapolated using a power law to calculate density changes at high doses. The changes in RI and length are subsequently calculated using the Lorentz–Lorenz relation and an established empirical equation, respectively. The effects of both the change in the RI and length contraction on OFS are modeled for both low and high doses using FIMMWAVE, a commercially available vectorial mode solver. An in-depth understanding of how radiation affects OFS may reveal various potential OFS applications in several types of radiation environments, such as nuclear reactors or in space.

42 ENGINEERING↗