Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “OPACITY”

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 721 records · Page 40

Steamworlds: Atmospheric Structure and Critical Mass of Planets Accreting Icy Pebbles

In the core accretion model, gas-giant planets first form a solid core, which then accretes gas from a protoplanetary disk when the core exceeds a critical mass. Here, we model the atmosphere of a core that grows by accreting ice-rich pebbles. The ice fraction of pebbles evaporates in warm regions of the atmosphere, saturating it with water vapor. Excess water precipitates to lower altitudes. Beneath an outer radiative region, the atmosphere is convective, following a moist adiabat in saturated regions due to water condensation and precipitation. Atmospheric mass, density, and temperature increase with core mass. For nominal model parameters, planets with core masses (ice + rock) between 0.08 and 0.16 Earth masses have surface temperatures between 273 and 647 K and form an ocean. In more massive planets, water exists as a supercritical convecting fluid mixed with gas from the disk. Typically, the core mass reaches a maximum (the critical mass) as a function of the total mass when the core is 2–5 Earth masses. The critical mass depends in a complicated way on pebble size, mass flux, and dust opacity due to the occasional appearance of multiple core-mass maxima. The core mass for an atmosphere of 50% hydrogen and helium may be a more robust indicator of the onset of gas accretion. This mass is typically 1–3 Earth masses for pebbles that are 50% ice by mass, increasing with opacity and pebble flux and decreasing with pebble ice/rock ratio.

79 ASTRONOMY AND ASTROPHYSICS↗

SUPERLUMINOUS SUPERNOVA SN 2015bn IN THE NEBULAR PHASE: EVIDENCE FOR THE ENGINE-POWERED EXPLOSION OF A STRIPPED MASSIVE STAR

We present nebular-phase imaging and spectroscopy for the hydrogen-poor superluminous supernova (SLSN) SN 2015bn, at redshift z = 0.1136, spanning +250–400 days after maximum light. The light curve exhibits a steepening in the decline rate from 1.4 mag (100 days){sup −1} to 1.7 mag (100 days){sup −1}, suggestive of a significant decrease in the opacity. This change is accompanied by a transition from a blue continuum superposed with photospheric absorption lines to a nebular spectrum dominated by emission lines of oxygen, calcium, and magnesium. There are no obvious signatures of circumstellar interaction or large {sup 56}Ni mass. We show that the spectrum at +400 days is virtually identical to a number of energetic SNe Ic such as SN 1997dq, SN 2012au, and SN 1998bw, indicating similar core conditions and strengthening the link between “hypernovae”/long gamma-ray bursts and SLSNe. A single explosion mechanism may unify these events that span absolute magnitudes of −22 < M {sub B} < −17. Both the light curve and spectrum of SN 2015bn are consistent with an engine-driven explosion ejecting 7–30 M {sub ⊙} of oxygen-dominated ejecta (for reasonable choices in temperature and opacity). A strong and relatively narrow O i λ 7774 line, seen in a number of these energetic events but not in normal supernovae, may point to an inner shell that is the signature of a central engine.

79 ASTRONOMY AND ASTROPHYSICS↗

The impact of unresolved turbulence on the escape fraction of Lyman Continuum photons

We investigate the relation between the turbulent Mach number (M) and the escape fraction of Lyman continuum photons (f{sub esc}) in high-redshift galaxies. Approximating the turbulence as isothermal and isotropic, we show that the increase in the variance in column densities from M=1 to M=10 causes f{sub esc} to increase by ≈25%, and the increase from M=1 to M=20 causes f{sub esc} to increases by ≈50% for a medium with opacity τ≈1. At a fixed Mach number, the correction factor for escape fraction relative to a constant column density case scales exponentially with the opacity in the cell, which has a large impact for simulated star-forming regions. Furthermore, in simulations of isotropic turbulence with full atomic/ionic cooling and chemistry, the fraction of HI drops by a factor of ≈2.5 at M≈10 even when the mean temperature is ≈5×10{sup 3} K. If turbulence is unresolved, these effects together enhance f{sub esc} by a factor >3 at Mach numbers above 10. Such Mach numbers are common at high redshifts where vigorous turbulence is driven by supernovae, gravitational instabilities, and merger activity, as shown both by numerical simulations and observations. These results, if implemented in the current hydrodynamical cosmological simulations to account for unresolved turbulence, can boost the theoretical predictions of the Lyman Continuum photon escape fraction and further constrain the sources of reionization.

79 ASTRONOMY AND ASTROPHYSICS↗

Infrared Emission from Kilonovae: The Case of the Nearby Short Hard Burst GRB 160821B

In this paper, we present constraints on Ks-band emission from one of the nearest short hard gamma-ray bursts, GRB 160821B, at z = 0.16, at three epochs. We detect a red relativistic afterglow from the jetted emission in the first epoch but do not detect any excess kilonova emission in the second two epochs. We compare upper limits obtained with Keck I/MOSFIRE to multi-dimensional radiative transfer models of kilonovae, that employ composition-dependent nuclear heating and LTE opacities of heavy elements. We discuss eight models that combine toroidal dynamical ejecta and two types of wind and one model with dynamical ejecta only. We also discuss simple, empirical scaling laws of predicted emission as a function of ejecta mass and ejecta velocity. Our limits for GRB 160821B constrain the ejecta mass to be lower than 0.03 M ⊙ for velocities greater than 0.1 c. At the distance sensitivity range of advanced LIGO, similar ground-based observations would be sufficiently sensitive to the full range of predicted model emission including models with only dynamical ejecta. Finally, the color evolution of these models shows that I–K color spans 7–16 mag, which suggests that even relatively shallow infrared searches for kilonovae could be as constraining as optical searches.

79 ASTRONOMY AND ASTROPHYSICS↗

A new, coupled transport-diffusion method for radiative transfer calculations

We derive and present a new frequency- and angle-integrated low-order system of equations designed to enhance the accuracy of a coupled, high-order (transport) solution of the thermal radiative transfer equations. In particular, our new low-order system is designed to use intensity-weighted opacities and anisotropic diffusion coefficients generated by a solution of the Implicit Monte Carlo (IMC) equations in order to predict the spatial dependence of the material temperature and radiation energies in the ensuing time cycle. The predicted temperature solution can then be exploited to generate appropriately time-centered opacities, specific heats, and Planck emission spectra for the upcoming IMC solution. Additionally, the relatively inexpensive solution of the low-order system can be iteratively solved to recommend an adaptive time step size before the IMC solution is computed. A test implementation has been implemented using existing software available from the Jayenne and Capsaicin projects at Los Alamos National Laboratory. We present initial results from a new driver code that has integrated these stochastic and deterministic software packages. (authors)

97 MATHEMATICS AND COMPUTING↗

Modelling Thomson Scattering in a Hydrogen Plasma at Stellar Interior Conditions Using the Hypernetted‐Chain Approach

Under the extreme conditions found in small stars, where electron degeneracy and Coulomb coupling are significant, accurate modeling of Thomson scattering is crucial for determining opacity, a primary quantity for stellar energy transport. We use hypernetted-chain calculations, incorporating quantum pseudopotentials and electron-exchange effects to obtain the electron–electron static structure factor to calculate the Thomson scattering transport cross-section for conditions prevailing in the interior of small stars. These results are compared to those from average-atom simulations and analytical calculations. Our findings support laboratory astrophysics experiments aimed at benchmarking opacity models for stellar interiors, particularly for red dwarf stars, and help to bridge theoretical models with observations.

ICF↗

Lorentzian dynamics and factorization beyond rationality

We investigate the emergence of topological defect lines in the conformal Regge limit of two-dimensional conformal field theory. We explain how a local operator can be factorized into a holomorphic and an anti-holomorphic defect operator connected through a topological defect line, and discuss implications on analyticity and Lorentzian dynamics including aspects of chaos. We derive a formula relating the infinite boost limit, which holographically encodes the “opacity” of bulk scattering, to the action of topological defect lines on local operators. Leveraging the unitary bound on the opacity and the positivity of fusion coefficients, we show that the spectral radii of a large class of topological defect lines are given by their loop expectation values. Factorization also gives a formula relating the local and defect operator algebras and fusion categorical data. We then review factorization in rational conformal field theory from a defect perspective, and examine irrational theories. On the orbifold branch of the c = 1 free boson theory, we find a unified description for the topological defect lines through which the twist fields are factorized; at irrational points, the twist fields factorize through “non-compact” topological defect lines which exhibit continuous defect operator spectra. Along the way, we initiate the development of a formalism to characterize non-compact topological defect lines.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A Variable Eddington Factor method for the 1-D grey radiative transfer equations with discontinuous Galerkin and mixed finite-element spatial differencing

The purpose of this paper is to present a Variable Eddington Factor (VEF) method for the 1-D grey radiative transfer equations that uses a lumped linear discontinuous Galerkin spatial discretization for the Sequations together with a constant-linear mixed finite-element discretization for the VEF moment and material temperature equations. The use of independent discretizations can be particularly useful for multiphysics applications such as radiation-hydrodynamics. The VEF method is quite old, but to our knowledge, this particular combination of differencing schemes has not been previously investigated for radiative transfer. We define the scheme and present computational results. As expected, the scheme exhibits second-order accuracy for the directionally-integrated intensity and material temperature, and behaves well in the thick diffusion limit. An important focus of this study is the treatment of the strong temperature dependence of the opacities and the spatial dependence of the opacities within each cell, which are not explicitly defined by the basic discretization schemes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Proper orthogonal decomposition based reduced-order modeling of flux-Limited gray thermal radiation

Here, in this work, a proper orthogonal decomposition (POD) based reduced-order model (ROM) is developed to solve gray, flux-limited thermal radiation diffusion. We focus on the variable opacity radiation penetration benchmark posed by Olson, Auer, and Hall. The T -3 relationship for opacity in conjunction with high-temperature radiation penetrating an initially cold material produces a strong thermal radiation shock. This class of problems is particularly challenging for standard POD-based reduced-order modeling due to the nonlinearities presented by 1) the T 4 source term, and 2) flux-limited diffusion operator. To address these challenges and develop a cost competitive ROM, we employ a “hyper-reduction” technique through discrete empirical interpolation (DEIM) and allow for adaptive reduced-order projections through principal interval decomposition (PID). Performance of the proposed methodology is quantified by comparing the cost savings and accuracy relative to a full-order computation. Reference solutions and snapshot data are obtained through a full-order calculation performed by the University of Chicago maintained astrophysics code, FLASH. For consistency and potential extensibility, the developed ROM is also implemented in FLASH. We find that in the initialization regime, where the thermal radiation wave is initially created by the warming of the material, this class of problems is highly reducible and suitable for POD-based ROMs. However, the strong convective nature of the wave propagation regime is less reducible and more challenging to create an efficient ROM.

42 ENGINEERING↗

Extracting a mixing parameter from 2D radiographic imaging of variable-density turbulent flow

We extract a suitably averaged fluctuating density from the two-dimensional radiographic image of a flow. The X-ray attenuation is given by the Beer–Lambert law which exponentially damps the incident beam intensity by a factor proportional to the density, opacity and thickness of the target. By making reasonable assumptions for the mean density, opacity and effective thickness of the target flow, we estimate the density fluctuation contribution to the attenuation. The extracted density fluctuations averaged across the thickness of the flow in the direction of the beam may be used to form the density–specific-volume correlation $b$. In a statistical description of variable-density turbulence, $b$ quantifies the degree of mixedness. The ability to extract a measure of mixedness from experimental data would be a powerful tool that could be used in the validation of mix models. The scheme proposed is tested for DNS data computed for variable density buoyancy-driven mixing. We quantify the deficits in the extracted value of $b$ due to target thickness, Atwood number and modeled signal noise. This analysis justifies using the proposed scheme to infer the mix parameter from thin targets at moderate to low Atwood numbers. Furthermore, to illustrate how the scheme might be used in a practical problem, we demonstrate its application to a radiographic image of counter-shear flow obtained from experiments at the National Ignition Facility.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Understanding parton evolution in matter from renormalization group analysis

We perform a renormalization group (RG) analysis of collinear hadron production in deep inelastic scattering on nuclei. We consider the limit where the parent parton energy E is large, while the medium opacity remains small. We identify the fixed order and leading enhanced medium contributions to the semi-inclusive cross sections and derive RG equations that resum multiple emissions near the endpoints of the splitting functions at first order in opacity. These evolution equations treat the same type of radiation enhancement in matter as the modified Dokshitzer-Gribov-Lipatov-Altarelli-Parisi approach, but differ in the way one regulates the collinear divergences. They provide a unique analytic insight into the problem of resummation and a faster and more efficient path to phenomenology. The new RG evolution framework is applied to study fragmentation in eA reactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Simulating x-ray Thomson scattering signals from high-density, millimetre-scale plasmas at the National Ignition Facility

We have developed a model for analysing x-ray Thomson scattering data from high-density, millimetre-scale inhomogeneous plasmas created during ultra-high pressure implosions at the National Ignition Facility in a spherically convergent geometry. The density weighting of the scattered signal and attenuation of the incident and scattered x-rays throughout the target are included using radial profiles of the density, opacity, ionization state, and temperature provided by radiation-hydrodynamics simulations. These simulations show that the scattered signal is strongly weighted toward the bulk of the shocked plasma and the Fermi degenerate material near the ablation front. We show that the scattered signal provides a good representation of the temperature of this highly nonuniform bulk plasma and can be determined to an accuracy of ca. 15% using typical data analysis techniques with simple 0D calculations. On the other hand, the mean ionization of the carbon in the bulk is underestimated. We suggest that this discrepancy is due to the convolution of scattering profiles from different regions of the target. Subsequently, we discuss modifications to the current platform to minimise the impact of inhomogeneities, as well as opacity, and also to enable probing of conditions more strongly weighted toward the compressed core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Design calculations for a xenon plasma x-ray shield to protect the NIF optical Thomson scattering diagnostic

An Optical Thomson Scattering (OTS) diagnostic is currently being developed for the National Ignition Facility (NIF) at Lawrence Livermore National Laboratory. This diagnostic is designed to make measurements of the hohlraum plasma parameters, such as the electron temperature and the density, during inertial confinement fusion (ICF) experiments. NIF ICF experiments present a very challenging environment for optical measurements; by their very nature, hohlraums produce intense soft x-ray emission, which can cause “blanking” (radiation induced opacity) of the radiation facing optical components. The soft x-ray fluence at the surface of the OTS blast shield, 60 cm from the hohlraum, is estimated to be ∼8 J cm{sup −2}. This is significantly above the expected threshold for the onset of “blanking” effects. A novel xenon plasma x-ray shield is proposed to protect the blast shield from x-rays and mitigate “blanking.” Estimates suggest that an areal density of 10{sup 19} cm{sup −2} Xe atoms will be sufficient to absorb 99.5% of the soft x-ray flux. Two potential designs for this shield are presented.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Characterization of Agfa Structurix series D4 and D3sc x-ray films in the 0.7–4.6 keV energy range

X-ray films remain a key asset for high-resolution x-ray spectral imaging in high-energy-density experiments conducted at the National Ignition Facility (NIF). The soft x-ray Opacity Spectrometer (OpSpec) fielded at the NIF has an elliptically shaped crystal design that measures x rays in the 900–2100 eV range and currently uses an image plate as the detecting medium. However, Agfa D4 and D3sc x-ray films’ higher spatial resolution provides increased spectral resolution to the data over the IP-TR image plates, driving the desire for regular use of x-ray film as a detecting medium. The calibration of Agfa D4 x-ray film for use in the OpSpec is communicated here. These calibration efforts are vital to the accuracy of the NIF opacity measurements and are conducted in a previously un-studied x-ray energy range under a new film development protocol required by NIF. The absolute response of Agfa D4 x-ray film from 705 to 4620 eV has been measured using the Nevada National Security Site Manson x-ray source. A broader range of energies was selected to compare results with previously published data. The measurements were taken using selected anodes, filters, and applied voltages to produce well-defined energy lines.

Dutra, E. C. (ORCID:0000000167649195)↗

Platform for probing radiation transport properties of hydrogen at conditions found in the deep interiors of red dwarfs

We describe an experimental concept at the National Ignition Facility for specifically tailored spherical implosions to compress hydrogen to extreme densities (up to ∼800× solid density, electron number density ne∼4×1025 cm−3) at moderate temperatures (T∼200 eV), i.e., to conditions, which are relevant to the interiors of red dwarf stars. The dense plasma will be probed by laser-generated x-ray radiation of different photon energy to determine the plasma opacity due to collisional (free–free) absorption and Thomson scattering. The obtained results will benchmark radiation transport models, which in the case for free–free absorption show strong deviations at conditions relevant to red dwarfs. This very first experimental test of free–free opacity models at these extreme states will help to constrain where inside those celestial objects energy transport is dominated by radiation or convection. Moreover, our study will inform models for other important processes in dense plasmas, which are based on electron–ion collisions, e.g., stopping of swift ions or electron–ion temperature relaxation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Statistical data analysis of x-ray spectroscopy data enabled by neural network accelerated Bayesian inference

Bayesian inference applied to x-ray spectroscopy data analysis enables uncertainty quantification necessary to rigorously test theoretical models. However, when comparing to data, detailed atomic physics and radiation transfer calculations of x-ray emission from non-uniform plasma conditions are typically too slow to be performed in line with statistical sampling methods, such as Markov Chain Monte Carlo sampling. Furthermore, differences in transition energies and x-ray opacities often make direct comparisons between simulated and measured spectra unreliable. Here, we present a spectral decomposition method that allows for corrections to line positions and bound–bound opacities to best fit experimental data, with the goal of providing quantitative feedback to improve the underlying theoretical models and guide future experiments. In this work, we use a neural network (NN) surrogate model to replace spectral calculations of isobaric hot-spots created in Kr-doped implosions at the National Ignition Facility. The NN was trained on calculations of x-ray spectra using an isobaric hot-spot model post-processed with Cretin, a multi-species atomic kinetics and radiation code. The speedup provided by the NN model to generate x-ray emission spectra enables statistical analysis of parameterized models with sufficient detail to accurately represent the physical system and extract the plasma parameters of interest.

47 OTHER INSTRUMENTATION↗

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↗

Resolving discrepancies in bang-time predictions for indirect-drive ICF experiments on the NIF: Insights from the Build-A-Hohlraum campaign

This study investigated discrepancies between measured and simulated x-ray drive in Indirect-Drive Inertial Confinement Fusion (ID-ICF) hohlraums at the National Ignition Facility. Despite advances in radiation-hydrodynamic simulations, a consistent “drive deficit” remains. Experimentally measured ID-ICF capsule bang-times are systematically 400–700 ps later than simulations predict. The Build-A-Hohlraum (BAH) campaign explored potential causes for this discrepancy by systematically varying hohlraum features, including laser entrance hole (LEH) windows, capsules, and gas fills. Overall, the agreement between simulated and experimental x-ray drive was found to be largely unaffected by these changes. The data allow us to exclude some hypotheses put forward to potentially explain the discrepancy. Errors in the local thermodynamic equilibrium (LTE) atomic modeling, errors in the modeling of LEH closure, and errors due to a lack of plasma species mix physics in simulations are shown to be inconsistent with our measurements. Instead, the data support the hypothesis that errors in NLTE emission modeling are a significant contributor to the discrepancy. X-ray emission in the 2–4 keV range is found to be approximately 30% lower than in simulations. This is accompanied by higher than predicted electron temperatures in the gold bubble region, pointing to errors in non-LTE modeling. Introducing an opacity multiplier of 0.87 on energy groups above 1.8 keV improves agreement with experimental data, reducing the bang-time discrepancy from 300 to 100 ps. These results underscore the need for refined NLTE opacity models to enhance the predictive power of hohlraum simulations.

Band emission↗