Aero-Thermal Characterization of Accelerating and Diffusing Passages Downstream of Rotating Detonation Combustors
Presentation from ASME Turbo Expo 2021
Braun, James↗
SEARCH · Engineering Papers
Results for “Flow cell”
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.
Aero-Thermal Characterization of Accelerating and Diffusing Passages Downstream of Rotating Detonation Combustors
Cycle benefits of rotating detonation engines show up to five percentage points of efficiency gain for low-pressure ratio engines. An optimal integration between the combustor and the turbine needs to be guaranteed to realize this potential gain. The rotating detonation combustor (RDC) exhausts transonic flow with shocks rotating at frequencies ranging from a few to tens of kilohertz depending on the number of present waves. Hence, the turbine design requires precise knowledge of the fluctuations and losses downstream of the combustor. This paper focuses on the quantification of fluctuations and losses for accelerating and diffusing passages. The analysis of the combustor is performed via reactive unsteady Reynolds Averaged Navier-Stokes (URANS) simulations. The unsteady RANS equations are solved via CFD++ from Metacomp with a one-step reaction mechanism for an H2-air mixture. The resolving of the boundary layer is achieved with a structured mesh of around 36 million cells. Inlet pressure of 10 bar and two different back pressures are applied to the combustor to model the interconnection with downstream turbines. Finally, we present and assess a methodology to reduce the computational time to model these passages ten times.
Braun, James↗
The importance of delayed neutron precursors in gamma dose calculations for activated primary heat exchanger components in the Molten Salt Breeder Reactor
The Molten Salt Reactor (MSR) Multiphysics Applications technical area in the Nuclear Energy Advanced Modeling and Simulation program has supported the development of 3D Monte Carlo models of the Molten Salt Breeder Reactor (MSBR) over the last couple of years. This MSBR model was previously run with the Shift Monte Carlo code to perform radiation shielding calculations in the reactor cell area. The MSBR is a 2250 MWth (1000 MWe) liquid-fueled molten salt reactor design developed at Oak Ridge National Laboratory in the 1970s. Determining the source terms from activated primary heat exchanger (HX) components is important because delayed and prompt neutron fluxes incident on these components affect the dose rate in the primary HX maintenance areas. This information can be used in the development of remote handling procedures required during shutdown for maintenance. A methodology has been developed and is proposed in this paper to quantify the activated source term from the primary HX components as a result of the movement of the delayed neutron precursors in flowing primary fuel salt through the primary HXs in the MSBR. The goal of this research is to evaluate the gamma dose rates in the maintenance hatches above the primary HX using the activated HX source terms. The study showed that the gamma dose rates are approximately two orders of magnitudes higher when accounting for the neutron activation from the movement of delayed neutron precursors through the HXs than when flowing fuel is not considered. Thus, the movement of delayed neutron precursors must be taken into account for accurately predicting the neutron activation of primary loop components.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
A nanoporous capacitive electrochemical ratchet for continuous ion separations
Directed ion transport in liquid electrolyte solutions underlies many phenomena in natural and industrial settings. While nature has evolved structures that drive continuous ion flow without Faradaic redox reactions, establishing this process in synthetic systems has been challenging. Here, in this study, we report an ion pump that drives aqueous ions against a force using a capacitive ratchet mechanism independent of redox reactions. Modulation of an electric potential between thin metallic layers on either face of a nanoporous alumina wafer immersed in solution results in persistent voltages and ionic currents. This occurs due to the nonlinear capacitive nature of electric double layers, whose repeated charging and discharging sustains a continuous ion flux. Using this approach, we demonstrate ratchet-driven electrodialysis that reaches a 50% decrease in the conductivity of the solution in a dilution cell. These ratchet-based ion pumps can enable continuous desalination and selective ion separation using an electrically powered device with no moving parts.
Kautz, Rylan [University of California, Irvine, CA↗
Modular Integrated System for Carbon-Neutral Methanol Synthesis Using Direct Air Capture and Carbon-Free Hydrogen Production
This study investigates the development and economic analysis of a modular integrated system for carbon-neutral methanol synthesis, leveraging direct air capture (DAC) and solid oxide electrolysis cells (SOEC) for carbon dioxide and hydrogen production, respectively. The proposed system integrates a novel building-based DAC process, functionalized solid sorbents, and low-energy SOEC technology, aiming to minimize operational and capital costs. A comparison between the base case system (1,000 t methanol/year) and a scaled-up model (14,758 t methanol/year) reveals significant improvements in efficiency and economic feasibility. The scaled-up system achieves a levelized cost of methanol (LCOM) of $740/t, a 7.5% reduction compared to that of conventional DAC-based systems, while utilizing existing building HVAC infrastructure for air handling. Detailed sensitivity analyses were conducted, evaluating the effects of plant capacity and air flow rate on the LCOM, demonstrating the scalability of the building-based DAC system. The cradle-to-gate life cycle analysis shows that the proposed process using renewable-sourced electricity achieves a 38% reduction in greenhouse gas (GHG) emission compared to reported values of green methanol production technologies that use a conventional DAC and a conventional methanol synthesis catalyst. When fossil-sourced electricity is used in the proposed process, it leads to about a 37.5% reduction in GHG emission in comparison to reported values for conventional methanol production technologies using steam methane reforming technology and fossil-sourced electricity.
alcohols↗
Autonomous Hydrogen Fueling Station
This project “Autonomous Hydrogen Fueling Station” covered the autonomous refueling with both gaseous hydrogen and liquid hydrogen. The part on gaseous hydrogen focused on the development of an autonomous robotic fueling arm that would couple to a fuel cell engine for hydrogen refueling without guidance from the forklift operator and budget period. Research was also covered for the robotic fueling with a commercial vehicle. The second phase of the project created the baseline for an autonomous liquid hydrogen transfer system that would minimize boil off losses by operating at thermodynamically efficient state points. For the development of the robotic fueling arm, testing was conducted to establish a baseline measurement of the accuracy and repeatability of a human operator positioning a lift truck in front of a dispenser. The goal was to establish the range of motion required for an autonomous fueling mechanism to mate a hydrogen nozzle with a receptacle on a fuel cell system installed in a forklift. The final design comprised a selective compliance articulated robot arm (SCARA)-type mechanism with two arms for horizontal motion and a ball screw for vertical movement and color and LIDAR cameras were used for marker identification and proximity awareness to guide the robotic arm to its target receptacle. Initial tests resulted in 199 out of 200 successful attempts at autonomous coupling of the dispensing coupler and a fuel cell engine, without hydrogen. The dispenser prototype was modified to include tubing for both hydrogen fuel and air purge lines, but subsequent tests were confounded by the shoulder motor over current errors which limited the robot from getting to the fully inserted position to achieve a positive latch. Robotic hydrogen refueling was successfully demonstrated over 1.5 hours of testing, Plug completed 29 successful latches with an average number of 4 sequential latches before failure. However, a robot capable of placing the nozzle with more force is required for higher reliability. For budget period two, a small scale (10 kg / transfer) automated control system was designed that would operate valves to control pressure and flow of liquid nitrogen between a source and receiving tank with an aim to minimize boil off losses by operating at the most thermodynamically efficient state points. Control system logic flow and a P&ID were developed prior to system safety characterization via HAZOP. A control narrative and system state points were defined. Delays in approval for a change of project objective and procurement issues precluded the construction and test of the final prototype system.
08 HYDROGEN↗
A real-time distributed solid oxide electrolysis cell (SOEC) model for cyber-physical simulation
System integration and dynamic operability between SOEC and balance-of-plant (BoP) components are major technical challenges before realizing rapid load following of SOEC systems. Cyber-physical simulation (CPS) is a leading-edge digital engineering approach and is regarded as the next step beyond Digital Twins. CPS approach can be used to research SOEC system integration and develop dynamic controls prior to actual pilot testing without using a real SOEC. To seamlessly couple with BoP hardware and access non-observable operational parameters (e.g., local temperature gradient) during transients, a distributed one-dimensional (1D) real-time SOEC model was developed. Its real-time execution was demonstrated for 20 to 640 nodes at the fixed time step of 5 ms. A higher excess air ratio enabled smaller local temperature gradients on SOEC solid materials and faster transients upon current density step change from 0.15 to 0.55 A cm -2 . During the transients, the magnitude of the peak temperature gradient nearly doubled in 10 s from -3.5 to -5.9 °C cm -1 . This represents a significant operating risk that can impact the dynamic operability of SOEC systems. In addition, the local temperature gradient was found to change directions on all nodes in SOEC solid materials, with the greatest impact on the upstream nodes. The SOEC model was also tested at the thermal neutral voltage using actual process air flow parameters as variable model inputs. Variable process air temperatures were found to induce alternating local temperature gradients on SOEC solid materials. These are new operational mechanisms for SOEC degradation relevant for load following operational modes yet distinct from previous reports. To mitigate these unfavorable features, the SOEC can be operated at voltages that are slightly (±20 mV) deviated from the thermal neutral voltage. Here, the corresponding net thermal energy change was less than 1.6% of the electric power consumption. This 1D real-time SOEC model established the basis of cyber-physical simulation of SOEC hybrid systems.
24 POWER TRANSMISSION AND DISTRIBUTION↗
Numerical Simulations of High Prandtl Number Liquid Jets Impinging on a Flat Plate
In this work, 3D simulations of oil jets impinging on a flat, heated wall are presented. The numerical setup uses the Volume of Fluid (VoF) method to model the two-phase flow. A careful grid definition across the liquid film, along with the use of the Conjugate Heat Transfer (CHT) approach allowed local heat transfer to be solved with fine resolution at the wall. Variations of liquid flow rate, liquid temperature and surface temperature allow to cover a wide range of local Reynolds and Prandtl numbers (226 < Re < 2850, 77 < Pr < 161). Resulting surface-averaged heat transfer compares very well with experimental measurements conducted in a previous study. In-depth analysis of the flow has identified expected features from the literature. In particular, the impact of jet axial velocity profiles on the heat transfer distribution in the stagnation zone was clearly stated. The increase in heat transfer when warming the liquid film was also reproduced and explained by a decrease in oil viscosity and an increase in film velocity. All those effects were taken into account in correlations for stagnation and local values of Nusselt number. A grid sensitivity study was also conducted, showing that if the grid solving the thermal boundary layer in the stagnation zone can be coarsened without impacting local and surface-averaged predictions of heat transfer, a minimum resolution (2 to 3 cells) within the thermal boundary layer is however required for an accurate prediction of heat transfer.
conjugate heat transfer↗
Strong‐Yet‐Ductile Eutectic Alloys Employing Cocoon‐Like Nanometer‐Sized Dislocation Cells
Abstract Eutectic alloys (EAs) with superior fluidity are known to be the easiest to cast into high‐quality ingots, making them the alloys of choice for making large‐sized structural parts. However, conventional EAs (CEAs) have never reached strength–ductility combinations on par with the best in other alloy categories. Via thermomechanical processing of cast Ni‐32.88wt%Fe‐9.53wt%Al CEAs, a cocoon‐like nano‐meshed (as fine as 26 nm) network of dislocations (CNN‐D) is produced via recovery annealing, through the rearrangement of cold‐work‐accumulated dislocations anchored by dense pre‐existing nanoprecipitates. In lieu of traditional plasticity mechanisms, such as TWIP and TRIP, the CNN‐D is particularly effective in eutectic lamellae with alternating phases, as it instigates nanometer‐spaced planar slip bands that not only dynamically refine the microstructure but also transmit from the FCC (face‐centered‐cubic) layers into the otherwise brittle B2 layers. These additional mechanisms for strengthening and strain hardening sustain stable tensile flow, resulting in a striking elevation of both strength and ductility to outrank not only all previous CEAs, but also the state of the art—additively manufactured eutectic high‐entropy alloys. The CNN‐D thus adds a novel microstructural strategy for performance enhancement, especially for compositionally complex alloys that increasingly make use of nanoprecipitates or local chemical order.
Chemistry↗
Mesh and model requirements for capturing deep-stall aerodynamics in low-Mach-number flows
Here, this work presents a comprehensive computational fluid dynamics investigation of the effects of grid resolution and turbulence-model choice for capturing the unsteady three-dimensional aerodynamic performance of NACA 0012 and 0021 airfoils, with specific focus on the deep-stall regime. At high angles of attack (α), wind turbine blades routinely experience vortex-induced vibrations, which can cause significant structural damages. Accurate predictions of post-stall aerodynamics can identify the frequencies at which such vibrations maybe triggered. In this context, the NACA 0012 airfoil simulations are conducted at a chord-based Reynolds number, Re c =2×10 6 , with the k-ω Shear-Stress Transport Reynolds-Averaged Navier-Stokes (RANS) and Improved Delayed Detached Eddy Simulation (IDDES) hybrid RANS-Large Eddy Simulation turbulence models. The effect of mesh resolution both in the wall-normal and spanwise directions is investigated. Only the IDDES model with a minimum spanwise resolution of 24 cells per chord length correctly predicts the aerodynamic forces. Spectral analysis shows the peak primary shedding frequency at α=30°, which signifies the end of the stall region. In the post-stall regime, both lift and drag frequencies drop asymptotically with increasing α. The Strouhal number, based on normalised chord length, remains nearly constant in this region. Based on this study, NACA 0021 airfoil runs are performed with IDDES for Re c =2.7×10 5 and 2.0×10 6 on the finest wall-normal mesh and three spanwise grids. Simulations conducted on the finer spanwise grids demonstrate grid independence and show good agreement with experiments. The effect of varying Rec on the airfoil frequency statistics is investigated. Additionally, comparison studies are presented to investigate the impact of airfoil thickness on the frequency content at Re c =2.0×10 6 . The results from the study provide guidance on the choice of mesh resolution with the IDDES model to accurately capture aerodynamic quantities for complex industrial applications.
17 WIND ENERGY↗
Imaginary admittance and charge transfer resistance correlate to the physiological status of Shewanella oneidensis cultures in real time
Monitoring microbial activity is essential for industrial and environmental applications to proceed efficiently. To minimize time and labor-intensive monitoring, a new paradigm is required for in-situ, real time analysis. Since bioconversion of organics is accomplished by microorganisms through the oxidation of feedstocks linked to the reduction of electron acceptors, microorganisms can be viewed as electrochemical catalysts. In this respect, cell membranes have an electrical potential, which is analogous to a conventional capacitor and linked dynamically to cellular activity. Here, in this study, we demonstrate the use of electrochemical impedance spectrometry (EIS) and cyclic voltammetry (CV) for monitoring microbial metabolic activity in real time, in-situ. The effect of organic electron donors as a function of concentration to the physiological status of strains of Shewanella oneidensis was determined. In this study, the pyomelanin overproducer (S. oneidensis ΔhmgA) and the pyomelanin deficient mutant (S. oneidensis ΔmelA) were chosen due to different surface electrochemical characteristics along with differences in oxygen utilization efficiency. CV, relative admittance, phase shift and permittivity changed with growth status and correlated with electron flow from organic carbon sources and terminal electron acceptor availability. This work offers a novel and inexpensive approach to real time monitoring with the advantage of abundant data.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Interpenetration and kinetic effects in converging, high-energy plasma jets
In this work, we report on numerical simulations of laser-driven convergent plasma fusion targets. These “inverted corona” fusion targets are useful for the study of counter-streaming and converging rarefied plasma flows, and previous experiments have demonstrated their potential as neutron sources. The scheme consists of a fuel layer lined along the interior surface of a hollow plastic shell that is laser-ablated and expands inward towards the target center. The plasma streams generated in these targets are initially nearly collisionless as they converge, leading to wide interaction length scales and long interaction time scales as the jets interpenetrate. Such kinetic effects impact mixing of constituent ions - a phenomenon not properly captured by single-fluid hydrodynamic simulations. Here we conduct numerical simulations using two different methods: (1) single-fluid simulations in HYDRA, and (2) kinetic-ion, fluid-electron hybrid particle-in-cell (PIC) simulations in the code Chicago. It is shown that the initially nearly collisionless plasma fronts interpenetrate deeply and lead to broader interaction regions in space and time resulting in significant beam-beam fusion. The two approaches make different, testable predictions for the effect of fuel-layer thickness on neutron yield.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Multi-Sourced Collaboration for the Production and Refining of Rare Elements and Critical Metals (Final Technical Report)
The project objective was to develop a feasible and cost-effective method for recovering rare earth elements (REEs) and critical materials (CMs) from coal and coal byproducts, resulting in high-purity individually separated REEs and CMs. The targeted REEs included Y, Pr, Nd, Gd, Dy, and Sm, with a purity of over 99.5%, while the CMs included Co, Mn, Ga, Sr, Li, Ni, Zn, and Ge, with a purity of over 90%. The project aimed to design a prototype facility capable of producing 1-3 tonnes/day of high-purity REO mixes. The work was divided into four designated circuits: 1) REE extraction and concentration, 2) REE separation and purification, 3) RE metal production, and 4) CM production. To achieve these goals, the project involved 11 tasks, including technology reviews, research, process flow diagram development, mass balance estimation, and preliminary technical-economic analysis. The project team included researchers from the University of Kentucky, University of Alabama and Virginia Tech as well as process specialists from Argonne National Laboratory. MP Materials provided technical support regarding rare earth markets and processing while Alliance Coal performed resource assessment. The project included a market analysis for Nd/Pr, Tb, Dy, Gd, Y, Co, Mn, Li, Sr, Ga, Ni, Zn, and Ge. These analyses provided insights into the supply and demand trends as well as historic and future projections of market price relative to purity requirements for these elements. Two coal resources were selected for the project: the West Kentucky No. 13 (Baker) Seam and an undisclosed lignite resource in the Illinois coal basin. The estimated quantities of REEs in these resources were calculated based on production samples and drilling data. It was estimated that there is adequate supply for an operation producing one metric ton daily of higher purity mixed rare earth oxides (MREO) for approximately 20 years at a site located in western Kentucky. In Circuit 1, project data was obtained from a pilot heap leach and REE concentration facility. It was concluded that the existing circuit, which generated a MREO concentrate, two types of CM mixed products, and Li- and Sr-containing waters, would be suitable feed for circuits 2-4. Data from the first-of-its-kind coal coarse refuse heap leach pilot pad played a crucial role in estimating reliable elemental concentrations of the pregnant leaching solution (PLS). The average total REE concentration in the PLS was found to be 28.6 ppm. In Circuit 2, several concepts were explored including a novel process referred to as solvent-assisted chromatography (SAC). This concept involved a novel columnar reactor that incorporated multiple mixer/settlers, thereby enabling the operation of counter-flowing aqueous and organic phases. Unfortunately, due to project time constraints, a complete fundamental modeling analysis could not be completed to fully evaluate the technology. Molten salt electrowinning was considered as an alternative for circuit 3 following circuit 2 purification circuit utilizing the novel SAC process. A mass and energy balance of Nd reduction to metal in a fluoride containing molten salt electrolyte was conducted. Comparisons were made with the current state of Asian molten salt electrorefining, and potential improvements in siphoning rare earth metals (REM) from the reactor were presented. A cost estimate was performed for the production of 1 tonne per day, which yielded a total of $2.29 million for the nine electrowinning (EW) cells required. The selected option for circuits 2 and 3 was a plasma distillation process, which initially separates rare earth elements (REEs) from other elements. This is followed by selective electrowinning in various ionic liquids. The selection was made on the basis of thermodynamic modeling and experimental data previously published by a project partner. The combination offers an innovative approach to integrated refining and RE metal production. For Circuit 4, an extensive literature review was conducted for the processing of the CMs. The ultimate decision was to utilize a combined plasma and ionic liquid process as well to produce individual high-purity concentrates of Zn, Ni, Co, Mn, and Mg. A separate flowsheet for Li and Sr was recommended, which would yield carbonates of these elements. Due to the lack of suitable experimental data at this time, a process recommendation could not be provided but several methods have been proposed for consideration. Lastly, a techno-economic analysis (TEA) was conducted to assess the effectiveness of the proposed process for further investigation. The TEA results revealed a capital expense (CapEx) of $737 million and an annual operational expense (OpEx) of $220 million. Due to the selected elements, the hypothetical heap leach pad can produce 1 metric tonne per day of REO equivalent, but a conscious decision was made to only treat targeted REEs, resulting in the production of 0.4 metric tonne of REM. An estimated annual revenue of $90.87 million was projected based on standard market pricing information provided by the funding agency. During the TEA, ten different modules were evaluated for costing purposes. The precipitation circuit was identified as the largest single operational expense, followed by the Mg/Mn process due to the amount of treated metal. In terms of capital expenditures, the heap leach process incurred the highest cost, followed by the Mg/Mn process. The scalability of the plasma process is a crucial consideration since the reactors cannot be scaled beyond the largest demonstrated size due to their reliance on surface area of the slag and vapor phase. The purity estimate for the REEs are generally 98%±2% to produce a metal. The purity level being lower than the project objective was due to the lack of specific experimental data needed to tighten the tolerance of the estimates. Based on literature and previous experience, the CMs are estimated as follows; Ga (95%+, metal), Sr (95%+, carbonate), Li (95%+, carbonate), Ni (98%±2%, metal), Zn (95%+, metal sponge), Ge (95%+, metal), Co (98%±2%, metal), and Mn (98%±2%, metal).
01 COAL, LIGNITE, AND PEAT↗
Atmospheric Pressure Plasma Deposition of TiO2: A Review
Atmospheric pressure plasma (APP) deposition techniques are useful today because of their simplicity and their time and cost savings, particularly for growth of oxide films. Among the oxide materials, titanium dioxide (TiO2) has a wide range of applications in electronics, solar cells, and photocatalysis, which has made it an extremely popular research topic for decades. Here, we provide an overview of non-thermal APP deposition techniques for TiO2 thin film, some historical background, and some very recent findings and developments. First, we define non-thermal plasma, and then we describe the advantages of APP deposition. In addition, we explain the importance of TiO2 and then describe briefly the three deposition techniques used to date. We also compare the structural, electronic, and optical properties of TiO2 films deposited by different APP methods. Lastly, we examine the status of current research related to the effects of such deposition parameters as plasma power, feed gas, bias voltage, gas flow rate, and substrate temperature on the deposition rate, crystal phase, and other film properties. The examples given cover the most common APP deposition techniques for TiO2 growth to understand their advantages for specific applications. In addition, we discuss the important challenges that APP deposition is facing in this rapidly growing field.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Enabling Combustion Science Simulations for Future Exascale Machines
Reacting flow simulations for combustion applications require extensive computing capabilities. Leveraging the AMReX library, the Pele suite of combustion simulation tools targets the largest supercomputers available and future exascale machines. We introduce PeleC, the compressible solver in the Pele suite, and detail its capabilities, including complex geometry representation, chemistry integration, and discretization. We present a comparison of development efforts using both OpenACC and AMReX's C++ performance portability framework for execution on multiple GPU architectures. We discuss relevant details that have allowed PeleC to achieve high performance and scalability. PeleC's performance characteristics are measured through relevant simulations on multiple supercomputers. The success of PeleC's design for exascale is exhibited through demonstration of a 160 billion cell simulation and weak scaling onto 100\% of Summit, an NVIDIA-based GPU supercomputer at Oak Ridge National Laboratory. Our results provide confidence that PeleC will enable future combustion science simulations with unprecedented fidelity.
combustion↗
Risk Assessment and Ventilation Modeling for Hydrogen Release in Vehicle Repair Garages
The availability of repair garage infrastructure for hydrogen fuel cell vehicles is becoming increasingly important for future industry growth. Ventilation requirements for hydrogen fuel cell vehicles can affect both retrofitted and purpose-built repair garages and the costs associated with these requirements can be significant. A hazard and operability study (HAZOP) was performed to identify key risk-significant scenarios related to hydrogen vehicles in a repair garage. Detailed simulations and modeling were performed using appropriate computational tools to estimate the location, behavior, and severity of hydrogen release based on key HAZOP scenarios. This work compares current fire code requirements to an alternate ventilation strategy to further reduce potentially hazardous conditions. Overall, the amount of flammable mass of hydrogen at any one time in the simulation is low compared to the total mass of hydrogen released, due to the low flow rate of a low pressure release. It is shown that position, direction, and velocity of ventilation have a significant impact on the amount of instantaneous flammable mass in the domain.
33 ADVANCED PROPULSION SYSTEMS↗
COMPUTATIONAL MODELING OF IGNITION AND PREMIXED FLAME PROPAGATION INITIATED BY A PRE-CHAMBER TURBULENT JET
Addressing the pressing need for reduced carbon emissions, Turbulent Jet Ignition (TJI) emerges as a promising technology for ultra-lean combustion, offering enhanced thermal efficiencies and minimized cyclic variability in spark-ignited engines. To facilitate rapid testing and integration of this technology, a robust computational modeling framework is crucial. This study delves into the predictive capabilities of computational models for main-chamber ignition and premixed flame propagation using a single-cycle TJI rig measured by Biswas et al. (Applied Thermal Engineering, vol 106, 2016). Employing an open-source compressible flow simulation solver with Large Eddy Simulation (LES) for turbulence modeling, the investigation integrates the conventional Laminar Finite Rate Chemistry (LFRC) model alongside the transported Probability Density Method (PDF) for turbulence-chemistry interaction. A fully-consistent Eulerian Monte-Carlo Fields (EMCF) method is utilized to approximate the transported PDF, while Interaction by Exchange with Mean is employed to close micro-mixing terms in stochastic differential equations. A reduced chemical reaction mechanism with 21 species and 84 reactions (DRM-19) is used for solving chemical kinetics, and a double Gaussian energy deposition model is used to approximate the spark ignition in the pre-chamber. An unstructured O-grid mesh with 0.3 million cells in the prechamber and 1 million cells in the main chamber is employed. Results are divided into two phases: pre-chamber initialization and full TJI simulations. Validation of the predicted pre-chamber flame propagation and the lean ignition in the main-chamber is carried out by using available experimental data. Under quiescent conditions, both the LFRC and transported PDF methods largely underestimate the flame speed and subsequent pressure growth in the pre-chamber. A linear momentum forcing technique is applied to investigate the impact of initial turbulence in the pre-chamber, demonstrating a notable influence on flame propagation. Fine-tuning of the forcing coefficient reproduces the sudden pressure growth observed in the experiment. The experimentally validated pre-chamber simulation serves as the initial condition for the full TJI simulations. It is found that the LFRC model fails to predict lean-ignition in the main-chamber, resulting in a misfiring event. Incorporation of turbulence-chemistry interaction using the transported PDF method substantially improves the prediction of the ignition event in the main-chamber, achieving fair qualitative agreement and quantitative validation of combustion parameters within ±10% of the reported experimental data. The rich simulation results consisting of a full set of statistical description of the thermo-chemical states enable us to gain deep insights into the ignition mechanisms in the main chamber, which is limited when done experimentally. A novel dual ignition phenomenon is revealed in the TJI rig for the first time. Initially, a primary ignition kernel is formed at a downstream location which eventually detaches from the main jet. As the jet momentum decreases, a secondary ignition event follows, this time at a more upstream location which eventually combines with the primary ignition kernel to form a single connected flame front. Investigation of these ignition sequences in chemical composition space reveal distinct differences between the two. The primary ignition event in the main-chamber is followed by a large concentration of active radicals from the pre-chamber jet, accelerating the chain-branching steps, characterizing what has been referred to as flame ignition. In contrast, the secondary ignition occurs in the absence of active radicals in the pre-chamber jet, hence characterized as jet ignition. Further analysis of the effect of pre-chamber jet characteristics on lean ignition in the main-chamber is conducted by setting up cases with different initial pressure ratios (por) between the two chambers, a non-dimensional parameter, ranging from 1.2 to 3.2. As the initial pressure ratio increases, jet momentum increases, with dual ignition observed in cases above por= 2.2. Case with por= 3.2 lead to misfiring. The effect of ignition sequence on global combustion characteristics of TJI is analyzed. Dual ignition events lead to non-monotonicity in combustion characteristics such as global reaction progress variable, flame penetration, and global heat release rate. In dual ignition events, although the rate of fuel consumption and global heat release rate is initially lower, the secondary ignition leads to a sudden increase in flame surface area, resulting in a sudden jump and promoting the overall performance of the TJI system.
42 ENGINEERING↗
Asynchronous distributed-memory task-parallel algorithm for compressible flows on unstructured 3D Eulerian grids
Here, we discuss the implementation of a finite element method, used to numerically solve the Euler equations of compressible flows, using an asynchronous runtime system (RTS). The algorithm is implemented for distributed-memory machines, using stationary unstructured 3D meshes, combining data-, and task-parallelism on top of the Charm++ RTS. Charm++’s execution model is asynchronous by default, allowing arbitrary overlap of computation and communication. Task-parallelism allows scheduling parts of an algorithm independently of, or dependent on, each other. Built-in automatic load balancing enables continuous redistribution of computational load by migration of work units based on real-time CPU load measurement. The RTS also features automatic checkpointing, fault tolerance, resilience against hardware failure, and supports power-, and energy-aware computation. We demonstrate scalability up to 25 x 10 9 cells at $\mathscr{O}$10 4 compute cores and the benefits of automatic load balancing for irregular workloads. The full source code with documentation is available at https://quinoacomputing.org.
42 ENGINEERING↗