Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical 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 307 records · Page 17

Prediction of pilot reserve attention capacity during air-to-air target tracking

Reserve attention capacity of a pilot was calculated using a pilot model that allocates exclusive model attention according to the ranking of task urgency functions whose variables are tracking error and error rate. The modeled task consisted of tracking a maneuvering target aircraft both vertically and horizontally, and when possible, performing a diverting side task which was simulated by the precise positioning of an electrical stylus and modeled as a task of constant urgency in the attention allocation algorithm. The urgency of the single loop vertical task is simply the magnitude of the vertical tracking error, while the multiloop horizontal task requires a nonlinear urgency measure of error and error rate terms. Comparison of model results with flight simulation data verified the computed model statistics of tracking error of both axes, lateral and longitudinal stick amplitude and rate, and side task episodes. Full data for the simulation tracking statistics as well as the explicit equations and structure of the urgency function multiaxis pilot model are presented.

Onstott, E. D.↗

Turbulent burning velocity of lean premixed hydrogen/air flames at engine conditions: Effects of turbulence intensity and length scale

For turbulent lean premixed hydrogen flames with strong thermodiffusively instabilities, most previous studies have focused on the influence of turbulence intensity, whereas the role of turbulence length scale is less well understood. Here, this study addresses this gap by conducting direct numerical simulations (DNS) of statistically planar turbulent premixed flames for a lean (ϕ=0.35) hydrogen/air mixture under independently varied turbulence intensity (u') and length scale (l T ) at engine-relevant thermodynamics conditions. Results show that as u' increases, the flame front becomes increasingly wrinkled, forming smaller cellular structures. In contrast, l T variations do not significantly alter the size of these structures. For the turbulent burning velocity (S T ), the normalized S T (i.e., S T /S L , where S L is the laminar flame speed) increases linearly with u', driven by both enhanced flame surface wrinkling (i.e., increased A T /A L ) and enhanced local burning rate (i.e., increased I 0 ). However, increasing l T reduces I 0 , despite a continued increase in A T /A L , resulting in only a marginal increase in S T /S L . To reveal the underlying mechanisms, especially the decreasing trend of I 0 with l T , local flame dynamics analyses are performed. It is found that as l T increases, the interaction between thermodiffusive effects and turbulence weakens due to the reduced tangential strain rate, while the flame curvature remains largely unchanged. This suppresses local reactivity enhancement and thus decreases I 0 , In contrast, an increase in u' enhances the interaction by amplifying both curvature fluctuation and tangential strain rate, leading to increased local reactivity (increased I 0 ). Finally, based on the DNS data, several new scaling models are proposed for the three global properties, S T /S L , A T /A L , and I 0 , and showed improvements compared to existing models. These findings provide new insights into the flame-turbulence interactions in thermodiffusively unstable hydrogen flames. The DNS dataset is also useful for the development of turbulent combustion models applicable to practical engine simulations.

Engine-relevant condition↗

Lattice-scale variations in viscosity are correlated with solution structure at mineral-water interfaces

At solid-liquid interfaces, the viscosity increases markedly from the bulk due to the collective interactions of ions and water molecules, influencing phenomena relevant to nanofluidics, colloidal dynamics, and electrochemistry. Here, in this study, we investigated dissipative forces at the boehmite-water interface using 3D atomic force microscopy. We observed an increase in interfacial solution viscosity, η, by 10-100-fold as the nanoprobe approached the surface in normal direction, with up to four oscillatory features showing average peaks of η/η bulk = 44-71. Moreover, the viscosity showed sub-nanometer variations within 0.5 nm from the interface, templated by the underlying crystal lattice and correlated with the interfacial solution structure. Beyond a near-wall region of approximately 1.2 nm, the dissipative response was comparable to that in bulk solution. Molecular dynamics simulations, along with statistical mechanical analyses, provided details on hydrodynamic structures near the interface. Specifically, the lattice-dependent dissipative responses are correlated with extensive hydrogen bonding by interfacial water molecules, which increased friction, particularly along the [001] direction. These results demonstrate how solution viscosity at mineral-water interfaces is anisotropic and correlated with the local solution structure, providing insights into the dynamics of nanocrystal attachment.

Viscosity↗

A Comparison of PV Resource Modeling for Sizing Microgrid Components

Microgrid systems are being deployed with increased frequency to meet critical building loads during unplanned power outages. Determining the appropriate type and size of components that make up a microgrid can greatly affect its ability to meet this need. This paper presents a new tool to appropriately size PV, battery, and generator capacities to meet all site loads, given a resilience goal. It uses a statistical model to simulate microgrid behavior under a large range of solar resource conditions, providing resilience planners with confidence that critical loads will be met even under extreme weather conditions. We compare this model with one that uses typical weather data for sizing microgrids, and quantify the resilience effects of using these different methods. Our method produces more conservative generator capacities and fuel requirements in an outage, predicting up to 30% more fuel required in the most extreme case. Microgrid systems that are designed using more conservative estimates are more likely to continue to serve load during emergency outages under a large range of conditions.

Newman, Sarah F.↗

Ripening of Rh Nanoparticle Catalysts in Reverse Water–Gas Shift via a Data-Driven Model Combining Physics, Theory, and Experiment

Degradation via sintering is an ongoing challenge that impedes the broad commercial success of supported metallic nanoparticle catalysts. To mitigate degradation via informed catalyst design and process operations, here we aim to disambiguate the underlying mechanisms of sintering by combining theory and experiment in a quantitative framework. While mechanistic sintering models exist, they only model a single sintering pathway, even though multiple sintering mechanisms can occur simultaneously or dominate at different stages of the process. Data-driven machine learning models have emerged as a means to represent complex processes through data regression. However, machine learning models have very large data needs and lack mechanistic insights due to their black-box encoding. To develop an interpretive model of catalyst degradation via sintering, we constructed a hybrid model combining mechanistic “physics-based” models and data-driven methods to obtain both reliable predictions and mechanistic insights regarding experimentally observed sintering phenomena. Focusing on nanoparticle sintering in the Rh–TiO 2 catalyst for the reverse water–gas shift (RWGS) reaction, the hybrid model couples a mechanistic term for Ostwald ripening with energy values calculated via density functional theory (DFT) with a parametric, data-driven discrepancy function term for unmodeled mechanisms. The hybrid model is trained using Bayesian inference with data collected from small-angle X-ray scattering (SAXS) in situ experiments wherein average nanoparticle diameter versus time was measured at three relevant operating temperatures. The calibrated hybrid model results show that an Ostwald ripening-only model parameterized with fixed DFT energies does not fully capture the time and temperature dependence of the SAXS-observed sintering kinetics, and that an additional functional contribution, or DFT energy calibration, is required to reconcile simulation and experiment. Analysis of the hybrid-model error confirms that the hybrid model outperforms both the purely mechanistic and purely data-driven alternatives in terms of expected predictive accuracy for time-evolving average particle sizes. Furthermore, the results support the hypothesis that the Ostwald ripening mechanism is less important for explaining the sintering phenomena as operating temperature increases under an assumed fixed DFT parameterization. This could be explained in one of two ways: either latent, unmodeled sintering mechanisms dominate at higher temperatures, or the DFT uncertainty increases with temperature. The proposed modeling approach directly links theory to experiments and simulations via a statistical hybrid modeling framework and can be extended to other catalytic systems to improve predictive models and mechanistic understanding.

Bayesian hybrid modeling↗

A Machine‐Learning‐Assisted Stochastic Cloud Population Model as a Parameterization of Cumulus Convection

Abstract A machine‐learning‐assisted stochastic cloud population model is coupled with the Advanced Research Weather Research and Forecasting (WRF) model to represent fluctuations in the cloud‐base mass flux associated with the life cycles and interactions among cumulus convection cells. In this cloud population model, the size distribution and the associated cloud‐base mass flux of the convective cells are related to their previous state and to the change in the total convective area via a transition function. The convective area tendency in turn is assumed to depend on the cloud‐base mass flux that is resolved by the host WRF model. The transition function is represented by a single hidden‐layer neural network trained by the evolution of convective cell size distributions in a 1‐km grid‐spacing WRF simulation run over the Australian Monsoon region. At every grid point of the host model, the cloud population model predicts the cell size and cloud‐base mass flux distributions from which a random sample of cells is fed to an entraining parcel model that calculates precipitation as well as the associated liquid water potential temperature and total moisture tendencies. These tendencies are averaged over the cells and provided to the host model. Several regional simulations are performed over tropical and midlatitude domains to test this as a potential approach to scale‐aware parameterization. It is shown that such an approach could be a new promising path to simulating realistic precipitation statistics and propagation of precipitation associated with the Madden‐Julian Oscillation while maintaining realistic depictions of the diurnal cycle over both land and ocean.

54 ENVIRONMENTAL SCIENCES↗

The local wavenumber model for computation of turbulent mixing

We present an overview of the current status in the development of a two-point spectral closure model for turbulent flows, known as the local wavenumber (LWN) model. The model is envisioned as a practical option for applications requiring multi-physics simulations in which statistical hydrodynamics quantities such as Reynolds stresses, turbulent kinetic energy, and measures of mixing such as density-correlations and mix-width evolution, need to be captured with relatively high fidelity. In this review, we present the capabilities of the LWN model since it was first formulated in the early 1990s, for computations of increasing levels of complexity ranging from homogeneous isotropic turbulence, inhomogeneous and anisotropic single-fluid turbulence, to two-species mixing driven by buoyancy forces. The review concludes with a discussion of some of the more theoretical considerations that remain in the development of this model.

42 ENGINEERING↗

Quantal diffusion approach for multinucleon transfer processes in the 58,64 Ni + 208 Pb reactions: Toward the production of unknown neutron-rich nuclei

In recent years, substantial efforts have been made for the study of multinucleon transfer reactions at energies around the Coulomb barrier both experimentally and theoretically, aiming at the production of unknown neutron-rich heavy nuclei. It is crucial to provide reliable theoretical predictions based on microscopic theories with sufficient predictive power. Purpose: This article aims to clarify the applicability of the quantal diffusion approach based on the stochastic mean-field (SMF) theory for multinucleon transfer processes. Isotope production cross sections are evaluated for the reactions of 64 Ni + 208 Pb at E c.m. = 268 MeV and 58 Ni + 208 Pb at E c.m. = 270 MeV and are compared with available experimental data. Methods: Three-dimensional time-dependent Hartree-Fock (TDHF) calculations are carried out for a range of initial orbital angular momenta with Skyrme SLy4d functional. Quantal diffusion equations, derived based on the SMF theory, for variances and covariance of neutron and proton numbers of reaction products are solved, with microscopic drift and diffusion coefficients obtained from time evolution of occupied single-particle orbitals in TDHF. Secondary de-excitation processes, both particle evaporation and fission, are simulated by a statistical compound-nucleus de-excitation model, gemini++. Dynamics of a fast isospin equilibration process followed by a slow drift toward the mass symmetry are commonly observed, as expected. Various reaction outcomes are evaluated, including average mass and charge numbers of reaction products, total kinetic energy loss (TKEL), scattering angle, contact time, and production cross sections for primary and secondary products. By comparing with the experimental data, we find that SMF and TDHF quantitatively reproduce experimental data for few-nucleon-transfer channels around the average values. In contrast, for many-nucleon-transfer channels, we find that the SMF approach provides much better description of the experimentally measured isotopic distributions. The results underline the importance of beyond-mean-field effects, especially one-body (mean-field) fluctuations and correlations, in describing multinucleon transfer processes. Moreover, through a combined analysis of SMF with a statistical model, gemini++, we find a significant contribution of transfer-induced fission, which is consistent with the experimental observation. In some cases, the SMF approach overestimates the isotopic width, requiring further improvements of the theoretical description. Possible ways to improve the description are discussed. The SMF approach is designed to describe the quantum many-body problem according to an ensemble of mean-field trajectories, taking into account part of many-body correlations in the description. As it requires feasible computational costs comparable to the ordinary TDHF approach, together with further model improvements, it will be a promising tool in the search for optimal reaction conditions to produce yet-unknown neutron-rich heavy nuclei through the multinucleon transfer reaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Hydration Mechanisms in Nanoparticle Interaction and Surface Energetics

Work on the project advanced molecular level understanding, prediction and control of nanoscale hydration in salt solutions under electric stimuli or ionic patterning. Permeation of nanoporous electrodes at preset voltage underlies the function of ultracapacitors. Transitory regulation of wetting in nanoporous media by electric field spans an array of applications in materials, energy storage, and separation sciences. The physically related modulation of nanoparticle solubility by surface charges can significantly extend the range of the nanomaterial applications, improve processing techniques, and potentially alleviate environmental concerns. Pore/solution equilibria, and activated kinetics of liquid gating by nanoconfined electrolytes are challenging problems at the forefront of experimental and theoretical research. Addressing these problems from a molecular perspective required the development of state-of-the-art simulation algorithms in statistical mechanics to capture complex processes in open systems under electric control. Parallel studies of wetting and dispersibility of polar and ionizing particles aim to uncover predictive relations between electrowetting, chemical functionalization, and geometry of nanomaterial particles. By nonequilibrium dynamic modeling, we elucidated electrolyte flow in nanochanels and associated electrokinetic energy conversion. Research on the program provided training opportunities for the next generation of scientists in computational chemistry. Insights, and methods from the project implicate broad segment of researchers in energy and nanosciences, materials and surface chemistry.

74 ATOMIC AND MOLECULAR PHYSICS↗

Interferometric binary phase estimations

We propose an interferometric scheme where each photon returns one bit of the binary expansion of an unknown phase. It sets up a method for estimating the phase value at arbitrary uncertainty. This strategy is global, since it requires no prior information, and it achieves the Heisenberg bound independently of the output statistics. We provide simulations and a characterization of this architecture.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The Information Length Concept Applied to Plasma Turbulence

A methodology to study statistical properties of anomalous transport in fusion plasma is investigated. Three time traces generated by the full-f gyrokinetic code GKNET are analyzed for this purpose. The time traces consist of heat flux as a function of the radial position, which is studied in a novel manner using statistical methods. The simulation data exhibit transport processes with both medium and long correlation length along the radius. A typical example of a phenomenon with long correlation length is avalanches. In order to investigate the evolution of the turbulent state, two basic configurations are studied, one flux-driven and one gradient-driven with decaying turbulence. The information length concept in tandem with Boltzmann–Gibbs and Tsallis entropy is used in the investigation. It is found that the dynamical states in both flux-driven and gradient-driven cases are surprisingly similar, but the Tsallis entropy reveals differences between them. This indicates that the types of probability distribution function are nevertheless quite different since the higher moments are significantly different.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Specific Effect of Peculiar Velocities on Dark-energy Constraints from Type Ia Supernovae

Peculiar velocities of Type Ia supernova (SN Ia) host galaxies affect the dark-energy parameter constraints in a small but very specific way: the parameters are biased in a single direction in parameter space that is a priori knowable for a given SN Ia data set. We demonstrate the latter fact with a combination of inference from a cosmological N-body simulation with overwhelming statistics applied to the Pantheon SN Ia data set, then confirm it by simple quantitative arguments. Here, we quantify small modifications to the current analyses that would ensure that the effect of cosmological parameters is essentially guaranteed to be negligible.

79 ASTRONOMY AND ASTROPHYSICS↗

The importance of Aitken mode aerosol particles for cloud sustenance in the summertime high Arctic – a simulation study supported by observational data

The potential importance of Aitken mode particles (diameters ~ 25–80 nm) for stratiform mixed-phase clouds in the summertime high Arctic (>80° N) has been investigated using two large-eddy simulation models. We find that, in both models, Aitken mode particles significantly affect the simulated microphysical and radiative properties of the cloud and can help sustain the cloud when accumulation mode concentrations are low (< 10–20 cm –3 ), even when the particles have low hygroscopicity (hygroscopicity parameter – κ=0.1). However, the influence of the Aitken mode decreases if the overall liquid water content of the cloud is low, either due to a higher ice fraction or due to low radiative cooling rates. An analysis of the simulated supersaturation (ss) statistics shows that the ss frequently reaches 0.5 % and sometimes even exceeds 1 %, which confirms that Aitken mode particles can be activated. The modelling results are in qualitative agreement with observations of the Hoppel minimum obtained from four different expeditions in the high Arctic. Our findings highlight the importance of better understanding Aitken mode particle formation, chemical properties and emissions, particularly in clean environments such as the high Arctic.

54 ENVIRONMENTAL SCIENCES↗

Reliability of cosmic dust data from Pioneers 8 and 9.

A comprehensive study of the characteristics and capabilities of the Pioneer cosmic dust mission is presented to facilitate accurate astronomical adaptations of the data by independent researchers. The characteristics of the sensor and associated electronics as they relate to dynamic range, field of view, and penetration effects in the film are discussed, comparisons are made between flight and simulation data, and statistical analyses are made of the reliability of the flight data. It is shown that the measurements from the Pioneer experiment are highly reliable and provide a valuable contribution to man's knowledge of the meteoroid environment of the solar system.

Gruen, E.↗

Thermodynamic analysis of turbulent combustion in a spark ignition engine. Experimental evidence

A method independent of physical modeling assumptions is presented to analyze high speed flame photography and cylinder pressure measurements from a transparent piston spark ignition research engine. The method involves defining characteristic quantities of the phenomena of flame propagation and combustion, and estimating their values from the experimental information. Using only the pressure information, the mass fraction curves are examined. An empirical burning law is presented which simulates such curves. Statistical data for the characteristics delay and burning angles which show that cycle to cycle fractional variations are of the same order of magnitude for both angles are discussed. The enflamed and burnt mass fractions are compared as are the rates of entrainment and burning.

Beretta, G. P.↗

Optimal mapping of irregular finite element domains to parallel processors

Mapping the solution domain of n-finite elements into N-subdomains that may be processed in parallel by N-processors is an optimal one if the subdomain decomposition results in a well-balanced workload distribution among the processors. The problem is discussed in the context of irregular finite element domains as an important aspect of the efficient utilization of the capabilities of emerging multiprocessor computers. Finding the optimal mapping is an intractable combinatorial optimization problem, for which a satisfactory approximate solution is obtained here by analogy to a method used in statistical mechanics for simulating the annealing process in solids. The simulated annealing analogy and algorithm are described, and numerical results are given for mapping an irregular two-dimensional finite element domain containing a singularity onto the Hypercube computer.

Flower, J.↗

LMSS propagation modeling at Virginia Tech

Recent efforts in the modeling of land mobile satellite systems are reported. These include descriptions of a simple model for prediction of fading statistics, a propagation simulator, and results from studies using the simulator. Predictions are compared to available measured data.

Stutzman, Warren L.↗