Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Texture evolution”

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 181 records · Page 10

On the exceptionally high ductility of Mg–2Zn-0.3Ca-0.2Ce-0.1Mn alloy

A newly developed Mg–2Zn-0.3Ca-0.2Ce-0.1Mn (in weight %, ZXEM2000) wrought alloy exhibits an exceptionally high tensile ductility (29% ± 1.4% tensile elongation) along the rolling direction at room temperature. In this work, we investigated the fundamental plastic deformation mechanisms that can contribute to the high ductility of ZXEM2000. Scanning electron microscopy and transmission electron microscopy were carried out at selected strain levels to characterize the microstructural evolution. Compared with the deformed microstructures from the hot-rolled pure Mg and Mg–3Al–1Zn alloy at similar conditions, our results revealed that the high ductility of ZXEM2000 could be attributed to four key factors: weaker texture, fine grain size, reduced twinning, and increased cross-slip frequency of < c +a > dislocation. Furthermore, these direct observations offer a key insight into improving ductility of magnesium alloys for future alloy design and development.

36 MATERIALS SCIENCE↗

Part 1: Aspects of lithospheric evolution on Venus. Part 2: Thermal and collisional histories of chondrite parent bodies

The geological evolution of distinctly different kinds of solar system objects is addressed. Venus has been observed over the past decade by orbital radars on both American and Soviet spacecraft. These surface measurements provide clues to the structure and evolution of the lithosphere. The parent bodies of chondritic meteorites, thought to resemble asteroids, represent the other end of the size spectrum of terrestrial objects. Their early thermal and collisional histories may be constrained by the chemical and textural record preserved in meteorite samples. Impact craters on Venus have been observed by the Soviet Venera 15/16 spacecraft. A formalism is presented by which the size-frequency distribution of impact craters may be used to estimate upper bounds on the mean global rates of volcanic resurfacing and lithospheric recycling on that planet over the past several hundred million years. The impact crater density reported from Venera observations, if valid for the entire Venus surface, indicates a mean volcanic flux no greater than 2 cu km/y, corresponding to a maximum average rate of resurfacing of about 4 km/b.y. For the lowest estimated mean crater retention age of the surface of Venus imaged by Venera 15/16, the rate of lithospheric recycling on Venus does not exceed 1.5 sq km/y. Ordinary chondrite meteorites show textural and chemical patterns indicative of varying intensities of thermal metamorphism. The conventional onion-shell model, which envisions highly metamorphosed material in the core and less intensely heated rocks near the surface, predicts an inverse relation between peak temperature and cooking rate, but none has been observed. A metamorphosed-planetesimal model is devised to explain this discrepancy, whereby heating occurs in planetesimals a few kilometers in radius which then accrete to form 100-km-radius parent bodies. Cooling rates are then randomly controlled by burial depth. Thermal and collisional constraints are examined, and the model is found to be applicable only to highly insulating Al-26-rich planetesimals that remain closely aggregated upon accretion. An alternative model is presented, in which onion-shell parent bodies are collisionally fragmented during metamorphism and then gravitationally reassembled. If reassembly times are short, then cooling rates would be determined by burial depth in the reaccreted parent body. This model, unlike previous ones, can explain both coherent and incoherent cooling of Breccia clasts by collisions during or after metamorphism, respectively.

Grimm, Robert E.↗

Eocene-Oligocene Metamorphism, Fluid Flow and Deformation in the Ruby Mountains-East Humboldt Range Metamorphic Core Complex

Metamorphic core complexes are ubiquitous in collapsed orogens globally and play a primary role in crustal exhumation. Here, we investigate the metamorphic history of the Ruby Mountains-East Humboldt Range metamorphic core complex, Nevada, using petrochronology to understand how magmatism, metamorphism and deformation interact to modulate crustal rheology, and the timing of exhumation within the Sevier orogenic belt. Field study, microstructural analysis, thermobarometry and laser ablation split-stream monazite, titanite and allanite petrochronology of the mylonitic footwall were integrated to elucidate the pressure-temperature-time (P-T-t) evolution. Major-, trace-element and quartz-in-garnet thermobarometry show peak metamorphism occurred at 5.5–6 kbar and 600–650°C across the study area. Monazite and titanite U-(Th)-Pb petrochronology constrain this episode of metamorphism to 88–81 Ma with the structurally deepest samples yielding dates down to ca. 71 Ma. The preservation of peak metamorphic mineral assemblages indicates these rocks remained at near-peak temperature conditions into the Cenozoic. Cenozoic metamorphism occurred during a punctuated episode contemporaneous with the Great Basin ignimbrite flare up. All samples show a distinct 39–37 Ma population of monazite, titanite and allanite dates that correlate with the emplacement of quartz diorite and gabbro intrusions throughout the footwall. Accessory phases contain zoning textures in backscattered electron images that are consistent with fluid-mediated dissolution-reprecipitation reactions. The U-(Th)-Pb analyses from these zones yield 39–27 Ma dates interpreted as the timing of fluid flow. The results of this study elucidate the integrated importance of magmatism, melting and fluid flow in driving the coupled evolution of Late Cretaceous metamorphism within the Sevier hinterland and the subsequent Cenozoic metamorphic core complex development that exhumed the middle crust.

58 GEOSCIENCES↗

Lunar breccias, petrology, and earth planetary structure

Topics covered include: (1) petrologic studies of poikiloblastic textured rocks; (2) petrology of aluminous mare basalts in breccia 14063; (3) petrology of Apollo 15 breccia 15459; (4) high-alumina mare basalts; (5) some petrological aspects of imbrium stratigraphy; (6) petrology of lunar rocks and implication to lunar evolution; (7) the crystallization trends of spinels in Tertiary basalts from Rhum and Muck and their petrogenetic significance; (8) the geology and evolution of the Cayman Trench; (9) The petrochemistry of igneous rocks from the Cayman Trench and the Captains Bay Pluton, Unalaska Island and their relation to tectonic processes at plate margins; and (10) the oxide and silicate mineral chemistry of a Kimberlite from the Premier Mine with implications for the evolution of kimberlitic magma.

Ridley, W. I.↗

Accretion and Heating of Particles by Supersonic Planetesimals

Recently, nebular shock waves have become one of the leading candidates for explaining the presence of chondrules in primitive meteorites. While shocks have been shown to be capable of explaining many of the features of chondrules, a major problem with the theory is that the source of the shocks remains unidentified. Among the suggested sources of the shocks are bow shocks created by supersonic planetesimals in the nebula. Recently, we studied the structure of the shocks that would form around such supersonic planetesimals. We found that particles that encountered the shocks at distances greater than 2 planetesimal radii from the planetesimal would cool too quickly to form the textures observed in chondrules. While the region of the shock far away from the planetesimal may not allow chondrules to form, the region closer in has not been studied in detail. In this work we consider the dynamical and thermal evolution of particles that encounter supersonic planetesimals in this region.

Ciesla, F. J.↗

Dynamics and Formation of Antiferromagnetic Textures in MnBi 2 Te 4 Single Crystal

We report coherent X-ray imaging of antiferromagnetic (AFM) domains and domain walls in MnBi 2 Te 4 , an intrinsic AFM topological insulator. This technique enables direct visualization of domain morphology without reconstruction algorithms, allowing us to resolve antiphase domain walls as distinct dark lines arising from the A-type AFM structure. The wall width is determined to be 550(30) nm, in good agreement with earlier magnetic force microscopy results. The temperature dependence of the AFM order parameter extracted from our images closely follows previous neutron scattering data. Remarkably, however, we find a pronounced hysteresis in the evolution of domains and domain walls: upon cooling, dynamic reorganizations occur within a narrow ∼1 K interval below 𝑇 𝑁 , whereas upon warming, the domain configuration remains largely unchanged until AFM order disappears. These findings reveal a complex energy landscape in MnBi 2 Te 4 , governed by the interplay of exchange, anisotropy, and domain-wall energies, and underscore the critical role of AFM domain-wall dynamics in shaping its physical properties. These sharply defined and hysteretically evolving walls may provide a controllable AFM texture in MnBi 2 Te 4 , hinting at potential use in low-power spintronic devices based on domain-wall dynamics.

36 MATERIALS SCIENCE↗

Conditions of formation of pyroxene excentroradial chondrules

Understanding the exact mechanism of origin of chondrules and the attendant physical conditions remains a key factor in understanding the early evolution of the solar system. The present investigation is concerned with an interpretation of the range of cooling histories needed to make chondrules. Chondrule-like spherules were formed in dynamic crystallization experiments, using gas-mixing facilities. It is found that textures in pyroxene-rich chondrules are well reproduced in the laboratory by cooling melts of chondrule composition from just above the liquidus. The range of pyroxene dendrite widths is slightly greater for natural chondrules than for droplets cooled from 50 C/hr to greater than 3000 C/hr. A blanketing medium of variable thickness explains the cooling rates lower than expected for radiative cooling. Experimental results, textures indicating incomplete melting, and isotopic disequilibrium are compatible with reheating of primitive material in the early solar nebula.

Hewins, R. H.↗

Benchmarking Microscale Ductility Measurements (Final Report of the Project DE-NE0008799)

Conventional macroscale experimentation is generally considered to be straightforward with few limitations. Conversely, micro/nanoscale experimentation presents numerous challenges in loading device design, sample preparation and handling, as well as accurate understanding of grain size and local texture effects on recorded measurements. Despite these challenges, nanopillar compression, MEMs based micro-tension, and nanoindentation approaches have been able to provide fundamental contributions to the understanding of material behavior at small lengthscales. However, the overarching shortcoming of these micro/nanoscale experimentation approaches, is the inability to directly translate measurements evaluated at the nm and µm length scales (e.g., hardness) to macroscale tensile material behavior (i.e., elastic modulus, yield strength, and ductility). The objectives of the proposed study are, 1) to establish best practices for obtaining tensile microscale ductility measurements, and 2) to validate methodologies to for comparing microscale ductility measurements to macroscale ductility measurements. In order to achieve these objectives, a multi-lengthscale, multi-temperature testing protocol and simulation framework are executed first on copper as a model material to validate the following approach, and second on reactor grade Zircaloy-2. Experiments are conducted on specimens extracted from the same test piece to ensure nominally identical grain size and texture from specimens to specimen. Motivated by the need to isolate the contribution of size-effects on obtained mechanical property measurements, specimens are manufactured with thicknesses at the micro- (1-10 µm), meso- (10-100's µm), and macroscales (sub-sized ASTM E8). In-situ full-field deformation techniques (scanning electron microscopy (SEM) grid methods and optical DIC) are incorporated into testing at each specimen length-scale to capture plasticity localization and evolution. Experimental testing for all specimens is conducted at both room temperature and elevated temperatures to probe the role of thermal activation on plastic deformation accommodation processes. Simulation efforts focus on examining the mechanical behavior of microscale specimens using a finite element approach with explicitly resolved grain morphologies, and an embedded crystal plasticity model. The cost-efficient implementation method allows for the modeling of a statistically significant number of both real (i.e., digital twin) and generated microstructures to obtain an understanding of the interrelationships between specimen microstructure and geometric variables (grain size, texture, specimen geometry, etc.) on microscale mechanical behavior.

36 MATERIALS SCIENCE↗

Benthic Marine Cyanobacterial Mat Ecosystems: Biogeochemistry and Biomarkers

Cyanobacterial mats are complete ecosystems that can include processes of primary production, diagenesis and lithification. Light sustains oxygenic photosynthesis, which in turn provides energy, organic matter and oxygen to the community. Due to both absorption and scattering phenomena, incident light is transformed with depth in the mat, both in intensity and spectral composition. Mobile photo synthesizers optimize their position with respect to this light gradient. When photosynthesis ceases at night, the upper layers of the mat become reduced and sulfidic. Counteracting gradients of oxygen and sulfide combine to provide daily-contrasting environments separated on a scale of a few mm. The functional complexity of mats, coupled with the highly proximal and ordered spatial arrangement of biota, offers the potential for a staggering number of interactions. At a minimum, the products of each functional group of microorganisms affect the other groups both positively and negatively. For example, cyanobacteria generate organic matter (potential substrates) but also oxygen (a toxin for many anaerobes). Anaerobic activity recycles nutrients to the photosynthesizers but also generates potentially toxic sulfide. The combination of benefits and hazards of light, oxygen and sulfide promotes the allocation of the various essential mat processes between light and dark periods, and to various depths in the mat. Observations of mats have produced numerous surprises. For example, obligately anaerobic processes can occur in the presence of abundant oxygen, highly reduced gases are produced in the presence of abundant sulfate, meiofauna thrive at high sulfide concentrations, and the mats' constituent populations respond to environmental changes in complex ways. While photosynthetic bacteria dominate the biomass and productivity of the mat, nonphotosynthetic, anaerobic processes constitute the ultimate biological filter on the ecosystem's emergent biosignatures, including those sedimentary textures, organic compounds, and minerals that enter the fossil record. The ability of cyanobacterial mats to channel abundant solar energy into the creation and maintenance of complex structures and processes has created a multitude of consequences, both for sedimentation and for the early evolution of our biosphere.

DesMarais, David J.↗

Understanding process-microstructure-property relationships in laser powder bed fusion of non-spherical Ti-6Al-4V powder

Powder feedstock is a major cost driver in metal additive manufacturing (AM). Replacing the spherical powder with the cost-efficient non-spherical one can reduce the feedstock cost up to 50% and attract more interest to adopt AM in production and new alloy development. Here, for this paper, a comprehensive study was conducted to understand process-microstructure-property relationships in laser powder bed fusion of hydride-dehydride Ti-6Al-4V powder. We demonstrated that variation of laser scan speed had a significant impact on the grain structure, pore evolution and properties compared to laser power. Dynamic X-ray radiography showed that with decreasing scan speed at a constant laser power, a transition from conduction to keyhole mode laser processing occurred, in which a deeper melt pool at lower scan speed intensified texture. In other words, an increase in laser scan speed resulted in formation of the refined prior β grains with shape factor of ~5, lowering the anisotropy. Furthermore, the degree of variant selection was evaluated based on the analyzed texture as a function of laser power and scan speed. With increasing laser scan speed, the dominant α/α boundary type was altered from type 2 to 4 and the degree of variant selection was noticeably decreased. On the other hand, increasing laser power left the morphology of prior β grains, their size, and the dominant α/α boundary (type 4) unchanged, while the texture and anisotropy were intensified, and the degree of variant selection was slightly decreased. Finally, dependency of surface roughness and microhardness were discussed as a function of laser processing parameters.

variant selection↗

Magnetic domains and domain wall pinning in atomically thin CrBr 3 revealed by nanoscale imaging

The emergence of atomically thin van der Waals magnets provides a new platform for the studies of two-dimensional magnetism and its applications. However, the widely used measurement methods in recent studies cannot provide quantitative information of the magnetization nor achieve nanoscale spatial resolution. These capabilities are essential to explore the rich properties of magnetic domains and spin textures. Here, we employ cryogenic scanning magnetometry using a single-electron spin of a nitrogen-vacancy center in a diamond probe to unambiguously prove the existence of magnetic domains and study their dynamics in atomically thin CrBr 3 . By controlling the magnetic domain evolution as a function of magnetic field, we find that the pinning effect is a dominant coercivity mechanism and determine the magnetization of a CrBr 3 bilayer to be about 26 Bohr magnetons per square nanometer. The high spatial resolution of this technique enables imaging of magnetic domains and allows to locate the sites of defects that pin the domain walls and nucleate the reverse domains. Our work highlights scanning nitrogen-vacancy center magnetometry as a quantitative probe to explore nanoscale features in two-dimensional magnets.

36 MATERIALS SCIENCE↗

Classification of JERS-1 Image Mosaic of Central Africa Using A Supervised Multiscale Classifier of Texture Features

In this paper, a multiscale approach is introduced to classify the Japanese Research Satellite-1 (JERS-1) mosaic image over the Central African rainforest. A series of texture maps are generated from the 100 m mosaic image at various scales. Using a quadtree model and relating classes at each scale by a Markovian relationship, the multiscale images are classified from course to finer scale. The results are verified at various scales and the evolution of classification is monitored by calculating the error at each stage.

Saatchi, Sassan↗

The granulite suite: Impact melts and metamorphic breccias of the early lunar crust

The granulite suite consists of two major types of rocks. One is coarse-grained and poikilitic with many euhedral crystals of olivine and plagioclase. These characteristics indicate crystallization from a melt; the poikilitic granulites are impact melt breccias. The other group is finer-grained and granoblastic, with numerous triple junctions; the granoblastic granulites are metamorphic rocks. Compositional groups identified by Lindstrom and Lindstrom contain both textural types. Two pyroxene thermometry indicates that both groups equilibrated at 1000 to 1150 C. Calculations suggest that the granoblastic group, which has an average grain size of about 80 microns, was annealed for less than 6 x 10 exp 4 y at 1000 C, and for less than 2500 y at 1150 C. Similar equilibration temperatures suggest that both groups were physically associated after impact events produced the poikilitic melts. Granulitic impactites hold important information about the pre-Nectarian bombardment history of the Moon, and the composition and thermal evolution of the early lunar crust. Granulitic impactites are widely considered to be an important rock type in the lunar crust, but how they formed is poorly understood. Metal compositions and elevated concentrations of meteoritic siderophile elements suggest that most lunar granulites are impact breccias. Their occurrence as clasts in approximately 3.9 Ga breccias, and Ar-(40-39) ages greater than or = 4.2 Ga for some granulites show that they represent a component of the lunar crust which formed prior to the Nectarian cataclysm. Petrographic characteristics of lunar granulites indicate at least two endmember textural variants which apparently formed in fundamentally different ways. One type has granoblastic textures consisting of equant, polygonal to rounded grains, and abundant triple junctions with small dispersions around 120 degrees indicating a close approach to textural equilibrium. As suggested by many authors, granoblastic granulites probably formed by subsolidus annealing and recrystallization of fragmental or glassy protoliths. Examples of this type include 15418, 78155, and 79215. The other textural type consists of poikilitic to poikiloblastic rocks with euhedral to subhedral plagioclase and olivine enclosed by interstitial pyroxene. In some cases, the texture resembles that of an orthocumulate. Examples of this type include 60035, 67955, and 77017. Rounding of grain edges is common in poikilitic granulites, but the regular crystal shapes and widely dispersed dihedral angles show they are far from textural equilibrium. The textures of poikilitic granulites are more consistent with the formation of these rocks by crystallization from a melt than by subsolidus metamorphism. A few samples have been recognized with textural characteristics transitional between those of the granoblastic and poikiloblastic endmembers (e.g., 72559, 78527). Pyroxene compositions taken from the literature and determined for this study by electron microprobe were used to calculate equilibration temperatures. The Kretz Ca transfer (solvus) thermometer and the Lindsley and Anderson graphical method both give similar temperatures, which range from approximately 1000 to 1150 C. There is no apparent temperature difference between granoblastic and poikilitic varieties, but there is a hint in these data that the more ferroan varieties equilibrated to lower temperatures. Additional studies are in progress to test this possibility.

Cushing, J. A.↗

Myco-Architecture off Planet: Growing Surface Structures at Destination

Our work focused on filling major select key technical knowledge gaps, and technical aspects to be addressed in a Phase II proposal were identified such as the secretion of glues or plastics by the fungi or bacteria to form biocomposites. Specifically, we raised the technical TRL by assessing the growth of in-house mycelial-producing strains on potential food substrates and analyzed the advantages and disadvantages to their use off planet. Growth of the mycelia on sawdust and powdered nutrients including dried, powdered algae, was used as a baseline. We measured growth as a function of temperature, tested the relevant material properties of the mycelia products produced, and measured the mass of input and output volume for production per volume of material. Mycotecture was produced in a bag simulating the proposed mission implementation. Post-production such as heat treatment (as is done with terrestrial mycotecture)was assessed.The material properties of the dried, frozen and baked mycelial outputs included density, strength,thermal insulation and flame retardation. Tensile/compression testing on biocomposite material specimens was performed to obtain three-dimensional elastic constants, as well as to examine deformation and fracture behavior. These mechanical tests were performed on both dormant and activated samples to understand the evolution of the structural material. In addition, experimental data obtained from the mechanical testing was used to build a failure prediction model that accounts for material anisotropy. This aided in future structural design through a quantitative understanding of the mechanical limits of the material. To examine the texture, defects, fracture surfaces of the material, high-resolution microstructural imaging was utilized before and after the mechanical tests. This microstructural analysis informed us about the macro-structural behavior and influenced structural design. Thermal decomposition analysis was also performed tounderstand the thermal limits.

Biocomposites↗

Incorporating dynamic recrystallization into a crystal plasticity model for high-temperature deformation of Ti-6Al-4V

During hot deformation of (α + β) titanium alloys, the simultaneous action of strain and temperature in the (α + β) regime facilitates dynamic recovery, dynamic recrystallization (DRX), and phase transformations via non-equilibrium paths. DRX is manifested in the form of fine recrystallized α or β grains. Here, in the present study, a two-phase crystal plasticity finite element framework (CP-DRX) has been developed which incorporates DRX kinetics into the crystal plasticity (CP) model to predict the flow characteristics of Ti-6Al-4V alloys during thermo-mechanical processing. The CP slip system parameters, as well as elastic properties from both α and β phases of Ti, are calibrated for different strain rate conditions. An EBSD-informed two-phase microstructure representation has been utilized in the CP-DRX framework to explore the time-dependent evolution of DRX microstructure and crystal orientation for different strain rate conditions. The proposed CP-DRX can capture the evolution of crystal orientation and plastic flow stress–strain response of polycrystalline Ti-6Al-4V during the deformation process. Furthermore, the proposed model is able to capture the softening behavior, observed in average stress–strain response from experiments performed using a Gleeble thermomechanical simulator and predict the recrystallization texture.

(α+ β)-Titanium (Ti) Alloys↗

Evidence for a single-layer van der Waals multiferroic

Here we searched for signatures of type-II multiferroic order in ultrathin van der Waals materials, to determine if such a state, comprising coexisting electrical polarization and spin order with strong magnetoelectric effects, could persist down to a single-atomic layer. We reported the first observation of such a state in triangular lattice material NiI 2 . This system has a layered structure and built-in magnetic frustration that promotes a complex magnetic ground state with proper-screw spin helices. The spin helix pattern breaks inversion symmetry and induces electrical polarization via magnetoelectric coupling. We use circular dichroic Raman measurements to directly probe the magneto-chiral ground state and its electromagnon modes originating from dynamic magnetoelectric coupling. Using birefringence and second-harmonic generation measurements, we detect a highly anisotropic electronic state simultaneously breaking three-fold rotational and inversion symmetry. The evolution of the optical signatures as a function of temperature and layer number surprisingly revealed a multiferroic state that persists down to the ultrathin limit of monolayer NiI 2 . These observations establish NiI 2 and transition metal dihalides as a new platform for studying emergent multiferroic phenomena, chiral magnetic textures and ferroelectricity in the 2D limit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Temperature-dependent mechanical properties and crystal plasticity parameters for additively manufactured Haynes-214 alloy: Experiments and numerical modeling

Our experimental mechanical testing data demonstrated that the additively manufactured (AM) laser powder bed fusion (L-PBF) Haynes-214 alloy exhibits non-linear mechanical properties as the temperature rises from ambient to 870 °C. Crystal plasticity (CP) simulations provide an effective approach to gaining deeper insights into microstructure-property linkages under thermomechanical loading. This method can reduce the need for costly high-temperature mechanical testing while accounting for the effects of crystallographic texture and grain morphology on the mechanical behavior of AM materials. However, calibrating a CP model is time-consuming because individual simulations are computationally expensive and hundreds (or more) of iterations over parameter sets may be required. To address this issue, we have designed a machine learning-differential evolution (ML-DE) CP framework that can accurately interpolate the tensile properties of AM L-PBF Haynes-214 alloy across a wide temperature range from ambient to 870 °C, with minimal reliance on experimental data. The framework uses electron backscatter diffraction (EBSD) measurements to generate statistically equivalent microstructural volume elements to serve as inputs to the CP modeling framework. Stress–strain curves were generated from 1000 CP simulations, which serve as the training data set for the three ML regression algorithms explored: linear, extra-trees, and multi-layer perceptron. These three regression models were independently evaluated to compare their efficiency and identify the most suitable algorithm for the given problem. Results revealed that the extra-trees ML regressor outperforms the other models in both qualitative and quantitative aspects with an R 2 of 0.98. Subsequently, the differential evolution optimization approach is employed to calibrate the ML-based CP material parameters with experimental results obtained at various temperatures. Finally, temperature-dependent CP material parameters are formulated. The effectiveness and efficiency of the designed framework are validated through comparison with experimental results, demonstrating a high degree of agreement. These calibrated parametric constitutive equations enable further use of the CP model to study the deformation behavior of this alloy under a wide range of thermo-mechanical loading conditions.

36 MATERIALS SCIENCE↗

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↗