Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stochastic 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 1,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Computer simulations of two-pickup-ion instabilities in a cometary environment

This paper presents a one-dimensional electromagnetic hybrid computer simulation of distant cometary environments, in which protons and oxygen ions are injected continuously at the same time. It is demonstrated that, when, instead of a single component, both protons and oxygen ions are injected, the faster growing shorter-wavelength proton-proton right-hand resonant instability dominates the early time physics, leading to linear temporal growth of the fluctuating magnetic field energy density and pitch angle scattering of the protons (but causing no significant response by the nonresonant oxygen ions). At later times, as the injected oxygen ions accumulate sufficient density, the oxygen-proton right-hand resonant instability emerges, exhibiting exponential temporal growth to a relatively large peak amplitude, phase bunching, and magnetically trapping the oxygen ions. At times after peak amplitude, the oxygen-proton fluctuations stochastically accelerate some of the injected protons.

Gary, S. Peter↗
DOE OSTI2021

A deterministic verification strategy for electrostatic particle-in-cell algorithms in arbitrary spatial dimensions using the method of manufactured solutions

As simulations of kinetic plasmas continue to increase in scope and complexity, a rigorous and straightforward method for verifying particle-in-cell (PIC) implementations is necessary to ensure their correctness. Here, in this paper, we present a deterministic method for the rigorous verification of multidimensional, multispecies, electrostatic particle-in-cell codes based on the method of manufactured solutions. Specifically, we prove that rigorous verification is possible through the exclusive examination of errors of grid quantities (i.e., moments and/or fields), allowing for a very light-weight and non-intrusive implementation in existing PIC codes. This is a marked improvement over earlier PIC verification approaches (only demonstrated with one species in 1D-1V), which rely on the comparison of cumulative distribution functions, and are invasive on the code base, introduce additional stochastic noise, are significantly more computationally expensive, and lack rigorous convergence properties. Interestingly, we show that different grid quantities feature different rates of convergence with the number of particles and mesh size, impacting the verification process. These theoretical results are confirmed numerically with a multi-species 2D-2V particle-in-cell code, including the ability of the method to detect order reduction due to an incorrect implementation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2011

A Simple Postflight Measure of Postural Atania in Astronauts

Astronauts returning from space flight universally present with postural ataxia. Throughout the Space Shuttle Program, measurement of ataxia has concentrated on sway in the anterior-posterior (AP) plane. The current investigation, as a part of a larger functional study, concentrated on characterizing postural instability using dynamic stabilographic sway patterns in both the AP and medial-lateral (ML) planes. To accomplish this goal, six astronauts from short-duration (Shuttle) and three from long-duration (ISS) flights were required to recover from a simulated fall. Subjects with eyes open, wearing running shoes lay prone on the floor for 2 minutes and then quickly stood up, maintained a quiet stance for 3 minutes, arms relaxed along the side of the body, and feet comfortably placed on the force plate. Crewmembers were tested twice before flight, on landing day (Shuttle only), and 1, 6, and 30 days after flight. Anterior-posterior and ML center-of-pressure (COP) coordinates were calculated from the ground reaction forces collected at 500 Hz. The 3-minute quiet stance trial was broken into three 1-minute segments for stabilogram diffusion analysis. A mean sway speed (rate of change of COP displacement) was also calculated as an additional postural stability parameter. While there was considerable variation, most of crewmembers tested exhibited increased stochastic activity evidenced by larger short-term COP diffusion coefficients postflight in both the AP and ML planes, suggesting significant changes in postural control mechanisms, particularly control of lower limb muscle function. As expected, postural instability of ISS astronauts on the first day postflight was similar to that of Shuttle crewmembers on landing day. Recoveries of stochastic activity and mean sway speed to baseline levels were typically observed by the 30th day postflight for both long-duration and short-duration crewmembers. Dynamic postural stability characteristics obtained in this low-impact study complement the data measured with computerized dynamic posturography.

Reschke, M. F.↗
DOE OSTI2022

Binding kinetics of harmonically confined random walkers

Diffusion-mediated binding of molecules under the influence of discrete spatially confining potentials is a commonly encountered scenario in systems subjected to explicit fields or implicit fields arising from tethering restraints. Here, in this work, we derive analytical expressions for the mean binding time of two random walkers geometrically confined by means of two harmonic potentials in one- and two-dimensional systems, which show excellent agreement with Brownian dynamics simulations. As a demonstration of its utility, we use this theory to maximize the communication speed in existing DNA walkers, obtaining quantitative agreement with previously reported experimental findings. The analytical expressions derived in this paper are broadly applicable to diverse systems, providing ways to characterize communication processes and optimize the rate of signal propagation for sensing and computing applications at the nanoscale.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Novel Noise-Aware Classical Optimizer for Variational Quantum Algorithms

A key component of variational quantum algorithms (VQAs) is the choice of classical optimizer employed to update the parameterization of an ansatz. It is well recognized that quantum algorithms will, for the foreseeable future, necessarily be run on noisy devices with limited fidelities. Thus, the evaluation of an objective function (e.g., the guiding function in the quantum approximate optimization algorithm (QAOA) or the expectation of the electronic Hamiltonian in variational quantum eigensolver (VQE)) required by a classical optimizer is subject not only to stochastic error from estimating an expected value but also to error resulting from intermittent hardware noise. Model-based derivative-free optimization methods have emerged as popular choices of a classical optimizer in the noisy VQA setting, based on empirical studies. However, these optimization methods were not explicitly designed with the consideration of noise. In this work we adapt recent developments from the “noise-aware numerical optimization” literature to these commonly used derivative-free model-based methods. We introduce the key defining characteristics of these novel noise-aware derivative-free model-based methods that separate them from standard model-based methods. In conclusion, we study an implementation of such noise-aware derivative-free model-based methods and compare its performance on demonstrative VQA simulations to classical solvers packaged in scikit-quant.

classical optimizers↗
NASA NTRS1994

Design and evaluation of a Stochastic Optimal Feed-forward and Feedback Technology (SOFFT) flight control architecture

This paper describes the design and evaluation of a stochastic optimal feed-forward and feedback technology (SOFFT) control architecture with emphasis on the feed-forward controller design. The SOFFT approach allows the designer to independently design the feed-forward and feedback controllers to meet separate objectives and then integrate the two controllers. The feed-forward controller has been integrated with an existing high-angle-of-attack (high-alpha) feedback controller. The feed-forward controller includes a variable command model with parameters selected to satisfy level 1 flying qualities with a high-alpha adjustment to achieve desired agility guidelines, a nonlinear interpolation approach that scales entire matrices for approximation of the plant model, and equations for calculating feed-forward gains developed for perfect plant-model tracking. The SOFFT design was applied to a nonlinear batch simulation model of an F/A-18 aircraft modified for thrust vectoring. Simulation results show that agility guidelines are met and that the SOFFT controller filters undesired pilot-induced frequencies more effectively during a tracking task than a flight controller that has the same feedback control law but does not have the SOFFT feed-forward control.

Ostroff, Aaron J.↗
DOE OSTI2023

Probing fractional statistics in quantum simulators of spin liquid Hamiltonians

Recent advances in programmable quantum devices brought to the fore the intriguing possibility of using them to realise and investigate topological quantum spin liquid phases. This new and exciting direction brings about important research questions on how to probe and determine the presence of such exotic, highly entangled phases in a noisy quantum environment. One of the most promising tools is investigating the behaviour of the topological excitations, and in particular their fractional statistics. In this work we put forward a generic route to achieve this, and we illustrate it in the specific case of \mathbb{Z}_2 ℤ 2 topological spin liquids implemented with the aid of combinatorial gauge symmetry. We design a convenient architecture to study signatures of fractional statistics via quasiparticle interferometry, and we assess its robustness to diagonal and off-diagonal disorder, as well as to dephasing - effects that are generally pervasive in current quantum programmable devices. Interestingly, when turned on its head, our scheme provides a remarkably clear test of the “quantumness” of these devices, with robust signatures that crucially hinge on quantum coherence and quantum interference effects, and cannot be mimicked by classical stochastic processes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Effect of void positioning on the detonation sensitivity of a heterogeneous energetic material

We show although it is well-established that voids profoundly influence the initiation and reaction behaviors of heterogeneous energetic materials such as polymer-bonded explosives (PBX) and propellants, there has been little study of how void location in different constituents in the microstructures of such materials affect the macroscale behavior. Here, we use three-dimensional (3D) mesoscale simulations to study how void placement within the reactive grains versus the polymer binder influences the shock-to-detonation transition (SDT) in a polymer-bonded explosive. The material studied here has a microstructure comprised of 75% PETN (pentaerythritol tetranitrate) grains and 25% HTPB (hydroxyl-terminated polybutadiene) polymer binder by volume. Porosities up to 10% in the form of spherical voids distributed in both the grains and polymer are considered. An Arrhenius reactive burn relation is used to model the chemical kinetics of the PETN grains under shock loading, thereby resolving the heterogeneous detonation behavior of the PBX. The influence of void location on the shock initiation sensitivity of the material is quantitatively ranked by comparing the predicted run distance to detonation (RDD) for each sample. The analysis includes inherent quantification of uncertainties arising from the stochastic variations in the microstructure morphologies and void distributions by using statistically equivalent microstructure sample sets (SEMSS), leading to probabilistic formulations for the RDD as a function of shock pressure. The calculations reveal that the location of voids in the composite microstructure significantly affects the RDD. Specifically, voids exclusively within the grains cause the PBX to be more sensitive (having shorter RDD) than voids in the polymer binder. Unique probabilistic relationships are derived to map the probability of observing RDD for each void location material case, allowing for prediction of initiation behavior anywhere in the shock pressure – RDD space. These findings agree with trends reported in the literature.

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

Confinement time and ambipolar potential in a relativistic mirror-confined plasma

Advanced aneutronic fusion fuels such as proton-Boron 11 tend to require much higher temperatures than conventional fuels like deuterium–tritium. For electrons, the bulk plasma temperature can approach a substantial fraction of the rest mass. In a mirror confinement system, where the electrons are confined by an ambipolar potential of at least five electron temperatures, the tail electrons which can escape the potential are fully relativistic, which must be taken into account in calculating their confinement. In this paper, simple estimates are employed to extend the scaling of the confinement time into the relativistic regime. By asymptotically matching this scaling to known solutions in the non-relativistic limit, accurate forms for the confinement time (and thus, the ambipolar potential) are obtained. Furthermore, these forms are verified using finite-element-based Fokker–Planck simulations over a wide range of parameters. Comparing relativistic and nonrelativistic mirror-confined plasmas with the same ratio of confining potential |eΦ| to electron temperature T e and the same mirror ratio R, the net result is a decrease in the confinement time due to relativistic effects by a factor of S ≡ (1 + 15T e /8m e c 2 )/(1 + 2|eΦ|/m e c 2 )⁠.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2000

A Probabilistic Framework for the Validation and Certification of Computer Simulations

The paper presents a methodology for quantifying, propagating, and managing the uncertainty in the data required to initialize computer simulations of complex phenomena. The purpose of the methodology is to permit the quantitative assessment of a certification level to be associated with the predictions from the simulations, as well as the design of a data acquisition strategy to achieve a target level of certification. The value of a methodology that can address the above issues is obvious, specially in light of the trend in the availability of computational resources, as well as the trend in sensor technology. These two trends make it possible to probe physical phenomena both with physical sensors, as well as with complex models, at previously inconceivable levels. With these new abilities arises the need to develop the knowledge to integrate the information from sensors and computer simulations. This is achieved in the present work by tracing both activities back to a level of abstraction that highlights their commonalities, thus allowing them to be manipulated in a mathematically consistent fashion. In particular, the mathematical theory underlying computer simulations has long been associated with partial differential equations and functional analysis concepts such as Hilbert spares and orthogonal projections. By relying on a probabilistic framework for the modeling of data, a Hilbert space framework emerges that permits the modeling of coefficients in the governing equations as random variables, or equivalently, as elements in a Hilbert space. This permits the development of an approximation theory for probabilistic problems that parallels that of deterministic approximation theory. According to this formalism, the solution of the problem is identified by its projection on a basis in the Hilbert space of random variables, as opposed to more traditional techniques where the solution is approximated by its first or second-order statistics. The present representation, in addition to capturing significantly more information than the traditional approach, facilitates the linkage between different interacting stochastic systems as is typically observed in real-life situations.

Ghanem, Roger↗
NASA NTRS2010

Postflight Quiet Stance Stability of Astronauts Following Recovery From a Simulated Fall

INTRODUCTION: Astronauts returning from space flight universally present with postural ataxia. Throughout the Space Shuttle Program, measurement of ataxia has concentrated on sway in the anterior-posterior plane. Implementation of an interdisciplinary pre- and postflight study (Functional Task Test, FTT) designed to evaluate both astronaut postflight functional performance and related physiological changes has allowed the investigation of postural instability by characterizing dynamic stabilographic sway patterns. METHODS: Six astronauts from short-duration (Shuttle) and three from long-duration (ISS) flights were required to recover from a simulated fall. Subjects with eyes open, wearing running shoes lay prone on the floor for 2 minutes and then quickly stood up, maintained a quiet stance for 3 minutes, arms relaxed along the side of the body, and feet comfortably placed on the force plate. Crewmembers were tested twice before flight, on landing day (Shuttle only), and 1, 6, and 30 days after flight. Anterior-posterior (AP) and medial-lateral (ML) center-of-pressure (COP) coordinates were calculated from the ground reaction forces collected at 500 Hz. The 3-minute quiet stance trial was broken into three 1-minute segments for stabilogram diffusion analysis. A mean sway speed (rate of change of COP displacement) was also calculated as an additional postural stability parameter. RESULTS/CONCLUSION: While there was considerable variation, most of crewmembers tested exhibited increased stochastic activity evidenced by larger short-term COP diffusion coefficients postflight in both the AP and ML planes, suggesting significant changes in postural control mechanisms, particularly control of lower limb muscle function. As expected, postural instability of ISS astronauts on the first day postflight was similar to that of Shuttle crewmembers on landing day. Recoveries of stochastic activity and mean sway speed to baseline levels were typically observed by the 30th day postflight for both long-duration and short-duration crewmembers. Dynamic postural stability characteristics obtained in this low-impact study complement the data measured with computerized dynamic posturography.

Reschke, M. F.↗
NASA NTRSDate not supplied

A Practical Approach to Uncertainty Quantification Using Probability Boxes

To date, while the use of CFD for aerospace vehicle design and development is prevalent, the documentation of uncertainties associated with the simulations are rare. Instead, the current state-of-the-art relies heavily on the experience of the CFD practitioner to estimate the uncertainty associated with their simulations through simple sensitivity studies or subject matter expertise. This practice will have to be replaced with a formal uncertainty quantification (UQ) process if CFD is to play an expanded role in the research and engineering design community, test and evaluation community, and ultimately certification for flight. Accounting for uncertainties in a formal manner is a tedious process. Moreover, the typical CFD practitioner is not likely to be familiar with formal UQ methods. These factors have prevented the adoption of UQ methods in the engineering design and development cycle. This presentation will outline a credible approach to UQ using Probability Boxes that is straightforward to apply, and can readily be automated using existing UQ tool sets such as the DAKOTA packaged developed at Sandia. The added expense incurred when moving away from a deterministic CFD process to a stochastic one that captures uncertainties to enable risk-informed decision making will be discussed, as well as effective ways to reduce the computational costs.

Uncertainty Quantification↗
NASA NTRS2017

Development of a Detailed Surface Chemistry Framework in DSMC

Many of the current direct simulation Monte Carlo (DSMC) codes still employ only simple surface catalysis models. These include only basic mechanisms such as dissociation, recombination, and exchange reactions, without any provision for adsorption and finite rate kinetics. Incorporating finite rate chemistry at the surface is increasingly becoming a necessity for various applications such as high speed re-entry flows over thermal protection systems (TPS), micro-electro-mechanical systems (MEMS), surface catalysis, etc. In the recent years, relatively few works have examined finite-rate surface reaction modeling using the DSMC method.In this work, a generalized finite-rate surface chemistry framework incorporating a comprehensive list of reaction mechanisms is developed and implemented into the DSMC solver SPARTA. The various mechanisms include adsorption, desorption, Langmuir-Hinshelwood (LH), Eley-Rideal (ER), Collision Induced (CI), condensation, sublimation, etc. The approach is to stochastically model the various competing reactions occurring on a set of active sites. Both gas-surface (e.g., ER, CI) and pure-surface (e.g., LH, desorption) reaction mechanisms are incorporated. The reaction mechanisms could also be catalytic or surface altering based on the participation of the bulk-phase species (e.g., bulk carbon atoms). Marschall and MacLean developed a general formulation in which multiple phases and surface sites are used and we adopt a similar convention in the current work. Microscopic parameters of reaction probabilities (for gas-surface reactions) and frequencies (for pure-surface reactions) that are require for DSMC are computed from the surface properties and macroscopic parameters such as rate constants, sticking coefficients, etc. The energy and angular distributions of the products are decided based on the reaction type and input parameters. Thus, the user has the capability to model various surface reactions via user-specified reaction rate constants, surface properties and parameters.

Gas↗
DOE OSTI2022

Model-Free Measurement of Local Entropy Production and Extractable Work in Active Matter

Time-reversal symmetry breaking and entropy production are universal features of nonequilibrium phenomena. Despite its importance in the physics of active and living systems, the entropy production of systems with many degrees of freedom has remained of little practical significance because the high dimensionality of their state space makes it difficult to measure. Here, in this work, we introduce a local measure of entropy production and a numerical protocol to estimate it. We establish a connection between the entropy production and extractability of work in a given region of the system and show how this quantity depends crucially on the degrees of freedom being tracked. We validate our approach in theory, simulation, and experiments by considering systems of active Brownian particles undergoing motility-induced phase separation, as well as active Brownian particles and E.coli in a rectifying device in which the time-reversal asymmetry of the particle dynamics couples to spatial asymmetry to reveal its effects on a macroscopic scale.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Distributed Stochastic Optimization of a Neural Representation Network for Time-Space Tomography Reconstruction

4D time-space reconstruction of dynamic events or deforming objects using X-ray computed tomography (CT) is an important inverse problem in non-destructive evaluation. Conventional back-projection based reconstruction methods assume that the object remains static for the duration of several tens or hundreds of X-ray projection measurement images (reconstruction of consecutive limited-angle CT scans). However, this is an unrealistic assumption for many in-situ experiments that causes spurious artifacts and inaccurate morphological reconstructions of the object. To solve this problem, we propose to perform a 4D time-space reconstruction using a distributed implicit neural representation (DINR) network that is trained using a novel distributed stochastic training algorithm. Our DINR network learns to reconstruct the object at its output by iterative optimization of its network parameters such that the measured projection images best match the output of the CT forward measurement model. Here, we use a forward measurement model that is a function of the DINR outputs at a sparsely sampled set of continuous valued 4D object coordinates. Unlike previous neural representation architectures that forward and back propagate through dense voxel grids that sample the object's entire time-space coordinates, we only propagate through the DINR at a small subset of object coordinates in each iteration resulting in an order-of-magnitude reduction in memory and compute for training. DINR leverages distributed computation across several compute nodes and GPUs to produce high-fidelity 4D time-space reconstructions. We use both simulated parallel-beam and experimental cone-beam X-ray CT datasets to demonstrate the superior performance of our approach.

36 MATERIALS SCIENCE↗
DOE OSTI2021

ELECTRON COOLING USING A PULSED AND DITHERING BEAM FROM AN ELECTROSTATIC ELECTRON COOLER

Electron cooling continues to be an invaluable technique to reduce and maintain the emittance in hadron storage rings such as the EIC and EICC where stochastic cooling is inefficient and radiative cooling is negligible. Extending the energy range of electron coolers beyond what is feasible with the conventional, electrostatic approach necessitates the use of RF fields for acceleration and, thus, a bunched electron beam. To experimentally investigate how the relative time structure of the two beams affects the cooling properties, we have set up a pulsed-beam cooling device by adding a synchronized pulsing circuit to the conventional electron source of the CSRm cooler at Institute of Modern Physics. Using both constant and modulated synchronization between electron pulses and ion bunches, we have measured the effects of the electron bunch length and longitudinal ion focusing strength on the temporal evolution of the longitudinal and transverse ion beam profile and demonstrated the detrimental effect of timing jitter as predicted by space-charge theory and simulations. Our experiment suggests a need for further investigations of specific aspects of bunched cooling such as synchro-betatron coupling and phase dithering effects when using short electron bunches to cool longer ion bunches. However, given the comparatively long IBS lifetime of higher-energy proton storage rings like the EIC, slow dithering could potentially provide an option to save cost on the electron cooler linac.

Wang, Haipeng↗
NASA NTRS1986

Two-phase flow

An experimental program to characterize the spray from candidate nozzles for icing-cloud simulation is discussed. One canidate nozzle, which is currently used for icing research, has been characterized for flow and drop size. The median-volume diameter (MVD) from this air-assist nozzle is compared with correlations in the literature. The new experimental spray facility is discussed, and the drop-size instruments are discussed in detail. Since there is no absolute standard for drop-size measurements and there are other limitations, such as drop -size range and velocity range, several instruments are used and results are compared. A two-phase model was developed at Pennsylvania State University. The model uses the k-epsilon model of turbulence in the continous phase. Three methods for treating the discrete phase are used: (1) a locally homogeneous flow (LHF) model, (2) a deterministic separated flow (DSF) model, and (3) a stochastic separated flow (SSF) model. In the LHF model both phases have the same velocity and temperature at each point. The DSF model provides interphase transport but ignores the effects of turbulent fluctuations. In the SSF model the drops interact with turbulent eddies whose properties are determined by the k-epsilon turbulence model. The two-phase flow model has been extended to include the effects of evaporation and combustion.

Tacina, Robert R.↗
DOE OSTI2023

Nonlocal suppression of Biermann battery magnetic-field generation for arbitrary atomic numbers and magnetization

The Biermann battery term of magnetohydrodynamics (MHD) generates a magnetic field where electron density gradients and electron temperature gradients are perpendicular to one another. Kinetic simulations and experiments have shown that the rate of magnetic-field generation is lower than Biermann when the electron mean free path becomes comparable to or greater than the temperature gradient scale length, known as the nonlocal regime. We investigate the nonlocal suppression of the Biermann term using simplified Fokker–Planck simulations covering a wide range of parameters. We provide the first fit for nonlocal Biermann suppression that has physically accurate behavior for small and large values of a suitable nonlocality parameter, valid for an arbitrary atomic number, and that includes the effect of magnetization on nonlocality. Lastly, the fit is intended to provide an approximate method to account for reduced magnetic-field generation in MHD codes and theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗