Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stochastic programming”

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 397 records · Page 22

A REDUCED ORDER MODELING APPROACH TO PROBABILISTIC CREEP-DAMAGE PREDICTIONS IN FINITE ELEMENT ANALYSIS

This paper introduces a computationally efficient Reduced Order Modeling (ROM) approach for the probabilistic prediction of creep-damage failure. Component-level probabilistic simulations are needed to assess the reliability and safety of high-temperature components. Full-scale probabilistic creep-damage modeling in finite element (FE) approach is computationally expensive requiring many hundreds of simulations to replicate the uncertainty of component failure. To that end, ROM is proposed to minimize the elevated computational cost while controlling the loss of accuracy. It is proposed that full-scale probabilistic simulations can be completed in 1D at a reduced cost, the extremum conditions extracted, and those conditions applied for lower-cost 2D/3D probabilistic simulations of components that capture the mean and uncertainty of failure. The probabilistic Sine-hyperbolic (Sinh) model is selected which in previous work was calibrated to alloy 304 stainless steel. The Sinh model includes probability density functions (pdfs) for test condition (stress and temperature), initial damage (i.e. microstructure), and material properties uncertainty. The Sinh model is programmed into ANSYS finite element software using the USERCREEP.F material subroutine. First, the Sinh model and FE code are subject to verification and validation to ensure the accuracy of the simulations. Numerous Monte Carlo simulations are executed in a 1D model to generate probabilistic creep deformation, damage, and rupture data. This data is analyzed and the probabilistic parameters corresponding to extreme creep response are extracted. The ROM concept is applied where only the extreme conditions are applied in the 2D probabilistic prediction of a component. The probabilistic predictions between the 1D and 2D geometry is compared to assess ROM for creep. The accuracy of the probabilistic prediction employing the ROM approach will potentially reduce the time and cost of simulating complex engineering systems. Future studies will introduce multi-stage Sinh, stochasticity, and spatial uncertainty for improved prediction.

36 MATERIALS SCIENCE↗

An effective description of charge diffusion and energy transport in a charged plasma from holography

We discuss the physics of sound propagation and charge diffusion in a plasma with non-vanishing charge density. Our analysis culminates the program initiated in to construct an open effective field theory of low-lying modes of the stress tensor and charge current in such plasmas. We model the plasma holographically as a Reissner-Nordström-AdS d+1 black hole, and study linearized fluctuations of longitudinally polarized scalar gravitons and photons in this background. We demonstrate that the perturbations can be decoupled and repackaged into the dynamics of two designer scalars, whose gravitational coupling is modulated by a non-trivial dilatonic factor. The holographic analysis allows us to isolate the phonon mode from the charge diffusion mode, and identify the combination of currents that corresponds to each of them. We use these results to obtain the real-time Gaussian effective action, which includes both the retarded response and the associated stochastic (Hawking) fluctuations, accurate to quartic order in gradients.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A NASTRAN-based computer program for structural dynamic analysis of Horizontal Axis Wind Turbines

This paper describes a computer program developed for structural dynamic analysis of horizontal axis wind turbines (HAWT's). It is based on the finite element method through its reliance on NASTRAN for the development of mass, stiffness, and damping matrices of the tower end rotor, which are treated in NASTRAN as separate structures. The tower is modeled in a stationary frame and the rotor in one rotating at a constant angular velocity. The two structures are subsequently joined together (external to NASTRAN) using a time-dependent transformation consistent with the hub configuration. Aerodynamic loads are computed with an established flow model based on strip theory. Aeroelastic effects are included by incorporating the local velocity and twisting deformation of the blade in the load computation. The turbulent nature of the wind, both in space and time, is modeled by adding in stochastic wind increments. The resulting equations of motion are solved in the time domain using the implicit Newmark-Beta integrator. Preliminary comparisons with data from the Boeing/NASA MOD2 HAWT indicate that the code is capable of accurately and efficiently predicting the response of HAWT's driven by turbulent winds.

Lobitz, Don W.↗

Thermal neutron scattering measurements and modeling of yttrium-hydrides for high temperature moderator applications

In this work, thermal neutron scattering measurements of yttrium hydride from 5 to 1,200 K were conducted to determine the change in inelastic scattering as a function of temperature and to probe for anharmonic effects. Additionally, measurements on samples of YH x from x = 1.62 to 1.90 were done to determine the effects of varying hydrogen concentration on the inelastic spectra. Changes in temperature affected the inelastic spectra in unanticipated ways, indicating that there are anharmonic effects, whereas hydrogen concentration did not significantly affect the inelastic spectra. These measurements were compared against the ENDF/B-VIII.0 thermal scattering files of YH 2 , as well as new thermal scattering files created by using the stochastic temperature-dependent effective potential (s-TDEP). Both libraries were found to be in good agreement with the experimental data at lower temperatures. At higher temperatures, the s-TDEP method is better at predicting the experimentally observed softening of phonon modes in acoustic and optical regions.

36 MATERIALS SCIENCE↗

Fuel-injector/air-swirl characterization

The objectives of this program are to establish an experimental data base documenting the behavior of gas turbine engine fuel injector sprays as the spray interacts with the swirling gas flow existing in the combustor dome, and to conduct an assessment of the validity of current analytical techniques for predicting fuel spray behavior. Emphasis is placed on the acquisition of data using injector/swirler components which closely resemble components currently in use in advanced aircraft gas turbine engines, conducting tests under conditions that closely simulate or closely approximate those developed in actual combustors, and conducting a well-controlled experimental effort which will comprise using a combination of low-risk experiments and experiments requiring the use of state-of-the-art diagnostic instrumentation. Analysis of the data is to be conducted using an existing, TEACH-type code which employs a stochastic analysis of the motion of the dispersed phase in the turbulent continuum flow field.

Mcvey, J. B.↗

Data From: Simulating bioclogging effects on dynamic riverbed permeability and infiltration, Water Resources Research

We collected a time series of Russian River infiltration rates to parameterize stochastic model development of infiltration conditions as a function of bioclogging. The time series of infiltration were collected from the Russian River Riverbank Filtration site located in Sonoma County California. Infiltration datasets are shown in units of m/day and were obtained using a seepage meter in 2012. To address the combined effects of bioclogging and disconnection on infiltration, we developed numerical representations of bioclogging processes based on these datasets using a within a one-dimensional, variably saturated flow model representing losing-connected and losing-disconnected rivers. All models and bioclogging formulations were used to create synthetic test cases for bioclogging.This research was supported by the Jane Lewis Fellowship from the University of California, Berkeley, the Sonoma County Water Agency (SCWA), the Roy G. Post Foundation Scholarship, the U.S. Department of Energy, Office of Science Graduate Student Research (SCGSR) Program, U.S. Department of Energy, Office of Science, Office of Biological and Environmental Research under award DE-AC02-05CH11231, and the UFZ-Helmholtz Centre for Environmental Research, Leipzig, Germany.

54 ENVIRONMENTAL SCIENCES↗

A stochastic bioburden model for spacecraft sterilization.

Development of a stochastic model of the probability distribution for the random variable representing the number of microorganisms on a surface as a function of time. The first basic principle associated with bioburden estimation is that viable particles are removed from surfaces. The second notion important to the analysis is that microorganisms in environments and on surfaces occur in clumps. The last basic principle relating to bioburden modeling is that viable particles are deposited on a surface. The bioburden on a spacecraft is determined by the amount and kind of control exercised on the spacecraft assembly location, the shedding characteristics of the individuals in the vicinity of the spacecraft, its orientation, the geographical location in which the assembly takes place, and the steps in the assembly procedure. The model presented has many of the features which are desirable for its use in the spacecraft sterilization programs currently being planned by NASA.

Roark, A. L.↗

Improved Cross Section Generation Capability of Griffin in FY22

The Griffin code is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. The code includes a variety of deterministic steady-state transport solvers for fixed source, k-eigenvalue, adjoint, and subcritical multiplication as well as transient solvers for spatial dynamics with the improved quasi-static method. Griffin uses cross section data in the ISOXML format generated from external deterministic or Monte Carlo cross section generation codes. In recent years, the MC 2 -3 modules have been added to Griffin for fast reactor cross section generation, and the self-shielding application programming interface (SSAPI) was implemented in the ISOXML module for thermal reactor cross section generation. The on-the-fly slowing down method and double-heterogeneity treatment have been implemented to SSAPI and verified against particulate fuel-bearing graphite-moderated thermal reactor problems with high accuracy. This year, work has been focused on improving the cross section generation capability of ISOXML and streamlining the cross section generation procedures. In addition, the form function data were added to ISOXML in order to support the pin power reconstruction capability that was newly implemented in Griffin in this fiscal year. To facilitate the cross section generation using MC 2 -3 and SSAPI in Griffin, the cross section generation workflows have been set up for both fast and thermal spectrum reactors. The MOOSE action system tool was devised for fast spectrum problems, and the MOOSE stochastic tool was adopted to the branch calculation procedure for thermal spectrum problems. Meanwhile, to ensure the accuracy of group-constants, the thermal up-scattering kernel calculator accounting for resonance scattering was implemented in ISOXML, demonstrating the accurate computation of a Doppler-broadened scattering kernel of any Legendre order within a reasonable timescale. Other aspects of ISOXML, such as deletion solver and data, documentation, ISOXML file management, and the interface for Mixture, were improved as well.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Cost and Throughput Analysis for the NASA Ames Arc Jet Modernization Program

NASA Ames Center is currently evaluating alternatives to modernize the Arc Jet Complex, a critical part of testing for NASA’s planetary missions. NASA’s Arc Jet Complex facilities “are used to simulate the aerothermodynamic heating that a spacecraft endures throughout hypersonic atmospheric entry, and to test candidate thermal protection system (TPS) materials and systems. “Because planetary mission schedules often have tight windows due to planetary alignment constraints, a small increase in schedule could result in a two-year delay. Such a delay could increase the cost of a $1billionmissionbyhundreds of millions of dollars due to project personnel pay and clean room storage. To avoid these costs, the authors support NASA Ames in evaluating return on investment (ROI) and effectiveness of alternatives for modernizing the complex. The first input into the ROI is the deconstruction and construction cost estimates, which are developed using independent research on highly specialized subsystems, vendor quotes, and Unified Facilities Criteria (UFC), depending on the facility and work package. One of the measures of effectiveness is throughput analysis of the test bays, as a main goal of the modernization is to increase the number of possible test runs per year. This analysis is conducted via a probabilistic simulation and accounts for a variety of stochastic factors that influence the sequence of test runs, such as the facility availability; test complexity; the need to pause to assess test results; test failure; and the possibility of a system failure. The methodologies for both these analyses are discussed, along with the challenges presented due to the unique nature of the highly specialized test equipment.

Jennifer Scheel↗

Enabling Catalog Simulations of Transient and Variable Sources Based on LSST Cadence Strategies

The American Astronomical Society. All rights reserved.. The Large Synoptic Survey Telescope (LSST) project will conduct a 10 year multi-band survey starting in 2022. Observing strategies for this survey are being actively investigated, and the science capabilities can be best forecasted on the basis of simulated strategies from the LSST Operations Simulator (OpSim). This paper describes the way in which OpSim simulates a stochastic realization of the sequence of LSST pointings over the survey duration, and is based on a model of the observatory (including telescope) and historical data of observational conditions. OpSim outputs contain a record of each simulated pointing of the survey along with a complete characterization of the pointing in terms of observing conditions, and some useful quantities derived from the characteristics of the pointing. Thus, each record can be efficiently used to derive the properties of observations of all astrophysical sources found in that pointing. However, in order to obtain the time series of observations (light curves) of a set of sources, it is often more convenient to compute all observations of an astrophysical source, and iterate over sources. In this document, we describe the open source python package OpSimSummary, which allows for a convenient reordering. The objectives of this package are to provide users with an Application Programming Interface for accessing all such observations and summarizing this information in the form of intermediate data products usable by third party software such as SNANA, thereby also bridging the gap between official LSST products and preexisting simulation codes.

79 ASTRONOMY AND ASTROPHYSICS↗

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.↗

Inter-Agency Consultative Group for Space Science (IACG): Handbook of Missions and Payloads

The ACE spacecraft design is based on the Charge Composition Explorer (CCE) built by Johns Hopkins University (JHU) and the Applied Physics Lab (APL) for the AMPTE program. ACE is designed as a spinning spacecraft with its spin axis aligned to the Earth-Sun axis. The ACE launch weight will be approx. 633 kg, including 105 kg of scientific instruments and 184 kg of propellant. Using a Delta-class expendable launch vehicle, ACE will be launched into an L1 libration point (240 R(sub e)) orbit. Telemetry will be 6.7 kbps average, using tape recorder storage with daily readout to DSN. The experiment power requirement is approximately 76 W nominal and 96 W peak. The prime objective of the ACE mission is: (1) to determine accurate elemental and isotropic abundances including solar matter, local interstellar matter and local galactic matter; (2) to study the origin of elements and evolutionary processing in galactic nucleosynthesis, galactic evolution, origin and evolution of the solar system; (3) to study coronal formation and solar-wind acceleration processes; and (4) to study particle acceleration and transport, including coronal shock acceleration, stochastic flare acceleration, interplanetary shock acceleration, and interstellar acceleration and propagation. To accomplish this objective, ACE will perform comprehensive and coordinated determinations of the elemental and isotopic composition of energetic nuclei accelerated on the Sun, in interplanetary space, and from galactic sources. These observations will span five decades in energy, from solar wind to galactic cosmic ray energies, and will cover the element range from H-1 to Zr-40. Comparison of these samples of matter will be used to study the origin and subsequent evolution of both solar system and galactic material by isolating the effects of fundamental processes that include nucleosynthesis, charged and neutral particle separation, bulk plasma acceleration, and the acceleration of suprathermal and high-energy particles.

Source record↗

235-F GoldSim Fate and Transport Model: Uncertainty Quantification

Building 235-F was configured with two missions in mind: Actinide Billet Line (ABL) and the fabrication of Pu-238 oxide for space program applications. ABL produced Np-237 billets for use in SRS reactors, whereas the design process, fabrication, and examination of Pu-238 oxide powder occurred within the following areas, respectively: Plutonium Experimental Facility (PEF), Plutonium Fuel Form (PuFF), and Old Metallography Lab (OML). By 1990 production ceased and by 2006 de-inventory occurred; however, assays have shown significant holdup remains within ABL and PuFF. As a result, 235-F is a Category 2 nuclear facility, with plans to undergo deactivation and decommission (D and D) via In-Situ Disposal (ISD). The purpose of this project is to ensure United States Environmental Protection Agency (USEPA) groundwater radiation maximum contaminant level (MCL) and dosage standards are met during the D and D of 235-F by quantifying uncertainty through probabilistic modeling and evaluation of various ISD alternatives. GoldSim is a dynamic modeling software package with a graphical, object-oriented interface capable of capturing the influence of complex system input variability on probabilistic system outcomes. A GoldSim stochastic fate and transport model for 235-F was developed and matched with a PORFLOW deterministic model to simulate probabilistic release and flow of radionuclides from ABL and PuFF into the vadose zone, the Upper Three Runs (UTR) Aquifer, and UTR Creek. The GoldSim model was used to probabilistically evaluate four ISD scenarios against USEPA groundwater radiation MCLs and dosage standards. The deterministic 235-F GoldSim fate and transport model continues to be refined to match the results of the PORFLOW deterministic model to ensure the model accurately represents radionuclide movement through the groundwater system. The stochastic variables that are utilized within the GoldSim model are founded on the most current data; a conservative perspective is taken where needed. Alignment with the PORFLOW deterministic model, coupled with input stochastic variability, allows the probabilistic 235-F GoldSim model to capture the conservative breadth of possible outcomes for radionuclide fate within this particular system. This ensures that the USEPA MCLs and dosage limits hold even in the worst case scenarios.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Overview of NASARTI (NASA Radiation Track Image) Program: Highlights of the Model Improvement and the New Results

This presentation summarizes several years of research done by the co-authors developing the NASARTI (NASA Radiation Track Image) program and supporting it with scientific data. The goal of the program is to support NASA mission to achieve a safe space travel for humans despite the perils of space radiation. The program focuses on selected topics in radiation biology that were deemed important throughout this period of time, both for the NASA human space flight program and to academic radiation research. Besides scientific support to develop strategies protecting humans against an exposure to deep space radiation during space missions, and understanding health effects from space radiation on astronauts, other important ramifications of the ionizing radiation were studied with the applicability to greater human needs: understanding the origins of cancer, the impact on human genome, and the application of computer technology to biological research addressing the health of general population. The models under NASARTI project include: the general properties of ionizing radiation, such as particular track structure, the effects of radiation on human DNA, visualization and the statistical properties of DSBs (DNA double-strand breaks), DNA damage and repair pathways models and cell phenotypes, chromosomal aberrations, microscopy data analysis and the application to human tissue damage and cancer models. The development of the GUI and the interactive website, as deliverables to NASA operations teams and tools for a broader research community, is discussed. Most recent findings in the area of chromosomal aberrations and the application of the stochastic track structure are also presented.

Ponomarev, Artem L.↗

Novel machine-learning method for spin classification of neutron resonances

The performance of nuclear reactors and other nuclear systems depends on a precise understanding of the neutron interaction cross sections for materials used in these systems. These cross sections exhibit resonant structure whose shape is determined in part by the angular-momentum quantum numbers of the resonances. The correct assignment of the quantum numbers of neutron resonances is, therefore, paramount. In this project, we apply machine learning to automate the quantum number assignments using only the resonances' energies and widths and not relying on detailed transmission or capture measurements. The classifier used for quantum number assignment is trained using stochastically generated resonance sequences whose distributions mimic those of real data. Here we explore the use of several physics-motivated features for training our classifier. These features amount to out-of-distribution tests of a given resonance's widths and resonance-pair spacings. We pay special attention to situations where either capture widths cannot be trusted for classification purposes or where there is insufficient information to classify resonances by the total spin J. We demonstrate the efficacy of our classification approach using simulated and actual 52 Cr resonance data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

GRACE Final Technical Report

The GRACE project (Grid that is Risk-Aware for Clean Electricity) was a five-year research initiative funded by the U.S. Department of Energy's Advanced Research Projects Agency- Energy (ARPA-E) under the PERFORM program. Led by Duke University's Nicholas School of the Environment, with contributions from The Ohio State University, North Carolina State University, Dartmouth College, and Pacific Northwest National Laboratory, the project addressed a fundamental challenge in grid management: conventional software plans for a single most-likely outcome and relies on reserves as a buffer, leaving utilities poorly equipped for the growing variability introduced by renewable energy. GRACE demonstrated a better approach: explicitly representing thousands of plausible future scenarios and choosing operating schedules that perform well across all of them.

14 SOLAR ENERGY↗

A stochastic optimal feedforward and feedback control methodology for superagility

A new control design methodology is developed: Stochastic Optimal Feedforward and Feedback Technology (SOFFT). Traditional design techniques optimize a single cost function (which expresses the design objectives) to obtain both the feedforward and feedback control laws. This approach places conflicting demands on the control law such as fast tracking versus noise atttenuation/disturbance rejection. In the SOFFT approach, two cost functions are defined. The feedforward control law is designed to optimize one cost function, the feedback optimizes the other. By separating the design objectives and decoupling the feedforward and feedback design processes, both objectives can be achieved fully. A new measure of command tracking performance, Z-plots, is also developed. By analyzing these plots at off-nominal conditions, the sensitivity or robustness of the system in tracking commands can be predicted. Z-plots provide an important tool for designing robust control systems. The Variable-Gain SOFFT methodology was used to design a flight control system for the F/A-18 aircraft. It is shown that SOFFT can be used to expand the operating regime and provide greater performance (flying/handling qualities) throughout the extended flight regime. This work was performed under the NASA SBIR program. ICS plans to market the software developed as a new module in its commercial CACSD software package: ACET.

Halyo, Nesim↗

Benchmark for Fuel Shuffling and Depletion for Pebble-Bed Reactors

Pebble bed reactors have specific operational characteristics when their fuel-cycle and fueling operations are considered. They are specifically distinguished by other type of nuclear reactor designs by their online fuel recycling scheme, where the fuel elements that have not yet reached discharge burnup can be reloaded and recycled continuously during normal operation. The fuel in a pebble bed reactor is not stationary and stochastically moves through the core once or several times during its lifetime, which allows them to operate without requiring a large excess reactivity hold for the burnup. However, this characteristic of pebble bed reactors introduces challenges in simulation, as each pebble can take many different trajectories through the core, its composition depends on the details of the irradiation history that is unique to its aggregated path through the core. For predicting the safety performance characteristics, such as source term, maximum fuel temperatures and fuel failure rates, etc., it is important to accurately incorporate the movement of pebbles through the core during their lifetime in a multi-physics simulation together with other phenomena. The equilibrium core analysis for pebble bed reactors are performed with multi-physics tools including fuel depletion in a multi pass reload coupled to the fuel movement. Currently, there are only a few legacy multi-physics simulation tools that can implement the pebble flow characteristics and perform equilibrium core analysis for pebble bed reactors. However, there are development efforts on-going under Department of Energy's Nuclear Energy Advanced Modelling and Simulation program and also in private industry for including these capabilities into their modelling and simulation tools. Any new development in the modelling and simulation tools needs to be validated by using tools such as experiments, analytical solutions or code-to-code benchmarks. In this work, a code-to-code benchmark for the equilibrium core analysis capability of pebble bed reactors was developed. Multiple cases were identified to capture different fuel cycle strategies that can be used in PBRs. The results of each case are presented in terms of overall equilibrium core characteristics: the discharge burnup; spatial burnup distribution; spatial isotopic distributions; axial and radial neutron flux distributions and power history of fuel elements per pass through the core for both a prototypical pebble bed High Temperature Gas-cooled Reactor and a prototypical pebble bed Fluoride-salt cooled High temperature Reactor.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗