Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reduce order 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,027 records · Page 57

Recycling of Lead Pastes from Spent Lead–Acid Batteries: Thermodynamic Constraints for Desulphurization

Lead–acid batteries are important to modern society because of their wide usage and low cost. The primary source for production of new lead–acid batteries is from recycling spent lead–acid batteries. In spent lead–acid batteries, lead is primarily present as lead pastes. In lead pastes, the dominant component is lead sulfate (PbSO 4 , mineral name anglesite) and lead oxide sulfate (PbO•PbSO 4 , mineral name lanarkite), which accounts for more than 60% of lead pastes. In the recycling process for lead–acid batteries, the desulphurization of lead sulfate is the key part to the overall process. In this work, the thermodynamic constraints for desulphurization via the hydrometallurgical route for recycling lead pastes are presented. The thermodynamic constraints are established according to the thermodynamic model that is applicable and important to recycling of lead pastes via hydrometallurgical routes in high ionic strength solutions that are expected to be in industrial processes. The thermodynamic database is based on the Pitzer equations for calculations of activity coefficients of aqueous species. The desulphurization of lead sulfates represented by PbSO 4 can be achieved through the following routes. (1) conversion to lead oxalate in oxalate-bearing solutions; (2) conversion to lead monoxide in alkaline solutions; and (3) conversion to lead carbonate in carbonate solutions. Among the above three routes, the conversion to lead oxalate is environmentally friendly and has a strong thermodynamic driving force. Oxalate-bearing solutions such as oxalic acid and potassium oxalate solutions will provide high activities of oxalate that are many orders of magnitude higher than those required for conversion of anglesite or lanarkite to lead oxalate, in accordance with the thermodynamic model established for the oxalate system. An additional advantage of the oxalate conversion route is that no additional reductant is needed to reduce lead dioxide to lead oxide or lead sulfate, as there is a strong thermodynamic force to convert lead dioxide directly to lead oxalate. As lanarkite is an important sulfate-bearing phase in lead pastes, this study evaluates the solubility constant for lanarkite regarding the following reaction, based on the solubility data, PbO•PbSO 4 + 2H + ⇌ 2Pb 2+ + SO 4 2– + H 2 O(l).

25 ENERGY STORAGE↗

Model predictive control of mixing controlled compression ignition operation for low reactivity fuels

Using gasoline or other low reactivity fuels with a pilot injection or port fuel injection in a compression ignition engine has shown great potential in reducing NOx emissions while keeping high thermal efficiency compared to diesel. However, excessive combustion noise is caused by a high maximum pressure rise rate in the cylinder due to the higher fractions of premixed charge of the low-reactivity fuel. This noise can result in structural damage to engine components and as such, combustion noise limits the range of the operating parameters and makes the control of such engines challenging. In this study, a simulation environment was built up in MATLAB/Simulink leveraging a physics-based zero-dimension combustion model to capture the in-cylinder pressure time traces as well as metrics relevant to thermal efficiency and combustion noise. Here, in order to also facilitate the control of emissions, machine learning models were investigated to capture NOx emissions. A kernel-based extreme learning machine (K-ELM) performed best and had a coefficient of correlation (R-squared) of 0.998. The combustion and NOx emission models are valid for not only conventional gasoline fuel but also oxygenated alternative fuel blends at three different pilot injection strategies. In order to track key combustion metrics while keeping noise and emissions within constraints, a model predictive control (MPC) was applied for a compression ignition engine operating with a range of potential fuels and fuel injection strategies. The MPC is validated under different scenarios, including a load step change, fuel type change, and injection strategy change, with proportional–integral (PI) control as the baseline. The simulation results show that MPC reduces about 26% of ringing intensity in the transient process and 17% at the steady state for E30. Generally, MPC can optimize the overall performance through modifying the main injection timing, pilot fuel mass, and exhaust gas recirculation (EGR) fraction.

42 ENGINEERING↗

Benchmarking MCNP models to support gamma irradiation experiments

Seismic isolators and damping devices have been used for 30+ years to reduce seismic forces and deployed in 10,000+ structures, including bridges, buildings, and mission-critical infrastructure. Studies have shown that seismic isolation of nuclear power plants can reduce the seismic risk by several orders of magnitude and substantially reduce the overnight capital cost. Advanced reactor developers are considering seismic isolation as an integral design feature in their reactor design. It can either be deployed at the building level, isolating the entire reactor building, or at the component level, isolating individual pieces of safety-class equipment, such as a reactor vessel or a steam generator. For equipment isolation, located close to the reactor vessel, the seismic isolators and/or damping devices may be exposed to gamma and neutron radiation, with a possible effect on their mechanical properties which would affect the response of the isolated equipment under earthquake shaking. The Department of Energy, as part of its Nuclear Energy University Program (NEUP), in collaboration with the Idaho National Laboratory (INL), is funding a project at the University at Buffalo (UB) to investigate the effect of gamma radiation on the mechanical properties of various seismic protective devices. (When isolating equipment, the seismic protective devices will be located outside the reactor vessel where the neutron exposure will be insignificant.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

An In Situ , Automated High-Explosives Aging Method Utilizing Two-Dimensional Gas Chromatography–Mass Spectrometry

Understanding chemical changes that occur in high explosives as they age is of great importance to the safe employment and storage of these compounds. Traditional methods of aging high explosives even under accelerated aging conditions are time intensive with durations on the order of months to years. The nature of traditional aging analyses reduces each sample to a snapshot data point often separated widely in time, requiring many assumptions as to how the degradation products develop. Further complicating matters, several analytical techniques are typically employed for each sample analysis in order to ascertain an entire picture of the decomposition pathways. To address these shortcomings with existing methods, a new method of accelerated aging of high explosives utilizing comprehensive two-dimensional gas chromatography coupled to high-resolution mass spectrometry (GC × GC-HRMS) was developed using 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (CL-20) as a model compound for method development. This in situ automated method reduces the time scale of aging to a matter of hours using the inlet of the GC × GC as the aging vessel. GC × GC in combination with HRMS allowed for the collection of both evolved gases and other decomposition products produced during the entire aging process in real time with HRMS providing far greater certainty in identification of explosives aging products. Additionally, this method allowed for a higher throughput of samples with greatly simplified sample preparation. Chemometric analysis of the GC × GC-HRMS data set via the alteration analysis (ALA) enabled discovery of statistically significant chemical changes providing insight into the variation of decomposition pathways with varying aging temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling cross-beam energy transfer with sector ray tracing

Ray-based cross-beam energy transfer (CBET) models are an essential feature of the radiation-hydrodynamic codes used to simulate inertial confinement fusion implosions, but full 3D ray-based CBET calculations can have a prohibitively high computational cost. Sector ray tracing can be used to reduce the cost by orders of magnitude in cases where the coronal plasma and laser drive can be approximated as spherically symmetric. An extension of sector ray tracing (section ray tracing) can be used to relax the assumption of a spherically symmetric drive while still retaining most of the computational savings of a sector ray trace. We discuss the foundations of sector and section ray tracing and compare them to full ray tracing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Climatology of Cloud-Top Radiative Cooling in Marine Shallow Clouds

A one-year’s worth of near-global marine shallow single-layer cloud top radiative cooling (CTRC) is derived from a radiative transfer model with inputs from the satellite cloud retrievals and reanalysis sounding. The mean cloud top radiative flux divergence is 61 Wm -2 , decomposed into the longwave and shortwave components of 73 and -11 W m -2 , respectively. Equatorward of 30°N/S, the CTRC is largely a reflection of free-atmospheric specific humidity distribution: a dry atmosphere enhances CTRC by reducing downward thermal radiation. Consequently, the cooling minimizes in the “wet” tropics and maximizes in the “dry” eastern subtropics. Poleward of 30°N/S, the CTRC decreases slightly due to the colder clouds that emit less effectively. Furthermore, the CTRC exhibits distinctive seasonal cycles with stronger cooling in the winter and has amplitudes of order 10–20 Wm -2 in stratocumulus-rich regions. The datasets were used to train a machine-learning model that substantially speeds up the retrieval.

54 ENVIRONMENTAL SCIENCES↗

Slow manifold reduction for plasma science

The classical Chapman–Enskog procedure admits a substantial geometrical generalization known as slow manifold reduction. This generalization provides a paradigm for deriving and understanding most reduced models in plasma physics that are based on controlled approximations applied to problems with multiple timescales. In this Review we develop the theory of slow manifold reduction with a plasma physics audience in mind. In particular we illustrate (a) how the slow manifold concept may be used to understand breakdown of a reduced model over sufficiently-long time intervals, and (b) how a discrete-time analogue of slow manifold theory provides a useful framework for developing implicit integrators for temporally-stiff plasma models. For readers with more advanced mathematical training we also use slow manifold reduction to explain the phenomenon of inheritance of Hamiltonian structure in dissipation-free reduced plasma models. Various facets of the theory are illustrated in the context of the Abraham–Lorentz model of a single charged particle experiencing its own radiation drag. As a culminating example we derive the slow manifold underlying kinetic quasineutral plasma dynamics up to first-order in perturbation theory. This first-order result incorporates several physical effects associated with small deviations from exact charge neutrality that lead to slow drift away from predictions based on the leading-order approximation n e = Z i n i .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

PSCF+: An Extended and Improved Open-Source Software Package for Polymer Self-Consistent Field Calculations

This work introduces PSCF+, a freely available open‑source software tool for modeling how block copolymers self‑assemble into ordered nanoscale structures. The software allows researchers to study a wide range of polymer models and interaction types, making it possible to directly compare theoretical predictions with molecular and simulation‑based studies in a consistent way. PSCF+ incorporates several recently developed numerical methods that significantly reduce computer memory requirements and speed up calculations, enabling efficient use of modern GPU hardware. Here, it also provides automated tools to calculate free‑energy landscapes and phase boundaries, which are essential for understanding and predicting material behavior. Together, these advances make PSCF+ a powerful and flexible platform for studying polymer self‑assembly and support the design of new polymeric materials for scientific and technological applications.

Chemical calculations↗

Tuning perovskite nanocrystal superlattices for superradiance in the presence of disorder

The cooperative emission of interacting nanocrystals is an exciting topic fueled by recent reports of superfluorescence and superradiance in assemblies of perovskite nanocubes. Several studies estimated that coherent coupling is localized to a small fraction of nanocrystals (10 -7 –10 -3 ) within the assembly, raising questions about the origins of localization and ways to overcome it. In this work, we examine single-excitation superradiance by calculating radiative decays and the distribution of superradiant wave function in two-dimensional CsPbBr 3 nanocube superlattices. The calculations reveal that the energy disorder caused by size distribution and large interparticle separations reduces radiative coupling and leads to the excitation localization, with the energy disorder being the dominant factor. The single-excitation model clearly predicts that, in the pursuit of cooperative effects, having identical nanocubes in the superlattice is more important than achieving a perfect spatial order. The monolayers of large CsPbBr 3 nanocubes (L NC = 10–20 nm) are proposed as model systems for experimental tests of superradiance under conditions of non-negligible size dispersion, while small nanocubes (L NC = 5–10 nm) are preferred for realizing the Dicke state under ideal conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Twisted bilayer graphene. IV. Exact insulator ground states and phase diagram

Here, we derive the exact insulator ground states of the projected Hamiltonian of magic-angle twisted bilayer graphene (TBG) flat bands with Coulomb interactions in various limits, and study the perturbations away from these limits. We define the (first) chiral limit where the AA stacking hopping is zero, and a flat limit with exactly flat bands. In the chiral-flat limit, the TBG Hamiltonian has a U(4) × U(4) symmetry, and we find that the exact ground states at integer filling – 4 ≤ $\textit{ν}$ ≤ 4 relative to charge neutrality are Chern insulators of Chern numbers $ν_C$ = 4 – |$\textit{ν}$ |, 2 – |$\textit{ν}$ |, $\cdots$, |$\textit{ν}$ | – 4, all of which are degenerate. This confirms recent experiments where Chern insulators are found to be competitive low-energy states of TBG. When the chiral-flat limit is reduced to the nonchiral-flat limit which has a U(4) symmetry, we find $\textit{ν}$ = 0 , ± 2 has exact ground states of Chern number 0, while $\textit{ν}$ = ±1, ± 3 has perturbative ground states of Chern number $ν_C$ = ± 1, which are U(4) ferromagnetic. In the chiral-nonflat limit with a different U(4) symmetry, different Chern number states are degenerate up to second-order perturbations. In the realistic nonchiral-nonflat case, we find that the perturbative insulator states with Chern number $ν_C$ = 0 (0 < |$ν_C$| < 4 – |$\textit{ν}$ |) at integer fillings $\textit{ν}$ are fully (partially) intervalley coherent, while the insulator states with Chern number |$ν_C$| = 4 – |$\textit{ν}$ | are valley polarized. However, for 0 < |$ν_C$| ≤ 4 – |$\textit{ν}$ |, the fully intervalley coherent states are highly competitive (0.005 meV/electron higher). At nonzero magnetic field |$\textit{B}$ | > 0, a first-order phase transition for $\textit{ν}$ = ± 1, ± 2 from Chern number $ν_C$ = sgn ($\textit{νB}$)(2 – |$\textit{ν}$ |) to $ν_C$ = sgn ($\textit{νB}$)(4 – |$\textit{ν}$ |) is expected, which agrees with recent experimental observations. Lastly, the TBG Hamiltonian reduces into an extended Hubbard model in the stabilizer code limit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Economic Storage Size Optimization for Electric Vehicle Extreme-Fast Charging Stations

En-route charging infrastructure for electric vehicles is critical to support transportation needs. These charging stations are likely to have high loads and especially sharp peak loads given fast charging capabilities needed to meet transportation schedules. In order to reduce both strain on distribution grid infrastructure and charging station operational costs, many stations are likely to employ behind the meter storage. This paper demonstrates a behind the meter storage sizing optimization that employs an open-source agent-based vehicle behavior model (BEAM) to determine the best sizing across many scenarios. This optimization and analysis is novel in that it examines how storage size impacts not only charging station cost and peak load, but also vehicle queue times. The optimization is also applied across a wide analysis region with sufficient diversity and numbers to provide novel statistical analysis of optimal sizes.

Aka, Julius↗

Design of a Single Layer Metamaterial for Pressure Offloading of Transtibial Amputees

While using a prosthesis, transtibial amputees may experience pain and discomfort brought on by large pressure gradients at the interface between the residual limb and prosthetic socket. Current prosthetic interface solutions attempt to alleviate these pressure gradients by using soft homogenous liners to reduce and distribute pressures. This research investigates an additively manufactured metamaterial inlay with a tailored mechanical response to reduce peak pressure gradients around the limb. The inlay uses a hyperelastic behaving metamaterial (US10244818) comprised of triangular pattern unit cells which can be 3D printed with walls of various thicknesses controlled by draft angles. The hyperelastic material properties are modeled using a Yeoh 3rd order model. The 3rd order coefficients can be adjusted and optimized, which corresponds to a change in the unit cell wall thickness to create an inlay that can meet the unique offloading needs of an amputee. Finite element analysis simulations evaluated the pressure gradient reduction from: 1) a common homogenous silicone liner, 2) a prosthetist's inlay prescription that utilizes three variations of the metamaterial, and 3) a metamaterial solution with optimized Yeoh 3rd order coefficients. When compared to a traditional homogenous silicone liner for two unique limb loading scenarios, the prosthetist prescribed inlay and the optimized material inlay can achieve equal or greater pressure gradient reduction capabilities. These results show the potential feasibility of implementing this metamaterial as a method of personalized medicine for transtibial amputees by creating a customizable interface solution to meet the unique performance needs of an individual patient.

60 APPLIED LIFE SCIENCES↗

Models and Strategies for Optimal Demand Side Management in the Chemical Industries

Deregulation and the increase of renewable electricity generation from wind and solar photovoltaics have transformed the U.S. electricity market. Economic and environmental benefits notwithstanding, the presence of renewables has increased variability and uncertainty on the supply side of the grid. Managing demand, rather than generation – a strategy referred to as “demand response (DR)” – is an attractive approach for mitigating this imbalance. DR efforts aim to reduce electricity usage during peak demand times, lessening stress on the grid. Industrial users are particularly attractive entities for DR participation since they present large, localized loads that can provide significant relief on grid demand and –unlike other large loads, such as buildings – are minimally dependent on human needs and preferences. In this project, we accomplished three main objectives. (1) We developed data-driven low-order DR scheduling-relevant dynamic models of chemical processes. Concurrently, we studied the formulation and solution of the associated optimal DR production scheduling problems. (a) A prototype air separation unit (ASU) model was used to generate simulated operating data for initial modeling efforts, which enabled the later use of industrial data for data-driven modeling. (b) We utilized Hammerstein-Wiener (HW) and Finite Step Response (FSR) models to represent nonlinear plant dynamics. (c) The HW models were linearized using exact linearization so they could potentially be embedded in power system models, which are formulated as mixed integer linear programs (MILPs). (d) We solved DR optimization problems under uncertainty and found that even naïve predictions of electricity price and product demand led to significant cost savings benefits. (2) Our DR scheduling optimization problem formulations are amenable to real-time solution. (a) We utilized Lagrangian Relaxation (LR) to efficiently solve the optimization problem by decoupling subproblems linked by complicating constraints. (b) We have achieved computation times for the 3-day DR scheduling problem of an ASU as low as 1.88 minutes. (3) Our representations of the DR behavior of chemical process as grid-level batteries were embedded in power system models. (a) For a small-scale grid, we found that incorporating the dynamics of the chemical plant in the optimal power flow calculations resulted in better resource management leading to up to 15% and 46% cost reduction for the grid and chemical plant operations, respectively, during periods of power line congestion. We have published several works dedicated to modeling and solving DR optimization problems from the user side. These were published in top peer-reviewed journals and are summarized in this report. The most recent work (and papers in preparation) considers DR scheduling from the grid side. Future efforts will consider networked plants (e.g., air separation units operating on a common pipeline) for DR participation, which is expected to amplify the capabilities of industrial DR participants to perform load-shifting. Our consideration of uncertainty in DR has inspired future directions in this area as well: we plan to develop multistage methods to fully account for the effects of uncertainty in DR scheduling.

24 POWER TRANSMISSION AND DISTRIBUTION↗

The PI3K/mTOR inhibitor Gedatolisib eliminates dormant breast cancer cells in organotypic culture, but fails to prevent metastasis in preclinical settings

Dormant, disseminated tumor cells (DTCs) are thought to be the source of breast cancer metastases several years or even decades after initial treatment. To date, a selective therapy that leads to their elimination has not been discovered. While dormant DTCs resist chemotherapy, evidence suggests that this resistance is driven not by their lack of proliferation, but by their engagement of the surrounding microenvironment, via integrin-β1-mediated interactions. Because integrin-β1-targeted agents have not been translated readily to the clinic, signaling nodes downstream of integrin-β1 could serve as attractive therapeutic targets in order to sensitize dormant DTCs to therapy. By probing a number of kinases downstream of integrin-β1, we determined that PI3K inhibition with either a tool compounds or a compound (PF-05212384; aka Gedatolisib) in clinical trials robustly sensitizes quiescent breast tumor cells seeded in organotypic bone marrow cultures to chemotherapy. These results motivated the preclinical study of whether Gedatolisib—with or without genotoxic therapy—would reduce DTC burden and prevent metastases. Despite promising results in organotypic culture, Gedatolisib failed to reduce DTC burden or delay, reduce or prevent metastasis in murine models of either triple-negative or estrogen receptor-positive breast cancer dissemination and metastasis. This result held true whether analyzing Gedatolisib on its own (vs. vehicle-treated animals) or in combination with dose-dense doxorubicin and cyclophosphamide (vs. animals treated only with dose-dense chemotherapies). These data suggest that PI3K is not the node downstream of integrin-β1 that confers chemotherapeutic resistance to DTCs. More broadly, they cast doubt on the strategy to target PI3K in order to eliminate DTCs and prevent breast cancer metastasis.

60 APPLIED LIFE SCIENCES↗

Field demonstration and implementation analysis of model predictive control in an office HVAC system

Model Predictive Control (MPC) is a promising technique to address growing needs for heating, ventilation, and air-conditioning (HVAC) systems to operate more efficiently and with greater flexibility. However, due to a number of factors, including the required implementation expertise, lack of high quality data, and a risk-adverse industry, MPC has yet to gain widespread adoption. While many previous studies have shown the advantages of MPC, few analyzed the implementation effort and associated practical challenges. In addition, previous work has developed an open-source, Modelica-based tool-chain that automatically generates optimal control, parameter estimation, and state estimation problems aimed at facilitating MPC implementation. Therefore, this study demonstrates usage of this tool-chain to implement MPC in a real office building, discusses practical challenges of implementing MPC, and estimates the implementation effort associated with various tasks in order to inform the development of future workflows and serve as an initial benchmark for their impact on reducing implementation effort. This study finds that the implemented MPC saves approximately 40% of HVAC energy over the existing control during a two-month trial period and that tasks related to data collection and controller deployment activities can each require as much effort as model generation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Measurement of free chlorine levels in water using potentiometric responses of biofilms and applications for monitoring and managing the quality of potable water

We report residual free chlorine is not monitored continuously at scale in drinking water distribution systems because existing real-time sensor technologies require frequent maintenance, cleaning, and calibration, which makes these products too costly to be used throughout a distribution system. As a result, current measurement approaches require manual sampling, which is not feasible for the consistent monitoring of free chlorine because chlorine concentrations vary significantly throughout pipeline distribution and over time and space. This research presents an alternative and cost-effective method of predicting free chlorine levels in drinking water using graphite electrodes coated with naturally grown microbial biofilms. This Microbial Potentiometric Sensor (MPS) array was installed in a Continuously Mixed Batch Reactor (CMBR), and drinking water containing variable free chlorine concentrations. The chlorine concentrations were introduced in a controlled manner, and the MPS signals were monitored over time. MPS signals were measured from the change in Open Circuit Potential (OCP) across the MPS array in real-time. An empirically derived relationship between the normalized change in OCP and free chlorine was established by fitting individual and average MPS data to a decaying exponential growth function in order to predict free chlorine levels. The results show that free chlorine can be predicted with reasonable accuracy, with model validation showing an average absolute error of ±0.09 ppm below 1.1 ppm and ±0.30 ppm between 1.1 and 2.7 ppm. However, the accuracy of predictions was reduced at higher free chlorine levels. The researchers conclude that MPS systems may benefit drinking water distribution systems by measuring free chlorine. These advantages of the MPS are especially pronounced in the developing world because this system is inexpensive and does not require routine maintenance or cleaning. The system relies on a naturally forming and regenerating biofilm and an inexpensive potentiometric meter to produce stable measurements.

54 ENVIRONMENTAL SCIENCES↗

Radiation Enhanced Anion Diffusion in Chromia

Cr 2 O 3 is a ubiquitous passivating film on alloys that is vital to their stability and performance in hostile environments. In this study, we directly observe and measure how anion transport through the oxide is affected by 400 keV Ar 2+ ion irradiation as a function of temperature (30–500 °C) and dose (0.33–2.00 dpa) in single-crystal Cr 2 O 3 films using embedded 18 O isotopic tracers and atom probe tomography. Diffusion coefficients are experimentally determined and compared to a chemical rate theory model. Experimental and theoretical results are broadly in agreement, showing anion diffusion increases by at least 4–5 orders of magnitude upon irradiation with a weak dependence on temperature and a stronger dose-driven sink-strength dependence. These results reveal that radiation could significantly reduce the protectiveness of Cr 2 O 3 films when passivity relies upon limited anion diffusivity, even at relatively modest levels of irradiation and low temperatures.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Microkinetic insights into the role of catalyst and water activity on the nucleation, growth, and dissolution during COF-5 synthesis

The chemical pathway for synthesizing covalent organic frameworks (COFs) involves a complex medley of reaction sequences over a rippling energy landscape that cannot be adequately described using existing theories. Even with the development of state-of-the-art experimental and computational tools, identifying primary mechanisms of nucleation and growth of COFs remains elusive. Other than empirically, little is known about how the catalyst composition and water activity affect the kinetics of the reaction pathway. In this work, for the first time, we employ time-resolved in situ Fourier transform infrared spectroscopy (FT-IR) coupled with a six-parameter microkinetic model consisting of ~10 million reactions and over 20 000 species. The integrated approach elucidates previously unrecognized roles of catalyst pK a on COF yield and water on growth rate and size distribution. COF crystalline yield increases with decreasing pK a of the catalysts, whereas the effect of water is to reduce the growth rate of COF and broaden the size distribution. The microkinetic model reproduces the experimental data and quantitatively predicts the role of synthesis conditions such as temperature, catalyst, and precursor concentration on the nucleation and growth rates. Furthermore, the model also validates the second-order reaction mechanism of COF-5 and predicts the activation barriers for classical and non-classical growth of COF-5 crystals. The microkinetic model developed here is generalizable to different COFs and other multicomponent systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗