Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nonlinear 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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Modeling of two-storey precast school building using Ruaumoko 2D program

The long-distant earthquake loading from Sumatra and Java Island had caused some slight damages to precast and reinforced concrete buildings in West Malaysia such as cracks on wall panels, columns and beams. Subsequently, the safety of existing precast concrete building is needed to be analyzed because these buildings were designed using BS 8110 which did not include the seismic loading in the design. Thus, this paper emphasizes on the seismic performance and dynamic behavior of precast school building constructed in Malaysia under three selected past earthquakes excitations ; El Centro 1940 North-South, El Centro East-West components and San Fernando 1971 using RUAUMOKO 2D program. This program is fully utilized by using prototype precast school model and dynamic non-linear time history analysis. From the results, it can be concluded that two-storey precast school building has experienced severe damage and partial collapse especially at beam-column joint under San Fernando and El Centro North-South Earthquake as its exceeds the allowable inter-storey drift and displacement as specified in Eurocode 8. The San Fernando earthquake has produced a massive destruction to the precast building under viscous damping, ξ = 5% and this building has generated maximum building displacement of 435mm, maximum building drift of 0.68% and maximum bending moment at 8458kNm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Center for Extended Magnetohydrodynamic Modeling (CEMM), Final Technical Report for UW-Madison Engineering Physics

The Center for Extended Magnetohydrodynamic Modeling (CEMM) focused on the further development, verification, validation, and application of two scientific computation programs, the NIMROD and M3D-C1 codes, for modeling macroscopic dynamics in magnetized plasma. Both NIMROD and M3D-C1 solve the mathematical equations that model two-fluid 3D magnetohydrodynamic (MHD). They are designed and optimized to be efficient and accurate for toroidal confinement physics calculations. The effort improved the fidelity of the physics models, the realism of the treatment of the surrounding vessel, coils, and other structures, the accuracy and efficiency of the algorithms, and the ability to scale calculations to thousands of processors in order to take advantage of leading-edge computers. The primary application areas of Resonant Magnetic Perturbations (RMP), Edge Harmonic Oscillations (EHO), Edge Localized Modes (ELMs), Disruption Studies, Sawteeth and Stationary States, and Kinetic MHD and Neoclassical Tearing Modes (NTM) are of concern for the international ITER experiment, which is under construction in France. The effort by the UW-Madison Engineering Physics component of CEMM improved the numerical methods in the NIMROD code with respect to modeling MHD interchange and two-fluid effects. It also improved the algorithm for solving large algebraic systems to enhance NIMROD's computational performance when running the code in parallel on increasing numbers of processor cores. The modeling efforts included a study of the nonlinear evolution of ballooning, and the findings were supported by analytical theory and computation. Numerical MHD modeling also examined the influences of two-fluid effects on ELMs and RMPs, where "two-fluid" refers to separation of electron dynamics from ion dynamics. Another related aspect considered the effects of a localized layer of electrical current density near the edge of the plasma confinement region. Two-fluid effects on resonant reconnecting instabilities were verified and applied to computations of sawtooth and drift-tearing dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Process Simulation of Cold Pressing and Sintering of Armstrong CP-Ti Powders

A computational methodology is presented for the process simulation of cold pressing and sintering of Armstrong CP-Ti powders. Since the powder consolidation is governed by specific pressure-dependent constitutive equations, solution algorithms were developed for the ABAQUS user material subroutine, UMAT, for computing the plastic strain increments based on an implicit integration of the nonlinear yield function, flow rule, and hardening equations. Sintering was simulated using a model based on diffusional creep using the user subroutine CREEP. The initial mesh, stress, and density for the simulation of sintering were obtained from the results of the cold pressing simulation, minimizing the errors from decoupling the cold pressing and sintering simulations. Numerical simulation results are presented for the cold compaction followed by a sintering step of the Ti powders. The numerical simulation results for the relative density were compared to those measured from experiments before and after sintering, showing that the relative density can be accurately predicted. Notice: This manuscript has been authored by UT-Battelle, LLC, under Contract No. DE-AC05-00OR22725 with the U.S. Department of Energy. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, world-wide license to publish or reproduce the published form of this manuscript, or allow others to do so, for United States Government purposes. ACKNOWLEDGEMENTS This research was sponsored by the U.S. DOE, and carried out at ORNL, under Contract DE-AC05-00OR22725 with UT-Battelle, LLC. This research was sponsored by the U.S. DOE, EERE Industrial Technology Program Office under CPS Agreement # 17881.

36 MATERIALS SCIENCE↗
DOE OSTI2018

AOI [3]: Smart Refractory Sensor Systems for Wireless Monitoring of Temperature, Health, and Degradation of Slagging Gasifiers

The objective of the work was to develop refractory “smart bricks”, which would contain embedded temperature, strain/stress, and spallation sensors throughout the volume of highchromia (-Cr 2 O 3 ) refractory brick. The proposed work included work to interconnect the sensors to the reactor exterior, where the sensor signals may be processed by low-power electronics and transmitted wirelessly to a central processing hub. The data processing and wireless transmitter hardware was specifically designed to be isolated (with low power consumption) and to be adaptable to future implementation of energy-harvesting strategies for extended life. Finally, the collected data was incorporated into a model to estimate refractory degradation, a technique that could help monitor the health of the refractory in real-time. The long-term goal of this program was to demonstrate high-temperature, wireless sensor arrays for in situ three-dimensional (3-D) refractory monitoring or mapping for slagging gasification systems. The research was in collaboration with HarbisonWalker International (HWI) Technology Center in West Mifflin, PA. HWI is a leading developer and manufacturer of ceramic refractory products for high-temperature applications. The work completed focused on the following areas: 1) Investigation of the chemical stability, microstructural evolution, grain growth kinetics, degree of homogeneity (quantitative image analysis), and electrical properties of refractory oxide-silicide composites at temperatures between 750-1450ºC; 2) Fabrication of silicide-alumina composite and oxide thermocouples and thermistor preforms and the development of techniques to embed them into high-chromia refractory bricks to form “smart bricks”; 3) Utilization of commercial off-the-shelf discrete components to prototype circuits for interfacing between smart brick sensors and the wireless sensor network. The prototypes were then used to design an integrated circuit for thermistor, thermocouple, and capacitive-based smart brick sensor interfacing; 4) Interfacing of the smart bricks with embedded sensors with wireless motes thus yielding a complete signal chain. This end-to-end data collection system was tested on a furnace heated to 1350°C; 5) Development of a slag penetration model and a nonlinear unknown input filter for the data from the embedded sensors for estimating temperature and extent of slag penetration.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Approximating Trajectory Constraints With Machine Learning – Microgrid Islanding With Frequency Constraints

Here, we introduce deep earning aided constraint encoding to tackle the frequency-constraint microgrid scheduling problem. The nonlinear function between system operating condition and frequency nadir is approximated by using a neural network, which admits an exact mixed-integer formulation (MIP). This formulation is then integrated with the scheduling problem to encode the frequency constraint. With the stronger representation power of the neural network, the resulting commands can ensure adequate frequency response in a realistic setting in addition to islanding success. The proposed method is validated on a modified 33-node system. Successful islanding with a secure response is simulated under the scheduled commands using a detailed three-phase model in Simulink. The advantages of our model are particularly remarkable when the inertia emulation functions from wind turbine generators are considered.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2015

Light Water Reactor Sustainability Program Advanced Seismic Soil Structure Modeling

Risk calculations should focus on providing best estimate results, and associated insights, for evaluation and decision-making. Specifically, seismic probabilistic risk assessments (SPRAs) are intended to provide best estimates of the various combinations of structural and equipment failures that can lead to a seismic induced core damage event. However, in some instances the current SPRA approach has large uncertainties, and potentially masks other important events (for instance, it was not the seismic motions that caused the Fukushima core melt events, but the tsunami ingress into the facility). SPRA’s are performed by convolving the seismic hazard (this is the estimate of all likely damaging earthquakes at the site of interest) with the seismic fragility (the conditional probability of failure of a structure, system, or component given the occurrence of earthquake ground motion). In this calculation, there are three main pieces to seismic risk quantification, 1) seismic hazard and nuclear power plants (NPPs) response to the hazard, 2) fragility or capacity of structures, systems and components (SSC), and 3) systems analysis. Two areas where NLSSI effects may be important in SPRA calculations are, 1) when calculating in-structure response at the area of interest, and 2) calculation of seismic fragilities (current fragility calculations assume a lognormal distribution for probability of failure of components). Some important effects when using NLSSI in the SPRA calculation process include, 1) gapping and sliding, 2) inclined seismic waves coupled with gapping and sliding of foundations atop soil, 3) inclined seismic waves coupled with gapping and sliding of deeply embedded structures, 4) soil dilatancy, 5) soil liquefaction, 6) surface waves, 7) buoyancy, 8) concrete cracking and 9) seismic isolation The focus of the research task presented here-in is on implementation of NLSSI into the SPRA calculation process when calculating in-structure response at the area of interest. The specific nonlinear soil behavior included in the NLSSI calculation presented in this report is gapping and sliding. Other NLSSI effects are not included in the calculation. The results presented in this report document initial model runs in the linear and nonlinear analysis process. Final comparisons between traditional and advanced SPRA will be presented in the September 30th deliverable.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Size effects on thermo-mechanical performance of U-10Mo monolithic fuel plates

Monolithic fuel is a fuel form that is considered for the conversion of high performance research reactors. This plate-type fuel consists of a high density U-Mo fuel in monolithic form that is sandwiched between zirconium diffusion barriers, and encapsulated in an aluminum cladding. To date, large number of plates have been irradiated with satisfactory perforamce. The program is now moving into the qualification phase, a predecessor to the timely conversion of the target reactors. It must be shown that the fuel system meets the safety standards and performs well in reactor. The requirement to satisfactory irradiation performance under normal operating conditions is primarily demonstrated by a successful testing. Since each reactor employs distinct fuel plate geometries for various consideration with unique plate design features and attributes, a single “generic” plate geometry capturing all of the extremities is not achievable. Furthermore, testing all these geometric and irradiation parameters on a large size plate is not practical. Therefore, a smaller, “down-scaled” versions of fuel plates, are often employed for experimental purposes. This limitation consequently requires much more cautious performance evaluations, as thermal and mechanical response of a plate with certain geometry may not be representative for a plate with a different geometry. To investigate if plate size has any effects on irradiation performance, the plates with various geometric dimensions were parametrically evaluated. In particular, length and width of the plates were varied between the bounding values. Temperature, deformation, stress values were comparatively evaluated. The results have indicated that effects of geometric ratios and plate size variations in length and width directions are insignificant. However, wider plates could become more prone to a warping-type deformation, if there are nonlinearities.

42 ENGINEERING↗
DOE OSTI2019

Size effects on thermo-mechanical performance of U-10Mo monolithic fuel plates

Monolithic fuel is a fuel form that is considered for the conversion of high performance research reactors. This plate-type fuel consists of a high density U-Mo fuel in monolithic form that is sandwiched between zirconium diffusion barriers, and encapsulated in an aluminum cladding. To date, large number of plates have been irradiated with satisfactory perforamce. The program is now moving into the qualification phase, a predecessor to the timely conversion of the target reactors. It must be shown that the fuel system meets the safety standards and performs well in reactor. The requirement to satisfactory irradiation performance under normal operating conditions is primarily demonstrated by a successful testing. Since each reactor employs distinct fuel plate geometries for various consideration with unique plate design features and attributes, a single “generic” plate geometry capturing all of the extremities is not achievable. Furthermore, testing all these geometric and irradiation parameters on a large size plate is not practical. Therefore, a smaller, “down-scaled” versions of fuel plates, are often employed for experimental purposes. This limitation consequently requires much more cautious performance evaluations, as thermal and mechanical response of a plate with certain geometry may not be representative for a plate with a different geometry. To investigate if plate size has any effects on irradiation performance, the plates with various geometric dimensions were parametrically evaluated. In particular, length and width of the plates were varied between the bounding values. Temperature, deformation, stress values were comparatively evaluated. The results have indicated that effects of geometric ratios and plate size variations in length and width directions are insignificant. However, wider plates could become more prone to a warping-type deformation, if there are nonlinearities.

42 ENGINEERING↗
DOE OSTI2023

MODELING A SPENT NUCLEAR FUEL CASK SEISMIC TEST

The US Department of Energy Spent Fuel and Waste Science and Technology (SFWST) program is planning to conduct a series of full-scale shake table tests to simulate hypothetical earthquake conditions and record the response of surrogate spent nuclear fuel (SNF) assemblies in a canister dry storage system mockup. The shake table motions will represent a range of hypothetical earthquake conditions at hypothetical locations in the continental US to generally define the range of mechanical loads that SNF can be expected to experience during extended dry storage periods. The earthquake conditions will represent seismic hazards in the 2,000-to-20,000-year return period range. The test will use instrumented pressurized water reactor fuel assemblies with surrogate mass inside the fuel rods instead of radioactive fuel pellets. Dummy assemblies with similar mass and dimensions to actual fuel assemblies will occupy the rest of the fuel assembly locations in a SNF canister that holds 32 fuel assemblies. The canister will be located inside a mockup vertical concrete cask. Instrumentation will record the motion of the major components of the complex dynamic system, and strain gauges will be used to record the cladding strain at select locations. Preparations for the test require modeling predictions to identify the range of response of the system and to help select specific earthquake cases to be simulated on the shake table from a large set of potential cases. This paper describes the pretest nonlinear finite element modeling efforts that have been completed to date, including cask system level modeling and fuel assembly modeling in LS-DYNA. The cask system level models are critical for anticipating sliding or tipping of an unanchored cask during the test. The fuel assembly model is needed to estimate the range of cladding strain response and fuel assembly structural response to be expected during the test and cladding strain measurements will be key metrics for model validation and the development of modeling best practices after the test is concluded.

Klymyshyn, Nicholas A.↗
DOE OSTI2010

Presto 4.18 user's guide.

Presto is a Lagrangian, three-dimensional explicit, transient dynamics code that is used to analyze solids subjected to large, suddenly applied loads. The code is designed for a parallel computing environment and for problems with large deformations, nonlinear material behavior, and contact. Presto also has a versatile element library that incorporates both continuum elements and structural elements. This user's guide describes the input for Presto that gives users access to all the current functionality in the code. The environment in which Presto is built allows it to be coupled with other engineering analysis codes. Using a concept called scope, the input structure reflects the fact that Presto can be used in a coupled environment. The user's guide describes how scope is implemented from the outermost to the innermost scopes. Within a given scope, the descriptions of input commands are grouped based on functionality of the code. For example, all material input command lines are described in a chapter of the user's guide for all the material models that can be used in Presto.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2018

Load Mitigation on Floating Offshore Wind Turbines with Advanced Controls and Tuned Mass Dampers

General Electric, the National Renewable Energy Laboratory (NREL), the University of Massachusetts Amherst (UMass), and Glosten have recently completed a US Department of Energy (DOE)-funded research program to study technologies for mitigating loads on floating offshore wind turbines through the use of advanced turbine controls and tuned mass dampers (TMDs). The analysis was based upon the Glosten PelaStar tension leg platform (TLP) with GE Haliade 150 turbine, a system developed in a previous front end engineering design (FEED) study funded by the Energy Technology Institute (ETI) in the UK. The platform was designed for the WaveHub wave energy research site, with a mean water depth of 59-m. Loads were analyzed by running time-domain simulations in four 50-year return period (50-YRP) ultimate load state (ULS) conditions and 77 fatigue load state (FLS) environmental conditions. In 50-YRP conditions advanced controls are not active. The influence of TMDs on ULS loads have been reported previously (Park et al. [2]). In FLS conditions advanced controls and TMDs afford dramatic reductions in fatigue damage, offering the potential of significant savings in tower structural requirements. Simulations in turbine idling conditions were run in OrcaFlex, and simulations in operating conditions were run in FAST. Simulations were run with a baseline turbine controller, representative of the current state of the art, and an advanced controller developed by NREL to use collective and individual blade pitch control to maintain rotor speed and reduce tower loads. UMass developed a number of TMD types, with varying system configurations, including passive nonlinear dampers and semi-actively controlled dampers with an inverse velocity groundhook control algorithm. Loads and accelerations in FLS conditions were evaluated on the basis of damage equivalent loads (DELs), and fatigue damage was computed by Miner's summations of stress cycles at the tower base. To study sensitivity to water depth, loads were analyzed at both the 59-m WaveHub depth and a more commercially realistic depth of 100 m. TMDs reduce fatigue damage at the tower-column interface flange by up to 52% in 59-m water depth and up to 28% in 100 m water depth. Advanced controls reduce fatigue damage at the tower-column flange by up to 22% in 59-m water depth and up to 40% in 100 m water depth. The most effective load-mitigation strategy is combining advanced controls with TMDs. This strategy affords a 71% reduction in fatigue damage in both 59-m and 100-m water depths.

17 WIND ENERGY↗
DOE OSTI2019

Redox transistors for neuromorphic computing

Efficiency bottlenecks inherent to conventional computing in executing neural algorithms have spurred the development of novel devices capable of “in-memory” computing. Commonly known as “memristors,” a variety of device concepts including conducting bridge, vacancy filament, phase change, and other types have been proposed as promising elements in artificial neural networks for executing inference and learning algorithms. In this article, we review the recent advances in memristor technology for neuromorphic computing and discuss strategies for addressing the most significant performance challenges, including nonlinearity, high read/write currents, and endurance. As an alternative to two-terminal memristors, we introduce the three-terminal electrochemical memory based on the redox transistor (RT), which uses a gate to tune the redox state of the channel. Decoupling the “read” and “write” operations using a third terminal and storage of information as a charge-compensated redox reaction in the bulk of the transistor enables high-density information storage. These properties enable low-energy operation without compromising analog performance and nonvolatility. Finally, we discuss the RT operating mechanisms using organic and inorganic materials, approaches for array integration, and prospects for achieving the device density and switching speeds necessary to make electrochemical memory competitive with established digital technology.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Novel, Middle and Long Wave Infrared Laser Sources For Accelerator Applications

The main accomplishment during Years 2020-2024 of the DOE Accelerator Stewardship Program is proof-of-principle demonstration of amplification of an ultrashort 10 µm seed pulse in a multiatmosphere CO 2 laser optically pumped by nanosecond pulses of either a ~4.3 µm Fe:ZnSe laser or a 2.8 µm electro-optically (EO) Q-switched Er: Laser. This has been accomplished at the UCLA Neptune Laboratory by using about 100 mJ, ~100 ns long pulses generated by the pump laser built for this purpose by the UAB team. Implementation of EO Q-switching regime in the solid-state Er:laser system made possible nanosecond synchronization between the pump pulse and the 0.5 ps seed pulse produced by a commercial Ti:sapphire laser based mid-IR OPA/DFG. The injection mode-locking regime realized in a compact Mult atmosphere 3 Hz CO 2 -He regenerative amplifier enabled generation of a train of ultrashort pulses separated by the cavity round-trip time of ~3.2 ns. The measured gain lifetime of ~1 µs at a total pressure of ~12-15 atm indicates that lifetime of the upper laser level is a critical limiting factor in designing high-pressure CO 2 amplifiers. A single picosecond 10 mm pulse with an energy 0.1-0.2 mJ was extracted from the amplifier’s cavity using an ultrafast Ge semiconductor switch whose reflectivity was photo-induced by a femtosecond 800 nm pulse. For direct amplification regime, nonlinear absorption in Ge at 10 mm at GW/cm 2 intensities limited the extracted energy. Therefore, the generated energy could be further increased by minimizing losses in the switch and optimizing the cavity mode size. Based on the obtained data in collaboration with BNL colleagues, we have designed a high-pressure CO 2 amplifier to be used as a preamplifier for the future experiments at the ATF BNL. The latter work is the subject of funding availability.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2010

Presto 4.16 user s guide.

Presto is a three-dimensional transient dynamics code with a versatile element library, nonlinear material models, large deformation capabilities, and contact. It is built on the SIERRA Framework [1, 2]. SIERRA provides a data management framework in a parallel computing environment that allows the addition of capabilities in a modular fashion. Contact capabilities are parallel and scalable. The Presto 4.16 User's Guide provides information about the functionality in Presto and the command structure required to access this functionality in a user input file. This document is divided into chapters based primarily on functionality. For example, the command structure related to the use of various element types is grouped in one chapter; descriptions of material models are grouped in another chapter. The input and usage of Presto is similar to that of the code Adagio [3]. Adagio is a three-dimensional quasi-static code with a versatile element library, nonlinear material models, large deformation capabilities, and contact. Adagio, like Presto, is built on the SIERRA Framework [1]. Contact capabilities for Adagio are also parallel and scalable. A significant feature of Adagio is that it offers a multilevel, nonlinear iterative solver. Because of the similarities in input and usage between Presto and Adagio, the user's guides for the two codes are structured in the same manner and share common material. (Once you have mastered the input structure for one code, it will be easy to master the syntax structure for the other code.) To maintain the commonality between the two user's guides, we have used a variety of techniques. For example, references to Adagio may be found in the Presto user's guide and vice versa, and the chapter order across the two guides is the same. On the other hand, each of the two user's guides is expressly tailored to the features of the specific code and documents the particular functionality for that code. For example, though both Presto and Adagio have contact functionality, the content of the chapter on contact in the two guides differs. Important references for both Adagio and Presto are given in the references section at the end of this chapter. Adagio was preceded by the codes JAC and JAS3D; JAC is described in Reference 4; JAS3D is described in Reference 5. Presto was preceded by the code Pronto3D. Pronto3D is described in References 6 and 7. Some of the fundamental nonlinear technology used by both Presto and Adagio are described in References 8, 9, and 10. Currently, both Presto and Adagio use the Exodus II database and the XDMF database; Exodus II is more commonly used than XDMF. (Other options may be added in the future.) The Exodus II database format is described in Reference 11, and the XDMF database format is described in Reference 12. Important information about contact is provided in the reference document for ACME [13]. ACME is a third-party library for contact. One of the key concepts for the command structure in the input file is a concept referred to as scope. A detailed explanation of scope is provided in Section 1.2. Most of the command lines in Chapter 2 are related to a certain scope rather than to some particular functionality.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2010

Adagio 4.16 user s guide.

Adagio is a three-dimensional, implicit solid mechanics code with a versatile element library, nonlinear material models, and capabilities for modeling large deformation and contact. Adagio is a parallel code, and its nonlinear solver and contact capabilities enable scalable solutions of large problems. It is built on the SIERRA Framework [1, 2]. SIERRA provides a data management framework in a parallel computing environment that allows the addition of capabilities in a modular fashion. The Adagio 4.16 User's Guide provides information about the functionality in Adagio and the command structure required to access this functionality in a user input file. This document is divided into chapters based primarily on functionality. For example, the command structure related to the use of various element types is grouped in one chapter; descriptions of material models are grouped in another chapter. The input and usage of Adagio is similar to that of the code Presto [3]. Presto, like Adagio, is a solid mechanics code built on the SIERRA Framework. The primary difference between the two codes is that Presto uses explicit time integration for transient dynamics analysis, whereas Adagio is an implicit code. Because of the similarities in input and usage between Adagio and Presto, the user's guides for the two codes are structured in the same manner and share common material. (Once you have mastered the input structure for one code, it will be easy to master the syntax structure for the other code.) To maintain the commonality between the two user's guides, we have used a variety of techniques. For example, references to Presto may be found in the Adagio user's guide and vice versa, and the chapter order across the two guides is the same. On the other hand, each of the two user's guides is expressly tailored to the features of the specific code and documents the particular functionality for that code. For example, though both Presto and Adagio have contact functionality, the content of the chapter on contact in the two guides differs. Important references for both Adagio and Presto are given in the references section at the end of this chapter. Adagio was preceded by the codes JAC and JAS3D; JAC is described in Reference 4; JAS3D is described in Reference 5. Presto was preceded by the code Pronto3D. Pronto3D is described in References 6 and 7. Some of the fundamental nonlinear technology used by both Presto and Adagio are described in References 8, 9, and 10. Currently, both Presto and Adagio use the Exodus II database and the XDMF database; Exodus II is more commonly used than XDMF. (Other options may be added in the future.) The Exodus II database format is described in Reference 11, and the XDMF database format is described in Reference 12. Important information about contact is provided in the reference document for ACME [13]. ACME is a third-party library for contact. One of the key concepts for the command structure in the input file is a concept referred to as scope. A detailed explanation of scope is provided in Section 1.2. Most of the command lines in Chapter 2 are related to a certain scope rather than to some particular functionality.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Bayesian Monte Carlo Evaluation of Imperfect (n, 233 U) Data and Model

Conventional nuclear data evaluation methods using generalized linear least squares make the following assumptions: prior and posterior probability distribution functions (PDFs) of all model parameters and data are normal (Gaussian); the linear approximation is sufficiently accurate to minimize the cost function (even for nonlinear models); the model (e.g., of neutron cross section) and experimental data (including covariance data) are without defect and prior PDFs of parameters and measured data are known perfectly. Neglect of covariance between model parameters and measured data in conventional evaluations contributes to imperfections. These assumptions are inherent to the generalized linear least squares minimization method commonly used for resolved resonance region neutron cross section evaluations but are often not justified due to the presence of non-normal PDFs, nonlinear models (e.g., R-matrix formalism), and inherent imperfections in data and models (e.g., imperfect covariance data). Here, these assumptions are removed in a mathematical framework of Bayes’ theorem, which is implemented using the Metropolis-Hastings Monte Carlo method. Most importantly, new parameters are introduced to parameterize discrepancies between the theoretical model and measured data to quantify judgement about discrepancies or imperfections in a reproducible manner. An evaluation of 233U in the eV region using the ENDF-B/VIII.0 library and transmission data (Guber et al.) is presented, and posterior parameters are compared to those obtained by conventional evaluation methods. This example illustrates the effects of removing the most harmful assumption: that of model-data perfection.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

A Hybrid AI/ML and Computational Mechanics Based Approach for Time-Series State and Fatigue Life Estimation of Nuclear Reactor Components

Environmental fatigue modeling is a complex problem due to multiple failure modes and their intermixing. The failure modes are function of various underlying causes in addition to the corrosive effect of reactor coolant environment. Some of the major causes are time-dependence of material associated with cyclic loading, load sequence effect associated with random/variable amplitude loading, effect of strain amplitude and rates, effect of varying temperature (along both temporal and spatial directions) and the effect of mean strain and stress. The nonlinear intermixing of failure modes associated with above mentioned causing parameters makes the environmental fatigue modeling is a challenging task. Because of this challenge, fatigue is traditionally being modeled based on experimental data. However, test based empirical approach often requires hundreds of fatigue tests to model the above-mentioned intermixing failure causes even for a single material system. The problem is further exaggerated for reactor component made from multi-material systems such as made from both carbon and stainless-steel base metals and their similar and dissimilar metal welds. With the difficulty of conducting hundreds of fatigue tests to capture the above-mentioned intermixing failure causes, fatigue modeling approaches often depends on empirical models based on limited available test data such as available through ASME code and NUREG 6909. However, these limited test-data-based models may not be enough to accurately predict the life of reactor components. Accurate prediction of life of reactor component would become a necessity, particularly when the license of the reactors to be extended for long-term-operation (LTO) that is for well beyond its original design life of 40 years. The requirement of extending the license of reactor under LTO requires hundreds of fatigue tests to be conducted to understand the mechanism associated with the above-mentioned interdependent failure causes. However, conducting large number of fatigue tests is not a feasibility due to the cost involved. To address this issues Argonne National Laboratory (ANL) with the sponsorship of DOE Light Water Reactor Sustainability (LWRS) program trying to develop a hybrid predictive modeling approach. This is based on limited experiment-data, Artificial-intelligence (AI) – Machine-Learning (ML) - Deep-Learning (DL) based techniques and Multiphysics-computational-mechanics based modeling tools. The hybrid approach not-only can improve the accuracy of the existing stress analysis and fatigue modeling approach but also can reduce the over-dependency on test-based approach. Towards this goal following are some of the major contributions based on ANL’s FY-20 environmental fatigue modeling activities: 1) A cyclic plasticity material model database for 82/182 dissimilar metal weld, which can be readily shared with US nuclear industry and regulatory agency on request. 2) A well validated analytical modeling methodology to perform cycle-by-cycle stress prediction under both constant amplitude fatigue loading and variable amplitude fatigue loading (with load-sequence effect). 3) An AI/ML/DL based methodology to predict unmeasurable cyclic strain based on other available sensor signals. This type of approach can be used for estimating strain in real reactor components from other sensor readings. 4) An AI/ML based approach to improve the US capability on environmental fatigue testing. This is by improving ANL’s existing environmental fatigue testing capacity to conduct ASME required strain-controlled tests (by controlling strain amplitudes and its rate), while not measuring the strain (due to the difficulty of placing an extensometer in a narrow autoclave in a PWR-water-test system). 5) A simulation and experiment based probabilistic modeling methodology for time-series fatigue state and life estimation of reactor metal such as dissimilar metal weld.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2011

Near Field Modeling Of Spe1 Experiment And Prediction Of The Second Source Physics Experiments (SPE2)

Motion along joints and fractures in the rock has been proposed as one of the sources of near-source shear wave generation, and demonstrating the validity of this hypothesis is a focal scientific objective of the source physics experimental campaign in the Climax Stock granitic outcrop. A modeling effort has been undertaken by LLNL to complement the experimental campaign, and over the long term provide a validated computation capability for the nuclear explosion monitoring community. The approach involves performing the near-field nonlinear modeling with hydrodynamic codes (e.g., GEODYN, GEODYN-L), and the far-field seismic propagation with an elastic wave propagation code (e.g., WPP). the codes will be coupled together to provide a comprehensive source-to-sensor modeling capability. The technical approach involves pre-test predictions of each of the SPE experiments using their state of the art modeling capabilities, followed by code improvements to alleviate deficiencies identified in the pre-test predictions. This spiral development cycle wherein simulations are used to guide experimental design and the data from the experiment used to improve the models is the most effective approach to enable a transition from the descriptive phenomenological models in current use to the predictive, hybrid physics models needed for a science-based modeling capability for nuclear explosion monitoring. The objective of this report is to describe initial results of non-linear motion predictions of the first two SPE shots in the Climax Stock: a 220-lb shot at a depth of 180 ft (SPE No.1), and a 2570-lb shot at a depth of 150 ft (SPE No.2). The simulations were performed using the LLNL ensemble granite model, a model developed to match velocity and displacement attenuation from HARDHAT, PILE DRIVER, and SHOAL, as well as Russian and French nuclear test data in granitic rocks. This model represents the state of the art modeling capabilities as they existed when the SPE campaign was launched in 2010, and the simulation results presented here will establish a baseline that will be used for gauging progress as planned modeling improvements are implemented during the remainder of the SPE program. The initial simulations were performed under 2D axisymmetric conditions assuming the geologic medium to be a homogeneous half space. However, logging data obtained from the emplacement hole reveal two major faults that intersect the borehole at two different depth intervals (NSTec report, 2011) and four major joint sets. To evaluate the effect of these discrete structures on the wave forms generated they have performed 2D and 3D analysis with a Lagrangian hydrocode, GEODYN-L that shares the same material models with GEODYN but can explicitly take joints and fault into consideration. They discuss results obtained using these two different approaches in this report.

58 GEOSCIENCES↗