Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stochastic optimization”

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,477 records · Page 82

Selective Shielding Bone Marrow: An Approach to Create Personal Protective Equipment for Gamma Radiation

Many casualties of the Hiroshima and Nagasaki atomic bombs, and many of the firefighters who first responded to the Chernobyl nuclear power plant accident, became ill with acute radiation syndrome (ARS). Most of the people who do not recover from ARS will die within a few weeks to a few months after exposure, with the primary cause of death being the destruction of the body's bone marrow (BM). This is exemplified by the Chernobyl firefighters, the majority of whom received whole body doses between 0.8 and 10 Gy, and the United Nations Scientific Committee on the Effects of Atomic Radiation (UNSCEAR) final report on the disaster concluded that the underlying cause of death of all 28 firefighters who succumbed to ARS was from bone marrow failure. BM is comprised of hematopoietic stem cells (HSCs) which are responsible for the constant renewal of blood cells. Due to their high rate of proliferation, HSCs are especially vulnerable to ionizing radiation, but are endowed with remarkable regenerative potential. Owing to their central role in blood production, lethal irradiation of HSCs leads to death from severe anemia, infection and internal bleeding. This relationship between high doses of radiation and HSC apoptosis has led to the use of HSC transplantation as a life-saving intervention in cases of acute exposure. This is demonstrable in medical practice, where thousands of individuals have undergone supra-lethal Total Body Irradiation (TBI) for purposes of cancer therapy, and were rescued by well-matched BM transplantation. Indeed, life-threatening damage may be reversed by BM transplantation in individuals receiving doses of radiation of up to 10 Gy. At doses > 10 Gy, damage to gastrointestinal (GI) tissue may become a limiting factor in survival. Remarkably, in the common procedure of BM transplantation, the number of HSCs extracted from a single active BM site containing less than 5% of the donor's BM tissue is sufficient to support the complete reconstitution of the HSC compartment in a lethally irradiated recipient. This capacity of BM to expand and replenish is due to the high regenerative potential of the HSCs it harbors. In contrast to mature blood cells, which are dispersed throughout the body, HSCs are confined to the bones, allowing for effective targeted shielding. This has been confirmed in several animal models, where sparing the HSC-rich area of a subject receiving otherwise full-body irradiation is sufficient to support hematopoietic functions and allow survival. Approximately 48% of the human body's active BM is contained within the lumbar vertebrae and the pelvic region. This high concentration of BM makes the pelvic region an ideal area to shield for preventing the serious hematopoietic effects of radiation exposure. Stem-cell rich tissues like active BM are also more sensitive to the stochastic effects (i.e. cancer) of radiation exposure. The high amount of BM in the pelvic region, combined with the presence of the sensitive gastrointestinal system and female gonads in this same area of the body allows for the possibility of significantly reducing stochastic effects of radiation by limiting radiation exposure to this vital area. Thus, shielding this region holds great promise for both acute and protracted exposures. Thus, a specialized radiation shielding device for the protection of active BM concentrations is presented here. This belt-like selective shielding device focuses on the protection of BM that is present in the pelvis and hip bones, protecting the medullary volume from which BM is commonly extracted for transplantation (i.e. the iliac bones) while allowing relatively unhindered movement of the wearer. In order to optimize the use of shielding materials towards the protection of active BM, the shielding uniquely brings into account the natural shielding properties of human tissue, with thickness being inversely related to the thickness and radiation attenuation of the underlying tissue at each point surrounding the area being protected. The findings resulting from both experimental and simulation testing of this selective shielding device are presented here.

61 RADIATION PROTECTION AND DOSIMETRY↗

The Role of Concentrating Solar Power Toward High Renewable Energy Penetrated Power Systems

Achieving high renewable energy penetrated power systems requires considerable operational flexibility to hedge the variability and uncertainty of variable renewable energy (VRE) generation. Compared with VRE sources, concentrating solar power (CSP) is an emerging controllable renewable generation technique that utilizes solar thermal power to generate electricity. The operational dispatchability of CSP would contribute to the power system transition toward high renewable penetration. In this paper, we explore how the generation portfolio will change toward high renewable energy penetrations, how much cost is involved, and what role CSP will play in realizing a high renewable energy penetrated power system. This study relies on a stochastic two-stage generation and transmission expansion planning model with CSP plants. The model captures the uncertainty and variability of renewable generation and the flexibility limits of thermal plants. With the target of achieving a renewable-dominated minimum-cost system with an expected renewable energy penetration level, the investments of both generation and transmission facilities are optimized. A case study on IEEE test systems with renewable technology cost data in 2050 is performed to analyze the value of CSP toward high renewable energy penetrated power systems.

14 SOLAR ENERGY↗

A stochastic scan strategy for grain structure control in complex geometries using electron beam powder bed fusion

Spatial control of microstructure within a three-dimensional component has been a dream of materials scientists for centuries. However, limitations in traditional manufacturing processes prevent detailed control over the distribution of microstructures in a single part. Here, we demonstrate the ability to control grain structure and crystallographic texture during metal additive manufacturing for arbitrary cross-sections of a practical size, with profound implications for the design and optimization of next-generation products. The key to this advance is a new geometry agnostic scan path algorithm that manipulates the spatial distribution of solidification conditions. Utilizing a fundamental understanding of solidification dynamics and a model of the heat transfer during processing, we have designed this algorithm to manipulate the natural competition between epitaxial dendrite growth and grain nucleation. With this algorithm, we successfully controlled the grain structure of Ni-based superalloy IN718 in the shape of the Mona Lisa.

36 MATERIALS SCIENCE↗

THE BENEFITS OF VLBI ASTROMETRY TO PULSAR TIMING ARRAY SEARCHES FOR GRAVITATIONAL RADIATION

Precision astrometry is an integral component of successful pulsar timing campaigns. Astrometric parameters are commonly derived by fitting them as parameters of a timing model to a series of pulse times of arrival (TOAs). TOAs measured to microsecond precision over spans of several years can yield position measurements with sub-milliarcsecond precision. However, timing-based astrometry can become biased if a pulsar displays any red spin noise or a red signal produced by the stochastic gravitational wave background. We investigate how noise of different spectral types is absorbed by timing models, leading to significant estimation biases in the astrometric parameters. We find that commonly used techniques for fitting timing models in the presence of red noise (Cholesky whitening) prevent the absorption of noise into the timing model remarkably well if the time baseline of observations exceeds several years, but are inadequate for dealing with shorter pulsar data sets. Independent of timing, pulsar-optimized very long baseline interferometry (VLBI) is capable of providing position estimates precise to the sub-milliarcsecond levels needed for high-precision timing. In order to make VLBI astrometric parameters useful in pulsar timing models, the transformation between the International Celestial Reference Frame (ICRF) and the dynamical solar system ephemeris used for pulsar timing must be constrained to within a few microarcseconds. We compute a transformation between the ICRF and pulsar timing frames and quantitatively discuss how the transformation will improve in coming years. We find that incorporating VLBI astrometry into the timing models of pulsars for which only a couple of years of timing data exist will lead to more realistic assessments of red spin noise and could enhance the amplitude of gravitational wave signatures in post-fit timing residuals by factors of 20 or more.

79 ASTRONOMY AND ASTROPHYSICS↗

Mixing state evolution of agglomerating particles in an aerosol chamber: Comparison of measurements and particle-resolved simulations

This article introduces a validation study of the stochastic particle-resolved aerosol model PartMC with experimental data from an aerosol chamber experiment. For the experiment, a scanning mobility particle sizer and a single-particle soot photometer were used to monitor the aerosol mixing state evolution of two initially externally-mixed aerosol populations of ammonium sulfate and black carbon particles undergoing agglomeration. We applied an efficient optimization algorithm (ProSRS) to determine several unconstrained simulation parameters and were able to successfully reproduce number concentrations and size distributions of mixed particles that formed by agglomeration. The PartMC modeling method in conjunction with the optimization procedure provides a tool for detailed comparisons of chamber experiments and modeling, where aerosol mixing state is the focus of investigation.

54 ENVIRONMENTAL SCIENCES↗

Post-Disaster Microgrid Formation for Enhanced Distribution System Resilience

This paper proposes a deep reinforcement learning (DRL) based approach for post-disaster critical load restoration in active distribution systems to form microgrids through network reconfiguration to minimize critical load curtailments. Distribution networks are represented as graph networks, and optimal network configurations with microgrids are obtained by searching for the optimal spanning forest. The constraints to the research question being explored are the radial topology and power balance. Unlike existing analytical and population-based approaches, which necessitate the repetition of entire analyses and computation for each outage scenario to find the optimal spanning forest, the proposed approach, once properly trained, can quickly determine the optimal, or near-optimal, spanning forest even when outage scenarios change. When multiple lines fail in the system, the proposed approach forms microgrids with distributed energy resources in active distribution systems to reduce critical load curtailment. The proposed DRL-based model learns the action-value function using the REINFORCE algorithm, which is a model-free reinforcement learning technique based on stochastic policy gradients. A case study was conducted on a 33-node distribution test system, demonstrating the effectiveness of the proposed approach for post-disaster critical load restoration.

active distribution systems↗

Structure–Property Relationships of Recycled Lithium-Ion Battery Cathodes: Microstructure Optimization Using Virtual Materials Testing

The increasing demand for sustainable battery technologies requires effective recycling strategies for end-of-life lithium-ion battery cathodes. In this study, virtual materials testing, a well-established framework for modeling conventionally manufactured NMC-based cathodes, is applied to partially recycled cathodes. To this end, virtual cathodes consisting of mixtures of pristine and recycled NMC particles are utilized to systematically analyze structure–property relationships depending on mixing ratios and different spatial arrangement strategies. For this purpose, a stochastic 3D model is developed that is capable of generating virtual cathodes with arbitrary volume fractions of active materials and mixing ratios of pristine and recycled NMC particles. Particularly, the stochastic 3D model can mimic the different size distributions of pristine and recycled particles that are observed in image data. Additionally, the model allows the structuring of pristine and recycled NMC either uniformly mixed or layer-wise arranged, mimicking single- and dual-layer cathodes. Subsequently, a systematic computational analysis is conducted to assess the influence of increasing active material ratios of recycled particles, ranging from 0 % to 100 %, while maintaining a constant overall active material volume fraction. The impact of particle mixing on cathode performance is evaluated by examining transport-relevant geometrical descriptors and effective properties, such as geodesic tortuosity, specific surface area, and tortuosity factor.

25 ENERGY STORAGE↗

Optimization and simulation of flight control laws under parameter uncertainty and external disturbances

Several tasks pertinent to flight control in parameter uncertainty and wind-gust loading were successfully completed. Identification algorithms for extracting stability and control derivatives from flight data taking gust loading into account were developed. They were verified by simulation and evaluated throughly on actual flight data taken on a Lockheed Jet Star flying in turbulence. In particular the need for automatically generated dither-like inputs was studied. Criteria for performance evaluation using stochastic models were developed for gust alleviation as well as handling quantities. Algorithms for assessing degradation in performance due to parameter uncertainty were developed and evaluated using flight test data.

Source record↗

A finite state, finite memory minimum principle, part 2

The development of the theory of the finite - state, finite - memory (FSFM) stochastic control problem is discussed. The sufficiency of the FSFM minimum principle (which is in general only a necessary condition) was investigated. By introducing the notion of a signaling strategy as defined in the literature on games, conditions under which the FSFM minimum principle is sufficient were determined. This result explicitly interconnects the information structure of the FSFM problem with its optimality conditions. The min-H algorithm for the FSFM problem was studied. It is demonstrated that a version of the algorithm always converges to a particular type of local minimum termed a person - by - person extremal.

Sandell, N. R., Jr.↗

A framework for stochastic simulation of distribution practices for hotel reservations

The focus of this study is primarily on the Greek hotel industry. The objective is to design and develop a framework for stochastic simulation of reservation requests, reservation arrivals, cancellations and hotel occupancy with a planning horizon of a tourist season. In Greek hospitality industry there have been two competing policies for reservation planning process up to 2003: reservations coming directly from customers and a reservations management relying on tour operator(s). Recently the Internet along with other emerging technologies has offered the potential to disrupt enduring distribution arrangements. The focus of the study is on the choice of distribution intermediaries. We present an empirical model for the hotel reservation planning process that makes use of a symbolic simulation, Monte Carlo method, as, requests for reservations, cancellations, and arrival rates are all sources of uncertainty. We consider as a case study the problem of determining the optimal booking strategy for a medium size hotel in Skiathos Island, Greece. Probability distributions and parameters estimation result from the historical data available and by following suggestions made in the relevant literature. The results of this study may assist hotel managers define distribution strategies for hotel rooms and evaluate the performance of the reservations management system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

A stochastic perturbation method to generate inflow turbulence in large-eddy simulation models: Application to neutrally stratified atmospheric boundary layers

Despite the variety of existing methods, efficient generation of turbulent inflow conditions for large-eddy simulation (LES) models remains a challenging and active research area. Herein, we extend our previous research on the cell perturbation method, which uses a novel stochastic approach based upon finite amplitude perturbations of the potential temperature field applied within a region near the inflow boundaries of the LES domain [Muñoz-Esparza et al., “Bridging the transition from mesoscale to microscale turbulence in numerical weather prediction models,” Boundary-Layer Meteorol., 153, 409–440 (2014)]. The objective was twofold: (i) to identify the governing parameters of the method and their optimum values and (ii) to generalize the results over a broad range of atmospheric large-scale forcing conditions, U{sub g} = 5 − 25 m s{sup −1}, where U{sub g} is the geostrophic wind. We identified the perturbation Eckert number, Ec=U{sub g}{sup 2}/ρc{sub p}θ{sup ~}{sub pm}, to be the parameter governing the flow transition to turbulence in neutrally stratified boundary layers. Here, θ{sup ~}{sub pm} is the maximum perturbation amplitude applied, c{sub p} is the specific heat capacity at constant pressure, and ρ is the density. The optimal Eckert number was found for nonlinear perturbations allowed by Ec ≈ 0.16, which instigate formation of hairpin-like vortices that most rapidly transition to a developed turbulent state. Larger Ec numbers (linear small-amplitude perturbations) result in streaky structures requiring larger fetches to reach the quasi-equilibrium solution, while smaller Ec numbers lead to buoyancy dominated perturbations exhibiting difficulties for hairpin-like vortices to emerge. Cell perturbations with wavelengths within the inertial range of three-dimensional turbulence achieved identical quasi-equilibrium values of resolved turbulent kinetic energy, q, and Reynolds-shear stress, . In contrast, large-scale perturbations acting at the production range exhibited reduced levels of , due to the formation of coherent streamwise structures, while q was maintained, requiring larger fetches for the turbulent solution to stabilize. Additionally, the cell perturbation method was compared to a synthetic turbulence generator. The proposed stochastic approach provided at least the same efficiency in developing realistic turbulence, while accelerating the formation of large-scales associated with production of turbulent kinetic energy. Also, it is computationally inexpensive and does not require any turbulent information.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Stochastic response analysis, order reduction, and output feedback controllers for flexible spacecraft

Real disturbances and real sensors have finite bandwidths. The first objective of this paper is to incorporate this finiteness in the 'open-loop modal cost analysis' as applied to a flexible spacecraft. Analysis based on residue calculus shows that among other factors, significance of a mode depends on the power spectral density of disturbances and the response spectral density of sensors at the modal frequency. The second objective of this article is to compare performances of an optimal and a suboptimal output feedback controller, the latter based on 'minimum error excitation' of Kosut. Both the performances are found to be nearly the same, leading us to favor the latter technique because it entails only linear computations. Our final objective is to detect an instability due to truncated modes by representing them as a multiplicative and an additive perturbation in a nominal transfer function. In an example problem it is found that this procedure leads to a narrow range of permissible controller gains, and that it labels a wrong mode as a cause of instability. A free beam is used to illustrate the analysis in this work.

Hablani, H. B.↗

A hybrid inventory management system respondingto regular demand and surge demand

This paper proposes a hybrid policy for a stochastic inventory system facing regular demand and surge demand. The combination of two different demand patterns can be observed in many areas, such as healthcare inventory and humanitarian supply chain management. The surge demand has a lower arrival rate but higher demand volume per arrival. The solution approach proposed in this paper incorporates the level crossing method and mixed integer programming technique to optimize the hybrid inventory policy with both regular orders and emergency orders. The level crossing method is applied to obtain the equilibrium distributions of inventory levels under a given policy. The model is further transformed into a mixed integer program to identify an optimal hybrid policy. A sensitivity analysis is conducted to investigate the impact of parameters on the optimal inventory policy and minimum cost. Numerical results clearly show the benefit of using the proposed hybrid inventory model. The model and solution approach could help healthcare providers or humanitarian logistics providers in managing their emergency supplies in responding to surge demands.

99 GENERAL AND MISCELLANEOUS↗

A Large Deviations Analysis of Certain Qualitative Properties of Parallel Tempering and Infinite Swapping Algorithms

Parallel tempering, or replica exchange, is a popular method for simulating complex systems. The idea is to run parallel simulations at different temperatures, and at a given swap rate exchange configurations between the parallel simulations. From the perspective of large deviations it is optimal to let the swap rate tend to infinity and it is possible to construct a corresponding simulation scheme, known as infinite swapping. In this paper we propose a novel use of large deviations for empirical measures for a more detailed analysis of the infinite swapping limit in the setting of continuous time jump Markov processes. Using the large deviations rate function and associated stochastic control problems we consider a diagnostic based on temperature assignments, which can be easily computed during a simulation. Here we show that the convergence of this diagnostic to its a priori known limit is a necessary condition for the convergence of infinite swapping. Finally, the rate function is also used to investigate the impact of asymmetries in the underlying potential landscape, and where in the state space poor sampling is most likely to occur.

97 MATHEMATICS AND COMPUTING↗

Constructing rigorous and broad biosurveillance networks for detecting emerging zoonotic outbreaks

Determining optimal surveillance networks for an emerging pathogen is difficult since it is not known beforehand what the characteristics of a pathogen will be or where it will emerge. The resources for surveillance of infectious diseases in animals and wildlife are often limited and mathematical modeling can play a supporting role in examining a wide range of scenarios of pathogen spread. We demonstrate how a hierarchy of mathematical and statistical tools can be used in surveillance planning help guide successful surveillance and mitigation policies for a wide range of zoonotic pathogens. The model forecasts can help clarify the complexities of potential scenarios, and optimize biosurveillance programs for rapidly detecting infectious diseases. Using the highly pathogenic zoonotic H5N1 avian influenza 2006-2007 epidemic in Nigeria as an example, we determined the risk for infection for localized areas in an outbreak and designed biosurveillance stations that are effective for different pathogen strains and a range of possible outbreak locations. We created a general multi-scale, multi-host stochastic SEIR epidemiological network model, with both short and long-range movement, to simulate the spread of an infectious disease through Nigerian human, poultry, backyard duck, and wild bird populations. We chose parameter ranges specific to avian influenza (but not to a particular strain) and used a Latin hypercube sample experimental design to investigate epidemic predictions in a thousand simulations. We ranked the risk of local regions by the number of times they became infected in the ensemble of simulations. These spatial statistics were then complied into a potential risk map of infection. Finally, we validated the results with a known outbreak, using spatial analysis of all the simulation runs to show the progression matched closely with the observed location of the farms infected in the 2006-2007 epidemic.

59 BASIC BIOLOGICAL SCIENCES↗

Exploring Interferometry Diagnostics for Optical Stochastic Cooling at FAST/IOTA

Optical Stochastic Cooling (OSC) is an advanced beam-cooling technique that will precede the traditional stochastic cooling. This method leverages optical radiation and high-precision feedback to cool the particles more efficiently than traditional stochastic methods by more than three orders of magnitude. As such, it is an enabler for the development of next generation discovery science machines at the frontiers of energy and intensity. This paper focuses on the development of the second phase of OSC and improving extreme beam cooling technique in accelerators. Along with the assembly and use of an interferometer to provide diagnostics and optimization to OSC systems. The interferometer can allow for the observation and analysis of fringe patterns from a recreated simpler version in the laser room. By leveraging interferometric techniques, OSC can achieve higher efficiency, verification, stability and performance.

Teriba, Folashade↗

High fidelity multiphysics tightly coupled model for a lead cooled fast reactor concept and application to statistical calculation of hot channel factors

A tightly coupled multiphysics code system is established using the MOOSE framework for hot channel factor (HCF) evaluation on a Lead Fast Reactor (LFR) concept. The coupled system is driven by the Griffin multiphysics coupling capability under which the MOOSE Heat Transfer module and NekRS computational fluid dynamics solver are coupled for conjugate heat transfer using the Cardinal application. The coupled capability is demonstrated on an LFR assembly model based on materials and geometry of a prototypical lead-cooled fast reactor design by Westinghouse Electric Company, LLC. Moreover, the work integrates the Multiphysics Object Oriented Simulation Environment (MOOSE) Stochastic Tools Module (STM) to perform calculations for statistical analysis of HCF. Furthermore, the coupling strategy and workflow demonstrated in this paper is not only useful for predicting accurate hot channel factors for different kinds of advanced reactors but also for other engineering applications such as control rod worth assessment, generation of high-fidelity database for Artificial intelligence (AI)/machine learning (ML) training, design optimization and multi-resolution modeling.

Cardinal↗