Ablation Front Research Articles

Modeling ablation of laminated composites: A novel manual mesh moving finite element analysis procedure with ABAQUS

The recent improvements of the commercial, general purpose Finite Element Analysis (FEA) software have allowed them to be used for modeling ablation problems. ABAQUS for example, provides a user-subroutine UMESHMOTION along with the Arbitrary Lagrangian-Eulerian (ALE) adaptive remesh algorithm that enable the users to couple the heat transfer with the progressive surface recession (i.e., ablation). However, such a numerical capability is limited to model ablation problems when the ablation front (i.e., receding surface) is confined within a single material domain (e.g., homogenous material). For ablation problems when the ablation front proceeds from one material domain to another (e.g., laminated composite materials that consist of multiple laminate layers with different material orientations), such a numerical capability is insufficient, since the mesh is not allowed to flow from one material domain to another when the UMESHMOTION subroutine is used. In this paper, a novel computational procedure that sequentially performs the general heat transfer analysis and the general static analysis in ABAQUS, is proposed enabling the capability of ABAQUS for modeling ablation of laminated composites. The proposed procedure is verified by comparing the predictions of temperature and ablation histories of a two-dimensional isotropic panel (i.e., with single material domain) with those obtained using the traditional UMESHMOTION+ALE method. In addition, a case study of applying the proposed procedure to predict the thermal and ablation response of a laminated carbon fiber reinforced epoxy matrix (CFRP) composite panel subjected to a high-intensity and short-duration radiative heat flux is demonstrated.

New radiographic image processing tested on the simple and double-flux platform at OMEGA

Ablation fronts and shocks are two radiative/hydrodynamic features ubiquitous in inertial confinement fusion (ICF). A specially designed shock-tube experiment was tested on the OMEGA laser facility to observe that these two features evolve at once and to assess thermodynamical and radiative properties. It is a basic science experiment aimed at improving our understanding of shocked and ablated matter which is critical to ICF design. At all time, these two moving “interfaces” separate the tube into three distinct zones where matter is ablated, shocked, or unshocked. The simple-flux or double-flux experiments, respectively, one or two halfraum-plus-tube, have been thought up to observe and image these zones using x-ray and visible image diagnostic. The possibility of observing all three regions at once was instrumental in our new radiographic image processing used to remove the backlighter background otherwise detrimental to quantitative measurements. By so doing, after processing the radiographic images of the 15 shots accumulated during the 2013 and 2015 campaigns, a quantitative comparison between experiments and our radiative hydrocode simulations was made possible. One point of the principal Hugoniot of the aerogel used as a light material in the shock-tube was inferred from that comparison. Most surprisingly, rapid variations of relative-transmission in the ablated region were observed during radiographic irradiations while it remained constant in the shocked region. This effect might be attributed to the spectral distribution variability of the backlighter during the radiographic pulse. Numerically, that distribution is strongly dependent upon NLTE models and it could potentially be used as a means to discriminate among them.

Modeling of thermal response and ablation in laminated glass fiber reinforced polymer matrix composites due to lightning strike

Abstract Thermal response and ablation of laminated glass fiber reinforced polymer matrix composites subjected to lightning strike are studied. The associated nonlinear time-dependent heat transfer model includes specific features of lightning arcs observed in physical measurements such as lightning channel radius expansion, non-uniform lightning current density, and associated heat flux. Moving spatially and temporally non-uniform lightning-current-induced heat flux boundary and moving boundary due to material phase transition caused by rapid surface ablation are also included. To predict moving phase boundary in the laminated anisotropic composites, an element deletion method is developed and embedded into finite element analysis (FEA), which is performed using ABAQUS. The Umeshmotion + ALE method based on the user subroutine Umeshmotion and arbitrary Lagrangian–Eulerian (ALE) adaptive mesh technique is also used, when applicable (i.e., moving phase boundary is confined within a top layer of the composite laminate). Heat transfer analysis is performed for a non-conductive laminated glass fiber reinforced polymer matrix composite panel representing the SNL 100-00 wind turbine tip. Thermal response of the panel subjected to pulsed and continuing lightning currents at three different lightning protection levels, LPL I, LPL II, and LPL III, is studied. Temperature-dependent anisotropic thermal properties of the composite panel are included in the analysis. The FEA results include temperature distributions and ablation zone profiles. The results show the Umeshmotion + ALE method is sufficient for the pulsed lightning current at all three LPL levels since the moving phase boundary, i.e. the ablation front, is found to be confined within the top layer of the laminate. For the continuing lightning currents at all three LPL levels, the Umeshmotion + ALE method is not applicable since the moving phase boundary comes to rest at depths exceeding the thickness of the top layer of the composite laminate.

Combustion Waves and Fronts in Flows: Flames, Shocks, Detonations, Ablation Fronts and Explosion of Stars

<i>Combustion Waves and Fronts in Flows: Flames, Shocks, Detonations, Ablation Fronts and Explosion of Stars</i>

A dual-PIXE tomography setup for reconstruction of Germanium in ICF target

Inertial Confinement Fusion (ICF) is one type of fusion energy research which could initiate nuclear fusion reactions through heating and compressing thermonuclear fuel. Compared to a pure plastic target, Germanium doping into the CH ablator layer by Glow Discharge Polymer (GDP) technique can increase the ablation velocity and the standoff distance between the ablation front and laser-deposition region. During target fabrication process, quantitative doping of Ge should be accurately controlled. Particle Induced X-ray Emission Tomography (PIXE-T) can make not only quantification of the concentration, but also reconstruction of the spatial distribution of doped element. The Si (Li) detector for PIXE tomography technique had a disadvantage of low counting rate. To make up this deficiency, another detector of Si (Li) with the same configuration positioned at the opposite side with the same detective angle 135° have been implemented. Simultaneously acquired elemental maps of Ge obtained using two detectors may be different because of the X-ray absorption along the X-ray exit route in the target. In this paper, the X-ray detection efficiency is drastically improved by this dual-PIXE tomography system.

Mix and hydrodynamic instabilities on NIF

Several new platforms have been developed to experimentally measure hydrodynamic instabilities in all phases of indirect-drive, inertial confinement fusion implosions on National Ignition Facility. At the ablation front, instability growth of pre-imposed modulations was measured with a face-on, x-ray radiography platform in the linear regime using the Hydrodynamic Growth Radiography (HGR) platform. Modulation growth of "native roughness" modulations and engineering features (fill tubes and capsule support membranes) were measured in conditions relevant to layered DT implosions. A new experimental platform was developed to measure instability growth at the ablator-ice interface. In the deceleration phase of implosions, several experimental platforms were developed to measure both low-mode asymmetries and high-mode perturbations near peak compression with x-ray and nuclear techniques. In one innovative technique, the self-emission from the hot spot was enhanced with argon dopant to "self-backlight" the shell in-flight. To stabilize instability growth, new "adiabat-shaping" techniques were developed using the HGR platform and applied in layered DT implosions.

Effect of high-Z dopant on the laser-driven ablative Richtmyer–Meshkov instability

AbstractThe effects of high-Z dopant on the laser-driven ablative Richtmyer–Meshkov instability (RMI) are investigated by theoretical analysis and radiation hydrodynamics simulations. It is found that the oscillation amplitude of ablative RMI depends on the ablation velocity, the blow-off plasma velocity and the post-shock sound speed. Owing to enhancing the radiation at the plasma corona and increasing the radiation temperature at the ablation front, the high-Z dopant in plastic target can significantly increase the ablation velocity and the blow-off plasma velocity, leading to an increase in oscillation frequency and a reduction in oscillation amplitude of the ablative RMI. The high-Z dopant in plastic target is beneficial to reduce the seed of ablative Rayleigh–Taylor instability. These results are helpful for the design of direct drive inertial confinement fusion capsules.

P2 asymmetry of Au's M-band flux and its smoothing effect due to high-Z ablator dopants

X-ray drive asymmetry is one of the main seeds of low-mode implosion asymmetry that blocks further improvement of the nuclear performance of “high-foot” experiments on the National Ignition Facility [Miller et al., Nucl. Fusion 44, S228 (2004)]. More particularly, the P2 asymmetry of Au's M-band flux can also severely influence the implosion performance of ignition capsules [Li et al., Phys. Plasmas 23, 072705 (2016)]. Here we study the smoothing effect of mid- and/or high-Z dopants in ablator on Au's M-band flux asymmetries, by modeling and comparing the implosion processes of a Ge-doped ignition capsule and a Si-doped one driven by X-ray sources with P2 M-band flux asymmetry. As the results, (1) mid- or high-Z dopants absorb hard X-rays (M-band flux) and re-emit isotropically, which helps to smooth the asymmetric M-band flux arriving at the ablation front, therefore reducing the P2 asymmetries of the imploding shell and hot spot; (2) the smoothing effect of Ge-dopant is more remarkable than Si-dopant because its opacity in Au's M-band is higher than the latter's; and (3) placing the doped layer at a larger radius in ablator is more efficient. Applying this effect may not be a main measure to reduce the low-mode implosion asymmetry, but might be of significance in some critical situations such as inertial confinement fusion (ICF) experiments very near the performance cliffs of asymmetric X-ray drives.

First Measurements of Fuel-Ablator Interface Instability Growth in Inertial Confinement Fusion Implosions on the National Ignition Facility.

Direct measurements of hydrodynamic instability growth at the fuel-ablator interface in inertial confinement fusion (ICF) implosions are reported for the first time. These experiments investigate one of the degradation mechanisms behind the lower-than-expected performance of early ICF implosions on the National Ignition Facility. Face-on x-ray radiography is used to measure instability growth occurring between the deuterium-tritium fuel and the plastic ablator from well-characterized perturbations. This growth starts in two ways through separate experiments-either from a preimposed interface modulation or from ablation front feedthrough. These experiments are consistent with analytic modeling and radiation-hydrodynamic simulations, which say that a moderately unstable Atwood number and convergence effects are causing in-flight perturbation growth at the interface. The analysis suggests that feedthrough from outersurface perturbations dominates the interface perturbation growth at mode 60.

A hybrid-drive nonisobaric-ignition scheme for inertial confinement fusion

A new hybrid-drive (HD) nonisobaric ignition scheme of inertial confinement fusion (ICF) is proposed, in which a HD pressure to drive implosion dynamics increases via increasing density rather than temperature in the conventional indirect drive (ID) and direct drive (DD) approaches. In this HD (combination of ID and DD) scheme, an assembled target of a spherical hohlraum and a layered deuterium-tritium capsule inside is used. The ID lasers first drive the shock to perform a spherical symmetry implosion and produce a large-scale corona plasma. Then, the DD lasers, whose critical surface in ID corona plasma is far from the radiation ablation front, drive a supersonic electron thermal wave, which slows down to a high-pressure electron compression wave, like a snowplow, piling up the corona plasma into high density and forming a HD pressurized plateau with a large width. The HD pressure is several times the conventional ID and DD ablation pressure and launches an enhanced precursor shock and a continuous compression wave, which give rise to the HD capsule implosion dynamics in a large implosion velocity. The hydrodynamic instabilities at imploding capsule interfaces are suppressed, and the continuous HD compression wave provides main pdV work large enough to hotspot, resulting in the HD nonisobaric ignition. The ignition condition and target design based on this scheme are given theoretically and by numerical simulations. It shows that the novel scheme can significantly suppress implosion asymmetry and hydrodynamic instabilities of current isobaric hotspot ignition design, and a high-gain ICF is promising.

Development of MPS Method for Analyzing Melt Spreading Behavior and MCCI in Severe Accidents

Spreading of molten core (corium) on reactor containment vessel floor and molten corium-concrete interaction (MCCI) are important phenomena in the late phase of a severe accident for assessment of the containment integrity and managing the severe accident. The severe accident research at Waseda University has been advancing to show that simulations with moving particle semi-implicit (MPS) method (one of the particle methods) can greatly improve the analytical capability and mechanical understanding of the melt behavior in severe accidents. MPS models have been developed and verified regarding calculations of radiation and thermal field, solid-liquid phase transition, buoyancy, and temperature dependency of viscosity to simulate phenomena, such as spreading of corium, ablation of concrete by the corium, crust formation and cooling of the corium by top flooding. Validations have been conducted against experiments such as FARO L26S, ECOKATS-V1, Theofanous, and SPREAD for spreading, SURC-2, SURC-4, SWISS-1, and SWISS-2 for MCCI. These validations cover melt spreading behaviors and MCCI by mixture of molten oxides (including prototypic UO2-ZrO2), metals, and water. Generally, the analytical results show good agreement with the experiment with respect to the leading edge of spreading melt and ablation front history of concrete. The MPS results indicate that crust formation may play important roles in melt spreading and MCCI. There is a need to develop a code for two dimensional MCCI experiment simulation with MPS method as future study, which will be able to simulate anisotropic ablation of concrete.

Transient effects in unstable ablation fronts and mixing layers in HEDP

We report results obtained for two elementary unstable flow configurations relevant to high energy density physics: the ablation front instability and the Rayleigh–Taylor -instability induced mixing layer. These two flows are characterized by a transience of their perturbation dynamics. In the ablative flow case, this perturbation dynamics transience takes the form of finite-durations of successive linear-perturbation evolution phases until reaching regimes of decaying oscillations. This behaviour is observed in various regimes: weakly or strongly accelerated ablation fronts, irradiation asymmetries or initial external-surface defects, and is a result of the mean-flow unsteadiness and stretching. In the case of the Rayleigh–Taylor-instability induced mixing layer, perturbation dynamics transience manifests itself through the extinction of turbulence and mixing as the flow reaches a stable state made of two stably stratified layers of pure fluids separated by an unstratified mixing layer. A second feature, also due to compressibility, takes the form of an intense acoustic wave production, mainly localized in the heavy fluid. Finally, we point out that a systematic short-term linear-perturbation dynamics analysis should be undertaken within the framework of non-normal stability theory.

Experimental evidence of a bubble-merger regime for the Rayleigh-Taylor Instability at the ablation front

Under the Discovery Science program, the longer pulses and higher laser energies provided by the National Ignition Facility (NIF) have been harnessed to study, first time in indirect-drive, the highly nonlinear stage of the Rayleigh-Taylor Instability (RTI) at the ablation front. A planar plastic package with pre-imposed two-dimensional broadband modulations is accelerated for up to 12 ns by the x-ray drive of a gas-filled gold radiation cavity with a radiative temperature plateau at 175 eV. This extended tailored drive allows a distance traveled in excess of 1 mm for a 130 μm thick foil, a factor 3x larger than previously achieved on other laser facilities. As a consequence, we have measured the ablative RTI in transition from the weakly nonlinear stage up to the deep nonlinear stage for various initial conditions. A bubble merger regime has been observed and the ablative stabilization strength varied by changing the plastic dopant from iodine to germanium.

Indirect-drive ablative Richtmyer Meshkov node scaling

The ablation front Rayleigh Taylor hydroinstability growth dispersion curve for indirect-drive implosions has been shown to be dependent on the Richtmyer Meshkov growth during the first shock transit phase. In this paper, a simplified treatment of the first shock ablative Richtmyer-Meshkov (ARM) growth dispersion curve is used to extract differences in ablation front perturbation growth behavior as function of foot pulse shape and ablator material for comparing the merits of various ICF design option.

Hydrodynamic instability measurements in DT-layered ICF capsules using the layered-HGR platform

The first measurements of hydrodynamic instability growth at the fuel-ablator interface in an ICF implosion are reported. Previous instability measurements on the National Ignition Facility have used plastic capsules to measure ablation front Rayleigh-Taylor growth with the Hydro.-Growth Radiography (HGR) platform. These capsules substituted an additional thickness of plastic ablator material in place of the cryogenic layer of Deuterium- Tritium (DT) fuel. The present experiments are the first to include a DT ice layer, which enables measurements of the instability growth occurring at the fuel-ablator interface. Instability growth at the fuel-ablator interface is seeded differently in two independent NIF experiments. In the first case, a perturbation on the outside of the capsule feeds through and grows on the interface. Comparisons to an implosion without a fuel layer produce a measure of the fuel's modulation. In the second case, a modulation was directly machined on the inner ablator before the fuel layer was added. The measurement of growth in these two scenarios are compared to 2D rad-hydro modeling.

Three-dimensional simulations of low foot and high foot implosion experiments on the National Ignition Facility

In order to achieve the several hundred Gbar stagnation pressures necessary for inertial confinement fusion ignition, implosion experiments on the National Ignition Facility (NIF) [E. I. Moses et al., Phys. Plasmas 16, 041006 (2009)] require the compression of deuterium-tritium fuel layers by a convergence ratio as high as forty. Such high convergence implosions are subject to degradation by a range of perturbations, including the growth of small-scale defects due to hydrodynamic instabilities, as well as longer scale modulations due to radiation flux asymmetries in the enclosing hohlraum. Due to the broad range of scales involved, and also the genuinely three-dimensional (3D) character of the flow, accurately modeling NIF implosions remains at the edge of current simulation capabilities. This paper describes the current state of progress of 3D capsule-only simulations of NIF implosions aimed at accurately describing the performance of specific NIF experiments. Current simulations include the effects of hohlraum radiation asymmetries, capsule surface defects, the capsule support tent and fill tube, and use a grid resolution shown to be converged in companion two-dimensional simulations. The results of detailed simulations of low foot implosions from the National Ignition Campaign are contrasted against results for more recent high foot implosions. While the simulations suggest that low foot performance was dominated by ablation front instability growth, especially the defect seeded by the capsule support tent, high foot implosions appear to be dominated by hohlraum flux asymmetries, although the support tent still plays a significant role. For both implosion types, the simulations show reasonable, though not perfect, agreement with the data and suggest that a reliable predictive capability is developing to guide future implosions toward ignition.

Hydrodynamic instabilities and mix studies on NIF: predictions, observations, and a path forward

The goals of the Mix Campaign are to determine how mix affects performance, locate the “mix cliff”, locate the source of the mix, and develop mitigation methods that allow performance to be increased. We have used several different drive pulse shapes and capsule designs in the Mix Campaign, to understand sensitivity to drive peak power, level of coast, rise time to peak power, adiabat, and dopant level in the capsule. Ablator material mixing into the hot spot has been shown conclusively with x-ray spectroscopy. The observed neutron yield drops steeply when the hot spot mix mass becomes too large. The mix appears to be driven by ablation- front Rayleigh-Taylor instabilities. A high foot, higher adiabat drive has a more stable ablation front and has allowed the mix mass in the hot spot to be reduced significantly. Two recent high foot shots achieved neutron yields > 1015 and measured neutron yield over clean 1D simulation (YOC) > 50%, which was one of the central goals of the Mix Campaign.

Convergent ablation measurements with gas-filled rugby hohlraum on OMEGA

Convergent ablation experiments with gas-filled rugby hohlraum were performed for the first time on the OMEGA laser facility. A time resolved 1D streaked radiography of capsule implosion is acquired in the direction perpendicular to hohlraum axis, whereas a 2D gated radiography is acquired at the same time along the hohlraum axis on a x-ray framing camera. The implosion trajectory has been measured for various kinds of uniformly doped ablators, including germanium-doped and silicon-doped polymers (CH), at two different doping fraction (2% and 4% at.). Our experiments aimed also at measuring the implosion performance of laminated capsules. A laminated ablator is constituted by thin alternate layers of un-doped and doped CH. It has been previously shown in planar geometry that laminated ablators could mitigate Rayleigh Taylor growth at ablation front. Our results confirm that the implosion of a capsule constituted with a uniform or laminated ablator behaves similarly, in accordance with post-shot simulations performed with the CEA hydrocode FCI2.

Ablation driven by hot electrons in shock ignition

An analytical model for the ablation driven by hot electrons is developed. The hot electrons are assumed to carry on the totality of the absorbed laser energy. Efficient energy coupling requires to keep the critical surface sufficiently close to the ablation front. To achieve this goal for high laser intensities a short enough laser wavelength is required. Scaling laws for the ablation pressure and the other relevant magnitudes of the ablation cloud are found in terms of the laser and target parameters.

The Ignition Physics Campaign on NIF: Status and Progress

We have made significant progress in ICF implosion performance on NIF since the 2011 IFSA. Employing a 3-shock, high adiabat CH (“High-Foot”) design, total neutron yields have increased 10-fold to 6.3 x1015 (a yield of ∼ 17 kJ, which is greater than the energy invested in the DT fuel ∼ 12kJ). At that level, the yield from alpha self-heating is essentially equivalent to the compression yield, indicating that we are close to the alpha self-heating regime. Low adiabat, 4-shock High Density Carbon (HDC) capsules have been imploded in conventional gas-filled hohlraums, and employing a 6 ns, 2-shock pulse, HDC capsules were imploded in near-vacuum hohlraums with overall coupling ∼ 98%. Both the 4- and 2-shock HDC capsules had very low mix and high yield over simulated performance. Rugby holraums have demonstrated uniform x-ray drive with minimal Cross Beam Energy Transfer (CBET), and we have made good progress in measuring and modelling growth of ablation front hydro instabilities.