Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermodynamic equilibrium”

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 145 records · Page 8

Insights into the mechanisms controlling the residual corrosion rate of borosilicate glasses

Borosilicate glasses are widely used to confine high-level radioactive wastes. The lifetime of these materials could reach hundreds of thousands of years if leaching of the glass into groundwater enables the formation of a passivating gel layer. Even in this regime, the glass will never stop corroding as thermodynamic equilibrium between glass and solution cannot be achieved. Therefore, accurate predictions of glass durability including passivation, require a deep understanding of the mechanisms controlling the so-called residual rate. However, despite tremendous efforts, these mechanisms remain poorly understood. Here, focusing on the behavior of the soluble elements of the International Simple Glass (B, Na, and Ca), we show that the residual rate is controlled by the behavior of B, a glass former supposed to dissolve instantaneously when in contact with water and thus widely considered as an ideal tracer. We then demonstrate that B release is controlled by multiple processes highly dependent on the pH. At the beginning of the passivating layer formation, the hydrolysis of B-O-Si linkages is rate-limiting and has an activation energy of ~60 kJ mol –1 , a value slightly lower than that for breaking Si-O-Si linkages. Once the fraction of closed pores resulting from gel restructuring is high enough, then diffusion of both reactants (water molecules) and some products (mainly B aq , Ca aq ) through the growing gel layer becomes rate-limiting. Consequently, B and Ca accumulate in an inner layer referred to as the active zone, with potential feedback on the B-O-Si hydrolysis. A new paradigm, including B as a key element of the system, is proposed to develop a comprehensive model for the corrosion of borosilicate glass.

36 MATERIALS SCIENCE↗

Observations of grain-boundary phase transformations in an elemental metal

The theory of grain boundary (the interface between crystallites, GB) structure has a long history 1 and the concept of GBs undergoing phase transformations was proposed 50 years ago 2,3 . The underlying assumption was that multiple stable and metastable states exist for different GB orientations 4,5,6 . The terminology ‘complexion’ was recently proposed to distinguish between interfacial states that differ in any equilibrium thermodynamic property 7 . Different types of complexion and transitions between complexions have been characterized, mostly in binary or multicomponent systems 8,9,10,11,12,13,14,15,16,17,18,19 . Simulations have provided insight into the phase behaviour of interfaces and shown that GB transitions can occur in many material systems 20,21,22,23,24 . However, the direct experimental observation and transformation kinetics of GBs in an elemental metal have remained elusive. In this work we demonstrate atomic-scale GB phase coexistence and transformations at symmetric and asymmetric $[11\bar{1}]$ tilt GBs in elemental copper. Atomic-resolution imaging reveals the coexistence of two different structures at Σ19b GBs (where Σ19 is the density of coincident sites and b is a GB variant), in agreement with evolutionary GB structure search and clustering analysis 21,25,26 . We also use finite-temperature molecular dynamics simulations to explore the coexistence and transformation kinetics of these GB phases. Our results demonstrate how GB phases can be kinetically trapped, enabling atomic-scale room-temperature observations. Our work paves the way for atomic-scale in situ studies of metallic GB phase transformations, which were previously detected only indirectly 9,15,27,28,29 , through their influence on abnormal grain growth, non-Arrhenius-type diffusion or liquid metal embrittlement.

36 MATERIALS SCIENCE↗

Chiral assemblies of pinwheel superlattices on substrates

The unique topology and physics of chiral superlattices make their self-assembly from nanoparticles highly sought after yet challenging in regard to (meta)materials. In this report we show that tetrahedral gold nanoparticles can transform from a perovskite-like, low-density phase with corner-to-corner connections into pinwheel assemblies with corner-to-edge connections and denser packing. Whereas corner-sharing assemblies are achiral, pinwheel superlattices become strongly mirror asymmetric on solid substrates as demonstrated by chirality measures. Liquid-phase transmission electron microscopy and computational models show that van der Waals and electrostatic interactions between nanoparticles control thermodynamic equilibrium. Variable corner-to-edge connections among tetrahedra enable fine-tuning of chirality. The domains of the bilayer superlattices show strong chiroptical activity as identified by photon-induced near-field electron microscopy and finite-difference time-domain simulations. The simplicity and versatility of substrate-supported chiral superlattices facilitate the manufacture of metastructured coatings with unusual optical, mechanical and electronic characteristics.

42 ENGINEERING↗

Terahertz electric-field-driven dynamical multiferroicity in SrTiO 3

The emergence of collective order in matter is among the most fundamental and intriguing phenomena in physics. In recent years, the dynamical control and creation of novel ordered states of matter not accessible in thermodynamic equilibrium is receiving much attention. The theoretical concept of dynamical multiferroicity has been introduced to describe the emergence of magnetization due to time-dependent electric polarization in non-ferromagnetic materials. In simple terms, the coherent rotating motion of the ions in a crystal induces a magnetic moment along the axis of rotation. Here we provide experimental evidence of room-temperature magnetization in the archetypal paraelectric perovskite SrTiO 3 due to this mechanism. We resonantly drive the infrared-active soft phonon mode with an intense circularly polarized terahertz electric field and detect the time-resolved magneto-optical Kerr effect. A simple model, which includes two coupled nonlinear oscillators whose forces and couplings are derived with ab initio calculations using self-consistent phonon theory at a finite temperature, reproduces qualitatively our experimental observations. A quantitatively correct magnitude was obtained for the effect by also considering the phonon analogue of the reciprocal of the Einstein–de Haas effect, which is also called the Barnett effect, in which the total angular momentum from the phonon order is transferred to the electronic one. Our findings show a new path for the control of magnetism, for example, for ultrafast magnetic switches, by coherently controlling the lattice vibrations with light.

36 MATERIALS SCIENCE↗

Real-time observation of dynamic structure of liquid-vapor interface at nanometer resolution in electron irradiated sodium chloride crystals

The dynamics and structure of the liquid and vapor interface has remained elusive for decades due to the lack of an effective tool for directly visualization beyond micrometer resolution. Here, we designed a simple liquid-cell for encapsulating the liquid state of sodium for transmission electron microscopic (TEM) observation. The real-time dynamic structure of the liquid-vapor interface was imaged and videoed by TEM on the sample of electron irradiated sodium chloride (NaCl) crystals, a well-studied sample with low melting temperature and quantum super-shells of clusters. The nanometer resolution images exhibit the fine structures of the capillary waves, composed of first-time observed three zones of structures and features, i.e. flexible nanoscale fibers, nanoparticles/clusters, and a low-pressure area that sucks the nanoparticles from the liquid to the interface. Although the phenomenons were observed based on irradiated NaCl crystals, the similarities of the phenomenons to predictions suggest our real-time ovserved dynamic structure might be useful in validating long-debated theoretical models of the liquid-vapor interface, and enhancing our knowledge in understanding the non-equilibrium thermodynamics of the liquid-vapor interface to benefit future engineering designs in microfluidics.

36 MATERIALS SCIENCE↗

Pressure-induced suppression of charge density phases across the entire rare-earth tritellurides by optical spectroscopy

The rare-earth tritellurides (RTe 3 ) are a distinct class of 2D layered materials that recently gained significant attention due to hosting such quantum collective phenomena as superconductivity or charge density waves (CDWs). Many members of this van der Waals (vdW) family crystals exhibit CDW behavior at room temperature, i.e., RTe 3 compound where R = La, Ce, Pr, Nd, Sm, Gd, and Tb. Here, our systematic studies establish the CDW properties of RTe 3 when the vdW spacing/interaction strength between adjacent RTe 3 layers is engineered under extreme hydrostatic pressures. Using a non-destructive spectroscopy technique, pressure-dependent Raman studies first establish the pressure coefficients of phonon and CDW amplitude modes for a variety of RTe 3 materials, including LaTe 3 , CeTe 3 , PrTe 3 , NdTe 3 , SmTe 3 , GdTe 3 , and TbTe 3 . Results further show that the CDW phase is eventually suppressed at high pressures when the interlayer spacing is reduced and interaction strength is increased. Comparison between different RTe 3 materials shows that LaTe 3 with the largest thermodynamic equilibrium interlayer spacing (smallest chemical pressure) exhibits the most stable CDW phases at high pressures. In contrast, CDW phases in late RTe 3 systems with the largest internal chemical pressures are suppressed easily with applied pressure. Overall results provide comprehensive insights into the CDW response of the entire RTe 3 series under extreme pressures, offering an understanding of CDW formation/engineering in a unique class of vdW RTe 3 material systems.

36 MATERIALS SCIENCE↗

Multi-scale computational screening and mechanistic insights of cyclic amines as solvents for improved lignocellulosic biomass processing

A computational screening workflow for the efficient deconstruction of cellulose, lignin and hemicellulose fractions of lignocellulosic biomass using cyclic amines as solvents. Lignocellulosic biomass is a promising feedstock for production of affordable fuels and chemicals from renewable resources. Effective solubilization and subsequent deconstruction of its cellulose, hemicellulose, and lignin fractions is essential for the viability of future biorefineries. This study used quantum chemistry-based equilibrium thermodynamics methods to evaluate the potential of 650 cyclic amines to solubilize cellulose, hemicellulose, and lignin. The activity coefficients of solvent - biopolymer interactions were predicted using the COSMO-RS (COnductor-like Screening MOdel for Real Solvents) method and used to identify cyclic amines that can efficiently dissolve and extract selective fractions of biopolymers during biomass pretreatment. Among the 650 cyclic amines, 1-piperazineethanmaine was predicted to be an effective solvent for extracting all three polymers and was experimentally shown to achieve the highest lignin removal (97.1%). Non-covalent interaction, reduced density gradient and quantum chemical calculations were performed to elucidate the dissolution mechanism of lignin, cellulose and hemicellulose and gain further molecular level insights into the interactions between the cyclic amines and biomass polymers that promote efficient solubilization and extraction. These analyses indicated that 1-piperazineethanmaine and 1-methylimidazole make noncovalent van der Waals, electrostatic interactions and hydrogen bonding with lignin, leading to enhanced lignin removal, while the strong intramolecular hydrogen bonding interactions in cellulose and hemicellulose result in weaker solvent-biopolymer interactions. Overall, the computational approach provided an efficient method for identifying cyclic amines tailored for optimal biomass pretreatment and resulted in the identification of a potential new class of solvents for effective biomass pretreatment.

Kumar, Nikhil↗

Aluminosilicate colloidal gels: from the early age to the precipitation of zeolites

Aluminosilicate hydrogels are often considered to be precursors for the crystallisation of zeolites carried out under hydrothermal conditions. The preparation of mechanically homogeneous aluminosilicate gels enables the study of these materials through bulk rheology and observation of the aging dynamics until the precipitation of crystalline zeolites. The first part of this study deals with the establishment of ternary state diagrams, in order to identify the range of chemical formulations that enable preparation of single-phase homogeneous gels. Then, by studying the viscoelastic moduli during the gelation reaction, and by yielding the gel under large deformation, we propose an empirical law considering the partial order of reaction on each chemical element, to predict the gelation time according to the chemical formulation. The scaling behavior of the elastic properties of this colloidal gel shows a transition from a strong link behavior to a weak link regime. Long term aging results in the shrinkage of the gel, accompanied by syneresis of interstitial liquid at the surface. Zeolites precipitate through crystallisation by a particle attachment mechanism, when thermodynamic equilibrium is reached. The stoichiometry of the precipitated zeolites is not only consistent with the concentration of the remaining species in the supernatant but, surprisingly, it is also very close to the partial order of the reaction of the chemical elements involved in the determination of the critical gel point. This indicates a strong correlation between the morphology of the soft amorphous gel network that is formed at an early age and those of the final solid precipitated crystals.

36 MATERIALS SCIENCE↗

Bandgap engineering of SrZrS 3 chalcogenide perovskite via substitutional doping for photovoltaic applications: a first-principles DFT study

Strontium zirconium sulfide (SrZrS 3 ) has garnered significant attention for photovoltaic (PV) applications due to its excellent optoelectronic properties, high chemical and moisture stability, and non-toxicity. However, the bandgaps of both the α- and β-phases lie outside the optimum ranges for both single-junction solar cells (SJSCs) and tandem solar cells (TSCs), thereby limiting their applications in PV technologies. In this study, we employed hybrid density functional theory to engineer the band gaps of α- and β-SrZrS 3 through substitutional doping. We considered three dopants (Hf, Sn, and Ti) at the Zr site of SrZrS 3 with varying doping concentrations and investigated their effects on the structural, electronic, and optical properties of the materials. We found that Sn and Ti doping effectively lowers the band gaps of both α- and β-SrZrS 3 , whereas Hf doping increases them. For x values up to 0.25, the band gaps of the α-SrZr 1−x Sn x S 3 and α-SrZr 1−x Ti x S 3 are within the optimum range for SJSCs, and those of β-SrZr 1−x Sn x S 3 and β-SrZr 1−x Ti x S 3 lie within the optimum range for Si/perovskite as well as perovskite/perovskite TSCs. The three dopants exhibited significant effects on the optical properties of both α- and β-SrZrS 3 , including the absorption coefficient, energy-loss functions, reflectivity, and refractivity spectra. Thermodynamic stability analysis revealed that for both phases, SrZr 1−x Hf x S 3 can be synthesized via exothermic processes, whereas the formation of SrZr 1−x Ti x S 3 and SrZr1−xSn x S 3 is endothermic and hence, not thermodynamically favorable. Further analysis showed that SrZr 1−x Ti x S 3 in both α- and β-phases are stable under thermodynamic equilibrium conditions, whereas SrZr 1−x Sn x S 3 is prone to dissociation into ternary phases (SrZrS 3 and SrSnS 3 ), especially at higher doping concentrations. These results show that Ti doping is effective in tuning the band gaps of α- and β-SrZrS 3 toward the optimal values for PV applications.

14 SOLAR ENERGY↗

CF + excitation in the interstellar medium

The detection of CF + in interstellar clouds potentially allows astronomers to infer the elemental fluorine abundance and the ionization fraction in ultraviolet-illuminated molecular gas. Because local thermodynamic equilibrium (LTE) conditions are hardly fulfilled in the interstellar medium (ISM), the accurate determination of the CF + abundance requires one to model its non-LTE excitation via both radiative and collisional processes. Here, we report quantum calculations of rate coefficients for the rotational excitation of CF + in collisions with para- and ortho-H 2 (for temperatures up to 150 K). As an application, we present non-LTE excitation models that reveal population inversion in physical conditions typical of ISM photodissociation regions (PDRs). We successfully applied these models to fit the CF + emission lines previously observed toward the Orion Bar and Horsehead PDRs. The radiative transfer models achieved with these new rate coefficients allow the use of CF + as a powerful probe to study molecular clouds exposed to strong stellar radiation fields.

79 ASTRONOMY AND ASTROPHYSICS↗

The chemical composition of globular clusters in the Local Group

We present detailed chemical abundance measurements for 45 globular clusters (GCs) associated with galaxies in (and, in one case, beyond) the Local Group. The measurements are based on new high-resolution integrated-light spectra of GCs in the galaxies NGC 185, NGC 205, M 31, M 33, and NGC 2403, combined with reanalysis of previously published observations of GCs in the Fornax dSph, WLM, NGC 147, NGC 6822, and the Milky Way. The GCs cover the range -2.8 < [Fe/H] < -0.1 and we determined abundances for Fe, Na, Mg, Si, Ca, Sc, Ti, Cr, Mn, Ni, Cu, Zn, Zr, Ba, and Eu. Corrections for non local thermodynamic equilibrium effects are included for Na, Mg, Ca, Ti, Mn, Fe, Ni, and Ba, building on a recently developed procedure. For several of the galaxies, our measurements provide the first quantitative constraints on the detailed composition of their metal-poor stellar populations. Overall, the GCs in different galaxies exhibit remarkably uniform abundance patterns of the α, iron-peak, and neutron-capture elements, with a dispersion of less than 0.1 dex in [α/Fe] for the full sample. There is a hint that GCs in dwarf galaxies are slightly less α-enhanced (by ~0.04 dex on average) than those in larger galaxies. One GC in M 33 (HM33-B) resembles the most metal-rich GCs in the Fornax dSph (Fornax 4) and NGC 6822 (SC7) by having α-element abundances closer to scaled-solar values, possibly hinting at an accretion origin. A principal components analysis shows that the α-element abundances strongly correlate with those of Na, Sc, Ni, and Zn. Several GCs with [Fe/H] < -1.5 are deficient in Mg compared to other α-elements. We find no GCs with strongly enhanced r-process abundances as reported for metal-poor stars in some ultra-faint dwarfs and the Magellanic Clouds. The similarity of the abundance patterns for metal-poor GCs in different environments points to similar early enrichment histories and only allow for minor variations in the initial mass function.

79 ASTRONOMY AND ASTROPHYSICS↗

Rotational state-to-state transition rate coefficients for H 2 O + H 2 O collisions at nonequilibrium conditions

Aims.The goal is to develop a database of rate coefficients for rotational state-to-state transitions in H 2 O + H 2 O collisions that is suitable for the modeling of energy transfer in nonequilibrium conditions, in which the distribution of rotational states of H 2 O deviates from local thermodynamic equilibrium. Methods.A two-temperature model was employed that assumed that although there is no equilibrium between all possible degrees of freedom in the system, the translational and rotational degrees of freedom can be expected to achieve their own equilibria independently, and that they can be approximately characterized by Boltzmann distributions at two different temperatures,T kin andT rot . Results.Upon introducing our new parameterization of the collisional rates, taking into account their dependence on bothT kin andT rot , we find a change of up to 20% in the H 2 O rotational level populations for both ortho and para-H 2 O for the part of the cometary coma where the nonequilibrium regime occurs.

Astronomy & Astrophysics↗

The R-Process Alliance: Hunting for gold in the near-UV spectrum of 2MASS J05383296–5904280

Over the past few years, the R-Process Alliance (RPA) has successfully carried out a search for stars that are highly enhanced in elements produced via the rapid neutron-capture (r-) process. In particular, the RPA has identified a number of relatively bright, highly r-process-enhanced (r-II) stars, suitable for observations with the Hubble Space Telescope (HST), facilitating abundance derivation of elements such as gold (Au) and cadmium (Cd). This paper presents the detailed abundances derived for the metal-poor ([Fe/H] = −2.55) highly r-process-enhanced ([Eu/Fe] = +1.29) r-II star 2MASS J05383296–5904280. One-dimensional local thermodynamic equilibrium (LTE) elemental abundances were derived via equivalent width and spectral synthesis using high-resolution high signal-to-noise near-UV HST/STIS and optical Magellan/MIKE spectra. Abundances were determined for 43 elements, including 26 neutron-capture elements. In particular, abundances of the rarely studied elements Nb, Mo, Cd, Lu, Os, Pt, and Au are derived from the HST spectrum. These results, combined with RPA near-UV observations of two additional r-II stars, increase the number of Cd abundances derived for r-process-enriched stars from seven to ten and Au abundances from four to seven. A large star-to-star scatter is detected for both of these elements, highlighting the need for more detections enabling further investigations, specifically into possible non-LTE effects.

chemically peculiar stars↗

The R -Process Alliance: Exploring the cosmic scatter among ten r -process sites with stellar abundances

Context. The astrophysical origin of the rapid neutron-capture process (r-process), responsible for producing roughly half of the elements heavier than iron, remains uncertain. Detailed chemical signatures from the oldest, most metal-poor stars, which act as fossil records of the earliest nucleosynthesis events, can be used to identify the dominant r-process sites. Aims. We present a homogeneous chemical abundance analysis of ten r-process element-enhanced stars. These old and metal-poor stars are strongly enriched in r-process elements with minimal contamination from other nucleosynthetic sources. By focusing on this chemically pure sample, we aim to investigate intrinsic variations in the r-process abundance patterns and explore their implications for the nature and potential diversity of r-process sites. Methods. We performed a detailed chemical abundance analysis of high-resolution, high-signal-to-noise spectra. For each star, we inspected over 1400 individual absorption lines using a combination of equivalent width measurements and spectral synthesis. The analysis was conducted under the assumption of 1D local thermodynamic equilibrium and employing the MOOG radiative transfer code. Results. We derived abundances for 54 chemical species, including 29 neutron-capture (n-capture) elements, covering the full mass range of the r-process abundance pattern. A kinematic analysis reveals that stars likely originated from ten kinematically distinct systems. Based on this assumption, we used the sample to probe the maximum variation expected from ten independent r-process nucleosynthesis events and computed the intrinsic dispersion of each element relative to Zr and Eu for the light and heavy r-process elements, respectively. This exercise resulted in a remarkably low cosmic scatter across the ten r-process sites enriching these stars; for the rare earth and third peak elements, for example, we find σ [La/Eu] = 0.08 and σ [Os/Eu] = 0.11 dex, while the scatter between light and heavy elements, σ [Zr/Eu] , is slightly higher at 0.18 dex. Conclusions. The elemental abundance patterns across the ten independent r-process sites show remarkably small cosmic dispersions. This minimal dispersion suggests a high degree of uniformity in r-process yields across diverse astrophysical environments.

Astronomy and AstroPhysics↗

Modeling radiative-shocks created by laser–cluster interactions

Radiative-shocks induced by laser–cluster interactions are modeled using radiation-hydrodynamic simulations. A good agreement—in both shock velocity and density profiles—is obtained between experiment and simulations, indicating that non-local thermodynamic equilibrium (NLTE) radiative effects are important in the experimental regime examined, particularly at early times (≤30 ns) due to the elevated temperatures (≥35 eV). The enhanced NLTE radiative emission causes the shock to be reduced in amplitude, increased in width, and reduced in propagation velocity, while the amplitude of the radiative precursor is increased. As the density and temperature conditions are relatively modest, this potentially has important implications for the scalings that are used in laboratory–astrophysics to transform between laboratory and astrophysical scales, which do not hold for non-LTE systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Deep learning for NLTE spectral opacities

Computer simulations of high energy density science experiments are computationally challenging, consisting of multiple physics calculations including radiation transport, hydrodynamics, atomic physics, nuclear reactions, laser–plasma interactions, and more. To simulate inertial confinement fusion (ICF) experiments at high fidelity, each of these physics calculations should be as detailed as possible. However, this quickly becomes too computationally expensive even for modern supercomputers, and thus many simplifying assumptions are made to reduce the required computational time. Much of the research has focused on acceleration techniques for the various packages in multiphysics codes. In this work, we explore a novel method for accelerating physics packages via machine learning. The non-local thermodynamic equilibrium (NLTE) package is one of the most expensive calculations in the simulations of indirect drive inertial confinement fusion, taking several tens of percent of the total wall clock time. We explore the use of machine learning to accelerate this package, by essentially replacing the physics calculation with a deep neural network that has been trained to emulate the physics code. Overall, we demonstrate the feasibility of this approach on a simple problem and perform a side-by-side comparison of the physics calculation and the neural network inline in an ICF Hohlraum simulation. We show that the neural network achieves a 10× speed up in NLTE computational time while achieving good agreement with the physics code for several quantities of interest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multiphase phase-field approach for solid–solid phase transformations via propagating interfacial phase in HMX

In this study, a thermodynamically consistent multiphase phase-field theory has been formulated to describe temperature-induced solid–solid phase transformations via an interfacial phase. This theory satisfies all thermodynamic equilibrium and stability conditions. It allows us to formulate the analytical solutions of the interface profile, energy, width, and velocity for each of the propagating solid–melt and solid–solid interfaces. The formation of the intermediate melt during solid–solid phase transformation has been studied hundreds of degrees below melting temperature in the HMX energetic crystal. The effect of the penalizing term in the phase-field model and two nanoscale material parameters on the appearance and degree of disordering of intermediate melt has been analyzed for the non-equilibrium interface. The developed approach is applicable to phase transitions in geological, pharmaceutical, ferroelectric, colloidal, and superhard materials where phase transformations occur via the intermediate phase.

Roy, Arunabha M. (ORCID:000000034790726X)↗

Near-field models and simulations of pellet ablation in tokamaks

Here, detailed numerical studies of the ablation of a single neon pellet in the plasma disruption mitigation parameter space have been performed. Simulations were carried out using FronTier, a hydrodynamic/low magnetic Reynolds number MHD code with explicit tracking of material interfaces. FronTier's physics models resolve the pellet surface ablation and the formation of a dense, cold cloud of ablated material, the deposition of energy from hot plasma electrons, the expansion of the ablation cloud along magnetic field lines, and the radiation losses. A local thermodynamic equilibrium model based on Saha equations has been used to resolve atomic processes in the cloud and Redlich–Kwong corrections to the ideal gas equation of state for cold and dense gases, which have been used near the pellet surface. The FronTier pellet code is a next generation of the code described. It has been validated against the semi-analytic improved neutral gas shielding model in the 1D spherically symmetric approximation. Main results include quantification of the influence of atomic processes and Redlich–Kwong corrections on the pellet ablation in spherically symmetric approximation and verification of analytic scaling laws in a broad range of pellet and plasma parameters. Using axially symmetric MHD simulations, the properties of ablation channels and the reduction of pellet ablation rates in magnetic fields of increasing strength have been studied. While the main emphasis has been given to neon pellets for the plasma disruption mitigation, selected results on deuterium fueling pellets have also been presented.

97 MATHEMATICS AND COMPUTING↗