Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reproducibility”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Forced and Unforced Variability of Twentieth Century North American Droughts and Pluvials

Research on the forcing of drought and pluvial events over North America is dominated by general circulation model experiments that often have operational limitations (e.g., computational expense, ability to simulate relevant processes, etc). We use a statistically based modeling approach to investigate sea surface temperature (SST) forcing of the twentieth century pluvial (1905-1917) and drought (1932-1939, 1948-1957, 1998-2002) events. A principal component (PC) analysis of Palmer Drought Severity Index (PDSI) from the North American Drought Atlas separates the drought variability into five leading modes accounting for 62% of the underlying variance. Over the full period spanning these events (1900-2005), the first three PCs significantly correlate with SSTs in the equatorial Pacific (PC 1), North Pacific (PC 2), and North Atlantic (PC 3), with spatial patterns (as defined by the empirical orthogonal functions) consistent with our understanding of North American drought responses to SST forcing. We use a large ensemble statistical modeling approach to determine how successfully we can reproduce these drought/pluvial events using these three modes of variability. Using Pacific forcing only (PCs 1-2), we are able to reproduce the 1948-1957 drought and 1905-1917 pluvial above a 95% random noise threshold in over 90% of the ensemble members; the addition of Atlantic forcing (PCs 1-2-3) provides only marginal improvement. For the 1998-2002 drought, Pacific forcing reproduces the drought above noise in over 65% of the ensemble members, with the addition of Atlantic forcing increasing the number passing to over 80%. The severity of the drought, however, is underestimated in the ensemble median, suggesting this drought intensity can only be achieved through internal variability or other processes. Pacific only forcing does a poor job of reproducing the 1932-1939 drought pattern in the ensemble median, and less than one third of ensemble members exceed the noise threshold (28%). Inclusion of Atlantic forcing improves the ensemble median drought pattern and nearly doubles the number of ensemble members passing the noise threshold (52%). Even with the inclusion of Atlantic forcing, the intensity of the simulated 1932-1939 drought is muted, and the drought itself extends too far into the southwest and southern Great Plains. To an even greater extent than the 1998-2002 drought, these results suggest much of the variance in the 1932-1939 drought is dependent on processes other than SST forcing. This study highlights the importance of internal noise and non SST processes for hydroclimatic variability over North America, complementing existing research using general circulation models.

Cook, Benjamin I.↗
NASA NTRS2013

Supersonic Retropulsion Surface Preparation of Carbon Fiber Reinforced Epoxy Composites for Adhesive Bonding

Surface preparation is widely recognized as a key step to producing robust and predictable bonds in a precise and reproducible manner. Standard surface preparation techniques, including grit blasting, manual abrasion, and peel ply, can lack precision and reproducibility, which can lead to variation in surface properties and subsequent bonding performance. The use of a laser to ablate composite surface resin can provide an efficient, precise, and reproducible means of preparing composite surfaces for adhesive bonding. Advantages include elimination of physical waste (i.e., grit media and sacrificial peel ply layers that ultimately require disposal), reduction in process variability due to increased precision (e.g. increased reproducibility), and automation of surface preparation, all of which improve reliability and process control. This paper describes a Nd:YAG laser surface preparation technique for composite substrates and the mechanical performance and failure modes of bonded laminates thus prepared. Additionally, bonded specimens were aged in a hot, wet environment for approximately one year and subsequently mechanically tested. The results of a one year hygrothermal aging study will be presented.

Palmieri, Frank L.↗
NASA NTRS2013

Stellar Laboratories: New GeV and Ge VI Oscillator Strengths and their Validation in the Hot White Dwarf RE0503-289

State-of-the-art spectral analysis of hot stars by means of non-LTE model-atmosphere techniques has arrived at a high level of sophistication. The analysis of high-resolution and high-S/N spectra, however, is strongly restricted by the lack of reliable atomic data for highly ionized species from intermediate-mass metals to trans-iron elements. Especially data for the latter has only been sparsely calculated. Many of their lines are identified in spectra of extremely hot, hydrogen-deficient post-AGB stars. A reliable determination of their abundances establishes crucial constraints for AGB nucleosynthesis simulations and, thus, for stellar evolutionary theory. Aims. In a previous analysis of the UV spectrum of RE 0503-289, spectral lines of highly ionized Ga, Ge, As, Se, Kr, Mo, Sn, Te, I, and Xe were identified. Individual abundance determinations are hampered by the lack of reliable oscillator strengths. Most of these identified lines stem from Ge V. In addition, we identified Ge VI lines for the first time. We calculated Ge V and Ge VI oscillator strengths in order to reproduce the observed spectrum. Methods. We newly calculated Ge V and Ge VI oscillator strengths to consider their radiative and collisional bound-bound transitions in detail in our non-LTE stellar-atmosphere models for the analysis of the Ge IV-VI spectrum exhibited in high-resolution and high-S/N FUV (FUSE) and UV (ORFEUS/BEFS, IUE) observations of RE 0503-289. Results. In the UV spectrum of RE 0503-289, we identify four Ge IV, 37 Ge V, and seven Ge VI lines. Most of these lines are identified for the first time in any star. We can reproduce almost all Ge IV, GeV, and Ge VI lines in the observed spectrum of RE 0503-289 (T(sub eff) = 70 kK, log g = 7.5) at log Ge = -3.8 +/- 0.3 (mass fraction, about 650 times solar). The Ge IV/V/VI ionization equilibrium, that is a very sensitive T(sub eff) indicator, is reproduced well. Conclusions. Reliable measurements and calculations of atomic data are a prerequisite for stellar-atmosphere modeling. Our oscillator-strength calculations have allowed, for the first time, Ge V and Ge VI lines to be successfully reproduced in a white dwarf s (RE 0503-289) spectrum and to determine its photospheric Ge abundance.

Rauch, T.↗
NASA NTRS2016

Spatial Sampling of Weather Data for Regional Crop Yield Simulations

Field-scale crop models are increasingly applied at spatio-temporal scales that range from regions to the globe and from decades up to 100 years. Sufficiently detailed data to capture the prevailing spatio-temporal heterogeneity in weather, soil, and management conditions as needed by crop models are rarely available. Effective sampling may overcome the problem of missing data but has rarely been investigated. In this study the effect of sampling weather data has been evaluated for simulating yields of winter wheat in a region in Germany over a 30-year period (1982-2011) using 12 process-based crop models. A stratified sampling was applied to compare the effect of different sizes of spatially sampled weather data (10, 30, 50, 100, 500, 1000 and full coverage of 34,078 sampling points) on simulated wheat yields. Stratified sampling was further compared with random sampling. Possible interactions between sample size and crop model were evaluated. The results showed differences in simulated yields among crop models but all models reproduced well the pattern of the stratification. Importantly, the regional mean of simulated yields based on full coverage could already be reproduced by a small sample of 10 points. This was also true for reproducing the temporal variability in simulated yields but more sampling points (about 100) were required to accurately reproduce spatial yield variability. The number of sampling points can be smaller when a stratified sampling is applied as compared to a random sampling. However, differences between crop models were observed including some interaction between the effect of sampling on simulated yields and the model used. We concluded that stratified sampling can considerably reduce the number of required simulations. But, differences between crop models must be considered as the choice for a specific model can have larger effects on simulated yields than the sampling strategy. Assessing the impact of sampling soil and crop management data for regional simulations of crop yields is still needed.

upscaling↗
NASA NTRS2017

Global Gridded Crop Model Evaluation: Benchmarking, Skills, Deficiencies and Implications.

Crop models are increasingly used to simulate crop yields at the global scale, but so far there is no general framework on how to assess model performance. Here we evaluate the simulation results of 14 global gridded crop modeling groups that have contributed historic crop yield simulations for maize, wheat, rice and soybean to the Global Gridded Crop Model Intercomparison (GGCMI) of the Agricultural Model Intercomparison and Improvement Project (AgMIP). Simulation results are compared to reference data at global, national and grid cell scales and we evaluate model performance with respect to time series correlation, spatial correlation and mean bias. We find that global gridded crop models (GGCMs) show mixed skill in reproducing time series correlations or spatial patterns at the different spatial scales. Generally, maize, wheat and soybean simulations of many GGCMs are capable of reproducing larger parts of observed temporal variability (time series correlation coefficients (r) of up to 0.888 for maize, 0.673 for wheat and 0.643 for soybean at the global scale) but rice yield variability cannot be well reproduced by most models. Yield variability can be well reproduced for most major producing countries by many GGCMs and for all countries by at least some. A comparison with gridded yield data and a statistical analysis of the effects of weather variability on yield variability shows that the ensemble of GGCMs can explain more of the yield variability than an ensemble of regression models for maize and soybean, but not for wheat and rice. We identify future research needs in global gridded crop modeling and for all individual crop modeling groups. In the absence of a purely observation-based benchmark for model evaluation, we propose that the best performing crop model per crop and region establishes the benchmark for all others, and modelers are encouraged to investigate how crop model performance can be increased. We make our evaluation system accessible to all crop modelers so that other modeling groups can also test their model performance against the reference data and the GGCMI benchmark.

wheat↗
NASA NTRSDate not supplied

GC/MS Method Development for Separating Lunar Volatile Ice Simulant Headspace Gases

Various investigators propose the lunar surface contains widely distributed volatiles, especially water- like species, i.e. OH and H2O. Surface volatiles are theorized to exist as a hydrated regolith layer, concentrated in extremely cold polar permanently shadowed regions (PSR), and/or solar wind implantation reservoirs in lunar glasses. The proposed sources of lunar surface volatiles range from cometary impacts, solar wind, or a supply present during moon formation. Future Artemis missions aim to collect and return the samples containing volatiles collected near lunar polar craters or PSRs. We, as advanced curation scientists, are responsible for developing techniques and methodologies for preserving returned sample integrity as much as possible. Pristine volatile-bearing samples are invaluable to the scientific community seeking to unravel the history of the solar system. Realistically, a sample will experience alteration during collection, transportation back to earth, and storage. The Planetary Exploration and Astromaterials Research Lab (PEARL) seeks to understand temperature and pressure effects on high-fidelity volatile-containing regolith simulants, the foundation for the future of cold curation. This abstract outlines the separation, identification, and quantification of headspace gases over volatile ice feed stock material using gas chromatography/mass spectrometry (GC/MS). Preliminary objectives concentrated on sample handling, reproducibility, and understanding the elution characteristics for each analyte. Initial GC/MS method development experiments utilized diluted static headspace sample preparation. Diluted samples were used because sampling headspace gases directly from a vial containing liquid analyte resulted in overloading of the column and detector. Overloading is evident based on chromatogram peak shapes and instrument contamination, or carry over, between experiments. A mixture of three alcohols were used for a majority of the sample handling and reproducibility studies. Reproducibility was tested via multiple users, calibration curves, and check standards. Stock solutions of condensed lunar volatile analytes included methanol, ammonia in methanol, hydrogen sulfide in water, and an equal volume mixture of methanol, ethanol, and isopropanol. Current samples use room air as the headspace sample matrix, however future experiments will incorporate an inert purge gas, such as argon or nitrogen. Three mL of each analyte solution were capped in separate 20 mL crimp top GC vials. Dilutions were carried out by removing an aliquot of headspace gases with a calibrated 1 mL gastight syringe and immediately transferring to a 20 mL capped crimp top vial. The GC/MS is a Thermo Fisher Trace 1310/ISQ 7000 with a TriPlus RSH autosampler and split/splitless injector module. The experiments outlined in this abstract use the following hardware: a 2.5 mL gastight headspace syringe tool, 1 mm ID x 78.5 mm length ultra-inert straight injection liner, and a TG-BondQ 30 m × 0.32 mm × 10 μm column. Various parameters, such as hardware selection and the temperature, pressure, and split ratio set points, continue to evolve as the overall experiment is refined. Diluted headspace chromatograms were collected for the individual stock solutions. Retention times, peak shapes, and mass spectra were evaluated and added to the data processing method for each molecule of interest. Figure 1 shows the total ion chromatograms for the three major lunar volatile simulant stock solutions: methanol, 7 N ammonia in methanol, and 0.4% hydrogen sulfide in water. Tailing peak shapes for ammonia (2.98 min rt) and water (4.06 min rt) indicate the molecules are not properly eluting from the selected column with the current separation method. Additionally, hydrogen sulfide and ammonia have overlapping peak windows, which could impact quantification. Ongoing experiments aim to address the peak shape and overlapping via the separation method and hardware selection. Sample preparation reproducibility experiments used stock solution containing equal volumes of a non- interactive mixture of methanol, ethanol, and isopropanol. Mass spectrum ion traces were used to identify and quantify all three alcohols. Peaks were automatically detected, identified, and integrated through the mass spectra detection and processing parameters. Calibration response curves and check standards were used to evaluate the validity of the sample preparation procedure. Figure 2 shows the methanol chromatogram peak area versus total headspace dilution volume transferred from the alcohol mixture vial. The calibration response curves and check standards validate sample preparation procedure. Continuing data analysis efforts are working towards correlating the peak area and instrument response factor to the headspace analyte concentration and condensed phase composition. Static headspace gas chromatography theory relies on Dalton’s law, Raoult’s law, Henry’s Law, and the Kolb and Ettre equation to associate peak area to the analyte composition in a non-ideal solution. Equation 1 is a simplified expression derived from the aforementioned theories. Future experiments involve liquid injections of the individual stock solutions, liquid and headspace analysis of various stock solution combinations, and the addition of regolith simulants to the mixtures. Temperature is another variable expected to affect reaction rates and will be explored.

Cecilia L. Amick↗
DOE OSTI2026

Toward the Neutrino Discovery Platform: An Auditable, Uncertainty-Bearing Toolchain for MINERvA Open-Data Cross-Section Analysis

The Neutrino Discovery Platform (NDP) aims to accelerate DUNE-era science by making the neutrino program's existing datasets analyzable through fast, reproducible, and auditable workflows. We report a working version of two of its layers, data curation and agentic orchestration, built and tested end to end on MINERvA open data. The guiding lesson throughout is that a cross section is a measurement, and not just a plotted shape, only if it carries a defensible systematic-uncertainty budget, a trustworthy unfolding, and a reproducible record. Using a single medium-energy playlist pair from the MINERvA open-data release (about $2.05\times10^{17}$ protons on target of data), we first reproduced the shapes of two published charged-current inclusive $\nu_\mu$ measurements through a complete extraction ladder: selection, background subtraction, D'Agostini unfolding, efficiency correction, and flux normalization. These shape-level reproductions ran and tracked the published results, but they lacked the systematic-uncertainty machinery that defines a MINERvA cross section. To supply it, we vendored and built the MINERvA Analysis Toolkit and developed a many-universe systematic-uncertainty tool that produces a portable covariance artifact, a parallel event-loop runner, and a per-run auditability harness. Validated against a published covariance release, the toolchain reproduces the released statistical, flux, and muon-energy-scale terms and shows that they account for roughly 63\% of the total variance, with the remainder unreleased. Using this same infrastructure, we then performed a measurement of our own design, the hadronic recoil-energy distribution of low-energy ($E_\nu<2.5$~GeV) charged-current inclusive events, and found data/simulation shape agreement of $\chi^2/\mathrm{ndf}=1.26$. Together these results show that the platform supports original physics and not only reproductions.

Breaux, Auto [Tulane U. (main)]↗
DOE OSTI2022

Anomalous structural recovery in the near glass transition range in a polymer glass: Data revisited in light of temperature variability in vacuum oven‐based experiments*

Abstract There is considerable interest in data reported by Cangialosi et al. [Cangialosi, D., Boucher, V. M., Alegría, A., & Colmenero, J. (2013). Physical Review Letters , 111 (9), 095701] that purportedly showed unusual irregular responses in the isothermal structural recovery behavior of glasses because the observation challenges the general view of the dynamics being smooth functions of time as observed in extensive dilatometric experiments by Kovacs [Kovacs, A. J. (1964). In Fortschritte der Hochpolymeren‐Forschung , 3 , 394–507. Springer, Berlin, Heidelberg.] in the 1960s. The reported data show an apparent two mechanism structural (enthalpy) recovery at aging temperatures ranging from 6 to 17 K below which is also controversial as it could not be reproduced in work from Koh and Simon [Koh, Y. P., & Simon, S. L. (2013). Macromolecules , 46(14), 5815–5821.] where a second plateau observed in the enthalpy loss curve was not reproduced in nearly 1‐year‐long experiments at 15 K below . Therefore, it is important to determine what might be possible reasons for the apparent non‐smooth relaxations. Here, we examine the possibility that the poor temperature control of a typical vacuum oven could explain the anomalous results. We used the Tool–Narayanaswamy–Moynihan (TNM) model of structural recovery to calculate the effect of typical vacuum oven temperature variation on the structural recovery of polystyrene and find that we can reproduce the reported experimental results. The issue of temperature control using a vacuum oven is discussed, and data from thermocouple measurements of temperature at various locations inside a vacuum oven are shown.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Isotopic analysis of sub-nanogram neodymium loads using new ATONA™ amplifiers

The long-term gains of the ATONA™ amplifiers are shown to have less than 1 ppm variability. Repeat measurements of the JNdi-1 reference material demonstrate the ability of the ATONA™ amplification technology to make measurements of 143 Nd/ 144 Nd ratios with 23 ppm external reproducibility on 100-pg loads. The effect of increasing integration time on analytical reproducibility is also displayed as increasing integration time from 10 to 30s reduced the external measurement uncertainty from 37 to 23 ppm. These measurements represent an improvement of more than a factor of 3 in external measurement reproducibility relative to previously published 143 Nd/ 144 Nd measurements of 100-pg loads. Furthermore, this new technology will allow for the measurement of smaller samples for precise isotope ratios and open new avenues of research in the geological and nuclear forensic communities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Logarithmic accuracy of parton showers: a fixed-order study

We formulate some first fundamental elements of an approach for assessing the logarithmic accuracy of parton-shower algorithms based on two broad criteria: their ability to reproduce the singularity structure of multi-parton matrix elements, and their ability to reproduce logarithmic resummation results. We illustrate our approach by considering properties of two transverse-momentum ordered final-state showers, examining features up to second order in the strong coupling. In particular we identify regions where they fail to reproduce the known singular limits of matrix elements. Here, the characteristics of the shower that are responsible for this also affect the logarithmic resummation accuracies of the shower, both in terms of leading (double) logarithms at subleading N C and next-to-leading (single) logarithms at leading N C .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Development and validation of an automatic method to detect the latest contracting viable left ventricular segments to assist guide CRT therapy from gated SPECT myocardial perfusion imaging

Objectives: The purpose of this study is to use ECG-gated SPECT MPI to detect the latest contracting viable left ventricular (LV) segments to help guide the LV probe placement used in CRT therapy and to validate segment selection against the visual integration method by experts. Methods: For each patient, the resting ECG-gated SPECT MPI short-axis images were sampled in 3D to generate a polar map of the perfusion distribution used to determine LV myocardial viability, and to measure LV synchronicity using our phase analysis tool. In the visual integration method, two experts visually interpreted the LV viability and mechanical dyssynchrony from the short-axis images and polar maps of viability and phase, to determine the latest contracting viable segments using the 17-segment model. In the automatic method, the apical segments, septal segments, and segments with more than 50% scar were excluded as these are not candidates for CRT LV probe placement. Amongst the remaining viable segments, the segments, whose phase angles were within 10° of the latest phase angle (the most delayed contracting segment), were identified for potential CRT LV probe placement and ranked based on the phase angles of the segments. Both methods were tested in 36 pre-CRT patients who underwent ECG-gated SPECT MPI. The accuracy was determined as the percent agreement between the visual integration and automatic methods. The automatic method was performed by a second independent operator to evaluate the inter-operator processing reproducibility. Results: In all the 36 patients, the LV lead positions of the 1st choices recommended by the automatic and visual integration methods were in the same segments in 35 patients, which achieved an agreement rate of 97.2%. In the inter-operator reproducibility test, the LV lead positions of the 1st choices recommended by the two operators were in the same segments in 25 patients, and were in the adjacent segments in 7 patients, which achieved an overall agreement of 88.8%. Conclusions: An automatic method has been developed to detect the latest contracting viable LV segments to help guide the LV probe placement used in CRT therapy. The retrospective clinical study with 36 patients suggests that this method has high agreement against the visual integration method by experts and good inter-operator reproducibility. Consequently, this method is promising to be a clinical tool to recommend the CRT LV lead positions.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2024

Computational analysis of flame initiation, quenching, and re-ignition in a prechamber natural gas engine under varying EGR-dilution levels

The on-road natural-gas (NG) fueled transportation relies on stoichiometric spark-ignition engines for the advantages of simple after-treatment system despite the efficiency penalty relative to lean-burn combustion strategies. Exhaust gas recirculation (EGR) has the potential to reduce this efficiency gap at low to moderate loads without the need for complex lean-exhaust aftertreatment systems. However, EGR dilution leads to reduced combustion stability and increased cycle-to-cycle variability. A promising technology that has the potential to achieve reliable operation under diluted conditions is the prechamber ignition (or turbulent jet ignition) which uses chemically active turbulent jets generated from combustion inside a prechamber to initiate, stabilize and accelerate combustion of the mixture inside the main chamber. The present work focusses on developing a RANS-based CFD approach to accurately reproduce in-cylinder phenomena in a stoichiometric NG prechamber-assisted heavy-duty engine without relying on complex combustion models that account for turbulence-chemistry interactions. This is necessary because reactive prechamber jets at high EGR dilution tend to extinguish while emerging into the main chamber, which is followed by a phase of re-ignition — a phenomenon that conventional G-equation or well-stirred reactor combustion models cannot reproduce. With addition of a damping multiplier to the well-stirred reactor model, the predictions are seen to show good agreement with experimental pressure evolution and combustion images acquired from a single cylinder Cummins N-14 optical diesel engine retrofitted with a prechamber ignition system. Model predictions of local heat release in the flame and temperature evolution inside the flame are used to investigate combustion dynamics in the prechamber and the main chamber. It is seen that the well-stirred reactor model with the inclusion of damping is able to reproduce the temporary reduction in heat release within the flame, which can be considered equivalent to quenching of jets, and the subsequent re-ignition of the flame inside the main chamber. The delay between quenching and re-ignition depends on the amount of dilution, as explained by an illustration of flame evolution in a Borghi diagram.

Prechamber ignition↗
DOE OSTI2021

Magnetoelastic transition and magnetocaloric effect in induction melted Fe 100-x Rh x bulk alloys with x = 50, 51

Magnetoelastic transitions (METs) in bulk in nearly equiatomic Fe-Rh alloys produced by arc melting may show poor reproducibility related to insufficient chemical homogeneity and presence of impurity phases in variable concentrations. To better understand the synthesis conditions that reliably yield bulk FeRh materials with reproducible MET characteristics, Fe 100-x Rh x alloys with x = 50, 50.5 and 51 at. % were prepared by induction melting and thermal annealing under identical conditions. The fabricated samples were cut into several slices, followed by characterization of METs in each of the slices using isothermal and isofield magnetization measurements, differential scanning calorimetry, and direct measurements of the magnetocaloric effect. All of the slices exhibit METs between the AFM and FM states, but the transitions are abrupt with nearly the same change of magnetization, ΔM, when x = 50.5 and 51, whereas for the x = 50 alloy the transition spreads over a wide temperature interval and ΔM may fluctuate by as much as 10 % from one specimen to another. A comparison of the magnetocaloric responses of x = 50 and 51 materials is presented. The clearly different effect of the magnetic field on the transition in both directions leads to significant differences in the reversibility and maximum values of the magnetic field-induced entropy and adiabatic temperature changes, as well as average hysteresis losses. In terms of reproducibility, our results suggest that induction melting is a more appropriate technique to prepare these binary alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Numerical study of the ignition behavior of a post-discharge kernel in a turbulent stratified crossflow

Ensuring robust ignition is critical for the operability of aeronautical gas-turbine combustors. For ignition to be successful, an important aspect is the ability of the hot gas generated by the spark discharge to initiate combustion reactions, leading to the formation of a self-sustained ignition kernel. This study focuses on this phenomena by performing simulations of kernel ignition in a crossflow configuration that was characterized experimentally. First, inert simulations are performed to identify numerical parameters correctly reproducing the kernel ejection from the ignition cavity, which is here modeled as a pulsed jet. In particular, the kernel diameter and the transit time of the kernel to the reacting mixture are matched with measurements. Considering stochastic perturbations of the ejection velocity of the ignition kernel, the variability of the kernel transit time is also reproduced by the simulations. Subsequently, simulations of a series of ignition sequences are performed with varying equivalence ratio of the fuel-air mixture in the crossflow. The numerical results are shown to reproduce the ignition failure that occurs for the leanest equivalence ratio ($φ = 0.6$). For higher equivalence ratios, the simulations are shown to capture the sensitivity of the ignition to the equivalence ratio, and the kernel successfully transitions into a propagating flame. Significant stochastic dispersion of the ignition strength is observed, which relates to the variability of the transit time of the kernel to the re- active mixture. An analysis of the structure of the ignition kernel also highlights the transition towards a self-propagating flame for successful ignition conditions.

42 ENGINEERING↗
DOE OSTI2024

Stochastic modal velocity field in rough-wall turbulence

Stochastically generated instantaneous velocity profiles are used to reproduce the outer region of rough-wall turbulent boundary layers in a range of Reynolds numbers extending from the wind tunnel to field conditions. Each profile consists in a sequence of steps, defined by the modal velocities and representing uniform momentum zones (UMZs), separated by velocity jumps representing the internal shear layers. Height-dependent UMZ is described by a minimal set of attributes: thickness, mid-height elevation, and streamwise (modal) and vertical velocities. These are informed by experimental observations and reproducing the statistical behaviour of rough-wall turbulence and attached eddy scaling, consistent with the corresponding experimental datasets. Sets of independently generated profiles are reorganized in the streamwise direction to form a spatially consistent modal velocity field, starting from any randomly selected profile. The operation allows one to stretch or compress the velocity field in space, increases the size of the domain and adjusts the size of the largest emerging structures to the Reynolds number of the simulated flow. By imposing the autocorrelation function of the modal velocity field to be anchored on the experimental measurements, we obtain a physically based spatial resolution, which is employed in the computation of the velocity spectrum, and second-order structure functions. The results reproduce the Kolmogorov inertial range extending from the UMZ and their attached-eddy vertical organization to the very-large-scale motions (VLSMs) introduced with the reordering process. The dynamic role of VLSM is confirmed in the –u'w' co-spectra and in their vertical derivative, representing a scale-dependent pressure gradient contribution.

42 ENGINEERING↗
DOE OSTI2026

Enhancing Sensitivity in Targeted Single-Cell Proteomics by Coupling a Dual Ion Funnel Interface with Triple Quadrupole Mass Spectrometer

Single-cell proteomics (SCP) has emerged as a powerful approach for understanding cellular heterogeneity and biological processes at unprecedented resolution. However, the extremely limited protein content of individual cells (femtogram to picogram levels) pushes current mass spectrometry instrumentation to its sensitivity limits, creating a critical analytical bottleneck. While selected reaction monitoring (SRM) using triple quadrupole (QqQ) instruments 1 offers advantages in sensitivity and reproducibility for targeted proteomics quantification, SRM still struggles with sensitivity for quantification of moderate- or low-abundance proteins from single-cell sample amounts. Here, we report the development and systematic evaluation of a dual ion funnel interface designed to address the sensitivity limitation by significantly enhancing ion transmission efficiency in commercial QqQ mass spectrometers. The dual ion funnel interface, composed of a curved S-funnel followed by a conventional ion funnel, improves ion transmission efficiency while reducing chemical noise through selective ion focusing. The performance of the dual ion funnel interface was systematically compared to standard interface on a TSQ Vantage platform across samples with different levels of complexity. The dual funnel interface demonstrated to provide up to 25-fold improvement in sensitivity across a wide range of protein concentrations in different biological matrices (low complex mouse macrophage and high complex human cells). Critically, enhanced sensitivity was accompanied by increased analytical reproducibility with lower coefficient of variations. Most importantly, the dual funnel interface enabled reliable quantification of low-abundance proteins that were barely detectable or not detected by the standard interface, extending analysis to single-cell equivalent amounts while maintaining excellent reproducibility. These results demonstrate that the dual funnel interface addresses the critical bottleneck in quantitative targeted proteomics, providing a technological foundation for ultrasensitive targeted SCP that requires both high sensitivity and robust quantitative performance.

Min, Sehong↗
DOE OSTI2020

Minimal Experimental Bias on the Hydrogen Bond Greatly Improves Ab Initio Molecular Dynamics Simulations of Water

Experiment Directed Simulations (EDS) is a method within a class of techniques seeking to improve molecular simulations by minimally biasing the system Hamiltonian to reproduce certain experimental observables. In a previous application of EDS to ab initio molecular dynamics (AIMD) simulation based on electronic density functional theory (DFT), the AIMD simulations of water were biased to reproduce its experimentally derived solvation structure. In particular, by solely biasing the O-O pair correlation functions, other structural and dynamical properties that were not biased were improved. In this work, the hypothesis is tested that directly biasing the OH pair correlation (and hence the H-O∙∙∙H hydrogen bonding), will provide an even better improvement of DFT-based water properties in AIMD simulations. The logic behind this hypothesis is that for most electronic DFT descriptions of water the hydrogen bonding is known to be deficient due to anomalous charge transfer and over polarization in the DFT. Using recent advances to the EDS learning algorithm, we thus train a minimal bias on AIMD water that reproduces the O-H radial distribution function derived from the highly accurate MB-pol model of water. Finally, it is then confirmed that biasing the O-H pair correlation alone can lead to improved AIMD water properties, with structural and dynamical properties in even closer to experiment than the previous EDS-AIMD model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗