Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Accounting”

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

Nucleation kinetics model for primary crystallization in Al–Y–Fe metallic glass

The high density of aluminum nanocrystals (>10 21 m -3 ) that develop during the primary crystallization in Al-based metallic glasses indicates a high nucleation rate (~10 18 m -3 s -1 ). Several studies have been advanced to account for the primary crystallization behavior, but none have been developed to completely describe the reaction kinetics. Recently, structural analysis by fluctuation electron microscopy has demonstrated the presence of the Al-like medium range order (MRO) regions as a spatial heterogeneity in as-spun Al 88 Y 7 Fe 5 metallic glass that is representative for the class of Al-based amorphous alloys that develop Al nanocrystals during primary crystallization. From the structural characterization, an MRO seeded nucleation configuration is established, whereby the Al nanocrystals are catalyzed by the MRO core to decrease the nucleation barrier. The MRO seeded nucleation model and the kinetic data from the delay time (τ) measurement provide a full accounting of the evolution of the Al nanocrystal density (N v ) during the primary crystallization under isothermal annealing treatments. Moreover, the calculated values of the steady state nucleation rates (J ss ) predicted by the nucleation model agree with the experimental results. Moreover, the model satisfies constraints on the structural, thermodynamic, and kinetic parameters, such as the critical nucleus size, the interface energy, and the volume-free energy driving force that are essential for a fully self-consistent nucleation kinetics analysis. In conclusion, the nucleation kinetics model can be applied more broadly to materials that are characterized by the presence of spatial heterogeneities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theoretical investigation of braking of tearing mode rotation by resistive walls in ITER

The locking of the 2/1 tearing mode to the resistive wall in the ITER tokamak (15 MA inductive scenario 2) is investigated theoretically using a cylindrical asymptotic matching model. The model takes into account the fact that ITER plasmas will effectively be surrounded by two walls; the inner blanket module layer with a time constant of about 23 ms, and the outer vacuum vessel with a time constant of about 380 ms. The model also takes cognizance of the fact that neither the blanket module layer nor the vacuum vessel can be accurately described as “thin” walls (in the ordinarily accepted sense). The model incorporates changes in both the plasma poloidal and the toroidal angular velocity profiles, in response to the electromagnetic braking torque that develops at the rational surface, because it turns out that neoclassical poloidal flow-damping is not strong enough to completely suppress changes in the poloidal velocity. Finally, the model accurately calculates changes in the poloidal and toroidal plasma angular velocity profiles by evolving the full angular equations of motion, taking the electromagnetic braking torque, plasma inertia, plasma viscosity, and poloidal flow-damping into account. The time required for the 2/1 tearing mode to grow from a small amplitude to a sufficient one to lock to the walls is found to be about 3.5 s. As a result, the critical full radial island width at which wall locking is triggered is found to be about 9% of the plasma minor radius.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modification of Lie's transform perturbation theory for charged particle motion in a magnetic field

It is pointed out that the conventional Lie transform perturbation theory for the guiding center motion of charged particles in a magnetic field needs to be modified for ordering inconsistency. There are two reasons. First, the ordering difference between the temporal variation of gyrophase and that of the other phase space coordinates needs to be taken into account. Second, it is also important to note that the parametric limit of the derivative of a function is not equivalent to the derivative of the limit function. When these facts are taken into account, the near identity transformation rule for one form related to the Lagrangian is modified. With the modified near identity transformation rule, the drift motion of charged particles can be described in the first order, instead of the second order and beyond through a tedious expansion process as in the conventional formulation. Here, this resolves the discrepancy between the direct and Lie transform treatments in the Lagrangian perturbation theory for charged particle motion in a magnetic field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Shock state distributions in porous tantalum and characterization with multipoint velocimetry

Heterogenous materials under shock compression can be expected to reach different shock states throughout the material according to local differences in microstructure and the history of wave propagation. Here, in this work, a compact, multiple-beam focusing optic assembly is used with high-speed velocimetry to interrogate the shock response of porous tantalum films prepared through thermal-spray deposition. The distribution of particle velocities across a shocked interface is compared to results obtained using a set of defocused interferometric beams that sampled the shock response over larger areas. The two methods produced velocity distributions along the shock plateau with the same mean, while a larger variance was measured with narrower beams. The finding was replicated using three-dimensional, mesoscopically resolved hydrodynamics simulations of solid tantalum with a pore structure mimicking statistical attributes of the material and accounting for radial divergence of the beams, with agreement across several impact velocities. Accounting for pore morphology in the simulations was found to be necessary for replicating the rise time of the shock plateau. The validated simulations were then used to show that while the average velocity along the shock plateau could be determined accurately with only a few interferometric beams, accurately determining the width of the velocity distribution, which here was approximately Gaussian, required a beam dimension much smaller than the spatial correlation lengthscale of the velocity field, here by a factor of ∼30×, with implications for the study of other porous materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stationary Bragg reflection of laser light in inhomogeneous absorbing plasmas inside inertial confinement fusion Hohlraums

Laser-produced plasma in inertial confinement fusion (ICF) Hohlraums are marked with density non-uniformity whose length scale can go down to micrometers. This scale is of the order of the laser wavelength. The WKB approximation, which is classically used in radiation-hydrodynamic codes to compute the laser trajectory, cannot correctly take into account such small-scale inhomogeneity of the plasma. Going beyond this approximation, we predict a novel mechanism for the laser reflection. Here we show that an electromagnetic plane wave with wave number k resonates with the k B =2 k Fourier component of a multimode perturbation of the background density and generates a reflected wave. It is the first time that this reflection is considered for stationary inhomogeneous ICF plasmas, and the energy absorption is taken into account. This mechanism, which is a form of Bragg reflection, can occur away from the critical surface and generate a drift of the location of the laser absorption. Furthermore, this absorption will be periodically modulated with a k B wave number. The stationary Bragg reflection can explain ongoing discrepancies between experimental and numerical data about laser trajectory and absorption in ICF Hohlraums.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Shaded fraction and backtracking in single-axis trackers on rolling terrain

A generalized closed-form equation for the shaded collector fraction in solar arrays on rolling or undulating terrain is provided for single-axis tracking and fixed-tilt systems. The equation accounts for different rotation angles between the shaded and shading trackers, cross-axis slope between the two trackers, and offset between the collector plane and axis of rotation. The validity of the equation is demonstrated through comparison with numerical ray-tracing simulations and remaining minor sources of error are quantified. Additionally, a simple procedure to determine backtracking rotations for each row in an array installed on the rolling terrain (varying in the direction perpendicular to the tracker axes) is provided. The backtracking equation accounts for a desired shaded fraction (including complete shade avoidance) as well as an axis-collector offset. Test cases are provided to facilitate implementation of these equations.

14 SOLAR ENERGY↗

Monte Carlo toolkit for designing and validating step-range-filter spectrometer designs

Here, we present a Monte Carlo toolkit for validating step range filter (SRF) spectrometer designs. Geant4 is used to transport charged particles through the SRF filters to generate synthetic SRF data that include realistic CR-39 effects. Synthetic SRF spectra generated by this method inherently account for instrument response and allow for the quantification of SRF performance before shots. The usefulness of this toolkit is demonstrated through its application to a number of problems. A new broadband SRF for the ∼10 MeV wide 3He3He proton spectrum is validated, and an analysis method for analyzing 3He3He-p SRF data that accounts for instrument response is put forth. In addition, an SRF design for the compact recoil-proton spectrometer (CRS) on the Z-machine is validated. Finally, a new calibration technique for the DD-p SRF is proposed and validated.

Johnson, T. M. (ORCID:0000000193032949)↗

Particle balance of deuterium during deuterium shattered pellet injection shutdown in DIII-D

A particle balance analysis was conducted during a deuterium (D 2 ) shattered pellet injection-induced plasma shutdown on the DIII-D tokamak to determine why less than 20% of the pellet material is assimilated into the core plasma by the mid-current quench (CQ). Initially, most of the D 2 is injected as frozen shards and ionized upon entering the vessel. During the thermal quench, ionized particles move to the divertors and subsequently to the center post (CP) walls, where they rapidly recycle and partially accumulate as neutrals without assimilating into the core plasma. In contrast, the particle flux to the outer midplane walls is negligible, despite being accompanied by hot plasma with electron temperatures exceeding 100 eV. During mid-CQ, volume recombination effects, although not large enough to impact overall particle balance, were significant enough to require accounting for accurate interpretation of fast-framing camera D-alpha signals and the estimation of the CP wall particle flux. In addition, toroidal asymmetries, observed in measurements of toroidal electron density perturbations and the phase of magnetohydrodynamic modes, are present throughout the shutdown and can account for a discrepancy in the assimilation rate for up to 50% of the observed D 2 particle inventory. These sources and sinks of particles and fluxes were identified using absolutely calibrated D-alpha brightness and Langmuir probes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hydrodynamic expansion and near-infrared absorption of x-ray heated aluminum plasmas

We use x-ray pulses from dense argon plasmas at the Z Machine (Sandia National Laboratories) to generate hypersonic aluminum plasmas akin to material ejecta during proposed planetary defense missions, fusion reactor wall excursions, and other high-energy density processes. Near-infrared absorption is used to diagnose the controlled expansion of the plasmas through cylindrical cavities following their generation from x-ray heating of solid aluminum 7075 alloy. The data are compared to multidimensional radiation hydrodynamics simulations utilizing the ALEGRA multiphysics code, accounting for the dynamics of radiation scattering, material phase change, plasma expansion, thermal re-irradiation, and interactions with the cavity and with the infrared beams. To allow for accurate simulation, density functional theory is used to apply the Hagen–Rubens relation for the far-infrared and is adjoined with a detailed configuration accounting model using the Propaceos code, producing opacities spanning 10 −1 –10 4 eV photon energy for aluminum 7075 alloy, and in comparison with pure aluminum. The model is found to agree with experimental data in the higher-fluence regime when the Hagen–Rubens relation is applied. The ejected material, which is observed to travel up to 55 km/s, is comprised of a strongly ionized, non-LTE plasma front at ∼10 eV temperature followed by a weakly ionized LTE gas at higher density. The present findings lend some confidence to the broad-range equation of state and infrared opacity models for weakly ionized aluminum plasmas while demonstrating an approach to their future refinement, with potential application to astrophysical plasmas and other extreme processes.

Adiabatic process↗

Hot-spot model for inertial confinement fusion implosions with an applied magnetic field

Imposing a magnetic field on inertial confinement fusion implosions magnetizes the electrons in the compressed fuel; this suppresses thermal losses, which increases temperature and fusion yield. Indirect-drive experiments at the National Ignition Facility with 12 and 26 T applied magnetic fields demonstrate up to 40% increase in temperature, 3× increase in fusion yield, and indicate that magnetization alters the radial temperature profile [Moody et al., Phys. Rev. Lett. 129, 195002 (2022); Lahmann et al., APS DPP (2022)]. In this work, we develop a semi-analytic hot-spot model, which accounts for the two-dimensional (2D) Braginskii anisotropic heat flow due to an applied axial magnetic field. First, we show that hot-spot magnetization alters the radial temperature profile, increasing the central peakedness, which is most pronounced for moderately magnetized implosions (with 8–14 T applied field), compared to both unmagnetized (with no applied field) and highly magnetized (with 26 T or higher applied field) implosions. This model explains the trend in the experimental data, which finds a similarly altered temperature profile in the 12 T experiment. Next, we derive the hot-spot model for gas-filled (Symcap) implosions, accounting for the effects of magnetization on the thermal conduction and in changing the radial temperature (and density) profiles. Using this model, we compute predicted central temperature amplification and yield enhancement scaling with the applied magnetic field. The central temperature fits the experimental data accurately, and the discrepancy in the yield suggests a systematic (independent of applied field) degradation, such as mix, and additional degradation in the reference unmagnetized shot, such as reduced laser drive, increased implosion asymmetry, or the magnetic field suppressing ablator mixing into the hot-spot.

Alpha particles↗

Extended molecular eigenmodes treatment of dipole–dipole NMR relaxation in real fluids

Traditional models of NMR relaxation fail to account for the complex, multi-exponential behavior of the autocorrelation function in realistic systems characterized by soft-interactions and molecules that are chemically and physically complex. Here, in this study, we describe the relative diffusion of the spin dipoles by means of a Fokker–Planck equation that includes an interaction potential of mean force to account for the response of the physical/chemical environment around the dipoles. By numerically solving the Fokker–Planck equation for the diffusion propagator, we estimate dipole–dipole NMR relaxation for like- and unlike-spin systems via its eigenmode solution. We test the model against molecular simulations of diffusing dipoles with harmonic potentials and also validate using experimental longitudinal relaxation data from real systems, including Gd(III)–aqua and Gd(III)–DO3A–butrol complexes, the latter being an important MRI contrast agent. Using this novel approach, we predict both the inner- and outer-shell contributions to the relaxivity rates with excellent accuracy at frequencies relevant to MRI. We also show that, under the appropriate assumptions, our framework naturally recovers the Bloembergen–Purcell–Pound, the Solomon–Bloembergen–Morgan, and the Hwang–Freed models. Our implementation is general and publicly available for application to a broad range of systems.

Pinheiro dos Santos, Thiago J. [Rice Univ., Housto↗

Revisiting the Anisotropic Complex Refractive Indices of Sodium Nitrate for Interpretation of the Reflectance Spectrum of Pressed Pellets

Reflectance spectroscopy is notoriously confounding in that the measured spectral response is highly dependent upon the morphology of the sample. Fortunately, all such perturbations are neatly encoded by the complex refractive index of the sample. Herein, we seek to quantitatively model the measured infrared reflectance spectrum of a specularly flat pressed pellet sample of the birefringent compound, sodium nitrate. Single crystals of sodium nitrate were synthesized via a slow evaporation process and spectroscopically analyzed using polarization-dependent infrared single-angle reflectance spectroscopy. The anisotropic complex refractive index was measured from 7500 to 300 cm-1 (1.33 to 33.33 µm). The deduced anisotropic optical constants were found to be consistent with those previously reported. Once measured and validated, the optical constants were applied to model the pressed pellet reflectance spectrum. It was evident that an average of the anisotropic refractive indices was insufficient to account for the measured pellet reflectance. In order to account for contributions of all possible microcrystalline orientations within the pellet, the Python package PYELLI was used to calculate a basis set of orientation-dependent reflection spectra spanning the distinct ? and ? Euler rotations of the uniaxial crystal. When the population of orientations was allowed to vary freely in a spectral fit analysis, the fit-deduced orientations were tightly clustered along f = 45º, hinting at residual anisotropy in the pressed pellet sample. Conversely, an equally valid spectral fit (with marginally worse fit metric) was obtained when the population was constrained to an isotropic distribution of orientations. Subsequent non-zero cross-polarization reflectance measurements likewise suggested anisotropy in the pellet. However, both grazing incidence wide-angle x-ray scattering and scanning electron microscopy measurements revealed that the microcrystal orientations at the surface of the pressed pellet sample were isotropically distributed (and that the average crystallite size was larger than ?/10). Application of the measured complex refractive indices for modeling the reflectance spectrum of the pressed pellet, and rectification of these seemingly contradictory observations will be discussed.

Wilhelm, Michael J.↗

Optimal binning of correlated measurements

Experimental measurements are commonly represented on a discrete grid, requiring a balance between granularity and statistical noise. Two strategies have traditionally been used to improve such representations: selecting an appropriate bin width to control discretization error and applying kernel-based smoothing to suppress fluctuations. Despite their shared goal, these approaches have largely developed independently, without a unified statistical description of how discretization and correlation jointly determine measurement precision. Here, we extend the discussion of optimal interval averaging to a correlation-aware setting by Gaussian process regression, which explicitly accounts for correlations among neighboring bins. Starting from first principles, we derive the mean-squared error of discretized measurements and obtain closed-form asymptotic expressions for the optimal bin width and correlation length. When recast in reduced variables, the theory reveals distinct universal scaling laws governing the error in the correlation-free and correlation-controlled regimes. Characterized by intrinsically smooth intensity profiles and counting-based statistics, neutron scattering measurements are well suited for demonstrating the enhanced error contraction enabled by inter-bin correlations. We show that such improvement is achievable over the experimentally accessible Q-range and across multiple instruments and material systems. These results show that explicitly accounting for correlations systematically reshapes the limits of precision in discretized, noise-limited measurements. More broadly, the framework provides a transferable statistical foundation for optimizing data representation, inference, and experimental design across the physical and data sciences.

Tung, Chi-Huan [ORNL] (ORCID:0000000221972074)↗

On “incorrectnes” of the rate equation for photon density of semiconductor lasers

The authors of works [1, 2] published in Quantum Electronics assert that the rate equation for photon density with a term taking into account the contribution of spontaneous emission to the laser mode is erroneous (i.e., the spontaneous term should be excluded from consideration) and criticise works [3, 4]. In the present work, we analyse previously published papers and show that the necessity of taking into account the contribution of spontaneous emission to the laser mode follows from the quantum-mechanical rate equation for semiconductor lasers. The results obtained in [1, 2] are analysed. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

High concentration of ultrafine particles in the Amazon free troposphere produced by organic new particle formation

The large concentrations of ultrafine particles consistently observed at high altitudes over the tropics represent one of the world’s largest aerosol reservoirs, which may be providing a globally important source of cloud condensation nuclei. However, the sources and chemical processes contributing to the formation of these particles remain unclear. Here we investigate new particle formation (NPF) mechanisms in the Amazon free troposphere by integrating insights from laboratory measurements, chemical transport modeling, and field measurements. To account for organic NPF, we develop a comprehensive model representation of the temperature-dependent formation chemistry and thermodynamics of extremely low volatility organic compounds as well as their roles in NPF processes. We find that pure-organic NPF driven by natural biogenic emissions dominates in the uppermost troposphere above 13 km and accounts for 65 to 83% of the column total NPF rate under relatively pristine conditions, while ternary NPF involving organics and sulfuric acid dominates between 8 and 13 km. The large organic NPF rates at high altitudes mainly result from decreased volatility of organics and increased NPF efficiency at low temperatures, somewhat counterbalanced by a reduced chemical formation rate of extremely low volatility organic compounds. These findings imply a key role of naturally occurring organic NPF in high-altitude preindustrial environments and will help better quantify anthropogenic aerosol forcing from preindustrial times to the present day.

54 ENVIRONMENTAL SCIENCES↗

Intracellular pathways for lignin catabolism in white-rot fungi

White-rot fungi, the most efficient organisms at breaking down lignin from wood in Nature, utilize lignin degradation products as a carbon source. 15 Lignin is a biopolymer found in plant cell walls that accounts for 30% of the organic carbon in the biosphere. White-rot fungi (WRF) are considered the most efficient organisms at degrading lignin in Nature. While lignin depolymerization by WRF has been exhaustively studied, the possibility that WRF are able to utilize lignin as a carbon source is still a matter of controversy. Here we employ 13C-labeling and systems biology approaches to demonstrate that two WRF, Trametes versicolor and Gelatoporia subvermispora, funnel lignin-derived aromatic compounds into central carbon metabolism via intracellular catabolic pathways. These results provide insights into global carbon cycling in 20 soil ecosystems, and furthermore establishes a foundation for employing WRF in simultaneous lignin depolymerization and bioconversion to bioproducts – a key step towards enabling a sustainable bioeconomy.

Del Cerro, Carlos↗

Shipping regulations lead to large reduction in cloud perturbations

Global shipping accounts for 13% of global emissions of SO2, which, once oxidized to sulfate aerosol, acts to cool the planet both directly by scattering sunlight and indirectly by increasing the albedo of clouds. This cooling due to sulfate aerosol offsets some of the warming effect of greenhouse gasses and is the largest uncertainty in determining the change in the Earth’s radiative balance by human activity. Ship tracks—the visible manifestation of the indirect of effect of ship emissions on clouds as quasi-linear features—have long provided an opportunity to quantify these effects. However, they have been arduous to catalog and typically studied only in particular regions for short periods of time. Using a machine-learning algorithm to automate their detection we catalog more than 1 million ship tracks to provide a global climatology. We use this to investigate the effect of stringent fuel regulations introduced by the International Maritime Organization in 2020 on their global prevalence since then, while accounting for the disruption in global commerce caused by COVID-19. We find a marked, but clearly nonlinear, decline in ship tracks globally: An 80% reduction in SO x emissions causes only a 25% reduction in the number of tracks detected.

54 ENVIRONMENTAL SCIENCES↗

Lagrangian large eddy simulations via physics-informed machine learning

High-Reynolds number homogeneous isotropic turbulence (HIT) is fully described within the Navier–Stokes (NS) equations, which are notoriously difficult to solve numerically. Engineers, interested primarily in describing turbulence at a reduced range of resolved scales, have designed heuristics, known as large eddy simulation (LES). LES is described in terms of the temporally evolving Eulerian velocity field defined over a spatial grid with the mean-spacing correspondent to the resolved scale. This classic Eulerian LES depends on assumptions about effects of subgrid scales on the resolved scales. Here, we take an alternative approach and design LES heuristics stated in terms of Lagrangian particles moving with the flow. Our Lagrangian LES, thus L-LES, is described by equations generalizing the weakly compressible smoothed particle hydrodynamics formulation with extended parametric and functional freedom, which is then resolved via Machine Learning training on Lagrangian data from direct numerical simulations of the NS equations. The L-LES model includes physics-informed parameterization and functional form, by combining physics-based parameters and physics-inspired Neural Networks to describe the evolution of turbulence within the resolved range of scales. The subgrid-scale contributions are modeled separately with physical constraints to account for the effects from unresolved scales. We build the resulting model under the differentiable programming framework to facilitate efficient training. We experiment with loss functions of different types, including physics-informed ones accounting for statistics of Lagrangian particles. We show that our L-LES model is capable of reproducing Eulerian and unique Lagrangian turbulence structures and statistics over a range of turbulent Mach numbers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗