Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “unstructured”

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 415 records · Page 23

SIERRA Low Mach Module: Fuego User Manual (V.5.10)

The SIERRA Low Mach Module: Fuego, henceforth referred to as Fuego, is the key element of the ASC fire environment simulation project. The fire environment simulation project is directed at characterizing both open large-scale pool fires and building enclosure fires. Fuego represents the turbulent, buoyantly-driven incompressible flow, heat transfer, mass transfer, combustion, soot, and absorption coefficient model portion of the simulation software. Using MPMD coupling, Scefire and Nalu handle the participating-media thermal radiation mechanics. This project is an integral part of the SIERRA multi-mechanics software development project. Fuego depends heavily upon the core architecture developments provided by SIERRA for massively parallel computing, solution adaptivity, and mechanics coupling on unstructured grids.

97 MATHEMATICS AND COMPUTING↗

An Investigation of Variables Affecting Plutonium Hydriding (20210946DI): Surface Phase Transitions and Reactions

Tusas was originally developed for phase field simulation of solidification (see [33, 6, 32]) and evolution of polycrystalline solid mechanics (see [16]). Tusas is particularly developed for PDEs compatible with structured or unstructured Lagrange (nodal) finite element discretizations and explicit (Euler) or implicit (Euler, Trapezoid, BDF2) temporal discretizations.

97 MATHEMATICS AND COMPUTING↗

Enhancing Risk Analysis of Accidental Release Using CFD Modeling and Machine Learning

This project report presents the results from the AU31 project that is entitled “Enhancing Risk Analysis of Accidental Release Using CFD Modeling and Machine Learning” performed between FY2021 and FY2022. The technical objective of this project is to enhance risk analysis of accidental release of radioactive materials and provide novel risk forecasting and capabilities for assessing unplanned hazardous chemical releases for the U.S. Department of Energy (DOE). The approach includes a cross-platform finite element analysis of airborne chemical release, machine learning (ML) for simulation and forecasting, and a big data approach to data management, transformation and analysis capable of handling petabyte-scale unstructured and structured data. The outcome of these findings were assembled and accessible as a final output via an intuitive web hosted geospatial-based user interface for demonstration and general dissemination. These new results also provide a mechanism of combining this knowledge in a risk analysis framework tool. The geospatial analysis and output can utilize 5 years of qualified meteorological data to illustrate projected plume footprints that are generated from computational aspects of this study. Specifically, we show the development of a geospatial risk analysis capabilities to respond to and plan for accidental release of radioactive materials and to provide a risk forecasting tool for assessing hazardous chemical releases for the DOE. The focus of this effort is germane to common practices in risk analysis; and it is highly relevant to the timely and transparent release of information related to materials at risk including hazardous materials and chemical vapors and rapid assessment of the potential impacts on the workforce at active sites across the DOE complex.

42 ENGINEERING↗

CrossLink: General Overview [Slides]

Problem: Traditional mesh generation approaches are labor intensive and have limited robustness when applied to parametric design exploration and optimization of complex geometries. While automatic mesh generation approaches exist, they tend to generate tetrahedral or mixed-hybrid meshes which are generally unsuitable for physics applications with strong shock waves, thin boundary layers, and strong gradients. In addition, simulations sizes in the billions of cells are becoming more common with traditional mesh generation methods quickly reaching scalability limits. Solution: CrossLink offers a topology-based mesh generation approach with unstructured block-filling methods and a scalable mesh generation engine. In addition, CrossLink incorporates a python-based API for seamless workflow integration and robust repeatability of the geometry handling and mesh generation process. This makes it ideal for parametric design study and optimization of complex geometries. Finally, future versions of CrossLink will offer a parametric mesh capability that optimizes a high-order mesh and enables reconstruction of the final mesh in memory by the physics solver.

97 MATHEMATICS AND COMPUTING↗

Improving the representation of isopycnal mixing in E3SM (Final Report)

The DOE-supported Energy Exascale Earth System Model (E3SM) is a major attempt by the Department of Energy to develop a new Earth System Model using unstructured grids, which allow for the model to concentrate resolution where it is most needed and avoid numerical artifacts experienced by previous Earth System Models (ESMs) using regular grids in which the lateral resolution becomes extremely fine as the resolution shrinks at the poles. This effort required rewriting many of the algorithms previously developed for regular grids. One algorithm that was problematic in the first version of E3SM was the representation of mixing due to turbulent ocean eddies with scales smaller than the model grid. These mesoscale eddies are the primary way by which tracers are stirred along density surfaces in the ocean interior. In particular, the process of isopycnal mixing exchanges fresher, more oxygenated waters from polar regions with saltier, nutrient-rich and oxygen poor waters from tropical regions. Previous work in our group has shown that this mechanism is vitally important for bringing oxygen into poorly ventilated tropical regions (Gnanadesikan et al., 2012, 2013) and plays a significant role in determining the uptake of anthropogenic carbon dioxide (Gnanadesikan et al., 2015). However, in E3SMv1 this process was turned off, as turning it on caused the model to become unstable. Additionally, the rate of isopycnal mixing is directly proportional to a mixing coefficient $A_{Redi}$ whose value varies from less than 400 m 2 /s to 2000 m 2 /s across contemporary climate models. The main thrusts of this proposal therefore were to: 1.) Improve the numerics of isopycnal mixing in E3SM; 2.) Identify ways in which climate models are sensitive to the isopycnal mixing coefficient; and 3.) Develop better representations of the isopycnal mixing coefficient.

54 ENVIRONMENTAL SCIENCES↗

Truchas-Flow

The truchas flow algorithm is designed for simulating the behavior of incompressible (but not necessarily constant density) fluids on unstructured meshes while allowing for solid and void phases to serve as moving boundaries. The interface between fluids, solids and void is tracked via a geometric volume of fluid approach. This allows for a single fluid formulation of the governing equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Demonstration of acoustic monitoring for structural health of microreactors: Through use of neural networks and resonant ultrasound spectroscopy

Nuclear microreactors prioritize modularity and portability and are intended to be a cost-effective technology for non-conventional nuclear markets. As such, the development of microreactors into a safe and feasible solution for energy security applications will necessitate the development of non-destructive technologies to monitor the integrity of inaccessible reactors components during operation. This demonstration applies linear and nonlinear acoustic techniques, in combination with machine learning, to detect and classify mechanical changes (stress and damage) in a test article which are broadly representative of potential operating challenges within a functioning microreactor. All necessary data has been collected for this demonstration, with minor experimental issues identified that can be addressed in follow-on work. Motivated by the expected conditions within a functioning microreactor, we have demonstrated our monitoring techniques on a core-block-like test article using an unstructured excitation source that approximates the noisy acoustic environment expected during reactor operation. At all stress states mechanically applied to the test article, a machine learning model using an artificial neural network was able to classify with 100% accuracy whether a 3D laser vibrometry point measurement was made on an intact or artificially defective test article. Further, model predictions about whether the defect interface was rough or smooth were 95% accurate, indicating the ability of acoustic techniques to recover defect characteristics. Resonant ultrasound spectroscopy (RUS) was also applied to the dataset to provide further quantitative insights about material properties. RUS analysis was ultimately hampered by several minor experimental and data issues, limiting results to certain cases for this demonstration. Last, analysis using nonlinear RUS exhibited sensitivity to changing levels of applied stresses for each intact and defective state. As presented in this demonstration, acoustic monitoring exhibits sensitivity to stress changes, which are of concern due to high thermal gradients expected during startup and operation. Further, our techniques distinguish between measurements made on intact and damaged test articles.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Athena-I Modeling with MCNP6.3 [Poster]

This objective of this project is to develop a model of the Athena-I experiment by creating unstructured mesh (UM) geometries in CUBIT and to perform calculations of the models in Monte Carlo N-Particle (MCNP) transport code version 6.3 (MCNP6.3).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Oktavian Modeling and Analysis with MCNP6.3 [Poster]

Project Goals: Apply MCNP6 variance reduction (VR) techniques to improve the Oktavian benchmark experiment calculations; Verify VR methods and electron transport with MCNP6.3 unstructured mesh (UM) geometry. All MCNP calculations in this poster are for mode n, p, e.

42 ENGINEERING↗

Tidal Resource Gaps Analysis Technical Report

The tidal resource gaps project was created to address a growing body of evidence that models underpredict tidal current speeds compared to measurements at a number of the top-ranking tidal energy sites. In response, this project compared opportunistic tidal power density measurements from 16 tidal energy hot spot sites with estimates from resource assessment data to identify discrepancies. To improve the accuracy of resource estimates from model data, updated data from eight improved model simulations were obtained. Model improvements included grid refinement, domain coupling, and the use of unstructured or nested grids. New resource estimates were calculated from the updated models, and these data were used to update the tidal hot spots-a list of promising tidal energy sites around the United States.

16 TIDAL AND WAVE POWER↗

Niowave Neutron Source Converter: Lead-Bismuth-Eutectic (LBE) Windowless Target Design and Evaluation

Los Alamos National Laboratory (LANL) is working with Niowave on the design and evaluation of their lead-bismuth-eutectic (LBE) windowless target (i.e., neutron source converter). Niowave plans to use 200 kW electron beam at 40 MeV beam energy to produce neutrons by photonuclear reaction with LBE. Then, the neutrons undergo fission at the surrounding uranium target assembly (UTA) to produce Molybdenum 99 (Mo-99) as a fission product, which eventually decays to Technetium-99 (Tc-99m). Tc-99m is one of important radioisotopes that is used for medical diagnostics. LANL conducted 3D multiphysics analysis for the Niowave neutron converter design and provided design assessment in thermohydraulic aspects. LANL conducted radiation transport calculations using Monte-Carlo N-Particle (MCNP) code with unstructured meshing scheme. The 3D volumetric heating profiles in the LBE and Stainless-Steel (SS) housing were imported into multiphase computational fluid dynamics (CFD) to obtain 3D temperature profiles of LBE and SS through conjugate heat transfer (CHT) analysis. The key findings are: LBE film thickness at the center of the beam is approximately 1.6 cm with a maximum LBE velocity of approximately 1.8 m/s, which is below a 2 m/s limit to avoid erosion issues on supporting structures; Heat deposition in the LBE peaks at ~1 cm depth from the LBE free-surface because of the forward interactions of electron, photon, and neutron with LBE; LBE maximum temperature is ~360 °C, which is below LBE evaporation initiative temperature, ~450 °C; LBE-SS interface temperature is ~350 °C, which is below the safety thermal limit to prevent severe corrosion on SS. The results indicate that Niowave’s neutron converter design satisfies both hydraulic and thermal criteria for safe operation. By virtue of such computational analysis, Niowave can move toward establishing an experimental setup to experimentally test their LBE neutron converter. The following sections describe the detailed work done by LANL.

43 PARTICLE ACCELERATORS↗

Multi-Fidelity Scheme for Accelerating 4D Finite Element Analysis for Mircoreactors

Next-generation microreactors are currently being designed to be operated terrestrial and extraterrestrial for remote surface power production. These systems will provide an alternative source of carbon-free energy that is versatile and can be utilized for various applications. These applications of microreactors have prompted the usage of high-fidelity unstructured finite element (FE) based approaches to provide time-dependent solutions for multiple physics fields. Computing these high-fidelity (4-D) solutions for multiple design iterations, physics fields, and transient events requires an immense computational cost. These high-fidelity solutions are computational expensive and the cost can be reduced through the implementation of less accurate low-fidelity solutions. Unlike the current fleet of commercial nuclear reactors, these next-gen systems present challenges due to the material and physical limitations required. Due to these constraints, the brute force technique of parameterizing important system characteristics determines whether a design meets the project objective. A system such as a nuclear reactor could have thousands of design parameters that affect system performance. In order to analyze the entire parameter space of such a complex system would require millions of CPU hours and countless design iterations. This costly approach is not practical due to regulatory and budget limitations. In this proposal, an approach that utilizes hybrid high-fidelity and low-fidelity FE models to reduce the computational cost of evaluating these applications in 4-D will be presented. The accuracy of the high-fidelity model and the computational efficiency of the low-fidelity model are taken advantage of to produce a solution that closely resembles the full order high-fidelity solution. By utilizing an unconverged coarse FE mesh, operation limits such as temperature, structural loading, etc., can be evaluated in an accelerated fashion and then can be spatially interpolated onto a finer mesh. The resulting error arising from the coarse mesh can be mitigated with a discrepancy function that actively quantifies and corrects the error in the coarse solution. This discrepancy function can be periodically updated with high-fidelity calculations across the temporal domain thus requiring less iterations on the fine mesh. As a result, the computational cost can be reduced for the evaluation and design iteration of microreactors. The proposal for this research contains three sections: Section 2 provides a literature review, Section 3 outlines the methodology for the proposed multi-fidelity scheme, and Section 4 displays preliminary results of the proposed research.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Enhance BISON Metal Fuel Transient Models and Enable Assessment Based on Out-of-Pile Transient Experiments

In this report, the activities and achievements made by Argonne National Laboratory for the Nuclear Energy Advanced Modeling and Simulation (NEAMS) BISON code metallic fuel validation and verification (V&V) project in FY2023, with a focus on the capability development and assessment case development for the out-of-pile transient experiments, are summarized. The liquid cladding penetration model was improved based on a comprehensive set of fuel behavior test apparatus (FBTA) test results. The new model was implemented into BISON along with a simple model to account for fuel melting effects. The metallic fuel meshing capability was also upgraded to allow local refinement using unstructured triangular mesh. Based on the previous success in developing the BISON-FIPD integration powered BISON assessment case, this database enabled BISON evaluation and V&V capability was expanded to the out-of-pile transient experiments by introducing BISON-OPTD integration. A series of initial out-of-pile transient assessment cases were implemented into BISON to provide powerful platforms to evaluate current BISON’s capabilities in simulating metallic fuel performance during out-of-pile transients. The established capabilities will also benefit the calibration and V&V of future advanced mechanistic models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Two MCNP Models for Computational Performance Benchmarking

The purpose of this report is to describe two non-trivial MCNP models and execution configurations that are suitable to characterizing computational performance. To that end, this report uses a constructive solid geometry (CSG) representation of the Oak Ridge National Laboratory (ORNL) Pool Critical Assembly (PCA) [1–3] to perform a k-eigenvalue calculation and an unstructured mesh (UM) representation of the International Commission on Radiological Protection (ICRP) publication 145 (ICRP145) male human phantom [4, 5] to perform a fixed-source calculation assuming that a 1 MeV photon source is distributed throughout the phantom’s liver.

97 MATHEMATICS AND COMPUTING↗

Heating, DPA, and He production in the Proton Beam Window Assembly for the Second Target Station

This report details the analysis of the heating rates, dpa values, and He production in the Proton Beam Window (PBW), the assembly around the PBW, and the shielding above it. Neutronics calculations use a detailed Unstructured Mesh (UM) converted from the solid CAD model of the PBW design of June 2022. The rest of the geometry is taken from the CSG master model (as of June 2022). For this analysis, a proton beam profile with a footprint of ≈90 cm 2 has been used. This proton profile was created by a Ring to Second Target magnet array configuration using octupole magnets. Maximum values in the PBW (aluminum) are 5.4 J/cc/pulse heating rate, 1.1 dpa per year, and 660 appm He-production per year. These results within ~ 10% with the previous calculations from June 2020. Average heating and dpa rates have been calculated for each component in the PBW assembly and will be used to evaluate and design the cooling needs and make decisions on removable components for the PBW. It is desirable to optimize and simplify the PBW design by only replacing components that are at the end of their lifetime. Additional analysis for the energy deposition in the PBW is performed for two super-Gaussian beam profiles with a footprint of 60 cm 2 and 30 cm 2 respectively. Also, an estimate of the PBW lifetime is performed for two realistic quadrupole beam profiles.

43 PARTICLE ACCELERATORS↗

The International Conference on Surface Plasmon Photonics 10

The International Conference on Surface Plasmon Photonics (SPP Conference) is widely regarded as the premier global conference in nanoscale photonics. The SPP10 meeting was held in the United States for the first time, at Rice University in Houston, Texas, on May 21-26, 2023. This conference was an outstanding venue for communicating research breakthroughs and accomplishments in the represented fields, as well as for visibility, networking, and recruiting, for students, postdoctoral fellows, and faculty throughout the US and the world as they attended this conference. The SPP10 conference hosted sessions on fundamental science, innovations, and applications in modern photonics with high potential for commercialization. Plenary sessions, invited talks from leading scientists in the field, and contributed talks from the larger photonics community on highly multidisciplinary topics of Quantum Plasmonics and Photonics; Energy Harvesting and Plasmonic Chemistry; Active Photonics and Devices; Ultrafast and Nonlinear Phenomena; and Machine Learning/Artificial Intelligence. The meeting schedule includes 3 plenary speakers, 31 Invited talks and 58 contributed talks. There were three well attended poster sessions which included lively discussions. Lunch time in the University Colleges allowed unstructured time for additional networking and discussions. The attendees represented the spectrum of endeavor in this field coming from academia, industry, and government laboratories, both U.S. and foreign scientists, senior researchers, young investigators, and students. The co-chairs of the conference were Prof. Naomi Halas and Prof. Peter Nordlander. The grant from DoE supported student’s and early career scientist’s registration and housing on campus.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Quantum Search Approaches to Sampling-Based Motion Planning

In this paper, we present a novel formulation of traditional sampling-based motion planners as database-oracle structures that can be solved via quantum search algorithms. We consider two complementary scenarios: for simpler sparse environments, we formulate the Quantum Full Path Search Algorithm (q-FPS), which creates a superposition of full random path solutions, manipulates probability amplitudes with Quantum Amplitude Amplification (QAA), and quantum measures a single obstacle free full path solution. For dense unstructured environments, we formulate the Quantum Rapidly Exploring Random Tree algorithm, q-RRT, that creates quantum superpositions of possible parent-child connections, manipulates probability amplitudes with QAA, and quantum measures a single reachable state, which is added to a tree. As performance depends on the number of oracle calls and the probability of measuring good quantum states, we quantify how these errors factor into the probabilistic completeness properties of the algorithm. We then numerically estimate the expected number of database solutions to provide an approximation of the optimal number of oracle calls in the algorithm. We compare the q-RRT algorithm with a classical implementation and verify quadratic run-time speedup in the largest connected component of a 2D dense random lattice. We conclude by evaluating a proposed approach to limit the expected number of database solutions and thus limit the optimal number of oracle calls to a given number.

97 MATHEMATICS AND COMPUTING↗

Processing MCNP Elemental Edit Outputs

The Monte Carlo N-Particle (MCNP) transport code version 6 (also known as MCNP6) has the capability for tracking particles on unstructured mesh (UM) geometry models embedded into constructive solid geometry (CSG) cells. A UM geometry is a collection of elements representing a solid geometry. The first step of MCNP UM modeling is using other software packages to create a finite element mesh representation of a solid 3D geometry. Computer-aided design (CAD) or computer-aided manufacturing (CAM) software is typically used to create a solid geometry model, which is later imported into mesh generation software to create a UM model. The MCNP UM feature was originally designed for models generated by the Abaqus/CAE software. The MCNP code version 6.0 and later can process UM models formatted as Abaqus input files. MCNP can process a UM model consisting of several different element types including linear tetrahedral or hexahedral elements and calculate quantities of interest such as flux and energy deposition at elements. An MCNP UM simulation provides high-fidelity elemental edit (i.e., tally) outputs, which can be further used in multiphysics calculations. The MCNP UM feature was used for multiphysics simulations where quantities of interest calculated by MCNP are used as inputs for heat transfer calculations in Abaqus. MCNP6.3 can produce two types of elemental edit output (EEOUT) file formats: ASCII and HDF5. An EEOUT file type must be requested on an EMBED card while output type (flux or energy deposition) must be requested on an EMBEE card. We wrote Python3 scripts to extract energy deposition values in an ASCII or HDF5 EEOUT file and compute a heat flux profile for an Abaqus heat transfer calculation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗