Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock loading”

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 235 records · Page 13

Unveiling the origins of work-hardening enhancement and mechanical instability in laser shock peened titanium

Laser shock-peening (LSP) produces complex gradients in microstructures and residual stresses, which enables fabrication of metallic engineering components with superior mechanical properties. In this work, we use non-destructive and spatially resolved high-energy X-ray diffraction (HE-XRD) to investigate the influences of LSP on microstructure, surface topography, and residual stress for high purity titanium plates. LSP is found to produce large compressive in-plane residual stresses near the peened surface which monotonically decay to zero about 2.5 mm below the surface. These properties were also tracked during in-situ tensile loading, allowing stress partitioning and work-hardening rates to be measured as a function of depth. The surface region is found to have the highest work-hardening rate and remain mechanically stable up to sample failure. In addition, a crystal rigid rotation along a transversal axis of about 27 & nbsp; was observed at large applied strain near the surface region, which is attributed to the formation of denser dispersed shear bands. Furthermore, a sudden drop in tensile stress was observed at a depth of similar to 160 mu m from the LSP surface, indicating the formation of a localized mechanical instability zone (LMIZ). The origin of this LMIZ could be attributed to the increase in shear stress induced by a load transfer from LSP surface due to the formation of denser dispersed shear bands. These observations indicate a localized transition in deformation behavior has occurred from dispersed shear banding to LMIZ and homogeneous deformation.

36 MATERIALS SCIENCE↗

Shear Strength and Release of Large Cryogenic Pellets from the Barrel of a Shattered Pellet Injector for Disruption Mitigation

Shattered pellet injection (SPI) has been chosen as the baseline disruption mitigation system on ITER due to its ability to rapidly inject material deep into the plasma to greatly increase the plasma density and radiate the thermal energy. SPI utilizes a mechanical punch or high-pressure gas to release and accelerate a pellet that has been cryogenically desublimated in the barrel of a pipe gun. Various material injection combinations could possibly be implemented during different phases of a disruption event to radiate plasma energy, reduce electromagnetic loads on machine components, avoid the formation of runaway electrons, or to dissipate runaway electrons that form. Each injection phase could possibly utilize combinations of deuterium, neon, or argon. In this paper we outline experimental measurements of pellet material shear strength at SPI operating temperatures to understand the force needed to release SPI pellets. Deuterium, neon, argon, and deuterium-neon mixture pellets with diameters of 8.5, 12.5, and 15.7 mm are formed at a range of relevant gas pressures and temperatures and dislodged from the cold zone with a slow-moving piston driven by a motor. The slow-moving piston is kept above the triple point temperature of the material while the pellet is forming, then cooled to below the triple point temperature before contacting the pellet to minimize any thermal conduction to the pellet. The piston incorporates a load cell to measure the force applied when the pellet breaks away from the cold zone in the barrel.The ability of the gas and punch methods to exceed the shear strength of the studied pellet materials for release has been analyzed. High pressure gas delivered by fast opening valves produce pressure shock to the pellet due to supersonic expansion of the propellant gas. Pressure (and therefore, force) oscillations are present due to transverse density propagation throughout the breech volume. Mechanical punches deliver an impact force through a high-kinetic energy impact. The effect of the mechanical shock on the pellet has been explored and will be presented in this paper. Scaling to larger ITER size SPI pellets will be described.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Modified Void Nucleation and Growth Model (MNAG) for Damage Evolution in BCC Ta

A void coalescence term was proposed as an addition to the original void nucleation and growth (NAG) model to accurately describe void evolution under dynamic loading. The new model, termed as modified void nucleation and growth model (MNAG model), incorporated analytic equations to explicitly account for the evolution of the void number density and the void volume fraction (damage) during void nucleation, growth, as well as the coalescence stage. The parameters in the MNAG model were fitted to molecular dynamics (MD) shock data for single-crystal and nanocrystalline Ta, and the corresponding nucleation, growth, and coalescence rates were extracted. The results suggested that void nucleation, growth, and coalescence rates were dependent on the orientation as well as grain size. Compared to other models, such as NAG, Cocks–Ashby, Tepla, and Tonks, which were only able to reproduce early or later stage damage evolution, the MNAG model was able to reproduce all stages associated with nucleation, growth, and coalescence. The MNAG model could provide the basis for hydrodynamic simulations to improve the fidelity of the damage nucleation and evolution in 3-D microstructures.

36 MATERIALS SCIENCE↗

Prediction of Probabilistic Shock Initiation Thresholds of Energetic Materials Through Evolution of Thermal-Mechanical Dissipation and Reactive Heating

The ignition threshold of an energetic material (EM) quantifies the macroscopic conditions for the onset of self-sustaining chemical reactions. The threshold is an important theoretical and practical measure of material attributes that relate to safety and reliability. Historically, the thresholds are measured experimentally. In this work, we present a new Lagrangian computational framework for establishing the probabilistic ignition thresholds of heterogeneous EM out of the evolutions of coupled mechanical-thermal-chemical processes using mesoscale simulations. Furthermore, the simulations explicitly account for microstructural heterogeneities, constituent properties, and interfacial processes and capture processes responsible for the development of material damage and the formation of hotspots in which chemical reactions initiate. The specific mechanisms tracked include viscoelasticity, viscoplasticity, fracture, post-fracture contact, frictional heating, heat conduction, reactive chemical heating, gaseous product generation, and convective heat transfer. To determine the ignition threshold, the minimum macroscopic loading required to achieve self-sustaining chemical reactions with a rate of reactive heat generation exceeding the rate of heat loss due to conduction and other dissipative mechanisms is determined. Probabilistic quantification of the processes and the thresholds are obtained via the use of statistically equivalent microstructure sample sets (SEMSS). The predictions are in agreement with available experimental data.

36 MATERIALS SCIENCE↗

Anisotropic strength behavior of single-crystal TATB

High-rate strength behavior plays an important role in the shock initiation of high explosives, with plastic deformation serving to localize heat into hot spots and as a mechanochemical means to enhance reactivity. Recent simulations predict that detonation-like shocks produce highly reactive nanoscale shear bands in the layered crystalline explosive TATB (1,3,5-triamino-2,4,6-trinitrobenzene), but the thresholds leading to this response are poorly understood. We utilize molecular dynamics to simulate the high-rate compressive stress–strain response of TATB, with a focus on understanding flow behavior. Here, the dependence of strength on pressure and loading axis (crystal orientation) is explored. The deformation mechanisms fall broadly into two categories, with compression along crystal layers activating a buckling/twinning mode and compression normal to the layers producing nanoscale shear bands. Despite the complexity of the underlying mechanisms, the crystal exhibits relatively straightforward stress–strain curves. Most of the crystal orientations studied show rapid strain softening following the onset of yielding, which settles to a steady flow state. Trajectories are analyzed using five metrics for local states and structural order, but most of these metrics yield similar distributions for these deformation mechanisms. On the other hand, a recently proposed measure of intramolecular strain energy is found to most cleanly distinguish between these mechanisms, while also providing a plausible connection with mechanochemically accelerated decomposition kinetics. Localization of intramolecular strain energy is found to depend strongly on crystal orientation and pressure.

36 MATERIALS SCIENCE↗

Computing virtual dark-field X-ray microscopy images of complex discrete dislocation structures from large-scale molecular dynamics simulations

Dark-field X-ray microscopy (DFXM) is a novel diffraction-based imaging technique that non-destructively maps the local deformation from crystalline defects in bulk materials. While studies have demonstrated that DFXM can spatially map 3D defect geometries, it is still challenging to interpret DFXM images of the high-dislocation-density systems relevant to macroscopic crystal plasticity. This work develops a scalable forward model to calculate virtual DFXM images for complex discrete dislocation structure(s) (DDS) obtained from atomistic simulations. Our new DDS-DFXM model integrates a non-singular formulation for calculating the local strain from the DDS and an efficient geometrical optics algorithm for computing the DFXM image from the strain field. We apply the model to complex DDS obtained from a large-scale mol­ecular dynamics simulation of compressive loading on single-crystal silicon. Simulated DFXM images exhibit prominent contrast for dislocation features between the multiple slip systems, demonstrating the potential of DFXM to resolve features from dislocation multiplication. In conclusion, the integrated DDS-DFXM model provides a toolbox for DFXM experimental design and image interpretation in the context of bulk crystal plasticity for a range of measurements across shock plasticity and the broader materials science community.

X-ray imaging↗

Development of pebble-based extruded carbon rods for extreme plasma heat flux environments

This work presents first experiments toward the development of continuously renewable (extrudable) pebble-based carbon rods for use as plasma-facing components in extreme steady-state plasma flux environments. The primary envisioned application of this work is a first wall that can survive long-term in future magnetic fusion power reactors while also improving recovery of the reactor fuel (tritium and deuterium atoms). Bench tests applying extreme steady-state front-surface heat loads of up to 50 MW/m 2 are presented. Continuous pebble rod front-surface recession and intact pebble recovery are successfully demonstrated, at a rate of order 0.2 cm/s. Numerical simulations of the pebble rod front-surface recession are able to match observations reasonably well, indicating that the recession mechanism can be understood as occurring due to pebble thermal expansion and resulting shock and cracking of the inter-pebble binder. Tests of the pebble rod extrusion demonstrate that friction between the rods and the stainless-steel extrusion channel is tolerably low (<50 N for the expected channel length) over a wide range of temperatures. Front-surface outgassing rates below 1000 Torr L/s/m 2 are achieved, believed to be sufficiently low for use in magnetic fusion reactors. In conclusion, initial parametric scans over pebble rod size and binder fraction to vary front-surface recession rates are presented.

36 MATERIALS SCIENCE↗

Examination of the cerium α-ε phase transition under dynamic loading with x-ray diffraction

Recent examination of the cerium Hugoniot with pyrometry and x-ray diffraction (XRD) has revealed a narrow solid-liquid two-phase region. It has been suggested that nonequilibrium melting may be occurring along the Ce Hugoniot, with either melt kinetics or a sluggish α-ε transition impeding the transition. In particular, the kinetics of the α-ε is unknown and the location of the phase boundary is in dispute. Static measurements suggest a nearly vertical phase boundary that intersects the Hugoniot at 6–7 GPa. This lies in direct conflict with dynamic measurements along the Hugoniot observing α-Ce through incipient melt. This work presents dynamic experiments using XRD to examine the behavior of the α-ε phase transition. The results show that the α-ε phase transition occurs through a tetragonal distortion, with the transition beginning at temperatures below the solid Hugoniot. Following the initial deviation from an ideal fcc structure, the c/a ratio is found to gradually increase with no steady value observed in the ε phase within the range of these experiments (below 17 GPa). Multiple diffraction patterns captured during the peak stress state show no significant change in c/a ratio prior to uniaxial release, upon which Ce reverts to an fcc structure. Furthermore, the results indicate that the α-ε transition occurs rapidly, both on loading and release. An examination of the c/a ratio with increasing temperatures suggests 11.5 GPa as a lower bound for the location of the α-ε-liquid triple point.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

X-ray imaging of liquid-liquid Mg cavitation

The coupling of dynamic loading platforms with synchrotron light sources have enabled experiments examining void formation in optically opaque materials. Much like spallation in a solid, cavitation occurs when rarefaction waves collide within a liquid. The resulting deformation leads to a complex flow field with numerous voids. It has been suggested that compaction of such a structure may be a significant source of ejecta (i.e. shallow bubble collapse). Experiments imaging supported and unsupported shocks interacting with a free surface in shock melted Mg with synchrotron radiation are presented. Here, the results are used to examine the structure of a cavitation plane that forms in the unsupported case and provide an estimate of the porosity within this region.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Novel PBX formulations containing thermally-expandable microspheres for on-demand control of explosive behavior

In this work, we present the formulation and analysis of inert Plastic-Bonded eXplosive (PBX)-surrogates loaded with thermally expandable microspheres (TEMs). TEMs consist of a thermoplastic acrylonitrile shell (10-50 µm diameter) encapsulating an inert low boiling hydrocarbon. Upon heating, the TEMs expand as the shell softens while the hydrocarbon gasifies, increasing the internal pressure and expanding the particle by as much as 120 vol%. We hypothesize that TEM expansion within a PBX will introduce tunable changes in local density and porosity, and ultimately the shock sensitivity. This paper focuses on microstructural analysis of surrogate PBX-TEM formulations, both before and after thermally- induced TEM expansion. The formulations were analyzed by scanning electron microscopy (SEM), ultra-small angle neutron scattering (USANS), and X-ray Computed Tomography (x-ray CT) under various thermal conditions. Evidence suggests that TEM expansion within a PBX is a viable method to impart microstructural changes that are predicted to have a measurable, on-demand effect on the initiation sensitivity. Future work includes development of PBX-TEM formulations with HMX and shock initiation experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Intense cyclic heating effects on thermo-fracture and thermal shock of solid tungsten and open-cell tungsten foam

Here we investigate here the effects of transient (cyclic) arc-jet plasma and laser heating on fracture behavior of W-foam and solid tungsten. The two key parameters that control the foam thermomechanical response are its density and mean cell size expressed in Pores Per Inch (PPI). Tungsten foam samples were fabricated with Chemical Vapor Deposition (CVD) with variety of PPI and relative density. These were tested under two types of qualitatively different conditions: (1) high-enthalpy arc-jet, and (2) high-power cyclic laser heating. None of the foam samples showed macroscopic through-thickness cracks. However, distributed micro-cracks were observed on ligaments and their triple junctions. Under the same loading conditions, W-foam and solid tungsten showed similar crack network pattern and characteristic length-scale. However, Crack Opening Displacement (COD) was twice as large in solid W as compared to W-foam. Foam samples that have been previously exposed to a low-pressure helium plasma showed significant changes in their surface forming nano-texture fuzz which was removed by subsequent testing in the arc-jet. Extensive fracture and re-crystallization were observed in the thin solid W disk that was fully-constrained from expansion. Thicker and fully-constrained solid W disks did not display recrystallization, grain growth, and extensive cracking. However, thicker disks that were free to expand showed some recrystallization and extensive through-thickness cracks due to less effective cooling and thus higher temperatures. Laser beam testing showed no visible damage formation at 0.19 GW/m 2 and 0.38 GW/m 2 for both low-density (23%) and high-density (43%) foams at low pulses (100-1000). Micro-cracks were observed after 10,000 pulses at 0.19 GW/m 2 in both foams, and in low-density foam after 100,000 at 0.38 GW/m 2 . The nature of thermomechanical damage in W-foam exposed to extreme power (GW/m 2 ) short-duration laser pulses was found to be qualitatively similar to that of high power (MW/m 2 ) long-duration arc-jet.

36 MATERIALS SCIENCE↗

Novel fabrication tools for dynamic compression targets with engineered voids using photolithography methods

Mesoscale imperfections, such as pores and voids, can strongly modify the properties and the mechanical response of materials under extreme conditions. Tracking the material response and microstructure evolution during void collapse is crucial for understanding its performance. In particular, imperfections in the ablator materials, such as voids, can limit the efficiency of the fusion reaction and ultimately hinder ignition. To characterize how voids influence the response of materials during dynamic loading and seed hydrodynamic instabilities, we, in this paper, have developed a tailored fabrication procedure for designer targets with voids at specific locations. Our procedure uses SU-8 as a proxy for the ablator materials and hollow silica microspheres as a proxy for voids and pores. By using photolithography to design the targets’ geometry, we demonstrate precise and highly reproducible placement of a single void within the sample, which is key for a detailed understanding of its behavior under shock compression. This fabrication technique will benefit high-repetition rate experiments at x-ray and laser facilities. Insight from shock compression experiments will provide benchmarks for the next generation of microphysics modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multifunctional Utilization of Pitch-Coated Carbon Fibers in Lithium-Based Rechargeable Batteries

This work reports carbon fibers as an electrode material and current collector for dual-carbon and lithium-ion batteries. Here, fully graphitic and semigraphitic carbon fibers undergo anion intercalation beyond 4.5 V versus Li/Li + . Symmetric dual-carbon full cells using a pitch-coated fully graphitic carbon fiber mat as both cathode and anode delivers an energy-density of 276 and 322 Wh kg -1 at 342 W kg -1 power density in the voltage range of 3.0–5.0 and 3.0–5.2 V, respectively. On the other hand, nongraphitic carbon fibers do not exhibit anion intercalation up to 5.2 V and can be used as a current collector. They also possess a larger number of Li + storage sites in their randomly oriented microstructure when used as an anode. A lithium-ion full cell with double carbon-coated C-LiFePO 4 loaded on nongraphitic carbon fiber as a cathode and pitch-coated nongraphitic carbon fiber as an anode exhibits 202.6 and 75.2 Wh kg -1 energy density at power densities of 46.75 W kg -1 and 11.7 kW kg -1 , respectively in 2.0–3.5 V range. This pitch-coated carbon fiber-based battery configuration eliminates the need for metal foils and costly fluorinated binders, lowers overall weight of the cell, and is capable of sustaining mechanical stress and thermal shock.

PF6− intercalation↗

Virtual texture analysis to investigate the deformation mechanisms in metal microstructures at the atomic scale

Understanding the deformation behavior of metallic materials at high strain rates requires the characterization of plasticity contributors such as twins, phase transformed regions, and dislocations. However, predicting the contributions from phase transformation and twinning relies on a complete understanding of the selection of variants for various loading orientations and the evolution of their volume fractions. This manuscript presents a new virtual texture (VirTex) analysis approach to characterize phase transformation and twinning variants in deformed microstructures generated using molecular dynamics (MD) simulations. Furthermore, the VirTex method involves the construction of a rotation matrix to calculate the angle/axis pairs and misorientation angles for each atom in the microstructure. Any changes in the orientation angle from angle/axis pairs and/or structure types are analyzed to determine the nucleation and evolution of variants in the microstructure. The study uses shock deformed single-crystal Fe, Ta, and Cu to analyze the variant selections for phase transformation or twinning or both in BCC and FCC systems. In addition, the VirTex analysis is able to characterize the phase transformation and twinning variants in nanocrystalline Fe and Ta microstructures. Besides characterizing variants, orientation mapping also provides an accelerated and on-the-fly approach for quantifying twin fractions in MD microstructures.

36 MATERIALS SCIENCE↗

Diffuse scattering from dynamically compressed single-crystal zirconium following the pressure-induced α → ω phase transition

The prototypical α → ω phase transition in zirconium is an ideal test bed for our understanding of polymorphism under extreme loading conditions. After half a century of study, a consensus had emerged that the transition is realized via one of two distinct displacive mechanisms, depending on the nature of the compression path. However, recent dynamic-compression experiments equipped with diffraction diagnostics performed in the past few years have revealed new transition mechanisms, demonstrating that our understanding of the underlying atomistic dynamics and transition kinetics is in fact far from complete. We present classical molecular dynamics simulations of the α → ω phase transition in single-crystal zirconium shock compressed along the [0001] axis using a machine-learning-class potential. The transition is predicted to proceed primarily via a modified version of the two-stage Usikov-Zilberstein mechanism, whereby the high-pressure ω phase heterogeneously nucleates at boundaries between grains of an intermediate β phase. We further observe the fomentation of atomistic disorder at the junctions between β grains, leading to the formation of highly defective interstitial material between the ω grains. We directly compare synthetic x-ray diffraction patterns generated from our simulations with those obtained using femtosecond diffraction in recent dynamic-compression experiments, and show that the simulations produce the same unique, anisotropic diffuse scattering signal unlike any previously seen from an elemental metal. Our simulations suggest that the diffuse signal arises from a combination of thermal diffuse scattering, nanoparticlelike scattering from residual kinetically stabilized α and β grains, and scattering from interstitial defective structures. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗

Deformation and Delamination of Dynamically Bulging Bilayered Films

We employ finite element simulations to study the dynamic bulging of films under shock tube conditions and compare these responses with quasi-static bulging. Both scenarios subject the film to time-varying pressure. Two film types are considered: uniform material films and bilayers with a relatively soft interior layer and a hard exterior layer. The interface between the materials in the bilayer cases is modeled as a cohesive zone that follows a bilinear normal traction-separation cohesive law. Additionally, bilayer films with preexisting cracks in the interior layer are analyzed. Parametric studies are also conducted to explore the effects of loading rate and specimen size. The mechanical response of all materials is assumed to be strain-rate independent, focusing solely on the inertial effects in the response of the films under dynamic loading, and the absence of such effects under quasi-static loading. The results indicate that quasi-static bulging exhibits a predictable plate bending like deformation. In contrast, depending on the loading rate and specimen size, dynamic bulging can give rise to elastic waves, and a sequence of deformation processes including initial uniform acceleration, followed by bending-like deformation, and finally through-thickness biaxial stretching. Dynamic loading also leads to larger bulge formation and greater deformation than quasi-static loading. A key finding is the role of preexisting cracks in driving delamination under dynamic loading, which is absent in quasi-static loading. This work enhances our understanding of dynamic bulging and motivates further research on characterizing delamination in layered films using dynamic bulge tests.

36 MATERIALS SCIENCE↗

Titanium alloy response sensitivity to variations in spectral reconstructions of National Ignition Facility xenon line-emission x-ray sources

Thermomechanical shock experiments on the National Ignition Facility (NIF) aim to study high strain rate dynamic material response. In such experiments, the NIF laser is used to generate high fluence x-ray emission sources, which irradiate material samples of interest. Under sufficiently high x-ray energy deposition, thermomechanical impulses are generated in the materials. While it is known that the characteristics of x-ray generated impulses vary as a function of incident x-ray spectra, it remains unclear how spectral assumptions and uncertainties in NIF spectral reconstructions affect our interpretation of impulsive loading. Here, in this paper, we simulate the response of a standard titanium alloy baseline sample to synthetic analytically derived and measured NIF xenon line-emission x-ray sources with a radiation hydrodynamics code. We vary the source spectral characteristics based on different source reconstruction techniques to understand the resulting variation in baseline sample response and compare the simulated response with experimental results. We find that the response is highly sensitive to assumptions made about the spectral contents and that knowledge of spectral uncertainties bounds our understanding of the resulting material response. The results of this effort help to extend our ability to use baseline material samples to extract quantitative properties from x-ray experiments on the NIF.

Alloys↗

Physics Guided Simulation of Electrostatic Discharge: Technical Report

Triboelectrically-charged objects may create threshold sparks, electrostatic discharge (ESD) events, to equilibrate charge between themselves and other relatively charged objects. ESD events exhibit many complex physical phenomena. They are a nexus of several fields of physics with disparate characteristic scales: plasma physics, chemical kinetics, hydrodynamics, circuit models, etc. These scales can span many orders of magnitude from the varied collisions thermalizing information within a plasma on the $\mathcal{O}(fs/ps)$ to the physical size of the plasma channel on the $\mathcal{O}(100µm)$, to the speed of a nonlinear hydrodynamic wave propagating at $\mathcal{O}(µm, ns)$. These threshold ESD events may occur in situations of programmatic importance, delivering energy and power profiles to a “victim load” generating deleterious consequences. To predict and mitigate these consequences we must answer questions about the spark’s energy budget: how much energy goes into producing the spark channel; how much gets radiated away; how much energy is advected away into the hydrodynamics; and how much energy is delivered to a victim load. An ESD simulation toolset has been created and evolved in order to answer these questions. An appropriate, physics guided implementation for simulation can be done by gaining insight into its constituent physics and leveraging that intuition to choose a suitable numerical operator. We examine in detail the chemical kinetics, circuit discharge, and hydrodynamics to deter mine dominant regimes, values, timescales, and interactions to uncover the underlying physical dynamics. We also examine and propose model reduction schemes for high-dimensional chemical kinetics. We use past and current work with experimentally validated and theoretically-verified hydrodynamics to calculate applicability limits of the non-ionizing strong shock limit. We quantify the energy budget from a hydrodynamic perspective and demonstrate that a significant fraction of the stored energy is “earmarked” for hydrodynamic advection as an energy terminus. Lastly, we combine the constituent physics of an ESD event (chemical kinetics, circuit model, and hydrodynamics) into a cohesive, actionable toolset and obtain promising results from an isothermal test case. We then propose a viable, modular evolution of the ESD toolset based upon the performed examination of the physics uncovering dominant physical scales and the stiffness of the compositional differential system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗