SEARCH · Engineering Papers

Results for “diffusion couples”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Diffusive crack growth at a bimaterial interface

The diffusional microcrack growth behavior in a bimaterial system is investigated with an aim at estimating service life of advanced ceramic composites under creep-rupture conditions. The crack is assumed to grow via a coupled surface and grain-boundary diffusion under steady state conditions. The tensile stress distribution along the interface ahead of the moving crack tip is solved, and it is found that a new length parameter exists as a scaling factor for which the solution becomes identical to the single phase case when plotted on the nondimensional physical plane. In contrast to the elastic stress solution which shows singularity at the tip, together with oscillatory character away from the tip, the creep stresses have a peak value away from the tip due to a wedging effect and interfacial sliding eliminates stress oscillation resulting in a decoupling between mode I and mode II loading. This solution ties the far-field loading parameter to the crack tip conditions in terms of the unknown crack velocity to give a specific V-K(1) relationship. It is shown that an exponent of 12 in the conventional crack growth power law emerges at the higher applied stress range.

Chuang, Tze-Jer↗

DOE OSTI2022

Elucidating the impact of the ionomer equivalent weight on a platinum group metal-free PEMFC cathode via oxygen limiting current

Leveraging the interactions between ionomer and catalyst can increase the performance of proton exchange membrane fuel cells. The impacts of the equivalent weight (EW) of perfluorosulfonic acid–based ionomers on the platinum group metal-free electrode structure and fuel cell performance have not been fully explored. Four membrane electrode assemblies (MEAs) were prepared by using a commercial Fe–N–C catalyst, two perfluorosulfonic acid ionomers with different EWs, that is, Aquivion 720 (A720) and Nafion 1100 (N1100), and two ionomer-to-catalyst (I/C) ratios. The four MEAs were characterized to understand the impact of the ionomer EW and content on the capacitance, proton conductivity, and mass transport on the cathode. The mass transport resistance was measured for the first time using a new oxygen reduction reaction limiting current method enabling to couple the effects of oxygen diffusion with liquid water generation. Low EW ionomer combined with a moderate I/C results in improved performance due to its enhanced proton conductivity. However, when used at high I/C, it can cause severe water flooding at high current density due to the enhanced liquid water uptake, especially at high relative humidity, resulting in lower catalyst utilization and higher mass transport resistance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Using a coupled dispersion model to estimate depletion of a tritium oxide plume by a forest

Tritium processing facilities may release tritium oxide (HTO) to the atmosphere which poses potential health risks to exposed co-located workers and to offsite individuals. Most radiological consequence analyses determine HTO dose by applying Gaussian plume models to simulate the transport of HTO. Within these models, deposition velocity is used to assess the sum of all deposition processes acting on the plume. While this may account for vegetative and soil uptake or respiration processes, it may currently lack inclusion of the complex interactions within heterogeneous forested environments. In this complex morphology, dispersion patterns are significantly altered by changing flow regimes above and below the forest canopy and by the transfer of plume material across the canopy boundary. To determine the effects of a heterogeneous forest canopy on an airborne HTO plume, a Gaussian plume model coupled with an advection-diffusion plume model was applied to estimate transport in the free atmosphere above the forest and within the forest canopy and understory. During 2012, wind speed and wind direction measurements taken at 5 heights, ranging from 2-m to 28-m, on an instrumented meteorological tower located in a loblolly pine forest at the Department of Energy (DOE) Savannah River Site (SRS), near Aiken, SC. From these measurements, model predictions were made over a full spectrum of meteorological conditions. Deposition and resuspension velocities were calculated based on the model-predicted flux of plume material across the top of the forest canopy. Additionally, net deposition velocity of the plume material was calculated as the difference between the deposition and resuspension velocities. The 1st and 5th percentile net deposition velocities were estimated to be 0.7 cm s-1 and 1.2 cm s-1, respectively.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

DOE OSTI2022

MEUMAPPS (C++ Version)

Many materials, metal alloys in particular, have features on the on micrometer or nanometer scale that have a large impact on the properties of the material. These features are known as the microstructure of the material. Understanding why and how the microstructure forms in a material is of fundamental scientific interest as well as of significant technological interest. The capability to predict microstructure evolution in a material allows the intentional design of microstructures and hence the intentional design of material properties. The phase-field method is one of the leading methods for predicting microstructure evolution. One of the most significant problems for phase-field models is their computational expense. Even limited phase-field simulations can easily require thousands of CPU core-hours to complete, which significantly limits their use. This code provides both a general framework for creating scalable, GPU-accelerated phase-field model applications as well as several applications themselves. The code is capable of using hundreds of GPUs efficiently, which greatly reduces the time required to perform simulations. The code is written with an emphasis on performance portability, that is the ability for the code to run efficiently on a number of different computing architectures without modification of the source code. The performance portability of this code is primarily enabled through the use of two libraries, Kokkos (performance portable data structures and execution patterns) and heFFTe (performance portable distributed 3D fast Fourier transforms). The code consists of a core library, applications, and tests. The core library includes shared functionality between applications. This includes interfaces with fast Fourier transform (FFT) libraries such as heFFTe, data structures based on Kokkos, file input and output capabilities, and a solver for infinitesimal strain mechanical equilibrium problems. Five applications are included in the code. The flagship application is the MEUMAPPS-SS application, which implements the Kim-Kim-Suzuki phase-field model for precipitation for an arbitrary number of phases and components in a metal alloy. Five simpler applications are also included that solve the Eshelby inclusion problem, Allen-Cahn equation, the coupled Allen-Cahn and diffusion equations, and the Cahn-Hilliard equation. The code includes two applications to solve the Cahn-Hilliard equation, one with constant-step-size first-order time integration and the second with adaptive high-order time integration.

DeWitt, Stephen [Oak Ridge National Lab. (ORNL), O↗

NASA NTRS1979

Chemically frozen multicomponent boundary layer theory of salt and/or ash deposition rates from combustion gases

There is increased interest in, and concern about, deposition and corrosion phenomena in combustion systems containing inorganic condensible vapors and particles (salts, ash). To meet the need for a computationally tractable deposition rate theory general enough to embrace multielement/component situations of current and future gas turbine and magnetogasdynamic interest, a multicomponent chemically 'frozen' boundary layer (CFBL) deposition theory is presented and its applicability to the special case of Na2SO4 deposition from seeded laboratory burner combustion products is demonstrated. The coupled effects of Fick (concentration) diffusion and Soret (thermal) diffusion are included, along with explicit corrections for effects of variable properties and free stream turbulence. The present formulation is sufficiently general to include the transport of particles provided they are small enough to be formally treated as heavy molecules. Quantitative criteria developed to delineate the domain of validity of CFBL-rate theory suggest considerable practical promise for the present framework, which is characterized by relatively modest demands for new input information and computer time.

Rosner, D. E.↗

NASA NTRS1991

An investigation of voids formation mechanisms and their effects on freeze and thaw processes of lithium and lithium fluoride

The mechanisms of void formation during the cooldown and freezing of lithium coolant within the primary loop of SP-100 type systems are investigated. These mechanisms are: (1) homogeneous nucleation; (2) heterogeneous nucleation; (3) normal segregation of helium gas dissolved in liquid lithium; and (4) shrinkage of lithium during freezing. To evaluate the void formation potential due to segregation, a numerical scheme that couples the freezing and mass diffusion processes in both the solid and liquid regions is developed. The results indicated that the formation of He bubbles is unlikely by either homogeneous or heterogeneous nucleation during the cooldown process. However, homogeneous nucleation of He bubbles following the segregation of dissolved He in liquid lithium ahead of the solid-liquid interface is likely to occur. Results also show that total volume of He void is insignificant when compared to that of shrinkage voids. In viewing this, the subsequent research focuses on the effects of shrinkage void forming during freezing of lithium on subsequent thaw processes are investigated using a numerical scheme that is based on a single (solid/liquid) cell approach. The cases of lithium-fluoride are also investigated to show the effect of larger volume shrinkage upon freezing on the freeze and thaw processes. Results show that a void forming at the wall appreciably reduces the solid-liquid interface velocity, during both freeze and thaw, and causes a substantial rise in the wall temperature during thaw. However, in the case of Li, the maximum wall temperature was much lower than the melting temperature of PWC-11, which is used as the structure material in the SP-100 system. Hence, it is included that a formation of hot spots is unlikely during the startup or restart of the SP-100 system.

El-Genk, Mohamed S.↗

NASA NTRS1991

Mechanisms of voids formation during cooldown and freezing of lithium in SP-100 type systems

The mechanisms of void formation during the cooldown and freezing of lithium coolant within the primary loop of SP-100 type systems are investigated. These mechanisms are: (a) homogeneous nucleation, (b) heterogeneous nucleation, (c) normal segregation of helium gas dissolved in liquid lithium, and (d) shrinkage of lithium during freezing. To evaluate the void formation potential due to segregation, a numerical scheme that couples the freezing and mass diffusion processes in both the solid and liquid regions is developed. The results indicated that the formation of He bubbles is unlikely by either homogeneous or heterogeneous nucleation during the cooldown process. However, homogeneous nucleation of He bubbles following the segregation of dissolved He in liquid Lithium ahead of the solid-liquid interface is likely to occur. Results also show that total volume of He void is insignificant when compared to that of shrinkage voids.

Yang, Jae Y.↗

NASA NTRS1992

Solar-wind/magnetospheric dynamos: MHD-scale collective entry of the solar wind energy, momentum and mass into the magnetosphere

A quasi open MHD (Magnetohydrodynamic) scale anomalous transport controlled boundary layer model is proposed, where the MHD collective behavior of magnetofluids (direct dynamo effect, anomalous viscous interaction and anomalous diffusion of the mass and the magnetic field) plays the main role in the conversion of the Solar Wind (SW) kinetic and magnetic energy into electromagnetic energy in the Magnetosphere (MSp). The so called direct and indirect dynamo effects are based on inductive and purely dissipative energy conversion, respectively. The self organization ability of vector fields in turbulent magnetofluids implies an inductive response of the plasma, which leads to the direct dynamo effect. The direct dynamo effect describes the direct formation of localized field aligned currents and the transverse Alfven waves and provides a source for MHD scale anomalous diffusivity and viscosity. The SW/MSp coupling depends on the dynamo efficiency.

Song, Yan↗

NASA NTRS2003

Unified Models of Turbulence and Nonlinear Wave Evolution in the Extended Solar Corona and Solar Wind

The PI (Cranmer) and Co-I (A. van Ballegooijen) made significant progress toward the goal of building a "unified model" of the dominant physical processes responsible for the acceleration of the solar wind. The approach outlined in the original proposal comprised two complementary pieces: (1) to further investigate individual physical processes under realistic coronal and solar wind conditions, and (2) to extract the dominant physical effects from simulations and apply them to a one-dimensional and time-independent model of plasma heating and acceleration. The accomplishments in the report period are thus divided into these two categories: 1a. Focused Study of Kinetic MHD Turbulence. We have developed a model of magnetohydrodynamic (MHD) turbulence in the extended solar corona that contains the effects of collisionless dissipation and anisotropic particle heating. A turbulent cascade is one possible way of generating small-scale fluctuations (easy to dissipate/heat) from a pre-existing population of low-frequency Alfven waves (difficult to dissipate/heat). We modeled the cascade as a combination of advection and diffusion in wavenumber space. The dominant spectral transfer occurs in the direction perpendicular to the background magnetic field. As expected from earlier models, this leads to a highly anisotropic fluctuation spectrum with a rapidly decaying tail in the parallel wavenumber direction. The wave power that decays to high enough frequencies to become ion cyclotron resonant depends on the relative strengths of advection and diffusion in the cascade. For the most realistic values of these parameters, though, there is insufficient power to heat protons and heavy ions. The dominant oblique waves undergo Landau damping, which implies strong parallel electron heating. We thus investigated the nonlinear evolution of the electron velocity distributions (VDFs) into parallel beams and discrete phase-space holes (similar to those seen in the terrestrial magnetosphere) which are an alternate means of heating protons via stochastic interactions similar to particle-particle collisions. 1b. Focused Study of the Multi-Mode Detailed Balance Formalism. The PI began to explore the feasibility of using the "weak turbulence," or detailed-balance theory of Tsytovich, Melrose, and others to encompass the relevant physics of the solar wind. This study did not go far, however, because if the "strong" MHD turbulence discussed above is a dominant player in the wind's acceleration region, this formalism is inherently not applicable to the corona. We will continue to study the various published approaches to the weak turbulence formalism, especially with an eye on ways to parameterize nonlinear wave reflection rates. 2. Building the Unified Model Code Architecture. We have begun developing the computational model of a time-steady open flux tube in the extended corona. The model will be "unified" in the sense that it will include (simultaneously for the first time) as many of the various proposed physical processes as possible, all on equal footing. To retain this generality, we have formulated the problem in two interconnected parts: a completely kinetic model for the particles, using the Monte Carlo approach, and a finite-difference approach for the self-consistent fluctuation spectra. The two codes are run sequentially and iteratively until complete consistency is achieved. The current version of the Monte Carlo code incorporates gravity, the zero-current electric field, magnetic mirroring, and collisions. The fluctuation code incorporates WKJ3 wave action conservation and the cascade/dissipation processes discussed above. The codes are being run for various test problems with known solutions. Planned additions to the codes include prescriptions for nonlinear wave steepening, kinetic velocity-space diffusion, and multi-mode coupling (including reflection and refraction).

Cranmer, Steven R.↗

DOE OSTI2024

Elucidating the Link Between Alkali Metal Ions and Reaction-Transport Mechanisms in Cathode Electrodes for Alkali-ion Batteries

Our long-term goal is to improve the reliability of electrode materials and their ability to transport and store various metal ions for electrochemical energy storage applications. The main objective of this work was to investigate the intrinsic relationship between the role of alkali metal ions and electrochemically driven mechanical stability and kinetic properties of battery materials. The overall question was “What is the role of alkali metal ions on the electrochemical and mechanical behavior of cathode electrodes? Our guiding hypothesis was that intercalation of larger alkali metal ions (Na and K) inevitably alters the coupled transport-reaction processes during battery operation in organic electrolytes, leading to more intensive chemo-mechanical instabilities in cathode electrodes, resulting in rapid capacity fade. To validate the hypothesis, we experimentally characterized the reaction-transport processes and governing forces driving the instability of electrode materials in different alkali metal-ion environments. The project had three main tasks. The first one was to investigate intercalation-induced strains and associated stress generation, and their impact on structural deformations in composite cathode electrodes. The second task focused on identifying potential-dependent dynamic changes in the electrode-electrolyte interface in alkali metal ion batteries. The last task was focused on determining how larger alkali metal ions with slower diffusivity affect the transport-mechanics coupling at faster scan rates, compared to smaller ions with faster diffusivity in electrodes. We shortly provided the outcome of each task in the accomplishment section. This project produced 10 peer-reviewed publications (9 research papers and one review manuscript) and supported two Ph.D. students, who graduated from Oklahoma State University.

25 ENERGY STORAGE↗

DOE OSTI2025

Mechanistic understanding of carbon mineralization in fracture systems using microfluidics

Carbon mineralization in mafic and ultramafic rocks presents an opportunity for permanent carbon storage in the Earth's subsurface. However, due to their lower permeability, pre-existing fracture networks are key for mineralization to occur. Therefore, to fully develop this technology, a mechanistic understanding of the mineralization behavior in fractures with the consideration of hydrodynamic components is required. We use high-pressure microfluidics to investigate key mechanisms influencing dissolution–precipitation in a fracture network. The experiments were conducted in micromodels made of natural rocks with a comb-shaped flow channel to mimic a fracture network. This enabled studying the effect of injection rate on coupled dissolution–precipitation in advection and diffusion-dominated flow paths. We used gypsum carbonation as an analog reaction to allow for realistic experimental time frames due to its rapid reaction kinetics. The experimental work is coupled with high-fidelity numerical simulations to enhance our understanding of the parameters affecting the mineralization reaction. Our results demonstrate the importance of flow rate on the rate and nature of the gypsum carbonation reaction revealing that higher flow rates enable deeper penetration of the mineral precipitation front into the dead-end channels. This is an important finding since for sustained mineralization in a fracture network, precipitation in dead-ends while still allowing for flowing fractures is critical. Detailed characterization of the precipitates showed that lower flow rates led to porous and loose precipitates in the form of aragonite while higher flow rates mimicked supersaturation behavior leading to the formation of calcite. The reactive transport simulations further demonstrated the significance of flow velocity in advection-dominated channels to influence the efficiency of carbon mineralization in diffusion-dominated channels, potentially clogging of dead-end channels. These findings highlight the need for coupling chemical, mechanical, and hydrodynamic processes to evaluate the nature and extent of carbon mineralization in fractured media critical for permanent storage in mafic and ultramafic formations. This research further highlights the need for more investigation in potential subsurface fracture generation techniques to aid carbon mineralization.

25 ENERGY STORAGE↗

NASA NTRS1987

Directional solidification in an imperfect furnace

The upward directional solidification of a binary alloy in which the rejected solute is more dense than the solvent is considered. The system is stably stratified both thermally and solutally. If heat is lost through the sidewalls, horizontal thermal gradients drive a steady double-diffusive convection. Solidification cells with small aspect ratios (width/length) A are considered, and asymptotic solutions to the steady double-diffusive system are obtained as A tends to 0 under the assumption heat losses are also small. When the thermal Grashof number is large, buoyancy boundary layers are at the sidewalls. The shape of the crystal-melt interface is determined by the coupled effects of local thermodynamic equilibrium, diffusion and convection.

Brattkus, K.↗

NASA NTRS1993

Precombination Cloud Collapse and Baryonic Dark Matter

A simple spherical model of dense baryon clouds in the hot big bang 'strongly nonlinear primordial isocurvature baryon fluctuations' is reviewed and used to describe the dependence of cloud behavior on the model parameters, baryon mass, and initial over-density. Gravitational collapse of clouds before and during recombination is considered including radiation diffusion and trapping, remnant type and mass, and effects on linear large-scale fluctuation modes. Sufficiently dense clouds collapse early into black holes with a minimum mass of approx. 1 solar mass, which behave dynamically like collisionless cold dark matter. Clouds below a critical over-density, however, delay collapse until recombination, remaining until then dynamically coupled to the radiation like ordinary diffuse baryons, and possibly producing remnants of other kinds and lower mass. The mean density in either type of baryonic remnant is unconstrained by observed element abundances. However, mixed or unmixed spatial variations in abundance may survive in the diffuse baryon and produce observable departures from standard predictions.

Hogan, Craig J.↗

DOE OSTI2021

Spatial profiles of photon chemical potential in near-field thermophotovoltaic cells

Emitted photons stemming from the radiative recombination of electron–hole pairs carry chemical potential in radiative energy converters. This luminescent effect can substantially alter the local net photogeneration in near-field thermophotovoltaic cells. Several assumptions involving the luminescent effect are commonly made in modeling photovoltaic devices; in particular, the photon chemical potential is assumed to be zero or a constant prescribed by the bias voltage. The significance of photon chemical potential depends upon the emitter temperature, the semiconductor properties, and the injection level. Hence, these assumptions are questionable in near-field thermophotovoltaic devices with nanoscale separation distances between the emitter and the cell. In the present work, an iterative solver that combines fluctuational electrodynamics with the drift-diffusion model is developed to tackle the coupled photon and charge transport problem, enabling the determination of the spatial profile of photon chemical potential beyond the detailed balance approach. Here, the difference between the results obtained by allowing the photon chemical potential to vary spatially and by assuming a constant value demonstrates the limitations of the conventional approaches. This study is critically important for performance evaluation of near-field thermophotovoltaic systems with nanoscale vacuum gaps.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2023

On the Impact of Mechanics on Electrochemistry of Lithium-Ion Battery Anodes

Abstract Models exploring electrochemistry-mechanics coupling in liquid electrolyte lithium-ion battery anodes have traditionally incorporated stress impact on thermodynamics, bulk diffusive transport, and fracture, while stress-kinetics coupling is more explored in the context of all solid-state batteries. Here, we showcase the existence of strong link between active particle surface pressure and reaction kinetics affecting performance even in liquid electrolyte systems. Traction-free and immobile particle surface mechanical boundary conditions are used to delineate the varying pressure magnitudes in graphite host during cycling. Both tensile and compressive stresses are generated in traction-free case, while a fixed surface subjects the entire particle to a compression state. Pressure magnitudes are nearly two to three orders of magnitude higher for the latter resulting in significant depression of open circuit potential and improvement of exchange current densities compared to stress-free state. The results demonstrate the need for incorporating stress-kinetics linkage in models and provide a rationale for putting battery electrodes under compression to improve kinetics.

25 ENERGY STORAGE↗

DOE OSTI2021

Atomistic observation on diffusion-mediated friction between single-asperity contacts

Super-lubricity, a sliding regime in which friction almost vanishes, has been well-documented with the advent of atomic force microscopy (AFM) and molecular dynamics (MD) simulation due to its enormous implications for manufacturing micro-/nano-electromechanical systems (MEMS/NEMS). Nevertheless, the atomistic mechanisms for super-lubricious behaviors are still elusive, primarily due to a lack of direct observation of interface at atomic resolution during frictional process. Here by using in situ high-resolution transmission electron microscopy (HRTEM) coupled with AFM, we report diffusion-mediated formation of a loosely-packed interfacial layer between two metallic asperities at infinitesimal normal forces, which is responsible for the sliding regime, super-lubricity exhibiting ultra-low friction forces and continuous sliding. By contrast, the loosely-packed interfacial layer vanishes along with the super-lubricity, leading to typical stick-slip friction. Additionally, in situ TEM observation and MD simulations reveals the crucial role that atom diffusion plays in atomic friction, and provides new insights into the fundamental mechanisms of super-lubricity.

42 ENGINEERING↗

DOE OSTI2021

Multiphysics and Multiscale Modeling of Coupled Transport of Chloride Ions in Concrete

Chloride ions (Cl−)-induced corrosion is one of the main degradation mechanisms in reinforced concrete (RC) structures. In most situations, the degradation initiates with the transport of Cl− from the surface of the concrete towards the reinforcing steel. The accumulation of Cl− at the steel-concrete interface could initiate reinforcement corrosion once a threshold Cl− concentration is achieved. An accurate numerical model of the Cl− transport in concrete is required to predict the corrosion initiation in RC structures. However, existing numerical models lack a representation of the heterogenous concrete microstructure resulting from the varying environmental conditions and the indirect effect of time dependent temperature and relative humidity (RH) on the water adsorption and Cl− binding isotherms. In this study, a numerical model is developed to study the coupled transport of Cl− with heat, RH and oxygen (O2) into the concrete. The modeling of the concrete microstructure is performed using the Virtual Cement and Concrete Testing Laboratory (VCCTL) code developed by the U.S. National Institute of Standards and Technology (NIST). The concept of equivalent maturation time is utilized to eliminate the limitation of simulating concrete microstructure using VCCTL in specific environmental conditions such as adiabatic. Thus, a time-dependent concrete microstructure, which depends on the hydration reactions coupled with the temperature and RH of the environment, is achieved to study the Cl− transport. Additionally, Cl− binding isotherms, which are a function of the pH of the concrete pore solution, Cl− concentration, and weight fraction of mono-sulfate aluminate (AFm) and calcium-silicate-hydrate (C-S-H), obtained from an experimental study by the same authors are utilized to account for the Cl− binding of cement hydration products. The temperature dependent RH diffusion was considered to account for the transport of Cl− with moisture transport. The temperature and RH diffusion in the concrete domain, composite theory, and Cl− binding and water adsorption isotherms are used in combination, to estimate the ensuing Cl− diffusion field within the concrete. The coupled transport process of heat, RH, Cl−, and O2 is implemented in the Multiphysics Object-Oriented Simulation Environment (MOOSE) developed by the U.S. Idaho National Laboratory (INL). The model was verified and validated using data from multiple experimental studies with different concrete mixture proportions, curing durations, and environmental conditions. Additionally, a sensitivity analysis was performed to identify that the water-to-cement (w/c) ratio, the exposure duration, the boundary conditions: temperature, RH, surface Cl− concentration, Cl− diffusion coefficient in the capillary water, and the critical RH are the important parameters that govern the Cl− transport in RC structures. In a case study, the capabilities of the developed numerical model are demonstrated by studying the complex 2D diffusion of Cl− in a RC beam located in two different climatic regions: warm and humid weather in Galveston, Texas, and cold and dry weather in North Minnesota, Minnesota, subjected to time varying temperature, RH, and surface Cl− concentrations.

composite theory↗

DOE OSTI2021

Hygro-thermo-mechanical modeling of partially saturated air-entrained concrete containing dissolved salt and exposed to freeze-thaw cycles

Highlights: • A model for cyclic freezing of air-entrained concrete containing salt is proposed. • The freeze-thaw hysteresis, also including the effect of dissolved salt, is modeled. • The diffusive and convective flux of salt is coupled to the freeze-thaw processes. • The freeze-thaw response of air-entrained concrete containing salt is reproduced. • Liquid absorption from reservoirs increases with an increased concentration of salt. In cold regions, understanding the freeze-thaw behavior of air-entrained concrete is important for designing durable structures and assessing the remaining service life of existing structures. This study presents a hygro-thermo-mechanical multiphase model that describes the cyclic freeze-thaw behavior of partially saturated air-entrained concrete containing dissolved salt. An equilibrium and a non-equilibrium approach are adopted to model the ice formation, including the freeze-thaw hysteresis, inside the porous network. The model also considers the diffusive and convective transport of the dissolved salt coupled to the freeze-thaw processes. Two examples are presented to verify and highlight the capabilities of the model. The first example shows that the model is capable of reproducing the experimentally observed mechanical response of specimens containing NaC1-solutions of different concentrations. In the second example, a larger absorption of liquid from an external reservoir is obtained with an increasing salt concentration in the reservoir, which is consistent with experimental observations.

36 MATERIALS SCIENCE↗