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 1,207 records · Page 67

A unified radiative magnetohydrodynamics code for lightning-like discharge simulations

A two-dimensional Eulerian finite difference code is developed for solving the non-ideal magnetohydrodynamic (MHD) equations including the effects of self-consistent magnetic field, thermal conduction, resistivity, gravity, and radiation transfer, which when combined with specified pulse current models and plasma equations of state, can be used as a unified lightning return stroke solver. The differential equations are written in the covariant form in the cylindrical geometry and kept in the conservative form which enables some high-accuracy shock capturing schemes to be equipped in the lightning channel configuration naturally. In this code, the 5-order weighted essentially non-oscillatory scheme combined with Lax-Friedrichs flux splitting method is introduced for computing the convection terms of the MHD equations. The 3-order total variation diminishing Runge-Kutta integral operator is also equipped to keep the time-space accuracy of consistency. The numerical algorithms for non-ideal terms, e.g., artificial viscosity, resistivity, and thermal conduction, are introduced in the code via operator splitting method. This code assumes the radiation is in local thermodynamic equilibrium with plasma components and the flux limited diffusion algorithm with grey opacities is implemented for computing the radiation transfer. The transport coefficients and equation of state in this code are obtained from detailed particle population distribution calculation, which makes the numerical model is self-consistent. This code is systematically validated via the Sedov blast solutions and then used for lightning return stroke simulations with the peak current being 20 kA, 30 kA, and 40 kA, respectively. The results show that this numerical model consistent with observations and previous numerical results. The population distribution evolution and energy conservation problems are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Calibration of a gated flat field spectrometer as a function of x-ray intensity

We present an experimental determination of the response of a gated flat-field spectrometer at the Shenguang-II laser facility. X-rays were emitted from a target that was heated by laser beams and then were divided into different intensities with a step aluminum filter and collected by a spectrometer. The transmission of the filter was calibrated using the Beijing Synchrotron Radiation Facility. The response characteristics of the spectrometer were determined by comparing the counts recorded by the spectrometer with the relative intensities of the x-rays transmitted through the step aluminum filter. The response characteristics were used to correct the transmission from two shots of an opacity experiment using the same samples. The transmissions from the two shots are consistent with corrections, but discrepant without corrections.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Recent advances in theoretical and numerical studies of wire array Z-pinch in the IAPCM

Fast Z-pinch has produced the most powerful X-ray radiation source in laboratory and also shows the possibility to drive inertial confinement fusion (ICF). Recent advances in wire-array Z-pinch researches at the Institute of Applied Physics and Computational Mathematics are presented in this paper. A typical wire array Z-pinch process has three phases: wire plasma formation and ablation, implosion and the MRT instability development, stagnation and radiation. A mass injection model with azimuthal modulation coefficient is used to describe the wire initiation, and the dynamics of ablated plasmas of wire-array Z-pinches in (r, θ) geometry is numerically studied. In the implosion phase, a two-dimensional(r, z) three temperature radiation MHD code MARED has been developed to investigate the development of the Magneto-Rayleigh-Taylor(MRT) instability. We also analyze the implosion modes of nested wire-array and find that the inner wire-array is hardly affected before the impaction of the outer wire-array. While the plasma accelerated to high speed in the implosion stage stagnates on the axis, abundant x-ray radiation is produced. The energy spectrum of the radiation and the production mechanism are investigated. The computational x-ray pulse shows a reasonable agreement with the experimental result. We also suggest that using alloyed wire-arrays can increase multi-keV K-shell yield by decreasing the opacity of K-shell lines. In addition, we use a detailed circuit model to study the energy coupling between the generator and the Z-pinch implosion. Recently, we are concentrating on the problems of Z-pinch driven ICF, such as dynamic hohlraum and capsule implosions. Our numerical investigations on the interaction of wire-array Z-pinches on foam convertors show qualitative agreements with experimental results on the “Qiangguang I” facility. An integrated two-dimensional simulation of dynamic hohlraum driven capsule implosion provides us the physical insights of wire-array plasma acceleration, shock generation and production, hohlraum formation, radiation ablation and fuel compression.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Heat loads to divertor nearby components from secondary radiation evolved during plasma instabilities

A fundamental issue in tokamak operation related to power exhaust during plasma instabilities is the understanding of heat and particle transport from the core plasma into the scrape-off layer and to plasma-facing materials. During abnormal and disruptive operation in tokamaks, radiation transport processes play a critical role in divertor/edge-generated plasma dynamics and are very important in determining overall lifetimes of the divertor and nearby components. This is equivalent to or greater than the effect of the direct impact of escaped core plasma on the divertor plate. We have developed and implemented comprehensive enhanced physical and numerical models in the upgraded HEIGHTS package for simulating detailed photon and particle transport in the evolved edge plasma during various instabilities. The paper describes details of a newly developed 3D Monte Carlo radiation transport model, including optimization methods of generated plasma opacities in the full range of expected photon spectra. Response of the ITER divertor's nearby surfaces due to radiation from the divertor-developed plasma was simulated by using actual full 3D reactor design and magnetic configurations. We analyzed in detail the radiation emission spectra and compared the emission of both carbon and tungsten as divertor plate materials. The integrated 3D simulation predicted unexpectedly high damage risk to the open stainless steel legs of the dome structure in the current ITER design from the intense radiation during a disruption on the tungsten divertor plate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Conditions for supersonic bent Marshak waves

Supersonic radiation diffusion approximation is an useful method to study the radiation transportation. Considering the 2-d Marshak theory, and an invariable source temperature, conditions for supersonic radiation diffusion are proved to be coincident with that for radiant flux domination in the early time when √(ε)x{sub f}/L≪1. However, they are even tighter than conditions for radiant flux domination in the late time when √(ε)x{sub f}/L≫1, and can be expressed as M>4(1+ε/3)/3 and τ>1. A large Mach number requires the high temperature, while the large optical depth requires the low temperature. Only when the source temperature is in a proper region the supersonic diffusion conditions can be satisfied. Assuming a power-low (in temperature and density) opacity and internal energy, for a given density, the supersonic diffusion regions are given theoretically. The 2-d Marshak theory is proved to be able to bound the supersonic diffusion conditions in both high and low temperature regions, however, the 1-d theory only bounds it in low temperature region. Taking SiO{sub 2} and the Au, for example, these supersonic regions are shown numerically.

36 MATERIALS SCIENCE↗

A semi-analytic model of magnetized liner inertial fusion

Presented is a semi-analytic model of magnetized liner inertial fusion (MagLIF). This model accounts for several key aspects of MagLIF, including: (1) preheat of the fuel (optionally via laser absorption); (2) pulsed-power-driven liner implosion; (3) liner compressibility with an analytic equation of state, artificial viscosity, internal magnetic pressure, and ohmic heating; (4) adiabatic compression and heating of the fuel; (5) radiative losses and fuel opacity; (6) magnetic flux compression with Nernst thermoelectric losses; (7) magnetized electron and ion thermal conduction losses; (8) end losses; (9) enhanced losses due to prescribed dopant concentrations and contaminant mix; (10) deuterium-deuterium and deuterium-tritium primary fusion reactions for arbitrary deuterium to tritium fuel ratios; and (11) magnetized α-particle fuel heating. We show that this simplified model, with its transparent and accessible physics, can be used to reproduce the general 1D behavior presented throughout the original MagLIF paper [S. A. Slutz et al., Phys. Plasmas 17, 056303 (2010)]. We also discuss some important physics insights gained as a result of developing this model, such as the dependence of radiative loss rates on the radial fraction of the fuel that is preheated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation study of 3–5 keV x-ray conversion efficiency from Ar K-shell vs. Ag L-shell targets on the National Ignition Facility laser

Tailored, high-flux, multi-keV x-ray sources are desirable for studying x-ray interactions with matter for various civilian, space and military applications. For this study, we focus on designing an efficient laser-driven non-local thermodynamic equilibrium 3–5 keV x-ray source from photon-energy-matched Ar K-shell and Ag L-shell targets at sub-critical densities (∼n{sub c}/10) to ensure supersonic, volumetric laser heating with minimal losses to kinetic energy, thermal x rays and laser-plasma instabilities. Using HYDRA, a multi-dimensional, arbitrary Lagrangian-Eulerian, radiation-hydrodynamics code, we performed a parameter study by varying initial target density and laser parameters for each material using conditions readily achievable on the National Ignition Facility (NIF) laser. We employ a model, benchmarked against Kr data collected on the NIF, that uses flux-limited Lee-More thermal conductivity and multi-group implicit Monte-Carlo photonics with non-local thermodynamic equilibrium, detailed super-configuration accounting opacities from CRETIN, an atomic-kinetics code. While the highest power laser configurations produced the largest x-ray yields, we report that the peak simulated laser to 3–5 keV x-ray conversion efficiencies of 17.7% and 36.4% for Ar and Ag, respectively, occurred at lower powers between ∼100–150 TW. For identical initial target densities and laser illumination, the Ag L-shell is observed to have ≳10× higher emissivity per ion per deposited laser energy than the Ar K-shell. Although such low-density Ag targets have not yet been demonstrated, simulations of targets fabricated using atomic layer deposition of Ag on silica aerogels (∼20% by atomic fraction) suggest similar performance to atomically pure metal foams and that either fabrication technique may be worth pursuing for an efficient 3–5 keV x-ray source on NIF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mimicking Martian dust: An in-vacuum dust deposition system for testing the ultraviolet sensors on the Curiosity rover

We have designed and developed an in-vacuum dust deposition system specifically conceived to simulate and study the effect of accumulation of Martian dust on the electronic instruments of scientific planetary exploration missions. We have used this device to characterize the dust effect on the UV sensor of the Rover Environmental Monitoring Station in the Mars science Laboratory mission of NASA in similar conditions to those found on Mars surface. The UV sensor includes six photodiodes for measuring the radiation in all UV wavelengths (direct incidence and reflected); it is placed on the body of Curiosity rover and it is severely affected by the dust deposited on it. Our experimental setup can help to estimate the duration of reliable reading of this instrument during operation. We have used an analogous of the Martian dust in chemical composition (magnetic species), color, and density, which has been characterized by X-ray spectroscopy. To ensure a Brownian motion of the dust during its fall and a homogeneous coverage on the instrumentation, the operating conditions of the vacuum vessel, determined by partial pressures and temperature, have to be modified to account for the different gravities of Mars with respect to Earth. We propose that our designed device and operational protocol can be of interest to test optoelectronic instrumentation affected by the opacity of dust, as can be the degradation of UV photodiodes in planetary exploration.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Extended application of Kohn-Sham first-principles molecular dynamics method with plane wave approximation at high energy—From cold materials to hot dense plasmas

An extended first-principles molecular dynamics (FPMD) method based on Kohn-Sham scheme is proposed to elevate the temperature limit of the FPMD method in the calculation of dense plasmas. The extended method treats the wave functions of high energy electrons as plane waves analytically and thus expands the application of the FPMD method to the region of hot dense plasmas without suffering from the formidable computational costs. In addition, the extended method inherits the high accuracy of the Kohn-Sham scheme and keeps the information of electronic structures. This gives an edge to the extended method in the calculation of mixtures of plasmas composed of heterogeneous ions, high-Z dense plasmas, lowering of ionization potentials, X-ray absorption/emission spectra, and opacities, which are of particular interest to astrophysics, inertial confinement fusion engineering, and laboratory astrophysics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical simulations of the ablative Rayleigh-Taylor instability in planar inertial-confinement-fusion targets using the FastRad3D code

The ablative Rayleigh-Taylor (RT) instability is a central issue in the performance of laser-accelerated inertial-confinement-fusion targets. Historically, the accurate numerical simulation of this instability has been a challenging task for many radiation hydrodynamics codes, particularly when it comes to capturing the ablatively stabilized region of the linear dispersion spectrum and modeling ab initio perturbations. Here, we present recent results from two-dimensional numerical simulations of the ablative RT instability in planar laser-ablated foils that were performed using the Eulerian code FastRad3D. Our study considers polystyrene, (cryogenic) deuterium-tritium, and beryllium target materials, quarter- and third-micron laser light, and low and high laser intensities. An initial single-mode surface perturbation is modeled in our simulations as a small modulation to the target mass density and the ablative RT growth-rate is calculated from the time history of areal-mass variations once the target reaches a steady-state acceleration. By performing a sequence of such simulations with different perturbation wavelengths, we generate a discrete dispersion spectrum for each of our examples and find that in all cases the linear RT growth-rate γ is well described by an expression of the form γ=α [kg/(1+ϵ kLm)]1/2−βkVa, where k is the perturbation wavenumber, g is the acceleration of the target, Lm is the minimum density scale-length, Va is the ablation velocity, and ϵ is either one or zero. The dimensionless coefficients α and β in the above formula depend on the particular target and laser parameters and are determined from two-dimensional simulation results through the use of a nonlinear curve-fitting procedure. While our findings are generally consistent with those of Betti et al. (Phys. Plasmas 5, 1446 (1998)), the ablative RT growth-rates predicted in this investigation are somewhat smaller than the values previously reported for the same target and laser parameters. It is speculated that differences in the equation-of-state and opacity models are largely responsible for the discrepancy. Resolution of this issue awaits the development of better experimental diagnostics capable of measuring small-wavelength (5–20 μm) perturbation growth due to the ablative RT instability in the linear regime.

Bates, J. W. (ORCID:0000000188087240)↗

Monochromatic backlighting of direct-drive cryogenic DT implosions on OMEGA

Backlighting is a powerful technique to observe the flow of cold and dense material in highenergy- density-plasma experiments. High-performance, direct-drive cryogenic deuterium-tritium (DT) implosions are a challenging backlighting configuration because of the low opacity of the DT shell, the high shell velocity, the small size of the stagnating shell, and the very bright selfemission of the hot core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simultaneous diagnosis of radial profiles and mix in NIF ignition-scale implosions via X-ray spectroscopy

In a National Ignition Facility implosion, hydrodynamic instabilities may cause the cold material from the imploding shell to be injected into the hot-spot (hot-spot mix), enhancing the radiative and conductive losses, which in turn may lead to a quenching of the ignition process. The bound-bound features of the spectrum emitted by high-Z ablator dopants that get mixed into the hot-spot have been previously used to infer the total amount of mixed mass; however, the typical errorbars are larger than the maximum tolerable mix. We present in this paper an improved 2D model for mix spectroscopy which can be used to retrieve information on both the amount of mixed mass and the full imploded plasma profile. By performing radiation transfer and simultaneously fitting all of the features exhibited by the spectra, we are able to constrain self-consistently the effect of the opacity of the external layers of the target on the emission, thus improving the accuracy of the inferred mixed mass. The model's predictive capabilities are first validated by fitting simulated spectra arising from fully characterized hydrodynamic simulations, and then, the model is applied to previously published experimental results, providing values of mix mass in agreement with previous estimates. Finally, we show that the new self consistent procedure leads to better constrained estimates of mix and also provides insight into the sensitivity of the hot-spot spectroscopy to the spatial properties of the imploded capsule, such as the in-flight aspect ratio of the cold fuel surrounding the hotspot.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A review on ab initio studies of static, transport, and optical properties of polystyrene under extreme conditions for inertial confinement fusion applications

Polystyrene (CH), commonly known as “plastic,” has been one of the widely used ablator materials for capsule designs in inertial confinement fusion (ICF). Knowing its precise properties under high-energy-density conditions is crucial to understanding and designing ICF implosions through radiation–hydrodynamic simulations. For this purpose, systematic ab initio studies on the static, transport, and optical properties of CH, in a wide range of density and temperature conditions (ρ= 0.1 to 100 g/cm 3 and T = 10 3 to 4 × 10 6 K), have been conducted using quantum molecular dynamics (QMD) simulations based on the density functional theory. We have built several wide-ranging, self-consistent material-properties tables for CH, such as the first-principles equation of state (FPEOS), the QMD-based thermal conductivity (Κ QMD ) and ionization, and the first-principles opacity table (FPOT). This paper is devoted to providing a review on (1) what results were obtained from these systematic ab initio studies; (2) how these self-consistent results were compared with both traditional plasma-physics models and available experiments; and (3) how these first-principles–based properties of polystyrene affect the predictions of ICF target performance, through both 1-D and 2-D radiation–hydrodynamic simulations. In the warm dense regime, our ab initio results, which can significantly differ from predictions of traditional plasma-physics models, compared favorably with experiments. When incorporated into hydrocodes for ICF simulations, these first-principles material properties of CH have produced significant differences over traditional models in predicting 1-D/2-D target performance of ICF implosions on OMEGA and direct-drive–ignition designs for the National Ignition Facility. Lastly, we will discuss the implications of these studies on the current small-margin ICF target designs using a CH ablator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Magnetic reconnection in the low solar chromosphere with a more realistic radiative cooling model

Magnetic reconnection is the most likely mechanism responsible for the high temperature events that are observed in strongly magnetized locations around the temperature minimum in the low solar chromosphere. This work improves upon our previous work [Ni et al., Astrophys. J. 852, 95 (2018)] by using a more realistic radiative cooling model computed from the OPACITY project and the CHIANTI database. We find that the rate of ionization of the neutral component of the plasma is still faster than recombination within the current sheet region. For low β plasmas, the ionized and neutral fluid flows are well-coupled throughout the reconnection region resembling the single-fluid Sweet-Parker model dynamics. Decoupling of the ion and neutral inflows appears in the higher β case with β 0 =1.46, which leads to a reconnection rate about three times faster than the rate predicted by the Sweet-Parker model. In all cases, the plasma temperature increases with time inside the current sheet, and the maximum value is above 2×10 4 K when the reconnection magnetic field strength is greater than 500 G. While the more realistic radiative cooling model does not result in qualitative changes of the characteristics of magnetic reconnection, it is necessary for studying the variations of the plasma temperature and ionization fraction inside current sheets in strongly magnetized regions of the low solar atmosphere. It is also important for studying energy conversion during the magnetic reconnection process when the hydrogen-dominated plasma approaches full ionization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Absolute Hugoniot measurements from a spherically convergent shock using x-ray radiography

Here, the canonical high pressure equation of state measurement is to induce a shock wave in the sample material and measure two mechanical properties of the shocked material or shock wave. For accurate measurements, the experiment is normally designed to generate a planar shock which is as steady as possible in space and time, and a single state is measured. A converging shock strengthens as it propagates, so a range of shock pressures is induced in a single experiment. However, equation of state measurements must then account for spatial and temporal gradients. We have used x-ray radiography of spherically converging shocks to determine states along the shock Hugoniot. The radius-time history of the shock, and thus its speed, was measured by radiographing the position of the shock front as a function of time using an x-ray streak camera. The density profile of the shock was then inferred from the x-ray transmission at each instant of time. Simultaneous measurement of the density at the shock front and the shock speed determines an absolute mechanical Hugoniot state. The density profile was reconstructed using the known, unshocked density which strongly constrains the density jump at the shock front. The radiographic configuration and streak camera behavior were treated in detail to reduce systematic errors. Measurements were performed on the Omega and National Ignition Facility lasers, using a hohlraum to induce a spatially uniform drive over the outside of a solid, spherical sample and a laser-heated thermal plasma as an x-ray source for radiography. Absolute shock Hugoniot measurements were demonstrated for carbon-containing samples of different composition and initial density, up to temperatures at which K-shell ionization reduced the opacity behind the shock. Here we present the experimental method using measurements of polystyrene as an example.

47 OTHER INSTRUMENTATION↗

Experimental considerations to observe two ionizing fronts in systems with a sharp absorption edge

This paper examines the experimental requirements to observe two shock fronts driven by a single x-ray source in systems with a sharp absorption edge. Here, we consider systems where the peak of the x-ray radiation drive coincides with the K-edge of the carbon, which occurs at a photon energy of 284 eV, causing photons to be deposited in two regions. The low-energy photons (E < 284 eV) penetrate further and drive the main shock, while the higher-energy photons (E > 284 eV) are absorbed in the ablated plasma. These higher-energy photons create an ionization front, which then produces a second shock, termed an edge-shock. Using a different radiation-hydrodynamics code and different opacity and equation of state tables, we replicate the previous work and build upon them to explore the conditions required to form the edge shock. We find that having the material K-edge coincide with the spectral domain of the radiation source is necessary but not sufficient on its own to drive the edge-shock.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Using multiple x-ray emission images of inertially confined implosions to identify spatial variations and estimate confinement volumes (invited)

We describe two methods to analyze multiple x-ray images of a small, self-emitting object, and we apply these methods to the stagnating hotspots in inertial confinement fusion experiments. The first method, the common integrated profile, can be employed to assess and quantify spatial variations in opacity. It is both a simple assessment of consistency and a sophisticated measurement of variations in a region that is otherwise difficult to observe. Second, we present a method to estimate volumes of highly asymmetric objects using multiple images of x-ray emission. The method is based on image intensities and does not require any explicit assumption of symmetry

47 OTHER INSTRUMENTATION↗

Implementing time resolved electron temperature capability at the NIF using a streak camera

A new capability on NIF has been implemented to measure the temperature of x-ray emitting sources. Although it is designed primarily for Inertial Confinement Fusion (ICF), it can be used for any hot emitting source that is well modeled. The electron temperature (T e ) of the hot spot within the core of imploded inertial confinement fusion capsules is an effective indicator of implosion performance. Currently, there are spatially and temporally integrated T e inferences using image plates. A temporally resolved measurement of Te will help elucidate the mechanisms for hot spot heating and cooling such as conduction to fuel, alpha-heating, mix and radiative losses. To determine the temporally resolved T e of hot spots, specific filters are added to an existing x-ray streak camera “Streaked Polar Instrumentation for Diagnosing Energetic Radiation (SPIDER)” to probe the emission spectrum during the x-ray burn history of implosions on the National Ignition Facility (NIF). One of the difficulties in inferring the hot spot temperature is the attenuation of the emission due to opacity from the shell and fuel. Therefore, a series of increasingly thick titanium filters were implemented to isolate the emission in specific energy regions that are sensitive to temperatures above 3keV while not significantly influenced by the shell/fuel attenuation. Additionally, a relatively thin zinc filter was used to measure the contribution of colder emission sources. Since the signal levels of the emission through the thicker filters are relatively poor, a dual slit (aperture) was designed to increase detected signal at the higher end of the spectrum. Herein, the design of the filters and slit are described, an overview of the solving technique is provided, and initial electron temperature results are reported.

47 OTHER INSTRUMENTATION↗