Search NASA⌕ Search

SEARCH · Search NASA

Results for “Arrhenius equation”

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 55 records · Page 3

Non-intrusive data-driven model reduction for differential–algebraic equations derived from lifting transformations

In this paper we present a non-intrusive data-driven approach for model reduction of nonlinear systems. The approach considers the particular case of nonlinear partial differential equations (PDEs) that form systems of partial differential–algebraic equations (PDAEs) when lifted to polynomial form. Such systems arise, for example, when the governing equations include Arrhenius reaction terms (e.g., in reacting flow models) and thermodynamic terms (e.g., the Helmholtz free energy terms in a phase-field solidification model). Using the known structured form of the lifted algebraic equations, the approach computes the reduced operators for the algebraic equations explicitly, using straightforward linear algebra operations on the basis matrices. The reduced operators for the differential equations are inferred from lifted snapshot data using operator inference, which solves a linear least squares regression problem. The approach is illustrated for the nonlinear model of solidification of a pure material. The lifting transformations reformulate the solidification PDEs as a system of PDAEs that have cubic structure. The operators of the lifted system for this solidification example have affine dependence on key process parameters, permitting us to learn a parametric reduced model with operator inference. Numerical experiments show the effectiveness of the resulting reduced models in capturing key aspects of the solidification dynamics.

42 ENGINEERING↗

Numerical modeling of chemical vapor deposition (CVD) in a horizontal reactor

In the present numerical prediction of the deposition rate of silicon from silane in a CVD process, the conservation equations for mass, momentum, energy, and chemical species are solved on a staggered grid using the SIMPLE algorithm, while the rate of chemical reactions in the gas phase and on the susceptor surface is obtained from an Arrhenius rate equation. Predicted deposition rates as a function of position along the susceptor with and without the gas phase chemical reaction are compared with the available experimental and numerical data; agreement is excellent except at the leading edge of the susceptor, where the deposition rate is overpredicted.

Sheikholeslami, M. Z.↗

Modeling of outgassing and matrix decomposition in carbon-phenolic composites

A new release rate equation to model the phase change of water to steam in composite materials was derived from the theory of molecular diffusion and equilibrium moisture concentration. The new model is dependent on internal pressure, the microstructure of the voids and channels in the composite materials, and the diffusion properties of the matrix material. Hence, it is more fundamental and accurate than the empirical Arrhenius rate equation currently in use. The model was mathematically formalized and integrated into the thermostructural analysis code CHAR. Parametric studies on variation of several parameters have been done. Comparisons to Arrhenius and straight-line models show that the new model produces physically realistic results under all conditions.

Mcmanus, Hugh L.↗

Permeation Rate Equations for Hydrogen and Deuterium in a Palladium-Silver Alloy

Mass transfer of a gas through a selective, solid membrane is an effective method for separation of desired species. This selective permeability is evident in the flow of hydrogen and the isotope deuterium through palladium-silver metal alloy media. In this study, based upon Sieverts’s law and the Arrhenius diffusion equation, an empirical correlation was developed to determine the steady-state permeation rate R, dependence on media temperature T, gas supply pressure p(sub S), and gas backpressure p(sub B) on the lower pressure side. Because of an extensive range of experimental conditions and complete reporting of raw data, the research by Ackerman and Koskinas was used as a source for data allowing empirical equation fitting. Unfortunately, those authors reported best-fit equations that poorly represented their own results. To improve the modeling of the original data and demonstrate the quality of the measurements, the current study develops improved hydrogen and deuterium permeation rate equations: P = 4.22×10(exp −6) A[exp(−704/T)](sq. root p(sub S) − sq. root p(sub B))/t for hydrogen P = 2.12×10(exp −6) A[exp(−468/T )](sq. root p(sub S) − sq. root p(sub B))/t for deuterium for values of cross-sectional area A (sq.cm), medium thickness t (cm), pressure p (psia), and temperature T (Kelvin), giving a permeation rate in mole/minute. These equations model permeation rate data more closely than do several other existing literature sources.

Smith, Phillip J.↗

A Model for the Oxidation of Carbon Silicon Carbide Composite Structures

A mathematical theory and an accompanying numerical scheme have been developed for predicting the oxidation behavior of carbon silicon carbide (C/SiC) composite structures. The theory is derived from the mechanics of the flow of ideal gases through a porous solid. The result of the theoretical formulation is a set of two coupled nonlinear differential equations written in terms of the oxidant and oxide partial pressures. The differential equations are solved simultaneously to obtain the partial vapor pressures of the oxidant and oxides as a function of the spatial location and time. The local rate of carbon oxidation is determined using the map of the local oxidant partial vapor pressure along with the Arrhenius rate equation. The nonlinear differential equations are cast into matrix equations by applying the Bubnov-Galerkin weighted residual method, allowing for the solution of the differential equations numerically. The numerical method is demonstrated by utilizing the method to model the carbon oxidation and weight loss behavior of C/SiC specimens during thermogravimetric experiments. The numerical method is used to study the physics of carbon oxidation in carbon silicon carbide composites.

Sullivan, Roy M.↗

Oxidation Behavior of Carbon Fiber-Reinforced Composites

OXIMAP is a numerical (FEA-based) solution tool capable of calculating the carbon fiber and fiber coating oxidation patterns within any arbitrarily shaped carbon silicon carbide composite structure as a function of time, temperature, and the environmental oxygen partial pressure. The mathematical formulation is derived from the mechanics of the flow of ideal gases through a chemically reacting, porous solid. The result of the formulation is a set of two coupled, non-linear differential equations written in terms of the oxidant and oxide partial pressures. The differential equations are solved simultaneously to obtain the partial vapor pressures of the oxidant and oxides as a function of the spatial location and time. The local rate of carbon oxidation is determined at each time step using the map of the local oxidant partial vapor pressure along with the Arrhenius rate equation. The non-linear differential equations are cast into matrix equations by applying the Bubnov-Galerkin weighted residual finite element method, allowing for the solution of the differential equations numerically.

Sullivan, Roy M.↗

A Model for the Oxidation of C/SiC Composite Structures

A mathematical theory and an accompanying numerical scheme have been developed for predicting the oxidation behavior of C/SiC composite structures. The theory is derived from the mechanics of the flow of ideal gases through a porous solid. Within the mathematical formulation, two diffusion mechanisms are possible: (1) the relative diffusion of one species with respect to the mixture, which is concentration gradient driven and (2) the diffusion associated with the average velocity of the gas mixture, which is total gas pressure gradient driven. The result of the theoretical formulation is a set of two coupled nonlinear differential equations written in terms of the oxidant and oxide partial pressures. The differential equations must be solved simultaneously to obtain the partial vapor pressures of the oxidant and oxides as a function of space and time. The local rate of carbon oxidation is determined as a function of space and time using the map of the local oxidant partial vapor pressure along with the Arrhenius rate equation. The nonlinear differential equations are cast into matrix equations by applying the Bubnov-Galerkin weighted residual method, allowing for the solution of the differential equations numerically. The end result is a numerical scheme capable of determining the variation of the local carbon oxidation rates as a function of space and time for any arbitrary C/SiC composite structures.

Sullivan, Roy M.↗

The influence of temperature on area-specific impedance and capacity of Li-ion cells with nickel-containing positive electrodes

The temperature-dependent behaviors of five nickel-containing positive electrodes (NCA, NMC811, NMC622, NMC532, and NMC111) in lithium-ion batteries are investigated using an electrochemical protocol involving rate studies, mild aging (~100 cycles), and hybrid pulse power characterization (HPPC). Tests are conducted using coin-cells with graphite negative electrodes at -20 °C, 0 °C, 20 °C, and 40 °C. Three techniques are compared for determining the area-specific impedance (ASI): i) fits to the rate study average voltages, ii) fitting to the entire voltage curves using a regularization scheme, and iii) HPPC. When fit to an Arrhenius-type equation, all methods yield similar apparent activation energies (±2 kJ/mol) for the impedance, which range from -20 to -31 kJ/mol for the electrodes. Impedance growth increases with temperature but remains at less than 0.2% per cycle for most electrodes and temperatures. NCA and NMC811 are the exceptions, which yield 0.5% and 1.5% increases in ASI per cycle, respectively, at 40 °C. For cells with the same electrodes, the capacities are similar at 20 and 40 °C but reduce at lower temperatures, with up to a 52% reduction at -20 °C and 2C. The fade in energy of the cells during C/3 cycling is attributed to decreasing capacity as opposed to increasing ASI.

25 ENERGY STORAGE↗

On the Arrhenius-like behavior of conductivity during flash sintering of 3 mol% yttria stabilized zirconia ceramics

Here, we discuss the Arrhenius behavior of electrical conductivity σ during flash sintering of 3 mol% yttria stabilized zirconia (3YSZ). In situ x-ray diffraction is used to determine sample temperature. Sintering contribution to σ is excluded by investigating the flash event on a dense ceramic. We show that total conductivity follows an Arrhenius-like equation in both dense and green samples, even during Stage II of flash. The non-linearity often verified during flash sintering of 3YSZ is therefore related to the furnace temperature instead of the sample temperature when building the ln(σ) vs. 1/T plots. Furthermore, we verified a change in the activation energy for conduction prior to the ignition of the flash event and discussed the possible mechanisms. For instance, the decrease in activation energy from Stage I to II in the dense sample is attributed to a contribution from electronic carriers.

36 MATERIALS SCIENCE↗

Late Stage Domain Coarsening Dynamics of Lamellar Block Copolymers

Block copolymer (BCP) phase separation dynamics can be expected to differ significantly from that of the polymer blends due to the constraint of chain connectivity. BCP phase separation dynamics has been studied theoretically, but there has been little experimental evidence to confirm the BCP domain growth scaling laws put forward theoretically. Here, we investigate the dynamics of late-stage lamellar BCP domain coarsening and find that the scaling exponent for lamellar domain growth is ≈1/6 (0.17), irrespective of the annealing temperature, a value close to the scaling exponent of 0.2 predicted by theoretical studies. Furthermore, we show that the prefactors in the domain coarsening equation show Arrhenius dependence on temperature, indicating that the BCP domain growth dynamics is Arrhenius over the temperature range investigated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Assessing models of force-dependent unbinding rates via infrequent metadynamics

Protein–ligand interactions are crucial for a wide range of physiological processes. Many cellular functions result in these non-covalent “bonds” being mechanically strained, and this can be integral to proper cellular function. Broadly, two classes of force dependence have been observed—slip bonds, where the unbinding rate increases, and catch bonds, where the unbinding rate decreases. Despite much theoretical work, we cannot predict for which protein–ligand pairs, pulling coordinates, and forces a particular rate dependence will appear. Here, we assess the ability of MD simulations combined with enhanced sampling techniques to probe the force dependence of unbinding rates. We show that the infrequent metadynamics technique correctly produces both catch and slip bonding kinetics for model potentials. We then apply it to the well-studied case of a buckyball in a hydrophobic cavity, which appears to exhibit an ideal slip bond. Finally, we compute the force-dependent unbinding rate of biotin–streptavidin. Here, the complex nature of the unbinding process causes the infrequent metadynamics method to begin to break down due to the presence of unbinding intermediates, despite the use of a previously optimized sampling coordinate. Allowing for this limitation, a combination of kinetic and free energy computations predicts an overall slip bond for larger forces consistent with prior experimental results although there are substantial deviations at small forces that require further investigation. Furthermore, this work demonstrates the promise of predicting force-dependent unbinding rates using enhanced sampling MD techniques while also revealing the methodological barriers that must be overcome to tackle more complex targets in the future.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A New Analytical Expression for Estimating the Adiabatic Temperature Rise in Lithium-Ion Batteries During High-Power Pulses

A coupled, thermal-electrochemical model is used to explain why the Joule heating assumption (i.e., I 2 R) does not provide a good representation of the temperature rise during high power pulses in lithium-ion batteries, even in cases where the reversible heat generation can be neglected. The poor agreement occurs because the internal resistance changes during the pulse due to the opposing effects of mass transport in the electrolyte, which raises the resistance, and heat generation (temperature rise), which lowers the resistance. These insights are used to propose a new analytical expression for predicting the temperature rise during adiabatic pulses with limited experimental and physical data. The expression accounts for thermal effects using a Taylor series expansion of an Arrhenius-type equation. It accounts for transport effects using an approximate solution to the one-dimensional diffusion equation. The new expression is shown to accurately estimate the simulated, adiabatic temperature rise across a range of loadings (1 to 4 mAh cm -2 ) and C-rates (1C to 10C) for cells containing LiNi 0.5 Mn 0.3 Co 0.2 O 2 positive electrodes and graphite negative electrodes. In conclusion, it is also shown to accurately estimate the experimental, adiabatic temperature rise measured for cells with LiFePO 4 positive electrodes and mesocarbon negative electrodes with minimal changes to the fitted parameters.

25 ENERGY STORAGE↗

Flue Gas Desulfurization (FGD) Wastewater Treatment Using Polybenzimidazole (PBI) Hollow Fiber (HF) Membranes

Polybenzimidazole (PBI) hollow fiber membranes were used to treat flue gas desulfurization (FGD) wastewater (WW) from a coal fired power plant. Membranes were tested using both single salt solutions and real FGD WW. The PBI membranes showed >99% rejection for single salt solutions of NaCl, MgCl2, CaSO4, and CaCl2 at approximately 2000 PPM (parts per million). The membranes also showed >97% rejection for FGD WW concentrations ranging from 6900 to 14,400 PPM total dissolved solids (TDS). The pH of the FGD WW was adjusted between 3.97–8.20, and there was an optimal pH between 5.31 and 7.80 where %rejection reached a maximum of >99%. The membranes were able to operate stably up to 50 °C, nearly doubling the water flux as compared to room temperature, and while maintaining >98% salt rejection.

activation energy↗

Yielding and deformation behavior of the single crystal nickel-base superalloy PWA 1480

Interrupted tensile tests were conducted to fixed plastic strain levels in 100 ordered single crystals of the nickel based superalloy PWA 1480. Testing was done in the range of 20 to 1093 C, at strain rate of 0.5 and 50%/min. The yield strength was constant from 20 to 760 C, above which the strength dropped rapidly and became a stong function of strain rate. The high temperature data were represented very well by an Arrhenius type equation, which resulted in three distinct temperature regimes. The deformation substructures were grouped in the same three regimes, indicating that there was a fundamental relationship between the deformation mechanisms and activation energies. Models of the yielding process were considered, and it was found that no currently available model was fully applicable to this alloy. It was also demonstrated that the initial deformation mechanism (during yielding) was frequently different from that which would be inferred by examining specimens which were tested to failure.

Milligan, W. W., Jr.↗

Preliminary study: Moisture-polymer interaction. Stuby objectives

The problems associated with mathematically modeling water-module interaction phenomena, including sorption and desorption, diffusion, and permeation are discussed. With reliable analytical models, an extensive materials data base, and solar radiation surface meteorological observations (SOLMET) weather data, predicting module lifetimes in realistic environments can become a practical reality. The status of the present techniques of simulating the various transport mechanisms was reported. The Dent model (a modified Brunauer-Emmet-Teller) approach represented polyvinyl butyral (PVB) sorption data. A 100-layer material model and Fick's diffusion model gave diffusivity values exhibiting adequate agreement with those measured for PVB. Diffusivity of PVB is concentration dependent, decreasing as the water content in PVB increases. The temperature dependence of diffusion in PVB is well modeled by the Arrhenius rate equation. Equilibrium conductivity and leakage current data are well represented by Hearle's model for bulk ionic conductivity. A nodal network analysis using the Systems Improved Numerical Differencing Analyzer (SINDA) Thermal Analyzer gave reasonable correlation with measurable data. It is concluded that realistic lifetime predictions seem to be feasible.

Wen, L. C.↗

Yielding and deformation behavior of the single crystal superalloy PWA 1480

Interrupted tensile tests were conducted to fixed plastic strain levels in 100 ordered single crystals of the nickel based superalloy PWA 1480. Testing was done in the range of 20 to 1093 C, at strain rate of 0.5 and 50 percent/min. The yield strength was constant from 20 to 760 C, above which the strength dropped rapidly and became a strong function of strain rate. The high temperature data were represented very well by an Arrhenius type equation, which resulted in three distinct temperature regimes. The deformation substructures were grouped in the same three regimes, indicating that there was a fundamental relationship between the deformation mechanisms and activation energies. Models of the yielding process were considered, and it was found that no currently available model was fully applicable to this alloy. It was also demonstrated that the initial deformation mechanism (during yielding) was frequently different from that which would be inferred by examining specimens which were tested to failure.

Milligan, Walter W.↗

Evaluation of transport conditions during physical vapor transport growth of opto-electronic crystals

Transport conditions were evaluated during the vapor phase growth of mercurous chloride crystals in a closed tube. Experimentally observed growth rates were much smaller than those calculated by the Hertz-Knudsen (H-K) equation. The Arrhenius behavior of growth rate with the temperature was used to derive the sticking coefficient. A one-dimensional diffusion model was used to calculate the total mass flux and was compared with the condensing flux. It was predicted that growth occurred in the convecto-diffusive range.

Singh, N. B.↗

A classical equation that accounts for observations of non-Arrhenius and cryogenic grain boundary migration

Observations of microstructural coarsening at cryogenic temperatures, as well as numerous simulations of grain boundary motion that show faster migration at low temperature than at high temperature, have been troubling because they do not follow the expected Arrhenius behavior. This work demonstrates that classical equations, that are not simplified, account for all these oddities and demonstrate that non-Arrhenius behavior can emerge from thermally activated processes. According to this classical model, this occurs when the intrinsic barrier energies of the processes become small, allowing activation at cryogenic temperatures. Additional thermal energy then allows the low energy process to proceed in reverse, so increasing temperature only serves to frustrate the forward motion. This classical form is shown to reconcile and describe a variety of diverse grain boundary migration observations.

36 MATERIALS SCIENCE↗