Numerical Prediction of Boundary Layer Flow Around An Immersed Plate in Carbopol Yield Stress Fluid
2022 Society for Industrial and Applied Mathematics Annual Meeting, Pittsburgh, PA, July 11-15, 2022
Konan, Arthur↗
SEARCH · Engineering Papers
Results for “Yield Prediction”
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.
Investigating strangeness enhancement with multiplicity in pp collisions using angular correlations
A study of strange hadron production associated with hard scattering processes and with the underlying event is conducted to investigate the origin of the enhanced production of strange hadrons in small collision systems characterised by large charged-particle multiplicities. For this purpose, the production of the single-strange meson ${\textrm{K}}_{\textrm{S}}^0$ and the double-strange baryon Ξ ± is measured, in each event, in the azimuthal direction of the highest-p T particle (“trigger” particle), related to hard scattering processes, and in the direction transverse to it in azimuth, associated with the underlying event, in pp collisions at $\sqrt{s}$ = 5.02 TeV and $\sqrt{s}$ = 13 TeV using the ALICE detector at the LHC. The per-trigger yields of ${\textrm{K}}_{\textrm{S}}^0$ and Ξ ± are dominated by the transverse-to-leading production (i.e., in the direction transverse to the trigger particle), whose contribution relative to the toward-leading production is observed to increase with the event charged-particle multiplicity. The transverse-to-leading and the toward-leading Ξ ± /${\textrm{K}}_{\textrm{S}}^0$ yield ratios increase with the multiplicity of charged particles, suggesting that strangeness enhancement with multiplicity is associated with both hard scattering processes and the underlying event. The relative production of Ξ ± with respect to ${\textrm{K}}_{\textrm{S}}^0$ is higher in transverse-to-leading processes over the whole multiplicity interval covered by the measurement. The ${\textrm{K}}_{\textrm{S}}^0$ and Ξ ± per-trigger yields and yield ratios are compared with predictions of three different phenomenological models, namely Pythia8.2 with the Monash tune, Pythia8.2 with ropes and EPOS LHC. The comparison shows that none of them can quantitatively describe either the transverse-to-leading or the toward-leading yields of ${\textrm{K}}_{\textrm{S}}^0$ and Ξ ± .
hadron-hadron scattering↗
Robust verification of stochastic simulation codes
We introduce a robust verification tool for computational codes, which we call Stochastic Robust Extrapolation based Error Quantification (StREEQ). Unlike the prevalent Grid Convergence Index (GCI) [1] method, our approach is suitable for both stochastic and deterministic computational codes and is generalizable to any number of discretization variables. Building on ideas introduced in the Robust Verification [2] approach, we estimate the converged solution and orders of convergence with uncertainty using multiple fits of a discretization error model. In contrast to Robust Verification, we perform these fits to many bootstrap samples yielding a larger set of predictions with smoother statistics. Here, bootstrap resampling is performed on the lack-of-fit errors for deterministic code responses, and directly on the noisy data set for stochastic responses. This approach lends a degree of robustness to the overall results, capable of yielding precise verification results for sufficiently resolved data sets, and appropriately expanding the uncertainty when the data set does not support a precise result. For stochastic responses, a credibility assessment is also performed to give the analyst an indication of the trustworthiness of the results. Furthermore, this approach is suitable for both code and solution verification, and is particularly useful for solution verification of high-consequence simulations..
97 MATHEMATICS AND COMPUTING↗
Watching a hydroperoxyalkyl radical (•QOOH) dissociate
A prototypical hydroperoxyalkyl radical (•QOOH) intermediate, transiently formed in the oxidation of volatile organic compounds, was directly observed through its infrared fingerprint and energy-dependent unimolecular decay to hydroxyl radical and cyclic ether products. Direct time-domain measurements of •QOOH unimolecular dissociation rates over a wide range of energies were found to be in accord with those predicted theoretically using state-of-the-art electronic structure characterizations of the transition state barrier region. Unimolecular decay was enhanced by substantial heavy-atom tunneling involving O-O elongation and C-C-O angle contraction along the reaction pathway. Master equation modeling yielded a fully a priori prediction of the pressure-dependent thermal unimolecular dissociation rates for the •QOOH intermediate—again increased by heavy-atom tunneling—which are required for global models of atmospheric and combustion chemistry.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Investigation of post-injection strategies for diesel engine Catalyst Heating Operation using a vapor-liquid-equilibrium-based spray model
Most multidimensional engine simulations spend much time solving for non-equilibrium spray dynamics (atomization, collision, vaporization). However, their accuracy is limited by significant grid dependency, and the need for extensive calibration. This is critical for modeling cold-start diesel fuel post injections, which occur at low temperatures and pressures, far from typical model validation ranges. At the same time, resolving micron-scale spray phenomena would render full Eulerian multiphase calculations prohibitive. In this study, an improved phase equilibrium based approach was implemented and assessed for simulating diesel catalyst heating operation strategies. A phase equilibrium solver based on the model by Yue and Reitz [1] was implemented: a fully multiphase CFD solver is employed with an engineering-size engine grid, and fuel injection is modeled using the standard Lagrangian parcels approach. Mass and energy from the liquid parcels are released to the Eulerian multiphase mixture according to an equilibrium-based liquid jet model. An improved phase equilibrium solver was developed to handle large real-gas mixtures such as those from accurate chemical kinetics mechanisms. The liquid-jet model was improved such that momentum transfer to the Eulerian solver better reproduces the physical spray jet structure. Validation of liquid/vapor penetration predictions showed that the model yields accurate results with very limited tuning and low sensitivity to the few calibration constants. In-cylinder simulations of diesel catalyst heating operation strategies showed that capturing spray structure is paramount when short, transient injection pulses and low temperatures are present. Furthermore, the EP model provides improved predictions of post-injection spray structure and ignitability, while conventional spray modeling does not capture the increase of liquid penetration during the expansion stroke. Finally, the only important EP model calibration constant, Cliq, does not affect momentum transfer, but it changes the local charge cooling distribution through the local energy transfer, which makes it candidate to additional research. The results confirm that non-equilibrium spray processes do not need to be resolved in engineering simulations of high-pressure diesel sprays.
42 ENGINEERING↗
Conceptual Design of a Thermal–to-14-MeV Neutron Conversion Device for Use in the Advanced Test Reactor
Here, this paper details the conceptual design of a thermal-to-14-MeV neutron converter consisting of a mixture of lithium and deuterium in a blanket material. Such a device operates a two-step reaction, first generating tritons via thermal neutron absorption in the tritium breeding material, and in the second step, high-energy neutrons are produced either via deuterium-tritium fusion reaction or with tritium reacting with lithium. A thermal-to-14-MeV neutron converter significantly hardens the neutron spectrum by virtually removing thermal neutrons and adding a high-energy 14 MeV component to the neutron spectrum. While similar concepts have been previously proposed and tested in other reactors, the unique characteristics of the Advanced Test Reactor (ATR), namely its important thermal flux (up to $10^{15}$ n.cm$^{-2}$.s$^{-1}$) makes it markedly attractive for obtaining a very large fast neutron flux, usable for irradiation studies under neutron flux conditions prototypical of fusion reactors. The paper provides a description of a new computational scheme developed for handling the coupled neutron-triton transport mechanism using the Geant4 toolkit. Resulting neutron spectra and high-energy neutron yields are summarized for different irradiation positions and potential neutron breeder materials. Maximum predicted thermal-to-14-MeV neutron yields are on the order of $2.10^{-4}$, which is consistent with previous studies found in the literature. Thus, when placed inside the ATR, such a neutron converter will be providing the largest high-energy neutron source available for activation and irradiation studies of materials foreseen for use in fusion reactors. Future steps will involve qualifying the computational scheme using the ATR critical facility using activation foil measurements.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Improvements to the macroscopic-microscopic approach of nuclear fission
The well-established macroscopic-microscopic (mac-mic) description of nuclear fission enables the prediction of fission-fragment yields for a broad range of fissioning systems. In this work, we present several key enhancements to this approach. We improve upon the microscopic sector of nuclear potential-energy surfaces by magnifying the resolution of the Lipkin-Nogami equations and strengthening the Strutinsky procedure, thus reducing spurious effects from the continuum. We further present a novel deterministic method for calculating fission dynamics under the assumption of strongly damped nucleonic motion. Our technique directly determines the evolution of the scissioned shape distribution according to the number of random-walk steps rather than the statistical accumulation of fission events. We show that our new technique is equivalent to the Metropolis random walk pioneered over the past decade by Randrup and colleagues. It further improves upon it because we remove the need for altering the nuclear landscape via a biased potential. With our final improvement, we calculate fission fragments mass and charge distributions using particle number projection, which affords the simultaneous calculation of both mass- and charge-yield distributions. Fission fragments are thus calculated from the quantum-mechanical A -body states of the potential-energy surface rather than from the collective mass asymmetry variable α g of the finite-range liquid-drop model used in past work. We highlight the success of our enhancements by predicting the odd-even staggering and the charge polarization for the neutron-induced fission of 233 U and 235 U .
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
A Simplified Method for Predicting Shaker Voltage in IMMATs
Impedance Matched Multi-Axis Tests (IMMATs) can replicate in-service vibration induced stress more accurately than single axis shaker table tests as they can better match a part’s operational boundary conditions and excite it in multiple degrees of freedom simultaneously. Here, the shakers used in IMMATs are less powerful than shaker tables, so shaker force limits can be exceeded during tests if they are not placed adequately for the desired environment. The ability to predict shaker voltage and force before performing a test is, therefore, helpful in selecting shaker locations so that their limits are not exceeded. In this study, electrodynamic shakers were modeled as discrete electromechanical systems, and the shaker parameters were chosen to match experimentally obtained acceleration/voltage frequency response functions (FRFs). These models were coupled to a finite element model of the device under test (DUT) via dynamic substructuring, and the substructured model was demonstrated to accurately predict shaker voltage as well as the error in reproducing the environment at multiple accelerometer locations. A simple method called the FRF Multiplication method, in which the FRF of the substructured system is approximated as the product of two separate FRFs of the shaker and DUT respectively, was proposed and applied to the same system, yielding similar voltage and error predictions to those obtained using substructuring. Simple case studies were presented to explore the applicability of the proposed method, and it was demonstrated to have similar accuracy to the substructuring method in a range of cases. Additionally, we showed that while it was not possible to derive a unique model of the shakers from acceleration/voltage FRFs alone, the models that could be obtained were sufficient to predict test error almost perfectly and shaker voltage with less than 40 percent error.
42 ENGINEERING↗
Improved deep learning prediction of antigen–antibody interactions
Identifying antibodies that neutralize specific antigens is crucial for developing effective immunotherapies, but this task remains challenging for many target antigens. The rise of deep learning–based computational approaches presents a promising avenue to address this challenge. Here, we assess the performance of a deep learning approach through two benchmark tests aimed at predicting antibodies for the receptor-binding domain of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein. Three different strategies for constructing input sequence alignments are employed for predicting structural models of antigen–antibody complexes. In our initial testing set, which comprises known experimental structures, these strategies collectively yield a significant top-ranked prediction for 61% of cases and a success rate of 47%. Notably, one strategy that utilizes the sequences of known antigen binders outperforms the other two, achieving a precision of 90% in a subsequent test set of ~1,000 antibodies, balanced between true and control antibodies for the antigen, albeit with a lower recall of 25%. Our results underscore the potential of integrating deep learning methods with single B cell sequencing techniques to enhance the prediction accuracy of antigen–antibody interactions.
Science & Technology - Other Topics↗
Development of a deep potential model for F and CF 2 etching of Si and SiO 2
An understanding of plasma-surface interactions at increasingly smaller scales is invaluable for the development of novel technologies and processing techniques. Molecular dynamics (MD) simulations can provide insights into atomic-scale interactions, though they are restricted by the availability of interatomic potentials. Machine learning methods, such as Deep Potential Molecular Dynamics (DeepMD), provide a systematic framework for the development of accurate and flexible ab initio-based models. In this work, we develop DeepMD models for the ion-enhanced etching of Si and SiO 2 by F and CF 2 radicals. We employ an active learning process to expand the data set on which the model is trained and demonstrate its effect on the model accuracy. The DeepMD results are compared to data from classical MD simulations and experiments. Physical sputtering yields of SiO 2 by Ar + ions show good agreement with previous simulation results using conventional interatomic potentials, though the predicted depth profiles are different. Etching yields are calculated as a function of ion energy and neutral to ion flux ratio for the Ar + ion-enhanced etching of SiO 2 and Si by F atoms, as well as for etching of SiO 2 by CF 2 radicals, showing reasonable agreement with experimental data. Finally, an ion-enhanced surface kinetic model is fitted to the DeepMD etch yields, and the fitted parameters are compared to quantities computed directly from DeepMD simulations. This study illustrates how molecular dynamics simulations using machine learning potentials can provide an accurate model of etching processes relevant to device manufacturing.
Kounis-Melas, Andreas [Princeton Univ., NJ (United↗
Predicting solid state material platforms for quantum technologies
Semiconductor materials provide a compelling platform for quantum technologies (QT). However, identifying promising material hosts among the plethora of candidates is a major challenge. Therefore, we have developed a framework for the automated discovery of semiconductor platforms for QT using material informatics and machine learning methods. Different approaches were implemented to label data for training the supervised machine learning (ML) algorithms logistic regression, decision trees, random forests and gradient boosting. We find that an empirical approach relying exclusively on findings from the literature yields a clear separation between predicted suitable and unsuitable candidates. In contrast to expectations from the literature focusing on band gap and ionic character as important properties for QT compatibility, the ML methods highlight features related to symmetry and crystal structure, including bond length, orientation and radial distribution, as influential when predicting a material as suitable for QT.
36 MATERIALS SCIENCE↗
Uncertainty-aware particle segmentation for electron microscopy at varied length scales
Electron microscopy is indispensable for examining the morphology and composition of solid materials at the sub-micron scale. To study the powder samples that are widely used in materials development, scanning electron microscopes (SEMs) are increasingly used at the laboratory scale to generate large datasets with hundreds of images. Parsing these images to identify distinct particles and determine their morphology requires careful analysis, and automating this process remains challenging. In this work, we enhance the Mask R-CNN architecture to develop a method for automated segmentation of particles in SEM images. We address several challenges inherent to measurements, such as image blur and particle agglomeration. Moreover, our method accounts for prediction uncertainty when such issues prevent accurate segmentation of a particle. Recognizing that disparate length scales are often present in large datasets, we use this framework to create two models that are separately trained to handle images obtained at low or high magnification. By testing these models on a variety of inorganic samples, our approach to particle segmentation surpasses an established automated segmentation method and yields comparable results to the predictions of three domain experts, revealing comparable accuracy while requiring a fraction of the time. These findings highlight the potential of deep learning in advancing autonomous workflows for materials characterization.
36 MATERIALS SCIENCE↗
Accurate and uncertainty-aware multi-task prediction of HEA properties using prior-guided deep Gaussian processes
Surrogate modeling techniques have become indispensable in accelerating the discovery and optimization of high-entropy alloys (HEAs), especially when integrating computational predictions with sparse experimental observations. This study systematically evaluates the training and testing performance of four prominent surrogate models—conventional Gaussian processes (cGP), Deep Gaussian processes (DGP), encoder-decoder neural networks for multi-output regression and eXtreme Gradient Boosting (XGBoost)—applied to a hybrid dataset of experimental and computational properties of the 8-component HEA system Al-Co-Cr-Cu-Fe-Mn-Ni-V. We specifically assess their capabilities in predicting correlated material properties, including yield strength, hardness, modulus, ultimate tensile strength, elongation, and average hardness under dynamic/quasi-static conditions, alongside auxiliary computational properties. The comparison highlights the strengths of hierarchical deep modeling approaches in handling heteroscedastic, heterotopic, and incomplete data commonly encountered in materials science. Our findings illustrate that combined surrogate models such as DGPs infused with machine-learned priors outperform other surrogates by effectively capturing inter-property correlations and by assimilating prior knowledge. This enhanced predictive accuracy positions the combined surrogate models as powerful tools for robust and data-efficient materials design.
36 MATERIALS SCIENCE↗
Automated Calibration for Rapid Optical Spectroscopy Sensor Development for Online Monitoring
An automated platform has been developed to assist researchers in the rapid development of optical spectroscopy sensors to quantify species from spectral data. This platform performs calibration and validation measurements simultaneously. Real-time, in situ monitoring of complex systems through optical spectroscopy has been shown to be a useful tool; however, building calibration models requires development time, which can be a limiting factor in the case of radiological or otherwise hazardous systems. While calibration time can be reduced through optimized design of experiments, this study approached the challenge differently through automation. The ATLAS (Automated Transient Learning for Applied Sensors) platform used pneumatic control of stock solutions to cycle flow profiles through desired calibration concentrations for multivariate model construction. Additionally, the transients between desired concentrations based on flow calculations were used as validation measurements to understand model predictive capabilities. This automated approach yielded an incredible 76% reduction in model development time and a 60% reduction in sample volume versus estimated manual sample preparation and static measurements. The ATLAS system was demonstrated on two systems: a three-lanthanide system with Pr/Nd/Ho representing a use case with significant overlap or interference between analyte signatures and an alternate system containing Pr/Nd/Ni to demonstrate a use case in which broad-band corrosion species signatures interfered with more distinct lanthanide absorbance profiles. Both systems resulted in strong model prediction performance (RMSEP < 9%). Lastly, ATLAS was demonstrated as a tool to simulate process monitoring scenarios (e.g., column separation) in which models can be further optimized to account for day-to-day changes as necessary (e.g., baseline correction). Ultimately, ATLAS offers a vital tool to rapidly screen monitoring methods, investigate sensor fusion, and explore more complex systems (i.e., larger numbers of species).
47 OTHER INSTRUMENTATION↗
Genomic prediction of regional-scale performance in switchgrass ( Panicum virgatum ) by accounting for genotype-by-environment variation and yield surrogate traits
Switchgrass is a potential crop for bioenergy or carbon capture schemes, but further yield improvements through selective breeding are needed to encourage commercialization. To identify promising switchgrass germplasm for future breeding efforts, we conducted multisite and multitrait genomic prediction with a diversity panel of 630 genotypes from 4 switchgrass subpopulations (Gulf, Midwest, Coastal, and Texas), which were measured for spaced plant biomass yield across 10 sites. Our study focused on the use of genomic prediction to share information among traits and environments. Specifically, we evaluated the predictive ability of cross-validation (CV) schemes using only genetic data and the training set (cross-validation 1: CV1), a subset of the sites (cross-validation 2: CV2), and/or with 2 yield surrogates (flowering time and fall plant height). We found that genotype-by-environment interactions were largely due to the north–south distribution of sites. The genetic correlations between the yield surrogates and the biomass yield were generally positive (mean height r = 0.85; mean flowering time r = 0.45) and did not vary due to subpopulation or growing region (North, Middle, or South). Genomic prediction models had CV predictive abilities of –0.02 for individuals using only genetic data (CV1), but 0.55, 0.69, 0.76, 0.81, and 0.84 for individuals with biomass performance data from 1, 2, 3, 4, and 5 sites included in the training data (CV2), respectively. To simulate a resource-limited breeding program, we determined the predictive ability of models provided with the following: 1 site observation of flowering time (0.39); 1 site observation of flowering time and fall height (0.51); 1 site observation of fall height (0.52); 1 site observation of biomass (0.55); and 5 site observations of biomass yield (0.84). The ability to share information at a regional scale is very encouraging, but further research is required to accurately translate spaced plant biomass to commercial-scale sward biomass performance.
09 BIOMASS FUELS↗
Charge radii of exotic neon and magnesium isotopes
We compute the charge radii and ground-state energies of even-mass neon and magnesium isotopes from neutron number N = 8 to the dripline. Our calculations are based on nucleon-nucleon and three-nucleon potentials from chiral effective field theory that include Δ isobars. These potentials yield an accurate saturation point and symmetry energy of nuclear matter. We use the coupled-cluster method and start from an axially symmetric reference state. Binding energies and two-neutron separation energies largely agree with data, and the dripline in neon is accurate. The computed charge radii are accurate for many isotopes where data exist. Finer details, such as isotope shifts, however, are not accurately reproduced. These chiral potentials indicate a subshell closure at N = 14 for the radii (but not for two-neutron separation energies) and a decrease in charge radii at N = 8 (observed in neon and predicted for magnesium). Furthermore, they yield a continued increase of charge radii as neutrons are added beyond N = 14 yet underestimate the large increase at N = 20 in magnesium.
20 ≤ A ≤ 38↗
Non-linear relationships between daily temperature extremes and US agricultural yields uncovered by global gridded meteorological datasets
Global agricultural commodity markets are highly integrated among major producers. Prices are driven by aggregate supply rather than what happens in individual countries in isolation. Furthermore, estimating the effects of weather-induced shocks on production, trade patterns and prices hence requires a globally representative weather data set. Recently, two data sets that provide daily or hourly records, GMFD and ERA5-Land, became available. Starting with the US, a data rich region, we formally test whether these global data sets are as good as more fine-scaled country-specific data in explaining yields and whether they estimate similar response functions. While GMFD and ERA5-Land have lower predictive skill for US corn and soybeans yields than the fine-scaled PRISM data, they still correctly uncover the underlying non-linear temperature relationship. All specifications using daily temperature extremes under any of the weather data sets outperform models that use a quadratic in average temperature. Correctly capturing the effect of daily extremes has a larger effect than the choice of weather data. In a second step, focusing on Sub Saharan Africa, a data sparse region, we confirm that GMFD and ERA5-Land have superior predictive power to CRU, a global weather data set previously employed for modeling climate effects in the region.
54 ENVIRONMENTAL SCIENCES↗
SO(3)-invariance of informed-graph-based deep neural network for anisotropic elastoplastic materials
This work examines the frame-invariance (and the lack thereof) exhibited in simulated anisotropic elasto-plastic responses generated from supervised machine learning of classical multi-layer and informed-graph-based neural networks, and proposes different remedies to fix this drawback. The inherent hierarchical relations among physical quantities and state variables in an elasto-plasticity model are first represented as informed, directed graphs, where three variations of the graph are tested. While feed-forward neural networks are used to train path-independent constitutive relations (e.g., elasticity), recurrent neural networks are used to replicate responses that depends on the deformation history, i.e. or path dependent. In dealing with the objectivity deficiency, we use the spectral form to represent tensors and, subsequently, three metrics, the Euclidean distance between the Euler Angles, the distance from the identity matrix, and geodesic on the unit sphere in Lie algebra, can be employed to constitute objective functions for the supervised machine learning. In this, the aim is to minimize the measured distance between the true and the predicted 3D rotation entities. Following this, we conduct numerical experiments on how these metrics, which are theoretically equivalent, may lead to differences in the efficiency of the supervised machine learning as well as the accuracy and robustness of the resultant models. Neural network models trained with tensors represented in component form for a given Cartesian coordinate system are used as a benchmark. Our numerical tests show that, even given the same amount of information and data, the quality of the anisotropic elasto-plasticity model is highly sensitive to the way tensors are represented and measured. The results reveal that using a loss function based on geodesic on the unit sphere in Lie algebra together with an informed, directed graph yield significantly more accurate rotation prediction than the other tested approaches.
42 ENGINEERING↗