Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model code”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Fuel performance analysis of Cr-coated Zircaloy-4 cladding during a prototypical LOCA event using BISON

Here, deformation and failure of chromium (Cr) coated Zircaloy-4 (Zry-4) were studied in loss-of-coolant accident (LOCA) conditions using the BISON fuel performance code. The BISON validation model simulating Halden research reactor experiments was extended to include Cr coatings and higher rod internal pressures to simulate high-burnup fuel. The transient simulations show Cr coatings help relieve stress in the Zry-4 substrate during the transient, delaying the onset of high-temperature creep, which leads to ballooning and bursting of the cladding. A nominal Cr coating thickness of 30 µm delays clad failure by 26 seconds and increases the clad burst temperature by 40 K. A parametric study showed that time to failure and burst temperature both increase with coating thickness, and Cr-coated cladding offers burst resistance under a wide range of rod internal pressures simulating high-burnup fuel. Results indicate that a thin Cr coating provides resistance against ballooning and bursting of cladding during LOCA events.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Developing a robust strength model using physically-informed genetic programming

The strength of materials is influenced by a range of external conditions, such as temperature and deformation rate. Consequently, materials that demonstrate substantial variations in their mechanical behavior due to fluctuations in temperature and strain rate require complex strength models to accurately predict material performance in real-world applications. To predict such complex behavior, a robust and flexible strength model is necessary. In this work, we utilize genetic programming-based symbolic regression (GPSR) to develop data-driven strength models that accurately represent the measured stress–strain responses of tin across a wide range of strain, strain rate and temperature regimes. The GPSR models are constrained by physically-informed conditions, which leads to significant improvement in extrapolation. The best model is integrated into a multi-physics code to perform Taylor impact simulations, validating the model’s accuracy and robustness. In conclusion, the model predictions showed excellent agreement with experimental results, particularly when compared to predictions using traditional strength models.

Genetic programming↗
DOE OSTI2022

Thick target neutron yields from Beryllium, Carbon, Tungsten, and Lead targets irradiated by 26.7 MeV/nucleon 4 H e ions

The angular and energy distributions of secondary neutrons produced from 26.7 MeV/nucleon 4 He ions, stopping in thick Be, C, W and Pb targets are measured by the time of flight method. The GEANT4, PHITS, and FLUKA Monte Carlo simulation codes with different physics models are employed to simulate the neutron yields and the results are compared with the experimental data. For the thick Be and C targets, GEANT4 with INCL++ reproduces the experimental data reasonably well. Whereas PHITS with JQMD-2.0 overestimates the neutron yields at intermediate energies and FLUKA with PEANUT underestimated the neutron yields at intermediate energies. For W and Pb targets, FLUKA with PEANUT physics model agree with experimental data best. PHITS with JQMD-2.0 overestimates the neutron yield at intermediate energies at 0°, and GEANT4 with INCL++ underestimates the neutron yield at intermediate energies at three angles.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

LISE$^{++}_{cute}$, the latest generation of the LISE ++ package, to simulate rare isotope production with fragment-separators

The LISE ++ software for fragment separator simulations has undergone a major update. The package, widely used at rare isotope beam facilities, can be used to predict intensities and purities of rare isotope beams and for planning and running of experiments using in-flight separators. It is especially useful for radioactive beam production as its results can be quickly compared to on-line data. The LISE ++ package has been ported to the Qt-framework in order to support modern compilers and computing methods. The benefits include 64-bit operation and LISE ++ availability on three different platforms: Windows, MacOS and Linux. In addition, the porting provides the ability to take advantage of future computational improvements. The updated package is named LISE$^{++}_{cute}$ to indicate a major step forward from the previous Borland-based versions. In addition to porting to the new platform, new main features and modifications have been added, mostly devoted to improving models and implementing other codes involved in rare isotope production at FRIB. Finally, a summary of modifications completed to improve the functionality of the code are discussed in this work, as well as future plans.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Delivery drone route planning over a battery swapping network

Many enterprises invest on drone delivery research and development to drop off packages at consumers’ doorsteps in a matter of minutes. We study delivery drone route planning over a battery swapping network allowing farther reach by penetrating current battery capacity constraints. A mixed-integer nonlinear programming model is created to plan efficient drone routing over the swapping machines by minimizing the delivery lead time. We develop an exact solution method, evaluate its performance, and compare it with a straightforward nonlinear solver application. A case study highlights the applicability of the model. Data and source code to the solver are publicly shared.

drone battery swapping↗
DOE OSTI2021

Synthesis of magnesiowüstite nanocrystallites embedded in an amorphous silicate matrix via low energy multiple ion implantations

The synthesis process is presented for experimentally simulating modifications in cosmic dust grains using sequential ion implantations or irradiations followed by thermal annealing. Cosmic silicate dust analogues were prepared via implantation of 20–80 keV Fe - , Mg - , and O - ions into commercially available p-type silicon (100)wafers. The as-implanted analogues are amorphous with a Mg/(Fe+Mg) ratio of 0.5 tailored to match theoretical abundances in circumstellar dusts. Before the ion implantations were performed, Monte-Carlo-based ion-solid interaction codes were used to model the dynamic redistribution of the implanted atoms in the silicon substrate. 600 keV helium ion irradiation was performed on one of the samples before thermal annealing. Two samples were thermally annealed at a temperature appropriate for an M-class stellar wind, 1000 K, for 8.3 h in a vacuum chamber with a pressure of 1 x 10 -7 torr. The elemental depth profiles were extracted utilizing Rutherford Backscattering Spectrometry (RBS) in the samples before and after thermal annealing. X-ray diffraction (XRD)analysis was employed for the identification of various phases in crystalline minerals in the annealed analogues. Transmission electron microscopy (TEM) analysis was utilized to identify specific crystal structures. RBS analysis shows redistribution of the implanted Fe, Mg, and O after thermal annealing due to incorporation into the crystal structures for each sample type. XRD patterns along with TEM analysis showed nanocrystalline Mg and Fe oxides with possible incorporation of additional silicate minerals.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Evaluating wind speed and power forecasts for wind energy applications using an open-source and systematic validation framework

Building on the verification and validation work developed under the Second Wind Forecast Improvement Project, this work exhibits the value of a consistent procedure to evaluate wind power forecasts. We established an open-source Python code base tailored for wind speed and wind power forecast validation, WE-Validate. The code base can evaluate model forecasts with observations in a coherent manner. To demonstrate the systematic validation framework of WE-Validate, we designed and hosted a forecast evaluation benchmark exercise. We invited forecast providers in industry and academia to participate and submit forecasts for two case studies. We then evaluated the submissions with WE-Validate. Our findings suggest that ensemble means have reasonable skills in time series forecasting, whereas they are often inferior to single ensemble members in wind ramp forecasting. Adopting a voting scheme in ramp forecasting that allows ensemble members to detect ramps independently leads to satisfactory skill scores. Throughout this document, we also emphasize the importance of using statistically robust and resistant metrics as well as equitable skill scores in forecast evaluation.

17 WIND ENERGY↗
DOE OSTI2025

Evidence for an Early Formation of Serenitatis Basin at 4.25 Ga Shifts Lunar Chronology

Troctolite sample 76535, collected in Serenitatis basin during Apollo 17, formed at least 50 km deep, experienced maximum shock pressures of 6 GPa, and has a 40 Ar/ 39 Ar excavation age of 4.25 Ga. Previous work attributed 76535 to the South Pole-Aitken (SPA) basin, presumably dating the SPA-impact and constraining lunar bombardment history. Here we use the iSALE-2D shock-physics code and gravity inversion modeling to determine if instead the Serenitatis impact event excavated 76535. We find nearly 140,000 km 3 of material (∼2% of near-surface ejecta) matching the depth and pressure constraints of 76535 is displaced to the surface during crater collapse of a Serenitatis-like impact event. We conclude that the Serenitatis impact event possibly excavated 76535, redefining its age to 4.25 Ga, 300 My older than the consensus age based on Apollo 17 samples. This finding would provide an important anchor point where lunar chronology where bombardment flux is especially uncertain.

Bjonnes, Evan [Lawrence Livermore National Laborat↗
DOE OSTI2022

A Dataset of 3D Structural and Simulated Transport Properties of Complex Porous Media

Physical processes that occur within porous materials have wide-ranging applications including - but not limited to - carbon sequestration, battery technology, membranes, oil and gas, geothermal energy, nuclear waste disposal, water resource management. The equations that describe these physical processes have been studied extensively; however, approximating them numerically requires immense computational resources due to the complex behavior that arises from the geometrically-intricate solid boundary conditions in porous materials. Here, we introduce a new dataset of unprecedented scale and breadth, DRP-372: a catalog of 3D geometries, simulation results, and structural properties of samples hosted on the Digital Rocks Portal. The dataset includes 1736 flow and electrical simulation results on 217 samples, which required more than 500 core years of computation. This data can be used for many purposes, such as constructing empirical models, validating new simulation codes, and developing machine learning algorithms that closely match the extensive purely-physical simulation. This article offers a detailed description of the contents of the dataset including the data collection, simulation schemes, and data validation.

3D images↗
DOE OSTI2024

Wind Turbine Design Optimization for Hydrogen Production

To help meet the need for inexpensive green fuels, we are working on wind turbine design optimization specifically for hydrogen production. We have thus far achieved a 1.53% decrease in LCOH as compared to a turbine optimized for LCOE using the same code, design variables, and models. We accomplished this by optimizing some components of the wind turbine tower, rotor, and drivetrain design with hydrogen production and costs in the design loop.

hydrogen↗
DOE OSTI2025

Impact of T - and ρ -dependent decay rates and new (n, γ ) cross-sections on the s process in low-mass asymptotic giant branch stars

Aims. We study the impact of nuclear input related to weak-decay rates and neutron-capture reactions on predictions for the slow neutron-capture process (s process) in asymptotic giant branch (AGB) stars. We provide the first database of surface abundances and stellar yields of the isotopes heavier than iron from the Monash models. Methods. We ran nucleosynthesis calculations with the Monash post-processing code for seven stellar structure evolution models of low-mass AGB stars with three different sets of nuclear inputs. The reference set has constant decay rates and represents the set used in the previous Monash publications. The second set contains the temperature and density dependence of β decays and electron captures based on the default rates of nuclear NETwork GENerator (NETGEN). In the third set, we further update 92 neutron-capture rates based on re-evaluated experimental cross sections from the ASTrophysical Rate and rAw data Library. We compare and discuss the predictions of the sets relative to each other in terms of isotopic surface abundances and total stellar yields. We also compare the results to isotopic ratios measured in presolar stardust silicon carbide (SiC) grains from AGB stars. Results. The new sets of models result in a ∼66% solar s-process contribution to the p-nucleus 152 Gd, confirming that this isotope is predominantly made by the s process. The nuclear input updates result in predictions for the 80 Kr/ 82 Kr ratio in the He intershell and surface 64 Ni/ 58 Ni, 94 Mo/ 96 Mo, and 137 Ba/ 136 Ba ratios that are more consistent with the corresponding ratios measured in stardust; however, the new predicted 138 Ba/ 136 Ba ratios are higher than the typical values of the SiC grains. The W isotopic anomalies are in agreement with data from the analyses of other meteoritic inclusions. We confirm that the production of 176 Lu and 205 Pb is affected by too large uncertainties in their decay rates from NETGEN.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Laser-driven flyer plate impact: Computational studies guided by experiments

We present a computational approach using a multimaterial, arbitrary Lagrangian–Eulerian code termed ALE3D to model the nanosecond/micrometer dynamics of the launch of 0.5–4.5 km/s laser-driven metal flyer plates and the impact with stationary targets of Pyrex and fused silica glasses, and Lexan and Plexiglas polymers, producing pressures in the target in the range of 5–20 GPa. The simulations are compared to experimental results where the flyer velocity profile and the velocity profile imparted to the target material were measured with high-speed velocimetry. The experimental flyer launch by a high-intensity pulsed laser is modeled by depositing heat into a thin vaporizable layer under the flyer plate. This model produces a flyer plate that has not been exposed to the laser pulse, allowing us to compare the properties of the real flyer to a simulated ideal flyer. The simulations of target impact are generally in good agreement with the experiment except at the highest impact velocities where the shock release process in the experiment is slower than that in the simulation. The cause of this disagreement is attributed to an inadequate description of the shock viscosity during the nanosecond unloading process.

Stekovic, Svjetlana (ORCID:0000000163878070)↗
DOE OSTI2021

Cross-beam energy transfer saturation by ion trapping-induced detuning

The performance of direct-drive inertial confinement fusion implosions relies critically on the coupling of laser energy to the target plasma. Cross beam energy transfer (CBET), the resonant exchange of energy between intersecting laser beams mediated by ponderomotively driven ion-acoustic waves (IAW), inhibits this coupling by scattering light into unwanted directions. Furthermore, the variety of beam intersection angles and varying plasma conditions in an implosion results in IAWs with a range of phase velocities. Here we show that CBET saturates through a resonance detuning that depends on the IAW phase velocity and that results from trapping-induced modifications to the ion distribution functions. For smaller phase velocities, the modifications to the distribution functions can rapidly thermalize in the presence of mid-Z ions, leading to a blueshift in the resonant frequency. For larger phase velocities, the modifications can persist, leading to a redshift in the resonant frequency. Ultimately, these results may reveal pathways towards CBET mitigation and inform reduced models for radiation hydrodynamics codes to improve their predictive capability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Deceleration-stage Rayleigh–Taylor growth in a background magnetic field studied in cylindrical and Cartesian geometries

Experiments have identified the Rayleigh–Taylor (RT) instability as one of the greatest obstacles to achieving inertial confinement fusion. Consequently, mitigation strategies to reduce RT growth and fuel–ablator mixing in the hotspot during the deceleration phase of the implosion are of great interest. In this work, the effect of seed magnetic fields on deceleration-phase RT growth are studied in planar and cylindrical geometries under conditions relevant to the National Ignition Facility (NIF) and Omega experiments. The magnetohydrodynamic (MHD) and resistive-MHD capabilities of the FLASH code are used to model imploding cylinders and planar blast-wave-driven targets. Realistic target and laser parameters are presented that suggest the occurrence of morphological differences in late-time RT evolution in the cylindrical NIF case and a measurable difference in spike height of single-mode growth in the planar NIF case. The results of this study indicate the need for target designs to utilize an RT-unstable foam–foam interface in order to achieve sufficient magnetic field amplification to alter RT evolution. Benchmarked FLASH simulations are used to study these magnetic field effects in both resistive and ideal MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Automatic Differentiation in MetaPhysicL and Its Applications in MOOSE

Efficient solution via Newton’s method of nonlinear systems of equations requires an accurate representation of the Jacobian, corresponding to the derivatives of the component residual equations with respect to the degrees of freedom. In practice these systems of equations often arise from spatial discretization of partial differential equations used to model physical phenomena. These equations may involve domain motion or material equations that are complex functions of the systems’ degrees of freedom. Computing the Jacobian by hand in these situations is arduous and prone to error. Finite difference approximations of the Jacobian or its action are prone to truncation error, especially in multiphysics settings. Symbolic differentiation packages may be used, but often result in an excessive number of terms in realistic model scenarios. An alternative to symbolic and numerical differentiation is automatic differentiation (AD), which propagates derivatives with every elementary operation of a computer program, corresponding to continual application of the chain rule. Automatic differentiation offers the guarantee of an exact Jacobian at a relatively small overhead cost. In this work, we outline the adoption of AD in the Multiphysics Object Oriented Simulation Environment (MOOSE) via the MetaPhysicL package. We describe the application of MOOSE’s AD capability to several sets of physics that were previously infeasible to model via hand-coded or Jacobian-free simulation techniques, including arbitrary Lagrangian-Eulerian and level-set simulations of laser melt pools, phase-field simulations with free energies provided through neural networks, and metallic nuclear fuel simulations that require inner Newton loop calculation of nonlinear material properties.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Impact of High-Reactivity Advanced Test Reactor Experiments on Photon Heating in Nearby Experiment Locations

The Advanced Test Reactor’s (ATR’s) distinctive ability to provide a wide range of irradiation conditions is attractive for programs pursuing fuel qualification experiments. These potentially high-fuel-load experiments are a relatively new development and produce unexplored effects on nearby experiments. Here, this paper explores how photon heating of such an experiment may affect other nearby experiment programs, ultimately serving to better inform decisions regarding experiment design and risks to programmatic goals. The MC21 (Monte Carlo for the 21st Century) code is used to model and study how gamma heat generation rates and axial effects impact different ATR positions. The results reveal that the proximity of a given experiment’s position to the high-fuel-load one can significantly alter that experiment’s expected axial profile.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Challenges and Opportunities Specific to Microreactors: Compact and Light-Weight Shielding During Operation and Transportation

The problem of designing a shield for a mobile microreactor has been addressed. The principal challenges encountered in this endeavor are the dual goals of increasing radiation protection and decreasing weight. It is recognized that these two goals are contradictory. Increasing shielding naturally increases weight. Conversely, decreasing weight naturally calls for decreasing material inventory, including possibly decreasing shield mass. The operational and regulatory drivers for shielding and transportation constraints are identified and summarized. The main constraints pertain to volume, mass, radiation dose, material feasibility and availability, and decommissioning expectations. The latter are limits on residual contamination and dose rates after the reactor is removed. Progress toward solving the challenge of shielding a transportable microreactor is made by considering a tungsten tetraboride (WB4) shell surrounding the reactor vessel. The shield and a generic microreactor were modeled using the MCNP code. It was determined that the presence of a 16-cm-thick WB4 shield reduces the dose equivalent rate during operations to 3.6 mrem/h just outside of a concrete vault and dome and soil berm surrounding the reactor, which is very near the occupational dose rate limits. Without the WB4 layer, the dose rate at the same location outside the berm would be 120 mrem/h. Most importantly, the presence of the WB4 shield substantially protects the surrounding concrete and soil from activation, thus decreasing the duration of administrative control for the site following removal of the reactor from 16 to 3.4 years. At that point, the site may be decommissioned and returned to unrestricted public access. A limitation on the incorporation of a WB4 shield as an intrinsic component of a microreactor is that it results in a design that is transportable only by using the largest road conveyances and that the load would require special permits for transit through the U.S. highway system. New directions for further improvements and possibly more effective solutions to the shielding problem of a mobile microreactor are discussed in the conclusion section.

Microreactor↗