Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer experiment”

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,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Standard Model prediction of the Bc lifetime

Applying an operator product expansion approach we update the Standard Model prediction of the B c lifetime from over 20 years ago. The non-perturbative velocity expansion is carried out up to third order in the relative velocity of the heavy quarks. The scheme dependence is studied using three different mass schemes for the b ¯ and c quarks, resulting in three different values consistent with each other and with experiment. Special focus has been laid on renormalon cancellation in the computation. Uncertainties resulting from scale dependence, neglecting the strange quark mass, non-perturbative matrix elements and parametric uncertainties are discussed in detail. The resulting uncertainties are still rather large compared to the experimental ones, and therefore do not allow for clear-cut conclusions concerning New Physics effects in the B c decay.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Margin to onset of nucleate boiling and flow instability studies for preliminary MITR design-demonstration element thermal-hydraulics

This study covers the onset of flow instability (OFI) preliminary results obtained from leveraging correlations, in addition to the preliminary thermal hydraulics results such as pressure, flow velocity, temperature, and oxide layer over the design demonstration experiment (DDE) for the Massachusetts Institute of Technology Reactor (MITR). Current computational fluid dynamics (CFD) models in fluid structure interaction (FSI) have added the capability of assessing margins to onset of nucleate boiling (ONB). This study initiates the capability to model the margin to OFI and ONB presented for the MITR. Such study is supportive of the United States High Performance Research Reactor (USHPRR) program. Previous studies provided preliminary thermal-hydraulic and mechanical analyses of the hydrodynamic effects in the MITR DDE under conservative approximations for plate power distribution. In conclusion, this study focuses on providing insights into the OFI future research direction optimizing the transport of thermal energy, mass-flow rates, flow-channel geometries, and boundary conditions.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Accurate simulation of surfaces and interfaces of ten FCC metals and steel using Lennard–Jones potentials

Abstract The earlier integration of validated Lennard–Jones (LJ) potentials for 8 fcc metals into materials and biomolecular force fields has advanced multiple research fields, for example, metal–electrolyte interfaces, recognition of biomolecules, colloidal assembly of metal nanostructures, alloys, and catalysis. Here we introduce 12-6 and 9-6 LJ parameters for classical all-atom simulations of 10 further fcc metals (Ac, Ca (α), Ce (γ), Es (β), Fe (γ), Ir, Rh, Sr (α), Th (α), Yb (β)) and stainless steel. The parameters reproduce lattice constants, surface energies, water interfacial energies, and interactions with (bio)organic molecules in 0.1 to 5% agreement with experiment, as well as qualitative mechanical properties under standard conditions. Deviations are reduced up to a factor of one hundred in comparison to earlier Lennard–Jones parameters, embedded atom models, and density functional theory. We also explain a quantitative correlation between atomization energies from experiments and surface energies that supports parameter development. The models are computationally very efficient and applicable to an exponential space of alloys. Compatibility with a wide range of force fields such as the Interface force field (IFF), AMBER, CHARMM, COMPASS, CVFF, DREIDING, OPLS-AA, and PCFF enables reliable simulations of nanostructures up to millions of atoms and microsecond time scales. User-friendly model building and input generation are available in the CHARMM-GUI Nanomaterial Modeler. As a limitation, deviations in mechanical properties vary and are comparable to DFT methods. We discuss the incorporation of reactivity and features of the electronic structure to expand the range of applications and further increase the accuracy.

Kanhaiya, Krishan (ORCID:0000000236222655)↗
DOE OSTI2024

ATLAS Data Analysis using a Parallel Workflow on Distributed Cloud-based Services with GPUs

A new type of parallel workflow is developed for the ATLAS experiment at the Large Hadron Collider, that makes use of distributed computing combined with a cloud-based infrastructure. This has been developed for a specific type of analysis using ATLAS data, one popularly referred to as Simulation-Based Inference (SBI). The JAX library is used for the parts of the workflow to compute gradients as well as accelerate program execution using just-in-time compilation, which becomes essential in a full SBI analysis and can also offer significant speed-ups in more traditional types of analysis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Shock tube simulations for the three-layer Richtmyer–Meshkov instability with single-mode and multimode perturbations

While the canonical two-component, single-mode Richtmyer–Meshkov instability (RMI) has been extensively studied, relatively less work has focused on the effects of an additional intermediate-density middle layer. This work investigates such three-material RMI configurations at two Atwood number scenarios using the ARES hydrodynamics code. After validation against previous experimental and computational studies, setups corresponding to recent three-layer shock tube experiments are simulated. Cases with both single-mode and multimode perturbations are studied to quantify mixing across the interface between the materials with highest and intermediate density. In particular, this work is able to comprehensibly examine differences between two- and three-dimensional setups for the single-mode and multimode problems. Observations from previous two-layer investigations still apply in the three-layer setup, but over the time horizons considered, there appears to be insufficient nonlinear mode coupling to create significant differences between two- and three-dimensional simulations following the first passage of a shock. Finally, additional reshock simulations have additional nonlinear growth that does result in expected differences between two- and three-dimensional cases in this three-layer setup, but significant differences do not manifest during the time horizon studied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Development of a data-driven neural network model for electron thermal transport in NSTX

A data-driven electron thermal transport neural network (ETT-NN) model, trained on TRANSP interpretative analysis results of National Spherical Torus Experiment (NSTX), was developed to enable faster and more accurate ETT computation for spherical tokamaks (STs). The model incorporates both convolutional NNs and recurrent NNs, allowing it to simultaneously account for the spatial and temporal non-localities and multi-scale features of turbulent transport, which have been considered only in a limited manner in conventional models. The model was validated through interpretative analysis and predictive simulations using Tokamak Reactor Integrated Automated Suite for Simulation and Computation, demonstrating relatively high accuracy. Additionally, parameter scans were performed on test discharges known to exhibit specific turbulent modes, such as microtearing mode, trapped electron mode, kinetic ballooning mode, and electron temperature gradient mode. The scanning results revealed that the ETT-NN model exhibits the same trends as those observed in conventional gyrokinetic simulations or theories, while also capturing the global nature of turbulent transport, indicating that the data-driven model accurately reflects the underlying physical characteristics. Furthermore, due to the dimensionless nature of the model, we can feasibly expand its applicability by incorporating data from other devices and uncovering the characteristics of ETT in STs in the future.

NSTX↗
DOE OSTI2024

0D, 1D, 2D, and 3D simulations of an idealized coaxial impedance-matched Marx generator

We have conducted 0D, 1D, 2D, and 3D simulations of an idealized coaxial impedance-matched Marx generator (IMG) []. The 0D calculations were conducted with a four-element circuit model; the 1D, 2D, and 3D calculations were conducted with highly resolved, fully electromagnetic representations. The IMG consists of 30 stages distributed axially and connected electrically in series. Each stage is powered by two bricks separated by 180° and connected electrically in parallel. Each brick comprises two opposite-polarity capacitors in series with a single switch. The bricks drive an internal impedance-matched coaxial transmission line terminated by a resistive load. The simulations neglect effects due to the switch-triggering circuit, the capacitor-charging circuit, external conducting boundaries, and reactive components of the load. We find dimensionality does not significantly affect the electrical power delivered by the IMG to its load: peak load powers estimated by the 0D, 1D, 2D, and 3D simulations agree to within 1%. The 3D calculations demonstrate that electromagnetic power radiated by the bricks, and axial gaps between stages, reduces the peak load power by less than ∼ 1 % . Each simulation assumes the load impedance is 34% above that at which the load power is maximized. Operating an IMG with such an overmatched load offers several advantages while decreasing the peak load power by only 2%. The 0D, 1D, 2D, and 3D models outlined herein could be adapted to assess computationally competing IMG designs, and conduct a variety of numerical IMG experiments, an IMG is constructed. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Pion-Induced Radiative Corrections to Neutron β Decay

In this work, we compute the electromagnetic corrections to neutron β decay using a low-energy hadronic effective field theory. We identify new radiative corrections arising from virtual pions that were missed in previous studies. The largest correction is a percent-level shift in the axial charge of the nucleon proportional to the electromagnetic part of the pion-mass splitting. Smaller corrections, comparable to anticipated experimental precision, impact the β-ν angular correlations and the β asymmetry. We comment on implications of our results for the comparison of the experimentally measured nucleon axial charge with first-principles computations using lattice QCD and on the potential of β decay experiments to constrain beyond-the-standard-model interactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Absence of High-Pressure Ground-State Reentrant Ferroelectricity in PbTiO 3

We study ferroelectricity in the classic perovskite ferroelectric PbTiO 3 to high pressures with density functional theory (DFT) and experimental diamond-anvil techniques. We use second harmonic generation spectroscopy to detect lack of inversion symmetry. Consistent with early understanding and experiments, we find that ferroelectricity disappears at moderate pressures. However, DFT computations show that the disappearance arises from the overtaking of zone boundary instabilities, and not from the squeezing out of the off-centering ferroelectric displacements with pressure, as previously thought. Moreover, at high pressures the distorted perovskite phases are metastable with respect to a new dense centrosymmetric postperovskite phase with P⁢2 1 /m symmetry and 8-coordinated Ti, which becomes stable at about 70 GPa.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Rounding Error Analysis of Mixed Precision Block Householder QR Algorithms

Although mixed precision arithmetic has recently garnered interest for training dense neural networks, many other applications could benefit from the speedups and lower storage cost if applied appropriately. The growing interest in employing mixed precision computations motivates the need for rounding error analysis that properly handles behavior from mixed precision arithmetic. We develop mixed precision variants of existing Householder QR algorithms and show error analyses supported by numerical experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Demonstration and Validation of Additively Manufactured Intensified Device for Enhanced Carbon Capture

In FY 2018 and 2019, ORNL personnel designed, prototyped, characterized, and tested the enhancement of amine-based carbon capture using additively manufactured intensified devices, i.e., devices combining multiple thermodynamic and kinetic processes into one unit operation. State of the art scrubbing technology often employs inter-stage cooling which requires multiple stages of single-purpose unit operations, leading to larger equipment size, higher capital cost, and sometimes less than optimal operating conditions for each piece of equipment. It was proposed that integrating multiple thermodynamic and kinetic operations into one unit could mitigate or negate these drawbacks. Heat and mass transfer studies demonstrated that an additively manufactured intensified device could significantly improve capture efficiency at bench scale. This project is a continuation of FEAA130 and aims at (1) demonstration and validation of enhanced CO 2 capture with 3D printed intensified devices (i.e., mass exchanger packing with internal cooling channels), for low-aqueous-solvent (LAS) based capture at the ORNL-constructed laboratory-scale (4.5-feet tall) column, namely Column A, with built-in cooling and (2) computational scoping on the conceptual design of a flexible and modular larger column, namely Column B, at ORNL for further demonstration of enhanced capture for aqueous solvent. To execute this multi-faceted project scope, an integrated project team from ORNL’s Manufacturing Science Division is tasked with applying capabilities in computational fluid dynamics (CFD), additive manufacturing, and absorber-scale demonstration/validation experiments.

42 ENGINEERING↗
DOE OSTI2022

LLNL Explosives Anisotropy Research

Lawrence Livermore National Laboratory scientists and engineers led a multi-institutional team in executing a series of high explosives tests that successfully demonstrated fundamental principles of anisotropy, a possible enabler for improved weapon and munition safety. Working under snowy and frigid conditions on Idaho’s Snake River Plain, a 13-member team from LLNL carried out 52 explosives shots over four days in mid-November at the Idaho National Laboratory’s (INL) National Security Test Range (NSTR) to complete the study. The broader anisotropy (ANISO) team included high explosives handlers and volunteers from INL, Los Alamos National Laboratory, Marine Raiders from the Marine Special Operations Command and members of the U.S. Special Operations Command. The purpose of the study was to explore theoretical methods of creating anisotropic explosives — explosives that perform differently depending on the direction the detonation wave moves through the explosive — by engineering certain physical features in the charges and obtaining basic data from testing. The work is part of an overall effort by the Lab to develop anisotropic explosives that could be used in munitions to reduce the severity and lethality of an unintended detonation without sacrificing performance. The data gathered from the study will be used to design and construct follow-on experiments at LLNL’s High Explosives Applications Facility (HEAF) and validate computer models for future anisotropic assemblies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Advancement of hybrid fluid-kinetic modeling for HEDP and ICF science

We report on the development progress of a hybrid fluid-kinetic code for simulating fluids and plasmas in a wide range of environments, such as laser–matter interactions, inertial confinement fusion, magnetic confinement fusion, and pulsed power. The suite of numerical tools under development utilizes heterogeneous computer architectures and leverages the benefits of particle–based simulation techniques. By working to combine the kinetic particle-in-cell (PIC) model with a particle-based fluid simulation technique, such as smoothed particle hydrodynamics, we are developing a flexible framework capable of accurately modeling complex flows within and between kinetic and fluid regimes. The TriForce code is under development as a C++ framework for parallel, 3D, particle-based, hybrid fluid-kinetic plasma simulations. The fluid half of TriForce will be based upon the meshless smoothed-particle-hydrodynamics (SPH) approach, well-suited for shear, mixing, and turbulence, whereas the kinetic half resembles a traditional particle-in-cell (PIC) code; other particle-based approaches to fluid modeling that do use a mesh are also possible to use and are under investigation. Maxwell’s electromagnetic field equations are solved either via explicit or implicit algorithms or approximated via resistive magnetohydrodynamics (MHD) using an Ohm’s law and resulting induction equation (extended MHD is under development). A primary goal of enabling direct comparisons, from the same code, between results from the variants of MHD and implicit electromagnetic solutions is to improve our fundamental understanding of systems with magnetic fields. The code is under development to recover results from both radiation-MHD and fully kinetic codes in those limits, and is continuing to be developed from other follow-on grants to operate in between where both descriptions may co-exist and interact. For certain applications, it is desired for a simulation to contain fluid ions and electrons as well as kinetic ions and electrons. Typically, it is too computationally intensive to model a full-scale ICF or HEDP experiment fully kinetically since many cycles are expended with very small time steps on modeling the fluid part of a material that is well treated by the fluid approximation. In this case, many traditional PIC particles can be replaced with a single fluid particle representing the thermal part of the distribution function, and there are fewer needed kinetic particles, which describe the non-thermal part and can be sub-cycled relative to the fluid particle advance. Furthermore, a pure fluid code may, depending on the problem, simply lack many physically important details that are beyond the scope of its reduced approximations and assumptions. In this report, we summarize the objectives achieved in the development of the collisional and kinetic half of the code, and the physics problems to which the code has been applied in the areas of advanced and innovative fusion concepts, pulsed power, and magneto-inertial fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Operational Evolution of FTS3: A DevOps Driven Approach to Elastic Operations

The File Transfer Service (FTS3) is a distributed data movement service developed at CERN and widely used to transfer data across the Worldwide LHC Computing Grid (WLCG). At Fermilab, FTS3 supports data transfers for multiple experiments, including Intensity Frontier experiments such as DUNE, enabling reliable data movement between WebDAV endpoints in Europe and the Americas.​ At CHEP 2021, we reported on the initial containerized deployment of FTS3 on OKD, the community Kubernetes distribution of Red Hat OpenShift. In this work, we present the subsequent evolution of this deployment, focusing on new operational capabilities introduced to improve scalability, robustness, and long-term maintainability.​ We describe the adoption of more secure and reproducible container build workflows, the integration of DevOps-driven operational practices, and enhancements in monitoring and automation. A key new result is the introduction of horizontal scaling and elastic resource management, allowing FTS3 components to dynamically adapt to workload variations while maintaining service reliability. We also discuss improvements in fault tolerance and operational procedures derived from production experience.​ Finally, we summarize lessons learned from operating FTS3 as a Kubernetes-native service and outline how these developments have improved the resilience and efficiency of data movement operations at Fermilab.

Munoz Flores, Victor Leopoldo [Fermilab]↗
DOE OSTI2020

Simulating Real Atmospheric Boundary Layers at Gray-Zone Resolutions: How Do Currently Available Turbulence Parameterizations Perform?

Recent computational and modeling advances have led a diverse modeling community to experiment with atmospheric boundary layer (ABL) simulations at subkilometer horizontal scales. Accurately parameterizing turbulence at these scales is a complex problem. The modeling solutions proposed to date are still in the development phase and remain largely unvalidated. This work assesses the performance of methods currently available in the Weather Research and Forecasting (WRF) model to represent ABL turbulence at a gray-zone grid spacing of 333 m. We consider three one-dimensional boundary layer parameterizations (MYNN, YSU and Shin-Hong) and coarse large-eddy simulations (LES). The reference dataset consists of five real-case simulations performed with WRF-LES nested down to 25 m. Results reveal that users should refrain from coarse LES and favor the scale-aware, Shin-Hong parameterization over traditional one-dimensional schemes. Overall, the spread in model performance is large for the cellular convection regime corresponding to the majority of our cases, with coarse LES overestimating turbulent energy across scales and YSU underestimating it and failing to reproduce its horizontal structure. Despite yielding the best results, the Shin-Hong scheme overestimates the effect of grid dependence on turbulent transport, highlighting the outstanding need for improved solutions to seamlessly parameterize turbulence across scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Performance of Closed Loop Venturi Aspirated Aeration System: Experimental Study and Numerical Analysis with Discrete Bubble Model

In wastewater treatment plants, aeration plays a significant role as it encourages aerobic respiration of microbes, which are necessary to break down carbonaceous matter in the waste stream. This process can account for the majority of energy use in wastewater treatment plants. The aeration process is also necessary in odor control in lagoons and in the aquaculture industry. Generally, the aeration process is accomplished with compressors or blowers which can be of low efficiency due to ideal gas laws. This study introduces the idea of increasing aeration efficiency by looping water from a reservoir through a piping network which includes a venturi aspirator at its inlet. For this purpose, an experimental study has been conducted in a laboratory setup with a pump which pulls water from a tank, passes it through a Venturi aspirator and a helical piping network intended to increase bubble residence time, before depositing it back into the bulk fluid tank. This same process is modeled computationally, using a discrete bubble method (DBM), with good agreement with experiments. The overall purpose here is to determine the optimal configurations for standard aeration efficiency (SAE) and the standard oxygen transfer rate (SOTR). A parametric study has been implemented using the DBM based on different hydraulic and flow parameters. The model is also used to predict the SAE of a hypothetical aeration system. Results indicate that it is possible to achieve SAE values in the range of surface aerators or submersed jet aerators using the proposed aeration system with less complex components, thereby decreasing overall costs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗