Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Parameter calibration”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

High‐Resolution National‐Scale Water Modeling Is Enhanced by Multiscale Differentiable Physics‐Informed Machine Learning

Abstract The National Water Model (NWM) is a key tool for flood forecasting, planning, and water management. Key challenges facing the NWM include calibration and parameter regionalization when confronted with big data. We present two novel versions of high‐resolution (∼37 km 2 ) differentiable models (a type of hybrid model): one with implicit, unit‐hydrograph‐style routing and another with explicit Muskingum‐Cunge routing in the river network. The former predicts streamflow at basin outlets whereas the latter presents a discretized product that seamlessly covers rivers in the conterminous United States (CONUS). Both versions use neural networks to provide a multiscale parameterization and process‐based equations to provide a structural backbone, which were trained simultaneously (“end‐to‐end”) on 2,807 basins across the CONUS and evaluated on 4,997 basins. Both versions show great potential to elevate future NWM performance for extensively calibrated as well as ungauged sites: the median daily Nash‐Sutcliffe efficiency of all 4,997 basins is improved to around 0.68 from 0.48 of NWM3.0. As they resolve spatial heterogeneity, both versions greatly improved simulations in the western CONUS and also in the Prairie Pothole Region, a long‐standing modeling challenge. The Muskingum‐Cunge version further improved performance for basins >10,000 km 2 . Overall, our results show how neural‐network‐based parameterizations can improve NWM performance for providing operational flood predictions while maintaining interpretability and multivariate outputs. The modeling system supports the Basic Model Interface (BMI), which allows seamless integration with the next‐generation NWM. We also provide a CONUS‐scale hydrologic data set for further evaluation and use.

Song, Yalan [Civil and Environmental Engineering T↗
DOE OSTI2020

Active- and transfer-learning applied to microscale-macroscale coupling to simulate viscoelastic flows

Active- and transfer-learning are applied to microscale dynamics of polymer flows for the multiscale discovery of effective constitutive approximations required in viscoelastic flow simulation. The result is macroscopic rheology directly connected to a microstructural model. Micro and macroscale simulations are adaptively coupled by means of Gaussian process regression (GPR) to run the expensive microscale computations only as necessary. This multiscale method is demonstrated with flows of a polymer solution as a model system. At the microscale level dissipative particle dynamics (DPD) is employed to model the fluid as a suspension of bead-spring micro-structures subjected to steady shear flow. The results yield the non-Newtonian viscosity and the first normal stress difference at strain rates as training data used in a GPR model. DPD parameters are calibrated with respect to experimental data for a real polymer solution. Compliance with these data requires adjustment of the DPD model's cutoff radius, which then becomes a function of the second invariant of the strain rate tensor. The FENE-P model is chosen for the macroscale description using the spectral element method (SEM) to simulate channel flow and flow past a circular cylinder. The DPD results at the lowest possible shear strain rate yield an estimate of the zero-shear rate viscosity, which allows the initiation of the macroscale flow by SEM as a Newtonian fluid. The resulting strain-rate field is surveyed to determine additional shear strain rate sampling points for the DPD system. This new information allows an initial fitting of parameters of the constitutive equation followed by new SEM simulations at the macroscale. Additionally, guided by active-learning GPR to select new sampling points, this process continues until convergence is achieved. The effectiveness of this new simulation paradigm for viscoelastic flows is tested with different macroscale operating conditions. The effective closure learned in the channel simulation is then transferred directly to the flow past a circular cylinder at low Reynolds number, where the results show that only two additional DPD simulations are required to achieve a satisfactory constitutive model. With an increase of the Reynolds number, the active-learning scheme automatically detects the inaccuracy of the learned constitutive model, and initiates additional DPD simulations for the extra data needed to once again close the microscale-macroscale coupled system. This new paradigm of active- and transfer-learning for multiscale modeling is readily applicable to other microscale-macroscale coupled simulations of complex fluids and other materials. Furthermore, the coupling between microscale and macroscale solvers can be seamlessly implemented with our open source multiscale universal interface (MUI) library.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1993

Streak camera based SLR receiver for two color atmospheric measurements

To realize accurate two-color differential measurements, an image digitizing system with variable spatial resolution was designed, built, and integrated to a photon-counting picosecond streak camera, yielding a temporal scan resolution better than 300 femtosecond/pixel. The streak camera is configured to operate with 3 spatial channels; two of these support green (532 nm) and uv (355 nm) while the third accommodates reference pulses (764 nm) for real-time calibration. Critical parameters affecting differential timing accuracy such as pulse width and shape, number of received photons, streak camera/imaging system nonlinearities, dynamic range, and noise characteristics were investigated to optimize the system for accurate differential delay measurements. The streak camera output image consists of three image fields, each field is 1024 pixels along the time axis and 16 pixels across the spatial axis. Each of the image fields may be independently positioned across the spatial axis. Two of the image fields are used for the two wavelengths used in the experiment; the third window measures the temporal separation of a pair of diode laser pulses which verify the streak camera sweep speed for each data frame. The sum of the 16 pixel intensities across each of the 1024 temporal positions for the three data windows is used to extract the three waveforms. The waveform data is processed using an iterative three-point running average filter (10 to 30 iterations are used) to remove high-frequency structure. The pulse pair separations are determined using the half-max and centroid type analysis. Rigorous experimental verification has demonstrated that this simplified process provides the best measurement accuracy. To calibrate the receiver system sweep, two laser pulses with precisely known temporal separation are scanned along the full length of the sweep axis. The experimental measurements are then modeled using polynomial regression to obtain a best fit to the data. Data aggregation using normal point approach has provided accurate data fitting techniques and is found to be much more convenient than using the full rate single shot data. The systematic errors from this model have been found to be less than 3 ps for normal points.

Varghese, Thomas K.↗
DOE OSTI2021

SWAT model application for evaluating agricultural conservation practice effectiveness in reducing phosphorous loss from the Western Lake Erie Basin

Lake Erie is threatened by eutrophication and harmful algal blooms due to excess nutrient loading from agricultural sources. Agricultural conservation practices (ACPs) have been developed and implemented to reduce nutrient losses but estimating ACP effectiveness is challenging. The Soil and Water Assessment Tool (SWAT) has been used to investigate ACP effectiveness for water quality improvement. Many SWAT applications have been developed by different investigators to evaluate ACP effectiveness for reducing nutrient, particularly phosphorus (P), loading in the agriculturally-dominated Western Lake Erie Basin (WLEB). Our objective is to establish what has been achieved by past modeling research and make suggestions for future applications and improvements. We synthesized the findings of 28 SWAT modeling studies within the WLEB. Models generally performed satisfactorily against accepted criteria for streamflow and sediment, but performance for P loads, like soluble reactive P, was mostly “unsatisfactory”. The “unsatisfactory” performance maybe due to imperfections and idealizations in model formulations and/or parameterization. Thus, simulations of P transport and transformation processes need improvement. In addition, model parameter selection is the key part of model set-up. Most SWAT modeling studies used default values during initial set-up, then performed calibration and validation. It was found that the calibrated P related parameter values varied widely across different studies, even within the same watershed with some values unrealistic for the study areas. The phenomena of different combinations of model parameters producing similar outputs indicates equifinality. Equifinality in the baseline model may impact results when ACPs are incorporated. Furthermore, the unrealistic values used in ACP assessment undermine the credibility of ACP effectiveness. Future model applications should try to re-examine the calibrated P parameters and make sure they are realistic for the study area as well as reduce equifinality by constraining the model with characterization of watershed conditions, better understanding of hydrologic processes, and parameter values based on real-world observations. As a result, future model applications should focus on improving P transport and transformation processes, using measured watershed characteristics for parameterization, and improving reflections of climate change, which could result in more accurate assessments of ACP effectiveness to meet targeted goals.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Manifold learning for coarse-graining atomistic simulations: Application to amorphous solids

In this work, we introduce a generalized machine learning framework to probabilistically parameterize upper-scale models in the form of nonlinear PDEs consistent with a continuum theory, based on coarse-grained atomistic simulation data of mechanical deformation and flow processes. The proposed framework utilizes a hypothesized coarse-graining methodology with manifold learning and surrogate-based optimization techniques. Coarse-grained high-dimensional data describing quantities of interest of the multiscale models are projected onto a nonlinear manifold whose geometric and topological structure is exploited for measuring behavioral discrepancies in the form of manifold distances. A surrogate model is constructed using Gaussian process regression to identify a mapping between stochastic parameters and distances. Derivative-free optimization is employed to adaptively identify a unique set of parameters of the upper-scale model capable of rapidly reproducing the system's behavior while maintaining consistency with coarse-grained atomic-level simulations. The proposed method is applied to learn the parameters of the shear transformation zone (STZ) theory of plasticity that describes plastic deformation in amorphous solids as well as coarse-graining parameters needed to translate between atomistic and continuum representations. We show that the methodology is able to successfully link coarse-grained microscale simulations to macroscale observables and achieve a high-level of parity between the models across scales.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Alkali-silica reaction expansion model for confined concrete with stress-dependency and casting direction anisotropy

Alkali-silica reaction (ASR) is a deleterious chemical reaction between alkali hydroxyl ions and types of silica found in some aggregates of concrete. Owners and regulators of nuclear power plants aim to ensure the safety of the concrete structures with optimal maintenance strategies. A new model was developed for predicting the expansion of concrete structures affected by alkali-silica reaction. The model includes a novel combination of existing models as an alkali-silica reaction advancement model, a casting direction anisotropic expansion model, a stress-dependent anisotropic expansion model, and a material property evolution model dependent on the degree of ASR expansion. The model parameters were calibrated based on existing literature data and data generated by previous efforts of this study. The calibrated model was then validated with the experiments carried out in previous efforts of this study. The model was shown to accurately predict the ASR-expansion of large-scale reinforced concrete specimens with confinement.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Krypton-81 dating of the deep Continental Intercalaire aquifer with implications for chlorine-36 dating

Deep groundwater samples from the Continental Intercalaire (CI) aquifer in the Northern Tunisian Sahara have been analyzed for noble gases 3 He, 4 He, Ne and 81 Kr, and for 14 C to better constrain the groundwater residence time of this large transboundary aquifer. Its significant radiogenic 4 He content and background-level 14 C both indicate water older than a few tens of thousands of years. Distinct helium concentrations and Ne/He ratios suggest different groundwater flow paths through the Tozeur and Kebili regions. 81 Kr is applied for the first time on the CI aquifer, providing more direct evidence for the presence of old groundwater with 81 Kr/Kr ratios at 63% – 15% of the atmospheric value, corresponding to residence times of 150 – 630 ka. These newly obtained 81 Kr ages enable a re-evaluation of the previously reported data from the Tozeur region in order to better constrain the input values of 36 Cl/Cl and Cl. The initial 36 Cl/Cl ratio is found to be 4 – 6 times larger than previously assumed, and the initial Cl concentration 6 – 15 times smaller. The results of 81 Kr dating and the recalibrated 36 Cl dating define an eastward age progression across > 500 km of flow path from the Algerian to Tunisian sections. Here, this study exemplifies the utility of 81 Kr not only as an independent and conservative chronometer, but also as a supplemental measurement to better calibrate the parameters required for other age tracers, including cosmogenic 36 Cl and radiogenic 4 He.

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

Third-body stabilization of supercritical CO 2 in CO oxidation: development and application of a ReaxFF force field for the CO/O/CO 2 system

Supercritical CO 2 (scCO 2 ) plays a crucial role as a solvent in separation processes, advanced power cycles, and materials processing. Nonetheless, the atomistic comprehension of how the dense scCO 2 matrix influences the fundamental reaction of carbon monoxide (CO) is still insufficiently explored. Experimental studies and molecular dynamics (MD) simulations frequently fail to detect the highly reactive, transient intermediates, such as atomic oxygen (O), that drive these reactions. Here, to address this issue, we have developed a novel ReaxFF reactive force field for the CO 2 /CO/O system. The force field parameters were calibrated using density functional theory and second-order Møller-Plesset calculations to model CO 2 crystal properties, intermolecular interactions, bond dissociation curves, and reaction energy barriers. The force field reproduces the cohesive energy of the CO 2 crystal, the pressure characteristics of bulk scCO 2 , the equation-of-state behavior over a wide pressure–density range, the pressure dependence of the C–O bond length under compression, and the structural properties of liquid and scCO 2 , as documented by experiments, ab-initio MD, and prominent non-reactive models. The force field was subsequently applied to study the CO + O → CO 2 reaction. In a dilute environment, the reaction is inefficient as the newly formed CO 2 rapidly dissociates due to excess kinetic and potential energy acquired from the exothermic reaction. Conversely, in a dense scCO 2 environment, the surrounding matrix acts as an efficient third body, stabilizing the emerging CO2 product via molecular collisions. Statistical analysis confirms an average excess energy dissipation of 133.9 ± 3.6 kcal/mol over 112.4 ± 17.9 ps. Kinetic energy decomposition reveals that ∼ 92% of the excess kinetic energy is stored in internal (rotational and vibrational) degrees of freedom. This ReaxFF force field establishes a mechanistic foundation for third-body stabilization in dense reactive environments.

Chowdhury, Emdadul Haque [Pennsylvania State Univ.↗
DOE OSTI2020

Nonlinear sparse Bayesian learning for physics-based models

This paper addresses the issue of overfitting while calibrating unknown parameters of over-parameterized physics-based models with noisy and incomplete observations. Here, a semi-analytical Bayesian framework of nonlinear sparse Bayesian learning (NSBL) is proposed to identify sparsity among model parameters during Bayesian inversion. NSBL offers significant advantages over machine learning algorithm of sparse Bayesian learning (SBL) for physics-based models, such as 1) the likelihood function or the posterior parameter distribution is not required to be Gaussian, and 2) prior parameter knowledge is incorporated into sparse learning (i.e. not all parameters are treated as questionable). NSBL employs the concept of automatic relevance determination (ARD) to facilitate sparsity among questionable parameters through parameterized prior distributions. The analytical tractability of NSBL is enabled by employing Gaussian ARD priors and by building a Gaussian mixture-model approximation of the posterior parameter distribution that excludes the contribution of ARD priors. Subsequently, type-II maximum likelihood is executed using Newton's method whereby the evidence and its gradient and Hessian information are computed in a semi-analytical fashion. We show numerically and analytically that SBL is a special case of NSBL for linear regression models. Subsequently, a linear regression example involving multimodality in both parameter posterior pdf and model evidence is considered to demonstrate the performance of NSBL in cases where SBL is inapplicable. Next, NSBL is applied to identify sparsity among the damping coefficients of a mass-spring-damper model of a shear building frame. These numerical studies demonstrate the robustness and efficiency of NSBL in alleviating overfitting during Bayesian inversion of nonlinear physics-based models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Impact of lead on an axisymmetric, single bead blown powder DED overhung geometry

Metal additive manufacturing can be utilized for the near net-shape manufacture of components in a layer wise technique. Traditionally, the buildup process is conducted in the direction of gravity or bottom to top in the vertical orientation; however, with the availability of commercial systems with additive heads or fixturing that can index from vertical, and the need to manufacture larger parts without modifying available systems, the limitation of manufacturing from bottom to top is removed enabling the buildup of larger components printed at an angle from vertical. The effect of gravity on the melt pool and as-deposited component quality when printing off vertical is unknown in literature, especially for axisymmetric components. This understanding is critical for the advancement of manufacturing components of increased size in 4 and 5-axis additive systems. This investigation utilizes blown powder directed energy deposition to evaluate the change in as-printed geometry when the start point is altered in relation to gravity while the part rotates to manufacture an axisymmetric component. The objectives of this work are to determine the impact of the deposition location on the geometric variability on axisymmetric components. This investigation tests the hypothesis that differences in layer height exist due to a change in catchment efficiency when the deposition location is moved to a tangent of the round geometry due to a change in the melt pool dynamics. It was found the ideal deposition location when printing at 26.6-degrees from vertical was at −90-degrees from the top center point of the component, along the tangent, where gravity was pushing the melt pool down, but the rotation of the part was pulling the deposited material towards the top center of the component. This work provides an understanding of layer height stability and catchment efficiency to guide print orientation strategy for high-aspect ratio components. It was found the −90-degree lead deposition had the best layer height stability at 0.6 % as compared to the top-center at 6.3 % and +90-degree deposition location at 9.2 % for the nominal programmed layer height. The change in layer height also effected the final diameter the greatest for the top center deposition location where the diameter diverged by 2.5 % during the overbuild condition and converged by 0.6 %. In conclusion, this finding will increase the manufacturing efficiency of axisymmetric components by increasing the passive stability of the printing process for the successful manufacture of parts without the need for perfectly calibrated manufacturing parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Stress localization investigation of additively manufactured GRCop-42 thin-wall structure

A full-field crystal plasticity (CP) framework is presented for the GRCop-42 alloy to study microscopic mechanical behavior and local stress heterogeneities. The microstructures of additively manufactured (AM) materials are often unique relative to conventionally processed materials, and the local thermal histories drive these differences during the build process. These thermal histories depend on the process parameters (laser power, scan speed, and scan strategy) and the part geometry. Prior research has shown that the mechanical properties of thin-walled structures can vary significantly with wall thickness due to changes in the thermal boundary conditions during manufacturing. It is, therefore, desirable to perform CP simulations based on the phenomenological constitutive model to predict the local mechanical responses induced by microstructural heterogeneities. This work generates representative microstructures based on experimentally collected grain information (i.e., texture) for grain scale stress analysis, and the material constitutive parameters are calibrated using the experimental mechanical testing data. Here, we specifically investigated the effect of crystallographic texture and grain morphologies on the size-dependent mechanical properties of AM GRCop-42. The selection of appropriate material properties for implementing an effective free surface boundary condition and the influence of adjacent buffer layers are also discussed. Analysis of local field results reveals a strong correlation between stress localization and the initial grain orientation. However, no significant relationship between the misorientation of the individual adjacent grains and the average misorientation is observed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Improving a Multilevel Turbulence Closure Model for a Shallow Lake in Comparison With Other 1-D Models

Lakes differ from lands in water availability, heat capacity, albedo, and roughness, which affect local surface-atmospheric interactions. This study modified a multilevel upper ocean model (UOM) for lake applications and evaluated its performance in Lake Taihu (China) with comprehensive measurements against three popular one-dimensional (1-D) lake models. These models were based on different concepts, including the self-similarity (FLake), the wind-driven eddy diffusion (LISSS), the k-ε turbulence closure (SIMSTRAT), and a simplified turbulence closure (UOM). The surface flux scheme in these models was unified to exclude the discrepancies in representing air-lake exchanges. All models in their default formulations presented obvious cold water temperature biases and largely underestimated the lake surface temperature (LST) diurnal range. For each model, these deficiencies were significantly reduced by incorporating new physics schemes or calibrated tunable parameters based on systematic sensitivity tests. The primary modifications for UOM included (1) a new scheme of decreased surface roughness lengths to better characterize the shallow lake, (2) a solar radiation penetration scheme with increased light extinction coefficient and surface absorption fraction to account for the high water turbidity, and (3) turbulent Prandtl number increased by a factor of 20 to reduce the turbulent vertical mixing. All other models were improved in these three aspects (roughness, extinction, and mixing) within their original formulations. Given these improvements, UOM showed superior performance to other models in capturing LST diurnal cycle and daily to seasonal variations, as well as summer-autumn vertical stratification changes. The new UOM is well suited for application in shallow lakes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Modeling root system growth around obstacles

State-of-the-Art models of Root System Architecture (RSA) do not allow simulating root growth around rigid obstacles. Yet, the presence of obstacles can be highly disruptive to the root system. We grew wheat seedlings in sealed petri dishes without obstacle and in custom 3D-printed rhizoboxes containing obstacles. Time-lapse photography was used to reconstruct the wheat root morphology network. We used the reconstructed wheat root network without obstacle to calibrate an RSA model implemented in the R-SWMS software. The root network with obstacles allowed calibrating the parameters of a new function that models the influence of rigid obstacles on wheat root growth. Experimental results show that the presence of a rigid obstacle does not affect the growth rate of the wheat root axes, but that it does influence the root trajectory after the main axis has passed the obstacle. The growth recovery time, i.e. the time for the main root axis to recover its geotropism-driven growth, is proportional to the time during which the main axis grows along the obstacle. Qualitative and quantitative comparisons between experimental and numerical results show that the proposed model successfully simulates wheat RSA growth around obstacles. Our results suggest that wheat roots follow patterns that could inspire the design of adaptive engineering flow networks.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Single Crystal Plasticity Model with Deformation Twinning for the High Rate Deformation of β-HMX

Deformation twinning is an active deformation mechanism during loading along certain crystallographic directions of β-cyclotetramethylene tetranitramine (β-HMX). In this work, a finite strain thermomechanical model is extended to include twinning as a deformation mechanism in addition to plastic slip. The stress is derived from the free energy expression including a term representing the equation of state. The crystal plasticity framework is used to divide the total strain into inelastic and elastic component, where the elastic part is used in the expression for the free energy. The twin systems are treated as pseudo slip systems and the shear rate on the twin systems is evaluated in terms of the twin resistance and the projection of the stress tensor. The model parameters were calibrated against a set of plate impact experiments performed on β-HMX by Dick et al. The remaining plate impact experiments are used to evaluate the predictive capability of the model. The quality of the model fits and predictions is discussed from a physical and modeling perspective. Particularly, the role of twin modeling on results is highlighted. Finally, the model is used to explore the relationship between the propensity for twinning and crystal orientation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Effect of surface defects on damage and failure of dynamically driven cylindrical shells

Here, in this work, expanding cylinder experiments are performed using a gas gun on the aluminum 6061 alloy to understand the dynamic deformation and failure behavior of this material in the presence of surface defects in the form of scratches. The experiments are then modeled with the modified Tepla (for tensile plasticity) model, which is equipped with the damage yield function of Gurson [“Continuum theory of ductile rupture by void nucleation and growth: Part I—Yield criteria and flow rules for porous ductile media,” (1977)] and Tvergaard and Needleman [Acta Metall. 32(1), 157–169 (1984)], over-stress viscosity [F. Addessio and J. Johnson, J. Appl. Phys. 74(3), 1640–1648 (1993)], and the underlying viscoplastic model of Preston et al. [J. Appl. Phys. 93(1), 211–220 (2003)] (Preston–Tonks–Wallace model) for the dense material. The Tepla model parameters were calibrated using the Bayesian approach against experimentally measured free-surface velocity data obtained from plate impact experiments on the aluminum 6061 alloy. The evolution of radial expansion and shapes of the deformed cylinder from the simulations on a pristine cylinder wall and a cylinder wall with a longitudinal surface defect shows excellent agreement with the corresponding experimental results. Finally, the role of surface defect size and cylinder wall geometry on the dynamic strength of the cylinder wall was also investigated through a series of simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Temperature-based reactive flow model for triaminotrinitrobenzene (TATB) plastic bonded explosives

A new reactive flow model is presented for triaminotrinitrobenzene (TATB)-based plastic bonded explosives, applicable to shock initiation and steady detonation problems of differing initial temperature. Temperature disequilibrium is assumed between unreacted explosive, material in the vicinity of compressed defects (called hot spots), and reaction products. The model incorporates temperature-dependent decomposition reaction rates. Particularly, Arrhenius model parameters were derived from quantum-based molecular dynamics simulations of TATB decomposition. Further, a model of detonation carbon aggregation is incorporated, describing the slow release of energy inherent to detonation in TATB-based materials. Model parameters were calibrated against gas gun shock initiation experiments and steady detonation rate stick tests. The predictive ability of the model in the shock initiation regime is tested against recent thin pulse experiments. The model is found to perform equally well in predicting the size-effect curve of ambient, cold, and hot rate sticks. The present work demonstrates the viability of incorporating results from subscale simulations into a continuum-scale reactive flow model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗