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,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Computational Fluid Dynamics Simulation of Gasoline Compression Ignition

Gasoline compression ignition (GCI) is a low temperature combustion (LTC) concept that has been gaining increasing interest over the recent years owing to its potential to achieve diesel-like thermal efficiencies with significantly reduced engine-out nitrogen oxides (NOx) and soot emissions compared to diesel engines. In this work, closed-cycle computational fluid dynamics (CFD) simulations are performed of this combustion mode using a sector mesh in an effort to understand effects of model settings on simulation results. One goal of this work is to provide recommendations for grid resolution, combustion model, chemical kinetic mechanism, and turbulence model to accurately capture experimental combustion characteristics. Grid resolutions ranging from 0.7 mm to 0.1 mm minimum cell sizes were evaluated in conjunction with both Reynolds averaged Navier-Stokes (RANS) and Large Eddy Simulation (LES) based turbulence models. Solution of chemical kinetics using the multi-zone approach is evaluated against the detailed approach of solving chemistry in every cell. The relatively small primary reference fuel (PRF) mechanism (48 species) used in this study is also evaluated against a larger 312-species gasoline mechanism. Based on these studies the following model settings are chosen keeping in mind both accuracy and computation costs – 0.175 mm minimum cell size grid, RANS turbulence model, 48-species PRF mechanism, and multi-zone chemistry solution with bin limits of 5 K in temperature and 0.05 in equivalence ratio. With these settings, the performance of the CFD model is evaluated against experimental results corresponding to a low load start of injection (SOI) timing sweep. The model is then exercised to investigate the effect of SOI on combustion phasing with constant intake valve closing (IVC) conditions and fueling over a range of SOI timings to isolate the impact of SOI on charge preparation and ignition. Simulation results indicate that there is an optimum SOI timing, in this case -30?aTDC (after top dead center), which results in the most stable combustion. Advancing injection with respect to this point leads to significant fuel mass burning in the colder squish region, leading to retarded phasing and ultimately misfire for SOI timings earlier than -42?aTDC. On the other hand, retarding injection beyond this optimum timing results in reduced residence time available for gasoline ignition kinetics, and also leads to retarded phasing, with misfire at SOI timings later than -15?aTDC.

Energy conversion↗
DOE OSTI2020

Driven-dissipative quantum mechanics on a lattice: Simulating a fermionic reservoir on a quantum computer

The driven-dissipative many-body problem remains one of the most challenging unsolved problems in quantum mechanics. The advent of quantum computers may provide a unique platform for efficiently simulating such driven-dissipative systems. But, there are many choices for how one can engineer the reservoir. One can simply employ ancilla qubits to act as a reservoir and then digitally simulate them via algorithmic cooling. A more attractive approach, which allows one to simulate an infinite reservoir, is to integrate out the bath degrees of freedom and describe the driven-dissipative system via a master equation, that can also be simulated on a quantum computer. In this work, we consider the particular case of noninteracting electrons on a lattice driven by an electric field and coupled to a fermionic thermostat. Then, we provide two different quantum circuits: the first one reconstructs the full dynamics of the system using Trotter steps, while the second one dissipatively prepares the final nonequilibrium steady state in a single step. We run both circuits on the IBM quantum experience. For circuit (i), we achieved up to five Trotter steps. When partial resets become available on quantum computers, we expect that the maximum simulation time can be significantly increased. Lastly, the methods developed here suggest generalizations that can be applied to simulating interacting driven-dissipative systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Determination of the Most Likely Wear Zones of Pipe Bends in the Pulsed Pneumatic Conveyor at the RT-1 Plant

The results of an experimental investigation and numerical modeling of the most likely high-wear zones in 90° pipeline bends in the pulsed pneumatic conveyor system at the RT-1 plant of the Mayak Industrial Association are presented. The experimental studies were conducted on a model of a pulsed pneumatic conveyor system; the numerical modeling of the gas-dynamic flow in the pipeline channel was performed in a model of a two-dimensional non-stationary motion of a compressible medium. It was found on the basis of the experimental results obtained using the model and numerical modeling in the homogeneous approximation that local regions with high wear of walls in the 90° pipe bends in the pulsed pneumatic conveyor do indeed exist and their experimentally found positions agree with the calculations.

42 ENGINEERING↗
DOE OSTI2019

Modeling of the Heat-Affected and Thermomechanically Affected Zones in a Ti-6Al-4V Inertia Friction Weld

Inertia friction welding has been used across the aerospace, automotive, and power-generation industries for the fabrication of complex axisymmetric components for over forty years. The process involves one axisymmetric piece being held stationary and another piece being brought into contact set to rotate about its axis of symmetry by a flywheel with the system under an applied load across the joint. Plasticization at the joint interface through the frictional heating sees the two pieces bond together. The titanium alloy Ti-6Al-4V has been widely studied for inertia welding applications. A successful selection of processing parameters (flywheel energy and mass, applied load) allows an inertia welding process which produces a very high-integrity weld, with a minimal heat-affected zone (HAZ) and thermomechanically affected zone (TMAZ), formed as a narrow band at the interface and extending further into the material. The width of this narrow band of heated material is dependent upon the process parameters used. A series of experimental inertia friction welds were performed using Ti-6Al-4V, and a finite element (FE) modeling framework was developed using the FE code Deform in order to predict the widths of the HAZ and TMAZ at the weld interface. The experimentally observed HAZ boundaries were correlated with the thermal fields from the FE model, while TMAZ boundaries were correlated with the Von Mises plastic strain fields.

42 ENGINEERING↗
DOE OSTI2018

Blowing polar skyrmion bubbles in oxide superlattices

Particle-like topological structures such as skyrmions and vortices have garnered ever-increasing interests due to their rich physical insights and potential broad applications in spintronics. Here, in this work, we discover the reversible switching between polar skyrmion bubbles and ordered vortex arrays in ferroelectric superlattices under an electric field, reminiscent of the Plateau-Raleigh instability in fluid mechanics. An electric field phase diagram is constructed, showing a wide stability window for the observed polar skyrmions. A “volcano”-like pontryagin density distribution is formed, indicating the formation of a smooth circular skyrmion. The topological charge Q at different applied field is calculated, verifying the field-driven topological transition between Q = 0 and Q = ±1 states. This study is a demonstration for the computational design of field-induced topological phase transitions, giving promise for the design of next-generation nanoelectronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Evaluation of enzymatic depolymerization of PET, PTT, and PBT polyesters

Millions of tons of waste polyester plastics, including polyethylene terephthalate (PET), polytrimethylene terephthalate (PTT), and polybutylene terephthalate (PBT), end up in the environment as soil and water contaminants. Recent advances in enzymatic polyester degradation have motivated researchers toward biorecycling of plastic wastes. Leaf-branch compost cutinase (LCC) enzymes have been proven to be effective in the biodegradation of PET. This study focuses on enzymatic depolymerization of waste PET, PTT, and PBT materials by using an ICCG variant of LCC (LCC ICCG ) produced from Escherichia coli BL21(DE3). The degradation efficiency of the polyesters was determined by the monomer terephthalic acid (TPA) released from the depolymerization reaction. It was found that the most efficient depolymerization was achieved for PET, followed by PTT and PBT. A kinetic model based on Langmuir adsorption and the Michaelis-Menten equation was developed to describe the enzymatic depolymerization of PET, PTT, and PBT with various enzyme and substrate loadings. The model simulation results revealed that the LCCICCG enzyme loading should be linearly increased as the work capacity of the polyester substrate increases. A specific enzyme loading of 0.91 mg/g PET is suggested to achieve 90% depolymerization of PET within three days. Furthermore, the experimental data and model simulation results can be used to help further engineer the enzyme and process to achieve a complete biodegradation of polyester wastes at large scale.

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

Modeling hydrodynamic and biomass pyrolysis effects of recycled product gases in a bubbling fluidized bed reactor

Fast pyrolysis of biomass in a fluidized bed reactor is typically conducted in a nitrogen gas environment. Recycling product gas can improve the economics of operating such a system by reducing reliance on pure process streams, but much less is known about how recycling pyrolysis product gas may affect fluidization behavior and pyrolysis kinetics. Therefore, gas effects in a fluidized bed biomass pyrolysis reactor were investigated using engineering correlations, low-order models, and CFD simulations for N 2 , H 2 , CO, CO 2 , and CH 4 carrier gas mixtures. Here, our findings reveal viscosity of a gas mixture can be significantly underestimated depending on the model and correlation. Furthermore, fluidization characteristics such as U mf and gas-solid convective heat transfer can be greatly affected by the gas properties. By utilizing H 2 as the fluidizing gas (instead of N 2 ), while maintaining a constant fluidization ratio (U s /U mf ), the bio-oil yields can be increased ~5%. This is due to the lower density H 2 producing similar hydrodynamics as N 2 at higher gas flow rates. These higher flow rates result in shorter gas residence times, and as a result, less secondary reactions that convert bio-oil to light gases and char. Model results also suggest that bio-oil yield is not significantly affected by the type of carrier gas used, with bio-oil yield varying by ~2% across different gas mixtures while maintaining constant flow rate. In conclusion, this indicates that recycled pyrolytic gas can be used as the carrier gas for biomass pyrolysis.

09 BIOMASS FUELS↗
DOE OSTI2024

A Finite Difference informed Random Walk solver for simulating radiation defect evolution in polycrystalline structures with strongly inhomogeneous diffusivity

Diffusivity of species and defects on grain boundaries is usually several orders of magnitude larger than that inside grains. Such strongly inhomogeneous diffusivity requires prohibitively high computational demands for modeling microstructural evolution. Here, this paper presents a highly-efficient numerical solver, combining the Finite Difference method and Random Walk model, designed for accurately modeling strongly inhomogeneous diffusion within polycrystalline structures. The proposed solver, termed Finite Difference informed Random Walk (FDiRW), integrates a customized Finite Difference (cFD) scheme tailored for fast diffusion along thin grain boundaries represented by a single-layer of nodes. Numerical experiments demonstrate that the FDiRW solver achieves an impressive efficiency gain of 1560x compared to traditional Finite Difference methods while maintaining accuracy, making it feasible for personal computer machines to handle diffusional systems with strongly inhomogeneous diffusivity across static polycrystalline microstructures. The model has been successfully applied to simulate radiation defect evolution, showcasing its scalability to engineering scales in both length and time dimensions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

NRAP-open-IAM: A flexible open-source integrated-assessment-model for geologic carbon storage risk assessment and management

Large-scale implementation of geologic carbon storage (GCS) to help reduce atmospheric greenhouse gas emissions requires stakeholder confidence that injected CO2 will remain contained and that potential subsurface environmental risks are acceptably small and manageable. The U.S. Department of Energy’s National Risk Assessment Partnership (NRAP) has developed an open-source integrated assessment model (NRAP-Open-IAM) to help address questions about a potential GCS site’s ability to effectively contain injected CO 2 and protect groundwater and other overlying environmentally sensitive receptors. NRAP-Open-IAM allows a user to: (1) incorporate relevant site geologic and injection scenario data; (2) characterize important site features and events;(3) couple fast prediction models of various system components of the engineered geologic system; and (4) execute stochastic, dynamic simulation of whole GCS system performance, leakage risk assessment, and uncertainty quantification. NRAP-Open-IAM is available on GitLab (https://gitlab.com/NRAP/OpenIAM), and is accompanied by multiple application examples and detailed user and developer guides.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Topological effects on separation of alkane isomers in metal–organic frameworks

Polymorphism in metal–organic frameworks (MOFs) means that the same chemical building blocks (nodes and linkers) can be used to construct isomeric MOFs with different topological networks. The choice of topology can substantially impact the pore network of the MOF, changing the sizes and shapes of the pores, which has implications for adsorption and separation applications. Here, we look at the influence of topology in 38 polymorphic MOFs on the separation of linear and branched C4–C6 alkane isomers, a separation of great importance to the petrochemical industry. We find that the MOF Cu 2 (1,4-benzenedicarboxylate) in nbo topology (nbo-Cu 2 BDC) has particularly high affinity for linear alkanes due to its small pore size, which excludes the branched isomers. Upon studying this MOF in further detail, we find that it can take either of two conformations: a cubic conformation, which is typical of nbo MOFs, and a unique star conformation that contains 1D triangular and hexagonal channels. The determination of which conformation the MOF will adopt depends on steric effects between the nodes and linkers.

42 ENGINEERING↗
DOE OSTI2013

Multi-GPU Implementation of a 3D Finite Difference Time Domain Earthquake Code on Heterogeneous Supercomputers

We have developed a highly scalable 3D Finite Difference GPU code for use in earthquake engineering and disaster management through regional petascale earthquake simulations. This MPI-CUDA code is based on a widely-used wave propagation code called AWP-ODC and restructured for high throughput and efficiency on a heterogeneous computing architecture. We present an effective communication reduction technique for leveraging GPUs with minimal PCI-e overhead, and a novel overlapping method to fully hide data communication latency between GPUs. The optimization concept used in this work can be extended to general stencil computing on a structured grid. The benchmarks demonstrated sustained 100 TFlops in single precision for 49 billion mesh points using 952 GPUs on the NCCS Titan Phase 5 system, which is a 77-fold speedup compared to the CPU version of the code. This multi-GPU implementation has been validated and used for a large-scale verification wave propagation simulation of Mw5.4 Chino Hills earthquake using 128 GPUs.

58 GEOSCIENCES↗
DOE OSTI2021

Effect of electrode rinse solutions on the electrodialysis of concentrated salts

The data presented indicate that modest attention to the chemistry of the electrode rinse solution can yield as much as 56% improvement in process efficiency. When operating electrodialysis (ED) with salt concentrations greater than 10,000 mg/l (as NaCl) the amperage can be relatively high even at low voltages. To maximize the overall rate of ion transport from diluate to concentrate, there is a need to minimize resistances in the electrode cells since these can represent as much as 30% of the total resistances in the entire process and in this particular ED unit, represented resistance roughly equivalent to that within the membrane stack. This work describes the performance of a 10-cell pair, 200 cm 2 (0.02 m 2 ) per membrane, pilot ED unit operated in batch mode at 5 V potential. The feedstock to the stack was varied from 0.5% to 6% NaCl. Results from Volt-amp profiles (2.0–15 V) were used as the rationale for choosing the preferred electrode rinse solution, 90 g/l (kg/m 3 ) disodium sulfate at pH 12.5. The recommended solution of 30 g/l (kg/m 3 ) disodium sulfate at neutral pH was tested against stronger solutions (60, 90, and 120 g/l (kg/m 3 ) disodium sulfate) at neutral pH and with 1 g/l (kg/m 3 ) sodium hydroxide added to yield solutions at pH 12.5. There were stark differences in the performance of the ED unit between the solutions at pH 7 and those at pH 12.5. The slopes of the Volt-amp profiles improved indicating less resistance to ion flow at pH 12.5. Most of the difference can be attributed to the greater conductivity of sodium hydroxide compared to disodium sulfate, where even a modest addition of 1 g/l (kg/m 3 ) sodium hydroxide increased the conductivity of the solution between 8 and 20% depending on the concentration of disodium sulfate. Most notably, the initiation voltage was 2.39 ± 0.08 V with solutions at pH 7 and 1.99 V ± 0.04 V when the rinse solutions were around pH 12.5, a difference of around 0.42 V. Since the full batch runs were performed at 5.0 V, this 0.42 V differential represented a loss of 8% in process efficiency. A series of tests showed a shift in initiation voltage occurred between pH 11.5 and 12.5. Nernst equations for water electrolysis are coupled with a simple flux model that estimates the concentration of hydroxide and hydronium ions at the electrode surfaces. This model indicates that the potential demanded by the electrodes is directly linked to the current. Furthermore, there is a steep drop in the potential demanded by the anode when there is sufficient hydroxide available to generate oxygen from hydroxide as opposed to generation from water. The voltage demand is minimal at pH 12.5 and reflective of the reference standard voltage (just over 1.23 V predicted at 0.1 amps). However, the flux model predicts that the voltage demanded by the electrodes at pH 7 is 1.61 V at 0.1 amps (the lowest amperage measurable with the power pack). In conclusion, comparing this result to the reference standard (E 0 = -1.23 V) yields a satisfactory explanation for the observed difference below reference standard of 0.42 V.

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

Comparison of Energy Efficiency and Power Density in Pressure Retarded Osmosis and Reverse Electrodialysis

Pressure retarded osmosis (PRO) and reverse electrodialysis (RED) are emerging membrane-based technologies that can convert chemical energy in salinity gradients to useful work. The two processes have intrinsically different working principles: controlled mixing in PRO is achieved by water permeation across salt-rejecting membranes, whereas RED is driven by ion flux across charged membranes. This study compares the energy efficiency and power density performance of PRO and RED with simulated technologically available membranes for natural, anthropogenic, and engineered salinity gradients (seawater-river water, desalination brine-wastewater, and synthetic hypersaline solutions, respectively). The analysis shows that PRO can achieve both greater efficiencies (54-56%) and higher power densities (2.4-38 W/m(2)) than RED (18-38% and 0.77-1.2 W/m(2)). The superior efficiency is attributed to the ability of PRO membranes to more effectively utilize the salinity difference to drive water permeation and better suppress the detrimental leakage of salts. On the other hand, the low conductivity of currently available ion exchange membranes impedes RED ion flux and, thus, constrains the power density. Both technologies exhibit a trade-off between efficiency and power density: employing more permeable but less selective membranes can enhance the power density, but undesired entropy production due to uncontrolled mixing increases and some efficiency is sacrificed. When the concentration difference is increased (i.e., natural -> anthropogenic -> engineered salinity gradients), PRO osmotic pressure difference rises proportionally but not so for RED Nernst potential, which has logarithmic dependence on the solution concentration. Because of this inherently different characteristic, RED is unable to take advantage of larger salinity gradients, whereas PRO power density is considerably enhanced. Additionally, high solution concentrations suppress the Donnan exclusion effect of the charged RED membranes, severely reducing the permselectivity and diminishing the energy conversion efficiency. This study indicates that PRO is more suitable to extract energy from a range of salinity gradients, while significant advancements in ion exchange membranes are likely necessary for RED to be competitive with PRO.

Yip, NY↗
DOE OSTI2024

Optimal control of large quantum systems: assessing memory and runtime performance of GRAPE

Abstract Gradient Ascent Pulse Engineering (GRAPE) is a popular technique in quantum optimal control, and can be combined with automatic differentiation (AD) to facilitate on-the-fly evaluation of cost-function gradients. We illustrate that the convenience of AD comes at a significant memory cost due to the cumulative storage of a large number of states and propagators. For quantum systems of increasing Hilbert space size, this imposes a significant bottleneck. We revisit the strategy of hard-coding gradients in a scheme that fully avoids propagator storage and significantly reduces memory requirements. Separately, we present improvements to numerical state propagation to enhance runtime performance. We benchmark runtime and memory usage and compare this approach to AD-based implementations, with a focus on pushing towards larger Hilbert space sizes. The results confirm that the AD-free approach facilitates the application of optimal control for large quantum systems which would otherwise be difficult to tackle.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Geometrical model for non-zero {theta}{sub 13}

Based on Friedberg and Lee's geometric picture by which the tribimaximal Pontecorvo-Maki-Nakawaga-Sakata leptonic mixing matrix is constructed, namely, corresponding mixing angles correspond to the geometric angles among the sides of a cube. We suggest that the three realistic mixing angles, which slightly deviate from the values determined for the cube, are due to a viable deformation from the perfectly cubic shape. Taking the best-fitted results of {theta}{sub 12} and {theta}{sub 23} as inputs, we determine the central value of sin{sup 2}2{theta}{sub 13} should be 0.0238, with a relatively large error tolerance; this value lies in the range of measurement precision of the Daya Bay experiment and is consistent with recent results from the T2K Collaboration.

42 ENGINEERING↗
DOE OSTI2019

Realization of ordered magnetic skyrmions in thin films at ambient conditions

Magnetic skyrmions have captivated physicists due to their topological nature and novel physical properties. In addition, skyrmions hold significant promise for future information technologies. A key barrier to realizing skyrmion-based devices has been stabilizing these spin structures under ambient conditions. In this paper, we demonstrate that the tunable magnetic properties of amorphous Fe/Gd mulitlayers enable the formation of skyrmion lattices which are stable over a large temperature and magnetic field parameter space, including room temperature and zero magnetic field. These skyrmions, having a hybrid nature displaying both Bloch-type and Néel-type characteristics, are stabilized by dipolar interactions rather than Dzyaloshinskii-Moriya interactions, typically considered a requirement for the generation of skyrmions. Small angle neutron scattering (SANS) was used in combination with soft x-ray microscopy to provide a unique, multiscale probe of the local and long-range order of these structures. The hexagonal lattice seen in SANS results from the hybrid skyrmion picture obtained with micromagnetic simulations. Furthermore, these results identify a pathway to engineer controllable skyrmion phases in thin film geometries which are stable at ambient conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Survey of Open-Source Tools for Transmission and Distribution Systems Research

This work presents a review of open-source electric power transmission and distribution systems analysis tools suitable for use by industry professionals and academic researchers. Due to the high complexity of the electric grid, there exist numerous tools and extensive research pertaining to nearly every aspect of the design, operation, and control of transmission and distribution networks. In addition to the commercial tools, a wide range of free, open-source tools, models, and data usable by the scientific community for related research have been developed by different organizations, including both international and US universities and national laboratories. However, due to the absence of a catalog of available tools, models and data, researchers often lack a knowledge of existing capabilities and may develop duplicative software and tools. Increasing awareness of these available resources seeks to accelerate their broader use, leading to more efficient and standardized grid analysis. This review paper (which is part of a larger survey effort that studied over 400 tools in the transmission, distribution, buildings, and electric vehicles space) outlines selected open-source resources that have been developed in power transmission and distribution systems research. It is anticipated that this work can serve as a guide for industry and academic researchers alike, ensuring that research efforts are well-channeled.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2018

Data-Driven Molecular Engineering of Solar-Powered Windows

Buildings are the centerpiece of modern living, with more than half of the world’s population now living in urban environments. This demographic evolution has led to building use be-coming the main drain of our energy resources. According to the US Energy Information Administration, in 2016, 40 percent of the total energy consumption in the United States came from building use. However, we could overcome this energy drain by embedding new environmental technologies into future cities to realize energy-sustainable buildings.

42 ENGINEERING↗