Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Advanced Codes”

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

Development and Validation of Physics-Based Sub-Models of High-Pressure Supercritical Fuel Injection at Diesel Conditions Insert

The project sought to develop and validate accurate, physics-based, numerical submodels and implement them in computational fluid dynamics (CFD) software codes to enable accurate predictions of diesel sprays at pressures higher than the current systems. The combined experimental-computational effort has four main objectives: (1) acquire spatially and temporally resolved scalar measurements of fuel-oxidizer mixing in the near field of the fuel jet, (2) develop software modules to calculate thermo-physical properties of diesel surrogates and their mixtures with oxidizer of interest in real-gas operating regime, (3) integrate the real-fluid model into a commercial and an open-source CFD code to simulate fuel-air mixing in high-pressure diesel sprays, and (4) assess the robustness, accuracy, and uncertainty of the integrated CFD solvers in computing diesel sprays. Project accomplishments are consistent with the original goals of the project. Although we acquired spatially and temporally resolved scalar measurements, they were limited to refractive index difference and not to thermodynamic properties such as temperature and equivalence ratio. However, we successfully developed the theoretical background needed to acquire such measurements at diesel conditions and demonstrated the capability of the RSD for scalar measurements in a canonical turbulent jet flow replicating features of the diesel jet. Additional work will be needed to implement these procedures at high-pressure, high-temperature engine-like conditions of the CPFR. Mainly, we did not anticipate the need to find the liquid boundary and its statistical features or the state relationships between refractive index and thermodynamic properties, assuming the widely known Dale-Gladstone equation would suffice. Still, efforts devoted to these tasks were essential to fulfil the goals of the project and serves as the foundation on which many of the conclusions and key contributions of this project are based. Further assessment and application of the developed CFD codes is recommended to fully realize the benefits of this research.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Modeling the Effects of Artificial Drainage on Agriculture-dominated Watersheds using a Fully Distributed Integrated Hydrology Model: Datasets, scripts, model files

This model-data archive supports the research paper that demonstrates the integration of agricultural drainage features—specifically, narrow engineered ditches and tile drains—into a fully distributed, basin-scale integrated surface-subsurface hydrology model (ISSHM), Amanzi-ATS. The model employs innovative computational meshes aligned with agricultural ditches and incorporates the physically based Hooghoudt's drainage equation to simulate tile drainage, offering a novel strategy that enhances the accuracy of hydrological simulations.The archived dataset includes input parameters, model configurations, and select simulation outputs for the Amanzi-ATS model that successfully captured the streamflow patterns in the Portage River Watershed as validated by USGS gauge readings. Jupyter notebook for the preparation of model inputs and post-processing of outputs are also included. The model's predictive performance achieved a normalized Kling-Gupta Efficiency (KGE) of 0.81, surpassing SWAT without the necessity for site-specific calibration.The Amanzi-ATS model presented in this modeL-data archive allows for numerical experiments to explore the shifts in the flow structure under different drainage scenarios. As a tool for advancing the understanding of distributed hydrological responses and nutrient cycling, this archived model provides valuable insights for researchers, modelers, and decision-makers involved in watershed management and environmental modeling.The Watershed Workflow package is implemented in Python3. The Jupyter notebooks can be executed through multiple open-source tools, for example, Anaconda Jupyter Lab, VS Studio Code, etc. Other data files include CSV and HDF5 files, which can be read through Python scripts. The input files for the ATS model, open-source integrated hydrology, and transport model, are in XML format and can be edited in any commonly used text editors.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Coupling Surface Flow with High-performance Subsurface Reactive Flow and Transport Code PFLOTRAN

Water exchange between the surface and subsurface is important for both water resource management and environmental protection. In this paper, we develop coupled surface and subsurface flow simulation capability in a parallel subsurface flow and reactive transport code PFLOTRAN. We sequentially couple the diffusion wave-based surface flow with the subsurface flow governedby the Richards equation in PFLOTRAN. These two flow domains are linked with a boundary condition switching method that ensures continuity of pressure and flux at the surface-subsurface interface. We verify the coupled code against other existing hydrologic models and observation data using a number of numerical experiments. The coupled hydrological model exhibits good performance in strong parallel scaling tests. The new coupled surface and subsurface simulator significantly advance community simulation capability towards improving integrated hydrologic and biogeochemical understanding of complex systems such as watersheds and river corridors. Keywords: Surface flow, Integrated hydrological modeling, Boundary condition switching, Parallel computing

Wu, Runjian↗
DOE OSTI2020

Radiation shielding analysis of the concrete hutch of a long beamline of the Advanced Photon Source Upgrade project

The radiation shielding analysis of the external concrete enclosure of a feature long beamline of the Advanced Photon Source Upgrade project is described in this work. The beamline has 2 undulators that can operate in two different configurations and the source terms are estimated for the worst configuration using the SRW code in the OASYS package as well as the STAC8 code. The shielding analysis is carried out using the STAC8 code, and for the case of pink beam, with the FLUKA code also. Effective dose rates on contact with the outside surface of the concrete walls for pink, reflected full spectrum and monochromatic beams are calculated. For the monochromatic beams, 2 sets of discrete bandwidths (BW) from the Double Crystal Monochromator (DCM) and the full reflectance of the Double Multilayer Monochromator (DMM) are used. Further, due to the size of the direct beam and the energy domain of interest in the direct and scattered fields, the importance of the equivalent dose to the lens of the eye is discussed with respect to the shielding calculations of such facilities. The dose rates outside the concrete hutch in the lateral directions calculated using STAC8 is about 50% higher than the FLUKA results for the case of pink beam but may be considered as a reasonable agreement due to the differences in the calculational methodologies. Those dose rates in the lateral and forward directions are higher for the pink beam to allow full occupancy from the radiation protection point of view. For the full spectrum after 3 reflections, the effective dose rates are less than 0.5 Sv/h in the lateral direction but is higher in the forward region below about 22° with respect to the beam direction. For the DCM beams, except for the scenario with one mirror and 0.1%BW, all dose rates are less than 0.5 Sv/h outside the concrete hutch in the lateral directions. Here, the higher harmonics are what contributes to the dose outside the shield and the bandwidths assumed to derive the photon flux thus becomes important. The dose rates with the DMM are generally lower compared to the 0.1% BW DCM beam but are higher than the results obtained with XOP calculated bandwidths.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Data on erosion and hydrogen fuel retention in Beryllium plasma-facing materials

ITER will use beryllium as a plasma-facing material in the main chamber, covering a total surface area of about 620 m 2 . Given the importance of beryllium erosion and co-deposition for tritium retention in ITER, significant efforts have been made to understand the behaviour of beryllium under fusion-relevant conditions with high particle and heat loads. This paper provides a comprehensive report on the state of knowledge of beryllium behaviour under fusion-relevant conditions: the erosion mechanisms and their consequences, beryllium migration in JET, fuel retention and dust generation. The paper reviews basic laboratory studies, advanced computer simulations and experience from laboratory plasma experiments in linear simulators of plasma–wall interactions and in controlled fusion devices using beryllium plasma-facing components. A critical assessment of analytical methods and simulation codes used in beryllium studies is given. The overall objective is to review the existing set of data with a broad literature survey and to identify gaps and research needs to broaden the database for ITER.

Erosion–deposition↗
DOE OSTI2023

On relaxations of the max k -cut problem formulations

Here, a tight continuous relaxation is a crucial factor in solving mixed integer formulations of many NP-hard combinatorial optimization problems. The (weighted) max k-cut problem is a fundamental combinatorial optimization problem with multiple notorious mixed integer optimization formulations. In this paper, we explore four existing mixed integer optimization formulations of the max k-cut problem. Specifically, we show that the continuous relaxation of a binary quadratic optimization formulation of the problem is: (i) stronger than the continuous relaxation of two mixed integer linear optimization formulations and (ii) at least as strong as the continuous relaxation of a mixed integer semidefinite optimization formulation. We also conduct a set of experiments on multiple sets of instances of the max k-cut problem using state-of-the-art solvers that empirically confirm the theoretical results in item (i). Furthermore, these numerical results illustrate the advances in the efficiency of global non-convex quadratic optimization solvers and more general mixed integer nonlinear optimization solvers. As a result, these solvers provide a promising option to solve combinatorial optimization problems. Our codes and data are available on GitHub.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Real-space solution to the electronic structure problem for nearly a million electrons

We report a Kohn–Sham density functional theory calculation of a system with more than 200 000 atoms and 800 000 electrons using a real-space high-order finite-difference method to investigate the electronic structure of large spherical silicon nanoclusters. Our system of choice was a 20 nm large spherical nanocluster with 202 617 silicon atoms and 13 836 hydrogen atoms used to passivate the dangling surface bonds. To speed up the convergence of the eigenspace, we utilized Chebyshev-filtered subspace iteration, and for sparse matrix–vector multiplications, we used blockwise Hilbert space-filling curves, implemented in the PARSEC code. For this calculation, we also replaced our orthonormalization + Rayleigh–Ritz step with a generalized eigenvalue problem step. We utilized all of the 8192 nodes (458 752 processors) on the Frontera machine at the Texas Advanced Computing Center. We achieved two Chebyshev-filtered subspace iterations, yielding a good approximation of the electronic density of states. Our work pushes the limits on the capabilities of the current electronic structure solvers to nearly 106 electrons and demonstrates the potential of the real-space approach to efficiently parallelize large calculations on modern high-performance computing platforms.

Chemistry↗
DOE OSTI2023

Volume 9: Critical Facility with Add-On Ion Beam

Completion of the originally planned critical facility has been considered as part of the High Flux Isotope Reactor (HFIR) Sustaining and Enhancing Neutron Science Initiative at Oak Ridge National Laboratory (ORNL). A working group of ORNL staff members was formed to develop the idea and to recommend one or more configurations to best support future HFIR operations and scientific capacities. HFIR was designed with a critical pool in the reactor bay that was projected to be outfitted as a critical facility. The primary purpose of the planned critical facility was to measure the subcritical worth of fresh fuel elements to support startup requirements. However, the Y-12 National Security Complex already housed a critical facility that was used for this purpose. This report presents the working group’s efforts, including determination of the proposed critical facility’s high-impact benefits. A generic low-power critical facility would be employed for reactor physics measurements, code and data validation, reactor operator and staff training, and education. This facility would be instrumental in supporting current HFIR operations, conversion of HFIR to low-enriched uranium, existing light water reactor operations, and advanced reactor development and deployment.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Modeling Heat Pipes with Non-Condensable Gases

The report showcases advancements in modeling heat pipes, considering both scenarios with and without non-condensable gases. To achieve this, two distinct modeling approaches are juxtaposed. Firstly, the effective conduction model, implemented in the MOOSE-based code Sockeye, is employed. Secondly, a first-of-a-kind two-phase Euler-Euler Computational Fluid Dynamics (CFD) model is developed using the STAR-CCM+ code. Both models undergo validation against experimental data, acknowledging the inherent uncertainties associated with each modeling assumption. Interestingly, the non-tuned CFD model surpasses the performance of the calibrated conduction model for heat pipes operating under both conditions: with and without non-condensable gases. It's worth noting, however, that the CFD models entail significantly longer runtimes compared to the conduction models. Nevertheless, the insights garnered from the CFD model shed invaluable light on the intricate operational dynamics of heat pipes. Future work involves broadening the validation scope of these models and continuing their development to enhance their utility as robust tools for heat pipe design and operational support.

42 ENGINEERING↗
DOE OSTI2022

High-Temperature Gas-Cooled Reactor Research Survey and Overview: Preliminary Data Platform Construction for the Nuclear Energy University Program

Since the U.S. Department of Energy Office of Nuclear Energy initiated the Nuclear Energy University Program (NEUP) in 2009, there are 29 NEUP projects focusing on high-temperature gas-cooled reactor (HTGR) research up to July 2022. The resultant research product, either experimental or computational, were published as final NEUP reports, journal articles and conference proceedings. However, these federally funded products have been scattered and sometimes cannot be easily accessed. To improve access to this valuable HTGR validation data and optimize the return on the significant investment made by the Department of Energy, the Advanced Reactor Technologies (ART) Gas-Cooled Reactor (GCR) program started a survey of completed and ongoing HTGR NEUP projects to develop a public-access database specific for HTGRs applications that can be used to retrieve computational fluid dynamics and system code validation data. This effort will help guide future NEUP-funded research, define new state of the ART Phenomena Identification and Ranking Table (PIRT), and promote the usage of this data in the codes validation matrices. This report provides an overview of the NEUP-funded HTGR-related research projects from Fiscal Year (FY) 2009–2021 and identifies validation knowledge gaps still existing in HTGR thermal-fluid research. A preliminary data platform has been developed for the 29 NEUP projects investigating HTGR thermal hydraulics, including their final reports as well as the available scientific publications. As an ultimate goal for this work, the ART-GCR program will create a central database at Idaho National Laboratory to identify, organize, and store these datasets generated by experimental investigations or computational models, experimental facility descriptions, and publicly-available academic products from the HTGR-related NEUP projects and provide future guidance for the storage and transmission of important project documentations for later NEUP projects as well.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Modeling of deuterium and carbon radiation transport in MAST-U tokamak advanced divertors

Modest effects of deuterium and carbon radiation opacity in the Super-X and snowflake divertor plasmas are predicted for MAST Upgrade tokamak with core plasma input power 2.5–5 MW and plasma current 1 MA. The radiation transport modeling is based on the SOLPS-EIRENE and UEDGE code divertor plasma predictions. Two radiation transport models are used: one is based on a full radiation transport equation implemented in the radiation transport and collisional-radiative code CRETIN (without feedback on the background plasma), and another is an internal self-consistent UEDGE model with ionization, recombination, and heating rates corrected for Ly α line trapping based on the escape probability model implemented in CRETIN. In MAST-U, the Super-X and snowflake divertor plasmas are predicted to reach detached regimes at lower upstream densities than the standard divertor, and the conclusion still holds with radiation transport effects included. At neutral densities m -3 , modest Ly α deuterium line trapping with optical depths 10–15 is predicted in the Super-X divertor. Divertor plasmas are optically thin to other Lyman and Balmer lines, as well as to strong C III and C IV lines that are responsible for most of divertor radiated power. Insignificant changes (within a few percent) to divertor deuterium ionization and recombination rates are found. Radiation fluxes on outer divertor target are modified within a factor of 2–3 when the radiation transport is accounted for, and a similar variation is found due to the line shape models that define the absorption and emission line profiles in the radiation transport modeling. The predicted Lyman and Balmer spectral intensities are significantly modified due to radiation trapping. A measurement of divertor radiation transport effects is discussed using the Ly β /Ba α line ratio. In the snowflake divertor configuration, divertor plasmas are found to be optically thin to Lyman series lines within a large range of parameter variations that include upstream density, divertor transport coefficients, and magnetic configurations. Finally, modest radiation transport effects are only found in a few cases with strongest divertor transport and magnetic configurations closest to the ideal snowflake configuration, however, the plasma background models that were used are yet to be validated with an experiment.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Smart Ventilation for Advanced California Homes

This project investigated smart ventilation approaches to minimize energy use for providing indoor air quality (IAQ) in high performance new California homes. Evaluation criteria included annual ventilation-related energy, peak energy and time-of-use savings, and the indoor air quality relative to a minimally code-compliant ventilation system. The simulations used CONTAM’s air flow and contaminant transport model, combined with the EnergyPlus building loads model. House types representing the default California Energy Code compliance homes were investigated for four California climate zones, covering a wide range of climate types. Both single and multi-zone smart ventilation controls were investigated. Contaminant sources included contaminants emitted continuously and varying with time, temperature and relative humidity, episodic emissions from occupant activities and outdoor particles. Single-zone ventilation controls that varied ventilation depending on outdoor temperatures were able to consistently save half of ventilation-related energy without compromising long-term IAQ. Ventilation strategies that tracked occupancy were less successful, because this work included generic contaminants with constant background emission rates. Energy performance for occupancy controls improved with a one-hour pre-occupancy flush out strategy. The addition of zoning ventilation controls did not offer significant IAQ to energy improvements compared to non-zonal versions of the same ventilation system type. The best controls had HVAC energy savings of 10-20%, with individual cases reaching up to 40% savings. However, these savings cannot be achieved without worsening personal exposures for at least one contaminant. A metric is needed to assess the competing changes in exposure to different contaminants in order to determine the net-health impacts of a control strategy. Controls that directly sensed contaminants and controlled them to acceptable levels showed that the California OEHHA limit for formaldehyde completely dominates system performance, with homes not able to meet the limit even with continuous operation of a fan sized to twice the current code minimum.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Dynamical logical qubits in the Bacon-Shor code

The Bacon-Shor code is a quantum error correcting subsystem code composed of weight-2 check operators that admits a single logical qubit, and has distance 𝑑 on a 𝑑×𝑑 square lattice. We show that when viewed as a Floquet code, by choosing an appropriate measurement schedule of the check operators, it can additionally host several dynamical logical qubits. Specifically, we identify a period-4 measurement schedule of the check operators that preserves logical information between the instantaneous stabilizer groups. Such a schedule not only measures the usual stabilizers of the Bacon-Shor code, but also measures and promotes gauge operators of the parent subsystem code to additional temporary stabilizers that protect the dynamical logical qubits against errors. We show that the code distance of these Floquet-Bacon-Shor codes scales as Θ⁢(𝑑/√𝑘) on an 𝑛=𝑑×𝑑 lattice with 𝑘 dynamical logical qubits, along with the logical qubit of the parent subsystem code. Unlike the usual Bacon-Shor code, the Floquet-Bacon-Shor code family introduced here can therefore saturate the subsystem bound 𝑘⁢𝑑=𝑂⁡(𝑛). Moreover, several errors are shown to be self-corrected purely by the measurement schedule itself. This work provides insights into the design space for dynamical codes and expands the known approaches for constructing Floquet codes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

CVEVOLVE

CVEvolve is an agentic AI system for autonomous algorithm discovery for scientific data processing. It creates workflows where large language model agents freely set up and configure development environments and evaluation harnesses, develop and improve data processing algorithms with designed exploration-exploitation balancing mechanisms, log history and findings in a structured database, and run holdout testing to ensure algorithm generalizability. CVEvolve offers a zero-code interface and does not require users to provide structured data and evaluation scripts.

Cherukara, MatthewJoseph [Argonne National Laborat↗
DOE OSTI2021

Summary of 2020 edge-plasma modeling for FNSFand plans for 2021

This year, following the study of liquid-wall options, our edge-plasma modeling has returned to analyzing issues associated with a solid tungsten-wall, Fusion Nuclear Science Facility (FNSF)-size device, primarily using the 2D UEDGE plasma/neutral transport code. The focus is to predict plasma/neutral heat and particle fluxes to various wall components and determine the concentration of injected and wall-sputtered impurities that reach the core plasma region. Our models have been advanced in a number of areas compared to the initial modeling done in the 2017 timeframe: this year we include each individual charge-state of impurity ions instead of assuming a fixed concentration of impurities relative to the deuterium-tritium (DT) fuel plasma, including the interaction between inner and outer divertor legs. Here edge fueling of the impurity and pumping in the private-flux region is incorporated. Increasingly detailed atomic-physics rate files for impurity radiation have been compared, as well as the role of molecules. The impact of spatial variation of the turbulence-driven radial plasma transport is considered to contribute to the understanding of the heat-flux-limit effects of choosing a double-null or single-null divertor configuration. Detailed plasma/neutral wall fluxes, densities, and temperatures have been provided to the ERO group for their more detailed studies of wall sputtering and near-surface impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Electric Vehicle Supply Equipment and Considerations for a Reasonable Rate of Return

This white paper identifies the following recommendations for states as they consider issues related to electric vehicle supply equipment program income under the National Electric Vehicle Infrastructure program and other programs covered by National Electric Vehicle Infrastructure Standards and Requirements in the U.S. Code of Federal Regulations (23 CFR 680): 1. Conduct value analysis for better outcomes 2. Promote pricing transparency

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Multigroup cross section generation capability in GRIFFIN

GRIFFIN is an advanced reactor multiphysics application built on the object-oriented simulation environment (MOOSE) and is jointly developed by Idaho National Laboratory and Argonne National Laboratory. The cross section application programming interface, originally developed for the PROTEUS code, has been integrated into GRIFFIN to prepare cross sections for thermal reactor applications with heterogeneous geometries. Additional improvements have been made by implementing an on-the-fly slowing down method, a double heterogeneity treatment capability, and updating the procedure to generate the fine multigroup library. The cross section preparation capability in GRIFFIN was verified for graphite-moderated TRISO fuel-based reactor benchmark problems: unit-cell problems of VHTR and EMPIRE micro reactor and HTTR assembly problems. Eigenvalues and multigroup cross sections of GRIFFIN agreed very well with those of the continuous-energy Monte Carlo code Serpent2 within 200 pcm in eigenvalue and 2% in cross sections. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Verification and Validation of High Explosive Reactive Burn Models Implemented in LANL's EAP and LAP Code Base

Reactive burn models represent a significant leap in high explosive (HE) modeling capability. The first generation of engineering models of HE detonation are called programmed burn models and they are largely based on the distance between a prescribed detonation point and each zone in a simulation. There have been many advancements to programmed burn models over the years and when the assumptions upon which they are based are met, a properly tuned programmed burn model can be highly accurate but if any of their assumptions is not met, as is the case for corner turning or weakly initiated HE burn, they will give the wrong answer. Reactive burn models represent an entirely new way of modeling HE burn. They use the local conditions of a zone – e.g. temperature, pressure or density – as calculated by a hydrocode to determine if and when the zone is going to detonate and if so, how rapidly. This difference opens up an entirely new set of capabilities for HE modeling. It makes it possible to accurately and predictively model phenomena like the effect of confinement and the formation of dead zones. Reactive burn models have seen sustained development effort at LANL for at least the last decade but several recent developments make it timely to transition reactive burn models from a research topic to a production tool. The main goal of this milestone is to facilitate and accelerate the adoption of reactive burn as a commonly available modeling option, with recommendations on the resolution that will be required and uncertainties associated with their modeling choices. To achieve this, we have performed verification, validation, and uncertainty quantification (UQ) assessments of AWSD and SURF/SURFplus in xRage and FLAG on a variety of different problems.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗