Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large system simulation”

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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Experimental observations of detached bow shock formation in the interaction of a laser-produced plasma with a magnetized obstacle

The magnetic field produced by planets with active dynamos, like the Earth, can exert sufficient pressure to oppose supersonic stellar wind plasmas, leading to the formation of a standing bow shock upstream of the magnetopause, or pressure-balance surface. Scaled laboratory experiments studying the interaction of an inflowing solar wind analog with a strong, external magnetic field are a promising new way to study magnetospheric physics and to complement existing models, although reaching regimes favorable for magnetized shock formation is experimentally challenging. This paper presents experimental evidence of the formation of a magnetized bow shock in the interaction of a supersonic, super-Alfvenic plasma with a strongly magnetized obstacle at the OMEGA laser facility. The solar wind analog is generated by the collision and subsequent expansion of two counterpropagating, laser-driven plasma plumes. The magnetized obstacle is a thin wire, driven with strong electrical currents. Hydrodynamic simulations using the FLASH code predict that the colliding plasma source meets the criteria for bow shock formation. Spatially resolved, optical Thomson scattering measures the electron number density, and optical emission lines provide a measurement of the plasma temperature, from which we infer the presence of a fast magnetosonic shock far upstream of the obstacle. Proton images provide a measure of large-scale features in the magnetic field topology, and reconstructed path-integrated magnetic field maps from these images suggest the formation of a bow shock upstream of the wire and as a transient magnetopause. Finally, we compare features in the reconstructed fields to two-dimensional MHD simulations of the system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Local ordering in Ge/Ge–Sn semiconductor alloy core/shell nanowires revealed by extended x-ray absorption fine structure (EXAFS)

Short-range atomic order in semiconductor alloys is a relatively unexplored topic that may promote design of new materials with unexpected properties. Here, local atomic ordering is investigated in Ge–Sn alloys, a group-IV system that is attractive for its enhanced optoelectronic properties achievable via a direct gap for Sn concentrations exceeding ≈10 at. %. The substantial misfit strain imposed on Ge–Sn thin films during growth on bulk Si or Ge substrates can induce defect formation; however, misfit strain can be accommodated by growing Ge–Sn alloy films on Ge nanowires, which effectively act as elastically compliant substrates. In this work, Ge core/Ge 1–x Sn x (x ≈ 0.1) shell nanowires were characterized with extended x-ray absorption fine structure (EXAFS) to elucidate their local atomic environment. Simultaneous fitting of high-quality EXAFS data collected at both the Ge K-edge and the Sn K-edge reveals a large (≈ 40%) deficiency of Sn in the first coordination shell around a Sn atom relative to a random alloy, thereby providing the first direct experimental evidence of significant short-range order in this semiconductor alloy system. Furthermore, comparison of path length data from the EXAFS measurements with density functional theory simulations provides alloy atomic structures consistent with this conclusion.

36 MATERIALS SCIENCE↗
DOE OSTI2024

DL_POLY Quantum 2.0: A modular general-purpose software for advanced path integral simulations

DL_POLY Quantum 2.0, a vastly expanded software based on DL_POLY Classic 1.10, is a highly parallelized computational suite written in FORTRAN77 with a modular structure for incorporating nuclear quantum effects into large-scale/long-time molecular dynamics simulations. This is achieved by presenting users with a wide selection of state-of-the-art dynamics methods that utilize the isomorphism between a classical ring polymer and Feynman’s path integral formalism of quantum mechanics. Here, the flexible and user-friendly input/output handling system allows the control of methodology, integration schemes, and thermostatting. DL_POLY Quantum is equipped with a module specifically assigned for calculating correlation functions and printing out the values for sought-after quantities, such as dipole moments and center-of-mass velocities, with packaged tools for calculating infrared absorption spectra and diffusion coefficients.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Roadmap for warm dense matter physics

This roadmap presents the state-of-the-art, current challenges and near future developments anticipated in the thriving field of warm dense matter (WDM) physics. Originating from strongly coupled plasma physics, high pressure physics and high energy density science, the WDM physics community has recently taken a giant leap forward. This is due to spectacular developments in laser technology, diagnostic capabilities, and computer simulation techniques. Only in the last decade has it become possible to perform accurate enough simulations & experiments to truly verify theoretical results as well as to reliably design experiments based on predictions. Consequently, this roadmap discusses recent developments of and contemporary challenges for theoretical methods and experimental techniques needed to describe, create and diagnose WDM. A large part of this roadmap is dedicated to specific WDM systems and applications in astrophysics, inertial confinement fusion and novel material synthesis.

dense astrophysical objects↗
DOE OSTI2024

Limits to forests-based mitigation in integrated assessment modelling: global potentials and impacts under constraining factors

Forests-based measures such as afforestation/reforestation (A/R) and reducing deforestation (RDF) are considered promising options to mitigate climate change, yet their mitigation potentials are limited by economic and biophysical factors that are largely uncertain. The range of mitigation potential estimates from integrated assessment models raises concerns about the capacity of land systems to provide realistic, cost-effective and permanent land-based mitigation. We use the Global Change Analysis Model to quantify the economic mitigation potential of forests-based measures by simulating a climate policy including a tax on greenhouse gas emissions from agriculture, forestry, and other land uses. In addition, we assess how constraining unused arable land (UAL) availability, forestland expansion rates, and global bioenergy demand may influence the forests-based mitigation potential by simulating scenarios with alternative combinations of constraints. Results show that the average forests-based mitigation potential in 2020–2050 increases from 738 MtCO 2 .yr -1 through a forestland increase of 86 Mha in the fully constrained scenario to 1394 MtCO 2 .yr -1 through a forestland increase of 146 Mha when all constraints are relaxed. Regional potentials in terms of A/R and RDF differ strongly between scenarios: unconstrained forest expansion rates mostly increase A/R potentials in northern regions (e.g., +120 MtCO 2 .yr -1 in North America); while unconstrained UAL conversion and low bioenergy demand mostly increase RDF potentials in tropical regions (e.g., +76 and +68 MtCO 2 .yr -1 in Southeast Asia, respectively). This study shows that forests-based mitigation is limited by many factors that constrain the rates of land use change across regions. These factors, often overlooked in modelling exercises, should be carefully addressed for understanding the role of forests in global climate mitigation and defining pledges towards the Paris Agreement.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Multi-Neutrino Entanglement and Correlations in Dense Neutrino Systems

The time-evolution of multi-neutrino entanglement and correlations are studied in two-flavor collective neutrino oscillations, relevant for dense neutrino environments, building upon previous works. Specifically, simulations performed of systems with up to 12 neutrinos using Quantinuum’s H1-1 20 qubit trapped-ion quantum computer are used to compute n-tangles, and two- and three- body correlations, probing beyond mean-field descriptions. n-tangle re-scalings are found to converge for large system sizes, signaling the presence of genuine multi-neutrino entanglement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Supo Thermal Model Development IV

This report describes the continuation of the Computational Fluid Dynamics (CFD) model of the Supo cooling system described in the reports, Supo Thermal Model Development, II, and III. In the previous reports, the solution natural convection of Supo was simulated using a 2D, axisymmetric model that altered the actual geometry. A 3D geometry of Supo will more accurately estimate the natural convection heat transfer coefficient (HTC) of the system since Supo contains several non-symmetric internal components that can only be modeled in 3D. This, in turn, will improve the natural convection HTC results for comparison purposes with past operational data and previous 2D CFD results. Unfortunately, the complexities of performing multiphase 3D CFD simulations has proven too large of a task for our current capabilities. Instead, this report will focus on the 3D model development for the 2.39 kW/L (31.3 kW) case and limitations and challenges associated with the 3D simulations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

2021 NIRT Mini Drill (After Action Report)

During the summer and fall of 2021, several functional area drills were held that focused on exercising Consequence Management’s (CM) ability to extract and use data from RadResponder for the purpose of answering intermediate-phase questions presented as technical inject requests for information (RFI) in Sandia National Laboratories (SNL) Consequence Management Operational System (COSMOS) software. The scenario chosen was that of Northern Lights 2016 (NL16) which was a large-scale nuclear power plant (NPP) release exercise in the state of Minnesota. The NL16 data was extracted from the Radiological Assessment and Monitoring System (RAMS) event where it was created and was reformatted for implanting to a new RadResponder event. Next, the beta-version of a laboratory sample data simulator was used to generate more sample data that was injected to the event. Five “mini-drills” were devised with each prompt defined by a data-based need. For each drill, a team of assessment and NARAC scientists worked the problem using the drill prompt and the available data in RadResponder. The teams held a kickoff meeting, had several days to work the problem, and then reported their results as well as observations in a hotwash. Several areas for improvement in both the software and process were identified during the course of these drills. This report will document the process of addressing each RFI and the discovered gaps in both software capability and methodology so that they can be considered for future development and investment by the CM and NIRT programs.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2023

Hybrid physics-based and machine learning tools for materials assessment

In this work, we develop novel physics-based and machine learning computational techniques to predict fundamental properties of metallic systems that affect radiation damage behavior in reactor structural materials. Oftentimes, atomistic predictions of engineering alloys simplify chemical compositions to a single element to reduce computational cost and complexity, introducing large sources of uncertainty and potentially missing important behavior. In addition, engineering alloys such as SS 316 are particularly challenging to simulate with first principles methods because of the additional degrees of freedom introduced by magnetic natures of the constituent elements, and very little data of this type exists within the literature. These novel methods aim to improve the qualitative prediction of radiation damage in engineering alloys by more accurately simulating their compositions, both by accelerating the computations and by developing novel analyses.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Evaluation of E3SM land model snow simulations over the western United States

Abstract. Seasonal snow has crucial impacts on climate, ecosystems, and humans, but it is vulnerable to global warming. The land component (ELM) of the Energy Exascale Earth System Model (E3SM) mechanistically simulates snow processes from accumulation, canopy interception, compaction, and snow aging to melt. Although high-quality field measurements, remote sensing snow products, and data assimilation products with high spatio-temporal resolution are available, there has been no systematic evaluation of the snow properties and phenology in ELM. This study comprehensively evaluates ELM snow simulations over the western United States at 0.125∘ resolution during 2001–2019 using the Snow Telemetry (SNOTEL) in situ networks, MODIS remote sensing products (i.e., MCD43 surface albedo product), the spatially and temporally complete (STC) snow-covered area and grain size (MODSCAG) and MODIS dust and radiative forcing in snow (MODDRFS) products (STC-MODSCAG/STC-MODDRFS), and the snow property inversion from remote sensing (SPIReS) product and two data assimilation products of snow water equivalent and snow depth – i.e., University of Arizona (UA) and SNOw Data Assimilation System (SNODAS). Overall the ELM simulations are consistent with the benchmarking datasets and reproduce the spatio-temporal patterns, interannual variability, and elevation gradients for different snow properties including snow cover fraction (fsno), surface albedo (αsur) over snow cover regions, snow water equivalent (SWE), and snow depth (Dsno). However, there are large biases of fsno with dense forest cover and αsur in the Rocky Mountains and Sierra Nevada in winter, compared to the MODIS products. There are large discrepancies of snow albedo, snow grain size, and light-absorbing particle-induced snow albedo reduction between ELM and the MODIS products, attributed to uncertainties in the aerosol forcing data, snow aging processes in ELM, and remote sensing retrievals. Against UA and SNODAS, ELM has a mean bias of −20.7 mm (−35.9 %) and −20.4 mm (−35.5 %), respectively, for spring, and −13.8 mm (−27.8 %) and −10.2 mm (−22.2 %), respectively, for winter. ELM shows a relatively high correlation with SNOTEL SWE, with mean correlation coefficients of 0.69 but negative mean biases of −122.7 mm. Compared to the snow phenology of STC-MODSCAG and SPIReS, ELM shows delayed snow accumulation onset dates by 17.3 and 12.4 d, earlier snow end dates by 35.5 and 26.8 d, and shorter snow durations by 52.9 and 39.5 d, respectively. This study underscores the need for diagnosing model biases and improving ELM representations of snow properties and snow phenology in mountainous areas for more credible simulation and future projection of mountain snowpack.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Initial Engineering Framework for Species Tracking in MSRs Involving Fuel Salt and Structures

This report implements a high-fidelity multiphysics modeling framework using the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program tools to track leaching and plating of materials between the fuel salt and structures in Molten Salt Reactors (MSRs). Specifically, this framework is demonstrated in modeling various processes within MSRs including hot and cold leg corrosion, redox potential driven corrosion, colloid precipitation and deposition, and noble metal plating and decay heat modeling. The model integrates neutronics, thermal-hydraulics, depletion, and thermochemistry to simulate the evolving chemical behavior of the fuel salt in interactions with structural materials. Initial findings suggest that hot-cold leg corrosion in MSR systems may become significant if the redox potential of the fuel salt is not sufficiently controlled. Additionally, the impact of noble metal decay heat deposition in large power reactors is explored. The resulting engineering framework for species tracking will be applied in the future toward specific validation work with ongoing experimental efforts focusing on both chloride and fluoride salt loop experiments with in-situ corrosion and redox control instrumentation.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Statistical Performance of Forced Oscillation Detectors in the Presence of Missing Measurements

In bulk power systems, measurement-based monitoring for large oscillations can help maintain system reliability. One of the challenges encountered in a recent field demonstration was the unavailability of measurements due to underlying measurement quality or communication problems. During the demonstration, the oscillation detector ignored a measurement location if even 10 seconds of data was missing. To extend the detector's ability to operate in these conditions, this paper evaluates the impact of three methods for addressing missing data. The strengths and weaknesses of each approach are evaluated using theoretical expressions for the probability of detection along with results from simulated data and publicly available field measurements. Based on these results, a suitable approach is identified that can extend the oscillation detector's performance when large segments of data are missing.

Follum, James D.↗
DOE OSTI2021

COVID 19 vaccine distribution solution to the last mile challenge: Experimental and simulation studies of ultra-low temperature refrigeration system

Most COVID-19 vaccines require ambient temperature control for transportation and storage. Both Pfizer and Moderna vaccines are based on mRNA and lipid nanoparticles requiring low temperature storage. The Pfizer vaccine requires ultra-low temperature storage (between -80 °C and -60 °C), while the Moderna vaccine requires -30 °C storage. Pfizer has designed a reusable package for transportation and storage that can keep the vaccine at the target temperature for 10 days. However, the last stage of distribution is quite challenging, especially for rural or suburban areas, where local towns, pharmacy chains and hospitals may not have the infrastructure required to store the vaccine. Also, the need for a large amount of ultra-low temperature refrigeration equipment in a short time period creates tremendous pressure on the equipment suppliers. In addition, there is limited data available to address ancillary challenges of the distribution framework for both transportation and storage stages. As such, there is a need for a quick, effective, secure, and safe solution to mitigate the challenges faced by vaccine distribution logistics. The study proposes an effective, secure, and safe ultra-low temperature refrigeration solution to resolve the vaccine distribution last mile challenge. Furthermore, the approach is to utilize commercially available products, such as refrigeration container units, and retrofit them to meet the vaccine storage temperature requirement. Both experimental and simulation studies are conducted to evaluate the technical merits of this solution with the ability to control temperature at -30 °C or -70 °C as part of the last mile supply chain for vaccine candidates.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Gear Test Assembly – Report Experimental Testing and Gear Analysis (FY2021 Midyear Report)

The Gear Test Assembly (GTA) has completed three experimental campaigns at the Mechanisms Engineering Test Loop (METL) facility. The most recent campaign tested Inconel 718 spur gears, tapered roller bearings of 52100 bearing steel with no heat treatment, and a lantern ring shaft seal with spring-loaded fasteners. Testing was performed in Test Vessel 1 (TV1) with 250°C sodium that had an oxide concentration of <10ppm. Testing began in February 2021 and was ended in March 2021 when a thrust bearing assembly failed, requiring shutdown and maintenance. The third experimental campaign completed 1,568 simulated fuel assembly maneuvers before the thrust bearing failure. The GTA was removed from TV1 and cleaned in the Carbonation System. The GTA is currently being prepared for the next experimental campaign. Nondestructive examination methods have been developed at Argonne to monitor the health of the GTA spur gears. Eddy current testing (ECT) and ultrasonic testing (UT) were performed prior to sodium testing and after each experimental campaign. The NDE has shown that the primary damage mechanism is mechanical wear on the gear tooth surface. Several larger nicks are observable on the gear faces, with additional large nicks present on the edges of the gear teeth. While this damage is observable, the overall health of the gears is adequate for more in-sodium testing. A total of 12,752 simulated fuel assembly maneuvers have been completed using this set of gears.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Modeling and Analysis of a Novel Offshore Binary Species Free-Floating Longline Macroalgal Farming System

The investigation of innovative macroalgal cultivation is important and needed to optimize farming operations, increase biomass production, reduce the impact on the ecosystem, and lower system and operational costs. However, most macroalgal farming systems (MFSs) are stationary, which need to occupy a substantial coastal area, require extensive investment in farm infrastructure, and cost high fertilizer and anchoring expenses. This study aims to model, analyze, and support a novel binary species free-floating longline macroalgal cultivation concept. The expected outcomes could provide a basis for the design and application of the novel MFS to improve biomass production, decrease costs, and reduce the impact on the local ecosystem. In this paper, Saccharina latissima and Nereocystis luetkeana were modeled and validated, and coupled with longline to simulate the binary species MFS free float in various growth periods and associated locations along the US west coast. Further, the numerical predictions indicated the possibility of failure on the longline and breakage at the kelp holdfasts is low. However, the large forces due to an instantaneous change in dynamic loads caused by loss of hydrostatic buoyancy when the longline stretches out of the water would damage the kelps. Buoy-longline contact interactions could damage the buoy, resulting in the loss of the system by sinking. Furthermore, the kelp-longline and kelp-kelp entanglements could potentially cause kelp damage.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Simulation of Liquid Lithium Divertor Geometry Using SOLPS-ITER

Here, plasma-facing component (PFC) geometries are evaluated for a Fusion Nuclear Science Facility (FNSF) design to find solutions that are compatible with flowing liquid lithium components while satisfying requirements from the perspective of plasma and neutral particle transport. Flowing liquid metal (LM) divertor systems offer important advantages due to continual regeneration of the surface material, but may require modifications from standard optimized tokamak divertor designs. Scrape Off Layer Plasma Simulator for the International Thermonuclear Experimental Reactor (SOLPS-ITER) simulations are used to compare standard vertical target, open, and balanced baffled geometries. It is found that the open geometry offers little control of neutral particles and has narrow operating windows where the divertor fluxes can be lowered to acceptable levels, while maintaining the required upstream density. The addition of neutral baffles allows greater control over the divertor conditions with less impact on the upstream density. The divertor neutral pressure in the baffled geometry is slightly higher than the open divertor, but not as large with a vertical target. A first-pass coupling of the SOLPS solution with sheath and implantation/erosion models to compute the lithium emission shows that the lithium is largely confined to the near-surface region due to the small ionization mean free path and large main ion flow toward the surface. The lithium source was varied by three orders of magnitude. At the highest rate, when the sourced lithium approaches the level of the main ion recycling flux, a high core lithium concentration (~10%) is observed with only moderate power dissipation to affect the divertor plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Rotation‐Controlled Diurnal Evolution of Uranus' Asymmetric Bow Shock at Equinox

Uranus possesses the most extreme magnetic and rotational geometry in the solar system, resulting in a uniquely dynamic and asymmetric interaction between its magnetosphere and the solar wind. Here we investigate the diurnal evolution of the Uranian bow shock (BS) at equinox using global multifluid magnetohydrodynamic simulations constrained by Voyager 2 observations. The BS structure and variability are quantified using the stand‐off distance, terminator distance, flaring parameter, and BS parameter, enabling a systematic assessment of its global asymmetry over one planetary rotation. Our results show that the BS exhibits pronounced diurnal asymmetry that is primarily controlled by planetary rotation. Even under steady upstream solar wind conditions, the BS undergoes periodic expansion and contraction, reflecting rotation‐driven reconfiguration of the planetary magnetospheric topology. This diurnal modulation represents a remarkably large asymmetry in planetary BS geometry within the solar system. These results establish planetary rotation as the intrinsic primary driver of Uranus' BS dynamics at equinox and provide essential constraints for the design and formulation of future space missions to the ice giants. These findings also offer a unique framework for understanding shock physics in the abundant population of ice‐giant exoplanets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Preventing Electrolyte Decomposition on a Ca Metal Electrode Interface Using an Artificial Solid‐Electrolyte Interphase

Abstract Calcium ion batteries are gaining attention as alternatives to lithium‐ion technology because they offer comparable properties at reduced cost and improved safety. However, progress has been limited because of the inability to efficiently and reversibly plate and strip Ca metal anodes in organic electrolytes. Moreover, the inorganic components of the solid‐electrolyte interphase (SEI) that form via decomposition of the electrolyte often do not allow for the diffusion of Ca ions. In this work, an approach combining density functional theory and ab initio molecular dynamics (AIMD) simulations is utilized to show that the use of a preformed artificial SEI layer of amorphous can potentially prevent electrolyte decomposition. First, Ca is shown to be able to intercalate into an amorphous layer (up to Ca 1.5 Al 2 O 3 ) and diffuse through on a reasonable time scale. Through calculation of the density of states, the system is found to remain insulating up to the equilibrium stoichiometry. Finally, AIMD simulations with a realistic organic electrolyte environment are used to show that this calcinated layer completely prevents the decomposition of solvent molecules. This approach can provide a route to efficient rechargeable Ca ion batteries, paving the way for cheap large‐scale energy storage.

Young, Joshua↗