Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulation techniques”

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Reconfigurable metamaterial absorber as an optical switch based on organic-graphene control

In this paper, we have suggested an effective optical switch with the cross fractal shape. It contains the main cross element that connected to four parasitic elements by four organic’s material junctions. This absorber is modified as a tunable element by using the graphene layer which can be considered for compensating the physical distortion. Each organic and parasitic load is placed over a graphene layer. These graphene layers separated with a SiO{sub 2} thin film. We have implemented a SiO{sub 2} layer as a spacer between the graphene and metal ground layer. The organic materials are popular for their various resistances value. We have applied a bulky dielectric slab with two different conductivity values for high and low resistance (HRS and LRS). This absorber works at the range of 90–100 THz for mid-infrared spectrum and the bandwidth of 4%. In the HRS mode the reflection value is about − 35 up to − 25 dB. When the organic material switched to LRS mode, this value drastically increased to − 5 dB. The designed absorber has both the tunable and switchable attributes. We have simulated this switch by the finite integral technique method with CST microwave studio.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Resonant anomaly detection with multiple reference datasets

An important class of techniques for resonant anomaly detection in high energy physics builds models that can distinguish between reference and target datasets, where only the latter has appreciable signal. Such techniques, including Classification Without Labels (CWoLa) and Simulation Assisted Likelihood-free Anomaly Detection (SALAD) rely on a single reference dataset. They cannot take advantage of commonly available multiple datasets and thus cannot fully exploit available information. In this work, we propose generalizations of CWoLa and SALAD for settings where multiple reference datasets are available, building on weak supervision techniques. We demonstrate improved performance in a number of settings with realistic and synthetic data. As an added benefit, our generalizations enable us to provide finite-sample guarantees, improving on existing asymptotic analyses.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Feedforward Control of Fuel Distribution on Advanced Dual-Fuel Engines with Varying Intake Valve Closing Timings

This study examines the dynamics and control of an engine operated with late intake valve closure (LIVC) timings in a dual-fuel combustion mode. The engine features a fuel delivery system in which diesel is direct-injected and natural gas is port-injected. Despite the benefits of LIVC and dual-fuel strategy, combining these two techniques resulted in efficiency losses due to the variability of the combustion process across cylinders. The difference in power production across cylinders ranges from 9% at an IVC of 570°ATDC* to 38% at an IVC of 620 °ATDC and indicates an increasingly uneven fuel distribution as the intake valve remains open longer in the compression stroke. This paper describes an approach for controlling the amount of fuel injected into each cylinders’ port of an inline six- cylinder heavy-duty dual-fuel engine to minimize the variations in fuel distribution across cylinder. Since measuring the actual fuel reaching each cylinder is not a practical option on stock engines, this study aimed to produce a feedforward control technique. Using a model developed on an engine simulation software, the optimal fuel injection amounts at each cylinder port across 125 operating conditions (consisting of sweeps of five engine speeds, five engine loads, and five IVC timings) were determined by using PI controllers that adjust the amount of fuel injected at each port based on the amount of fuel reaching the cylinders. The results from the PI controller were expressed as a function of the operating condition (engine speed, engine load, and IVC timing); and consequently, could be reproduced without relying on a measurement of the amount of fuel reaching each cylinder. The feedforward control strategy, unlike the PI controller, is easily implementable on a stock engine and significantly reduces the variations in fuel distribution across cylinders.

Kassa, Mateos↗
DOE OSTI2022

A Moving Embedded Boundary Approach for the Compressible Navier-Stokes Equations in a Block-Structured Adaptive Refinement Framework

A computational technique has been developed to perform compressible flow simulations involving moving boundaries using an embedded boundary approach within the block-structured adaptive mesh refinement (SAMR) framework of AMReX [1], [91], [92]. We leverage the SAMR capability to obtain quantitatively accurate results whilst using robust, second-order finite volume schemes. A conservative, unsplit, cut-cell approach is utilized and a ghost-cell approach is developed for computing the flux on the moving, embedded boundary faces. A third-order least-squares formulation has been developed to compute the wall velocity gradients, and was found to significantly improve the performance of the solver in terms of the quantitative comparison of surface quantities such as the skin friction coefficient. Various test cases are performed to validate the method, and compared with analytical, experimental, and other numerical results in literature. Inviscid and viscous test cases are performed that span a wide regime of flow speeds - acoustic (harmonically pulsating sphere), smooth flows (expansion fan created by a receding piston) and flows with shocks (shock-cylinder interaction, shock-wedge interaction, pitching NACA 0012 airfoil and shock-cone interaction). A closed system with moving boundaries - an oscillating piston in a cylinder, showed that the percentage error in mass within the system decreases with refinement, demonstrating that the numerical scheme is conservative with grid refinement, but is not discretely conservative. Viscous test cases involve that of a horizontally moving cylinder at Re = 40, an inline oscillating cylinder at Re = 100, and a transversely oscillating cylinder at Re = 185. The judicious use of adaptive mesh refinement with appropriate refinement criteria to capture the regions of interest leads to well-resolved flow features, and good quantitative comparison is observed with the results available in literature.

adaptive refinement↗
DOE OSTI2025

Burnup Monitoring for Pebble Bed Reactor Systems

A pebble burnup monitoring system is a required component for domestic reactor safety and safeguards applications associated with pebble bed reactors (PBRs). One of the main requirements of a PBR burnup monitoring system is that it needs to be capable of rapid measurements to assess the burnup of each individual pebble to determine whether to recirculate it in the reactor or discard it as spent fuel. This report considers three different approaches for a burnup monitoring system for pebbles discharged from the reactor core in a pebble bed modular reactor-400: • passive gamma spectrometry measurement, • passive neutron coincidence measurement, and • active neutron counter based on the differential die-away technique. Conceptual designs have been created for each of these detectors, and preliminary analysis has been performed using Monte Carlo N-Particle and Oak Ridge Isotope Generation code simulations. The advantages and practical limitations (e.g., high radiation background) of each system were identified. Simulations suggest that each of the three measurement techniques can be successfully employed to distinguish between pebbles based on their number of passes through the core and to quantify the burnup of pebbles.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

A Scintillation-Fiber Based Cosmic-ray Muon Tomography System for Imaging Dry Storage Cask

Because of the growth of the nuclear power industry in the United States and the policy decision to ban reprocessing of commercial spent nuclear fuel, the spent fuel inventory at commercial reactor sites has been increasing. Spent nuclear fuel presents extra complexity of security and proliferation issues. These fuel assemblies are practically inaccessible for inspection purposes, as reopening a DSC would require special facilities and would be tremendously expensive. Cosmic-ray muon Tomography has been studied in the past for various applications and was demonstrated to be a good non-destructive assay technique for high-Z materials, such as spent nuclear fuel assemblies in dry storage casks (DSC). Compared to conventional imaging techniques, muon imaging has some advantages due to its negligible health risk and high penetrating power. In this paper, we report our simulation study on imaging a DSC using muon Tomography technique. The Point of Closest Approach (POCA) and the Maximum-Likelihood Expectation- Maximization (ML-EM) algorithms are implemented and evaluated. A scaled-down prototype muon tomography system is proposed and studied. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Developing Ultrahigh-Resolution E3SM Land Model for GPU Systems

Designing and refactoring complex scientific code, such as the E3SM land model (ELM), for new computing architectures is challenging. This paper presents design strategies and technical approaches to develop a data-oriented, GPU-ready ELM model using compiler directives (OpenACC/OpenMP). We first analyze the datatypes and processes in the original ELM code. Then we present design considerations for ultrahigh-resolution ELM (uELM) development for massive GPU systems. These techniques include the global data-oriented simulation workflow, domain partition, code porting and data copy, memory reduction, parallel loop restructure and flattening, and race condition detection. We implemented the first version of uELM using OpenACC targeting the NVidia GPUs in the Summit supercomputer at Oak Ridge National Laboratory. During the implementation, we developed a software tool (named SPEL) to facilitate code generation, verification, and performance tuning using these techniques. The first uELM implementation for Nvidia GPUs on Summit delivered promising results: 1) over 98% of the ELM code was automatically generated and tuned by scripts. Most ELM modules had better computational performances than the original ELM code for CPUs. The GPU-ready uELM is more scalable than the CPU code on fully-loaded Summit nodes. Example profiling results from several modules are also presented to illustrate the performance improvements and race condition detection. The lessons learned and toolkit developed in the study are also suitable for further uELM deployment using OpenMP on the first US exascale computer, Frontier, equipped with AMD CPUs and GPUs.

Schwartz, Peter↗
DOE OSTI2016

Material point method modeling in oil and gas reservoirs

A computer system and method of simulating the behavior of an oil and gas reservoir including changes in the margins of frangible solids. A system of equations including state equations such as momentum, and conservation laws such as mass conservation and volume fraction continuity, are defined and discretized for at least two phases in a modeled volume, one of which corresponds to frangible material. A material point model technique for numerically solving the system of discretized equations, to derive fluid flow at each of a plurality of mesh nodes in the modeled volume, and the velocity of at each of a plurality of particles representing the frangible material in the modeled volume. A time-splitting technique improves the computational efficiency of the simulation while maintaining accuracy on the deformation scale. The method can be applied to derive accurate upscaled model equations for larger volume scale simulations.

02 PETROLEUM↗
DOE OSTI2023

Fast and accurate reduced-order modeling of a MOOSE-based additive manufacturing model with operator learning

One predominant challenge in additive manufacturing (AM) is to achieve specific material properties by manipulating manufacturing process parameters during the runtime. Such manipulation tends to increase the computational load imposed on existing simulation tools employed in AM. The goal of the present work is to construct a fast and accurate reduced-order model (ROM) for an AM model developed within the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, ultimately reducing the time/cost of AM control and optimization processes. Our adoption of the operator learning (OL) approach enabled us to learn a family of differential equations produced by altering process variables in the laser’s Gaussian point heat source. More specifically, we used the Fourier neural operator (FNO) and deep operator network (DeepONet) to develop ROMs for time-dependent responses. Furthermore, we benchmarked the performance of these OL methods against a conventional deep neural network (DNN)-based ROM. Ultimately, we found that OL methods offer comparable performance and, in terms of accuracy and generalizability, even outperform DNN at predicting scalar model responses. The DNN-based ROM afforded the fastest training time. Furthermore, all the ROMs were faster than the original MOOSE model yet still provided accurate predictions. FNO had a smaller mean prediction error than DeepONet, with a larger variance for time-dependent responses. Unlike DNN, both FNO and DeepONet were able to simulate time series data without the need for dimensionality reduction techniques. Finally, the present work can help facilitate the AM optimization process by enabling faster execution of simulation tools while still preserving evaluation accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Optimization of laser-driven cylindrical implosions on the OMEGA laser

Laser-driven cylindrical implosions were conducted on the OMEGA laser as part of the laser-driven mini-MagLIF (Magnetized Liner Inertial Fusion) Campaign. Gated x-ray images were analyzed to infer shell trajectories and study the energy coupling in these implosions. Two-dimensional and three-dimensional HYDRA simulations were performed and post-processed to produce synthetic x-ray self-emission images for comparison. An analysis technique, which could be applied to both experimental and simulated x-ray images, was developed to characterize the shape and uniformity of the implosion. The analysis leads to a measurement of the average implosion velocity and axial implosion length, which can then be used to optimize the beam pointing and energy balance for future experiments. Discrepancies between simulation results and experiments allude to important physical processes that are not accounted for in the simulations. In 2-D simulations, the laser beam’s azimuthal angle of incidence is not included because the $\phi$-direction is not simulated, thus energy absorption is over-predicted. The 3-D simulation results are more consistent with the experiments, but the simulations do not include a calculation of crossbeam energy transfer or non-local thermal transport, which affect the energy coupled to the implosion. Finally, by appropriately adjusting the simulated energy balance and flux limit, the simulations can accurately model the experiments, which have achieved uniform implosions over a 700-μm-long region at velocities of approximately 200 km=s.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Characterizing Tradeoffs in Memory, Accuracy, and Speed for Chemistry Tabulation Techniques

Chemistry tabulation is a common approach in practical simulations of turbulent combustion at engineering scales. Linear interpolants have traditionally been used for accessing precomputed multidimensional tables but suffer from large memory requirements and discontinuous derivatives. Higher-degree interpolants address some of these restrictions but are similarly limited to relatively low-dimensional tabulation. Artificial neural networks (ANNs) can be used to overcome these limitations but cannot guarantee the same accuracy as interpolants and introduce challenges in reproducibility and reliable training. These challenges are enhanced as the physics complexity to be represented within the tabulation increases. Here, we assess the efficiency, accuracy, and memory requirements of Lagrange polynomials, tensor product B-splines, and ANNs as tabulation strategies. We analyze results in the context of nonadiabatic flamelet modeling where higher dimension counts are necessary. While ANNs do not require structuring of data, providing benefits for complex physics representation, interpolation approaches often rely on some structuring of the table. Interpolation using structured table inputs that are not directly related to the variables transported in a simulation can incur additional query costs. This is demonstrated in the present implementation of heat losses. We show that ANNs, despite being difficult to train and reproduce, can be advantageous for high-dimensional, unstructured datasets relevant to nonadiabatic flamelet models. Furthermore we demonstrate that Lagrange polynomials show significant speedup for similar accuracy compared to B-splines.

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

Strong Lensing Parameter Estimation on Ground-Based Imaging Data Using Simulation-Based Inference

Current ground-based cosmological surveys, such as the Dark Energy Survey (DES), are predicted to discover thousands of galaxy-scale strong lenses, while future surveys, such as the Vera Rubin Observatory Legacy Survey of Space and Time (LSST) will increase that number by 1-2 orders of magnitude. The large number of strong lenses discoverable in future surveys will make strong lensing a highly competitive and complementary cosmic probe. To leverage the increased statistical power of the lenses that will be discovered through upcoming surveys, automated lens analysis techniques are necessary. We present two Simulation-Based Inference (SBI) approaches for lens parameter estimation of galaxy-galaxy lenses. We demonstrate the successful application of Neural Posterior Estimation (NPE) to automate the inference of a 12-parameter lens mass model for DES-like ground-based imaging data. We compare our NPE constraints to a Bayesian Neural Network (BNN) and find that it outperforms the BNN, producing posterior distributions that are for the most part both more accurate and more precise; in particular, several source-light model parameters are systematically biased in the BNN implementation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Raman spectroscopy and electrical properties of InAs nanowires with local oxidation enabled by substrate micro-trenches and laser irradiation

The thermal gradients along indium arsenide nanowires were engineered by a combination of fabricated micro-trenches in the supporting substrate and focused laser irradiation. This allowed local spatial control of thermally activated oxidation reactions of the nanowire on the scale of the diffraction limit. The locality of the oxidation was detected by micro-Raman mapping, and the results were found to be consistent with numerical simulations of the temperature profile. Applying the technique to nanowires in electrical devices the locally oxidized nanowires remained conducting with a lower conductance as expected for an effectively thinner conducting core.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Notes on modeling and simulation

These notes present a high-level overview of how modeling and simulation are carried out by practitioners. The discussion is of a general nature; no specific techniques are examined but the activities associated with all modeling and simulation approaches are briefly addressed. There is also a discussion of validation and verification and, at the end, a section on why modeling and simulation are useful.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Modeling and Simulation of Austenitic Welds and Coarse-grained Specimens: Part II

The Pacific Northwest National Laboratory (PNNL) is conducting confirmatory research for the U.S. Nuclear Regulatory Commission (NRC) to evaluate commercially available nondestructive examination (NDE) modeling and simulation software used in the nuclear industry. Simulation results from ultrasonic testing (UT) models can inform the design and qualification of inspection techniques and help interpret inspection results. CIVA is a modeling and simulation package developed by the French Alternative Energies and Atomic Energy Commission (CEA). CIVA was selected for this study because it is readily available and has been used for NDE in the US commercial nuclear power industry. This report is focused on completing the efforts initiated in the previous PNNL report to evaluate UT modeling and simulation performance, reliability, and accuracy in relation to common inservice inspection (ISI) scenarios in nuclear power plants (NPP). This work will be used to provide guidance when establishing methods to perform and evaluate simulations for more standardized model implementation, simulation analysis, and interpretation of results.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Domain Adaptation for Measurements of Strong Gravitational Lenses

Upcoming surveys are predicted to discover galaxy-scale strong lenses on the magnitude of 105, making deep learning methods necessary in lensing data analysis. Currently, there is insufficient real lensing data to train deep learning algorithms, but training only on simulated data results in poor performance on real data. Domain adaptation can bridge the gap between simulated and real datasets. We adopt domain adaptation on the estimation of Einstein radius in simulated galaxy-scale gravitational lensing images. We evaluate two domain adaptation techniques - domain adversarial neural networks (DANN) and maximum mean discrepancy (MMD). We train on a source domain of simulated lenses and apply it to a target domain with emulation of DES survey conditions. We show that both domain adaptation techniques can significantly improve the model performance on the more complex target domain datasets. Our results show the potential of using domain adaptation to perform analysis on future survey data with a deep neural network trained on simulated data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Recovering missing CFD data for high-order discretizations using deep neural networks and dynamics learning

Data I/O poses a significant bottleneck in large-scale CFD simulations; thus, practitioners would like to significantly reduce the number of times the solution is saved to disk, yet retain the ability to recover any field quantity (at any time instance) a posteriori . The objective of this work is therefore to accurately recover missing CFD data a posteriori at any time instance, given that the solution has been written to disk at only a relatively small number of time instances. We consider in particular high-order discretizations (e.g., discontinuous Galerkin), as such techniques are becoming increasingly popular for the simulation of highly separated flows. To satisfy this objective, this work proposes a methodology consisting of two stages: 1) dimensionality reduction and 2) dynamics learning. For dimensionality reduction, we propose a novel hierarchical approach. First, the method reduces the number of degrees of freedom within each element of the high-order discretization by applying autoencoders from deep learning. Second, the methodology applies principal component analysis to compress the global vector of encodings. This leads to a low-dimensional state, which associates with a nonlinear embedding of the original CFD data. For dynamics learning, we propose to apply regression techniques (e.g., kernel methods) to learn the discrete-time velocity characterizing the time evolution of this low-dimensional state. In conclusion, a numerical example on a large-scale CFD example characterized by nearly 13 million degrees of freedom illustrates the suitability of the proposed method in an industrial setting.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Constant-depth circuits for dynamic simulations of materials on quantum computers

Abstract Dynamic simulation of materials is a promising application for near-term quantum computers. Current algorithms for Hamiltonian simulation, however, produce circuits that grow in depth with increasing simulation time, limiting feasible simulations to short-time dynamics. Here, we present a method for generating circuits that are constant in depth with increasing simulation time for a specific subset of one-dimensional (1D) materials Hamiltonians, thereby enabling simulations out to arbitrarily long times. Furthermore, by removing the effective limit on the number of feasibly simulatable time-steps, the constant-depth circuits enable Trotter error to be made negligibly small by allowing simulations to be broken into arbitrarily many time-steps. For an N -spin system, the constant-depth circuit contains only $\mathcal {O}(N^{2})$ O ( N 2 ) CNOT gates. Such compact circuits enable us to successfully execute long-time dynamic simulation of ubiquitous models, such as the transverse field Ising and XY models, on current quantum hardware for systems of up to 5 qubits without the need for complex error mitigation techniques. Aside from enabling long-time dynamic simulations with minimal Trotter error for a specific subset of 1D Hamiltonians, our constant-depth circuits can advance materials simulations on quantum computers more broadly in a number of indirect ways.

Materials simulation↗