Search NASA⌕ Search

SEARCH · Search NASA

Results for “Equilibrium calculations”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Standardising the “Gregory method” for calculating equilibrium climate sensitivity

The equilibrium climate sensitivity (ECS) – the equilibrium global mean temperature response to a doubling of atmospheric CO 2 – is a high-profile metric for quantifying the Earth system's response to human-induced climate change. A widely applied approach to estimating the ECS is the “Gregory method” (Gregory et al., 2004), which uses an ordinary least squares (OLS) regression between the net radiative flux, N, and surface air temperature anomalies, ΔT, from a 150 year experiment in which atmospheric CO 2 concentrations are quadrupled. The ECS is determined by extrapolating the linear fit to N=0, i.e. the ΔT-intercept, indicating the point at which the system is back in equilibrium. This method has been used to compare ECS estimates across the CMIP5 and CMIP6 ensembles and will likely be a key diagnostic for CMIP7. Despite its widespread application, there is little consistency or transparency between studies in how the climate model data is processed prior to the regression, leading to potential discrepancies in ECS estimates. We identify 32 alternative data processing pathways, varying by differences in global mean weighting, net radiative flux variable, anomaly calculation method, and linear regression fit. Using 44 CMIP6 models, we systematically assess the impact of these choices on ECS estimates and calculate uncertainty ranges using two bootstrap approaches. While the inter-model ECS range is insensitive to the data processing pathway, individual outlier models exhibit notable differences. Approximating a model's native grid cell area (if irregular) with cosine of the latitude can decrease the ECS by 11 %, the choice of N-variable can change the ECS by 6 %, and some anomaly calculation methods can introduce spurious temporal correlations in the processed data. Beyond data processing choices, we also evaluate an alternative linear regression method – total least squares (TLS) – which has a more statistically robust basis than OLS. However, for consistency with previous literature, and given TLS may reduce the ECS compared to OLS (by up to 24 %), thereby making a known bias in the Gregory method worse, we do not feel there is sufficient clarity to recommend a transition to TLS in all cases. To improve reproducibility and comparability in future studies, we recommend a standardised Gregory method: weighting the global mean by cell area, using the top of the atmosphere (as opposed to the top of model) N-variable, and calculating anomalies by first applying a rolling average to the preindustrial control timeseries then subtracting from the raw CO 2 quadrupling experiment. This approach accounts for model drift while reducing noise in the data to best meet the pre-conditions of the linear regression. While CMIP6 results of the multi-model mean ECS appear insensitive to these processing choices, similar assumptions may not hold for CMIP7, underscoring the need for standardised data preparation in future climate sensitivity assessments.

Geosciences↗
DOE OSTI2025

Reduced Model of Ionization Lag in Intense Laser-Produced Plasmas

A physics-based empirical formula is derived to predict the ionization lag in underdense plasmas generated by an intense laser. Time-dependent nonlocal thermodynamic equilibrium calculations demonstrate significantly delayed ionization, due to rapid changes in plasma conditions, which critically impacts plasma properties such as opacity, emissivity, and heat transport. The reduced model, based on these calculations, enables the estimation of ionization lag without requiring in-depth knowledge of nonlocal thermodynamic equilibrium modeling. Furthermore, modeling reveals that the two-step ionization process—collisional excitation followed by photoionization—plays a crucial role in this ionization delay, with collisional excitation setting the timescale for ionization. Simulations across a range of elements, from beryllium to germanium, demonstrate that ionization lag is a widespread phenomenon, underscoring the importance of incorporating such processes into ionization modeling in radiation hydrodynamic simulations for various laser-plasma experiments.

direct drive↗
DOE OSTI2024

Steady-state properties of multi-orbital systems using quantum Monte Carlo

A precise dynamical characterization of quantum impurity models with multiple interacting orbitals is challenging. In quantum Monte Carlo methods, this is embodied by sign problems. A dynamical sign problem makes it exponentially difficult to simulate long times. A multi-orbital sign problem generally results in a prohibitive computational cost for systems with multiple impurity degrees of freedom even in static equilibrium calculations. Here, we present a numerically exact inchworm method that simultaneously alleviates both sign problems, enabling simulation of multi-orbital systems directly in the equilibrium or nonequilibrium steady-state. The method combines ideas from the recently developed steady-state inchworm Monte Carlo framework [Erpenbeck et al., Phys. Rev. Lett. 130, 186301 (2023)] with other ideas from the equilibrium multi-orbital inchworm algorithm [Eidelstein et al., Phys. Rev. Lett. 124, 206405 (2020)]. We verify our method by comparison with analytical limits and numerical results from previous methods.

Chemistry↗
DOE OSTI2024

Charge self-consistent density functional theory plus ghost rotationally invariant slave-boson theory for correlated materials

We present a charge self-consistent density functional theory combined with the ghost rotationally invariant slave-boson (DFT+gRISB) formalism for studying correlated materials. Here, this method is applied to SrVO 3 and NiO, representing prototypical correlated metals and charge-transfer insulators. For SrVO 3 , we demonstrate that DFT+gRISB yields an accurate equilibrium volume and effective mass close to experimentally observed values. Regarding NiO, DFT+gRISB enables the simultaneous description of charge-transfer and Mott-Hubbard bands, significantly enhancing the accuracy of the original DFT+RISB approach. Furthermore, the calculated equilibrium volume and spectral function reasonably agree with experimental observations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Acceleration of Thermochemistry Solves in MOOSE and Pronghorn

This work focuses on the development and implementation of strategies to accelerate thermochemical calculations within MOOSE-based multiphysics simulations, particularly for applications in MSRs. We highlight the inherent complexity of nuclear materials, which require a multiscale approach to accurately model their behavior across various physical domains, including mechanical, chemical, and thermal phenomena. Thermochemical equilibrium calculations are crucial for predicting material properties and enhancing the fidelity of these simulations. The integration of Thermochimica, a Gibbs energy minimizer, into MOOSE allows for the direct minimization of Gibbs energy at every point on the mesh. However, the computational cost of such integration is significant. To address this, we explored acceleration strategies such as multi-threading support and the use of a thermodynamic ValueCache to reduce redundant calculations. Additionally, we investigated modifications to Thermochimica to enable phase constraints and improve its coupling with phase-field models, which are essential for simulating microstructural evolution and corrosion in MSR. These efforts aim to optimize the computational efficiency and accuracy of multiphysics simulations, thereby supporting the development of reliable and efficient nuclear materials for next-generation reactor technologies.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Self-consistent equilibrium and transport simulations for NSTX-U plasmas enhanced via machine learning surrogate models

The Control-Oriented Transport SIMulator (COTSIM) is an advanced equilibrium and transport code designed for simulating tokamak discharges at computational speeds suitable for control applications. COTSIM’s modular framework enables users to select models that balance accuracy with speed according to specific needs, allowing the code to operate from fast to faster-than-real-time performance levels. This work presents recent enhancements to COTSIM’s predictive accuracy for NSTX-U scenarios, achieved by integrating neural-network-based surrogate models and self-consistent equilibrium calculations. To improve source deposition predictions, a surrogate model for NUBEAM has been incorporated. Additionally, a surrogate model for the Multi-Mode Module (MMM) now supports predictions of anomalous thermal, momentum, and particle diffusivities—key factors for modeling the evolution of temperature and rotation. Each surrogate model was specifically trained for the NSTX-U operational regime to enhance COTSIM’s accuracy while maintaining computational efficiency. Moreover, COTSIM now couples fixed-boundary equilibrium solvers with its transport solvers, enabling self-consistent predictions of plasma profiles and equilibrium evolution over the discharge. Simulation results demonstrate strong agreement between COTSIM and TRANSP predictions for NSTX-U discharges. These substantial advancements expand COTSIM’s utility in model-based control applications for NSTX-U. Potential applications include simultaneous optimization of equilibrium and transport scenarios, integration into digital twins, real-time profile estimation (e.g., temperature and rotation) from limited or noisy measurements, and advanced feedback-based scenario control.

Equilibrium and transport modeling↗
DOE OSTI2025

Catalytic hydrogenation of HMF to BHMF over copper catalysts

2,5-Bis(hydroxymethyl)furan (BHMF) is a bio-derived building block for polyester production, obtained via the hydrogenation of 5-hydroxymethylfurfural (HMF). First-principles thermodynamic equilibrium calculations indicate that this reaction is not thermodynamically limited under relevant conditions (e.g., 100 °C and high H 2 partial pressure). In this work, crude HMF was employed as the feedstock for BHMF synthesis. Initially, acidic impurities and humins were removed from unrefined HMF through filtration using a packed bed of γ-alumina. A comprehensive study of the filtration process is presented, including filtration kinetics, breakthrough curve analysis, and mathematical modeling. The purified HMF was subsequently hydrogenated over a 10 wt% CuZrO 2 catalyst, using ethanol as the reaction solvent. Batch reactions were first performed for collection of kinetic data to guide the transition to continuous flow operation. Kinetic data was collected in a fixed bed reactor at varying contact time, time on stream, temperature, and HMF concentration. This data was used to develop a kinetic model for HMF hydrogenation. Maximum BHMF production rates were achieved at 130 °C, accompanied by minor formation of byproducts from BHMF ring-opening reactions. The BHMF selectivity was 100 % at 100 °C although with lower reaction rates. Furthermore, catalyst stability tests revealed a loss of up to 50 % in catalytic activity within the first 24 h, likely due to the adsorption of HMF-derived oligomers that are not easily removed by filtration.

Crude HMF filtration↗
DOE OSTI2026

Assessing the Accuracy of Property Model Predictions for Cost Optimization of Desalination Technologies

Accurate modeling of seawater thermophysical and thermodynamic properties is critical for optimizing desalination processes. This study compares three seawater property models, a Reaktoro multicomponent model, the thermophysical seawater properties library from the Massachusetts Institute of Technology, and a simplified sodium chloride model, in the context of levelized cost of water (LCOW) minimization for reverse osmosis (RO) and mechanical vapor compression systems. Process simulations and cost optimizations reveal that although all three models yield comparable LCOW and specific energy consumption (SEC) estimates under baseline conditions, deviations among their predictions increase with salinity. Relative differences in LCOW and SEC reach up to 6% and 8%, respectively. RO results show greater variability due to differences in osmotic pressure predictions, which affect pressure constraints at high recoveries. Computational performance varies substantially; specifically, Reaktoro simulations are up to 28 times slower than empirical models due to their detailed equilibrium calculations. These results suggest that empirical models offer acceptable accuracy for routine desalination process design, while Reaktoro provides advantages in scenarios requiring detailed speciation, such as scaling or pH adjustment studies. These findings underscore the importance of selecting appropriate property models based on the modeling objective of desalination applications and motivate future work integrating thermodynamic rigor with empirical efficiency.

Physical and chemical properties↗
DOE OSTI2025

Modeling phase transformations in Mn-rich disordered rocksalt cathodes with machine-learning interatomic potentials

Mn-rich disordered rocksalt (DRX) cathode materials exhibit a phase transformation from a disordered to a partially disordered spinel-like structure (δ-phase) during electrochemical cycling. Here, in this computational study, we use charge-informed molecular dynamics with a fine-tuned CHGNet foundation potential to investigate the phase transformation in LixMn 0.8 Ti 0.1 O 1.9 F 0.1 . Our results indicate that transition metal migration occurs and reorders to form the spinel-like ordering in an FCC anion framework. The transformed structure contains a higher concentration of nontransition metal (0-TM) face-sharing channels, which are known to improve Li transport kinetics. Analysis of the Mn valence distribution suggests that the appearance of tetrahedral Mn 2+ is a consequence of spinel-like ordering, rather than the trigger for cation migration as previously suggested. Calculated equilibrium intercalation voltage profiles demonstrate that the δ-phase, unlike the ordered spinel, exhibits solid-solution signatures at low voltage. A higher Li capacity is obtained than in the DRX phase. This study provides atomic insights into solid-state phase transformation and its relation to experimental electrochemistry, highlighting the potential of machine-learning interatomic potentials for understanding complex oxide materials.

Zhong, Peichen [University of California, Berkeley↗
DOE OSTI2024

Bubble Mass Transport Measurement in the Large-Scale Test Loops at Oak Ridge National Laboratory

Molten salts are complex fluids that incorporate multi-phase behavior depending on the chemistry and the physical properties of the entrained components. These components include the carrier salt, the actinide fuel, and fission, activation, and corrosion products, the concentrations of which depends on the burnup history of the salt. Radionuclide transport from the salt into the cover gas / off-gas system depends on volatility as predicted by thermochemistry, but data from the Molten Salt Reactor Experiment (MSRE) conducted in the late 1960s suggest that bubble formation and transport are also important. Anomalously high fractions of noble metals were found in the off-gas system and were attributed to transportation with parent salt aerosols and their association with rising noble gas bubbles. The prediction of such phenomena requires coupled neutronic, thermal hydraulic, and chemical equilibrium calculations, the framework of which is being developed within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program supported by the US Department of Energy (DOE) Office of Nuclear Energy (NE). However, separate effects tests and models of experiments using tools such as SAM and Thermochimica can guide model developers through important processes. Conversely, model development guides the choice of experiments and systems to provide data that are relevant for validation. This report describes several experiments that tracked gas transport in molten salts, ranging from small-scale systems to large-scale loops. Gas transportation in a molten salt, LiCl-KCl, has been studied using the shadowgraph technique. Sensors such as residual gas analysis, Raman spectroscopy, and laser-induced breakdown spectroscopy have been tested for off-gas measurements. These data were used to interpret how the gases move through the upstream salt / cover gas and the interface between them. Differential pressure measurements were able to detect individual gas bubbles as they popped at the liquid–gas interface. Salt aerosols were collected on a cascade impactor. Their formation was also observed directly via the shadowgraph method, and most of these aerosols were launched ballistically into the cover gas. Fine mists could also be observed. Convection currents through the salt were visualized and can be used to calculate the thermophysical properties of the salt itself. The apparatuses described in this report and in a previous work (McFarlane et al. 2023) have been commissioned and are available for use in making further measurements of salt/surrogate fission product behavior. Novel approaches using neutron imaging are planned for the study of fluoride salts, which cannot be contained in quartz, so shadowgraph visualization is not available. Bubble transport in convective flow and in a slow-moving salt column are planned. The mobile laser-induced breakdown spectroscopy (LIBS) system is available for several applications, including iodine capture in a molten hydroxide scrubber, H2 transport though molten salts to complement Raman analysis, and online tracking of salt aerosol generation, transport, and deposition. This report summarizes the findings from FY24 and the plans for FY25. The work completes milestone M2AT-24OR0702013 of the DOE-NE Advanced Reactor Technology, Molten Salt Reactor Campaign, DOE-NE-5.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Emissions Characterization for Ammonia Fuel Blends in an Enclosed Swirl-Stabilized Diffusion Burner

This study investigates ammonia flames using the enclosed Sydney swirl burner (ESSB), focusing on detailed global emissions measurements. Emissions were measured via Fourier-transform-infrared (FTIR) spectrometer, utilizing a heated, long-path absorption cell. Hot, wet measurements of pertinent species were collected, and concentrations were determined via lineshape fitting in conjunction with the HITRAN database. Results show that partially cracked ammonia compositions yield lower emissions compared to pure NH3/H2, and small amounts of NH3 addition to CH4-containing fuel blends exhibit high levels of NOx and CO, with measurable HCN and unburnt CH4. Equilibrium calculations suggest trade-offs between chemical timescale, heat loss, and mixing. Future work will explore emissions sampling procedures and expand analysis using chemical reactor network modelling.

ammonia combustion↗
DOE OSTI2024

An assessment of ocean alkalinity enhancement using aqueous hydroxides: kinetics, efficiency, and precipitation thresholds

Abstract. Ocean alkalinity enhancement (OAE) is a promising approach to marine carbon dioxide removal (mCDR) that leverages the large surface area and carbon storage capacity of the oceans to sequester atmospheric CO2 as dissolved bicarbonate (HCO3-). One OAE method involves the conversion of salt in seawater into aqueous alkalinity (NaOH), which is returned to the ocean. The resulting increase in seawater pH and alkalinity causes a shift in dissolved inorganic carbon (DIC) speciation toward carbonate and a decrease in the surface ocean pCO2. The shift in the pCO2 results in enhanced uptake of atmospheric CO2 by the seawater due to gas exchange. In this study, we systematically test the efficiency of CO2 uptake in seawater treated with NaOH at aquarium (15 L) and tank (6000 L) scales to establish operational boundaries for safety and efficiency in advance of scaling up to field experiments. CO2 equilibration occurred on the order of weeks to months, depending on circulation, air forcing, and air bubbling conditions within the test tanks. An increase of ∼0.7–0.9 mol DIC per mol added alkalinity (in the form of NaOH) was observed through analysis of seawater bottle samples and pH sensor data, consistent with the value expected given the values of the carbonate system equilibrium calculations for the range of salinities and temperatures tested. Mineral precipitation occurred when the bulk seawater pH exceeded 10.0 and Ωaragonite exceeded 30.0. This precipitation was dominated by Mg(OH)2 over hours to 1 d before shifting to CaCO3,aragonite precipitation. These data, combined with models of the dilution and advection of alkaline plumes, will allow the estimation of the amount of carbon dioxide removal expected from OAE pilot studies. Future experiments should better approximate field conditions including sediment interactions, biological activity, ocean circulation, air–sea gas exchange rates, and mixing zone dynamics.

Ringham, Mallory C. (ORCID:0000000348020185)↗
DOE OSTI2025

Design of novel refractory equiatomic multi-principal elemental alloys based on Mo-Nb-Ti system for Gen IV reactor applications

Excellent irradiation damage resistance demonstrated by multi-principal elemental alloys (MPEAs) has sparked significant interest among researchers, prompting exploration into their vast compositional space, to validate their suitability for nuclear applications. A combined approach of thermodynamic and empirical parameters calculations alongside CALPHAD (CALculation of PHAse Diagrams) for phase formation predictions enable high-throughput material selection for sophisticated applications like nuclear, overcoming laborious and time-consuming experiments. Key thermodynamic and empirical parameters for eight novel equiatomic MPEAs, based on seven low thermal neutron cross section refractory elements, for predicting phase formation were calculated, and equilibrium and non-equilibrium simulations in CALPHAD were employed to comprehensively model the systems. Pseudo binary phase diagram simulations showed that Zr, V or equiatomic CrV additions to the base MoNbTi alloy (MoNbTi-Zr, MoNbTi-V and MoNbTi-CrV alloys) favor the formation of isomorphous body-centered cubic (BCC) phase at high temperatures, while Cr, Al, equiatomic ZrV, or equiatomic CrAl additions (MoNbTi-Cr, MoNbTi-Al, MoNbTi-ZrV or MoNbTi-CrAl alloys) limit the solubility of them. Equilibrium CALPHAD simulations at 750 oC were consistent with XRD results on MoNbTi, MoNbTiZr and MoNbTiCr alloys, and partially for others. Notably, elemental segregation observed in the backscattered electron (BSE) scanning electron microscopy (SEM) images of the alloys was accurately simulated through non-equilibrium Scheil solidification calculations in CALPHAD, further verified by experiments. The precipitation of TiCr2 Laves phase in Cr containing MoNbTiCr and MoNbTiCrAl was accurately predicted while discrepancies were noted in MoNbTiCrV. The equilibrium simulations also provided insights into phase compositions at specific temperatures offering a pathway for tailoring the desired microstructure and properties of these systems. Empirical parameters calculations successfully predicted random solid solution in the base MoNbTi alloy, and with an exception in MoNbTiV and MoNbTiAl, predicted intermetallic precipitation in the rest, especially, Laves phase precipitation in Cr containing alloys.

36 - MATERIALS SCIENCE↗
DOE OSTI2024

Improved Fuel Cycle Capability of Griffin for Fast Reactor Applications

Griffin is a MOOSE based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory under the Department of Energy Office of Nuclear Energy Nuclear Energy Advanced Modeling and Simulation Program. This fiscal year, the fuel cycle capability has been significantly extended by improving the assembly shuffling option to allow flexible fuel reloading in the multi-cycle depletion calculation and incorporating decay between cycles. An equilibrium core calculation capability was also implemented to find an equilibrium core. Additionally, an enrichment search capability was added to determine the enrichment condition that allows a core to reach an equilibrium cycle with the end-of-cycle k-effective meeting a user-specified target value. The updated fuel management capability has been extensively tested using the three-dimensional ABTR problem with different batch schemes, exhibiting reasonable solutions in terms of manual shuffling, equilibrium cycle, and enrichment search calculations. The cross section generation workflow capability for fast reactors was further verified to produce microscopic cross sections as well as Griffin core inputs for an ring-heterogeneous configuration. Rigorous verification tests using the ABTR problem demonstrated that the RH core calculations, with cross sections and Griffin inputs generated from the cross section workflow of Griffin, produced accurate solution for fast reactor problems. In addition, an option to convert delay neutron parameter data generated from MC 2 -3 in the DLAYXS format into XML format was added to support transient calculations using MC 2 -3-generated data.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Equilibrium Isotope Fractionation with LANL Thermochemical Code Magpie

This report demonstrates calculations of equilibrium isotope fractionation with thermochemical code magpie developed and maintained under the ASC-PEM-HE program at LANL. All the necessary background is provided, and the results of our calculations are compared with literature data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Transient Modeling and Simulation of a Generic Stable Salt Reactor

A SAM system-level model of a generic stable salt reactor has been developed to investigate thermal-hydraulic behavior and safety performance under steady and transient conditions. The model integrates information generated from a reactor physics analysis using PROTEUS and PERSENT, and a computation fluid dynamics (CFD) analysis using STAR-CCM+. A loose, iterative coupling scheme between PROTEUS and SAM is implemented to calculate the equilibrium power and temperature distributions in the steady-state critical core condition. The converged steady-state model is then used in PERSENT to calculate the four reactivity feedback temperature coefficients (Doppler, fuel density, coolant density, and core radial expansion) and kinetic parameters that are needed in SAM to model the temperature feedback effects in transient simulations. Within the fully enclosed liquid fuel pins, natural convection is the dominant heat transfer mechanism. The STAR-CCM+ model of the fuel pin considers conjugate heat transfer from the liquid fuel salt to the pin cladding and external reactor coolant. The CFD results of the axial and radial temperature profiles are used to empirically determine an effective fuel salt thermal conductivity in the SAM fuel pin model so that the temperatures predicted by the SAM model match as closely as possible the CFD results. In the central region of the fuel pin, the effective thermal conductivity is as high as similar to 60 times the physical fuel salt thermal conductivity. The whole-plant SAM model is then used to simulate an unprotected station blackout transient. The results of this simulation showed that the large negative fuel axial expansion reactivity feedback reduces fission power to similar to 2.4% nominal power. The core is cooled by natural circulation, which removes heat in the core to the emergency heat removal system, and ultimately, to the ambient. However, peak fuel salt and cladding temperatures can potentially reach as high as 1500 K, albeit briefly, if the shutdown mechanism fails to operate.

stable salt reactor; transient simulations; system↗
DOE OSTI2025

Designing a validation experiment for radio frequency condensation

Abstract Theoretical studies have suggested that nonlinear effects can lead to ‘radio frequency (RF) condensation’, where an initially broad current profile can coalesce in islands when they reach sufficient width. In suitable conditions, RF condensation can ‘self-focus’ the driven current to the center of an island, improving stabilization efficiency and reducing control complexity. In unsuitable conditions, the effect can prematurely deplete the RF energy before it reaches the island center, impairing stabilization. It is predicted that the RF condensation effect can significantly impact reactor-scale tokamaks. This paper presents a set of simulations investigating the conditions under which RF condensation might be encountered in present-day tokamaks. For concreteness, the calculations use equilibrium reconstructions for two shots from DIII-D and AUG. The Current Condensation Amid Magnetic Islands (OCCAMI) simulation code has been used for this investigation. The code takes as its input a numerically specified axisymmetric EFIT equilibrium solution, and it perturbatively constructs a 3D field with an island embedded at the appropriate rational surface. In the OCCAMI code, the GENRAY code is used for ray tracing and for calculating the power deposition along a ray trajectory, and GENRAY is coupled self-consistently to a solution of the thermal diffusion equation in the island. The simulation results described in the paper illuminate the conditions required for experimental validation of the theory of RF condensation. The simulations also provide an explanation of why the effect was not noticed in experiments prior to the publication of theoretical papers on the subject.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗