Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “materials simulations”

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 OSTI2025

Magnetically tuned metal-insulator transition in LaAlO 3 /SrTiO 3 nanowire arrays

A wide family of two-dimensional (2D) systems, including stripe-phase superconductors, sliding Luttinger liquids, and anisotropic 2D materials, can be modeled by an array of coupled one-dimensional (1D) electron channels or nanowire arrays. Here we report experiments in arrays of conducting nanowires with gate and field tunable interwire coupling, that are programmed at the LaAlO 3 /SrTiO 3 interface. We find a magnetically tuned metal-to-insulator transition in which the transverse resistance of the nanowire array increases by up to four orders of magnitude, which can be further tuned with a gate voltage. To explain this behavior, we develop a minimal model of a coupled two-wire system where a Wenzel-Kramers-Brillouin-based approach is used to estimate the transverse tunneling conductance. We demonstrate the existence of distinct conductance features and highlight the crucial role played by the field dependence of the interwire potential barrier on transport properties. Since our model makes minimal assumptions, we expect our predictions to hold for a wide class of coupled 1D systems. The nanowire arrays can serve as model systems to understand the origin of exotic behavior in correlated materials via analog quantum simulation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

singularity-opac

singularity-opac is a library for providing a unified interface for opacities, emissivities, and scattering cross-sections for materials and use in simulation codes. It is designed to be performance portable and run on CPU and GPU.

Dolence, Joshua↗
DOE OSTI2024

Annual Report on NCSP Technical Support task in BNL during FY24

This work reports on the investigations done during FY 24 in collaboration with two 2024 Summer undergraduate interns through the DOE SULI and BNL SURP programs. These efforts were partially funded by the Nuclear Criticality Safety Program through the Technical Support Succession Plan task. The first project, developed by Ian Snider, focused on testing the impact of uncertainties in thermal cross sections for multiple materials in nuclear reactor simulations. The second project aimed implement concurrent into BNL’s machine-learning code to correct spin mis-assignments in neutron resonances, the Bayesian Resonance Reclassifier (BRR).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Heat Transfer in Void Generating Foam Decomposition: Further Development

Continued development of the additive conductivity material model, used to simulate changes in heat transfer that occurs in void generating foam decomposition, has resulted in an improved model and new features. The previous version of the model was calibrated against the Aria Bulk Fluid Element (BFE) solution and proposed a third-order polynomial correction term best captured the increased heat transfer due to voids in the foam. An investigation of the Fuego Conjugate Heat Transfer (CHT) and Aria BFE solutions at several geometries revealed the CHT solution and BFE solution had differing behavior across length scales, especially at smaller scales. Five calibration studies, using the Fuego CHT as the calibration data, were carried out with polynomial functions of 4-th, 3-rd, 2-nd, 1-st and 0-th orders to determine the best correction function that generalized well across length scales. Each polynomial function was calibrated/trained on six different sized geometries and then tested on three uniquely sized geometries. This study revealed that the 1-st order additive conductivity model performed the best. A new feature of void formation scaling was implemented to more realistically capture the heat transfer as voids are created. A scaling term was added to the model to activate the conductivity correction as decomposition progresses.

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

An Overview of the Design and Optimized Operation of Vanadium Redox Flow Batteries for Durations in the Range of 4–24 Hours

An extensive review of modeling approaches used to simulate vanadium redox flow battery (VRFB) performance is conducted in this study. Material development is reviewed, and opportunities for additional development identified. Various crossover mechanisms for the vanadium species are reviewed, and their effects on its state of charge and its state of health assessed. A stack design focusing on flow fields and an electrode design tailored to various flow fields are reviewed. An operational strategy that takes these parameters into account is reviewed for various operating envelopes, chosen based on end user preference in terms of minimizing capital cost or operation and maintenance cost. This work provides a framework for the design and operation of a VRFB for various grid services.

25 ENERGY STORAGE↗
DOE OSTI2022

Efficiency Studies of Fast Neutron Tracking Using MCNP

Fast neutron identification and spectroscopy is of great interest to nuclear physics experiments. Using the neutron elastic scattering, the fast neutron momentum can be measured. Wang and Morris introduced the theoretical concept that the initial fast neutron momentum can be derived from up to three consecutive elastic collisions between the neutron and the target, including the information of two consecutive recoil ion tracks and the vertex position of the third collision or two consecutive elastic collisions with the timing information. Here, we also include the additional possibility of measuring the deposited energies from the recoil ions. In this paper, we simulate the neutron elastic scattering using the Monte Carlo N-Particle Transport Code (MCNP) and study the corresponding neutron detection and tracking efficiency. The corresponding efficiency and the scattering distances are simulated with different target materials, especially natural silicon (92.23%28Si, 4.67%29Si, and 3.1%30Si) and helium-4 (4He). The timing of collision and the recoil ion energy are also investigated, which are important characters for the detector design. We also calculate the ion traveling range for different energies using the software, “The Stopping and Range of Ions in Matter (SRIM)”, showing that the ion track can be most conveniently observed in 4He unless sub-micron spatial resolution can be obtained in silicon.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Design and Optimization of a Hybrid Design for Quantum Transduction

This study presents the mechanical design and analysis of a quantum electro-optical transducer engineered to operate at millikelvin temperatures within a dilution refrigerator. The transducer enables bidirectional microwave-optical frequency conversion through a hybrid architecture that integrates a superconducting radiofrequency (SRF) cavity with an electro-optic optical cavity. Among several design options investigated, the configuration offering the best thermal and mechanical performance was selected, yielding a robust solution with reduced sensitivity to fabrication tolerances, improved heat dissipation, as well as alignment precision. The design ensures uniform temperature distribution, enabling higher laser pump powers and, thus, increased conversion efficiency, while maintaining mechanical stresses safely below the material yield strength. Electromagnetic simulations further validate the design, demonstrating enhanced coupling between the optical and microwave modes, as well as a broader tuning range achieved with smaller tuner displacements.

quantum computing↗
DOE OSTI2023

Implementation and experimental validation of nonlocal damage in a large-strain elasto-viscoplastic FFT-based framework for predicting ductile fracture in 3D polycrystalline materials

Ductile materials, such as metal alloys, can undergo substantial deformation before failure. Additionally, these materials are usually of polycrystalline composition and exhibit strongly anisotropic behavior at small length scales. Previously developed fast Fourier transform (FFT)-based models can model ductile fracture of isotropic materials or the elastic–plastic behavior of anisotropic polycrystalline materials; however, there remains a need to couple both capabilities. This work extends a large-strain FFT-based crystal plasticity model to simulate ductile fracture of polycrystalline materials. In this work, a triaxiality-based continuum damage mechanics (CDM) formulation is incorporated into a large-strain elasto-viscoplastic FFT (LS-EVPFFT) framework. The CDM formulation is augmented with an integral-based nonlocal regularization approach that correctly handles gas-phase material necessary to model unconstrained surfaces. To validate the damage-enabled LS-EVPFFT framework, mesoscale copper tensile coupons were machined using microwire electrical discharge machining and experimentally characterized using electron backscatter diffraction. In-situ optical digital image correlation was performed during uniaxial testing to provide a side-by-side comparison of the experimental and computational strain fields and stress–strain responses. The damage-enabled LS-EVPFFT framework can simulate the complete macroscopic stress–strain response of ductile polycrystals to failure. The model reproduces necking behavior that qualitatively agrees with experimental observations. By leveraging the relatively low computational cost of the damage-enabled LS-EVPFFT framework, the framework presented here allows the ductile fracture response of 3D polycrystalline materials to be tractably predicted.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Scoping Study to Develop a Computational Fluid Dynamics Based Model to Predict Radiological Materials Packaging Temperatures within a Generic Staging Building

The objective of this work is to perform computational fluid dynamics simulations of a ventilated radiological-material-package staging building to determine the effect of including or excluding a variety of physical effects and computational methods on both the predicted package temperatures and required computational resources. The generic building contains 640 drum-packages containing heat-generating, radiological material supported on four racks that are eight levels tall. Additionally, there is a forced ventilation system, lighting, and insulated walls. Computational models were constructed that included or excluded (a) shelving, (b) effects of unsteadiness, and (c) radiation heat transfer. Simulations with each drum modeled separately were compared to simpler simulations with sets of four drums represented by an equivalent box-package. These models employed between 106 to 107 elements. Steady state simulations predicted package temperatures that were within 0.1°C of simulations that included transient effects and required only one-eighth the computational resources. Calculations that excluded shelving predicted temperatures within 0.6°C of simulations that included shelving and required one-fourth the computational resources. Excluding radiation heat transfer systematically increased temperatures by around 1.5°C but reduced computational resources by a factor of four. The box-package model reduced the computational resources by a factor of 3, but systematically predicted higher temperatures by around 1.1°C. These results will be used to develop an efficient computational fluid dynamics model to assess the ability of different staging building designs to prevent the temperature of package components from exceeding specified limits.

42 ENGINEERING↗
DOE OSTI2021

Impact of Materials Properties on Higher-Temperature Engine Operation

Here, we examine the effects on materials temperatures and engine efficiency via simulations of engines operating at temperatures which exceed the thermal limits of today’s materials. Potential focus areas include high-speed, high-load operation (in the fuel-enrichment zone) as well as conditions of selective cooling at lower speeds and loads. We focus on a light-duty DISI and a heavy-duty CI engine using GT-Power. Temperature distributions within the head, block, piston, and valves were obtained from 3D FEA simulations coupled with 1D GT-Power representations of the engine’s gas flow and combustion regions. We use experimentally measured thermal properties of current commercial alloys for specific engine components, as well as candidate developmental alloys with improved temperature tolerance, to gauge the effects of materials properties on engine performance, particularly focusing on operating areas where materials which can withstand higher temperatures can enable intensified combustion conditions and increased specific output. The outcome of this work is guidance on materials selection targets and benefits as OEMs meet engine-performance needs in the coming decades, to include Class 8 line-haul freight vehicles with reduced emissions, as well as design of future near net zero-carbon light-duty range extenders for more rapid electrification of medium- and heavy-duty vehicles.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Nuclear magnetic resonance and molecular simulation study of H 2 and CH 4 adsorption onto shale and sandstone for hydrogen geological storage

Understanding pure H 2 and H 2 /CH 4 adsorption and diffusion in earth materials is one vital step toward a successful and safe H 2 storage in depleted gas reservoirs. Despite recent research efforts such understanding is far from complete. In this work we first use Nuclear Magnetic Resonance (NMR) experiments to study the NMR response of injected H 2 into Duvernay shale and Berea sandstone samples, representing materials in confining and storage zones. Then we use molecular simulations to investigate H 2 /CH 4 competitive adsorption and diffusion in kerogen, a common component of shale. Our results indicate that in shale there are two H 2 populations, i.e., free H 2 and adsorbed H 2 , that yield very distinct NMR responses. However, only free gas presents in sandstone that yields a H 2 NMR response similar to that of bulk H 2 . About 10 % of injected H 2 can be lost due to adsorption/desorption hysteresis in shale, and no H 2 loss (no hysteresis) is observed in sandstone. Here, our molecular simulation results support our NMR results that there are two H 2 populations in nanoporous materials (kerogen). The simulation results also indicate that CH 4 outcompetes H 2 in adsorption onto kerogen, due to stronger CH 4 -kerogen interactions than H 2 -kerogen interactions. Nevertheless, in a depleted gas reservoir with low CH 4 gas pressure, about ~30 % of residual CH 4 can be desorbed upon H 2 injection. The simulation results also predict that H 2 diffusion in porous kerogen is about one order of magnitude higher than that of CH 4 and CO 2 . This work provides an understanding of H 2 /CH 4 behaviors in deleted gas reservoirs upon H 2 injection and predictions of H 2 loss and CH 4 desorption in H 2 storage.

08 HYDROGEN↗
DOE OSTI2022

Simulation for the fault current limiting operation of REBCO CORC superconducting cables with different core materials

With high power density and low loss, CORC superconducting cables composed of RE- Ba 2 Cu 3 O 7–δ (REBCO) coated conductors are of great interest for power transmission applications. They can also effectively be designed to serve as fault current limiters, thanks to their sharp superconducting-to-normal transition. Coupled electromagnetic-thermal finite element (FE) simulations implemented in COMSOL Multiphysics were developed and validated against the published results for the fault current limiting performance of two CORC cables of REBCO coasted conductors wound on either copper or stainless steel cores. For improved accuracy, temperature dependencies of electrical and thermal properties of all component materials were considered in the simulations. The simulations were performed in the cross-sections of the cables for every individual layer of each REBCO tape to deliver comprehensive understanding of evolution of the current distribution and temperature rise in those layers. Furthermore, the studies can suggest approaches to optimize cable design for fault current limiting (FCL) applications by assessing the role of individual components. In the paper, possible computational errors, challenges and improvements are also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Analysis of Design Alternatives of Actively Cooled RF Window for MPEX

The Materials Plasma Exposure eXperiment (MPEX) has been designed as a linear plasma divertor simulator in order to address plasma material interaction (PMI) science for next-generation fusion devices. It will have the capability to test neutron irradiated samples with plasma fluxes of greater than 10 24 m -2 s -1 . It is expected to operate steady state for up to 10 6 s to consider PMI affects through reactor end of life. The conceptual design of MPEX was completed in 2019, with preliminary design having begun in 2020. The plasma source for MPEX is a helicon antenna, where the energized helical antenna sits outside of the vacuum in order to minimize impurities in the plasma. It is expected to receive up to 200 kW of continuous power, and so the antenna and the window must be actively cooled. The water-cooled copper antenna has been operated at full power on the Proto-MPEX device (which is a test facility to demonstrate the plasma source and heating systems). The water-cooled window, however, is a novel component that must meet numerous competing requirements. It requires a low dielectric loss to allow the Radio Frequency (RF) power to create the plasma within the vacuum boundary. It must be structurally robust to handle the significant heat flux from the plasma and any heat from dielectric coupling. It must be compatible with the coolant (preferably water). It requires a vacuum seal that minimizes impurities into the plasma and does not compromise the structural integrity of the window. Two window designs have been tested. Results from these tests, where temperatures are measured and heat fluxes inferred from infrared camera data, have been correlated with thermal-structural simulations. When these simulations are extrapolated to the full power steady-state heat fluxes that are expected in MPEX, the designs do not appear to have the necessary structural robustness. This study explores design alternatives for the MPEX helicon antenna window, presents analysis results for several of the alternatives, and shows a viable solution that satisfies the requirements for MPEX operation.

Tipton Jr, Joseph↗
DOE OSTI2022

Newton trust-region methods with primary variable switching for simulating high temperature multiphase porous media flow

Coupling multiphase flow with energy transport due to high temperature heat sources introduces significant new challenges since boiling and condensation processes can lead to dry-out conditions with subsequent re-wetting. The transition between two-phase and single-phase behavior can require changes to the primary dependent variables adding discontinuities as well as extending constitutive nonlinear relations to extreme physical conditions. Practical simulations of large-scale engineered domains lead to Jacobian systems with a very large number of unknowns that must be solved efficiently using iterative methods in parallel on high-performance computers. Performance assessment of potential nuclear repositories, carbon sequestration sites and geothermal reservoirs can require numerous Monte-Carlo simulations to explore uncertainty in material properties, boundary conditions, and failure scenarios. Due to the numerical challenges, standard NR iteration may not converge over the range of required simulations and require more sophisticated optimization method like trust-region. In this study, we use the open-source simulator PFLOTRAN for the important practical problem of the safety assessment of future nuclear waste repositories in the U.S. DOE geologic disposal safety assessment Framework. The simulator applies the PETSc parallel framework and a backward Euler, finite volume discretization. We demonstrate failure of the conventional NR method and the success of trust-region modifications to Newton’s method for a series of test problems of increasing complexity. Trust-region methods essentially modify the Newton step size and direction under some circumstances where the standard NR iteration can cause the solution to diverge or oscillate. Furthermore, we show how the Newton Trust-Region method can be adapted for Primary Variable Switching (PVS) when the multiphase state changes due to boiling or condensation. The simulations with high-temperature heat sources which led to extreme nonlinear processes with many state changes in the domain did not converge with NR, but they do complete successfully with the trust-region methods modified for PVS. This implementation effectively decreased weeks of simulation time needing manual adjustments to complete a simulation down to a day. Finally, we show the strong scalability of the methods on a single node and multiple nodes in an HPC cluster.

54 ENVIRONMENTAL SCIENCES↗