Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Results”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Overland flow numerical model prediction, Lower Triangle Region in East River Watershed, Colorado, 3 days

This data package contains numerical simulation results of surface flow variables such as flow velocity and water depth in Lower Triangle Region in East River Watershed, Colorado. The surface flow is a consequence of a high intensity rainfall event with a total duration of 3 days, available at a resolution of 10 minutes. The results are computed on triangular multiresolution meshes with resolutions ranging from 10 meter to 80 meter. The data package also contains a simulation on a uniform triangular mesh with a resolution of 10 meter. The simulations consider surface flow only and neglect subsurface flow, infiltration, and evapotranspiration. The purpose of the data is to assess the quality of a mesh refinement strategy.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A phase-field model for non-isothermal phase transformation and plasticity in polycrystalline yttria-stabilized tetragonal zirconia

Here we propose an elastoplastic phase-field (PF) model to investigate the mechanics of tetragonal-to-monoclinic phase transformation (TMPT) and elastoplastic deformation of polycrystalline yttria-stabilized tetragonal zirconia (YSTZ). A Landau polynomial with non-vanishing chemical energy at the equilibrium temperature is introduced to account for the actual formation energies of the phases. The effects of different grain orientations, latent heat, and temperature on TMPT and deformation mechanisms are considered. The suppressive transformation effects of the grain boundaries (GBs) is modeled using an inhomogeneous kinetic coefficient in the bulk and GBs. The simulation results for single crystals demonstrate the capability of the model to reproduce the orientation-dependent compressive deformation of YSTZ similar to atomistic simulations and micropillar experiments. The single crystal with [100] crystallographic orientation along the loading direction (SC[100]) displays both TMPT and plasticity, SC[101] experiences only phase transformation, while SC[001] undergoes only plastic yielding. The TMPT induced by compressive loading exhibits shape memory effect (SME) below the equilibrium transformation temperature and pseudoelasticity above it, while the critical transformation stress increases with increasing loading temperature. The irrecoverable plastic strain is found to trap a part of the monoclinic phase, which prevents a complete reverse transformation. The polycrystalline cases also display SME and PE at low and high temperatures, respectively. Due to the orientation differences between grains and the stress concentrations at geometric nonlinearities, plastic deformation occurs in polycrystalline YSTZ for an applied load less than the yield stress. The results suggest a possible limitation of plasticity and an improvement of the shape recovery of YSTZ if one can control the orientation of the grains and/or increase the density of stacking faults at the GBs during material processing.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A Model Intercomparison Study of Aerosol‐Cloud‐Turbulence Interactions in a Cloud Chamber: 1. Model Results

This study presents the first model intercomparison of aerosol‐cloud‐turbulence interactions in a controlled cloudy Rayleigh‐Bénard Convection chamber environment, utilizing the Pi Chamber at Michigan Technological University. We analyzed simulated cloud chamber‐averaged statistics of microphysics and thermodynamics in a warm‐phase, cloudy environment under steady‐state conditions at varying aerosol injection rates. Simulation results from seven distinct models (DNS, LES, and a 1D turbulence model) were compared. Our findings demonstrate that while all models qualitatively capture observed trends in droplet number concentration, mean radius, and droplet size distributions at both high and low aerosol injection rates, significant quantitative differences were observed. Notably, droplet number concentrations varied by over two orders of magnitude between models for the same injection rates, indicating sensitivities to the model treatments in droplet activation and removal and wall fluxes. Furthermore, inconsistencies in vertical relative humidity profiles and in achieving steady‐state liquid water content suggest the need for further investigation into the mechanisms driving these variations. Despite these discrepancies, the models generally reproduced consistent power‐law relationships between the microphysical variables. This model intercomparison underscores the importance of controlled cloud chamber experiments for validating and improving cloud microphysical parameterizations. Recommendations for future modeling studies are also highlighted, including constraining wall conditions and processes, investigating droplet/aerosol removal (including sidewall losses), and conducting simplified experiments to isolate specific processes contributing to model divergence and reduce model uncertainties.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

First attempt to quantify W7-X island divertor plasma by local experiment-model comparison

As a complement to our recent work, which focused on understanding the basic detachment physics and general experimental and numerical trends observed in the W7-X island divertor, this paper compares EMC3-Eirene simulation results with different local diagnostics, including IR cameras, Langmuir probes, H α -cameras, and Thomson scattering. The main purposes are to (1) justify the simulation setup in the previous work, (2) identify the application limitations of the current EMC3-Eirene model, (3) verify the consistency of different diagnostics, and (4) isolate the main geometric and physical effects that need to be prioritized in further developing the EMC3-Eirene code and improving diagnostic capabilities. It turns out that the current version of the EMC3-Eirene code (without drifts) is not yet able to quantitatively reproduce all selected local measurements simultaneously under the current experimental conditions (in particular, the existence of error fields). Nevertheless, it can be shown that within a reasonable range of variation in magnetic configuration, cross-field transport, and SOL plasma state in the modeling, a region of overlap between the numerical results and the local measurements can be established. More accurate model-experiment comparisons will require clarification of error fields, implementation of drifts in EMC3, and improvement of diagnostic capabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Critical nuclei at hetero-phase interfaces

Two-step nucleation, in which a metastable intermediate phase acts as a precursor for nucleating a thermodynamically stable phase, has been widely observed in many materials systems and solid-state reactions. Among the advantages of two-step nucleation is that the stable phase may nucleate heterogeneously at the hetero-phase interface between the original and the precursory phases. Although heterogeneous nucleation (HN) theories for homo-phase grain boundaries and inert surfaces are well established, our understanding of HN at reactive hetero-phase interfaces remains incomplete. This deficiency stems from the discontinuity of the chemical potential driving force across the hetero-phase interface, which profoundly affects the fundamental properties of the nucleus in a way that is not properly accounted for in existing models. In this work, we incorporate these effects to extend the classical nucleation theory to HN at hetero-phase interfaces. Our extended model demonstrates that the nucleus shape along the minimum energy path is strongly size-dependent, and this additional degree of freedom can result in the reduction of the critical nucleus volume and associated activation energy barrier by orders of magnitude relative to conventional predictions. The simulation results are used to construct a sensitivity map in the parameter space of interfacial energy and bulk driving force ratios, which quantifies the difference in nucleation barriers predicted by different models.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Heat pump water heater enhanced with phase change materials thermal energy storage: Modeling study

A promising solution to improve the first hour rating (FHR) of a heat pump water heater (HPWH) involves employing a secondary tank which contains phase change material (PCM) capsules. To better understand the influence of PCM thermal storage on the HPWH operational performance, a dynamic model was developed to simulate and analyze the behavior of a newly developed HPWH technology that incorporates PCM storage into a standard HPWH to optimize key parameters such as the uniform energy factor and FHR. Mathematical models of several key components of the proposed HPWH-PCM integrated thermal energy storage (TES) system, e.g., water heater tank, PCM TES tank, evaporator, compressor, and expansion valve, have been elaborated. Also, a model-based control co-simulation platform was developed to integrate a embed PCM storage HPWH dynamic model with a control model for better supporting control design, analysis, verification, and validation. The model accuracy has been validated through comparing simulation results with lab test results, with a mean average percentage error of <5.5% for most of the selected performance variables. In addition, using the developed co-simulation platform, the demand response control strategy was studied to evaluate the load flexibility of the combined HPWH-PCM storage system by shifting system power usage to outside of the 3.5 h peak load period.

25 ENERGY STORAGE↗
DOE OSTI2025

Stress accommodation in nanoscale dolan bridges designed for superconducting qubits

Josephson junctions are the principal circuit element in numerous superconducting quantum information devices and can be readily integrated into large-scale electronics. However, device integration at the wafer scale necessarily depends on having a reliable, high-fidelity, and high-yield fabrication method for creating Josephson junctions. When creating Al/AlO x based superconducting qubits, the standard Josephson junction fabrication method relies on a sub-micron suspended resist bridge, known as a Dolan bridge, which tends to be particularly fragile and can often times fracture during the resist development process, ultimately resulting in device failure. In this work, we demonstrate a unique Josephson junction lithography mask design that incorporates stress-relief channels. Our simulation results show that the addition of stress-relief channels reduces the lateral stress in the Dolan bridge by more than 70% for all the bridge geometries investigated. In practice, our novel mask design significantly increased the survivability of the bridge during device processing, resulting in 100% yield for over 100 Josephson junctions fabricated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Large-scale Multiphysics Simulations of Small Modular Reactors Operating in Natural Circulation

Thanks to the advancements in high-performance computing, advanced modeling and simulation have become crucial in driving the development and deployment of next-generation nuclear reactors, such as small modular reactors (SMRs). SMRs offer the promise of cost-effective baseload electricity production and improved safety, while addressing some of the challenges associated with large reactor designs, such as high capital costs and extended construction timelines. As part of the Exascale Computing Project, the large-scale multiphysics simulation of an entire SMR primary system has been achieved by combining computational fluid dynamics and neutronics. In addition to the successful demonstration of full-core SMR simulations, the current study integrated the impact of natural circulation into the system. Natural circulation is the primary mechanism driving coolant circulation in SMRs. The mass flow rate in the core depends on the core power, and a numerical model has been developed to predict it. The pressure drop caused by the helical coil steam generator was also accounted for by developing a pressure drop correlation based on high-fidelity large eddy simulation results, further improving prediction accuracy. In conclusion, the results of the study demonstrate that the implemented natural circulation model is effective in predicting the responses of SMR full-core multiphysics simulations.

ECP↗
DOE OSTI2020

SAM Enhancements and Model Developments for Molten-Salt-Fueled Reactors

To support the development and utilization of the SAM code for molten-salt-fueled reactor (MSR) safety analysis and licensing, an effort was devoted to enhancing code capabilities and developing reference models for the MSR primary loop. A reference standard problem of a prototypical reactor design is foundational to NRC to verify the adequacy of computer codes and evaluation models for a specific reactor type. A thermo-fluid model of the Molten Salt Reactor Experiment (MSRE) has been developed. The MSRE primary loop model consists of a 2-D core region and external core components in 1-D or 0-D. Both the steady-state and a transient scenario are simulated. The porous medium model is utilized for both the molten salt fluid channels and the moderator matrix in the MSRE core. The use of the porous medium model for the MSRE core is verified with the higher-fidelity simulation results using 1-D representations of the fluid channels and 3-D modeling of core structures. To further enhance SAM multi-scale simulation capabilities for MSRs, a mass transport model is developed and implemented in SAM multi-dimension flow module for modeling of species transport such as delayed neutron precursors in MSRs.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Effects of Different Closure Choices in Core-collapse Supernova Simulations

The two-moment method is widely used to approximate the full neutrino transport equation in core-collapse supernova (CCSN) simulations, and different closures lead to subtle differences in the simulation results. In this paper, we compare the effects of closure choices on various physical quantities in 1D and 2D time-dependent CCSN simulations with our multigroup radiation hydrodynamics code FORNAX. We find that choices of the third-order closure relations influence the time-dependent simulations only slightly. Choices of the second-order closure relation have larger consequences than choices of the third-order closure, but these are still small compared to the remaining variations due to ambiguities in some physical inputs such as the nuclear equation of state. We also find that deviations in Eddington factors are not monotonically related to deviations in physical quantities, which means that simply comparing the Eddington factors does not inform one concerning which closure is better.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Ion temperature anisotropy model with cross‐field drifts in the scrape‐off layer

Abstract An ion temperature anisotropy model including cross‐field drifts has been implemented in the 2D transport code UEDGE. Two ion temperature equations are solved in the model to give separate ion temperatures perpendicular and parallel to the magnetic field, which results in an additional mirror force for ion parallel transport and also leads to an additional curvature drift in the perpendicular direction. In our previous work, an ion temperature anisotropy model was implemented and applied to the full 2D version of UEDGE including recycled neutrals. The simulation results showed that the mirror force due to anisotropy reduces ion parallel flow velocity, however, the simulations were presented without cross‐field electric and magnetic drifts. Nevertheless, in present medium‐size tokamaks or future compact fusion devices, cross‐fields are important. In this work, the anisotropy model is extended to be compatible with the inclusion of cross‐field drifts. Preliminary results on comparison between the model solutions with and without cross‐field drifts are presented.

Zhao, Menglong↗
DOE OSTI2021

Molecular dynamics and experimental study on the adhesion mechanism of polyvinyl alcohol (PVA) fiber in alkali-activated slag/fly ash

This paper aims to study the adhesion mechanism of polyvinyl alcohol (PVA) fiber within alkali-activated slag/fly ash (AASF) matrix using molecular dynamics (MD) simulation in combination with systematic experimental characterization. The adhesion of PVA to C-(N-)A-S-H gel with different Ca/(Si+Al) and Al/Si ratios was modeled using MD simulation, with the related adsorption enthalpy calculated and the adhesion mechanism explored. The experimentally attained chemical bonding energy of PVA fiber in AASF coincides well with the simulation results. In both cases, the adhesion enhances primarily with increasing Ca/(Si+Al) ratio of C-(N-)A-S-H gel. Additionally, MD simulation indicates preferential element distributions of Ca around PVA molecule, which was confirmed experimentally by the detection of the Ca-rich C-(N-)A-S-H gel in the interfacial transition zone (ITZ). }This study provides further insights into the adhesion mechanism of PVA fiber to C-(N-)A-S-H gel formed in AASF, which is particularly valuable for the future development of PVA-based high-performance alkali-activated composites.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The June 2012 North American Derecho: A Testbed for Evaluating Regional and Global Climate Modeling Systems at Cloud‐Resolving Scales

Abstract In this paper, we introduce a testbed for evaluating and comparing climate modeling systems at cloud resolving scales using hindcasts of the June 2012 North American derecho. To demonstrate its utility for model intercomparison, the testbed is applied to two models: the regionally‐refined Simple Cloud‐Resolving E3SM Atmosphere Model (SCREAM) at 6.5, 3.25 and 1.625 km grid spacing and the Weather Research and Forecasting (WRF) model with 3.2 and 1.6 km grid spacing. We find the simulation results to be highly sensitive to the initial conditions (ICs), initialization time, and model configurations, with ICs from the Rapid Refresh producing the best simulation. Significant improvement is identified in both models as horizontal grid spacing is refined. While a propagation delay of approximately 2 hr is found in both models, SCREAM at 1.625 km simulates the observed bow echo structure of the derecho well and predicts strong surface gusts that exceed 30 m/s. In comparison, WRF has difficulty producing surface wind over 25 m/s, with wind gusts in WRF 42%–46% lower than in SCREAM. However, WRF has a lower bias in simulating cloud top temperature and extent, but overestimates precipitation intensity. Both models reproduce the observed outgoing longwave radiation spatial patterns well (Pearson correlation >0.88), but, compared with NEXRAD observations, simulate generally larger areas of composite radar reflectivity >40 dBZ and underestimate the precipitating area by ∼47%.

2012 North American derecho↗
DOE OSTI2020

Design and Analysis of the Staggered Double Grating Slow Wave Circuit for 263 GHz Sheet Beam TWT

The complete design process of a high-frequency circuit for a 263 GHz staggered double grating, sheet beam traveling wave tube is presented in this article. Here, this device, for use in a pulsed electron paramagnetic resonance spectrometer instrument, requires >30 dB gain and a 20 GHz bandwidth centered at 263 GHz. In addition to the standard Pierce gain analysis and point by point particle in cell (PIC) simulation, a novel fast pulse response analysis technique was utilized to determine the accurate synchronization voltage and test the ``hot” feature of the high-frequency circuit. Since this new method can evaluate the entire bandwidth in a single simulation, it can significantly reduce the computational resources normally required for the conventional PIC analysis while also illustrating the potential instabilities near the band edge. Good agreement between the fast pulse response and PIC simulation results validate the accuracy of this novel method. According to the analysis results and given the manufacturing limits of nano/micro-computer numerical control (CNC) machining, an appropriate SWS, coupler, and sever designs were chosen to complement the simulated design. Finally, a prototype circuit was manufactured and cold tested to verify the design and inspect the machining tolerances.

42 ENGINEERING↗
DOE OSTI2021

Monte Carlo Simulation of Background and Source Measurements with CSG and CAD Geometries

Detecting radioactive sources in an urban environment is difficult due to the large magnitude and variability of the background radiation. To support the search mission of the National Nuclear Security Administration, a project was undertaken to determine if first-principles modeling and simulation could be used to accurately predict the response of radiation detectors in urban environments. This study demonstrates that the simulated detector responses to photon radiation from both an isotopic source and the background compare well to benchmark-quality measurements in a large urban-like test environment. Simulation results using a traditional combinatorial solid geometry and a CAD model generated from LiDAR measurements gave very similar detector responses. With this validation, further simulations could be used to predict detector responses for various situations in real urban settings.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Cloud droplet diffusional growth in homogeneous isotropic turbulence: bin microphysics versus Lagrangian super-droplet simulations

The increase in the spectral width of an initially monodisperse population of cloud droplets in homogeneous isotropic turbulence is investigated by applying a finite-difference fluid flow model combined with either Eulerian bin microphysics or a Lagrangian particle-based scheme. The turbulence is forced applying a variant of the so-called linear forcing method that maintains the mean turbulent kinetic energy (TKE) and the TKE partitioning between velocity components. The latter is important for maintaining the quasi-steady forcing of the supersaturation fluctuations that drive the increase in the spectral width. We apply a large computational domain (64 3 m 3 ), one of the domains considered in Thomas et al. (2020). The simulations apply 1 m grid length and are in the spirit of the implicit large eddy simulation (ILES), that is, with small-scale dissipation provided by the model numerics. This is in contrast to the scaled-up direct numerical simulation (DNS) applied in Thomas et al. (2020). Two TKE intensities and three different droplet concentrations are considered. Analytic solutions derived in Sardina et al. (2015), valid for the case when the turbulence integral timescale is much larger than the droplet phase relaxation timescale, are used to guide the comparison between the two microphysics simulation techniques. The Lagrangian approach reproduces the scalings relatively well. Representing the spectral width increase in time is more challenging for the bin microphysics because appropriately high resolution in the bin space is needed. The bin width of 0.5 µm is only sufficient for the lowest droplet concentration (26 cm -3 ). For the highest droplet concentration (650 cm -3 ), an order of magnitude smaller bin size is barely sufficient. The scalings are not expected to be valid for the lowest droplet concentration and the high-TKE case, and the two microphysics schemes represent similar departures. Finally, because the fluid flow is the same for all simulations featuring either low or high TKE, one can compare point-by-point simulation results. Such a comparison shows very close temperature and water vapor point-by-point values across the computational domain and larger differences between simulated mean droplet radii and spectral width. The latter are explained by fundamental differences in the two simulation methodologies, numerical diffusion in the Eulerian bin approach and a relatively small number of Lagrangian particles that are used in the particle-based microphysics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A software‐defined networks ‐based measurement method of network traffic for 6G technologies

Abstract In software‐defined networks (SDN) for 6G technologies, the controller should accurately and efficiently measure the flow traffic in switches for traffic engineering. Fine‐grained flow measurement can more accurately describe the traffic in the network, but it also consumes much more resources. To reduce the overhead incurred in the measurement process and obtain the approximate fine‐grained measurements, we propose a novel lightweight measurement scheme that runs in the controller. The novel lightweight measurement architecture consists of two parts: coarse‐grained measurement and interpolation‐optimization. In the first part, based on the SDN architecture, we use the pull‐based random sampling method to quickly obtain the coarse‐grained measurement of flow traffic through OpenFlow protocol. In the second part, we insert some discrete values into the coarse‐grained measurement with the interpolation theory, then we optimize interpolation results until finding the optimal fine‐grained flow traffic measurement by utilizing the multiconstraint method. We verified the feasibility of the proposed measurement method, and simulation results show that the measurement error of the proposed method is under 25%, but the flow measurement overhead of this method occupies only 3.3% compared with that of the fine‐grained method.

Huo, Liuwei↗
DOE OSTI2023

On the Feasibility of Simulation-Driven Portfolio Scheduling for Cyberinfrastructure Runtime Systems

Runtime systems that automate the execution of applications on distributed cyberinfrastructures need to make scheduling decisions. Researchers have proposed many scheduling algorithms, but most of them are designed based on analytical models and assumptions that may not hold in practice. The literature is thus rife with algorithms that have been evaluated only within the scope of their underlying assumptions but whose practical effectiveness is unclear. It is thus difficult for developers to decide which algorithm to implement in their runtime systems.To obviate the above difficulty, we propose an approach by which the runtime system executes, throughout application execution, simulations of this very execution. Each simulation is for a different algorithm in a scheduling algorithm portfolio, and the best algorithm is selected based on simulation results. The main objective of this work is to evaluate the feasibility and potential merit of this portfolio scheduling approach, even in the presence of simulation inaccuracy, when compared to the traditional one-algorithm approach. We perform this evaluation via a case study in the context of scientific workflows. Our main finding is that portfolio scheduling can outperform the best one-algorithm approach even in the presence of relatively large simulation inaccuracies.

Casanova, Henri↗