Search NASA⌕ Search

SEARCH · Search NASA

Results for “Engine 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,153 records · Page 64

Resource-Optimized Fermionic Local-Hamiltonian Simulation on a Quantum Computer for Quantum Chemistry

The ability to simulate a fermionic system on a quantum computer is expected to revolutionize chemical engineering, materials design, nuclear physics, to name a few. Thus, optimizing the simulation circuits is of significance in harnessing the power of quantum computers. Here, we address this problem in two aspects. In the fault-tolerant regime, we optimize the R z and T gate counts along with the ancilla qubit counts required, assuming the use of a product-formula algorithm for implementation. We obtain a savings ratio of two in the gate counts and a savings ratio of eleven in the number of ancilla qubits required over the state of the art. In the pre-fault tolerant regime, we optimize the two-qubit gate counts, assuming the use of the variational quantum eigensolver (VQE) approach. Specific to the latter, we present a framework that enables bootstrapping the VQE progression towards the convergence of the ground-state energy of the fermionic system. This framework, based on perturbation theory, is capable of improving the energy estimate at each cycle of the VQE progression, by about a factor of three closer to the known ground-state energy compared to the standard VQE approach in the test-bed, classically-accessible system of the water molecule. The improved energy estimate in turn results in a commensurate level of savings of quantum resources, such as the number of qubits and quantum gates, required to be within a pre-specified tolerance from the known ground-state energy. We also explore a suite of generalized transformations of fermion to qubit operators and show that resource-requirement savings of up to more than 20 % , in small instances, is possible.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tabulated Fluid Properties Research Report

The Multiphysics Object-Oriented Simulation Environment (MOOSE) enables a wide range of advanced nuclear reactor simulations.[6] Under the guidance of MOOSE’s Thermal Hydraulics Team,I worked to expand the capabilities of Tabulated Fluid Properties (TFP) in the fluid properties module. The fluid properties module allows the user to determine a variety of fluid properties by interpolating points between tabulated data. I implemented the ability to use bilinear interpolation instead of bicubic interpolation for interpolating tabulated data. I also changed the method of variable set inversions to use a 2-dimensional Newton’s Method utility that I created. Variable set inversions are often done from (v,e) to (p,T), where v is specific volume, e is specific internal energy, p is pressure and T is temperature. New routines have also been added into TFP such that it can be used with more applications, such as the Navier Stokes and Thermal Hydraulics modules in MOOSE for Pronghorn[5] and RELAP-7[1] respectively. This work was spurred by interest from NASA in testing a Nuclear Thermal Propulsion (NTP) engine system. NTP engines have drastically different fluid properties throughout the engine and Tabulated Fluid Properties provides the flexibility needed to properly simulate and test these engines.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Multi-Federate Co-Convergence with HELICS

In co-simulation studies, convergence refers to the ability of different simulation tools involved to achieve a consistent and stable solution. Convergence is a critical aspect of co-simulation as it determines the accuracy of the results obtained from the simulation. Convergence in co-simulation studies depends on several factors, this includes system complexity, the accuracy of the models used, and the numerical methods employed by the simulation tools. It is essential to ensure that the coupling interfaces between the different tools are designed to allow for data exchange in a consistent and accurate manner. Hierarchical Engine for Large-scale Infrastructure Co-Simulation (HELICS) is an open-source co-simulation framework developed for the energy domain. This paper explores the convergence performance of a set of co-simulation use cases. We further explore the use of a co-convergence helper federate to help with co-simulation convergence. The convergence efficacy of several algorithms (both gradient-based and gradient-free) is tested against these use cases. Finally, the sensitivity of these algorithms to several factors, such as system scaling and others, is tested and detailed in this paper. Our results show that for a subset of use cases, the co-convergence helper federate is able to improve co-simulation convergence significantly.

co-convergence↗

The role of cool-flame fluctuations in high-pressure spray flames, studied using high-speed optical diagnostics and Large-Eddy Simulations

Here this work investigates the low- and high-temperature ignition and combustion processes, applied to the Engine Combustion Network Spray A flame, combining advanced optical diagnostics and large-eddy simulations (LES). Simultaneous high-speed (50 kHz) formaldehyde (CH 2 0) planar laser-induced fluorescence (PLIF) and line-of-sight OH* chemiluminescence imaging were used to measure the low- and high-temperature flame, during ignition as well as during quasi-steady combustion. While tracking the cool flame at the laser sheet plane, the present experimental setup allows detection of distinct ignition spots and dynamic fluctuations of the lift-off length over time, which overcomes limitations for flame tracking when using schlieren imaging. After significant development to improve LES prediction of the low-and high-temperature flame position, both during the ignition processes and quasi-steady combustion, the simulations were analyzed to gain understanding of the mixture variance and how this variance affects formation/consumption of CH 2 0. Analysis of the high-temperature ignition period shows that a key improvement in the LES is the ability to predict heterogeneous ignition sites, not only in the head of the jet, but in shear layers at the jet edge close to the position where flame lift-off eventually stabilizes. The LES analysis also shows concentrated pockets of CH 2 0, in the center of jet and at 20 mm downstream of the injector (in regions where the equivalence ratio is greater than 6), that are of similar length scale and frequency as the experiment (approximately 5–6 kHz). The periodic oscillation of CH 2 0 match the frequency of pressure waves generated during auto-ignition and reflected within the constant-volume vessel throughout injection. The ability of LES to capture the periodic appearance and destruction of CH 2 0 is particularly important because these structures travel downstream and become rich premixed flames that affect soot production.

42 ENGINEERING↗

Study of Plasma Liner Driven Magnetized Target Fusion Via Advanced Simulations

The feasibility of the plasma liner driven Magnetized Target Fusion (MTF) via terascale numerical simulations will be assessed. In the MTF concept, a plasma liner, formed by merging of a number (60 or more) of radial, highly supersonic plasma jets, implodes on the target in the form of two compact plasma toroids, and compresses it to conditions of the fusion ignition. By avoiding major difficulties associated with both the traditional laser driven inertial confinement fusion and solid liner driven MTF, the plasma liner driven MTF potentially provides a low-cost and fast R&D path towards the demonstration of practical fusion energy. High fidelity numerical simulations of full nonlinear models associated with the plasma liner MTF using state-of-art numerical algorithms and terascale computing are necessary in order to resolve uncertainties and provide guidance for future experiments. At Stony Brook University, we have developed unique computational capabilities that ideally suite the MTF problem. The FronTier code, developed in collaboration with BNL and LANL under DOE funding including SciDAC for the simulation of 3D multi-material hydro and MHD flows, has beenbenchmarked and used for fundamental and engineering problems in energy science applications. We have performed 3D simulations of converging supersonic plasma jets, their merger and the formation of the plasma liner, and a study of the corresponding oblique shock problem. We have studied the implosion of the plasma liner on the magnetized plasma target by resolving Rayleigh-Taylor instabilities in 2D and 3D and other relevant physics and estimate thermodynamic conditions of the target at the moment of maximum compression and the hydrodynamic efficiency of the method.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Performance Evaluation of Heterogeneous GPU Programming Frameworks for Hemodynamic Simulations

Preparing for the deployment of large scientific and engineering codes on upcoming exascale systems with GPU-dense nodes is made challenging by the unprecedented diversity of device architectures and heterogeneous programming models. In this work, we evaluate the process of porting a massively parallel, fluid dynamics code written in CUDA to SYCL, HIP, and Kokkos with a range of backends, using a combination of automated tools and manual tuning. We use a proxy application along with a custom performance model to inform the results and identify additional optimization strategies. At scale performance of the programming model implementations are evaluated on pre-production GPU node architectures for Frontier and Aurora, as well as on current NVIDIA device-based systems Summit and Polaris. Real-world workloads representing 3D blood flow calculations in complex vasculature are assessed. Our analysis highlights critical trade-offs between code performance, portability, and development time.

Martin, Aristotle↗

Suppressing simulation bias in multi-modal data using transfer learning

Abstract Many problems in science and engineering require making predictions based on few observations. To build a robust predictive model, these sparse data may need to be augmented with simulated data, especially when the design space is multi-dimensional. Simulations, however, often suffer from an inherent bias. Estimation of this bias may be poorly constrained not only because of data sparsity, but also because traditional predictive models fit only one type of observed outputs, such as scalars or images, instead of all available output data modalities, which might have been acquired and simulated at great cost. To break this limitation and open up the path for multi-modal calibration, we propose to combine a novel, transfer learning technique for suppressing the bias with recent developments in deep learning, which allow building predictive models with multi-modal outputs. First, we train an initial neural network model on simulated data to learn important correlations between different output modalities and between simulation inputs and outputs. Then, the model is partially retrained, or transfer learned, to fit the experiments; a method that has never been implemented in this type of architecture. Using fewer than 10 inertial confinement fusion experiments for training, transfer learning systematically improves the simulation predictions while a simple output calibration, which we design as a baseline, makes the predictions worse. We also offer extensive cross-validation with real and carefully designed synthetic data. The method described in this paper can be applied to a wide range of problems that require transferring knowledge from simulations to the domain of experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

CompLaB v1.0: a scalable pore-scale model for flow, biogeochemistry, microbial metabolism, and biofilm dynamics

Abstract. Microbial activity and chemical reactions in porous media depend on the local conditions at the pore scale and can involve complex feedback with fluid flow and mass transport. We present a modeling framework that quantitatively accounts for the interactions between the bio(geo)chemical and physical processes and that can integrate genome-scale microbial metabolic information into a dynamically changing, spatially explicit representation of environmental conditions. The model couples a lattice Boltzmann implementation of Navier–Stokes (flow) and advection–diffusion-reaction (mass conservation) equations. Reaction formulations can include both kinetic rate expressions and flux balance analysis, thereby integrating reactive transport modeling and systems biology. We also show that the use of surrogate models such as neural network representations of in silico cell models can speed up computations significantly, facilitating applications to complex environmental systems. Parallelization enables simulations that resolve heterogeneity at multiple scales, and a cellular automaton module provides additional capabilities to simulate biofilm dynamics. The code thus constitutes a platform suitable for a range of environmental, engineering and – potentially – medical applications, in particular ones that involve the simulation of microbial dynamics.

58 GEOSCIENCES↗

Analysis of ignition and stabilization modes in diesel spray flames using large eddy simulations and chemical explosive mode analysis

In this work, spray flames under diesel engine conditions are numerically investigated to understand the ignition and stabilization characteristics, based on state-of-the-art large eddy simulations (LES) and chemical explosive mode analysis (CEMA). In particular, CEMA is employed to identify such critical flame features as auto-ignition, premixed flame fronts, cool and strongly burning flames, and local extinction. In particular, a recently developed CEMA-based criterion for demarcating local combustion modes, that is, auto-ignition, diffusion-assisted ignition, and extinction, is applied to gain in-depth insights into the interactions of local chemical and transport processes. A series of one-dimensional non-premixed flamelets with varying air temperatures and scalar dissipation rates are first analyzed using CEMA with a focus on the role of each local mode in the transient ignition process. CEMA-based diagnostics are then applied to diesel spray flames at different ambient temperatures. The structure of local combustion modes is identified, while the contribution from each local mode to overall flame burning is quantified at different stages of the flame developing process. Effects of the ambient temperature on local combustion mode distribution, and subsequently on the ignition and stabilization mechanisms, are further examined. Results show that auto-ignition is the main mechanism for flame initiation and stabilization, while the role of diffusion is substantially enhanced as ambient temperature decreases.

assisted-ignition↗

FY22 Laboratory Directed Research and Development Annual Report

The Laboratory Directed Research and Development (LDRD) program yields foundational scientific research and development (R&D) essential to growing SRNL’s core competencies, in alignment with SRNL’s Strategic Plan to provide long-term benefits to the Department of Energy (DOE), the National Nuclear Security Administration (NNSA), and other customers and stakeholders. Five strategic goals are outlined in SRNL’s strategic plan: 1) Provide applied science and engineering for EM’s active clean-up sites and LM’s post closure management sites; 2) Provide science-based solutions for gaps identified in nonproliferation strategic vision and support the government in actives impacting national security; 3) Lead ST&E as the central technical authority for processing tritium loaded reservoirs and support production of plutonium pits; 4) Align science and energy security programs by focusing modern modeling, simulation, and data analytics tools on materials engineering and performance applications; 5) Build a workforce for the future.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Scale-up Testing of Advanced Polaris Membrane in CO2 Capture Technology

This final technical report describes work conducted by Membrane Technology and Research, Inc. (MTR) for the U.S. Department of Energy (DOE), National Energy Technology Lab (NETL) on the scale-up and testing of advanced Polaris™ membrane CO2 capture technology at the Technology Centre Mongstad (TCM) under award number DE-FE0031591. The work was performed from August 1, 2018 through January 31, 2023. The overall goal of this project was to design, build and operate an advanced Polaris membrane CO2 capture system at TCM. MTR was assisted in this project by Trimeric Corporation (Trimeric), an engineering design services company, the Carbon Capture Simulation for Industry Impact (CCSI2), a partnership among national laboratories, industry, and academic institutions, and the Technology Centre Mongstad (TCM), who provided the host site for the slipstream field test. This report details the work conducted to scale-up MTR’s second-generation (Gen-2) Polaris membrane and advanced planar membrane modules to a final form factor optimized for commercial use; validate their performance in an engineering-scale field test at TCM; and to show the potential of the MTR process to meet DOE CO2 capture targets from large source point emitters. Work for this project included membrane optimization and scale-up, advanced planar module design and fabrication, design and fabrication of an engineering-scale field test membrane skid, operation of the field test skid processing Residue Fluid Catalytic Cracker (RFCC) industrial flue gas at TCM, and a detailed techno-economic analysis (TEA) of the MTR membrane post-combustion process for CO2 capture. This project validated recent membrane technology advancements at the engineering-scale, moves the MTR advanced post-combustion capture technology to TRL-6, and mitigates risk in future Large Pilot or Demonstration scale-up activities.

20 FOSSIL-FUELED POWER PLANTS↗

Soft Robotic Glove System for Wearable Haptic Teleoperation - 18373

Managing nuclear waste is a difficult, dangerous, and costly task. Teleoperation of robotic operators is an important functionality for this application, allowing robotic devices to operate in dangerous environments while still leveraging the existing abilities and technical knowledge of human workers. One difficulty in implementing effective teleoperation systems is that human operators rely on more senses beside sight: in particular, touch. The sense of touch could provide significant benefits for teleoperation, preventing environmental mishaps that could cost millions of dollars. To help remedy this problem, we present a soft haptic feedback glove system. The primary functionality of this system is to provide tactile and force feedback on the user's hands. This was achieved using custom cylindrical pneumatic actuators. These actuators consist of plastic sheeting folded into a hollow-walled cylinder which is connected to a pressure source. When pressurized, the cylinder walls expand and the cylinder straightens, with the force being related to the pressure. The glove system has one of these actuators over each of the user's three main fingers (index, middle, and thumb), where pressurizing the actuator applies a reaction force to the finger of the user, mimicking the resistance of a grasped object. The expansion of these actuators also causes them to physically press against the user's fingers, giving them the feeling of contact. To demonstrate the usability of this system, we implemented a virtual environment in V-REP simulating the interactions between a robotic Jaco Arm and a floating ball. This simulation receives hand and figure positions from the haptic glove and used them to drive the state of the simulated Jaco Arm. In response to the interactions between the Jaco and the floating ball, the simulation returns commands to the glove, instructing it to pressurize or depressurize its actuators as appropriate, allowing the user to feel the ball in their hand. We tested users in their ability to prevent the ball from floating away using this system with and without the haptic feedback. Results show that the addition of haptic feedback significantly improved the ability of users to interact precisely with the environment. We also used the finger positions of the haptic glove to control a physical robotic hand, and demonstrated the system's usability for control of real robots, though the lack of sensors on the robot made providing the tactile feedback impractical. However, there is more to grasping and holding an object than just the feeling of contact and the resistance to closing one's hand. The weight of the object will apply forces to the user's other joints as they hold it. To this end, we added a haptic feedback device for the user's wrist, a technique that can be generalized to the user's other joints. This is made up of reverse Pneumatic Artificial Muscles (rPAMs) stretching across the wrist between mounts on the hand and forearm. These rPAM actuators are tubes of silicone wrapped in thread such that they extend when pressurized. This extension can apply torque on the user's wrist, mimicking the feeling of gravity, or allowing the simulation of additional contact forces (such as the user's hand impacting an object). We tested this wrist device independently using a simulated path-following experiment. Users showed a significant improvement in performance when aided by haptic feedback during sinusoidal path experiments, demonstrating the usability of this approach. We then integrated the two systems, further demonstrating the usability of this approach. This project represents the first step towards a more effective robotic teleoperation system, allowing skilled human workers to perform tasks more effectively while remaining out of danger. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Different singularities in the functions of extended kinetic theory at the origin of the yield stress in granular flows

We use previous results from discrete element simulations of simple shear flows of rigid, identical spheres in the collisional regime to show that the volume fraction-dependence of the stresses is singular at the shear rigidity. Here, we identify the shear rigidity, which is a decreasing function of the interparticle friction, as the maximum volume fraction beyond which a random collisional assembly of grains cannot be sheared without developing force chains that span the entire domain. In the framework of extended kinetic theory, i.e., kinetic theory that accounts for the decreasing in the collisional dissipation due to the breaking of molecular chaos at volume fractions larger than 0.49, we also show that the volume fraction-dependence of the correlation length (measure of the velocity correlation) is singular at random close packing, independent of the interparticle friction. The difference in the singularities ensures that the ratio of the shear stress to the pressure at shear rigidity is different from zero even in the case of frictionless spheres: we identify that with the yield stress ratio of granular materials, and we show that the theoretical predictions, once the different singularities are inserted into the functions of extended kinetic theory, are in excellent agreement with the results of numerical simulations.

42 ENGINEERING↗

Porflow Implementation of Vadose Zone Conceptual Model for Slit and Engineered Trenches in the E-Area Low Level Waste Facility Performance Assessment

A summary of the conceptual models to be used in representing slit (ST) and engineered trenches (ET) in PORFLOW (ACRi 2018) simulations as part of the performance assessment (PA) of the E-Area Low Level Waste Facility (ELLWF) is presented. Key details that are discussed include: model geometries, spatially dependent hydro-stratigraphic representations, model dimensionality, and boundary conditions (i.e., infiltration rates, cover overhangs, subsidence). Accounting for differences in the percent of non-crushable materials, eighteen unique models, defined by seven hydro-stratigraphic groupings, will be used to represent slit and engineered trenches.

42 ENGINEERING↗

Design, Flow Field, and Heat Transfer Characterization of the Conjugate Aero-Thermal Test Facility at NETL

A new aerothermal test facility was constructed for the purpose of studying film cooling performance in an environment that accurately simulates conjugate heat transfer characteristics that exist in engine operation. Here, this paper details the design of the facility and the plan for conducting steady state film cooling experiments to improve the understanding of conjugate heat transfer scaling from laboratory to engine conditions. The test facility consists of two separate low channels (hot gas/coolant) and each gas path has a flow conditioning section, a convergent nozzle and a test section/channel with viewports. Numerical simulations were conducted to predict flow field characteristics supporting the design of the flow loop facility. Preliminary experiments were conducted to characterize the flow field using velocity and temperature profile measurements. In addition, infrared (IR) thermography methods were developed to measure surface temperatures on the hot side of the test plate. The IR measurement methods including calibration of the IR camera is explained in detail. It was concluded that appropriate hot gas path flow conditioning could be achieved using a strainer-like tube, a perforated plate, and a honeycomb-mesh screen system upstream of the test section. Flow field measurements from preliminary experiments showed that the boundary layer profile follows the law of the wall.

47 OTHER INSTRUMENTATION↗

SU-F-18C-08: A Validation Study of a Commercially Available Software Package's Absorbed Dose Estimates in a Physical Phantom

Purpose: This study assesses the accuracy of the absorbed dose estimates from CT scans generated by Monte Carlo (MC) simulation using a commercially available radiation dose monitoring software program. Methods: Axial CT studies of an anthropomorphic abdomen phantom with dose bores at a central location and 4 peripheral locations were conducted using a fixed tube current at 120 kV. A 100 mm ion chamber and a 0.6 cc ion chamber calibrated at diagnostic energy levels were used to measure dose in the phantom at each of the 5 dose bore locations. Simulations using the software program's Monte Carlo engine were run using a mathematical model of the anthropomorphic phantom to determine conversion coefficients between the CTDIvol used for the study and the dose at the location of the dose bores. Simulations were conducted using both the software's generic CT beam model and a refined model generated using HVL and bow tie filter profile measurements made on the scanner used for the study. Results: Monte Carlo simulations completed using the generalized beam model differed from the measured conversion factors by an absolute value average of 13.0% and 13.8% for the 100 mm and 0.6 cc ion chamber studies, respectively. The MC simulations using the scanner specific beam model generated conversion coefficients that differed from the CTDIvol to measured dose conversion coefficients by an absolute value average of 7.3% and 7.8% for the 100 mm and 0.6 cc ion chamber cases, respectively. Conclusion: A scanner specific beam model used in MC simulations generates more accurate dose conversion coefficients in an anthropomorphic phantom than those generated with a generalized beam model. Agreement between measured conversion coefficients and simulated values were less than 20% for all positions using the universal beam model.

07 ISOTOPE AND RADIATION SOURCES↗