Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic parameters”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A Potts Model parameter study of particle size, Monte Carlo temperature, and “Particle-Assisted Abnormal Grain Growth”

A Potts Model was proposed to account for temperature and particle heterogeneity-dependent Zener Pinning. This was accomplished by relating the random fluctuations of grain boundary position allowed by the Potts Model switching probability at finite simulation temperature to experimentally observed grain growth stagnation behavior. We assume that these fluctuations arise from random fluctuations in the thermodynamic and kinetic properties of the grain boundaries and/or interfaces as they change with temperature. As an application of this model, the grain growth kinetics of U-10 wt. % Mo nuclear fuels were simulated with the input of microstructural images during different heat treatment processes. Simulated average grain growth behavior is in good agreement with experiments.

Frazier, William E.↗
DOE OSTI2020

Mapping the kinetic evolution of metastable grain boundaries under non-equilibrium processing

The kinetic evolution of a multiplicity of metastable grain boundaries (GBs) under fast driving conditions are studied by atomistic modeling. Assisted with an enhanced statistical analysis, the energetic evolution of GBs over a broad metastability-temperature space is mapped out, wherein two distinct regimes—an ageing regime and a rejuvenating regime—are retrieved with high fidelity. By comparing the results under various conditions (e.g. random perturbations, isothermal annealing, and fast heating/cooling), it is shown that such ageing/rejuvenating mechanism map is universal, irrespective of the actual stimuli used to elicit the metastable GBs. Here, the ageing/rejuvenating phenomena are demonstrated to stem from the energy imbalance of uphill climbing and downhill dropping during sequential transitions in the system's potential energy landscape. Without the necessity of introducing free parameters, such model can reconcile experimentally measured hardness variation of nanocrystalline metals subjected to femto-second laser irradiation, and it therefore provides a novel perspective on achieving a plurality of interfacial states and facilitating previously inaccessible property regimes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Dependence of Solar Wind Proton Temperature on the Polarization Properties of Alfvénic Fluctuations at Ion-kinetic Scales

In this work, we use fluctuating magnetic helicity to investigate the polarization properties of Alfvénic fluctuations at ion-kinetic scales in the solar wind as a function of β p , the ratio of proton thermal pressure to magnetic pressure, and θ vB , the angle between the proton flow and local mean magnetic field, B 0 . Using almost 15 yr of Wind observations, we separate the contributions to helicity from fluctuations with wavevectors, k, quasi-parallel and oblique to B 0 , finding that the helicity of Alfvénic fluctuations is consistent with predictions from linear Vlasov theory. This result suggests that the nonlinear turbulent fluctuations at these scales share at least some polarization properties with Alfvén waves. Furthermore, we also investigate the dependence of proton temperature in the β p –θ vB plane to probe for possible signatures of turbulent dissipation, finding that it correlates with θ vB . The proton temperature parallel to B 0 is higher in the parameter space where we measure the helicity of right-handed Alfvénic fluctuations, and the temperature perpendicular to B 0 is higher where we measure left-handed fluctuations. This finding is inconsistent with the general assumption that by sampling different θ vB in the solar wind we can analyze the dependence of the turbulence distribution on θ kB , the angle between k and B 0 . After ruling out both instrumental and expansion effects, we conclude that our results provide new evidence for the importance of local kinetic processes that depend on θ vB in determining proton temperature in the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Proton affinities of pertechnetate (TcO 4 - ) and perrhenate (ReO 4 - )

The anions pertechnetate, TcO 4 - , and perrhenate, ReO 4 - , exhibit very similar chemical and physical properties. Revealing and understanding disparities between them enhances fundamental understanding of both. Electrospray ionization generated the gas-phase proton bound dimer (TcO 4 - )(H + )(ReO 4 - ). Collision induced dissociation of the dimer yielded predominantly HTcO 4 - and ReO 4 - , which according to Cooks’ kinetic method indicates that the proton affinity (PA) of TcO 4 - is greater than that of Re O 4 - . Density functional theory computations agree with the experimental observation, providing PA[TcO 4 - ] = 300.1 kcal mol -1 and PA[ReO 4 - ] = 297.2 kcal mol -1 . Attempts to rationalize these relative PAs based on elementary molecular parameters such as atomic charges indicate that the entirety of bond formation and concomitant bond disruption needs to be considered to understand the energies associated with such protonation processes. Although in both the gas and solution phases, TcO 4 - is a stronger base than ReO 4 - , it is noted that the significance of even such qualitative accordance is tempered by the very different natures of the underlying phenomena.

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

Towards fully predictive gyrokinetic full- f simulations: validation and triangularity studies in TCV

Designing economical magnetic confinement fusion power plants motivates computational tools that can estimate plasma behavior from engineering parameters without direct reliance on experimental measurement of the plasma profiles. In this work, we present full-f global long-wavelength gyrokinetic simulations of edge and scrape-off layer turbulence in tokamaks that use only magnetic geometry, heating power, and particle inventory as inputs. Unlike many modeling approaches that employ free parameters fitted to experimental data, raising uncertainties when extrapolating to reactor scales. This approach directly simulates turbulence and resulting profiles through gyrokinetics without such empirical adjustments. This is achieved via an adaptive sourcing algorithm in Gkeyll that strictly controls energy injection and emulates particle sourcing due to neutral recycling. We show that the simulated kinetic profiles compare reasonably well with Thomson scattering and Langmuir probe data for Tokamak á Configuration Variable (TCV) discharge #65125, and that the simulations reproduce characteristic features such as blob transport and self-organized electric fields. Applying the same framework to study triangularity effects suggests mechanisms contributing to the improved confinement reported for negative triangularity (NT). Simulations of TCV discharges #65125 and #65130 indicate that NT increases the E x B flow shear (by about 20% in these cases), which correlates with reduced turbulent losses and a modest change in the distribution of power exhaust to the vessel wall. While the physical models contain approximations that can be refined in future work, the predictive capability demonstrated here, evolving multiple profile relaxation times with kinetic electron and ion models in hundreds of GPU hours, indicates the feasibility of using Gkeyll to support design studies of fusion devices.

Hoffmann, Antoine Cyril David [Princeton Plasma Ph↗
DOE OSTI2021

Analytic quasi-steady evolution of marginally unstable kinetic plasma instabilities in the presence of drag and scattering

The 1D bump-on-tail problem is studied in order to determine the influence of drag on quasi-steady solutions near marginal stability (1-γ d /γ L << 1) when effective collisions are much larger than the instability growth rate (ν>>γ). In this common tokamak regime, it is rigorously shown that the paradigmatic Berk–Breizman cubic equation for the nonlinear mode evolution reduces to a much simpler differential equation, dubbed the time-local cubic equation, which can be solved directly. It is found that in addition to increasing the saturation amplitude, drag introduces a shift in the apparent oscillation frequency by modulating the saturated wave envelope. Furthermore, excellent agreement is found between the analytic solution for the mode evolution and both the numerically integrated Berk–Breizman cubic equation and fully nonlinear 1D Vlasov simulations. Experimentally isolating the contribution of drag to the saturated mode amplitude for verification purposes is explored but complicated by the reality that the amount of drag cannot be varied independently of other key parameters in realistic scenarios. While the influence of drag is modest when the ratio of effective drag to effective scattering α/ν is very small, it can become substantial when α/ν ≳ 0.5, suggesting that drag should be accounted for in quantitative models of fast-ion-driven instabilities in fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Computational Prediction of Coiled–Coil Protein Gelation Dynamics and Structure

Protein hydrogels represent an important and growing biomaterial for a multitude of applications, including diagnostics and drug delivery. We have previously explored the ability to engineer the thermoresponsive supramolecular assembly of coiled–coil proteins into hydrogels with varying gelation properties, where we have defined important parameters in the coiled–coil hydrogel design. Using Rosetta energy scores and Poisson–Boltzmann electrostatic energies, we iterate a computational design strategy to predict the gelation of coiled–coil proteins while simultaneously exploring five new coiled–coil protein hydrogel sequences. Provided this library, we explore the impact of in silico energies on structure and gelation kinetics, where we also reveal a range of blue autofluorescence that enables hydrogel disassembly and recovery. As a result of this library, we identify the new coiled–coil hydrogel sequence, Q5, capable of gelation within 24 h at 4 °C, a more than 2-fold increase over that of our previous iteration Q2. The fast gelation time of Q5 enables the assessment of structural transition in real time using small-angle X-ray scattering (SAXS) that is correlated to coarse-grained and atomistic molecular dynamics simulations revealing the supramolecular assembling behavior of coiled–coils toward nanofiber assembly and gelation. This work represents the first system of hydrogels with predictable self-assembly, autofluorescent capability, and a molecular model of coiled–coil fiber formation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Assessment of the literature about Be-W mixed material layer formation in the fusion reactor environment

Abstract All plasma facing surfaces in a fusion reactor, whether initially pure or an alloy, will rapidly evolve into a mixed material due to plasma-induced erosion, migration and redeposition. Beryllium (Be) erosion from the main chamber, and its transport and deposition on to a tungsten (W) divertor results in the growth of mixed Be-W layers, which can evolve to form beryllides. These Be-W mixed materials exhibit generally less desirable properties than pure tungsten or pure beryllium, such as lower melting points. In order to better understand the parameter space for growth of these alloys, this paper reviews the literature on Be-W mixed material formation experiments—in magnetically confined fusion reactors, in linear plasma test stands, and during thin-film deposition—and on computational modeling of Be-W interactions, as well as briefly assesses the Be-W growth kinetics. We conclude that the following kinetic steps drive the material mixing: adsorption of the implanted/deposited ion on the metal surface; diffusion of the implanted/deposited ion from surface into the bulk, which is accelerated by defects; and loss of deposited material through erosion. Adsorption dominates (or prevents) material mixing in thin-film deposition experiments, whereas diffusion drives material mixing in plasma exposures due to the energetic ion implantation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Assessment of the literature about Be-W mixed material layer formation in the fusion reactor environment

All plasma facing surfaces in a fusion reactor, whether initially pure or an alloy, will rapidly evolve into a mixed material due to plasma-induced erosion, migration and redeposition. Beryllium (Be) erosion from the main chamber, and its transport and deposition on to a tungsten (W) divertor results in the growth of mixed Be-W layers, which can evolve to form beryllides. These Be-W mixed materials exhibit generally less desirable properties than pure tungsten or pure beryllium, such as lower melting points. In order to better understand the parameter space for growth of these alloys, this paper reviews the literature on Be-W mixed material formation experiments—in magnetically confined fusion reactors, in linear plasma test stands, and during thin-film deposition—and on computational modeling of Be-W interactions, as well as briefly assesses the Be-W growth kinetics. We conclude that the following kinetic steps drive the material mixing: adsorption of the implanted/deposited ion on the metal surface; diffusion of the implanted/deposited ion from surface into the bulk, which is accelerated by defects; and loss of deposited material through erosion. Adsorption dominates (or prevents) material mixing in thin-film deposition experiments, whereas diffusion drives material mixing in plasma exposures due to the energetic ion implantation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

MULTI-DIMENSIONAL MODELING OF THE CFR ENGINE FOR THE INVESTIGATION OF SI NATURAL GAS COMBUSTION AND CONTROLLED END-GAS AUTOIGNITION

Engine knock and misfire are barriers to pathways leading to high-efficiency Spark-Ignited (SI) Natural Gas engines. The general tendency to knock is highly dependent on engine operating conditions and the fuel reactivity. The problem is further complicated by low emission limits and the wide range of chemical reactivity in pipeline quality natural gas. Depending on the region and the source of the natural gas, its reactivity, described by its methane number (analogous to the octane number for liquid SI fuels) can span from 65 - 95. In order to realize diesel-like efficiencies, SI natural gas engines must be designed to operate at high BMEP near knock limits over a wide range of fuel reactivity. This requires a deep understanding regarding the combustion-engine interactions pertaining to flame propagation and end-gas autoignition (EGAI). However, EGAI, if controlled, provides an opportunity to increase SI natural gas engine efficiency by increasing combustion rate and the total burned fuel, mitigating the effects of the slow flame speeds of natural gas fuels which generally reduce BMEP and increase unburned hydrocarbon emissions. For this reason, in order to study EGAI phenomenon, the present work highlights multi-dimensional computational fluid dynamics (CFD) models of the Cooperative Fuel Research (CFR) engine. The CFR engine models are used to investigate fuel-engine interactions that lead to EGAI with natural gas, including effects of fuel reactivity, engine operating parameters, and exhaust gas recirculation (EGR). A Three-Pressure Analysis, performed with GT-Power, was used to estimate initial and boundary conditions for the three-dimensional CFD model. CONVERGE CFD v2.4 was used for the three-dimensional CFD modeling where the level set G-Equation model and SAGE detailed chemical kinetics solver were used. An assessment of the different modeling approaches is also provided to evaluate their limitations, advantages and disadvantages, and for which situations they are most applicable. Model validation was performed with experimental data taken with a CFR engine over varying compression ratio, CA50, EGR fraction, and IMEP and shows good agreement in Peak Cylinder Pressure (PCP), PCP crank angle, and the location of the 10%, 50%, and 90% mass fraction burned (CA10, CA50, and CA90, respectively). The models can predict the onset crank angle and pressure rise rate for light, medium, and heavy EGAI under a variety of fuel reactivities and engine operating conditions.

Bestel, Diego↗
DOE OSTI2022

Predicting Asteroid Material Properties from a DART-like Kinetic Impact

Abstract NASA’s Double Asteroid Redirection Test (DART) mission is the first full-scale test of the kinetic impactor method for asteroid deflection, in which a spacecraft intentionally impacts an asteroid to change its trajectory. DART represents an important first step for planetary defense technology demonstration, providing a realistic assessment of the effectiveness of the kinetic impact approach on a near-Earth asteroid. The momentum imparted to the asteroid is transferred from the impacting spacecraft and enhanced by the momentum of material ejected from the impact site. However, the magnitude of the ejecta contribution is dependent on the material properties of the target. These properties, such as strength and shear modulus, are unknown for the DART target asteroid, Dimorphos, as well as most asteroids since such properties are difficult to characterize remotely. This study examines how hydrocode simulations can be used to estimate material properties from information available post-impact, specifically the asteroid size and shape, the velocity and properties of the impacting spacecraft, and the final velocity change imparted to the asteroid. Across >300 three-dimensional simulations varying seven material parameters describing the asteroid, we found many combinations of properties could reproduce a particular asteroid velocity. Additional observations, such as asteroid mass or crater size, are required to further constrain properties like asteroid strength or outcomes like the momentum enhancement provided by impact ejecta. Our results demonstrate the vital importance of having as much knowledge as possible prior to an impact mission, with key material parameters being the asteroid’s mass, porosity, strength, and elastic properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Pyrolysis of high-density polyethylene: Degradation behaviors, kinetics, and product characteristics

Pyrolysis is a promising technology for converting plastic waste into valuable raw materials while offering a potential solution to the global plastic pollution crisis. In this study, the thermal pyrolysis of high-density polyethylene (HDPE) is investigated in a drop tube reactor under nearly isothermal conditions. The impact of reaction temperature and gas/volatile residence time on carbon conversion and product distribution is examined across a range of 500–900°C and 3.6–32.2s, respectively. Non-condensable gas products detected by online mass spectrometry are H 2 , CH 4 , C 2 H 4 , C 2 H 6 , C 3 H 6 , and C 3 H 8 . At elevated temperatures and prolonged residence time, H 2 yield reaches as high as 8.6 wt% of the initial HDPE mass due to intensified cracking reactions of C 2 –C 3 hydrocarbons and long-chain aliphatic compounds. Consequently, pyrolysis tars consist mainly of polycyclic aromatic hydrocarbons (PAHs) with 5–7 rings, accompanied by visible coke deposition within the reactor. HDPE decomposition to volatiles is an endothermic process and it is complete at a temperature between 492°C and 525°C, depending on the heating rate employed, from non-isothermal thermogravimetric analysis and differential scanning calorimetry (TGA-DSC) measurements. The thermal degradation of HDPE pellets follows the two-dimensional nucleation growth model for conversion levels up to 0.8 with an apparent activation energy of 259–270 kJ/mol and a pre-exponential factor of 4.83 × 10 17 –1.37 × 10 19 min -1 , determined from various isoconversional methods such as Flynn-Wall-Ozawa (FWO), Kissinger-Akahira-Sunose (KAS), and Starink, along with Criado's master plots. Further, these findings provide valuable insights into optimizing process parameters and refining reactor design for pyrolysis, which can be integrated with gasification and reforming processes to enhance hydrogen production on a larger scale.

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

Lepto-axiogenesis

We propose a baryogenenesis mechanism that uses a rotating condensate of a Peccei-Quinn (PQ) symmetry breaking field and the dimension-five operator that gives Majorana neutrino masses. The rotation induces charge asymmetries for the Higgs boson and for lepton chirality through sphaleron processes and Yukawa interactions. The dimension-five interaction transfers these asymmetries to the lepton asymmetry, which in turn is transferred into the baryon asymmetry through the electroweak sphaleron process. QCD axion dark matter can be simultaneously produced by dynamics of the same PQ field via kinetic misalignment or parametric resonance, favoring an axion decay constant ƒ a ≲ 10 10 GeV, or by conventional misalignment and contributions from strings and domain walls with ƒ a ~ 10 11 GeV. The size of the baryon asymmetry is tied to the mass of the PQ field. In simple supersymmetric theories, it is independent of UV parameters and predicts the supersymmtry breaking mass scale to be $ \mathcal{O} $(10 - 10 4 ) TeV, depending on the masses of the neutrinos and whether the condensate is thermalized during a radiation or matter dominated era. The high supersymmetry breaking mass scale may be free from cosmological and flavor/CP problems. We also construct a theory where TeV scale supersymmetry is possible. Parametric resonance may give warm axions, and the radial component of the PQ field may give signals in rare kaon decays from mixing with the Higgs and in dark radiation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Atomic Layer Deposition Niobium Nitride Films for High-Q Resonators

Niobium nitride (NbN) is a useful material for fabricating detectors because of its high critical temperature and relatively high kinetic inductance. In particular, NbN can be used to fabricate nanowire detectors and mm-wave transmission lines. When deposited, NbN is usually sputtered, leaving room for concern about uniformity at small thicknesses. In this work, we present atomic layer deposition niobium nitride (ALD NbN) as an alternative technique that allows for precision control of deposition parameters such as film thickness, stage temperature, and nitrogen composition. Atomic-scale control over film thickness admits a high degree of uniformity for films 4–30 nm thick; control over deposition temperature gives rise to growth rate changes, which can be used to optimize film thickness and critical temperature. In order to charac- terize ALD NbN in the radio-frequency regime, we construct single-layer micro- wave resonators and test their performance as a function of stage temperature and input power. ALD processes can admit high resonator quality factors, with ≥ 43% of resonators above Q i = 10 6 , which in turn increase detector multiplexing capabilities. i Furthermore, we find critical temperatures in the range of 7.5–10.9 K that vary as a function of cycle count and deposition temperature.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Kinetic theory of particle-in-cell simulation plasma and the ensemble averaging technique

Abstract We derive the kinetic theory of fluctuations in physically and numerically stable particle-in-cell (PIC) simulations of electrostatic plasmas. The starting point is the single-time correlation at the start of the simulation between the statistical fluctuations of the weighted densities of macroparticle centers in the plasma particle phase-space. The fluctuations are associated with different initial conditions, typically due to the random initial conditions (in velocity space) of the macroparticles/simulation plasma, assigned according to their initial distribution of probability. The single-time correlations at all time steps and in each spatial grid cell are then determined from the Laplace–Fourier transforms of the discretized Klimontovich-like equation for the macroparticles and Maxwell’s equations for the fields, as computed by modern PIC codes. We recover the expressions for the electrostatic field and the plasma particle density fluctuation autocorrelation spectra as well as the kinetic equations describing the average evolution of PIC-simulated plasma particles, first derived by Langdon (1970b Proc. 4th Conf. Numerical Simulation of Plasmas ) using a test macroparticle approach perturbing a discretized Vlasovian plasma and then averaging the obtained physical quantity over the initial macroparticle velocity distribution. We generalize and extend these results to the modern algorithms in PIC codes using arbitrary macroparticle weights. Analytical estimates of statistical fluctuation amplitudes are derived as a function of the plasma simulation parameters, using the central limit theorem in the limit of a large number of macroparticles per cell. The theory is then used to analyze the ensemble averaging technique of PIC simulations where statistical averages are performed over ensembles of PIC simulations, modeling the same plasma physics problem but using different statistical realizations of the initial distribution functions of the macroparticles. This method is illustrated by linear Landau damping uncovering (from noise, which is usually considered numerical) the physical fluctuations driven by a single small amplitude electrostatic wave perturbing a PIC simulation plasma in equilibrium.

fluctuations correlations↗
DOE OSTI2020

High-yield growth kinetics and spatial mapping of single-walled carbon nanotube forests at wafer scale

Emerging commercial applications of vertically aligned, single-walled carbon nanotube (SWCNT) “forests” require synthesis that minimizes nanotube diameter while maximizing number density across substrate areas exceeding centimeter scale. To address this need, we synthesized SWCNT forests on full silicon wafers with notable reproducibility and uniformity, and co-optimized growth for small diameters and high densities across large areas to access new territory in this 3D parameter space. We mapped the spatial uniformity of key structural features using Raman microscopy, synchrotron X-ray scattering, and Rutherford backscattering spectrometry. Low C 2 H 2 flux over sub-nm Fe/Mo catalysts produced small-diameter SWCNTs (2.1 nm) at high number densities (2.26 × 10 12 cm -2 ) on wafers up to 6 in. Although removing Mo resulted in larger SWCNT diameters and lower densities (<0.7 × 10 12 cm -2 ), mass conversion rates from C 2 H 2 to SWCNT product were high and remarkably invariant for catalyst compositions and densities (i.e., 47.7% or 1.30 × 10 6 % g-catalyst -1 on 4-in. wafers). These carbon conversion efficiencies far exceed typical benchtop reactors and are on par with the best reported literature values. Our detailed elucidation of correlations among structural characteristics within this resource-efficient process is expected to guide future scale-up efforts of SWCNT forest growth beyond wafer scale.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Leveraging Temperature-Dependent (Electro)Chemical Kinetics for High-Throughput Flow Battery Characterization

The library of redox-active organics that are potential candidates for electrochemical energy storage in flow batteries is exceedingly vast, necessitating high-throughput characterization of molecular lifetimes. Demonstrated extremely stable chemistries require accurate yet rapid cell cycling tests, a demand often frustrated by time-denominated capacity fade mechanisms. We have developed a high-throughput setup for elevated temperature cycling of redox flow batteries, providing a new dimension in characterization parameter space to explore. We utilize it to evaluate capacity fade rates of aqueous redox-active organic molecules, as functions of temperature. We demonstrate Arrhenius-like behavior in the temporal capacity fade rates of multiple flow battery electrolytes, permitting extrapolation to lower operating temperatures. Collectively, these results highlight the importance of accelerated decomposition protocols to expedite the screening process of candidate molecules for long lifetime flow batteries.

25 ENERGY STORAGE↗
DOE OSTI2020

Cloud based tool for analysis of chemical kinetic mechanisms

The size and complexity of chemical kinetic models used to simulate combustion of hydrocarbon fuels continue to increase with our improved understanding of the underlying physical mechanisms, along with the improved ability to utilize such models as the result of increased computational power and efficiency. As mechanisms grow beyond thousands of species and tens of thousands of reactions, it becomes increasingly difficult to manually check for errors and physical inconsistencies. In this work, we present several automated methods to check for such issues in the specification of chemical reaction models. First, we demonstrate how discontinuities in thermodynamic data can cause simulation difficulties manifested as long wall-clock times and convergence failures. To correct this type of problem, we describe an automated method for refitting thermodynamic parameters. We also outline several methods to check the timescales of reaction rate coefficients to ensure physical consistency. All the methods are made available through a web-based tool (https://combustiontools.llnl.gov) intended to aid mechanism developers and users to improve the accuracy and performance of fuel models used by the combustion community.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗