Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “core 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,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

The Dark Energy Survey Supernova Programme: Modelling Selection Efficiency and Observed Core-collapse Supernova Contamination

The analysis of current and future cosmological surveys of Type Ia supernovae (SNe Ia) at high redshift depends on the accuratephotometric classification of the SN events detected. Generating realistic simulations of photometric SN surveys constitutes anessential step for training and testing photometric classification algorithms, and for correcting biases introduced by selectioneffects and contamination arising from core-collapse SNe in the photometric SN Ia samples. We use published SN time-seriesspectrophotometric templates, rates, luminosity functions, and empirical relationships between SNe and their host galaxies toconstruct a framework for simulating photometric SN surveys. We present this framework in the context of the Dark EnergySurvey (DES) 5-yr photometric SN sample, comparing our simulations of DES with the observed DES transient populations.We demonstrate excellent agreement in many distributions, including Hubble residuals, between our simulations and data.We estimate the core collapse fraction expected in the DES SN sample after selection requirements are applied and beforephotometric classification. After testing different modelling choices and astrophysical assumptions underlying our simulation,we find that the predicted contamination varies from 7.2 to 11.7 per cent, with an average of 8.8 per cent and an r.m.s. of 1.1 percent. Our simulations are the first to reproduce the observed photometric SN and host galaxy properties in high-redshift surveyswithout fine-tuning the input parameters. The simulation methods presented here will be a critical component of the cosmologyanalysis of the DES photometric SN Ia sample: correcting for biases arising from contamination, and evaluating the associatedsystematic uncertainty.

M Vincenzi↗
DOE OSTI2014

A spectral transform dynamical core option within the Community Atmosphere Model (CAM4)

A spectral transform dynamical core with an 85 spectral truncation resolution (T85) within the Community Atmosphere Model (CAM), version 4, is evaluated within the recently released Community Earth System Model, version 1.0 (CESM) global climate model. The spectral dynamical core option provides a well-known base within the climate model community from which to assess climate behavior and statistics, and its relative computational efficiency for smaller computing platforms allows it to be extended to perform climate length simulations using high-resolution configurations in the near term. To establish the characteristics of the CAM4 T85, an ensemble of simulations covering the present day observational period using forced sea surface temperatures and prescribed sea-ice extent are evaluated. Overall, the T85 ensemble attributes and biases are similar to a companion ensemble of simulations using the one degree finite volume (FV1) dynamical core, relative to observed and model derived datasets. Notable improvements with T85 compared to FV1 include the representation of wintertime Arctic sea level pressure and summer precipitation over the Western Indian subcontinent. The mean and spatial patterns of the land surface temperature trends over the AMIP period are generally well simulated with the T85 ensemble relative to observations, however the model is not able to capture the extent nor magnitude of changes in temperature extremes over the boreal summer, where the changes are most dramatic. Biases in the wintertime Arctic surface temperature and annual mean surface stress fields persist with T85 as with the CAM3 version of T85.

atmosphere model analysis↗
DOE OSTI2019

Phonon Thermal Conductivity and Phase Equilibria of Fractal Bi–Sb Nanoparticles

The thermal conductivity component associated with lattice vibrations is one of the quantities determining the thermoelectric activity of a material. We have simulated the dependences of phase composition and the phonon component of the thermal conductivity associated with it on the shape of nanoparticles of a Bi–Sb alloy with an equiatomic composition and with core–shell configuration. The shape of a particle is simulated by a coefficient corresponding to the extent of deviation of the particle shape from spherical or by its fractal dimension. It is shown that mutual solubilities of components depend on the nanoparticle shape and on the mutual arrangement of coexisting phases, and the thermodynamic equilibrium position for particles with complex morphology corresponds to the homogeneous state. Homogenization of a nanoparticle reduces the phonon component of its thermal conductivity by 70–80%.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

Direct Numerical Simulation of Autoignition in a Jet in a Cross-Flow Configuration: ALCF-2 Early Science Program (Technical Report)

Autoignition in turbulent flows is a challenging fundamental problem due to the intricate coupling of different physical and chemical processes extending over multiple flow and chemistry scales. At the same time, the improved understanding and ability to predict autoignition in flows characterized by considerable fluctuations of velocity, composition, and temperature is essential for the development of novel low-emission concepts for power generation. The aim of this project is to study the fundamental aspects of autoignition in a fuel-air mixing device directly applicable to mixing ducts in gas turbines. The NEK5000-based code for low Mach number reactive flows is used to perform very large scale direct numerical simulations of autoignition of a diluted hydrogen jet ejected in a cross-flowing stream of hot turbulent air in a laboratory-scale configuration. We report on our experience running NEK5000 on the new BGQ system at ALCF mira during the early science period (ESP). First of all, the most efficient problem size per MPI-rank is obtained through core-level efficiency metric measured from the target simulation. Furthermore, the most efficient number of ranks is found through strong scaling experiments. Finally, low-level insight into the observed parallel efficiency is enabled through IBM’s HPC Toolkit libraries.

Abdilghanie, A.↗
DOE OSTI2013

The Mean Climate of the Community Atmosphere Model (CAM4) in Forced SST and Fully Coupled Experiments

The Community Atmosphere Model, version 4 (CAM4), was released as part of the Community Climate System Model, version 4 (CCSM4). The finite volume (FV) dynamical core is now the default because of its superior transport and conservation properties. Deep convection parameterization changes include a dilute plume calculation of convective available potential energy (CAPE) and the introduction of convective momentum transport (CMT). An additional cloud fraction calculation is now performed following macrophysical state updates to provide improved thermodynamic consistency. A freeze-drying modification is further made to the cloud fraction calculation in very dry environments (e.g., the Arctic), where cloud fraction and cloud water values were often inconsistent in CAM3. In CAM4 the FV dynamical core further degrades the excessive trade-wind simulation, but reduces zonal stress errors at higher latitudes. Plume dilution alleviates much of the midtropospheric tropical dry biases and reduces the persistent monsoon precipitation biases over the Arabian Peninsula and the southern Indian Ocean. CMT reduces much of the excessive trade-wind biases in eastern ocean basins. CAM4 shows a global reduction in cloud fraction compared to CAM3, primarily as a result of the freeze-drying and improved cloud fraction equilibrium modifications. Regional climate feature improvements include the propagation of stationary waves from the Pacific into midlatitudes and the seasonal frequency of Northern Hemisphere blocking events. A 1° versus 2° horizontal resolution of the FV dynamical core exhibits superior improvements in regional climate features of precipitation and surface stress. Finally, improvements in the fully coupled mean climate between CAM3 and CAM4 are also more substantial than in forced sea surface temperature (SST) simulations.

cloud fraction calculation↗
DOE OSTI2025

The ePIC Simulation Campaign Workflow on the Open Science Grid

The ePIC collaboration is realizing the first experiment of the future Electron-Ion Collider (EIC) at the Brookhaven National Laboratory that will allow for a precision study of the nucleons and the nucleus at the scale of sea quarks and gluons through the study of electron-proton/ion collisions. This paper will discuss the current workflow for running centralized simulation campaigns for ePIC on the Open Science Grid (OSG) infrastructure. This involves monthly releases of ePIC software and container deployments to CVMFS, generation of input datasets in HepMC format according to collaboration-defined policy, using Snakemake in CI/CD for validation and benchmarking, and submitting jobs to the OSG condor scheduler for opportunistic running on available resources. File transfers utilize XrootD, and Rucio is used for data management. The workflow is continuously refined to improve daily throughput (currently 50-100k core hours per day) and minimize job failures. Since May 2023, monthly simulation campaigns employing the workflow have cumulatively used over 20 million core hours on the OSG and produced over 350 TB of simulation data. The campaigns incorporate simulations for the broad science program of the EIC and are actively used for the detector and physics studies in preparation of the Technical Design Report (TDR).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2011

Computational Fluid Dynamics Demonstration of Rigid Bodies in Motion

The Design Analysis Branch (NE-Ml) at the Kennedy Space Center has not had the ability to accurately couple Rigid Body Dynamics (RBD) and Computational Fluid Dynamics (CFD). OVERFLOW-D is a flow solver that has been developed by NASA to have the capability to analyze and simulate dynamic motions with up to six Degrees of Freedom (6-DOF). Two simulations were prepared over the course of the internship to demonstrate 6DOF motion of rigid bodies under aerodynamic loading. The geometries in the simulations were based on a conceptual Space Launch System (SLS). The first simulation that was prepared and computed was the motion of a Solid Rocket Booster (SRB) as it separates from its core stage. To reduce computational time during the development of the simulation, only half of the physical domain with respect to the symmetry plane was simulated. Then a full solution was prepared and computed. The second simulation was a model of the SLS as it departs from a launch pad under a 20 knot crosswind. This simulation was reduced to Two Dimensions (2D) to reduce both preparation and computation time. By allowing 2-DOF for translations and 1-DOF for rotation, the simulation predicted unrealistic rotation. The simulation was then constrained to only allow translations.

Camarena, Ernesto↗
DOE OSTI2015

Transmission electron microscopy and ab initio calculations to relate interfacial intermixing and the magnetism of core/shell nanoparticles

Significant efforts towards understanding bi-magnetic core-shell nanoparticles are underway currently as they provide a pathway towards properties unavailable with single-phased systems. Recently, we have demonstrated that the magnetism of γ-Fe2O3/CoO core-shell nanoparticles, in particular, at high temperatures, originates essentially from an interfacial doped iron-oxide layer that is formed by the migration of Co{sup 2+} from the CoO shell into the surface layers of the γ-Fe2O3 core [Skoropata et al., Phys. Rev. B 89, 024410 (2014)]. To examine directly the nature of the intermixed layer, we have used high-resolution transmission electron microscopy (HRTEM) and first-principles calculations to examine the impact of the core-shell intermixing at the atomic level. By analyzing the HRTEM images and energy dispersive spectra, the level and nature of intermixing was confirmed, mainly as doping of Co into the octahedral site vacancies of γ-Fe2O3. The average Co doping depths for different processing temperatures (150 °C and 235 °C) were 0.56 nm and 0.78 nm (determined to within 5% through simulation), respectively, establishing that the amount of core-shell intermixing can be altered purposefully with an appropriate change in synthesis conditions. Through first-principles calculations, we find that the intermixing phase of γ-Fe2O3 with Co doping is ferromagnetic, with even higher magnetization as compared to that of pure γ-Fe2O3. In addition, we show that Co doping into different octahedral sites can cause different magnetizations. This was reflected in a change in overall nanoparticle magnetization, where we observed a 25% reduction in magnetization for the 235 °C versus the 150 °C sample, despite a thicker intermixed layer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Control And Optimization Modular Modeling Application For Nuclear Deployment

The purpose of the COMMAND code is to provide a flexible, scalable tool for use in developing, integrating, and testing the technologies necessary for achieving autonomous operations of advanced nuclear reactors. The code enables users to efficiently implement custom simulations and experiments by combining key methods from different software modules. These modules are focused on: modeling and simulation tools, such as nuclear simulation tools used for high-fidelity modeling (e.g., Reactor Excursion and Leak Analysis Program [RELAP5-3D] and Monte Carlo N-Particle [MCNP]); machine learning and optimization tools (e.g., anomaly detection and data-driven modeling techniques); advanced control in its digital, high-performance, and supervisory control forms (e.g., proportional integral derivative (PID) control and model predictive control (MPC); and integration with hardware through industrial communication protocols. To ensure flexibility and scalability, COMMAND was designed to be both modular—the software “pieces” all inherit from generic building blocks and can be combined and connected to create complicated simulations—and high performing—designed for parallel processing, enabling simulations and experiments to take advantage of multi-core computers, servers, and nodes. The code is written in the Python programming language due to the language's popularity, active community, and open-source and cross-platform nature. Maintaining consistency with other simulation tools used within the nuclear energy community, users implement simulations and experiments through text input files, which define components, parameters, connections, etc., through lines of text. Given that COMMAND is written in Python, these input files are native Python scripts, and so use the standard Python structure and formatting. This also enables users to take advantage of Python's extensive package library to develop custom capabilities for their specific use cases.

Faber, Jacob [Idaho National Laboratory (INL), Ida↗
DOE OSTI2023

Microstability of β ~ 1 tokamak equilibria

High-power-density tokamaks offer a potential solution to design cost-effective fusion devices. One way to achieve high power density is to operate at a high β value (the ratio of thermal to magnetic pressure), i.e. β~1 . However, a β~1 state may be unstable to various pressure- and current-driven instabilities or have unfavourable microstability properties. To explore these possibilities, we generate β~1 equilibria and investigate their stability. First, we demonstrate the generation of high-β equilibria with the computer code VMEC. We then analyse these equilibria to determine their stability against the infinite-n ideal-ballooning mode. We follow that by engaging in a detailed microstability study using the GS2 code, beginning with assessments of electrostatic ion-temperature-gradient and trapped election mode instabilities. We observe interesting behaviour for the high-β equilibria – stabilization of these modes through two distinct mechanisms – large negative local shear and reversal of electron precession drift. Finally, we perform electromagnetic gyrokinetic simulations and observe enhanced stability in the outer core of high-β equilibria and absence of kinetic ballooning modes in the negative-triangularity, high-β equilibria. The enhanced outer-core stability of high-β equilibria is different from their lower-β counterparts and offers an alternative, potentially favourable regime of tokamak operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Core-collapse contamination in photometric samples of Type Ia Supernovae

This is an exciting time for cosmology with type Ia supernovae (SNe Ia). The recentlyconcluded Dark Energy Survey SN programme (DES-SN) has obtained the largest anddeepest high-redshift cosmological SN Ia sample, and the Vera Rubin Observatory isexpected to observe at least one order of magnitude more SNe Ia in the next decade.In both these experiments, only a limited fraction (.10 per cent) of the SNe can bespectroscopically classified. This leaves us with large ‘photometric’ SN samples, withthe potential for significant contamination by core-collapse SNe that may bias SN Iacosmological measurements. This thesis demonstrates how this contamination can bemodelled and accounted for in current and future cosmological analyses.First, I present state-of-the-art simulations of the SN universe. These are designedto accurately model the population of SNe Ia, peculiar SNe Ia and core-collapse SNe, aswell as their host galaxies. To improve the diversity and quality of the simulated corecollapseSNe, I build a new library of core-collapse SN templates using spectroscopicand photometric (optical and near-ultraviolet) data of 67 core-collapse SNe from theliterature. I account for our incomplete knowledge of core-collapse SN properties bygenerating a set of SN simulations (rather than a single one), each exploring differentmodelling choices and template libraries. I then characterise selection effects in theDES-SN survey and incorporate them in the simulations, thus obtaining a series ofDES-like simulated SN samples that can be compared to the observed DES-SN data.The agreement between the simulations and data is excellent across many observed SNproperties, including Hubble residuals. These simulations are the first to reproduce theobserved photometric SN and host galaxy properties in high-redshift surveys with no fine-tuning of the input parameters.I use my simulation framework to train and test the performance of SuperNNova,a photometric SN classifier based on recurrent neural networks. I explore differenttraining and validation strategies and show that, across all the DES-SN simulationstested, SuperNNova reduces core-collapse SN contamination to 0.8–3.5 per cent. Ithen show that biases due to contamination on the equation-of-state of dark energy,w, are < 0:008 when using our reference SuperNNova model. This compares to anexpected statistical uncertainty on w from the DES-SN sample of 0:039, thus showing that contamination is not a limiting systematic for the cosmological analysis of theDES-SN sample.The results presented in this thesis are the foundation of the DES SN Ia cosmologicalanalysis; they also provide important implications for the future of SN cosmology,as they demonstrate that contamination is not expected to significantly degrade thecosmological figure of merit of the Rubin SN Ia analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Modernization and optimization of a legacy open-source CFD code for high-performance computing architectures

Legacy codes remain a crucial element of today's simulation-based engineering ecosystem due to the extensive validation process and investment in such software. The rapid evolution of high-performance computing architectures necessitates the modernization of these codes. One approach to modernization is a complete overhaul of the code. However, this could require extensive investments, such as rewriting in modern languages, new data constructs, etc., which will necessitate systematic verification and validation to re-establish the credibility of the computational models. The current study advocates using a more incremental approach and is a culmination of several modernization efforts of the legacy code MFIX, which is an open-source computational fluid dynamics code that has evolved over several decades, widely used in multiphase flows and still being developed by the National Energy Technology Laboratory. Two different modernization approaches,‘bottom-up’ and ‘top-down’, are illustrated. Here, preliminary results show up to 8.5x improvement at the selected kernel level with the first approach, and up to 50% improvement in total simulated time with the latter were achieved for the demonstration cases and target HPC systems employed.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2016

Space Launch System Booster Separation Aerodynamic Database Development and Uncertainty Quantification

The development of the aerodynamic database for the Space Launch System (SLS) booster separation environment has presented many challenges because of the complex physics of the ow around three independent bodies due to proximity e ects and jet inter- actions from the booster separation motors and the core stage engines. This aerodynamic environment is dicult to simulate in a wind tunnel experiment and also dicult to simu- late with computational uid dynamics. The database is further complicated by the high dimensionality of the independent variable space, which includes the orientation of the core stage, the relative positions and orientations of the solid rocket boosters, and the thrust lev- els of the various engines. Moreover, the clearance between the core stage and the boosters during the separation event is sensitive to the aerodynamic uncertainties of the database. This paper will present the development process for Version 3 of the SLS booster separa- tion aerodynamic database and the statistics-based uncertainty quanti cation process for the database.

Chan, David T.↗
DOE OSTI2014

Influence of shape, size and internal structure on magnetic properties of core-edge nanodots with perpendicular anisotropy

The properties of perpendicularly magnetized isolated nanodots different in shape, size, and internal structure are simulated by micromagnetic calculations. Investigated dots are magnetically uniform, or they are composed of a core and an edge characterized by different anisotropy—stronger or weaker than that of the core. Based on calculated hysteresis loops, we discuss in details the magnetization reversal processes, stability of magnetic structures, and spin configurations in the dots.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Advancing Neutrino Simulation Modeling with MARLEY: Insights from the NNSA-MSIIP Internship

Core-collapse supernovae are intense sources of tens-of-MeV neutrinos. However, there is no experimental data to validate the current event generator MARLEY (Model of Argon Reaction Low Energy Yields) that can model supernova neutrinos. To validate the model, I developed a new muon capture feature within the MARLEY simulation framework by coding key functions in C++, Python, and ROOT. I generated and analyzed one million simulated muon capture events and comparing the results to experimental data. To improve the model, I worked on optimizing the model’s parameters to improve its precision using statistical methods.

Wong, Baker [Fermilab]↗
DOE OSTI2026

Interfacial Inversion of Stealth Surfactants

Amphiphilic macromolecular surfactants segregate to liquid–liquid interfaces, thereby reducing the interfacial tension and free energy. Here, we investigated “stealth surfactants” in the form of core–shell bottlebrush polymers comprised of pH-responsive diblock copolymer side chains forming a hydrophilic core and a hydrophobic shell, enabling solubility in oil. At liquid–liquid interfaces, these polymers undergo a structural “inversion”, with hydrophilic blocks segregating into the aqueous phase and hydrophobic blocks residing in the oil phase. The reconfiguration kinetics and surfactant properties are influenced by multiple factors, including the molecular weights of the backbone and side chain components, the hydrophilic-to-hydrophobic balance of the side chains, and the pH of the aqueous phase. An observed nonmonotonic dependence of interfacial tension with time is attributed to a progressive structural inversion, where the projected area of the macromolecule onto the interface decreases. To validate this inversion hypothesis, interfacial properties were characterized by sum-frequency generation vibrational spectroscopy, which revealed configurational changes of the core–shell bottlebrush polymers at the fluid interface and revealed a pH-dependent interfacial coverage. Coarse-grained molecular dynamics simulations supported these experimental findings, showing that the pH-responsive core and hydrophobic shell assume a time-averaged configuration with orientations parallel and perpendicular to the plane of the interface, respectively. These findings open routes to design multistimuli-responsive polymeric surfactants and compatibilizers, expanding their potential applications in advanced interfacial systems.

Stealth surfactants↗
DOE OSTI2014

High pressure induced phase transition and superdiffusion in anomalous fluid confined in flexible nanopores

The behavior of a confined spherical symmetric anomalous fluid under high external pressure was studied with Molecular Dynamics simulations. The fluid is modeled by a core-softened potential with two characteristic length scales, which in bulk reproduces the dynamical, thermodynamical, and structural anomalous behavior observed for water and other anomalous fluids. Our findings show that this system has a superdiffusion regime for sufficient high pressure and low temperature. As well, our results indicate that this superdiffusive regime is strongly related with the fluid structural properties and the superdiffusion to diffusion transition is a first order phase transition. We show how the simulation time and statistics are important to obtain the correct dynamical behavior of the confined fluid. Our results are discussed on the basis of the two length scales.

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

Characterization of swift heavy ion irradiation damage in ceria

We have examined microstructural evolution in irradiated ceria (CeO2) using swift heavy ion irradiation, electron microscopy, and atomistic simulation. CeO2, a UO2 fuel surrogate, was irradiated with gold ions at an energy of 1 GeV to fluences up to 1x1014 ions/cm2. Transmission electron microscopy accompanied by electron energy loss spectroscopy showed that the ion tracks were of similar size at all fluences, and that there was no chemical change in the ion track core. Classical molecular dynamics simulations of thermal spikes in CeO2 with energy deposition of 12 and 36 keV/nm show damage consisting of isolated point defects at the lower energy and defect clusters at 36 keV/nm, with no amorphization at either energy. Inferences are drawn from modeling about density changes in the ion track and the formation of interstitial loops that shed light on features observed by electron microscopy of swift heavy ion irradiated ceria.

Yablinsky, Clarissa↗