Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dilute solution”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Glycolate Analysis in Tank 22: Developing and Testing Analytical Methods for the Savannah River Site Liquid Waste System

Researchers developed and tested a range of techniques to support a low mg/L Limit of Quantitation (LOQ) of glycolate in radioactive Tank 22 waste solution by Ion Chromatography (IC) and by proton NMR. For the IC method, Dionex OnGuard II cartridges were tested as a means of removing alkali earth and transition metals that can interfere during ion chromatography applications especially when analyzing for low-concentration, chelating analytes. Additionally, concentrations of nitrate in the raw Tank 22 sample (5,000 to 10,000 mg/L) were managed by using reasonable levels of sample dilution. The resulting IC performance quality was enhanced by improving the baseline, increasing sensitivity [Limit of Quantitation (LOQ) of 12 mg/L and Limit of Detection (LOD) of 4.0 mg/L], and resolving analytes into well-defined Gaussian peaks when using the Dionex OnGuard II H⁺ cartridges to remove matrix interferences. High concentrations of nitrate limited the performance of the IC method. To achieve an acceptable baseline, samples with high nitrate content require more dilution, resulting in higher detection limits. Tests on samples with higher concentrations of nitrate (Tank 30 and 32 supernate with approximately 150000 mg/L nitrate) suggested that IC detection limits for glycolate in these samples would be > 500 mg/L. Thus, an LOQ of 12 mg/L is not feasible by this IC method on evaporator feed samples. An alternative method of glycolate analysis, using proton nuclear magnetic resonance (H NMR), was developed by the research team. In initial tests, the H NMR technique provided reasonable quantitation of glycolate in Tank 22 conditions by direct observation of the liquid. The H NMR method may provide improved detection limits for solutions with higher nitrate concentrations. We recommend further development of this analysis for high nitrate LWS samples such as evaporator feed and evaporator drop tank content. A method for pretreatment of samples using crystalline silicotitanate (CST) was developed and tested. The objective of the pretreatment was to facilitate analysis of samples with higher levels of radioactivity. The pretreatment did not influence the glycolate concentration in solution, and the pretreatment is projected to reduce Cs 137 activity in a sample by a factor of 16,200. The double strike CST pretreatment was demonstrated and would allow milliliters of higher activity samples to be transferred from the Shielded Cells and handled in a containment unit for glycolate analysis. The various studies validated IC and H NMR methods for glycolate analysis, defined the range of applicability, and demonstrated key supporting analytical protocols. Based on the results, high quality glycolate analysis of Tank 22 is feasible down to approximately 12 mg/L, with the potential for broader applicability of the methods to other conditions in the Savannah River Site Liquid Waste System (LWS).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Measurement of Plutonium Ratios for CRM 137A Recertification

This report documents the analytical methodologies used by the Chemical and Isotopics Mass Spectrometry (CIMS) Group at the US Department of Energy’s Oak Ridge National Laboratory (ORNL) for plutonium isotopic ratio characterizations on purified ~1 mg units prepared by Lawrence Livermore National Laboratory (LLNL) from a mother solution created from an existing 250 mg unit of certified reference material (CRM) 137. The analytical processes described are designed to maintain traceability of the final certified attributes and to eliminate sources of systematic bias for the measurements described. The measurands for certification are the plutonium isotope ratios: 238 Pu/ 239 Pu, 240 Pu/ 239 Pu, 241 Pu/ 239 Pu, and 242 Pu/ 239 Pu. The preparation and analytical process included the following activities: Preparation of reference materials to be used as calibrants to establish traceability to the SI unit of mass for the isotopic analyses and quality controls. For the plutonium isotopic analyses, CRM 128 (a CRM comprising equal atoms of 239 Pu and 242 Pu) was used as the primary calibrant; Each CRM 137A unit was reconstituted using dilute high-purity HNO3. From each, seven replicates were distributed for the isotopic characterizations by thermal ionization mass spectrometry (TIMS) using the total evaporation (TE) method for isotope ratio measurements and MC-ICP-MS using a dynamic peak hooping method for secondary isotope ratio measurements; Due to time delay in the analyses from the last purification of the original CRM 137 mother solution by LLNL, an aliquot of the three provided plutonium units was purified using Eichrom TEVA separation methods to remove americium in-growth and uranium impurities; Minimally corrected data were reported from each analysis using reporting templates provided by the New Brunswick Laboratory Program Office (NBL PO) for the specific methods used for isotopic measurements. In addition, a copy of the instrument exports, Excel templates for isotopic data corrections, and the uncertainty budgets completed with the GUM Workbench software package were submitted to NBL PO.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2025

Diffusion behavior of lanthanide fission products in bcc Fe cladding: A first-principles study

Fuel-cladding chemical interaction poses significant challenges in nuclear reactors, where fission products generated from nuclear fuel interact with Fe-based cladding materials, potentially compromising their structural integrity. This study investigates the diffusion behavior of lanthanide fission products, Lanthanum (La), Cerium (Ce), Praseodymium (Pr), and Neodymium (Nd), within body-centered cubic (bcc) Fe cladding using the density functional theory, nudged elastic band method, and self-consistent mean field theory. Our results reveal significant vacancy binding energies, particularly with the 1st and 2nd nearest neighbors, which diminish beyond the 5th nearest neighbor, with La exhibiting the strongest binding affinity, followed by Nd, Ce, and Pr. The nudged elastic band calculations indicate significant high barriers for the dissociation of 1st nearest neighbor vacancy-solute pairs for all fission products. The tracer diffusion coefficients of these fission products were derived in an Arrhenius form, with a magnetic correction that accounts for the high-temperature paramagnetic state. The significant trapping effect of vacancies caused by a very dilute concentration of fission products reduces vacancy mobility, potentially leading to modifications in point defect supersaturation, void nucleation, and swelling under irradiation. These represent critical challenges for irradiated cladding materials. The tracer diffusion coefficients indicate that Nd diffuses the fastest, followed by La, Ce, and Pr. Furthermore, this study provides essential insights for understanding fission product transport in cladding materials and informs future design strategies to mitigate fuel-cladding chemical interaction, ultimately enhancing nuclear reactor safety and performance.

Diffusion↗
DOE OSTI2024

Control over Conformational Landscapes of Polypeptoids by Monomer Sequence Patterning

The ability to program chain conformation and structure through control over the monomer sequence of synthetic polymers has broad implications for next-generation material design. While related problems of protein-folding and de novo design have generated accurate predictions of 3D folded chain structures, generalization to synthetic polymers remains intractable due to the requirement of large structural databases and the intrinsically disordered nature of polymer building blocks. In this work, polypeptoids, a class of peptidomimetic synthetic polymers, are utilized to build a general workflow for the study of relationships between monomer sequence and dynamic 3D chain structure in solution. Furthermore, this work demonstrates how control over the monomer sequence can alter the conformational landscape of synthetic polymers to deviate dramatically from classical chain statistics. Specifically, the distribution of end-to-end distances, as measured by double electron-electron resonance spectroscopy in dilute solvent, is systematically skewed towards shorter distances with an increasing number of hydrophobes and further refined by hydrophobe arrangement in amphiphilic polypeptoid chains.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Isopentane Disproportionation in Lewis Acidic Chloroaluminate Ionic Liquid

Lewis acidic chloroaluminate ionic liquid diluted in dichloromethane catalyzes the disproportionation of alkanes, such as isopentane, via carbenium ions in a reaction readily initiated by carbenium ion precursors such as tert-butyl chloride (TBC). The carbenium ion-AlCl 4 – ion-pairs stabilized by the ionic liquid-dichloromethane solution are the key intermediates in two distinctive kinetic regimes, i.e., a transient regime (0–5 min) and a steady-state regime (after 5 min). The transient regime constitutes the majority of isopentane conversion and is governed by the initial carbenium ion concentration. In the steady-state regime, disproportionation occurs at a considerably lower rate, affected by the carbenium ion concentration, the concentration of the ionic liquid, and the reaction temperature. The formation of alkenes observed in the 1 H NMR spectra of the reacting substrates, along with the DFT calculations, suggests that deprotonation of carbenium ion-pairs reduces their concentration, decreasing, in turn, the reaction rate. Kinetic modeling indicates that the transient regime is significantly controlled by the hydride transfer (k HT ) and the deprotonation rate constants (k DP ), while the steady-state regime is additionally influenced by the alkene protonation rate constant (k –DP ). The overall activation energy of the reaction at the steady state, expressed as E a,steady-state regime = E a,HT – E a,DP + E a,–DP , was 54 kJ/mol. The reaction mechanism and the kinetics highlight the potential of Lewis acid-catalyzed conversions of hydrocarbons under remarkably mild conditions.

Deprotonation↗
DOE OSTI2026

Modification of ion-temperature-gradient turbulence by impurities in stellarator plasmas

Recent nonlinear gyrokinetic simulations have shown that impurities can strongly modify the turbulent heat flux in stellarator plasmas. Here, the ion-temperature-gradient (ITG) dispersion relation in a plasma containing impurities is analytically solved in certain limits and an expression for the modification of the ITG growth rate by impurities is derived. The analytical expression is the sum of three terms corresponding to three different physical causes (impurity density gradient, impurity temperature gradient and dilution) of the change in the growth rate. The scalings predicted analytically for the modification of the growth rate are shown to be reproduced by linear gyrokinetic simulations. The conditions for reduction or increase of the ITG growth by impurities are also correctly predicted by the analytical solution to the dispersion relation. Finally, a remarkable correlation is found between the analytical expression for the modification of the growth rate and the modification of the turbulent heat flux obtained from nonlinear gyrokinetic simulations.

Calvo, Iván [Research Centre for Energy, Environme↗
DOE OSTI2025

UO2 microstructural evolutions induced by Ni, Mo, and W dopants for intentional forensics

The concept of tagging nuclear fuel with a chemical barcode to enable forensics analysis across the nuclear fuel cycle is an area of active investigation, particularly to ensure fabrication viability without disrupting current fuel performance. This study explored the feasibility of using Ni, Mo, and W isotopic double-spikes as dopants in UO2 fuel from the perspective of fuel fabrication. Doped UO2 pellets were produced using conventional fuel fabrication processes, including powder mixing, sieving, pressing, and sintering in a reductive atmosphere. Two composition levels, 100 and 1000 ppm, were evaluated for each dopant element with isotopic double-spike configurations. For the Ni system, additional dopant concentrations of 250 and 500 ppm were produced with nonperturbed isotopic ratios. The results demonstrated that successful incorporation of Ni, Mo, and W double-spikes into UO2 pellets occurred with minimal shift in final density or dopant loss during pellet fabrication. Isotopic analysis confirmed the presence of the double-spike signature even when diluted with natural isotopic material in ratio of 1:5 in the fabrication process. Microstructural examinations revealed different impacts on grain size compared with undoped UO2. This study showed that Ni incorporation up to ∼500 ppm promoted moderate grain growth, whereas the Mo and W systems caused grain size reduction at all concentrations. Changes in the UO2 lattice parameter as a function of composition were detected exclusively for Ni up to 500 ppm, indicating that the Ni solid solution was the main factor for the observed grain growth. Insoluble (Mo and W) or supersaturated (Ni > 500 ppm) conditions produced grain size reduction. The Ni-doped pellets in the solution range resulted in a final microstructure within fuel specifications, demonstrating its potential benefits of employing complex dopant systems for potential nuclear forensic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Coupling nitrate capture with ammonia production through bifunctional redox-electrodes

Abstract Nitrate is a ubiquitous aqueous pollutant from agricultural and industrial activities. At the same time, conversion of nitrate to ammonia provides an attractive solution for the coupled environmental and energy challenge underlying the nitrogen cycle, by valorizing a pollutant to a carbon-free energy carrier and essential chemical feedstock. Mass transport limitations are a key obstacle to the efficient conversion of nitrate to ammonia from water streams, due to the dilute concentration of nitrate. Here, we develop bifunctional electrodes that couple a nitrate-selective redox-electrosorbent (polyaniline) with an electrocatalyst (cobalt oxide) for nitrate to ammonium conversion. We demonstrate the synergistic reactive separation of nitrate through solely electrochemical control. Electrochemically-reversible nitrate uptake greater than 70 mg/g can be achieved, with electronic structure calculations and spectroscopic measurements providing insight into the underlying role of hydrogen bonding for nitrate selectivity. Using agricultural tile drainage water containing dilute nitrate (0.27 mM), we demonstrate that the bifunctional electrode can achieve a 8-fold up-concentration of nitrate, a 24-fold enhancement of ammonium production rate (108.1 ug h −1 cm −2 ), and a >10-fold enhancement in energy efficiency when compared to direct electrocatalysis in the dilute stream. Our study provides a generalized strategy for a fully electrified reaction-separation pathway for modular nitrate remediation and ammonia production.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2025

Design and Optimization of a Hybrid Design for Quantum Transduction

This study presents the mechanical design and analysis of a quantum electro-optical transducer engineered to operate at millikelvin temperatures within a dilution refrigerator. The transducer enables bidirectional microwave-optical frequency conversion through a hybrid architecture that integrates a superconducting radiofrequency (SRF) cavity with an electro-optic optical cavity. Among several design options investigated, the configuration offering the best thermal and mechanical performance was selected, yielding a robust solution with reduced sensitivity to fabrication tolerances, improved heat dissipation, as well as alignment precision. The design ensures uniform temperature distribution, enabling higher laser pump powers and, thus, increased conversion efficiency, while maintaining mechanical stresses safely below the material yield strength. Electromagnetic simulations further validate the design, demonstrating enhanced coupling between the optical and microwave modes, as well as a broader tuning range achieved with smaller tuner displacements.

quantum computing↗
DOE OSTI2021

High pressure, high flow rate batch mixing apparatus for high throughput experiments

An automated, high pressure, high flow rate batch mixing apparatus has been designed and constructed for rapid, stable, and repeatable mixing of multiple gases and vapors. The apparatus operates as an intermittent batch mixer with cycles of topping off fresh mixture to maintain pressure in an accumulator tank until consumed in an experimental apparatus. At high duty cycles, the apparatus can also function at steady state. This style of mixing is suitable for experiments such as high repetition rate shock tubes and other devices with intermittent flow demands. It is compact and portable, facilitating use in locations such as synchrotron light sources. The entire apparatus is heated to permit the mixing of vapors from species with low volatilities. The apparatus is fully automated and runs for extended periods with the only intervention being to refresh reagent supplies. Here, the accuracy and repeatability of the apparatus were verified by periodic gas sampling and analysis with gas chromatography. Multi-component mixtures spanning a wide range of complexity, dilution, and volatility of constituents have been prepared. The compositions of the majority of the mixture were found to be stable over several filling cycles, repeatable, and with the proper calibration of set-point conditions, accurate. Challenges were encountered preparing a mixture from multi-component liquids, and potential solutions are discussed.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Numerical investigation of microbial quorum sensing under various flow conditions

Microorganisms efficiently coordinate phenotype expressions through a decision-making process known as quorum sensing (QS). We investigated QS amongst distinct, spatially distributed microbial aggregates under various flow conditions using a process-driven numerical model. Model simulations assess the conditions suitable for QS induction and quantify the importance of advective transport of signaling molecules. In addition, advection dilutes signaling molecules so that faster flow conditions require higher microbial densities, faster signal production rates, or higher sensitivities to signaling molecules to induce QS. However, autoinduction of signal production can substantially increase the transport distance of signaling molecules in both upstream and downstream directions. We present empirical approximations to the solutions of the advection–diffusion–reaction equation that describe the concentration profiles of signaling molecules for a wide range of flow and reaction rates. These empirical relationships, which predict the distribution of dissolved solutes along pore channels, allow to quantitatively estimate the effective communication distances amongst multiple microbial aggregates without further numerical simulations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Biphasic solvents for post-combustion CO 2 capture from natural gas flue Gas

Fossil fuel fired power plants are generally expected to remain one of the most significant global sources of electricity for decades to come. Consequently, carbon management technologies are needed to reduce or eliminate ongoing emissions from these sources. Amongst the many techniques for carbon capture, aqueous amine-based absorbents (monoethanolamine, MEA, in particular) are, presently, considered the leading technology for post-combustion CO 2 point-source capture. These technologies are nevertheless limited by their high capital and regeneration energy costs. Biphasic solvents have been identified as an attractive alternative to traditional MEA based absorbents due to their potential energy savings. Thus far, however, the research on biphasic solvents has largely focused on their performance in coal flue gas while more dilute natural gas flue gas applications have received relatively little attention. Here, this work examines the performances of two novel biphasic solvent blends, diethylenetriamine (DETA) and triethylenetetramine (TETA), in CO 2 capture from a natural gas flue gas simulant. Across several regeneration tests, both solvents achieved considerable energy savings over the benchmark MEA solution. Specifically, the energy consumption per mol CO 2 recovered for the DETA-based and TETA-based solvents was 46 % and 35 % less than that of the benchmark MEA solution, respectively. Molecular dynamics simulations were also performed to gain a deeper understanding of the phase separation phenomena that occur as a consequence of CO 2 absorption. These simulations indicated that phase change was driven by the strong interaction between the absorption products and water, while the degree of separation depended on the CO 2 loading.

Biphasic solvents↗
DOE OSTI2020

Analysis of the multicomponent relativistic Boltzmann equation for electron scattering in big bang nucleosynthesis

Big-bang nucleosynthesis (BBN) is valuable as a means to constrain the physics of the early universe and it is the only probe of the radiation-dominated epoch. A fundamental assumption in BBN is that the nuclear velocity distributions obey Maxwell-Boltzmann (MB) statistics as they do in stars. Specifically, the BBN epoch is characterized by a dilute baryon plasma for which the velocity distribution of nuclei is mainly determined by the dominant Coulomb elastic scattering with mildly relativistic electrons. One must therefore deduce the momentum distribution for reacting nuclei from the multicomponent relativistic Boltzmann equation. However, the full multicomponent relativistic Boltzmann equation has only recently been analyzed and its solution has only been worked out in special cases. Moreover, a variety of schemes have been proposed that introduce nonthermal components into the BBN environment which can alter the thermal distribution of reacting nuclei. Here, we construct the relativistic Boltzmann equation in the context of BBN. We also derive a Langevin model and perform relativistic Monte-Carlo simulations which clarify the baryon distribution during BBN and can be used to analyze any relaxation from a nonthermal injection. Further, we show by these analyses that the thermalization process leads to a nuclear distribution function that remains very close to MB statistics even during the most relativistic environment relevant to BBN. Hence, the predictions of standard BBN remain unchanged.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Process Intensification and Scale-Up of a Continuous Enzymatic Hydrolysis and Separation Process

Combining separate unit operations into one where the best of each part can be maximized is one of the benefits of process intensification. An example is the combination of lignocellulosic biomass enzymatic hydrolysis with the downstream solid-liquid separation step to produce clarified sugars ready for fermentation or catalytic upgrading. The productivity and endpoint yield of enzymatic hydrolysis both enjoy the benefits of reduced feedback inhibition through the continuous removal of sugars by incorporating separations into the reactor. Likewise, the efficiency of recovering clarified sugars from the enzymatic hydrolysis slurry can be enhanced by operating separations equipment at steady-state conditions simultaneously with the continuously fed hydrolysis process. A key parameter that enables greater processing capacity while also raising the risks of failure is the solids loading or concentration. Higher solids loading allows for smaller reactor vessels and results in clarified sugars of higher concentration; however, required pumping power increases, reactor agitation may become ineffective, and membrane flux suffers. Feedstock material attributes influenced by upstream pretreatment must also be scrutinized more carefully: dilute-acid pretreated and deacetylated-and-disc-refined feedstocks exhibit different characteristics that affect agitation and pumping. The authors invite you to further explore the process science enabling the scale-up of this technology from conceptual work at the bench to pilot-scale industrially-relevant equipment where the challenges and solutions of integration and process optimization are expounded upon.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Correlating inter-particle forces and particle shape to shear-induced aggregation/fragmentation and rheology for dilute anisotropic particle suspensions: A complementary study via capillary rheometry and in-situ small and ultra-small angle X-ray scattering

Predicting the stability of suspensions composed of anisotropic particles is difficult due to the complexity and interplay of hydrodynamic and colloidal forces. For example, it is not well understood how non-spherical particle shapes affect suspension-rheological responses. This study presents evidence showing how the stacking of boehmite platelets dispersed in aqueous solution leads to moderate shear thinning, even at low boehmite mass fractions. Wide angle X-ray, small angle X-ray, and ultra-small angle X-ray scattering were used to simultaneously monitor changes to the size and fractal dimensions of boehmite aggregates from 6 – 10,000 Å as the sample was recirculated through an in situ capillary rheometer. The latter also provided simultaneous measurements of the viscosity of the suspensions. Computational fluid dynamics modeling of the apparatus was used to provide a more rigorous analysis of the fluid flow. In addition to shear-induced aggregation/fragmentation correlated to a complicated balance between hydrodynamic and colloidal forces, the results indicate that orientationally-dependent interactions, arising from non-spherical boehmite aggregate sub-units, may explain the relatively large experimental viscosity when the hydrodynamic force is small compared to colloidal forces. Ultimately, the results from this study can provide a foundation for future work investigating particle shape effects on colloidal and hydrodynamics forces.

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

Convection-Induced Compositional Patterning at Grain Boundaries in Irradiated Alloys

We consider the stability of precipitates formed at grain boundaries (GBs) by radiation-induced segregation in dilute alloys subjected to irradiation. Here, the effects of grain size and misorientation of symmetric-tilt GBs are quantified using phase field modeling. A novel regime is identified where, at long times, GBs are decorated by precipitate patterns that resist coarsening. Maps of the chemical Péclet number indicate that arrested coarsening takes place when solute advection dominates over thermal diffusion right up to the precipitate/matrix interface, preventing interfacial local equilibrium and overriding capillary effects. This contrasts with liquid-solid mixtures where convection always accelerates coarsening.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Assessing Options for Particulate Measurement Tools to Advance Wood Heaters (CRADA Final Report)

For approximately 40 years, particulate matter (PM) from residential wood heaters has been certified in North America using an EPA dilution tunnel with a filter-based gravimetric measurement. While reasonably repeatable, this approach reports only a single value for the entire burn sequence and does not provide information on PM at different points within the sequence, which would support improved wood heater designs. PM monitors that report with high time resolution offer a potential solution, but their performance for various types of heaters and across heater operating conditions has not been extensively evaluated. In this 18-month collaborative project, the Hearth, Patio & Barbecue Association (HPBA) and Lawrence Berkeley National Laboratory (LBNL) tested five continuous PM instruments spanning three measurement principles—light-scattering, inertial mass (TEOM), and triboelectric—against gravimetric and size-resolved references on three residential heater types (cordwood, catalytic cordwood, and pellet) across startup, high-burn, and low-burn operating phases.

42 ENGINEERING↗
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↗