SEARCH · Engineering Papers

Results for “STOCHASTIC PROCESSES”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 973 records · Page 54

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Kink mechanism in Cu/Nb nanolaminates explored by $\mathcal{in}$ $\mathcal{situ}$ pillar compression

We report Nano metallic laminates (NMLs) exhibit different failure modes depending on the loading conditions due to their mechanical anisotropies. Kinking is a typical failure mode in many NMLs compressed along a layer-parallel direction. However, a detailed description of the microstructure evolution during kink band (KB) formation and an in-depth understanding of the formation mechanisms are lacking. In this work, the KB process is investigated in Cu/Nb NMLs by in situ micro pillar compression in the scanning electron microscope (SEM) along a layer-parallel direction. Post-mortem S/TEM and transmission Kikuchi diffraction (TKD) analyses show that kink banding leads to significant microstructure changes characterized by an accumulation of geometrically necessary dislocations (GNDs) and of tilt geometrically necessary boundaries (GNBs) near KB boundaries (KBBs). The distinct microstructure evolution implies that KB formation is facilitated by the inhomogeneous microstructures resulting in constrained deformation modes. Specifically, dislocations active on slip planes nearly parallel to the interfaces make a major contribution to kink evolution after the onset of kinking. Once layer-parallel slip systems are activated, preexisting lattice dislocations and dislocations nucleating from interfaces will accumulate as GNDs near KBBs via the stochastic storage of lattice dislocations that have certain Burgers vectors. GNDs can further transform into GNBs via cross-slip and climb driven processes near the KBB. Furthermore, GNBs near KBBs can grow by incorporating more GNDs or by coalescence to accommodate the KB evolution. We further hypothesize that microstructural perturbations and their ensuing stresses can initiate KB formation in Cu/Nb NMLs.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Reliability-informed end-of-use decision making for product sustainability using two-stage stochastic optimization

The concept of circular economy has been diffused in recent decades to promote economic growth that does not add to the burden on natural resource extraction. Re-X options (e.g., reuse, repair, refurbish, remanufacture, recycle) have been gradually adopted in the product development process and optimized to reduce or eliminate waste and pollution. Although manufacturing incorporating Re-X options can be more environmentally friendly, it involves more sources of uncertainty than traditional manufacturing since the end-of-use products can be collected from multiple origins with various quantities and qualities, and the market demand for both new and remanufactured products cannot be forecasted perfectly. Thus, there is a need to optimize the Re-X policy to alleviate the negative impacts of the higher uncertainty. One option is using the reliability information of new products to estimate the end-of-use conditions and applying multi-stage stochastic optimization to capture multiple demand scenarios. This paper develops a two-stage stochastic optimization model to optimize the quality thresholds for reuse, recycling, and remanufacturing options. Our objective is to minimize the total cost, energy consumption, and environmental impact of producing and providing warranty service for a product family. The model employs reliability information of product components to estimate the warranty service cost and the end-of-use conditions of the returned resources. A case study on a general product family is implemented to illustrate the efficacy of the optimization model. Finally, results show that the two-stage optimization can achieve cost and environmental impact reduction for a hybrid manufacturing and remanufacturing process.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1988

Creep life prediction based on stochastic model of microstructurally short crack growth

A nondimensional model of microstructurally short crack growth in creep is developed based on a detailed observation of the creep fracture process of 304 stainless steel. In order to deal with the scatter of small crack growth rate data caused by microstructural inhomogeneity, a random variable technique is used in the model. A cumulative probability of the crack length at an arbitary time, G(bar a, bar t), and that of the time when a crack reaches an arbitary length, F(bar t, bar a), are obtained numerically by means of a Monte Carlo method. G(bar a, bar t), and F(bar t, bar a) are the probabilities for a single crack. However, multiple cracks generally initiate on the surface of a smooth specimen from the early stage of creep life to the final stage. TAking into account the multiple crack initiations, the actual crack length distribution observed on the surface of a specimen is predicted by the combination of probabilities for a single crack. The prediction shows a fairly good agreement with the experimental result for creep of 304 stainless steel at 923 K. The probability of creep life is obtained from an assumption that creep fracture takes place when the longest crack reaches a critical length. The observed and predicted scatter of the life is fairly small for the specimens tested.

Kitamura, Takayuki↗

NASA NTRS1989

Creep life prediction based on stochastic model of microstructurally short crack growth

A nondimensional model of microstructurally short crack growth in creep is developed based on a detailed observation of the creep fracture process of 304 stainless steel. In order to deal with the scatter of small crack growth rate data caused by microstructural inhomogeneity, a random variable technique is used in the model. A cumulative probability of the crack length at an arbitrary time, G(bar a, bar t), and that of the time when a crack reaches an arbitrary length, F(bar t, bar a), are obtained numerically by means of a Monte Carlo method. G(bar a, bar t), and F(bar t, bar a) are the probabilities for a single crack. However, multiple cracks generally initiate on the surface of a smooth specimen from the early stage of creep life to the final stage. Taking into account the multiple crack initiations, the actual crack length distribution observed on the surface of a specimen is predicted by the combination of probabilities for a single crack. The prediction shows a fairly good agreement with the experimental result for creep of 304 stainless steel at 923 K. The probability of creep life is obtained from an assumption that creep fracture takes place when the longest crack reaches a critical length. The observed and predicted scatter of the life is fairly small for the specimens tested.

Kitamura, Takayuki↗

NASA NTRS1992

Static assignment of complex stochastic tasks using stochastic majorization

We consider the problem of statically assigning many tasks to a (smaller) system of homogeneous processors, where a task's structure is modeled as a branching process, and all tasks are assumed to have identical behavior. We show how the theory of majorization can be used to obtain a partial order among possible task assignments. Our results show that if the vector of numbers of tasks assigned to each processor under one mapping is majorized by that of another mapping, then the former mapping is better than the latter with respect to a large number of objective functions. In particular, we show how measurements of finishing time, resource utilization, and reliability are all captured by the theory. We also show how the theory may be applied to the problem of partitioning a pool of processors for distribution among parallelizable tasks.

Nicol, David↗

NASA NTRS2006

Dynamic Theory of Relativistic Electrons Stochastic Heating by Whistler Mode Waves with Application to the Earth Magnetosphere

In the Hamiltonian approach an electron motion in a coherent packet of the whistler mode waves propagating along the direction of an ambient magnetic field is studied. The physical processes by which these particles are accelerated to high energy are established. Equations governing a particle motion by group symmetries of the problem were transformed in to a closed pair of nonlinear difference equations. The solutions of these equations have shown there exists the energetic threshold below that the electron motion is regular, and when the initial energy is above the threshold an electron moves stochastically. It is proved that the upper boundary of particle stochastic heating is conditioned by intrinsic property of the particle chaotic motion. Particle energy spectra and pitch angle electron scattering are described by the Fokker-Planck-Kolmogorov equations. It is shown that significant pitch angle diffusion occurs for the Earth radiation belt electrons with energies from a few keV up to a few MeV.

Khazanov, G. V.↗

DOE OSTI2024

Progressive Hedging Decomposition for Solutions of Large-Scale Process Family Design Problems

In previous work, we have introduced a mathematical model for solving a discretized version of the process family design problem. This involves two sets of decision variables. One set selects which unit module designs are included in the process platform out of a candidate set of options; the other set determines which of these unit module designs are assigned to each variant. In this work, we exploit a parallelized Progressive Hedging (PH) algorithm to solve even larger scale design problems. PH is a well-known algorithm traditionally used to solve stochastic programming problems. While our problem is not a two-stage stochastic programming problem, the structure is similar, and it can be directly mapped to the PH approach, which we employ here to solve this deterministic optimization problem. We decompose our problem by process variant. We treat the platform unit module design variables as first-stage and the assignment of unit module designs to variants as second-stage, solving the problem using mpi-sppy. We demonstrate this approach on case studies of CC, water desalination, and refrigeration.

Stinchfield, Georgia↗

DOE OSTI2024

Towards provably efficient quantum algorithms for large-scale machine-learning models

Large machine learning models are revolutionary technologies of artificial intelligence whose bottlenecks include huge computational expenses, power, and time used both in the pre-training and fine-tuning process. In this work, we show that fault-tolerant quantum computing could possibly provide provably efficient resolutions for generic (stochastic) gradient descent algorithms, scaling as $\mathcal{O}$(T 2 x polylog($n$)), where n is the size of the models and T is the number of iterations in the training, as long as the models are both sufficiently dissipative and sparse, with small learning rates. Based on earlier efficient quantum algorithms for dissipative differential equations, we find and prove that similar algorithms work for (stochastic) gradient descent, the primary algorithm for machine learning. In practice, we benchmark instances of large machine learning models from 7 million to 103 million parameters. We find that, in the context of sparse training, a quantum enhancement is possible at the early stage of learning after model pruning, motivating a sparse parameter download and re-upload scheme. Our work shows solidly that fault-tolerant quantum algorithms could potentially contribute to most state-of-the-art, large-scale machine-learning problems.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Data-Driven Global Sensitivity Analysis of Three-Phase Distribution System with PVs

Global sensitivity analysis (GSA) of distribution system with respect to stochastic PV variations plays an important role in designing optimal voltage control schemes. This paper proposes a Kriging, i.e., Gaussian process modeling enabled data-driven GSA method. The key idea is to develop a surrogate model that captures the hidden global relationship between voltage and real and reactive power injections from the historical data. With the surrogate model, the Sobol index can be conveniently calculated to assess the global sensitivity of voltage to various power injection variations. Comparison results with other model-based GSA methods on the IEEE 37-bus feeder, such as the polynomial chaos expansion and the Monte Carlo approaches demonstrate that the proposed method can achieve accurate GSA outcomes while maintaining high computational efficiency.

distribution system analysis↗

NASA NTRS1995

The Elastic Response of the Earth to Interannual Variations in Antarctic Precipitation

Measurements of elastic displacements of the bedrock surrounding large ice sheets have been proposed as a means to detect mass changes in these ice sheets. However, accumulation of glacial mass on the ice sheets is a noisy process, subject to large spatial and temporal variations in precipitation. We simulated the response of the Antarctic continent to a stochastic model of interannual precipitation variations and found that interannual variations in the elastic response of the earth are large when compared to the long-term mean of displacements produced by an assumed average ice mass imbalance of 10%. If, as some scientists predict, Antarctic ice mass changes in the future become dramatic, the long-term signal should be large enough to be detected by a few years of geodetic measurements, despite climatic noise.

Hager, Bradford H.↗

DOE OSTI2021

Modeling and simulation of vascular tumors embedded in evolving capillary networks

Here in this work, we present a coupled 3D–1D model of solid tumor growth within a dynamically changing vascular network to facilitate realistic simulations of angiogenesis. Additionally, the model includes erosion of the extracellular matrix, interstitial flow, and coupled flow in blood vessels and tissue. We employ continuum mixture theory with stochastic Cahn–Hilliard type phase-field models of tumor growth. The interstitial flow is governed by a mesoscale version of Darcy’s law. The flow in the blood vessels is controlled by Poiseuille flow, and Starling’s law is applied to model the mass transfer in and out of blood vessels. The evolution of the network of blood vessels is orchestrated by the concentration of the tumor angiogenesis factors (TAFs); blood vessels grow towards the increasing TAFs concentrations. This process is not deterministic, allowing random growth of blood vessels and, therefore, due to the coupling of nutrients in tissue and vessels, makes the growth of tumors stochastic. We demonstrate the performance of the model by applying it to a variety of scenarios. Numerical experiments illustrate the flexibility of the model and its ability to generate satellite tumors. Simulations of the effects of angiogenesis on tumor growth are presented as well as sample-independent features of cancer.

3D–1D coupled blood flow models↗

DOE OSTI2025

Effect of pattern transfer process on roughness of block copolymer patterns from directed self-assembly

Block copolymer-directed self-assembly (DSA) remains promising for improving pattern quality and reducing the stochastic variations that challenge high numerical aperture extreme ultraviolet lithography. Equally critical is refining pattern transfer methods for accurately transferring the rectified DSA features to the underlying substrate. We compare two atomic layer deposition (ALD)-based techniques: sequential infiltration synthesis (SIS) and dry liftoff, applied to polystyrene-block-poly(methyl methacrylate) (PS-b-PMMA) DSA patterns. Both methods utilize aluminum oxide hard masks, with one synthesized through infiltration into the PMMA domains and the other through conformal ALD coating. High-resolution scanning electron micrographs were analyzed to measure line edge, width, and placement roughness for both the line (PMMA) and space (PS) features. Although both methods yielded similar overall 3σ rms roughness, they differed significantly in the frequency-dependent power spectral density (PSD) profiles. SIS reduced line placement roughness at length scales associated with the polymer pitch, but increased space width roughness at low frequencies, whereas dry liftoff mimicked the frequency content of the original guiding pattern. This study underscores the importance of PSD evaluation in selecting optimal pattern transfer strategies for specific applications.

Block copolymers↗

DOE OSTI2021

Stochastic microstructural analysis of ceramic matrix composites using a high-fidelity multiscale framework

This paper presents a high-fidelity multiscale framework to simulate the mechanical behavior of ceramic matrix composites (CMCs), accounting for the complex micromorphology captured using detailed material characterization. First, high-resolution micrographs are obtained for the specific carbon fiber silicon-carbide-nitride matrix (C/SiNC) CMC to characterize the variability of the architectural features and manufacturing-induced defects at the microscale. An image processing algorithm is then used to precisely estimate the size and distribution of all subscale features and defects from the micrographs. The information is then used to generate a three-dimensional stochastic representative volume element (SRVE) to reconstruct microscale constituents accounting for the variability. As a result, the generated SRVEs are simulated using the high-fidelity generalized method of cells (HFGMC) micromechanics theory to investigate the effects of defects on the elastic properties of C/SiNC CMCs.

Ceramic Matrix Composites↗

DOE OSTI2023

Stochastic Price Generation for Evaluating Wholesale Electricity Market Bidding Strategies

This work presents a novel method for generating electricity price scenarios from statistical properties of past electricity prices using a hybrid statistical and reduced-form stochastic model. Previous work in applying stochastic differential equations (SDE) to model electricity prices has focused on daily average prices. To extend stochastic price generation methods to hourly or sub-hourly pricing, we address several weaknesses in the state-of-the-art: (1) we replace the mean-reversion component of the SDE with an ARIMA process that is better able to characterize the daily and weekly trends; (2) we extend the price-spike, or jump process to account for conditional probabilities of price spikes occurring in consecutive time steps by replacing the traditional Poisson process for modeling jumps with a generalized point process model inspired by brain neuron models; and (3) we replace the traditional method of estimating spike intensity with empirical variance with a Markov process based on observed price spike intensity transitions. The method is demonstrated with electricity prices from the US ERCOT market and a use-case example is provided for bidding an energy storage unit into the day-ahead and real-time energy markets of ERCOT using stochastic optimization methods. Results show that the the synthetic price model out performs a (naive) persistence forecast model by resulting in 24% to 47% more in profits over 168 simulated days.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

NASA NTRS1995

High Precision and Real Time Tracking of Low Earth Orbiters With GPS: Case Studies With TOPEX/POSEIDON and EUVE

TOPEX/POSEIDON carries a dual-frequency 6 channel GPS receiver while EUVE has a 12 channel single frequency receiver. Flying at an altitude of 1334 km, TOPEX/POSEIDON performs precise ocean altimetry, which demands the highest possible accuracy in determining the radial orbit component in post-processing. Radial RMS accuracies of about 2 cm were realized using reduced dynamic tracking techniques. In this approach, orbit errors due to force are substantially reduced by exploiting the geometric strength of GPS to solve for a set of stochastic forces. On EUVE, the emphasis was on evaluating real time positioning techniques with a single frequency receiver. The capability for real time 3D accuracies of 15 m in the presence of Selective Availability was shown. This was validated by comparing to a post-processed differential GPS truth orbit believed accurate to about 1 m.!.

Tracking low-Earth orbiters↗

DOE OSTI2020

A Two-Stage Decomposition Approach for AC Optimal Power Flow

The alternating current optimal power flow (AC-OPF) problem is critical to power system operations and planning, but it is generally hard to solve due to its nonconvex and large-scale nature. Furthermore, this paper proposes a scalable decomposition approach in which the power network is decomposed into a master network and a number of subnetworks, where each network has its own AC-OPF subproblem. This formulates a two-stage optimization problem and requires only a small amount of communication between the master and subnetworks. The key contribution is a smoothing technique that renders the response of a subnetwork differentiable with respect to the input from the master problem, utilizing properties of the barrier problem formulation that naturally arises when subproblems are solved by a primal-dual interior-point algorithm. Consequently, existing efficient nonlinear programming solvers can be used for both the master problem and the subproblems. The advantage of this framework is that speedup can be obtained by processing the subnetworks in parallel, and it has convergence guarantees under reasonable assumptions. The formulation is readily extended to instances with stochastic subnetwork loads. Numerical results show favorable performance and illustrate the scalability of the algorithm which is able to solve instances with more than 11 million buses.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

Development of a three-dimensional APOLLO3 neutrons deterministic scheme for the CABRI reactor

CABRI is an experimental reactor to study the fuel behavior during reactivity injection transients. These transients being highly multiphysics, the development of suitable modeling and simulation tools to simulate them is important for the optimization of the tests and the control of the experimental conditions. This paper focuses on the development of an APOLLO3 deterministic core calculation tool dedicated to the CABRI transient analysis. It represents the first stage of the incremental process for the implementation of a multiphysics time-dependent modeling of the CABRI transient. The neutron calculation scheme is based on a classical two-step approach. The first step consists of a 281-energy group calculation flux with the TDT-MOC (Method Of Characteristics) solver for cross-section space and energy (23 groups) collapsing for the CABRI different assembly clusters. The bias on a 2D core neutron calculation due to the self-shielding calculation and collapsing on a restricted pattern are investigated thanks to a comparison with a direct full 2D calculation on a quarter of core. The second step relies on a pin-resolved transport 3D transport core calculation with the SN solver MINARET. A progressive numerical validation process is followed to quantify the calculation biases on reactivity and reaction rates at each step using reference calculations with the stochastic code TRIPOLI4. The next development stage toward a multiphysics scheme will be the implementation of the 3D-kinetics equation resolution and the coupling with a core thermal-hydraulics model. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRS1992

A review of transport theory

Ways in which energy change terms arise in the transport equation and how the various terms relate to the modes of energy exchange between the particles and plasma are shown. It is argued that the transport equation cannot be used to describe the initial acceleration of thermal particles by plasma shocks or relativistic shocks where the energetic particle speeds are never much greater than the flow speeds. In most other situations, it describes almost any acceleration process that can be caused by a moving plasma. It describes shock acceleration for both parallel shocks and oblique ones, and stochastic acceleration by the turbulent motion of the scatterers as well as by their motion across the magnetic field.

Jones, Frank C.↗