Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “three phase modeling”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Fierro

FIERRO is a parallel C++ code designed to simulate fluid mechanics, heat transfer, and solid mechanics in two- and three dimensional space. FIERRO is written to run on homogeneous (CPU) and heterogeneous (CPU+GPU) high performance computing machines. Fierro can aid a) modeling and design efforts that have historically relied on commercial implicit and explicit finite element codes, b) numerical methods research, c) manufacturing research, and d) computer science research. The code contains diverse numerical methods to solve the governing physics equations for both quasi-static and dynamic problems. Mathematical optimization solvers are coupled to the numerical methods to research topology and shape optimization that has application to additive manufacturing, and to create novel numerical approaches. Phase-field methods with micromechanical solvers are provided to simulate microstructure formation and evolution in manufacturing processes. The micromechanical solvers can also help research efforts create continuum-scale constitutive models for solids, as a function of the microstructure, in situ in a calculation or in a stand-alone manner. No physical data exists within the code.

Morgan, Nathaniel↗
DOE OSTI2020

Computational Fluid Dynamics Analysis of a 12 MW Circulating Fluidized Bed Rise

A three-dimensional (3D) two fluid model (TFM) model of the 12 MW circulating fluidized bed (CFB) riser at Chalmers University was developed using the open-source Multiphase Flow with Interphase eXchanges (MFiX) Software Suite. The hydrodynamic behavior of the riser was validated by comparing axial pressure profiles from experimental measurements with the results from the simulations. Under some of the operating conditions where there was no external circulation of solids, the pressure profiles from the simulation matched well with the experimental measurements. For the cases where there was external circulation of solids, the pressure profiles with the monodisperse solid phase modeling underpredicted the solids concentration in the freeboard region. With the inclusion of polydispersity in the model, the predicted pressure profiles agreed well with the experimental measurements. A two-dimensional (2D) model of the riser was used to perform simulations of coal combustion. The predicted flue gas outlet compositions of oxygen and carbon dioxide compared reasonably to experiments. The experiments also reported some concentration of carbon monoxide which was overpredicted in the simulations. Additional experimental data such as axial oxygen concentration profile and pressure profiles were also compared with the experiment. It was observed that the oxygen concentration in the bottom bed from the simulations was higher than in the experiments, which suggested that most of the oxidation reactions occurred in the bottom bed. The concentration of volatiles and char in the simulations was lower in the bottom bed than in the experiments, which suggest that there is not sufficient penetration of the solid phase species into the bed.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2025

A phase-field fracture formulation for generalized standard materials: The interplay between thermomechanics and damage

Accurately modeling fracture of ductile materials poses open challenges in the field of computational mechanics due to the multiphysics nature of their failure processes. Integrating the interplay between thermodynamics and damage into ductile fracture models is vital for predicting critical failure modes. Here, in this paper, we develop a versatile phase-field (PF) framework for modeling ductile fracture, taking into account finite-strain elasto-plasticity. The framework stems from a variational formulation of constitutive relations for generalized standard materials (GSMs), whose response is described by a Helmholtz free energy and a dissipation pseudo-potential. Its variational structure is based on a minimum principle for a functional that expresses the sum of power densities for reversible and irreversible processes. By minimizing this functional with a constraint on a von Mises yield function, we derive the evolution equation for the equivalent plastic strain and an associative flow rule. This constrained optimization problem is analytically solved for a wide class of thermo-viscoplasticity models. The key innovations of the current work include (i) a cubic plastic degradation function that accounts for a non-vanishing damage-dependent yield stress, (ii) closed-form expressions of the Helmholtz free energy and dissipation pseudo-potential for three thermo-viscoplasticity models, (iii) an extended Johnson–Cook plasticity model with a nonlinear hardening law, and (iv) a plastic work heat source that depends on the plastic degradation function and a variable Taylor–Quinney (TQ) coefficient. The capabilities of the proposed framework are tested with the aid of four ductile fracture problems, including the Sandia Fracture Challenge. In each of these problems, we examine the evolution of relevant field variables such as the PF order parameter, the equivalent plastic strain, the temperature, and the internal power dissipation density, in addition to the overall structural response quantified by the force–displacement curve. These numerical studies demonstrate that the proposed framework effectively represents ductile fracture, yielding computational results that exhibit good agreement with experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A separations and purification process for improving yields and meeting fuel contaminant specifications for high-octane gasoline produced from dimethyl-ether over a Cu/BEA catalyst

In this work, we have been developing a three-step conversion of biomass-derived syngas to methanol to dimethyl-ether (DME) to non-aromatic hydrocarbons for use as high-octane gasoline and sustainable aviation fuel. This process produces branched alkanes from DME using a Cu/BEA catalyst and is a promising alternative to other syngas conversion processes such as Fischer-Tropsch to linear alkanes and traditional ZSM-5 catalyzed methanol to aromatic gasoline. In this short article we describe some advances in our understanding related to separations and purification via the use of more detailed experimental speciation in an updated process model involving multiple phase equilibrium-based separation steps. Primary modeled reactor outlet constituents (and weight %) are: C3 and lighter hydrocarbon gases (11.1%), C4s (54.5%), H 2 (1.2%), CO 2 (2.9%), water (5.0%), unreacted DME (16.5%), methanol (2.3%), and C5+ hydrocarbons (6.4%). DME (the primary reactant) and H 2 recycle and reuse are important for the overall process efficiency, and the recycle of C4s is important to increase the C5+ yield via reactivation and homologation. Thus H 2 , C4s, and DME are targeted for recycle, while methanol and water need to be removed from the product to conform with fuel specifications. Model predictions from Aspen Plus using the NRTL-RK property method indicate a fuel composition with C5+ content of 97.1 wt%, with minor constituents: 2.4 wt% C4s, 0.3 wt% methanol, 0.1 wt% DME, 0.03 wt% water, and 0.01 wt% C3s. These ranges of minor components conform with fuel quality requirements, and the modeled product is amenable for unconstrained blending to boost gasoline octane ratings.

09 BIOMASS FUELS↗
DOE OSTI2021

The benefits of diligence: how precise are predicted gravitational wave spectra in models with phase transitions?

Models of particle physics that feature phase transitions typically provide predictions for stochastic gravitational wave signals at future detectors and such predictions are used to delineate portions of the model parameter space that can be constrained. The question is: how precise are such predictions? Uncertainties enter in the calculation of the macroscopic thermal parameters and the dynamics of the phase transition itself. We calculate such uncertainties with increasing levels of sophistication in treating the phase transition dynamics. Currently, the highest level of diligence corresponds to careful treatments of the source lifetime; mean bubble separation; going beyond the bag model approximation in solving the hydrodynamics equations and explicitly calculating the fraction of energy in the fluid from these equations rather than using a fit; and including fits for the energy lost to vorticity modes and reheating effects. The lowest level of diligence incorporates none of these effects. We compute the percolation and nucleation temperatures, the mean bubble separation, the fluid velocity, and ultimately the gravitational wave spectrum corresponding to the level of highest diligence for three explicit examples: SMEFT, a dark sector Higgs model, and the real singlet-extended Standard Model (xSM). In each model, we contrast different levels of diligence in the calculation and find that the difference in the final predicted signal can be several orders of magnitude. Our results indicate that calculating the gravitational wave spectrum for particle physics models and deducing precise constraints on the parameter space of such models continues to remain very much a work in progress and warrants care.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The High Luminosity LHC interaction region magnets towards series production

The High Luminosity Large Hadron Collider (HL-LHC) is the new flagship project of CERN. First endorsed in 2013 and approved in 2016, HL-LHC is an upgrade of the accelerator aiming to increase by a factor of ten the statistics of the LHC collisions at the horizon of 2035-2040. HL-LHC relies on cutting edge technologies: among them, large aperture superconducting magnets will replace the present hardware to allow a smaller beam size in two interaction points (IPs). The project involves the construction of about 150 magnets of six different types: the quadrupole triplet, two main dipoles and three orbit correctors. The triplet, manufactured at CERN and in the USA, will consist of 30 magnets based on Nb3Sn technology, with an operational peak field of 11.4 T. These will be the first quadrupole Nb3Sn magnets installed in a particle accelerator. The other five types of magnets, all relying on Nb-Ti technology, present non-trivial challenges in the design and construction; they will be manufactured as part of in-kind contribution under the responsibility of institutes in Japan, China, Spain, and Italy. The project is now in the phase of transition between qualification through short models and prototypes and the beginning of the series construction. In this paper we review the magnet requirements, the reasons for selecting the design, the technological challenges with respect to previous projects, and we summarize the steps that have been taken to validate the baseline.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Impact of moment-based, energy integrated neutrino transport on microphysics and ejecta in binary neutron star mergers

We present an extensive study of the effects of neutrino transport in three-dimensional general relativistic radiation hydrodynamics (GRHD) simulations of binary neutron star (BNS) mergers using our moment-based, energy-integrated neutrino radiation transport (M1) scheme. Here, we consider a total of eight BNS configurations, while varying equation of state models, mass ratios, and grid resolutions, for a total of 16 simulations. We find that M1 neutrino transport is crucial in modeling the local absorption of neutrinos and the deposition of lepton number throughout the medium. We provide an in-depth look at the effects of neutrinos on the fluid dynamics and luminosity during the late inspiral and postmerger phases, the properties of ejecta and outflow, and the postmerger nucleosynthesis. The simulations presented in this work comprise an extensive study of the combined effect of the equation of state and M1 neutrino transport in GRHD simulations of BNS mergers, and establish that the solution provided by our M1 scheme is robust across system properties.

150 ≤ A ≤ 18959 ≤ A ≤ 8990 ≤ A ≤ 149↗
DOE OSTI2025

Magnetic excitations from the hexagonal spin clusters in the 𝑆 = $\frac{1}{2}$ distorted honeycomb lattice antiferromagnet Cu 2 ⁢(pymca)⁢ 3 ⁢(ClO 4 )

Cu 2 ⁢(pymca) ⁢3 (ClO 4 ) (pymca: pyrimidine-2-carboxylate) consists of a slightly distorted honeycomb lattice of Cu 2+ spins, which shows no long-range magnetic order down to 0.6 K. A magnetization study revealed 1/3 and 2/3 plateau phases [A. Okutani et al., J. Phys. Soc. Jpn. 88, 013703 (2019)], which is not expected for regular honeycomb antiferromagnets. Inelastic neutron scattering experiments were performed using a powder sample to investigate the exchange interactions of this material. The spin excitations from the singlet ground state to the first three triplet states, predicted from the antiferromagnetic hexagonal spin cluster interacting with 3.9 meV, were observed. Using the exact diagonalization methods, the intercluster coupling was estimated from the excitation peak width to be about 20% of the intracluster interaction, which is consistent with the previously reported value. Finally, our exchange path model explains the anisotropic exchange interactions in the distorted honeycomb plane.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

CMIP6-based Multi-model Hydroclimate Projection over the Conterminous US, Version 1.1

This dataset presents a suite of high-resolution downscaled hydro-climate projections over the conterminous United States (CONUS) based on multiple selected Global Climate Models (GCMs) from the Coupled Models Intercomparison Project phase 6 (CMIP6). The CMIP6 GCMs are downscaled using either statistical (DBCCA) or dynamical (RegCM) downscaling approaches based on two meteorological reference datasets (Daymet and Livneh). Subsequently, the downscaled precipitation, temperature, and wind speed are employed to drive two calibrated hydrologic models (VIC and PRMS), enabling the simulation of projected future hydrologic responses across the CONUS. Each ensemble member covers the 1980-2019 baseline and 2020-2059 near-future periods under the high-end (SSP585) emission scenario. Moreover, utilizing only DBCCA and Daymet, the projections are further extended to the 2060-2099 far-future period and across three additional emission scenarios (SSP370, SSP245, and SSP126). This dataset is formulated to support the SECURE Water Act Section 9505 Assessment for the US Department of Energy (DOE) Water Power Technologies Office (WPTO). For further details on this dataset, please refer to Kao et al. (2022) and Rastogi et al. (2022).

13 HYDRO ENERGY↗
DOE OSTI2022

Removal of Non-Specifically Bound Proteins Using Rayleigh Waves Generated on ST-Quartz Substrates

Label-free biosensors are plagued by the issue of non-specific protein binding which negatively affects sensing parameters such as sensitivity, selectivity, and limit-of-detection. In the current work, we explore the possibility of using the Rayleigh waves in ST-Quartz devices to efficiently remove non-specifically bound proteins via acoustic streaming. A coupled-field finite element (FE) fluid structure interaction (FSI) model of a surface acoustic wave (SAW) device based on ST-Quartz substrate in contact with a liquid loading was first used to predict trends in forces related to SAW-induced acoustic streaming. Based on model predictions, it is found that the computed SAW body force is sufficient to overcome adhesive forces between particles and a surface while lift and drag forces prevent reattachment for a range of SAW frequencies. We further performed experiments to validate the model predictions and observe that the excitation of Rayleigh SAWs removed non-specifically bound (NSB) antigens and antibodies from sensing and non-sensing regions, while rinsing and blocking agents were ineffective. An amplified RF signal applied to the device input disrupted the specific interactions between antigens and their capture antibody as well. ST-quartz allows propagation of Rayleigh and leaky SH-SAW waves in orthogonal directions. Thus, the results reported here could allow integration of three important biosensor functions on a single chip, i.e., removal of non-specific binding, mixing, and sensing in the liquid phase.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Simulating exotic phases of matter with bond-directed interactions with arrays of Majorana-Cooper pair boxes

It is suggested that networks of Majorana-Cooper pair boxes connected by metallic nanowirescan simulate various exotic states of matter. In this simulations Majorana-Cooper boxes play therole of effective spins S=1/2 and the metallic connections generate the Kondo screening and theRuderman-Kittel-Kasuya-Yosida (RKKY) interaction. Depending on what prevails - whether it isthe Kondo effect or the RKKY exchange, one will have either an effective spin model or a Kondolattice. The list of exotic stets includes the famous hexagonal Kitaev model, a generalization of thismodel for a Kondo lattice and various spin models with three-spin interactions. A special emphasizeis made on the discussion of the Kondo lattice scenario.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

How to identify different new neutrino oscillation physics scenarios at DUNE

Next generation neutrino oscillation experiments are expected to measure the remaining oscillation parameters with very good precision. They will have unprecedented capabilities to search for new physics that modify oscillations. DUNE, with its broad band beam, good particle identification, and relatively high energies will provide an excellent environment to search for new physics. If deviations from the standard three-flavor oscillation picture are seen however, it is crucial to know which new physics scenario is found so that it can be verified elsewhere and theoretically understood. We investigate several benchmark new physics scenarios by looking at existing long-baseline accelerator neutrino data from NOvA and T2K and determine at what sensitivity DUNE can differentiate among them. We consider sterile neutrinos and both vector and scalar non-standard neutrino interactions, all with new complex phases, the latter of which could conceivably provide absolute neutrino mass scale information. We find that, in many interesting cases, DUNE will have good model discrimination. We also perform a new fit to NOvA and T2K data with scalar NSI.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Cascade of phase transitions in a planar Dirac material

Here we investigate a model of interacting Dirac fermions in 2 + 1 dimensions with M flavors and N colors having the U(M)×SU(N ) symmetry. In the large-N limit, we find that the U(M) symmetry is spontaneously broken in a variety of ways. In the vacuum, when the parity-breaking flavor-singlet mass is varied, the ground state undergoes a sequence of M first-order phase transitions, experiencing M + 1 phases characterized by symmetry breaking U(M)→U(M – k)×U(k) with k ϵ {0, 1, 2, · · · , M}, bearing a close resemblance to the vacuum structure of three-dimensional QCD. At finite temperature and chemical potential, a rich phase diagram with first and second-order phase transitions and tricritical points is observed. Also exotic phases with spontaneous symmetry breaking of the form as U(3)→U(1) 3 , U(4)→U(2)×U(1) 2 , and U(5)→U(2)2×U(1) exist. For a large flavor-singlet mass, the increase of the chemical potential μ brings about M consecutive first-order transitions that separate the low-μ phase diagram with vanishing fermion density from the high-μ region with a high fermion density.

1/N expansion↗
DOE OSTI2022

Verification and Validation of a Modified Numerical Algorithm for Simulation of Transient Unconfined Groundwater Flow

This report extends verification and validation testing of the NUFT package of codes (Nitao 2000a,b) to the US1P module of code designed for simulation of single-phase mixed or “coupled” saturated and unsaturated (or “variably saturated”) groundwater flow. Importantly, the verification and validation testing of US1P in this report includes performance evaluation of a modified numerical algorithm that was not included in prior quality assurance (QA) of the NUFT package by Carle et al. (2014) for Underground Test Area (UGTA) project activities, which are directed at assessment of radionuclide contamination in groundwater sourced from underground nuclear test locations at the Nevada National Security Site (NNSS). An immediate purpose of this report is to provide QA for ongoing large-scale, three-dimensional (3-D) groundwater flow modeling of transient water levels associated with longterm water supply pumping and underground nuclear testing at the NNSS (Jackson and Fenelon, 2018; Jackson et al., 2021). For UGTA activities, the QA refers to the standards of ASTM (1996) for verification and validation testing of groundwater modeling codes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Determining the Solubility Behavior of Kogarkoite in Simulated Nuclear Waste

Kogarkoite (Na 3 FSO 4 ) is a sparingly soluble fluoride–sulfate double salt that has been identified in high level nuclear waste sludge at the Hanford Site and, more recently, in sludge batch compilation samples at the Savannah River Site (SRS). Due to its complex dissolution behavior, which exhibits an inverse dependence on sodium ion activity, the presence of this mineral poses significant challenges to waste retrieval and processing. Incomplete dissolution during sludge washing can lead to the retention of fluoride and sulfate in the high-level waste feed, potentially causing the formation of corrosive, immiscible molten salt layers, known as "glass gall,” in vitrification melters. Current efforts to optimize flowsheet parameters and wash-water volumes are hindered by the absence of a commercially available, certified reference material, which prevents the accurate calibration of analytical methods and the verification of dissolution kinetics. To address this critical gap, this research focuses on the laboratory synthesis of pure Kogarkoite to serve as a standard for comprehensive solubility and washing performance testing. A coupled synthesis and simulant campaign was executed using an evaporative crystallization protocol designed to replicate the dynamic concentration effects observed in tank farm operations. Thirteen simulant matrices were prepared by dissolving systematically varied ratios of sodium fluoride (NaF) and sodium sulfate (Na 2 SO 4 ) in deionized water under three distinct caustic regimes: 0.0 g (control), 4.0 g (~1 M), and 12.0 g (~3 M) sodium hydroxide (NaOH). While thermodynamic equilibrium models suggest that high-caustic environments should favor the stability of the double salt7, results from this evaporative study at 25 0 C revealed a distinct kinetic divergence. Simulants with high hydroxide loading predominantly yielded large, blocky crystals of sodium sulfate decahydrate (Na 2 SO 4 .10H 2 O). Successful synthesis of pure Kogarkoite was achieved exclusively in specific NaOH-free compositional windows, where the precipitate manifested as fine, opaque granular aggregates. Ion chromatography (IC) analysis confirmed phase purity through the simultaneous stoichiometric depletion of both fluoride and sulfate from the supernatant. This successful synthesis establishes a reproducible route to generate bulk Kogarkoite, enabling the subsequent phase of quantitative dissolution testing using inhibited water to optimize sludge-batch assembly.

Sarker, Md Sharif [Florida International Univ. (FI↗
DOE OSTI2021

Parametric study of solid-solid translucent phase change materials in building windows

Thermal energy storage and solar radiation management are crucial to improve the sustainability and energy efficiency of buildings. Compared with the implementation of phase change materials (PCMs) in opaque components, the energy saving potential of incorporating PCMs in transparent glazing windows is much less studied and not well understood. Here we present a comprehensive parametric study of novel PCM windows for building energy saving with a focus on optimizing and quantitatively distinguishing the contributions from the optical and thermal properties of the PCM, which is particularly useful for the design of solid-solid PCM windows. We investigate a reference commercial office building using EnergyPlus by developing an equivalent model of our PCM window that is compatible with EnergyPlus's modeling capabilities. Compared with a clear-clear double-pane window, the integration of 3 mm solid-solid PCMs with optimal properties in warm, mixed, and cold climates can respectively save up to 17.2%, 14.0%, and 5.8% energy for the HVAC (heating, ventilation, and air conditioning) system, and 9.4%, 6.7%, and 3.2% energy for the whole building. We also demonstrate that these energy savings are most sensitive to the solar absorptance of PCMs for all three climates. The optimal transition temperature varies with climate and is related to the climate and solar radiation heat gain. Other issues are also briefly discussed, such as hysteresis, window orientations, and the effect of interior lighting. Finally, although the optimal PCM windows show energy saving performance comparable with low-emissivity windows, the PCM windows provide a unique advantage in terms of shifting HVAC loads which can provide benefits to the electrical grid.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

A field–particle correlation analysis of a perpendicular magnetized collisionless shock

Using the field–particle correlation technique, we examine the particle energization in a three-dimensional (one spatial dimension and two velocity dimensions; 1D-2V) continuum Vlasov–Maxwell simulation of a perpendicular magnetized collisionless shock. The combination of the field–particle correlation technique with the high-fidelity representation of the particle distribution function provided by a direct discretization of the Vlasov equation allows us to ascertain the details of the exchange of energy between the electromagnetic fields and the particles in phase space. We identify the velocity-space signatures of shock-drift acceleration of the ions and adiabatic heating of the electrons arising from the perpendicular collisionless shock by constructing a simplified model with the minimum ingredients necessary to produce the observed energization signatures in the self-consistent Vlasov–Maxwell simulation. We are thus able to completely characterize the energy transfer in the perpendicular collisionless shock considered here and provide predictions for the application of the field–particle correlation technique to spacecraft measurements of collisionless shocks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Numerical Study of Coupled Water and Vapor Flow, Heat Transfer, and Solute Transport in Variably-Saturated Deformable Soil During Freeze-Thaw Cycles

As climate change intensifies, soil water flow, heat transfer, and solute transport in the active, unfrozen zones within permafrost and seasonally frozen ground exhibit progressively more complex interactions that are difficult to elucidate with measurements alone. For example, frozen conditions impede water flow and solute transport in soil, while heat and mass transfer are significantly affected by high thermal inertia generated from water-ice phase change during the freeze-thaw cycle. To assist in understanding these subsurface processes, the current study presents a coupled two-dimensional model, which examines heat conduction-convection with water-ice phase change, soil water (liquid water and vapor) and groundwater flow, advective-dispersive solute transport with sorption, and soil deformation (frost heave and thaw settlement) in variably saturated soils subjected to freeze-thaw actions. This coupled multiphysics problem is numerically solved using the finite element method. The model's performance is first verified by comparison to a well-documented freezing test on unsaturated soil in a laboratory environment obtained from the literature. Then based on the proposed model, we quantify the impacts of freeze-thaw cycles on the distribution of temperature, water content, displacement history, and solute concentration in three distinct soil types, including sand, silt and clay textures. The influence of fluctuations in the air temperature, groundwater level, hydraulic conductivity, and solute transport parameters was also comparatively studied. The results show that (a) there is a significant bidirectional exchange between groundwater in the saturated zone and soil water in the vadose zone during freeze-thaw periods, and its magnitude increases with the combined influence of higher hydraulic conductivity and higher capillarity; (b) the rapid dewatering ahead of the freezing front causes local volume shrinkage within the non-frozen region when the freezing front propagates downward during the freezing stage and this volume shrinkage reduces the impact of frost heave due to ice formation. This gradually recovers when the thawed water replenishes the water loss zone during the thawing stage; and (c) the profiles of soil moisture, temperature, displacement, and solute concentration during freeze-thaw cycles are sensitive to the changes in amplitude and freeze-thaw period of the sinusoidal varying air temperature near the ground surface, hydraulic conductivity of soil texture, and the initial groundwater levels. Our modeling framework and simulation results highlight the need to account for coupled thermal-hydraulic-mechanical-chemical behaviors to better understand soil water and groundwater dynamics during freeze-thaw cycles and further help explain the observed changes in water cycles and landscape evolution in cold regions.

54 ENVIRONMENTAL SCIENCES↗