Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusion coefficient”

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

Optimizing energy yield of monolithic perovskite/silicon tandem solar cells in real-world Conditions: The impact of luminescent coupling

Efficient light management is key to maximizing power conversion efficiency (PCE) in monolithic perovskite/silicon tandem solar cells. Achieving peak efficiency requires closely matched current generation in all junctions, especially in integrated configurations. However, real-world conditions vary significantly due to factors such as sunlight spectrum, diffuse-to-direct sunlight ratio, angular distribution of light, subcell temperature coefficients, and ground reflection. This study introduces a comprehensive optical and device simulation to optimize perovskite/silicon tandem cells, considering experimental luminescent coupling (LC) efficiency and its dependence on working conditions, alongside variations in radiative recombination, effect of temperature on absorptivity spectra, and cloud cover. Our results show potential energy yield improvements of up to 1.4 % with LC, based on current perovskite radiative recombination records, and up to 4 % with direct bandgap materials. Although radiative recombination's dependence on excitation intensity reduces output power and requires thicker absorbers, LC compensates for these losses. LC also lowers the optimized bandgap for the perovskite top cell from 1.72 eV to 1.64–1.68 eV, or even lower in regions with redshifted irradiance. Additionally, optimization revealed that thinner silicon bottom cells require a higher perovskite top cell bandgap, impacting the balance between fabrication cost and cell stability.

14 SOLAR ENERGY↗
DOE OSTI2024

Multinuclear PFGSTE NMR description of 39 K, 23 Na, 7 Li, and 1 H specific activation energies governing diffusion in alkali nitrite solutions

While pulsed field gradient stimulated echo nuclear magnetic resonance (PFGSTE NMR) spectroscopy has found widespread use in the quantification of self-diffusivity for many NMR-active nuclei, extending this technique to uncommon nuclei with unfavorable NMR properties remains an active area of research. Potassium-39 ( 39 K) is an archetypical NMR nucleus exhibiting an unfavorable gyromagnetic ratio combined with a very low Larmor frequency. Despite these unfavorable properties, this work demonstrates that 39 K PFGSTE NMR experiments are possible in aqueous solutions of concentrated potassium nitrite. Further, analysis of the results indicates that 39 K NMR diffusometry is feasible when the nuclei exhibit spin–lattice and spin–spin relaxation coefficients on the order of 60–100 ms and 50–100 ms, respectively. The diffusivity of 39 K followed Arrhenius behavior, and comparative 23 Na, 7 Li, and 1 H PFGSTE NMR studies of equimolal sodium nitrite and lithium nitrite solutions led to correlations between the enthalpy of hydration with the activation energy governing self-diffusion of the cations and also of water. Realizing the feasibility of 39 K PFGSTE NMR spectroscopy has a widespread impact across energy sciences because potassium is a common alkali element in energy storage materials and other applications.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Integration and validation of some modules for modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures

Three modules are integrated into the built-in OpenFOAM rhoCentralFoam solver towards accurate and efficient modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures within the OpenFOAM 10.0 framework. The first module is the mixture-averaged diffusion model. The second module is the built-in OpenFOAM Lagrangian solver coupled with optimised droplet drag coefficient and convective heat transfer coefficient sub-models. The last module is a sparse stiff chemistry solver based on dynamic adaptive hybrid integration (AHI-S). The optimised droplet sub-models are first verified in correct implementation for subsequent simulations in this work. Further, they show good accuracy against experimental and analytical data in the modelling of ammonia droplet acceleration and cooling in the flowing and/or low-temperature air. The accuracy and efficiency gains related to the mixture-averaged diffusion model and the AHI-S chemistry solver are examined by simulating 1-D detonation propagation in ammonia droplet-free/laden ammoniaoxygen mixtures. Numerical results of detonation propagation speed, gaseous temperature, density, and species distributions around the induction zone show good agreement with experimental data and analytical solutions. Compared to the built-in OpenFOAM diffusion model, the mixture-averaged diffusion model provides different numerical predictions of pulsating instabilities in detonation propagation. It shows better accuracy in depicting the detonation structure within the droplet-free section attributed to improved multi-component diffusion modelling. Compared to the built-in OpenFOAM solver EulerImplicit (backward Euler), the AHI-S chemistry solver reduces the computational cost by around 50%. It achieves satisfactory accuracy in calculating detonation propagation speed within the droplet-free section with the optimal efficiency when the safety factor, β, equals 0.5.

42 ENGINEERING↗
DOE OSTI2025

Computational Design of Cost-Effective High-Entropy Thermal/Environmental Barrier Coatings

Developing cost-effective thermal/environmental barrier coatings (TEBC) requires balance among various properties including low thermal conductivity, matching coefficient of thermal expansion (CTE), high thermal stability, high fracture toughness, and high recession resistance while being affordable. Low oxygen diffusivity is desirable as it can slow down oxygen transport to reach the underlying bond coating and hence delay oxidation of the bond coating. This project aims to design low-cost high-performance TEBC based on high entropy rare earth disilicates to protect SiC-based ceramic matrix composites from chemical and thermal attack for better performance of components in the hot section of gas turbine engines. To accelerate the TEBC design, we utilize first-principles density functional theory to predict key properties including phase stability, CTE, lattice thermal conductivity, temperature-dependent elastic constants, and oxygen diffusivity. Alloying elements including Yb, Y, Er, Eu, Gd, Lu, La, and Ce are considered, and modeling prediction are compared with available experimental results.

coefficient of thermal expansion↗
DOE OSTI2025

Flux Hypothesis for Odd Transport Phenomena

Onsager's regression hypothesis makes a fundamental connection between macroscopic transport phenomena and the average relaxation of spontaneous microscopic fluctuations. This relaxation, however, is agnostic to odd transport phenomena, in which fluxes run orthogonal to the gradients driving them. To account for odd transport, we generalize the regression hypothesis, postulating that macroscopic linear constitutive laws are, on average, obeyed by microscopic fluctuations, whether they contribute to relaxation or not. From this "flux hypothesis," Green-Kubo and reciprocal relations follow, elucidating the separate roles of broken time-reversal and parity symmetries underlying various odd transport coefficients. As an application, we derive and verify the Green-Kubo relation for odd collective diffusion in chiral active matter, first in an analytically tractable model and subsequently through molecular dynamics simulations of concentrated active spinners.

Hargus, Cory↗
DOE OSTI2026

Statistical inference of anomalous thermal transport with uncertainty quantification for interpretive 2D SOL models

The critical task of inferring anomalous cross-field transport coefficients is addressed in simulations of boundary plasmas with fluid models. A workflow for parameter inference in the UEDGE fluid code is developed using Bayesian optimization with parallelized sampling and integrated uncertainty quantification. In this workflow, transport coefficients are inferred by maximizing their posterior probability distribution, which is generally multidimensional and non-Gaussian. Uncertainty quantification is integrated throughout the optimization within the Bayesian framework that combines diagnostic uncertainties and model limitations. As a concrete example, we infer the anomalous electron thermal diffusivity $\chi_\perp$ from an interpretive 2D model describing electron heat transport in the conduction-limited region with radiative power loss. The workflow is first benchmarked against synthetic data and then tested on H-, L-, and I-mode discharges to match their midplane temperature and divertor heat flux profiles. We demonstrate that the workflow efficiently infers diffusivity and its associated uncertainty, generating 2D profiles that match 1D measurements. Future efforts will focus on incorporating more complicated fluid models and analyzing transport coefficients inferred from a large database of experimental results.

Bayesian optimization↗
DOE OSTI2024

A high-order finite difference method for moving immersed domain boundaries and material interfaces

Here, we present a high-order sharp treatment of immersed moving domain boundaries and material interfaces, and apply it to the advection-diffusion equation in two and three dimensions. The spatial discretization combines dimension-split finite difference schemes with an immersed boundary treatment based on a weighted least-squares reconstruction of the solution, providing stable discretizations with up to sixth order accuracy for diffusion terms and third order accuracy for advection terms. The temporal discretization relies on a novel strategy for maintaining high-order temporal accuracy in problems with moving boundaries that minimizes implementation complexity and allows arbitrary explicit or diagonally-implicit Runge-Kutta schemes. The approach is broadly compatible with popular PDE-specialized Runge-Kutta time integrators, including low-storage, strong stability preserving, and diagonally implicit schemes. Through numerical experiments we demonstrate that the full discretization maintains high-order spatial and temporal accuracy in the presence of complex 3D geometries and for a range of boundary conditions, including Dirichlet, Neumann, and flux conditions with large jumps in coefficients.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A high-temperature Rutherford Backscattering Spectrometry apparatus for in situ material characterization

A new methodology for high-temperature Rutherford Backscattering Spectrometry (HT-RBS) has been developed to enable in situ material characterization at elevated temperatures. A 3.5 MeV proton beam penetrates a 10-µm-thick 316L stainless steel foil mounted on a graphite substrate, with backscattered signals detected using an HT-RBS system. Conventional semiconductor detectors, primarily based on silicon, suffer significant performance degradation at temperatures higher than ~ 60 °C due to increased leakage current and noise, leading to signal distortion and failure. Here, to preserve spectral quality, a 5 µm aluminum foil shields the detector from thermal radiation, allowing reliable operation up to 900 °C at the target. A rotatable shutter provides additional thermal isolation during data collection pauses. In situ measurements of areal density changes of 316L stainless steel were conducted to validate the technique, revealing consistency with the known thermal expansion coefficient. The method facilitates seamless switching between irradiation and analysis, enabling continuous studies. This approach supports in situ investigations of diffusion, void swelling, creep, and corrosion, offering a versatile tool for advanced materials research.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Assessment of integral models for non-Boussinesq lazy plumes using numerical simulations

Integral modelling of turbulent buoyant plumes is crucial for rapid predictions of plume characteristics. While the governing equations are typically derived using self-similarity and a Boussinesq approximation, these assumptions may not hold for plumes originating from finite-area sources with large density ratios. Here, this work evaluates the accuracy of integral-scale models for non-Boussinesq lazy plumes using high-fidelity numerical simulations of turbulent helium plumes. We analyse the plume kinematics by computing vertical fluxes, plume radius and radial profiles, establishing some disparities between common practice and physical accuracy. We identify how the definition of the plume radius changes the perception of the plume structure when the flow is not self-similar and derive a relationship between the flux-based and threshold-based definitions without requiring self-similarity. We then examine the plume dynamics by evaluating the source terms from the governing plume equations. Our results support neglecting diffusive and viscous effects but emphasise the importance of the mean pressure gradient, even in the self-similar regime. Two coefficients need to be modelled: the well-known entrainment coefficient and the lesser-known momentum correction coefficient, which is a correction required for the momentum equation to account for self-similar and slender approximations. The momentum correction coefficient is found to be approximately constant and slightly greater than the assumed value of 1. The standard entrainment coefficient models perform well up to a local Richardson number three times the asymptotic value but overpredict entrainment for larger Richardson numbers. We propose a correction using the known finite limit of entrainment at infinite Richardson number.

Meehan, Michael Alexander [Sandia National Laborat↗
DOE OSTI2025

MMMnet: A Neural Network Surrogate for Real-Time Transport Prediction Based on the Updated Multi-Mode Model

The Multi-Mode Model (MMM) is a physics-based anomalous transport model integrated into TRANSP for predicting electron and ion thermal transport, electron and impurity particle transport, and toroidal and poloidal momentum transport. While MMM provides valuable predictive capabilities, its computational cost, although manageable for standard simulations, is too high for real-time control applications. MMMnet, a neural network-based surrogate model, is developed to address this challenge by significantly reducing computation time while maintaining high accuracy. Trained on TRANSP simulations of DIII-D discharges, MMMnet incorporates an updated version of MMM (9.0.10) with enhanced physics, including isotopic effects, plasma shaping via effective magnetic shear, unified correlation lengths for ion-scale modes, and a new physics-based model for the electromagnetic electron temperature gradient mode. A key advancement is MMMnet’s ability to predict all six transport coefficients, providing a comprehensive representation of plasma transport dynamics. MMMnet achieves a two-order-of-magnitude speed improvement while maintaining strong correlation with MMM diffusivities, making it well-suited for real-time tokamak control and scenario optimization.

DIII-D↗
DOE OSTI2024

Sorption kinetics and stability of conventional adsorbents for mercury remediation

In-situ remediation of mercury at numerous contaminated sites worldwide is a challenging and costly endeavor due to the persistency of this contaminant. In this study, we evaluated eight commercially available sorbent media ranging from carbon-, clays- and silica-based materials (PBC– Biochar, eSorb – Sorbster, nsPAC – Powdered Activated Carbon, fsPAC – Powdered Activated Carbon with Mackinawite, F300 – Filtrasorb 300, Si-SH – Silica Thiol, eBind – RemBind, Q-Clay – Organoclay PM-199), for their effectiveness in sorbing mercury (Hg 2+ ) and mercury complexed with dissolved organic matter (Hg-DOM). Under the chosen experimental conditions of this study, results showed that in the absence of DOM, the kinetic rates of Hg 2+ sorption onto the evaluated sorbents were in the order of 0.31 min-1 (Si-SH) to 2.98 min -1 (nsPAC), whereas in the presence of DOM, the rates varied from 0.16 min-1 (F300) to 0.95 min -1 (nsPAC). The measured sorption capacity for Hg 2+ in the absence of DOM varied from 3.02 mg/g (Q-Clay) to 35.15 mg/g (Si-SH), whereas in the presence of DOM, calculated partition coefficient (KD) ranged from 69.7 mL/g (Q-Clay) to 41,510 mL/g (Si-SH). Furthermore, kinetic data suggest liquid film diffusion was the rate-limiting steps governing mercury sorption onto the studied media. Overall, the obtained study parameters (kinetics/isotherm) are particularly important in informing robust engineering designs for deployment of the vast majority of evaluated sorbents. Thus, sorbent-based strategies offer viable solutions for cost-effective cleanup of mercury at industrially contaminated sites.

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

Integrating Ultra-Coarse-Grained Protein Models into Accessible Workflows for Multiscale Molecular Dynamics

To capture protein conformational transitions using molecular dynamics (MD), several simulation resolutions covering different spatial and temporal scales are typically needed. All-atom (AA) simulations provide fine resolution, but are computationally infeasible for large systems over longer durations. Coarse-grained (CG) and ultra-coarse-grained (UCG) models have a lower resolution and computational cost while still being able to conserve essential protein features. Prior work on a Multiscale Machinelearned Modeling Infrastructure (MuMMI) combined both AA and CG simulations to study RAS-RAF protein interactions, leveraging CG models for longer time scales and using AA to investigate unusual conformations in greater detail. However, MuMMI is still resource-intensive, and this study aims to maximize exploration of the protein conformational space while reducing computational cost. In this paper, we build on prior work that integrates UCG models based on heterogeneous elastic network modeling (hENM) into the MuMMI workflow. We demonstrate that UCG models enable accurate sampling of protein conformations, focusing on simulating RAS-RAF protein interactions. Using higher-resolution CG Martini simulation data, we can automatically refine intramolecular interactions in UCG models. We present a scalable Python package that uses fluctuations observed in higher-resolution CG Martini simulations to estimate bond coefficients of the UCG model. We built novel machine learning-based backmapping methods to recover more detailed CG Martini structures from UCG structures, using diffusion models to learn the mapping between scales. Finally, we present UCG-mini-MuMMI, an accessible and less compute-intensive version of MuMMI as a resource for the scientific community. Incorporating UCG models into MD studies is applicable to a broad range of systems and proteins, and our study offers insights into the advantages and limitations of these methods.

Chemical structure↗
DOE OSTI2024

Distinguishing Surface and Bulk Reactivity: Concentration-Dependent Kinetics of Iodide Oxidation by Ozone in Microdroplets

Iodine oxidation reactions play an important role in environmental, biological, and industrial contexts. The multiphase reaction between aqueous iodide and ozone is of particular interest due to its prevalence in the marine atmosphere and unique reactivity at the air–water interface. Here, we explore the concentration dependence of the I – + O 3 reaction in levitated microdroplets under both acidic and basic conditions. To interpret the experimental kinetics, molecular simulations are used to benchmark a kinetic model, which enables insight into the reactivity of the interface, the nanometer-scale subsurface region, and the bulk interior of the droplet. For all experiments, a kinetic description of gas- and liquid-phase diffusion is critical to interpreting the results. We find that the surface dominates the iodide oxidation kinetics under concentrated and acidic conditions, with the reactive uptake coefficient approaching an upper limit of 10 –2 at pH 3. In contrast, reactions in the subsurface dominate under more dilute and alkaline conditions, with inhibition of the surface reaction at pH 12 and an uptake coefficient that is 10× smaller. In conclusion, the origin of a changing surface mechanism with pH is explored and compared to previous ozone-dependent measurements.

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

Mesoscale Linear Elastic Modeling and Homogenization of Marine Energy Composites

The design of fiber-reinforced composite (FRC)-based components for marine energy applications necessitates a fundamental understanding of material properties and the resulting geometry to predict long-term performance. In this work, we present a modeling workflow to predict linear elastic and diffusive bulk properties at the mesoscale for an idealized geometry based on knowledge of fiber and resin properties. A parametric study was performed to identify the key model input parameters that influence bulk properties. Furthermore, we demonstrate how bulk properties can be leveraged in high-fidelity image-based simulations, where imperfections in tow geometry and voids captured during X-ray computed tomography imaging are explicitly represented within the simulation. Bulk properties of interest include moduli, Poisson’s ratios, hygroscopic swelling, diffusivity, and moisture uptake, which are key parameters for characterizing FRC performance within marine environments. Modeling predictions agreed well with experimental data, except for estimating swelling coefficients, likely due to crack accumulation as a function of moisture uptake. The mesoscale modeling workflow ultimately highlights a versatile framework for understanding the influence of material and geometric properties, which can be leveraged to rapidly assess new FRC-based components.

computational mechanics↗
DOE OSTI2025

Explosive Soot Challenge (Final Report)

This project assembled a broad ensemble of modeling and experimentation tools to study the morphological and optical properties of detonation soots in explosive fireballs. A gram-scale hemispherical high explosive was studied in a low-pressure controlled environment using in-situ experimentation with diffusely illuminated visible absorption spectroscopy, particle sizing through light scattering techniques, and post-test collections with subsequent morphological analysis. Hydrocode modeling was performed to replicate the detonation flow observations, and subsequent aerosol kinetics models provided particle size distributions and extinction coefficients from the hydrocode results. Experimentally observed soot morphologies agreed with expectation from the literature - a bimodal distribution was found, brought upon by the particles growing to a size where their inertia and fluid wakes are non-negligible. The aerosol kinetics model did not replicate the observed bimodal size distribution for lack of a coagulation kernel to represent the behavior. To recover particulate optical properties, a spectrally resolved absorption spectroscopy method termed Spectral diffuse back-illuminated extinction imaging (SBI-EI) was developed and implemented on two explosive types. Inverting the absorption spectra using a Kramers-Kronig consistent method yielded the complex index of refraction for the soots produced by the explosives. This method resulted in an unrealistic index of refraction for one of the two explosives, and this is suspected to be due to the model neglecting scattering brought upon by the large particle sizes observed. In addition to the core work, three additional studies were performed in parallel. These investigated the impact of scattering on diffuse absorption spectroscopy, studied how soots oxidate and sublimate in a well-controlled shock tube, and laid the theoretical groundwork for a new collision kernel to replicate the bimodal size distribution from the observations. Summaries of these efforts are included at the end of this report.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2025

Microstructure characteristics of LPBF&HIP fabricated graded composite transition joint between ferritic steel and austenitic stainless steel

Graded composite transition joints (GCTJs) offer a promising alternative to conventional dissimilar metal welds (DMWs) by enabling smooth compositional and microstructural transitions. However, GCTJs fabricated solely through additive manufacturing (AM) face challenges such as heat accumulation, complex parameter control, and elemental segregation. In this study, we propose a novel approach that relies on AM to design a spatially graded structure in one alloy and then employs hot isostatic pressing (HIP) as a diffusion bonding method to join it with a second alloy. Here, this method combines the flexibility of AM with the powder net-shaping advantage of HIP. Specifically, a series of closely packed austenitic stainless steel 304 conical structures were printed using laser powder bed fusion (LPBF) and then combined with ferritic steel P91 powder via HIP. By using electron backscatter diffraction (EBSD), electron probe microanalysis (EPMA), and transmission electron microscopy (TEM) techniques, the microstructure characteristics of the GCTJ of 304&P91, especially the interdiffusion zone (IDZ), have been systematically investigated. The microstructure at the interface transitions from austenite-ferrite (A+F) to austenite-martensite-ferrite (A+M+F), and finally to martensite-ferrite (M+F) due to diffusion. Additionally, the diffusion width between 304 and P91 increases with the volume fraction of P91. This unique design also ensures a gradual transition in both hardness and thermal expansion coefficient from 304 to P91, thereby enabling a smooth gradient in functional properties. Overall, this study proposes a novel approach for fabricating GCTJs and contributes to advancing design concepts in the field of dissimilar metal joining.

Additively manufacturing (AM)↗
DOE OSTI2025

Self-Trapped-Exciton Radiative Recombination in β–Ga 2 O 3 : Impact of Two Concurrent Nonradiative Auger Processes

The peculiarities of radiative and nonradiative processes associated with self-trapped intrinsic eXcitons in the excited β-Ga 2 O 3 crystals are studied via time-resolved techniques of induced absorption, transient grating, and photoluminescence (PL) at room temperature. The excitation above the bandgap is produced by laser pulses with linear light polarization parallel and orthogonal in the (–201) and (001) planes. We elucidate that the nonradiative recombination rate occurring in the eXciton prevails over its radiative emission rate in a wide range of free carrier concentration composed of excited and equilibrium electrons. Hence, the nonradiative recombination has no effect on the strong anisotropy and the shape of the eXciton emission band. However, we find out that the conventional ABC model of electron effective lifetime is insufficient for explanation of the excitation dependences. Inclusion of two nonradiative Auger mechanisms in a modified ABC formula provides excellent agreement of these dependences. We conclude that the trap-assisted Auger process is in proportion to the free electron density with coefficient B = 1.1 × 10 –11 cm 3 /s and appears at low/intermediate excitation, while the triple-particle Auger process is in proportion to Δn 2 with coefficient C = 8 × 10 –30 cm 6 /s and appears at high excitation conditions. The transition between two Auger mechanisms is accompanied by a rise of the eXciton diffusivity in preferred crystallographic directions where the radiative PL intensity is maximal. The diffusion length LD in these directions can reach values ~300 nm, but, at high excitations, LD becomes limited by Auger lifetimes. These findings pave the way for the implementation of self-trapped eXcitons into specific optoelectronic devices.

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

Phase-field model of freeze casting

Directional solidification of water-based solutions has emerged as a versatile technique for templating hierarchical porous materials. However, the underlying mechanisms of pattern formation remain incompletely understood. In this work, we present a detailed derivation and analysis of a quantitative phase-field model for simulating this nonequilibrium process. The phase-field model extends the thin-interface formulation of dilute binary alloy solidification with antitrapping to incorporate the highly anisotropic energetic and kinetic properties of the partially faceted ice-water interface. This interface is faceted in the basal plane normal to the ⟨0001⟩ directions and atomically rough in other directions within the basal plane. On the basal plane, the model reproduces a linear or nonlinear relationship between the interface growth rate and the kinetic undercooling that can be linked to experimental measurements. In both cases, spontaneous parity breaking of the solidification front is observed when the preferred growth direction is aligned with the temperature gradient. This phenomenon leads to the formation of partially faceted ice lamellae that drift laterally in one of the ⟨0001⟩ directions. Here, we demonstrate that the drifting velocity of the ice lamellae is controlled by the kinetics on the basal plane and converges as the thickness of the diffuse solid-liquid interface decreases. Furthermore, we examine the effect of the form of the kinetic anisotropy, which is chosen here such that the inverse of the kinetic coefficient varies linearly from a finite value in the ⟨0001⟩ directions to zero in all other directions within the basal plane, consistent with the assumption that the interface grows in local thermodynamic equilibrium in this plane. Our results indicate that the drifting velocity of ice lamellae is not affected by the slope of this linear relation, and the radius and undercooling at the tip of an ice lamella converge at relatively small slope values. Consequently, the phase-field simulations remain quantitative with computationally tractable choices of both the interface thickness and the slope assumed in the form of the kinetic anisotropy.

Materials science↗