Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stochastic approximation”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

S-OPT: A Points Selection Algorithm for Hyper-Reduction in Reduced Order Models

While projection-based reduced order models can reduce the dimension of full order solutions, the resulting reduced models may still contain terms that scale with the full order dimension. Hyper-reduction techniques are sampling-based methods that further reduce this computational complexity by approximating such terms with a much smaller dimension. The goal of this work is to introduce the points selection algorithm developed by Shin and Xiu as a hyper-reduction method. The selection algorithm was originally proposed as a stochastic collocation method for uncertainty quantification. Since the algorithm aims at maximizing a quantity $\mathcal{S}$ that measures both the column orthogonality and the determinant, we refer to the algorithm as S-OPT. Numerical examples are provided to demonstrate the performance of S-OPT and to compare its performance with a gappy proper orthogonal decomposition (POD) algorithm. Here, we found that using the S-OPT algorithm is shown to predict the full order solutions with higher accuracy than gappy POD especially when the number of sampling points is small, although we note that S-OPT shows slow asymptotic convergence with respect to the number of samples for some applications, e.g., Lagrangian hydrodynamics.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2014

Simulations of DSB Yields and Radiation-induced Chromosomal Aberrations in Human Cells Based on the Stochastic Track Structure Induced by HZE Particles

The formation of double‐strand breaks (DSBs) and chromosomal aberrations (CAs) is of great importance in radiation research and, specifically, in space applications. We are presenting a new particle track and DNA damage model, in which the particle stochastic track structure is combined with the random walk (RW) structure of chromosomes in a cell nucleus. The motivation for this effort stems from the fact that the model with the RW chromosomes, NASARTI (NASA radiation track image) previously relied on amorphous track structure, while the stochastic track structure model RITRACKS (Relativistic Ion Tracks) was focused on more microscopic targets than the entire genome. We have combined chromosomes simulated by RWs with stochastic track structure, which uses nanoscopic dose calculations performed with the Monte‐Carlo simulation by RITRACKS in a voxelized space. The new simulations produce the number of DSBs as function of dose and particle fluence for high‐energy particles, including iron, carbon and protons, using voxels of 20 nm dimension. The combined model also calculates yields of radiation‐induced CAs and unrejoined chromosome breaks in normal and repair deficient cells. The joined computational model is calibrated using the relative frequencies and distributions of chromosomal aberrations reported in the literature. The model considers fractionated deposition of energy to approximate dose rates of the space flight environment. The joined model also predicts of the yields and sizes of translocations, dicentrics, rings, and more complex‐type aberrations formed in the G0/G1 cell cycle phase during the first cell division after irradiation. We found that the main advantage of the joined model is our ability to simulate small doses: 0.05‐0.5 Gy. At such low doses, the stochastic track structure proved to be indispensable, as the action of individual delta‐rays becomes more important.

Ponomarev, Artem↗
NASA NTRS2014

Simulations of DSB Yields and Radiation-induced Chromosomal Aberrations in Human Cells Based on the Stochastic Track Structure iIduced by HZE Particles

The formation of double‐strand breaks (DSBs) and chromosomal aberrations (CAs) is of great importance in radiation research and, specifically, in space applications. We are presenting a new particle track and DNA damage model, in which the particle stochastic track structure is combined with the random walk (RW) structure of chromosomes in a cell nucleus. The motivation for this effort stems from the fact that the model with the RW chromosomes, NASARTI (NASA radiation track image) previously relied on amorphous track structure, while the stochastic track structure model RITRACKS (Relativistic Ion Tracks) was focused on more microscopic targets than the entire genome. We have combined chromosomes simulated by RWs with stochastic track structure, which uses nanoscopic dose calculations performed with the Monte‐Carlo simulation by RITRACKS in a voxelized space. The new simulations produce the number of DSBs as function of dose and particle fluence for high‐energy particles, including iron, carbon and protons, using voxels of 20 nm dimension. The combined model also calculates yields of radiation‐induced CAs and unrejoined chromosome breaks in normal and repair deficient cells. The joined computational model is calibrated using the relative frequencies and distributions of chromosomal aberrations reported in the literature. The model considers fractionated deposition of energy to approximate dose rates of the space flight environment. The joined model also predicts of the yields and sizes of translocations, dicentrics, rings, and more complex‐type aberrations formed in the G0/G1 cell cycle phase during the first cell division after irradiation. We found that the main advantage of the joined model is our ability to simulate small doses: 0.05‐0.5 Gy. At such low doses, the stochastic track structure proved to be indispensable, as the action of individual delta‐rays becomes more important.

Ponomarev, Artem↗
DOE OSTI2021

High-throughput bend-strengths of ultra-small polysilicon MEMS components

The strength distribution of polysilicon bend specimens, approximately 10 μm in size, is measured using a high-throughput microelectromechanical system fabrication and testing method. The distribution is predicted from reference tests on tensile specimens and finite element analysis of the bend specimen geometry incorporated into a stochastic extreme-value strength framework. Agreement between experiment and prediction suggests that the ultra-small specimens may be at the limit of extreme-value scaling and contain only one strength-controlling flaw/specimen.

Cook, Robert F. (ORCID:0000000304228881)↗
DOE OSTI2021

Optimization models for integrated biorefinery operations

Variations of physical and chemical characteristics of biomass lead to an uneven flow of biomass in a biorefinery, which reduces equipment utilization and increases operational costs. Uncertainty of biomass supply and high processing costs increase the risk of investing in the US’s cellulosic biofuel industry. We propose a stochastic programming model to streamline processes within a biorefinery. A chance constraint models system’s reliability requirement that the reactor is operating at a high utilization rate given uncertain biomass moisture content, particle size distribution, and equipment failure. The model identifies operating conditions of equipment and inventory level to maintain a continuous flow of biomass to the reactor. Furthermore, the sample average approximation method approximates the chance constraint and a bisection search-based heuristic solves this approximation. A case study is developed using real-life data collected at Idaho National Laboratory’s biomass processing facility. An extensive computational analysis indicates that sequencing of biomass bales based on moisture level, increasing storage capacity, and managing particle size distribution, increases utilization of the reactor and reduces operational costs.

09 BIOMASS FUELS↗
NASA NTRS1990

Time development of a perturbed-spherical nucleus in a pure supercooled liquid. I - Power-law growth of morphological instabilities. II - Nonlinear development

The linear growth stage of the morphological instabilities developing at the liquid-solid interface during crystal growth from an axisymmetric spherical nucleus is analyzed. The corresponding growth rate parameters are calculated numerically, and it is shown that morphological instabilities for free growth evolve according to a power law in agreement with WKB results and contrary to an exponential law found in a quasi-stationary approximation. The time evolution of an arbitrary perturbed liquid-solid interface from a linear into the nonlinear stage is studied. The initial perturbations include the eigenfunctions for the linear problem, localized perturbations of the grain-boundary type, and a stochastic noise. It is shown that the perturbations grow and spread in a wavelike manner. The formation of a predendritic growth stage is characterized by establishment of constant values of tip radius and velocity.

Pines, V.↗
NASA NTRS1994

Optimal active control for Burgers equations

A method for active fluid flow control based on control theory is discussed. Dynamic programming and fixed point successive approximations are used to accommodate the nonlinear control problem. The long-term goal of this project is to establish an effective method applicable to complex flows such as turbulence and jets. However, in this report, the method is applied to stochastic Burgers equation as an intermediate step towards this goal. Numerical results are compared with those obtained by gradient search methods.

Ikeda, Yutaka↗
NASA NTRS2006

What does reflection from cloud sides tell us about vertical distribution of cloud droplets?

In order to accurately measure the interaction of clouds with aerosols, we have to resolve the vertical distribution of cloud droplet sizes and determine the temperature of glaciation for clean and polluted clouds. Knowledge of the droplet vertical profile is also essential for understanding precipitation. So far, all existing satellites either measure cloud microphysics only at cloud top (e.g., MODIS) or give a vertical profile of precipitation sized droplets (e.g., Cloudsat). What if one measures cloud microphysical properties in the vertical by retrieving them from the solar and infrared radiation reflected or emitted from cloud sides? This was the idea behind CLAIM-3D (A 3D - cloud aerosol interaction mission) recently proposed by NASA GSFC. This presentation will focus on the interpretation of the radiation reflected from cloud sides. In contrast to plane-parallel approximation, a conventional approach to all current operational retrievals, 3D radiative transfer will be used for interpreting the observed reflectances. As a proof of concept, we will show a few examples of radiation reflected from cloud fields generated by a simple stochastic cloud model with prescribed microphysics. Instead of fixed values of the retrieved effective radii, the probability density functions of droplet size distributions will serve as possible retrievals.

Marshak, A.↗
DOE OSTI2023

Real-Time Krylov Theory for Quantum Computing Algorithms

Quantum computers provide new avenues to access ground and excited state properties of systems otherwise difficult to simulate on classical hardware. New approaches using subspaces generated by real-time evolution have shown efficiency in extracting eigenstate information, but the full capabilities of such approaches are still not understood. In recent work, we developed the variational quantum phase estimation (VQPE) method, a compact and efficient real-time algorithm to extract eigenvalues on quantum hardware. Here we build on that work by theoretically and numerically exploring a generalized Krylov scheme where the Krylov subspace is constructed through a parametrized real-time evolution, which applies to the VQPE algorithm as well as others. We establish an error bound that justifies the fast convergence of our spectral approximation. We also derive how the overlap with high energy eigenstates becomes suppressed from real-time subspace diagonalization and we visualize the process that shows the signature phase cancellations at specific eigenenergies. We investigate various algorithm implementations and consider performance when stochasticity is added to the target Hamiltonian in the form of spectral statistics. To demonstrate the practicality of such real-time evolution, we discuss its application to fundamental problems in quantum computation such as electronic structure predictions for strongly correlated systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Solving the $k$-Sparse Eigenvalue Problem with Reinforcement Learning

We examine the possibility of using a reinforcement learning (RL) algorithm to solve large-scale eigenvalue problems in which the desired the eigenvector can be approximated by a sparse vector with at most k nonzero elements, where k is relatively small compare to the dimension of the matrix to be partially diagonalized. Here, this type of problem arises in applications in which the desired eigenvector exhibits localization properties and in large-scale eigenvalue computations in which the amount of computational resource is limited. When the positions of these nonzero elements can be determined, we can obtain the k-sparse approximation to the original problem by computing eigenvalues of a k × k submatrix extracted from k rows and columns of the original matrix. We review a previously developed greedy algorithm for incrementally probing the positions of the nonzero elements in a k-sparse approximate eigenvector and show that the greedy algorithm can be improved by using an RL method to refine the selection of k rows and columns of the original matrix. We describe how to represent states, actions, rewards and policies in an RL algorithm designed to solve the k-sparse eigenvalue problem and demonstrate the effectiveness of the RL algorithm on two examples originating from quantum many-body physics.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1986

Heat and mass transfer in flames

Heat- and mass-transfer processes in turbulent diffusion flames are discussed, considering turbulent mixing and the structure of single-phase flames, drop processes in spray flames, and nonluminous and luminous flame radiation. Interactions between turbulence and other phenomena are emphasized, concentrating on past work of the author and his associates. The conserved-scalar formalism, along with the laminar-flamelet approximation, is shown to provide reasonable estimates of the structure of gas flames, with modest levels of empiricism. Extending this approach to spray flames has highlighted the importance of drop/turbulence interactions; e.g., turbulent dispersion of drops, modification of turbulence by drops, etc. Stochastic methods being developed to treat these phenomena are yielding encouraging results.

Faeth, G. M.↗
NASA NTRS1991

Structural factoring approach for analyzing stochastic networks

The problem of finding the distribution of the shortest path length through a stochastic network is investigated. A general algorithm for determining the exact distribution of the shortest path length is developed based on the concept of conditional factoring, in which a directed, stochastic network is decomposed into an equivalent set of smaller, generally less complex subnetworks. Several network constructs are identified and exploited to reduce significantly the computational effort required to solve a network problem relative to complete enumeration. This algorithm can be applied to two important classes of stochastic path problems: determining the critical path distribution for acyclic networks and the exact two-terminal reliability for probabilistic networks. Computational experience with the algorithm was encouraging and allowed the exact solution of networks that have been previously analyzed only by approximation techniques.

Hayhurst, Kelly J.↗
DOE OSTI2020

LANL studies of ESD spark energy discharge through September 2019

Two concentric-sphere apparatuses were used to measure the energies deposited by electric sparks through air to calibrated loads. Both sets of spheres had clean, uninsulated metal surfaces. The 127-mm diameter pair of spheres had a gap of 10 mm and a capacitance of 52 pF. The maximum energy stored on these spheres was 12 mJ and the maximum energy delivered to the 50-mΩ load was 11 µJ, with or without inductance. This load energy amounted to 0.09% of the available input energy; the remainder of the input energy was expended in forming the spark. The second set of spheres was 330 mm in diameter with a gap of 14 mm and a capacitance of 240 pF. The maximum energy stored on these spheres was 71 mJ and the maximum energy delivered to the 50-mΩ load was 150 µJ without added inductance. This load energy was 0.27% of the input energy; again, the remainder of the input energy was expended in forming the spark. The inclusion of 250 nH and 500 nH inductances in the discharge path of the smaller spheres caused the current discharges to resonate, but had little effect on the maximum energies deposited in the loads. At the time of writing, inductances had not yet been added to the larger spheres. The magnitudes of the observed spark discharges were stochastic, so more than 1200 tests were required to develop the statistical bounds of the spark behavior, i.e., to estimate the worst case (highest likely) energy depositions in the loads. The action integrals of the current discharges $∫ i^2dt$ were found to be approximately independent of the load resistances for loads of less than ~1 Ω. Hence the energy deposition scaled with resistance as 207 µJ/Ω and 2.93 mJ/Ω for the two apparatuses.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS1991

Probabilistic Structural Analysis Methods (PSAM) for select space propulsion system components

The fourth year of technical developments on the Numerical Evaluation of Stochastic Structures Under Stress (NESSUS) system for Probabilistic Structural Analysis Methods is summarized. The effort focused on the continued expansion of the Probabilistic Finite Element Method (PFEM) code, the implementation of the Probabilistic Boundary Element Method (PBEM), and the implementation of the Probabilistic Approximate Methods (PAppM) code. The principal focus for the PFEM code is the addition of a multilevel structural dynamics capability. The strategy includes probabilistic loads, treatment of material, geometry uncertainty, and full probabilistic variables. Enhancements are included for the Fast Probability Integration (FPI) algorithms and the addition of Monte Carlo simulation as an alternate. Work on the expert system and boundary element developments continues. The enhanced capability in the computer codes is validated by applications to a turbine blade and to an oxidizer duct.

Source record↗
DOE OSTI2026

Near-Optimal Performance of Stochastic Model Predictive Control

Here, this article presents a regret analysis for stochastic model predictive control (SMPC) in linear systems with quadratic performance index and additive and multiplicative uncertainties. Under a finite support assumption, the problem can be cast as a finite-dimensional quadratic program, but the problem becomes quickly intractable as the problem size grows exponentially in the horizon length. SMPC aims to compute approximate solutions by solving a sequence of problems with truncated prediction horizons and committing the solution in a receding-horizon fashion. Although this approach is widely used in practice, its performance relative to the optimal solution is not well understood. This article reports for the first time a rigorous near-optimal performance guarantee of SMPC: under stabilizability and detectability conditions, the regret of SMPC is exponentially small in the prediction horizon length, allowing SMPC to achieve near-optimal performance at a substantially reduced computational expense.

93E20, 93B45↗
NASA NTRS1986

Unimolecular diffusion-mediated reactions with a nonrandom time-modulated absorbing barrier

A diffusion-reaction model with time-dependent reactivity is formulated and applied to unimolecular reactions. The model is solved exactly numerically and approximately analytically for the unreacted fraction as a function of time. It is shown that the approximate analytical solution is valid even when the system is far from equilibrium, and when the reactivity probability is more complicated than a square-wave function of time. A discussion is also given of an approach to problems of this type using a stochastically fluctuating reactivity, and the first-passage time for a particular example is derived.

Bashford, D.↗
NASA NTRS1991

Probabilistic structural analysis algorithm development for computational efficiency

The PSAM (Probabilistic Structural Analysis Methods) program is developing a probabilistic structural risk assessment capability for the SSME components. An advanced probabilistic structural analysis software system, NESSUS (Numerical Evaluation of Stochastic Structures Under Stress), is being developed as part of the PSAM effort to accurately simulate stochastic structures operating under severe random loading conditions. One of the challenges in developing the NESSUS system is the development of the probabilistic algorithms that provide both efficiency and accuracy. The main probability algorithms developed and implemented in the NESSUS system are efficient, but approximate in nature. In the last six years, the algorithms have improved very significantly.

Wu, Y.-T.↗
NASA NTRS1989

Flutter suppression control law synthesis for the Active Flexible Wing model

The Active Flexible Wing Project is a collaborative effort between the NASA Langley Research Center and Rockwell International. The objectives are the validation of methodologies associated with mathematical modeling, flutter suppression control law development and digital implementation of the control system for application to flexible aircraft. A flutter suppression control law synthesis for this project is described. The state-space mathematical model used for the synthesis included ten flexible modes, four control surface modes and rational function approximation of the doublet-lattice unsteady aerodynamics. The design steps involved developing the full-order optimal control laws, reducing the order of the control law, and optimizing the reduced-order control law in both the continuous and the discrete domains to minimize stochastic response. System robustness was improved using singular value constraints. An 8th order robust control law was designed to increase the symmetric flutter dynamic pressure by 100 percent. Preliminary results are provided and experiences gained are discussed.

Mukhopadhyay, Vivek↗