Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sensitivity Simulation”

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

An experimental and kinetic modeling study of cyclopentane and dimethyl ether blends

Cyclopentane is a suitable naphthene, or cycloalkane, in a palette for multi-component gasoline surrogate fuels due to its presence in market fuels and its relevance to alkyl substituted cyclopentanes also present. However, the previous oxidation studies of cyclopentane have primarily focused on neat mixtures. Blending cyclopentane with dimethyl ether in this work therefore serves to inform our understanding of, and improve predictive models for, multi-component mixtures. Here, the auto-ignition of cyclopentane/dimethyl ether blends was studied in a high-pressure shock tube and in a rapid compression machine. A wide range of temperatures (650 – 1350 K) and elevated pressures of 20 and 40 bar were studied at equivalence ratios of 0.5, 1.0 and 2.0 in air for two blending ratios (30/70 and 70/30 mole% cyclopentane/di-methyl ether mixtures). A detailed kinetic model for cyclopentane was revised to capture the measured ignition delay times and apparent heat release rates in this study. Literature ignition delay time, jet-stirred reactor, and laminar burning velocity measurements of neat cyclopentane were used as additional validation. Improvements to the kinetic model were based on recent literature studies related to sub-models including cyclopentene and cyclopentadiene which allowed the removal of previous local rate-constant optimizations. Low temperature reactivity of cyclopentane was found to be controlled by the branching ratio between concerted elimination of HO͘ 2 and the strained formation of Q&#775OOH radicals in agreement with previous studies. In this study, the low branching ratio of Q&#775OOH formation increases the influence of a competing consumption pathway for cyclopentyl-peroxy (CPTO͘ 2 J) radicals. The sensitivity of the simulated ignition delay times to the formation of cyclopentyl hydroperoxide (CPTO 2 H), from CPTO͘ 2 J and HO͘ 2 , is discussed. The current model is used to analyze the influence of dimethyl ether on the reactivity of cyclopentane in the context of previous literature studies of dimethyl ether binary blends with ethanol and toluene.

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

Machine-learning techniques for model-independent searches in dijet final states

Anomaly detection methods used in a recent search for new phenomena by CMS at the CERN LHC are presented. The methods use machine learning to detect anomalous jets produced in the decay of new massive particles without depending on a specific theory model. The effectiveness of these approaches in enhancing sensitivity to various simulated signal samples is studied and compared using data collected in proton–proton collisions at a center-of-mass energy of 13 TeV. In an example analysis, the capabilities of anomaly detection methods are further demonstrated by identifying large-radius jets consistent with Lorentz-boosted hadronically decaying top quarks in a model-agnostic framework.

CMS↗
DOE OSTI2020

An Assessment of the Drag Models in the Case of a Shock Interacting With a Fixed Bed of Point Particles

Abstract In this work, three-dimensional Euler–Lagrange (EL) point-particle simulations of a shock wave interacting with a fixed bed of particles are carried out. The results from the particle-resolved (PR) simulations are used to assess the performance of the point-particle drag models during short time scales. We demonstrate that in a one-way coupled regime, the point-particle simulations recover the dominant gas dynamic features of the flow and are in a good agreement with the exact Riemann solution of a shock traveling through a sudden area contraction. Although the PR simulations are inviscid, we show that a dissipative drag is necessary to predict the mean behavior of the gas. As a model for the inviscid shock-induced (SI) drag two different models are presented in lieu of the quasi-steady drag. Finally, two-way coupled simulations are performed at four different particle volume fractions {0.10, 0.15, 0.20, 0.25} and three different incident shock Mach numbers {1.22, 1.66, 3.0} and compared against the data from PR inviscid simulations. At a lower Mach number (1.22), averaged flow quantities from the two-way coupled simulations agree well with the PR simulations. As the Mach number increases, we observe that the discrepancies between the point-particle and the PR simulations grow. A sensitivity analysis of the drag models involved reveals a strong influence of the inviscid-unsteady force on the gas quantities especially in the case of a strong shock interacting with a dense bed of particles. The use of Mach correlation beyond the subcritical regime coupled with the model for volume fraction correction is identified as a probable cause for the additional drag.

Engineering↗
DOE OSTI2023

Changing temperature profiles and the risk of dengue outbreaks

As temperatures change worldwide, the pattern and competency of disease vectors will change, altering the global distribution of both the burden of infectious disease and the risk of the emergence of those diseases into new regions. To evaluate the risk of potential summer dengue outbreaks triggered by infected travelers under various climate scenarios, we develop an SEIR-type model, run numerical simulations, and conduct sensitivity analyses under a range of temperature profiles. Our model extends existing theoretical frameworks for studying dengue dynamics by introducing temperature dependence of two key parameters: the mosquito extrinsic incubation period and the lifespan of mosquitoes, which empirical data suggests are both highly temperature dependent. We find that changing temperature significantly alters dengue risk in an inverted U-shape, with temperatures in the range 27-31°C producing the highest risk. As temperatures increase beyond 31°C, the determinants of dengue risk begin to shift from mosquito biting rate and carrying capacity to the duration of the human infectious period, suggesting that changing temperatures not only alter dengue risk but also the potential efficacy of control measures. To illustrate the role of spatial and temporal temperature heterogeneity, we select five US cities where the primary dengue vector, the mosquito Aedes aegypti , has been observed, and which have had dengue cases in the past: Los Angeles, Houston, Miami, Brownsville, and Phoenix. Our analysis suggests that an increase of 3°C leads to an approximate doubling of the risk of dengue in Los Angeles and Houston, but a reduction of risk in Miami, Brownsville, and Phoenix due to extreme heat.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

On the Response of a Herschel–Bulkley Fluid Due to a Moving Plate

In this paper, we study the boundary-layer flow of a Herschel–Bulkley fluid due to a moving plate; this problem has been experimentally investigated by others, where the fluid was assumed to be Carbopol, which has similar properties to cement. The computational fluid dynamics finite volume method from the open-source toolbox/library OpenFOAM is used on structured quad grids to solve the mass and the linear momentum conservation equations using the solver “overInterDyMFoam” customized with non-Newtonian viscosity libraries. The governing equations are solved numerically by using regularization methods in the context of the overset meshing technique. The results indicate that there is a good comparison between the experimental data and the simulations. The boundary layer thicknesses are predicted within the uncertainties of the measurements. The simulations indicate strong sensitivities to the rheological properties of the fluid.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Expected Gamma-Ray Burst Detection Rates and Redshift Distributions for the BlackCAT CubeSat Mission

We report the results of an extensive set of simulations exploring the sensitivity of the BlackCAT CubeSat to long-duration gamma-ray bursts (GRBs). BlackCAT is a NASA APRA-funded CubeSat mission for the detection and real-time subarcminute localization of high-redshift (z ≳ 3.5) GRBs. Thanks to their luminous and long-lived afterglow emissions, GRBs are uniquely valuable probes of high-redshift star-forming galaxies and the intergalactic medium. In addition, each detected GRB with a known redshift serves to localize a region of high-redshift star formation in three dimensions, enabling deep follow-on searches for host galaxies and associated local and large-scale structures. We explore two distinct models for the GRB redshift distribution and luminosity function, both consistent with Swift observations. We find that, for either model, BlackCAT is expected to detect a mean of 42 bursts per year on orbit, with 6.7% to 10% of these at z > 3.5. BlackCAT bursts will be localized to an r 90 ≲ 55'' precision and reported to the community within seconds. Due to the mission orbit and pointing scheme, bursts will be located in the night sky and well placed for deep multiwavelength follow-up observations. BlackCAT is on schedule to achieve launch readiness in 2025.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Influences of an entrainment–mixing parameterization on numerical simulations of cumulus and stratocumulus clouds

Abstract. Different entrainment–mixing processes can occur in clouds; however, a homogeneous mixing mechanism is often implicitly assumed in most commonly used microphysics schemes. Here, we first present a new entrainment–mixing parameterization that uses the grid mean relative humidity without requiring the relative humidity of the entrained air. Then, the parameterization is implemented in a microphysics scheme in a large eddy simulation model, and sensitivity experiments are conducted to compare the new parameterization with the default homogeneous entrainment–mixing parameterization. The results indicate that the new entrainment–mixing parameterization has a larger impact on the number concentration, volume mean radius, and cloud optical depth in the stratocumulus case than in the cumulus case. This is because inhomogeneous and homogeneous mixing mechanisms dominate in the stratocumulus and cumulus cases, respectively, which is mainly due to the larger turbulence dissipation rate in the cumulus case. Because stratocumulus clouds break up during the dissipation stage to form cumulus clouds, the effects of this new entrainment–mixing parameterization during the stratocumulus dissipation stage are between those during the stratocumulus mature stage and the cumulus case. A large aerosol concentration can enhance the effects of this new entrainment–mixing parameterization by decreasing the cloud droplet size and evaporation timescale. The results of this new entrainment–mixing parameterization with grid mean relative humidity are validated by the use of a different entrainment–mixing parameterization that uses parameterized entrained air properties. This study sheds new light on the improvement of entrainment–mixing parameterizations in models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Distinct drivers of recent seasonal precipitation increase over Central Asia: roles of anthropogenic aerosols and greenhouse gases

Observational evidence reveals a pronounced wetting trend over Central Asia in recent decades, with the most substantial increases occurring during winter and summer. Yet the extent to which the drivers of these changes differ seasonally remains unknown. Here, we use single-forcing experiments from the Precipitation Driver and Response Model Intercomparison Project (PDRMIP) to examine the effects of various external forcings on winter and summer precipitation across Central Asia and to explore the physical mechanisms underlying seasonal precipitation changes. We find that greenhouse gas (GHG) forcing mainly increases winter precipitation by enhancing atmospheric moisture content through warming. In contrast, in summer, Asian sulfate aerosols enhance precipitation by modulating the westerly jet, which strengthens atmospheric moisture transport into the region. Asian black carbon exerts an opposing influence that partially offsets the sulfate-induced effect. Further attribution analysis based on CMIP6 simulations reinforces these sensitivity results and shows that GHG forcing is the primary driver of winter precipitation increases whereas anthropogenic aerosols dominate summer trends. Future CMIP6 projections suggest that under moderate- to high-emission scenarios, winter precipitation will continue to rise due to increasing GHG concentrations, while summer precipitation may decline across much of Central Asia as a result of reduced aerosol emissions following Asian clean air policies. These findings highlight a distinct seasonality in the drivers of recent precipitation increase and suggest a plausible divergence in future winter and summer precipitation trends.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Numerical coupling of aerosol emissions, dry removal, and turbulent mixing in the E3SM Atmosphere Model version 1 (EAMv1) – Part 2: A semi-discrete error analysis framework for assessing coupling schemes

Abstract. Part 1 (Wan et al., 2024) of this study discusses the motivation and empirical evaluation of a revision to the aerosol-related numerical process coupling in the atmosphere component of the Energy Exascale Earth System Model version 1 (EAMv1) to address the previously reported issue of strong sensitivity of the simulated dust aerosol lifetime and dry removal rate to the model's vertical resolution. This paper complements that empirical justification of the revised scheme with a mathematical justification leveraging a semi-discrete analysis framework for assessing the splitting error of process coupling methods. The framework distinguishes the error due to numerical splitting from the error due to the time integration method(s) used for each individual process. Such a distinction results in a framework that provides an intuitive understanding of the causes of the splitting error. The application of this framework to the dust life cycle in EAMv1 confirms (i) that the original EAMv1 scheme artificially strengthens the effect of dry removal processes and (ii) that the revised splitting reduces that artificial strengthening. While the error analysis framework is presented in the context of the dust life cycle in EAMv1, the framework can be broadly leveraged to evaluate process coupling schemes, both in other physical problems and for any number of processes. This framework will be particularly powerful when the various process implementations support a variety of time integration approaches. Whereas traditional local truncation error approaches require separate consideration of each combination of time integration methods, this framework enables evaluation of coupling schemes independent of particular time integration approaches for each process while still allowing for the incorporation of these specific time integration errors if so desired. The framework also explains how the splitting error terms result from (i) the integration of individual processes in isolation from other processes and (ii) the choices of input state and time step size for the isolated integration of processes. Such a perspective has the potential for the rapid development of alternative coupling approaches that utilize knowledge both about the desired accuracy and about the computational costs of individual processes.

58 GEOSCIENCES↗
DOE OSTI2020

Data Requirements for Application of Risk-Based Dynamic Contingency Analysis to Evaluate Hurricane Impact to Electrical Infrastructure in Puerto Rico

This paper presents a risk-based dynamic contingency analysis framework that was used to evaluate the hurricane impact to electrical infrastructure in Puerto Rico. PNNL developed a scalable risk-based framework for identifying high-voltage transmission resilience improvements by classifying and prioritizing high-risk power grid contingencies (system failures) under hurricane impact. The risk-based framework is founded on grid outage definitions with their associated probabilities of occurrence from hurricane events, in combination with an impact assessment derived from detailed dynamic cascading failure analysis. This paper focuses on a discussion around data requirements for transmission resilience planning for hurricane events, derived from the development of the risk-based framework and its application to Puerto Rico. This paper launches an important first step in encouraging the engineering community and power system industry to move towards establishing resilience planning as a routine practice. Since actual results for Puerto Rico contain sensitive information, sample simulation results will be used to illustrate the data requirements and risk-based dynamic cascading framework on the Puerto Rico power grid, as well as demonstrate the potential for such a simulation framework. The paper includes a discussion on the lessons learned, importance and need for improved datasets that are not usually considered in traditional power system planning. The paper will also elaborate on how the scalable simulation framework and datasets might be expanded to larger footprints and leveraged for modelling other types of natural disasters.

DCAT, Puerto Rico, hurricane, Power System Stabili↗
DOE OSTI2022

Predictive Contaminant Transport Simulation with the P2R Model for the Composite Analysis Inventory Sensitivity Case

The Plateau to River (P2R) Model is a groundwater flow and contaminant fate and transport (F&T) simulation model used to support remedial activities conducted by the Central Plateau Cleanup Company (CPCCo) at the Hanford Site in Washington State. The model simulates contaminants of concern within the saturated zone of the uppermost aquifer beneath the Central Plateau and downgradient to the Columbia River. The overall objective of the saturated zone modeling effort is to provide a basis for making informed remedial action decisions based on descriptions of current and expected future contaminant concentrations in groundwater at decision points within and downgradient of the Central Plateau of the Hanford Site. Specifically, the purpose of this environmental calculation file (ECF) is to describe an inventory sensitivity case of the CA base case. The inventory sensitivity case implements a change in the activity contribution from the vadose zone in the A Trenches Area model (for tritium (H-3) and iodine-129 (I-129)), BC Cribs and Trenches model (for I-129 and technetium-99 (Tc-99)), and the PUREX Area model (for H-3 and I-129). All other simulated inventories are identical to the CA base case. The simulation of F&T of contaminants reported in this case will support dose predictions as part of the Hanford Site CA.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Back-to-back relative-excess observable to identify the chiral magnetic effect

Background: The chiral magnetic effect (CME) is extensively studied in heavy-ion collisions at RHIC and LHC. In the commonly used reaction-plane (RP)-dependent, charge-dependent azimuthal correlator ($Δ_γ$), both the close and back-to-back pairs are included. Many backgrounds contribute to the close pairs (e.g., resonance decays, jet correlations), whereas the back-to-back pairs are relatively free of those backgrounds. Purpose: In order to reduce those backgrounds, we propose a new observable, which only focuses on the back-to-back pairs, namely, the relative back-to-back opposite-sign over same-sign pair excess ($r_{BB}$) as a function of the pair azimuthal orientation with respect to the RP ($φ_{BB}$). Methods: We use analytical calculations and toy model simulations to demonstrate the sensitivity of $r_{BB}(φ_{BB})$ to the CME and its insensitivity to backgrounds. Results: With finite CME, the $φ_{BB}$ distribution of $r_{BB}$ shows a clear characteristic modulation. Its sensitivity to background is significantly reduced compared to the previous $Δ_γ$ observable. The simulation results are consistent with our analytical calculations. Conclusions: In this work, our studies demonstrate that the $r_{BB}(φ_{BB})$ observable is sensitive to the CME signal and rather insensitive to the resonance backgrounds.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Predictive Contaminant Transport Simulation with P2R Model for the Composite Analysis Limited Source Sensitivity Case.

The Plateau to River (P2R) is a groundwater flow and contaminant fate and transport (F&T) simulation model used to support remedial activities conducted by the Central Plateau Cleanup Company at the Hanford Site in Washington State. Figure 1 illustrates the P2R Model extents, discretization, and boundary conditions. The P2R Model is utilized in the Composite Analysis (CA) for the Hanford Site as the computational engine for computing F&T predictions as described in CP-60406, Hanford Site Composite Analysis Technical Approach Description: Groundwater. The model simulates contaminants of concern within the saturated zone of the uppermost aquifer beneath the Central Plateau and downgradient to the Columbia River. CP-57037, Model Package Report for the Plateau to River Model Version 8.3 documents the current version of the P2R Model including a description of the conceptual site model, model development and calibration, and limitations to the model application. Simulations conducted to support the dose calculations required by the CA are documented in ECF-HANFORD-19-0119, Predictive Flow Simulation with the P2R Model for the Composite Analysis Base Case and ECF-HANFORD-19-0120, Contaminant Transport Simulation with the P2R Model for the Composite Analysis Base Case.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Predictive Contaminant Transport Simulation with the P2R Model for the Composite Analysis Recharge Sensitivity Case

The Plateau-to-River (P2R) model is a groundwater flow and contaminant fate and transport (F&T) simulation model used to support remedial activities conducted by the Central Plateau Cleanup Company (CPCCo) at the Hanford Site in south-central Washington State. The overall objective of the saturated zone modeling effort is to provide a basis for making informed remedial action decisions based on descriptions of current and expected future contaminant concentrations in groundwater at decision points within and downgradient of the Central Plateau of the Hanford Site. Specifically, the purpose of this environmental calculation file (ECF) is to describe a recharge sensitivity case of the CA base case. The recharge sensitivity case implements a change in the activity contribution from the vadose zone in the A Trenches Area model (for tritium [H-3] and iodine-129 [I-129]), BC Cribs and Trenches model (for I-129 and technetium-99 [Tc-99]), or the Plutonium Uranium Extraction (PUREX) Area model (for H-3 and I-129). All other simulated inventories are identical to the CA base case. The simulation of fate and transport of contaminants reported in this case will support dose predictions as part of the updated Hanford Site CA.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Thiocyanate Ions Form Antiparallel Populations at the Concentrated Electrolyte/Charged Surfactant Interface

Anions play significant roles in the separation of lanthanides and actinides. The molecular-scale details of how these anions behave at aqueous interfaces are not well understood, especially at high ionic strengths. Here, we describe the interfacial structure of thiocyanate anions at a soft charged interface up to 5 M bulk concentration with combined classical and phase-sensitive and molecular dynamics (MD) simulations. At low concentrations thiocyanate ions are mostly oriented with their sulfur end pointing toward the charged surfactants. In this work, the VSFG signal reaches a plateau at around 100 mM bulk concentration, followed by significant changes above 1 M. At high concentrations a new thiocyanate population emerges with their sulfur end pointing toward the bulk liquid. The –CN stretch frequency is different for up and down oriented SCN – ions, indicating different coordination environments. These results provide key molecular-level insights for the interfacial behavior of complex anions in highly concentrated solutions.

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

In Situ Velocity‐Strain Sensitivity Near the San Jacinto Fault Zone Analyzed Through Train Tremors

Abstract We utilize train tremors as P‐wave seismic sources to investigate velocity‐strain sensitivity near the San Jacinto Fault Zone. A dense nodal array deployed at the Piñon Flat Observatory is used to detect and identify repeating train energy emitted from a railway in the Coachella valley. We construct P‐wave correlation functions across the fault zone and estimate the spatially averaged dt / t versus strain sensitivity to be 6.25 × 10 4 . Through numerical simulations, we explore how the sensitivity decays exponentially with depth. The optimal solution reveals a subsurface sensitivity of 1.2 × 10 5 and a depth decay rate of 0.05 km −1 . This sensitivity aligns with previous findings but is toward the higher end, likely due to the fractured fault‐zone rocks. The depth decay rate, previously unreported, is notably smaller than assumed in empirical models. This raises the necessity of further investigations of this parameter, which is crucial to study stress and velocity variations at seismogenic depth.

Geology↗
DOE OSTI2026

A Framework for Parametric and Predictive Uncertainty Quantification in the E3SM Land Model: Assessing Site and Observable Generalizability

Quantifying parametric uncertainty using observations from individual sites provides a critical foundation for Earth system modeling, serving as a necessary first step before scaling up to regional or global applications. This study introduces a novel computational framework designed to enhance model predictability by reducing parametric uncertainty and assessing site and observable generalizability using various observational constraints. The framework integrates five components: Model Simulation, Statistical Emulation, Global Sensitivity Analysis (GSA), Model Calibration, and Model Prediction. Using the E3SM land model, we simulated site-level land-atmosphere carbon and energy fluxes from 2003 to 2007 across five evergreen needleleaf FLUXNET sites, perturbing 26 vegetation-related model parameters. Gaussian process emulators were employed to expedite GSA and model calibration. Four critical parameters that strongly influence selected land-atmosphere fluxes were identified by GSA. Bayesian approaches were used to infer parameter probability distributions leveraging synthetic data and FLUXNET observations. The results reveal that posterior parameter distributions vary significantly across different sites and observables within the same plant functional type. Probabilistic predictions indicate that parameters calibrated at one site can enhance predictive accuracy at other sites, although site heterogeneity may sometimes outweigh parametric uncertainty. Additionally, the probabilistic predictions demonstrate that calibration for one variable can also improve predictability for other variables, thereby maximizing predictive capabilities with limited observations. This framework provides a powerful approach for reducing parametric uncertainty in Earth system models and deepening our understanding of carbon dynamics and energy cycles. Its adaptability makes it a valuable tool for broader applications in Earth system modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Mesh objective stochastic simulations of quasibrittle fracture

Continuum finite element (FE) modeling of damage and failure of quasibrittle structures suffers from the spurious mesh sensitivity due to strain localization. Here this issue has been addressed for deterministic analysis through the development of localization limiters. Here this study proposes a mechanism-based model to mitigate the mesh sensitivity in stochastic FE simulations of quasibrittle fracture. The interest is placed on the analysis of large-size structures, where the mesh size is conveniently chosen to be larger than the width of the fracture process zone as well as the correlation length of the random fields of constitutive properties. The present model is formulated within the framework of continuum damage mechanics. Two localization parameters are introduced to describe the evolution of the damage pattern of each finite element. These parameters are used to guide the energy regularization of the constitutive law, as well as to formulate the mesh-dependent probability distributions of constitutive properties. Depending on the prevailing damage pattern, different energy regularization schemes and mesh dependence of the probability distribution functions are used in the constitutive law. The model is applied to simulate the stochastic failure behavior of quasibrittle structures of different geometries featuring different failure processes including damage initiation, localization, and propagation. It is shown that using fixed probability distribution functions of constitutive properties could lead to strong mesh dependence of the prediction of the mean and variance of the structural load capacity. The probability distribution functions of constitutive properties must be linked to the damage pattern, which may evolve during the failure process. Such a mechanism-based modeling of the probability distributions of constitutive properties is essential for mitigating the spurious mesh sensitivity in stochastic FE analysis of quasibrittle fracture.

42 ENGINEERING↗