Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Static code analysis”

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 109 records · Page 6

Gravity at the horizon: on relativistic effects, CMB-LSS correlations and ultra-large scales in Horndeski's theory

We address the impact of consistent modifications of gravity on the largest observable scales, focusing on relativistic effects in galaxy number counts and the cross-correlation between the matter large scale structure (LSS) distribution and the cosmic microwave background (CMB). Our analysis applies to a very broad class of general scalar-tensor theories encoded in the Horndeski Lagrangian and is fully consistent on linear scales, retaining the full dynamics of the scalar field and not assuming quasi-static evolution. As particular examples we consider self-accelerating Covariant Galileons, Brans-Dicke theory and parameterizations based on the effective field theory of dark energy, using the hi class code to address the impact of these models on relativistic corrections to LSS observables. We find that especially effects which involve integrals along the line of sight (lensing convergence, time delay and the integrated Sachs-Wolfe effect—ISW) can be considerably modified, and even lead to O(1000%) deviations from General Relativity in the case of the ISW effect for Galileon models, for which standard probes such as the growth function only vary by O(10%). These effects become dominant when correlating galaxy number counts at different redshifts and can lead to ∼ 50% deviations in the total signal that might be observable by future LSS surveys. Because of their integrated nature, these deep-redshift cross-correlations are sensitive to modifications of gravity even when probing eras much before dark energy domination. We further isolate the ISW effect using the cross-correlation between LSS and CMB temperature anisotropies and use current data to further constrain Horndeski models. Forthcoming large-volume galaxy surveys using multiple-tracers will search for all these effects, opening a new window to probe gravity and cosmic acceleration at the largest scales available in our universe.

79 ASTRONOMY AND ASTROPHYSICS↗

Evolutionary Drift Models for Moving Target Defense

One of the biggest challenges faced by cyber defenders is that attacks evolve more rapidly than our ability to recognize them. We propose a moving target defense concept in which the means of detection is set in motion. This is done by moving away from static signature-based detection and instead adopting biological modeling techniques that describe families of related sequences. We present here one example for how to apply evolutionary models to cyber sequences, and demonstrate the feasibility of this technique on analysis of a complex, evolving software project. Specifically, we applied sequence-based and profile-based evolutionary models and report the ability of these models to recognize highly volatile code regions. We found that different drift models reliably identify different types of evolutionarily related code regions. The impact is that these (and possibly other) evolutionary models could be used in a moving target defense in which the "signature" being used to detect sequence-based behaviors is not a fixed signature but one that can recognize new variants of a known family based on multiple evolutionary models.

Oehmen, Christopher S.↗

Integrated FHR technology development: Tritium management, materials testing, salt chemistry control, thermal hydraulics and neutronics, associated benchmarking and commercial basis

This report summarizes the results of a 3-year Integrated Research Project (IRP) sponsored by the U.S. Department of Energy to address major technical challenges in the development of the Fluoride-salt-cooled High-temperature Reactor (FHR). The IRP universities included the Massachusetts Institute of Technology (MIT), the University of California at Berkeley (UCB), the University of Wisconsin (UW), and the University of New Mexico (UNM). The FHR is a new reactor concept less than 20 years old that combines (1) a clean fluoride salt coolant, (2) the graphite-matrix coated-particle fuel originally developed for High- Temperature Gas-Cooled Reactors and (3) passive decay heat removal systems from sodium fast reactors. The base-line design is a pebble-bed FHR that delivers heat between 600 to 700°C with a 7 Li 2 BeF 4 (flibe) salt coolant to the power cycle or industrial users. The Executive Summary describes what was accomplished and summarizes major conclusions. The main report contains stand-alone sections in each major area of work. Highlights are summarized below. Tritium. The coolant salt produces tritium necessitating understanding and control of tritium. Experimental work included measuring hydrogen uptake on carbon, tritium uptake during irradiation of graphite in 700°C salt and the permeability of tritium through metals with and without coatings. An FHR systems model was developed to predict tritium behavior and enable design of tritium control systems. The leading candidate for tritium removal from the flibe salt is carbon bed outside the reactor core. Thermal hydraulics and Neutronics. A 10-meter high FHR integral effects test facility, using Dowtherm A®, has been built and used to simulate the integral response of a representative FHR to improve understanding of FHR thermal hydraulics and validate thermal hydraulics codes used for the FHR. Neutronic analysis included code-to-code benchmarking, multiphysics modeling, uncertainty analysis for nuclear data uncertainties and shielding analysis. Twisted-tube heat exchangers were examined that may provide significant cost reductions because of their improved heat-transfer performance with liquid salts. The impacts on safety and design caused by uncertainties in physical properties of the salt were assessed with recommendations on where added physical property measurements are needed. Work began on understanding radiative heat transfer that becomes important in salt systems at higher temperatures. Materials. Static corrosion tests were conducted on graphite and materials in 700°C salt -- both in the UW laboratory and in the MIT reactor. In addition to static corrosion testing, a circulating salt corrosion loop was built and is operating. Methods to measure salt properties including redox were developed. The presence of carbon increases corrosion. With good redox control, the evidence indicates stainless steel can be used as a material of construction. Benchmarking. An additional series of workshops where held to enable benchmarking and comparison of results between different investigators. Commercialization. The FHR delivers higher average-temperature (~650°C) heat to the power cycle or industry than other reactors. This enables power cycles with lower-cost heat storage for variable electricity to the grid while the reactor operates at base-load -- providing a major competitive advantage by maximizing revenue from electricity sales. In October, 2016, a new startup company named Kairos Power was formed, and opened its first office in Jack London Square, Oakland, California, in January 2017. Kairos Power exists largely due to the work that was accomplished by this IRP and the previous DOE-funded FHR IRP that performed research on FHR technology from 2012 to 2015. In September, 2018, as this IRP is concluding, Kairos Power was staffed at 60 full time employees, working to develop technology derived from this IRP. Many of these employees worked with the IRP before joining Kairos Power 1 (https://kairospower.com). Other students and researchers from this IRP have joined other startup companies, including the Terrapower Molten Chloride Fast Reactor and the Terrestrial Energy development efforts.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Performance Exploration Through Optimistic Static Program Annotations

Compilers are limited by the static information directly or indirectly encoded in the program. Low-level languages, such as C/C++, are considered problematic as their weak type system and relaxed memory semantic allows for various, sometimes non-obvious, behaviors. Since compilers have to preserve the program semantics for all program executions, the existence of exceptional behavior can prevent optimizations that the developer would consider valid and might expect. Analyses to guarantee the absence of disruptive and unlikely situations are consequently an indispensable part of an optimizing compiler. However, such analyses have to be approximative and limited in scope as global and exact solutions are infeasible for any non-trivial program. In this paper, we present an automated tool to measure the effectmissing static information has on the optimizations applied to a given program. The approach generates an optimistically optimized program version which, compared to the original, defines a performance gap that can be closed by better compiler analyses and selective static program annotations. Our evaluation on six already tuned proxy applications for high-performance codes shows speedups of up to 20.6%. This clearly indicates that static uncertainty limits performance. At the same time, we observed that compilers are often unable to utilize additional static information. Thus, manual annotation of all correct static information is therefore not only error prone but also mostly redundant.

Auto Tuning↗

Griffin Software Development Plan

Griffin is a MOOSE-based reactor physics application for advanced reactor multiphysics modeling and simulation. The application is developed in a consistent multiphysics environment with strong software quality assurance. Griffin inherited most of the capabilities of MAMMOTH/Rattlesnake and is adopting the capabilities from PROTEUS that are needed in the code. The toolset includes a variety of deterministic radiation transport solvers for fixed source, k-eigenvalue, ad-joint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. The code contains the cross-section preparation capabilities applicable to fast and thermal reactors, including TRISO-fueled reactors. Core management capabilities include core performance, fuel depletion and shuffling, equilibrium core calculation, pebble-bed reactor run-in and equilibrium core, molten-salt reactor delayed neutron precursor drift, and control rod and drum movement with cusping correction. This software development plan presents the current and future capabilities and features in Griffin for the design and analysis of non-light-water reactor systems in steady-state and transient conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

RVA: A Plugin for ParaView 3.14

RVA is a plugin developed for the 64-bit Windows version of the ParaView 3.14 visualization package. RVA is designed to provide support in the visualization and analysis of complex reservoirs being managed using multi-fluid EOR techniques. RVA, for Reservoir Visualization and Analysis, was developed at the University of Illinois at Urbana-Champaign, with contributions from the Illinois State Geological Survey, Department of Computer Science and National Center for Supercomputing Applications. RVA was designed to utilize and enhance the state-of-the-art visualization capabilities within ParaView, readily allowing joint visualization of geologic framework and reservoir fluid simulation model results. Particular emphasis was placed on enabling visualization and analysis of simulation results highlighting multiple fluid phases, multiple properties for each fluid phase (including flow lines), multiple geologic models and multiple time steps. Additional advanced functionality was provided through the development of custom code to implement data mining capabilities. The built-in functionality of ParaView provides the capacity to process and visualize data sets ranging from small models on local desktop systems to extremely large models created and stored on remote supercomputers. The RVA plugin that we developed and the associated User Manual provide improved functionality through new software tools, and instruction in the use of ParaView-RVA, targeted to petroleum engineers and geologists in industry and research. The RVA web site (http://rva.cs.illinois.edu) provides an overview of functions, and the development web site (https://github.com/shaffer1/RVA) provides ready access to the source code, compiled binaries, user manual, and a suite of demonstration data sets. Key functionality has been included to support a range of reservoirs visualization and analysis needs, including: sophisticated connectivity analysis, cross sections through simulation results between selected wells, simplified volumetric calculations, global vertical exaggeration adjustments, ingestion of UTChem simulation results, ingestion of Isatis geostatistical framework models, interrogation of joint geologic and reservoir modeling results, joint visualization and analysis of well history files, location-targeted visualization, advanced correlation analysis, visualization of flow paths, and creation of static images and animations highlighting targeted reservoir features.

Source record↗

RVA. 3-D Visualization and Analysis Software to Support Management of Oil and Gas Resources

A free software application, RVA, has been developed as a plugin to the US DOE-funded ParaView visualization package, to provide support in the visualization and analysis of complex reservoirs being managed using multi-fluid EOR techniques. RVA, for Reservoir Visualization and Analysis, was developed as an open-source plugin to the 64 bit Windows version of ParaView 3.14. RVA was developed at the University of Illinois at Urbana-Champaign, with contributions from the Illinois State Geological Survey, Department of Computer Science and National Center for Supercomputing Applications. RVA was designed to utilize and enhance the state-of-the-art visualization capabilities within ParaView, readily allowing joint visualization of geologic framework and reservoir fluid simulation model results. Particular emphasis was placed on enabling visualization and analysis of simulation results highlighting multiple fluid phases, multiple properties for each fluid phase (including flow lines), multiple geologic models and multiple time steps. Additional advanced functionality was provided through the development of custom code to implement data mining capabilities. The built-in functionality of ParaView provides the capacity to process and visualize data sets ranging from small models on local desktop systems to extremely large models created and stored on remote supercomputers. The RVA plugin that we developed and the associated User Manual provide improved functionality through new software tools, and instruction in the use of ParaView-RVA, targeted to petroleum engineers and geologists in industry and research. The RVA web site (http://rva.cs.illinois.edu) provides an overview of functions, and the development web site (https://github.com/shaffer1/RVA) provides ready access to the source code, compiled binaries, user manual, and a suite of demonstration data sets. Key functionality has been included to support a range of reservoirs visualization and analysis needs, including: sophisticated connectivity analysis, cross sections through simulation results between selected wells, simplified volumetric calculations, global vertical exaggeration adjustments, ingestion of UTChem simulation results, ingestion of Isatis geostatistical framework models, interrogation of joint geologic and reservoir modeling results, joint visualization and analysis of well history files, location-targeted visualization, advanced correlation analysis, visualization of flow paths, and creation of static images and animations highlighting targeted reservoir features.

97 MATHEMATICS AND COMPUTING↗

The class 0 protostar BHR71: Herschel observations and dust continuum models

We use Herschel spectrophotometry of BHR71, an embedded Class 0 protostar, to provide new constraints on its physical properties. We detect 645 (non-unique) spectral lines among all spatial pixels. At least 61 different spectral lines originate from the central region. A CO rotational diagram analysis shows four excitation temperature components, 43, 197, 397, and 1057 K. Low-J CO lines trace the outflow while the high-J CO lines are centered on the infrared source. The low-excitation emission lines of H{sub 2}O trace the large-scale outflow, while the high-excitation emission lines trace a small-scale distribution around the equatorial plane. We model the envelope structure using the dust radiative transfer code, hyperion, incorporating rotational collapse, an outer static envelope, outflow cavity, and disk. The evolution of a rotating collapsing envelope can be constrained by the far-infrared/millimeter spectral energy distribution along with the azimuthally averaged radial intensity profile, and the structure of the outflow cavity plays a critical role at shorter wavelengths. Emission at 20–40 μm requires a cavity with a constant-density inner region and a power-law density outer region. The best-fit model has an envelope mass of 19 M{sub ⊙} inside a radius of 0.315 pc and a central luminosity of 18.8 L{sub ⊙}. The time since collapse began is 24,630–44,000 years, most likely around 36,000 years. The corresponding mass infall rate in the envelope (1.2 × 10{sup −5} M{sub ⊙} yr{sup −1}) is comparable to the stellar mass accretion rate, while the mass-loss rate estimated from the CO outflow is 20% of the stellar mass accretion rate. We find no evidence for episodic accretion.

79 ASTRONOMY AND ASTROPHYSICS↗

RVA: A Plugin for ParaView 3.14

RVA is a plugin developed for the 64-bit Windows version of the ParaView 3.14 visualization package. RVA is designed to provide support in the visualization and analysis of complex reservoirs being managed using multi-fluid EOR techniques. RVA, for Reservoir Visualization and Analysis, was developed at the University of Illinois at Urbana-Champaign, with contributions from the Illinois State Geological Survey, Department of Computer Science and National Center for Supercomputing Applications. RVA was designed to utilize and enhance the state-of-the-art visualization capabilities within ParaView, readily allowing joint visualization of geologic framework and reservoir fluid simulation model results. Particular emphasis was placed on enabling visualization and analysis of simulation results highlighting multiple fluid phases, multiple properties for each fluid phase (including flow lines), multiple geologic models and multiple time steps. Additional advanced functionality was provided through the development of custom code to implement data mining capabilities. The built-in functionality of ParaView provides the capacity to process and visualize data sets ranging from small models on local desktop systems to extremely large models created and stored on remote supercomputers. The RVA plugin that we developed and the associated User Manual provide improved functionality through new software tools, and instruction in the use of ParaView-RVA, targeted to petroleum engineers and geologists in industry and research. The RVA web site (http://rva.cs.illinois.edu) provides an overview of functions, and the development web site (https://github.com/shaffer1/RVA) provides ready access to the source code, compiled binaries, user manual, and a suite of demonstration data sets. Key functionality has been included to support a range of reservoirs visualization and analysis needs, including: sophisticated connectivity analysis, cross sections through simulation results between selected wells, simplified volumetric calculations, global vertical exaggeration adjustments, ingestion of UTChem simulation results, ingestion of Isatis geostatistical framework models, interrogation of joint geologic and reservoir modeling results, joint visualization and analysis of well history files, location-targeted visualization, advanced correlation analysis, visualization of flow paths, and creation of static images and animations highlighting targeted reservoir features.

Source record↗

Application of Combined Sustained and Cyclic Loading Test Results to Alloy 617 Elevated Temperature Design Criteria

Alloy 617 is a reference structural material for very high temperature components of advanced-gas cooled reactors with outlet temperatures in the range of 900-950°C . In order for designers to be able to use Alloy 617 for these high temperature components, Alloy 617 has to be approved for use in Section III (the nuclear section) of the American Society of Mechanical Engineers (ASME) Boiler and Pressure Vessel Code. A plan has been developed to submit a draft code for Alloy 617 to ASME Section III by 2015. However, the current rules in Subsection NH for the evaluation of strain limits and creep-fatigue damage using simplified methods based on elastic analysis have been deemed inappropriate for Alloy 617 at temperatures above 1200°F (650°C). The rationale for this exclusion is that at higher temperatures it is not feasible to decouple plasticity and creep deformation, which is the basis for the current simplified rules. This temperature, 1200 °F, is well below the temperature range of interest for this material in High Temperature Gas Cooled Reactor (HTGR) applications. The only current alternative is, thus, a full inelastic analysis which requires sophisticated material models which have been formulated but not yet verified. To address this issue, proposed code rules have been developed which are based on the use of elastic-perfectly plastic (EPP) analysis methods and which are expected to be applicable to very high temperatures.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Cyclic Ramberg–Osgood Parameters of 316SS Weld Metal Under In-Air and PWR-Water Environments

We report at present, the available Ramberg–Osgood (R–O) parameters for different metals (e.g., in ASME code and other literatures) are static (generally based on a tensile curve). These static R–O parameters cannot accurately model the cyclic plasticity behavior. This work presents the cyclic R–O material hardening parameters for 316 stainless steel similar metal welds. The parameters were estimated under various conditions (in-air at room temperature, 300 °C in-air, and in-air at primary water conditions for a pressurized water reactor (PWR)). It is anticipated that the reported results would be useful for computational mechanics-based shakedown analysis and fatigue life estimation of PWR components.

42 ENGINEERING↗

COUPLING OF RELAP5-SCDAP MOD4.0 AND NEUTRONIC CODES

High-fidelity and accurate nuclear system codes play a key role in the design and analysis of complex nuclear power plants, which consist of multiple subsystems, such as the reactor core (and its fuel, burnable poisons, control elements, etc.), the reactor internal structures, the vessel, and the energy conversion subsystem and beyond to grid demand. Most commonly the interplay between these various subsystems is modeled using coupled codes, each of which represents one of the subsystems. And the most common direct coupling is that of thermal-hydraulics and neutronics codes. The subject of this paper is the coupling of codes that model not only thermal-hydraulics and neutronics, but also structural components damage. Furthermore, the neutronic component is not limited to the sole core solver. The coupled code system encompasses thermal-hydraulics, material performance of the fuel, neutronic solver, and neutronic data preparation. Thus, this paper presents a framework for coupling RELAP5/SCDAPSIM/MOD4.0 with a suite of neutron kinetics codes that includes NESTLE, DRAGON and a version of the ENDF library. The version of the RELAP5/SCDAPSIM/MOD4.0 code used in this work is one developed by Innovate System Software (ISS) as part of the international SCDAP Development and Training Program (SDTP) for best-estimate analysis to model reactor transients including severe accident phenomena. This RELAP5/SCDAPSIM/MOD4.0 code version is also capable of predicting nuclear fuel performance. It uses nodal power distributions to calculate mechanical and thermal parameters such as heat-up, oxidation and meltdown of fuel rods and control rods, the ballooning and rupture of fuel rod cladding, the release of fission products from fuel rods, and the disintegration of fuel rods into porous debris and molten material. On the neutronics side, this work uses the NESTLE and DRAGON codes. NESTLE is a multi-dimensional static and kinetic neutronic code developed at North Carolina State University. It solves up to four energy groups neutron diffusion equations utilizing the Nodal Expansion Method (NEM) in Cartesian or hexagonal geometry. The DRAGON code, developed at Ecole Polytechnique de Montreal, performs lattice physics calculations based on the neutron transport equation and is capable of using very fine energy group structures. In this work, we have developed a coupling approach to exchange data among the various modules. In the coupling process, the generated nuclear data (in fine multigroup energy structure) are collapsed down into two- or four-group energy structures for use in NESTLE. The neutron kinetics and thermal-hydraulics modules are coupled at each time step by using the cross-section data. The power distribution results of the neutronic calculations are transmitted to the thermal-hydraulics code. The spatial distribution of coolant density and the fuel-moderator temperature, which result from the thermal-hydraulic calculations, are transmitted back to the neutron kinetics codes and then the loop is closed using new neutronics results. Details of the actual data transfers will be described in the full length paper.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Uncertainty quantification and sensitivity analysis of a nuclear thermal propulsion reactor startup sequence

The research presented in this article describes progress in applying stochastic methods, uncertainty quantification, parametric studies, and variance-based sensitivity analysis (also known as Sobol sensitivity analysis) to a full-core model of a nuclear thermal propulsion (NTP) system simulated via the radiation transport code Griffin to simulate neutronics. Our goal is to develop a reduced-order (surrogate) model that can be rapidly sampled with perturbations to multiple input parameters. In this NTP system, reactivity and power feedback affect the rotation of control drums (CDs), which is itself controlled by a hybrid proportional-integral-derivative (PID) controller actuated by the power demand and reactivity feedback from the numerical model. This model uses reactor kinetic feedback (mean generation time [Λ] and effective delayed neutron fraction [ β eff ] from a transient Griffin simulation executed via Griffin’s improved quasi-static solver to provide the kinetic parameters) as inputs to functions that control the CD rotation angle. By investigating numerous stochastic approaches, we developed a dual-purpose surrogate model of the NTP system, using polynomial regression in the Multiphysics Object-Oriented Simulation Environment (MOOSE) Stochastic Tools Module (STM). The trained model can be rapidly sampled while simultaneously perturbing various input parameters, such as coefficients on the PID control or temperature (directly affecting the neutron cross section). The surrogate model delivers accurate (within 5%) results at speeds orders of magnitude faster (minutes, not days of computational time) than the base model. Once the surrogate model has been trained, distributions of the uncertain parameters can be changed at will to investigate the effects of perturbing multiple inputs as well as the effects of these inputs on the model output. For example, coefficients used in the PID control system may vary due to some type of physical interference, or uncertainty may exist in the temperature of the neutron cross sections in various regions of the reactor. A distribution can be placed on these parameters, and operational boundaries can be determined. The goal of this work is to support development of an advanced control system for operating CDs in a functioning NTP system. This work is a scoping study of the MOOSE STM.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

A Similarity-Based Analysis Tool for Porting OpenMP Applications. In: Facing the Multicore-Challenge III. Lecture Notes in Computer Science, vol 7686.

Exascale computers are expected to exhibit an unprecedented level of complexity, thereby posing significant challenges for porting applications to these new systems. One of the ways to support this transition is to create tools that allow their users to benefit from prior successful porting experiences. The key to such an approach is the manner in which we define source code similarity, and whether similar codes can be ported in the same way to a given system. In this paper, we propose a novel approach based on the notion of similarity that uses static and dynamic code features to check if two serial subroutines can be ported with the same OpenMP strategy. Our approach creates an annotated family distance tree based on the syntactic structure of subroutines, where subroutines that belong to the same syntactic family and share the similar code features have a greater potential to be optimized in the same way. We describe the design and implementation of a tool, based upon a compiler and performance tool, that is used to gather the data to build this porting planning tree. We then validate our approach by analyzing the similarity in subroutines of the serial version of the NAS benchmarks and comparing how they were ported in the OpenMP version of the suite.

Ding, Wei↗

Near Real-time Data Analysis of Core-Collapse Supernova Simulations With Bellerophon

We present an overview of a software system, Bellerophon, built to support a production-level HPC application called CHIMERA, which simulates core-collapse supernova events at the petascale. Developed over the last four years, Bellerophon enables CHIMERA s geographically dispersed team of collaborators to perform data analysis in near real-time. Its n-tier architecture provides an encapsulated, end-to-end software solution that enables the CHIMERA team to quickly and easily access highly customizable animated and static views of results from anywhere in the world via a web-deliverable, cross-platform desktop application. In addition, Bellerophon addresses software engineering tasks for the CHIMERA team by providing an automated mechanism for performing regression testing on a variety of supercomputing platforms. Elements of the team s workflow management needs are met with software tools that dynamically generate code repository statistics, access important online resources, and monitor the current status of several supercomputing resources.

Lingerfelt, Eric J↗

Modeling the contributions to acoustic nonlinearity from complex dislocation networks using 3D dislocation dynamics

Nonlinear ultrasonic parameters are highly sensitive to microstructural features that affect macroscale material behavior, providing a nondestructive means to characterize their evolution. Although dislocations are known to be a strong source of acoustic nonlinearity, establishing quantitative links between the acoustic nonlinearity parameter (β), measured via Second Harmonic Generation, and dislocation morphology—such as dislocation length and density—remains an open challenge. This work advances the numerical modeling of dislocation–β relationships using 3D dislocation dynamics (DD) simulations in two approaches: a “static” method computing strain and stress fields from dislocation configurations in the absence of external loading, and a “quasi-static” method to estimate β from the curvature of dislocation lines under applied load. First, the static method is combined with finite element analysis to investigate a recent assertion that heterogeneous initial strain fields can induce higher harmonic generation in a linear elastic medium; the present results do not corroborate this outcome. Then, the quasi-static method is applied to multiple-dislocation scenarios through parametric studies, revealing behaviors not predicted by analytical models, such as the competing interactions of edge and screw dislocations and the significant influence of applied stress on β. Finally, the simulations are used to model SHG experimental results and validate the hypothesis that β can decrease during plastic deformation, despite increasing dislocation density. As the DD code used here is open-source, it provides a practical platform for future investigation into microstructure–β relationships important to the interpretation of SHG results.

Materials science↗

Data and scripts associated with the manuscript "A Hydrogeophysical Framework to Assess Infiltration During a Simulated Ecosystem-Scale Flooding Experiment”

This package contains the data and scripts used in "A hydrogeophysical framework to assess infiltration during a simulated ecosystem-scale flooding experiment" (Adebayo et al., 2023). The Geophysics folder contains two subfolders containing Ground Penetrating Radar (GPR) and Electrical Resistivity Imaging (ERI) datasets acquired between August 23 and August 26, 2021. The GPR subfolder contains raw data (Tempest_Project.GPZ) and processed results using the Sensors and Software EKKO Project software. The ERI subfolder contains two additional subfolders named Static and Timelapse. The Static subfolder contains background resistivity datasets acquired before infiltration, both raw and inverted results using the ResIpy code. The Timelapse subfolder contains repeated resistivity datasets acquired during infiltration on two profiles, along with raw and inverted results using the ResIpy code. The Flow model folder contains ATS model simulation files. The Geophysics-informed subfolder contains simulation datasets for two profiles generated using information from the geophysics data. The National_database-informed subfolder contains simulation datasets for two profiles generated using information from a national soil database. Each subfolder contains initial simulation results, multiple ensembles, and final simulation files. The Sensitivity_analysis folder contains the observation files exported from the flow model, notebooks used in running the code, and sensitivity analysis results showing the sensitivity of different parameters over a 48-hour period. We used mutual information-based global sensitivity analysis for the sensitivity analysis, following Jiang et al., 2022. This study shows a novel approach that integrates geophysical data with hydrological modeling to understand subsurface flow during coastal flooding.

54 ENVIRONMENTAL SCIENCES↗

Nuclear reactor transient analysis via a quasi-static kinetics Monte Carlo method

The predictor-corrector quasi-static (PCQS) method is applied to the Monte Carlo (MC) calculation for reactor transient analysis. To solve the transient fixed-source problem of the PCQS method, fission source iteration is used and a linear approximation of fission source distributions during a macro-time step is introduced to provide delayed neutron source. The conventional particle-tracking procedure is modified to solve the transient fixed-source problem via MC calculation. The PCQS method with MC calculation is compared with the direct time-dependent method of characteristics (MOC) on a TWIGL two-group problem for verification of the computer code. Then, the results on a continuous-energy problem are presented.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗