Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mixture 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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Modeling the reflectance spectrum of Callisto 0.25 to 4.1 microns

Models employing Hapke's (1981) radiative transfer theory are presented for the 0.2-4.1 micron reflectance spectrum. A simultaneous intimate, as well as aereal, mixture solution of ice and dark material is envisioned, in order to satisfy both absorption band depths and reflectance levels. The models indicate that the surface's ice component is rather large-grained, and that the major spectral features above about 2.5 microns cannot be accounted for by the ice. Spectra obtained for the nonice material were similar to each other; their absorption features resemble hydrated silicates bearing both oxidation states of iron.

Calvin, Wendy M.↗
DOE OSTI2021

An automated multi-component gas adsorption system (MC GAS)

The knowledge gap on adsorption of complex mixtures in the literature relative to single component data represents a persistent obstacle to developing accurate process models for adsorption separations. The collection of mixed gas adsorption data is an imminent need for improved understanding of the behavior of adsorbent systems in these diverse adsorption applications. Current approaches to understanding mixture adsorption using predictive theories based on pure component adsorption experiments often fail to capture the behavior of more complex, non-ideal systems. In this work, we present an automated volumetric instrument for the measurement of mixed gas adsorption isotherms. This instrument was validated by comparison to other in-house instruments and data available in the literature, and the binary adsorption measurements were found to be thermodynamically consistent. The automation of this instrument allows for rapid collection of high-quality mixture adsorption data.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Catalytic co-pyrolysis of LDPE and PET with HZSM-5, H-beta, and HY: experiments and kinetic modelling

In this study, the catalytic pyrolysis of low-density polyethylene (LDPE), polyethylene terephthalate (PET) and their mixture (1 : 1 wt/wt) with three zeolite catalysts (HZSM-5, H-beta, HY) was investigated using a thermogravimetric analyzer (TGA) and a Pyroprobe® micro-reactor coupled to a gas chromatograph mass spectrometer (Py-GC/MS). The TGA results demonstrated that during pyrolysis at 10 °C min –1 , on average, zeolite catalysts decreased the maximum decomposition temperature by 149 °C for LDPE while only decreasing by 8 °C for PET. The derivative thermogravimetric (DTG) curve evidenced interactions when the two polymers were catalytically co-pyrolyzed for all the three catalysts. A lumped nth order reaction scheme was able to accurately model both non-catalytic and catalytic pyrolysis and co-pyrolysis by using the least squares fitting approach for determining the kinetic parameters. The kinetic model was able to model well the interaction effects observed during catalytic co-pyrolysis of LDPE and PET with HZSM-5, H-beta, and HY (Fit% Wt% > 96%, Fit% DTG > 94%). Py-GC/MS experiments for the catalytic fast pyrolysis of LDPE revealed that HZSM-5 resulted in the highest selectivity to aromatic hydrocarbons (31.6%) and HY resulted in the highest selectivity to gasoline range C 5 –C 10 paraffins and olefins (40.9%). Catalytic fast pyrolysis of PET showed high selectivity to benzene for all the catalysts (>43%) and that HZSM-5 resulted in the highest selectivity to polyaromatic hydrocarbons (24.7%). Finally, the catalytic fast co-pyrolysis of LDPE and PET revealed interaction effects for all the three catalysts evidenced by a positive synergy% for alkylated benzenes (3–142%) and polyaromatics (105–187%) with a concomitant negative synergy% for benzene (24–36%) and C 5 –C 10 paraffins and olefins (27–53%).

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

Extended-MHD simulations of disruption mitigation via massive gas injection in SPARC

Recent developments to the M3D-C1 code enable higher fidelity modeling of disruptions, and can be applied in the design verification of reactor-scale tokamaks. Among these new capabilities is a method to mesh conducting vessel structures such as coils and passive plates, packing of the toroidal mesh around gas injectors, as well as anisotropic resistivity inside the vessel structures. We present extended-magnetohydrodynamic (MHD) simulations of disruption mitigation via massive gas injection (MGI) in SPARC. The goal of this study is to inform the disruption mitigation layout of SPARC and aid in the design of an effective gas injector configuration. Fully three-dimensional simulations with M3D-C1 are carried out for various injector configurations with the primary goal of determining the effect of different MGI parameters on heat loads and vessel forces. The simulations include a model for impurity ionization, recombination, advection and radiation, as well as spatially resolved conducting structures around the plasma. A localized mixture of deuterium and neon with a small toroidal and poloidal width is injected in up to six locations. We demonstrate that M3D-C1 can model a rapid shutdown via MGI using narrow and more realistic gas plumes than in previous simulations. As a result of the q = 1 surface in the SPARC baseline case a sawtooth is observed early in the simulations. Despite the sawtooth and the onset of edge MHD instabilities, the impurity distribution remains localized around the injector locations, but enables a radiative shutdown of the plasma. We find that using the maximum of six gas injectors results in a lower peaking factor and leads to a more even distribution of radiation toroidally than using two injectors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1975

Energy transfer by radiation in non-grey atomic gases in isothermal and non-isothermal slabs

A multiband model for the absorption coefficient of atomic hydrogen-helium plasmas is constructed which includes continuum and line contributions. Emission from 28 stronger lines of 106 that have been screened is considered, of which 21 are from hydrogen and 7 belong to helium, with reabsorption due to line-line, line-continuum overlap accurately accounted for. The model is utilized in the computation of intensities and fluxes from shock-heated slabs of 85% H2-15% He mixtures for slab thicknesses from 1 to 30 cm, temperature from 10,000 to 20,000 K, and for different densities. In conjunction with the multiband model, simple numerical schemes have been devised which provide a quick and comprehensive way of computing radiative energy transfer in nonisothermal and nongrey gases.

Poon, P. T. Y.↗
DOE OSTI2021

Mathematical modeling of novel porous transport layer architectures for proton exchange membrane electrolysis cells

Thin foil based porous transport layers (PTLs) that contain highly structured pore arrays have shown promise as anode PTLs in proton exchange membrane electrolysis cells. These novel PTLs, fabricated with advanced manufacturing techniques, produce thin, tunable, multifunctional layers with reduced flow and interfacial resistances and high thermal and electric conductivities. To further optimize their design, it is important to understand their fundamental impact on the transport of protons, electrons, and liquid/vapor mixtures in the electrode. In this work, we develop a two-dimensional multiphysics model to simulate the coupled electrochemistry and multiphase transport in an electrolysis cell operated with the novel PTL architecture. The results show that larger pores improve access of water to the anode catalyst layer, which is beneficial for both the oxygen evolution reaction and membrane hydration. Larger pore sizes also improve oxygen gas transport from the catalyst layer, because generated oxygen gas is forced to travel in-plane through the anode catalyst layer until it reaches a pore opening that is connected to a channel. The discussed results confirm that the proposed thin foil based PTLs are fundamentally different from conventional PTLs, such as felts or layered meshes. The model developed in this work also provides generalizable insight into fundamental PEMEC phenomena, such as the competition between liquid and gas phase transport, membrane hydration and water management, and nonuniform electrochemical reactions, which are processes relevant to all PEMEC designs.

25 ENERGY STORAGE↗
DOE OSTI2022

Metabolic flux modeling of Gluconobacter oxydans enables improved production of bioleaching organic acids

Bioleaching with organic acids is an option for recovery of rare earth elements (REE) and or other valuable metals from consumer or industrial wastes. Biohydrometallurgy can play an economically and environmentally sustainable role in REE recycling, particularly if agriculture or food wastes are used as the carbon substrates for production of the bioleaching agents. A challenge for organic acid based REE bioleaching is the incomplete conversion of the carbon substrate to an optimal organic acid mixture with high REE extraction capability. To help meet this challenge a metabolic flux model was developed for Gluconobacter oxydans to identify options for improving sugar conversion to chelating organic acids. Here, model predictions were consistent with laboratory data for the microbial consumption of glucose and production of gluconic acid and suggested that providing growth conditions with relatively high carbon to nitrogen molar ratios could result in production of lower pH biolixiviant and thus improve REE bioleaching yields.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1995

The Dripping Handrail Model: Transient Chaos in Accretion Systems

We define and study a simple dynamical model for accretion systems, the "dripping handrail" (DHR). The time evolution of this spatially extended system is a mixture of periodic and apparently random (but actually deterministic) behavior. The nature of this mixture depends on the values of its physical parameters - the accretion rate, diffusion coefficient, and density threshold. The aperiodic component is a special kind of deterministic chaos called transient chaos. The model can simultaneously exhibit both the quasiperiodic oscillations (QPO) and very low frequency noise (VLFN) that characterize the power spectra of fluctuations of several classes of accretion systems in astronomy. For this reason, our model may be relevant to many such astrophysical systems, including binary stars with accretion onto a compact object - white dwarf, neutron star, or black hole - as well as active galactic nuclei. We describe the systematics of the DHR's temporal behavior, by exploring its physical parameter space using several diagnostics: power spectra, wavelet "scalegrams," and Lyapunov exponents. In addition, we note that for large accretion rates the DHR has periodic modes; the effective pulse shapes for these modes - evaluated by folding the time series at the known period - bear a resemblance to the similarly- determined shapes for some x-ray pulsars. The pulsing observed in some of these systems may be such periodic-mode accretion, and not due to pure rotation as in the standard pulsar model.

Young, Karl↗
DOE OSTI2021

LANL Critical Benchmark Comparison Study and Subsequent Revision for Cases Involving LEU and MIX

As part of an international collaboration within the DOE Nuclear Criticality Safety Program (NCSP), LANL is involved in a comparison study to quantify differences in k-effective results from neutron transport simulations of critical benchmark experiments. The DOE NCSP Mission and Vision details the activity in which the French Institut de Radioprotection et de Sûreté Nucléaire (IRSN) leads the study with LANL, ORNL, and LLNL to compare results of various neutron transport codes and nuclear data libraries for ICSBEP benchmarks held in common by the entities. The task statement from the DOE NCSP Five-Year Execution Plan is given as: "The proposal is for IRSN to lead a new intercomparison based on the MORET code with the latest JEFF-3.3 data and ENDF/B-VIII.0 data, when available, using their existing comprehensive selection of 2,714 benchmarks and collate their results together with those from LLNL (COG), LANL (MCNP) and ORNL (SCALE). Due to the large number of benchmarks involved, this effort is envisioned to take three years with an additional year for IRSN to complete a summary report. The benchmark development will be performed independently to minimize modeling errors through discovery and resolution of discrepant results. A summary report will be generated (led by IRSN) to document the results of this study." This report documents results obtained for revisions made to cases involving Intermediate Enriched Uranium (IEU) and a mixture of Pu and Uranium (MIX), with a focus on the changes made to LANL benchmarks modeled with MCNP6 using ENDF/B-VII.1 nuclear data that appeared to have discrepant results when compared with results of other codes. Feedback was used to pinpoint review of particular benchmark input files and to revise them when necessary. This report documents the results of review and revision of specific benchmarks highlighted as possibly discrepant in the comparison study. In addition, there is an effort tied to this work involving collaboration between LANL XCP and NCS Divisions in the development of a shared review/revision procedure and use of a new benchmark repository.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2014

NASA Numerical and Experimental Evaluation of UTRC Low Emissions Injector

Computational and experimental analyses of a PICS-Pilot-In-Can-Swirler technology injector, developed by United Technologies Research Center (UTRC) are presented. NASA has defined technology targets for near term (called "N+1", circa 2015), midterm ("N+2", circa 2020) and far term ("N+3", circa 2030) that specify realistic emissions and fuel efficiency goals for commercial aircraft. This injector has potential for application in an engine to meet the Pratt & Whitney N+3 supersonic cycle goals, or the subsonic N+2 engine cycle goals. Experimental methods were employed to investigate supersonic cruise points as well as select points of the subsonic cycle engine; cruise, approach, and idle with a slightly elevated inlet pressure. Experiments at NASA employed gas analysis and a suite of laser-based measurement techniques to characterize the combustor flow downstream from the PICS dump plane. Optical diagnostics employed for this work included Planar Laser-Induced Fluorescence of fuel for injector spray pattern and Spontaneous Raman Spectroscopy for relative species concentration of fuel and CO2. The work reported here used unheated (liquid) Jet-A fuel for all fuel circuits and cycle conditions. The initial tests performed by UTRC used vaporized Jet-A to simulate the expected supersonic cruise condition, which anticipated using fuel as a heat sink. Using the National Combustion Code a PICS-based combustor was modeled with liquid fuel at the supersonic cruise condition. All CFD models used a cubic non-linear k-epsilon turbulence wall functions model, and a semi-detailed Jet-A kinetic mechanism based on a surrogate fuel mixture. Two initial spray droplet size distribution and spray cone conditions were used: (1) an initial condition (Lefebvre) with an assumed Rosin-Rammler distribution, and 7 degree Solid Spray Cone; and (2) the Boundary Layer Stripping (BLS) primary atomization model giving the spray size distribution and directional properties. Contour and line plots are shown in comparison with experimental data (where this data is available) for flow velocities, fuel, and temperature distribution. The CFD results are consistent with experimental observations for fuel distribution and vaporization. Analysis of gas sample results, using a previously-developed NASA NOx correlation, indicates that for sea-level takeoff, the PICS configuration is predicted to deliver an EINOx value of about three for the targeted supersonic aircraft. Emissions results at supersonic cruise conditions show potential for meeting the NASA goals with liquid fuel.

Hicks, Yolanda R.↗
NASA NTRS2014

NASA Numerical and Experimental Evaluation of UTRC Low Emissions Injector

Computational and experimental analyses of a PICS-Pilot-In-Can-Swirler technology injector, developed by United Technologies Research Center (UTRC) are presented. NASA has defined technology targets for near term (called "N+1", circa 2015), midterm ("N+2", circa 2020) and far term ("N+3", circa 2030) that specify realistic emissions and fuel efficiency goals for commercial aircraft. This injector has potential for application in an engine to meet the Pratt & Whitney N+3 supersonic cycle goals, or the subsonic N+2 engine cycle goals. Experimental methods were employed to investigate supersonic cruise points as well as select points of the subsonic cycle engine; cruise, approach, and idle with a slightly elevated inlet pressure. Experiments at NASA employed gas analysis and a suite of laser-based measurement techniques to characterize the combustor flow downstream from the PICS dump plane. Optical diagnostics employed for this work included Planar Laser-Induced Fluorescence of fuel for injector spray pattern and Spontaneous Raman Spectroscopy for relative species concentration of fuel and CO2. The work reported here used unheated (liquid) Jet-A fuel for all fuel circuits and cycle conditions. The initial tests performed by UTRC used vaporized Jet-A to simulate the expected supersonic cruise condition, which anticipated using fuel as a heat sink. Using the National Combustion Code a PICS-based combustor was modeled with liquid fuel at the supersonic cruise condition. All CFD models used a cubic non-linear k-epsilon turbulence wall functions model, and a semi-detailed Jet-A kinetic mechanism based on a surrogate fuel mixture. Two initial spray droplet size distribution and spray cone conditions were used: 1) an initial condition (Lefebvre) with an assumed Rosin-Rammler distribution, and 7 degree Solid Spray Cone; and 2) the Boundary Layer Stripping (BLS) primary atomization model giving the spray size distribution and directional properties. Contour and line plots are shown in comparison with experimental data (where this data is available) for flow velocities, fuel, and temperature distribution. The CFD results are consistent with experimental observations for fuel distribution and vaporization. Analysis of gas sample results, using a previously-developed NASA NOx correlation, indicates that for sea-level takeoff, the PICS configuration is predicted to deliver an EINOx value of about 3 for the targeted supersonic aircraft. Emissions results at supersonic cruise conditions show potential for meeting the NASA goals with liquid fuel.

Hicks, Yolanda R.↗
DOE OSTI2021

Accelerating phase-field-based microstructure evolution predictions via surrogate models trained by machine learning methods

Abstract The phase-field method is a powerful and versatile computational approach for modeling the evolution of microstructures and associated properties for a wide variety of physical, chemical, and biological systems. However, existing high-fidelity phase-field models are inherently computationally expensive, requiring high-performance computing resources and sophisticated numerical integration schemes to achieve a useful degree of accuracy. In this paper, we present a computationally inexpensive, accurate, data-driven surrogate model that directly learns the microstructural evolution of targeted systems by combining phase-field and history-dependent machine-learning techniques. We integrate a statistically representative, low-dimensional description of the microstructure, obtained directly from phase-field simulations, with either a time-series multivariate adaptive regression splines autoregressive algorithm or a long short-term memory neural network. The neural-network-trained surrogate model shows the best performance and accurately predicts the nonlinear microstructure evolution of a two-phase mixture during spinodal decomposition in seconds, without the need for “on-the-fly” solutions of the phase-field equations of motion. We also show that the predictions from our machine-learned surrogate model can be fed directly as an input into a classical high-fidelity phase-field model in order to accelerate the high-fidelity phase-field simulations by leaping in time. Such machine-learned phase-field framework opens a promising path forward to use accelerated phase-field simulations for discovering, understanding, and predicting processing–microstructure–performance relationships.

36 MATERIALS SCIENCE↗
NASA NTRS2001

A Small Aircraft Transportation System (SATS) Demand Model

The Small Aircraft Transportation System (SATS) demand modeling is a tool that will be useful for decision-makers to analyze SATS demands in both airport and airspace. We constructed a series of models following the general top-down, modular principles in systems engineering. There are three principal models, SATS Airport Demand Model (SATS-ADM), SATS Flight Demand Model (SATS-FDM), and LMINET-SATS. SATS-ADM models SATS operations, by aircraft type, from the forecasts in fleet, configuration and performance, utilization, and traffic mixture. Given the SATS airport operations such as the ones generated by SATS-ADM, SATS-FDM constructs the SATS origin and destination (O&D) traffic flow based on the solution of the gravity model, from which it then generates SATS flights using the Monte Carlo simulation based on the departure time-of-day profile. LMINET-SATS, an extension of LMINET, models SATS demands at airspace and airport by all aircraft operations in US The models use parameters to provide the user with flexibility and ease of use to generate SATS demand for different scenarios. Several case studies are included to illustrate the use of the models, which are useful to identify the need for a new air traffic management system to cope with SATS.

Long, Dou↗
NASA NTRSDate not supplied

Cloud-Based Solutions for Monitoring Coastal Ecosystems and Prioritization of Restoration Efforts Across Belize

In recent years the availability of automated change detection algorithms in Google Earth Engine has permitted cloud-based processing of large time series of satellite imagery. Models such as the Continuous Change Detection and Classification (CCDC), CCDC-Spectral Mixture Analysis (CCDC-SMA), and Landsat-based Detection of Trends in Disturbance and Recovery (LandTrendr) allow users to exploit decades of Earth Observations (EO) , leveraging the Landsat archive and data from other sensors to detect disturbance in forest ecosystems. Despite the wide adoption of these methods, robust documentation and growing community of users, little research has explored their use in mangrove environments. Mangroves are dynamic environments subject to changes due to not only the natural migration of mudflats but also coastal erosion, urban expansion, aquaculture practices, etc. This work aims to identify best practices for the application of these models to identify and monitor changes in the Belizean mangroves, which experienced an estimated 5.4% decrease in extent between 1980 and 2017 (Cherrington et al. 2020). Partnering directly with the Belizean Forest Department, our team will develop a replicable, efficient methodology to annually update the country’s mangrove extent employing EO-based change detection. This collaboration supports Belize’s marine conservation and climate resilience efforts, such as its Blue Bond agreement with the Nature Conservancy, building on national processes to monitor these critical ecosystems.

Mangroves↗
DOE OSTI2021

Progress Towards a New Model Chemolithoautotrophic Host (CRADA Final Report)

Some of the largest scale chemical processes for the production of important commodity chemicals and fuels involve an H 2 -containing gas mixture as an intermediate [e.g. ammonia, urea, methanol, Fischer-Tropsch (F-T) diesel, and hydrotreatment of crude petroleum to refined fuels]. The successful development of a model H 2 -oxidizing chemoautotrophic host could expand the range of fuels and chemicals produced from H 2 and CO 2 intermediates in the chemical, oil, and gas sectors in the near term, as well as from renewable and waste-derived sources of these gases that are expected to greatly expand in coming years. Among non-photosynthetic bacteria that can utilize H 2 and CO 2 , Cupriavidus necator (C. necator, formerly Ralstonia eutropha), is the best studied. C. necator is an excellent microbial host for the production of a variety of chemicals because it grows extremely quickly to very high cell densities autotrophically on H 2 and CO 2 , is genetically tractable, and has the ability to accumulate polymers, such as polyhydroxybutyrate, at industrial levels. Despite having great potential as a platform bioproduction host, genetic tools are limited, making metabolic engineering of this organism slow and laborious. In this CRADA project, we developed a number of genetic tools for C. necator: 1) Improvement of C. necator genetic transformation efficiency 2) Integration of heterologous genes into the C. necator chromosome and development of promoter library 3) Development of graded RBSs to control heterologous protein expression 4) Use of RBSs to demonstrate fatty alcohol production in C. necator 5) Demonstration of CRISPR-Cas9 gene editing in C. necator

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1983

Lagrange thermodynamic potential and intrinsic variables for He-3 He-4 dilute solutions

For a two-fluid model of dilute solutions of He-3 in liquid He-4, a thermodynamic potential is constructed that provides a Lagrangian for deriving equations of motion by a variational procedure. This Lagrangian is defined for uniform velocity fields as a (negative) Legendre transform of total internal energy, and its primary independent variables, together with their thermodynamic conjugates, are identified. Here, similarities between relations in classical physics and quantum statistical mechanics serve as a guide for developing an alternate expression for this function that reveals its character as the difference between apparent kinetic energy and intrinsic internal energy. When the He-3 concentration in the mixtures tends to zero, this expression reduces to Zilsel's formula for the Lagrangian for pure liquid He-4. An investigation of properties of the intrinsic internal energy leads to the introduction of intrinsic chemical potentials along with other intrinsic variables for the mixtures. Explicit formulas for these variables are derived for a noninteracting elementary excitation model of the fluid. Using these formulas and others also derived from quantum statistical mechanics, another equivalent expression for the Lagrangian is generated.

Jackson, H. W.↗
DOE OSTI2022

Assessing the Role of the Support in Acetic Acid Hydrodeoxygenation Selectivity on Pt/TiO2

To produce drop-in quality biofuels following the thermochemical conversion of biomass via catalytic fast pyrolysis (CFP), hydrodeoxygenation (HDO) reactions can be performed to remove excess oxygen and create a more stable bio-oil product. HDO reactions involve co-feeding the CFP vapor-phase product and H2 gas over bi-functional catalysts, such as noble-metal catalysts supported on reducible metal oxides (e.g., Pt/TiO2). Recent model-compound studies, focusing on important classes of species in the CFP vapor mixture, have sought to determine the role of the various Pt/TiO2 actives sites (i.e., Pt-metal, TiO2-support, Pt-TiO2-interfacial sites) in producing the observed desired deoxygenation and undesired decarboxylation/decarbonylation products. One important class of compounds in the CFP vapor-phase product is carboxylic acids (e.g., acetic acid). Prior experimental work on acetic acid HDO (AA-HDO) found that Pt/C and Pt/TiO2 catalysts favored the formation of undesired C-C and desired C-O bond-dissociation products; however, the fundamental surface chemistry driving this shift in selectivity has not been established. This presentation employs atomic-scale modeling to determine, through comparisons of adsorption and reaction energetics, how Pt-metal, Pt-TiO2-interface, and interfacial-vacancy sites catalyze the competing reaction pathways for AA-HDO. Our analysis indicates that hydroxyl vacancies at the Pt-TiO2 interface are critical for lowering barriers for C-O bond-cleavage steps, such that they become favorable/competitive with respect to C-C bond-dissociation steps.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗
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↗