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 361 records · Page 20

Advanced Drying Cycle Simulator Thermal Response to Commercial Drying Conditions

The purpose of this report is to document updates on testing of the apparatus built to simulate commercial drying procedures for spent nuclear fuel at the Nuclear Energy Work Complex at Sandia National Laboratories. Validation of the extent of water removal in a dry spent nuclear fuel storage system based on drying procedures used at nuclear power plants is needed to close existing technical gaps. Operational conditions leading to incomplete drying may have potential impacts on the fuel, cladding, and other components in the system during subsequent storage and disposal. A general lack of data suitable for model validation of commercial nuclear canister drying processes necessitates well-designed investigations of drying process efficacy and water retention that incorporate relevant physics and well-controlled boundary conditions. This report documents testing updates for the Advanced Drying Cycle Simulator (ADCS). This apparatus was built to simulate commercial drying procedures and quantify the amount of residual water remaining in a pressurized water reactor (PWR) fuel assembly after drying. The ADCS was constructed with a prototypic 17×17 PWR fuel skeleton and waterproof heater rods to simulate decay heat. These waterproof heaters are the next generation design to heater rods developed and tested at Sandia National Laboratories in FY20. This report describes preliminary testing of the ADCS through measurement and analysis of the thermal response of the system to a subset of commercial drying conditions that exclude the introduction of water, namely simulated decay heats and pressures relevant to commercial drying. This test series, referred to as a “dry” test series in this report, spans three uniform waterproof heater rod powers (representing spent fuel decay heats), four helium fill pressures, and six vacuum levels. This test series was conducted to cover the range of expected ADCS testing conditions for upcoming “wet” testing, where water will be introduced and a simulated commercial drying cycle will be performed. The dry test conditions were derived from the commercial drying conditions seen in the High Burnup Demonstration and the vacuum drying conditions chosen for a smaller scale Dashpot Drying Apparatus tested at Sandia National Laboratories in FY22. For a given uniform power and pressure/vacuum level, the ADCS was operated at constant power and pressure and allowed to reach steady state conditions. The thermal data obtained from these tests were analyzed, and the results can inform computational models built to simulate commercial drying processes by providing baseline thermal data prior to the introduction of water. Following the preliminary dry tests, a test plan for the ADCS will be developed to implement a drying procedure that begins with the introduction of water to the system and is based on measurements from the drying process used for the High Burnup Demonstration Project. While applying power to the simulated fuel rods, this procedure is expected to consist of filling the ADCS vessel with water, draining the water with applied pressure and multiple helium blowdowns, evacuating additional water with a vacuum drying sequence at successively lower pressures, and backfilling the vessel with helium. Additional investigations are expected to feature failed fuel rod simulators with engineered cladding defects and guide tubes with obstructed dashpots to challenge the drying system with multiple water retention sites. The data from these investigations is expected to inform the efficacy of commercial drying operations through the quantification of residual water in a prototypic-length dry storage canister.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Increasing Al-Pair Abundance in SSZ-13 Zeolite via Zeolite Synthesis in the Presence of Alkaline Earth Metal Hydroxide Produces Hydrothermally Stable Co-, Cu- and Pd-SSZ-13 Materials

Replacing alkaline for alkaline-earth metal hydroxide in the synthesis gel during the synthesis of siliceous SSZ-13 zeolite (Si/Al~10) yields SSZ-13 with novel, advantageous properties. Its NH 4 -form ion-exchanges higher amount of isolated divalent M(II) ions than the conventional one: this is the consequence of an increased number of Al pairs in the structure induced by the +2 charge of Sr(II) cations in the synthesis gel that force two charge-compensating AlO 4 – motives to reside closer together. We characterize the +2 state of Co(II) ions in these materials with infra-red spectroscopy and X-ray absorption spectroscopy measurements and show their utility for NOx pollutant adsorption from ambient air: the ones derived from SSZ-13 with higher Al pair content contain more isolated cobalt(II) and, thus, perform better as ambient-air NOx adsorbers. Notably, Co(II)/SSZ-13 with an increased number of Al pairs is significantly more hydrothermally stable than its NaOH-derived analogue. Loading Pd(II) into Co-SSZ-13(Sr) produces an active NOx adsorber (PNA) material that can be used for NOx adsorption from simulated diesel engine exhaust. The critical issue for these applications is hydrothermal stability of Pd-zeolites. Pd/SSZ-13 synthesized in the presence of Sr(OH) 2 does not lose its PNA capacity after extremely harsh aging at 850 and 900 °C (10 h in 10% H 2 O/air flow) and loses only ~55% capacity after hydrothermal aging at 930 °C. This can be extended to other divalent metals for catalytic applications, such as copper: we show that Cu/SSZ-13 catalyst can survive hydrothermal aging at 920 °C without losing its catalytic properties, metal dispersion and crystalline structure. Thus, we provide a new, simple, and scalable strategy for making remarkably (hydro)thermally stable metal-zeolite materials/catalysts with a number of useful applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multi-dimensional modeling of non-equilibrium plasma for automotive applications

While spark-ignition (SI) engine technology is aggressively moving towards challenging (dilute and boosted) combustion regimes, advanced ignition technologies generating non-equilibrium types of plasma are being considered by the automotive industry as a potential replacement for the conventional spark-plug technology. However, there are currently no models that can describe the low-temperature plasma (LTP) ignition process in computational fluid dynamics (CFD) codes that are typically used in the multi-dimensional engine modeling community. A key question for the engine modelers that are trying to describe the non-equilibrium ignition physics concerns the plasma characteristics. A key challenge is also represented by the plasma formation timescale (nanoseconds) that can hardly be resolved within a full engine cycle simulation. This paper reports on the multi-dimensional modeling of LTP generated by a nanopulsed high-voltage discharge in a pin-to-pin electrode configuration, and evaluates the effects of ambient pressure on the post-discharge results. It is shown that a nanopulsed delivery can result in a LTP or a mode transition to an arc-like event depending on the mixture properties in the gap between the two electrodes. It is further shown that numerical predictions of the LTP-to-arc transition point as a function of ambient pressure along with qualitative distributions of important radicals and temperature closely match recent experiments with identical geometry and discharge characteristics. As a result, the high-fidelity modeling effort described in this paper can be used as a key tool to deliver in-depth understanding of non-equilibrium plasmas for automotive applications and lead to the future implementation of dedicated models for ignition technologies that are based on non-equilibrium plasma physics.

Scarcelli, Riccardo↗

Thin Silicon Solar Cells: A Path to 35% Shockley-Queisser Limits

Crystalline silicon technology is expected to remain the leading photovoltaic industry workhorse for decades. We present here the objectives and workplan of a recently launched project funded by the U.S. Department of Energy through the Foundational Program to Advance Cell Efficiency II (FPACE II), which aims at leading crystalline silicon to an efficiency breakthrough. The project will tackle fundamental approach of materials design, defect engineering, device simulations and materials growth and characterization. Among the main novelties, the implementation of carrier selective contacts made of wide bandgap material or stack of materials is investigated for improved passivation, carrier extraction and carrier transport. Based on an initial selection of candidate materials, preliminary experiments are conducted to verify the suitability of their critical parameters as well as preservation of the silicon substrate surface and bulk properties. The target materials include III-V and metal-oxide materials.

Ding, Laura↗

Hiearchical Engine for Large Scale Infrastructure Simulation

HELICS ls a new open-source, cyber-physlcal-energy co-simulation framework for electric power systems. HELICS Is designed to support very-large-scale (100,000+ federates) co­simulations with off-the-shelf power-system, communication, market, and end-use tools. Other key features Include cross platform operating system support, the integration of both eventdrlven (e.g., packetlzed communication) and time-series (e.g.,power flow) simulations, and the ability to co-Iterate among federates to ensure physical model convergence at each time step.

Top, Philip↗

Predicting Cycle-to-cycle Variations in a Spark-ignition Engine using Multi-cycle Large Eddy Simulation

Multi-cycle large eddy simulation (LES) of a spark-ignition engine has been performed to predict cycle-to-cycle variations (CCV). To describe spark-ignited premixed flame propagation, the front propagation formulation (FPF) method is combined with a laminar-to-turbulent flame transition model that describes the non-equilibrium flame development processes during the initial stage of the kernel growth. In the laminar-to-turbulent flame transition model, the effects of turbulence on the variations in the laminar-to-turbulent flame transition time are modeled based on the Kolmogorov energy cascade. 15 consecutive cycle LES of in-cylinder processes in a spark-ignition four-valve single cylinder engine is conducted. The predicted in-cylinder pressure is compared with the experimental measurements. The range of CCV observed in the experiment is reproduced to an acceptable degree. The predicted CCV is found to be closely related to the difference in the turbu- lence field near the spark plug around the time of spark discharging, which is consistent with the previous LES studies. The current results also suggest the importance of capturing the variations of the laminar-to-turbulent flame transition in predicting CCV.

Su, Yunde↗

Predicting Cycle-to-cycle Variations in a Spark-ignition Engine using Multi-cycle Large Eddy Simulation

Multi-cycle large eddy simulation (LES) of a spark-ignition engine has been performed to predict cycle-to-cycle variations (CCV). To describe spark-ignited premixed flame propagation, the front propagation formulation (FPF) method is combined with a laminar-to-turbulent flame transition model that describes the non-equilibrium flame development processes during the initial stage of the kernel growth. In the laminar-to-turbulent flame transition model, the effects of turbulence on the variations in the laminar-to-turbulent flame transition time are modeled based on the Kolmogorov energy cascade. 15 consecutive cycle LES of in-cylinder processes in a spark-ignition four-valve single cylinder engine is conducted. The predicted in-cylinder pressure is compared with the experimental measurements. The range of CCV observed in the experiment is reproduced to an acceptable degree. The predicted CCV is found to be closely related to the difference in the turbulence field near the spark plug around the time of spark discharging, which is consistent with the previous LES studies. The current results also suggest the importance of capturing the variations of the laminar-to-turbulent flame transition in predicting CCV.

Yu, Yunde↗

A REDUCED ORDER MODELING APPROACH TO PROBABILISTIC CREEP-DAMAGE PREDICTIONS IN FINITE ELEMENT ANALYSIS

This paper introduces a computationally efficient Reduced Order Modeling (ROM) approach for the probabilistic prediction of creep-damage failure. Component-level probabilistic simulations are needed to assess the reliability and safety of high-temperature components. Full-scale probabilistic creep-damage modeling in finite element (FE) approach is computationally expensive requiring many hundreds of simulations to replicate the uncertainty of component failure. To that end, ROM is proposed to minimize the elevated computational cost while controlling the loss of accuracy. It is proposed that full-scale probabilistic simulations can be completed in 1D at a reduced cost, the extremum conditions extracted, and those conditions applied for lower-cost 2D/3D probabilistic simulations of components that capture the mean and uncertainty of failure. The probabilistic Sine-hyperbolic (Sinh) model is selected which in previous work was calibrated to alloy 304 stainless steel. The Sinh model includes probability density functions (pdfs) for test condition (stress and temperature), initial damage (i.e. microstructure), and material properties uncertainty. The Sinh model is programmed into ANSYS finite element software using the USERCREEP.F material subroutine. First, the Sinh model and FE code are subject to verification and validation to ensure the accuracy of the simulations. Numerous Monte Carlo simulations are executed in a 1D model to generate probabilistic creep deformation, damage, and rupture data. This data is analyzed and the probabilistic parameters corresponding to extreme creep response are extracted. The ROM concept is applied where only the extreme conditions are applied in the 2D probabilistic prediction of a component. The probabilistic predictions between the 1D and 2D geometry is compared to assess ROM for creep. The accuracy of the probabilistic prediction employing the ROM approach will potentially reduce the time and cost of simulating complex engineering systems. Future studies will introduce multi-stage Sinh, stochasticity, and spatial uncertainty for improved prediction.

36 MATERIALS SCIENCE↗

Effective Lanthanide Diffusivity through U-Zr Metallic Fuel

Fuel-cladding chemical interaction (FCCI) is one of the main concerns for the performance of U-Zr metallic fuels. Although FCCI results in complex phase transformations and morphological changes in both the fuel and cladding sides of the interface, the formation of a brittle layer in the cladding side of the interface (sometimes referred to as wastage) has one of the most significant impacts on the cladding mechanical integrity. This brittle layer is associated with the formation of intermetallic compounds between the cladding constituents Fe, Cr and the lanthanide fission product species. The formation of this cladding wastage layer has been the focus of recent efforts to develop a mechanistic modeling framework, including production of lanthanides in the fuel, transport through the fuel to the fuel-cladding interface, and phase transformation to form intermetallic phases. To calculate the rate of transport of lanthanides through the fuel, a multi-scale computational approach has been used. Atomistic calculations have been performed to determine the diffusivities of Nd, the most prevalent lanthanide species, through the solid fuel matrix and along pore surfaces. These diffusivities have been used as input parameters for a mesoscale model to determine an effective diffusion coefficient that accounts for the formation of porosity, its interconnection, and infiltration with bond sodium from the fuel-cladding gap. The effective diffusivity has been used in engineering-scale simulations using the BISON fuel performance code, and BISON simulations have been compared to representative validation cases from the EBR-II reactor.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Predicting Instability and the Effect of Wind Loading on Single-Axis Trackers

As PV modules continue to trend toward larger, thinner, and more flexible forms they grow more susceptible to damage from dynamic wind loading. As a result, understanding the impact of wind on PV systems, particularly when mounted on solar-tracking hardware, and identifying robust, stable array layouts and stow strategies is becoming increasingly important for the PV community. In our ongoing DuraMAT project, we are developing an open-source software package, PVade (PV aerodynamic design engineering), to simulate the cascading fluid-structure interaction that occurs within solar-tracking arrays to enable researchers to test hardware, layout, and tracker control changes, leading to enhanced stability and a reduction in wind-driven damage. We will give an overview of the PVade software, highlighting recent user-interface and algorithm developments, before presenting the latest outcomes from our ongoing validation campaign in which we analyze and compare with experimental data obtained from a DuraMAT 1 project. From there, we will present simulated results from a larger, multi-row array and highlight relationships between varying tracker angles and stability as measured by different experimentally validated metrics.

fluid structure interaction↗

Multi-scale modeling of wastage layer formation in metallic fuel cladding

Fuel-cladding chemical interaction (FCCI) is a major concern for U-Zr metallic fuels' performance, primarily due to the formation of a brittle layer (wastage) in the cladding. This brittle layer, resulting from intermetallic compounds between cladding constituents Fe, Cr, and lanthanide fission products, significantly impacts the cladding's mechanical integrity. Recent efforts focus on developing a mechanistic modeling framework to understand lanthanide production, transport to the fuel-cladding interface, and phase transformation to intermetallic phases. A multi-scale computational approach has been used to calculate lanthanide transport rates, with atomistic calculations determining Nd diffusivities through the solid fuel matrix and along pore surfaces. These diffusivities inform a mesoscale model to determine an effective diffusion coefficient, accounting for porosity and infiltration with bond sodium. This effective diffusivity is used in engineering-scale simulations via the BISON fuel performance code, which has been validated against EBR-II and FFTF reactor experiments.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

LES FOR SIMULATING THE GAS EXCHANGE PROCESS IN A SPARK IGNITION ENGINE

The gas exchange process is known to be a significant source of cyclic variability in Internal Combustion Engines (ICE). Traditionally, Large Eddy Simulations (LES) are expected to capture these cycle-to-cycle variations. This paper reports a numerical effort to establish best practices for capturing cyclic variability with LES tools in a Transparent Combustion Chamber (TCC) spark ignition engine. The main intention is to examine the sensitivity of cycle averaged mean and Root Mean Square (RMS) flow fields and Proper Orthogonal Decomposition (POD) modes to different computational hardware, adaptive mesh refinement (AMR) and LES sub-grid scale (SGS) models, since these aspects have received little attention in the past couple of decades. This study also examines the effect of near-wall resolution on the predicted wall shear stresses. LES is pursued with commercially available CONVERGE code. Two different SGS models are tested, a one-equation eddy viscosity model and dynamic structure model. The results seem to indicate that both mean and RMS fields without any SGS model are not much different than those with LES models, either one-equation eddy viscosity or dynamic structure model. Computational hardware results in subtle quantitative differences, especially in RMS distributions. The influence of AMR on both mean and RMS fields is negligible. The predicted shear stresses near the liner walls is also found to be relatively insensitive to near-wall resolution except in the valve curtain region.

Ameen, Muhsin M↗

Coupling subsurface and above-surface models for optimizing the design of borefields and district heating and cooling systems in the presence of varying water-table depth. In: Proceedings, 46th Workshop on Geothermal Reservoir Engineering

Dynamic energy simulation is important for the design and sizing of district heating and cooling systems with geothermal heat exchange. Current modeling approaches in building and district energy simulation tools typically consider heat conduction through the ground between boreholes, without flow of groundwater. On the other hand, detailed simulation tools for subsurface heat and mass transfer exist, but these fall short in simulating above-surface energy systems. To support the design and operation of such systems, we have developed a coupled model including a software package for building and district energy simulation, and software for detailed heat and mass transfer in the ground. For the first, we use the open-source Modelica Buildings Library, which includes dynamic simulation models for building and district energy and control systems. For the heat and mass transfer in the soil, we use the TOUGH simulator. TOUGH can model heat and multi-phase, multi-component mass transport for a variety of fluid systems, as well as chemical reactions, in fractured porous media. In previous work, we described the coupling of these software packages, including how time-dependent boundary conditions for the borehole walls are synchronized for use in Modelica and TOUGH. We verified that the coupled Modelica/TOUGH code produced consistent results with the original Modelica code for an idealized problem in which heat transfer was purely by conduction in a uniform geologic medium. Here, we examine less idealized problems for which TOUGH’s advanced capabilities for modeling fluid flow are required. The first problem has a shallow vadose zone, and the second problem has a thicker vadose zone with a water-table depth that varies in time, which requires a fine vertical grid discretization for the TOUGH model.

Doughtry, Christine↗

Application of High Performance Computing for Simulating Cycle-to-Cycle Variation in Dual-Fuel Combustion Engines

Interest in operational cost reduction is driving engine manufacturers to consider lower-cost fuel substitution in heavy-duty diesel engines. These dual-fuel (DF) engines could be operated either in diesel-only mode or operated with premixed natural gas (NG) ignited by a pilot flame of compression-ignited direct-injected diesel fuel. One promising application is that of large-bore, medium-speed engines such as those used in locomotives. With realistic natural gas substitution levels in the fleet of locomotives currently in service, such fuel substitution could result in billions of dollars of savings annually in the US alone. However, under certain conditions, dual-fuel operation can result in increased cycle-to-cycle variability (CCV) during combustion, resulting in variations in cylinder pressure and work extraction. In certain situations, the CCV of dual-fuel operation can be notably higher than that of diesel-only combustion under similar operating conditions. Excessive CCV can limit the NG substitution rate and operating range of a dual-fuel engine by increasing emissions and reducing engine stability, reliability and fuel efficiency via incomplete natural-gas combustion. Running multiple engine cycles in series to simulate CCV can be quite time consuming. Hence innovative modelling techniques and large computing resources are needed to investigate the factors affecting CCV in dual-fuel engines. This paper discusses the use of the High Performance Computing resource Titan, at the Oak Ridge Leadership Computing Facility at Oak Ridge National Laboratory, to investigate cycle-to-cycle combustion variability of a dual-fuel engine. The CONVERGE CFD software was used to simulate multiple, parallel single cycles of dual-fuel combustion with perturbed operating parameters and boundary conditions. These perturbations are imposed according to a sparse grids sampling of the parameter space. The sampling scheme chosen is similar to a design of experiments method but uses functions designed to minimize the number of samples required to achieve a desired degree of accuracy. The perturbed input parameters were connected to engine performance output over a large parameter space. Additionally, the same methodology was used to perform a set of diesel-only simulations at similar operating conditions to assess the numerical capability of capturing relatively low combustion variation. Lastly, this technique is expected to be useful to understand and predict combustion stability limits (knock and misfire) which are important for designing dual-fuel engines with improved operability and reliability.

33 ADVANCED PROPULSION SYSTEMS↗

Characterisation of two-stage ignition in diesel engine-relevant thermochemical conditions using direct numerical simulation

With the goal of providing a more detailed fundamental understanding of ignition processes in diesel engines, this study reports analysis of a direct numerical simulation (DNS) database. In the DNS, a pseudo turbulent mixing layer of dimethyl ether (DME) at 400 K and air at 900 K is simulated at a pressure of 40 atmospheres. At these conditions, DME exhibits a two-stage ignition and resides within the negative temperature coefficient (NTC) regime of ignition delay times, similar to diesel fuel. The analysis reveals a complex ignition process with several novel features. Autoignition occurs as a distributed, two-stage event. The high-temperature stage of ignition establishes edge flames that have a hybrid premixed/autoignition flame structure similar to that previously observed for lifted laminar flames at similar thermochemical conditions. In conclusion, a combustion mode analysis based on key radical species illustrates the multi-stage and multi-mode nature of the ignition process and highlights the substantial modelling challenge presented by diesel combustion.

33 ADVANCED PROPULSION SYSTEMS↗

Scalable Data Management of the Uintah Simulation Framework for Next-Generation Engineering Problems with Radiation

The need to scale next-generation industrial engineering problems to the largest computational platforms presents unique challenges. This paper focuses on data management related problems faced by the Uintah simulation framework at a production scale of 260K processes. Uintah provides a highly scalable asynchronous many-task runtime system, which in this work is used for the modeling of a 1000 megawatt electric (MWe) ultra-supercritical (USC) coal boiler. At 260K processes, we faced both parallel I/O and visualization related challenges, e.g., the default file-per-process I/O approach of Uintah did not scale on Mira. In this paper we present a simple to implement, restructuring based parallel I/O technique. We impose a restructuring step that alters the distribution of data among processes. The goal is to distribute the dataset such that each process holds a larger chunk of data, which is then written to a file independently. This approach finds a middle ground between two of the most common parallel I/O schemes–file per process I/O and shared file I/O–in terms of both the total number of generated files, and the extent of communication involved during the data aggregation phase. To address scalability issues when visualizing the simulation data, we developed a lightweight renderer using OSPRay, which allows scientists to visualize the data interactively at high quality and make production movies. Finally, this work presents a highly efficient and scalable radiation model based on the sweeping method, which significantly outperforms previous approaches in Uintah, like discrete ordinates. The integrated approach allowed the USC boiler problem to run on 260K CPU cores on Mira.

97 MATHEMATICS AND COMPUTING↗