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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,063 records · Page 59

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Effect of the volume fraction gradient on the phase interaction force model for disperse two-phase flows

In this work, the effects of the particle volume fraction gradient on fluid-particle interactions are studied. The phase interaction force is decomposed into three terms. For the first term, namely the symmetrized force density, we present theoretical reasoning and numerical evidence to assume that it is independent of the particle volume fraction gradient. The second term is the particle volume fraction gradient times a newly introduced diffusion stress. The third term is the divergence of the particle-fluid-particle (PFP) stress. If this assumption of independence of the particle volume fraction gradient for the first term can be verified, to the first order of the ratio of the mean distance between particles to the macroscopic lengthscale, all three terms can be studied and modeled in flows with uniform particle distributions. Models thus obtained are applicable to statistically inhomogeneous flows, with the second and third terms accounting for statistical inhomogeneity. To verify this assumption, numerical simulations of flows passing fixed arrays of particles are performed. Both uniform and nonuniform particle volume fractions are studied and compared for disperse multiphase flows with the particle Reynolds numbers ranging from 1 to 100, and particle volume fraction ranging from 1% to 26% in statistically steady states. It is found that the symmetrized force (first) term can be well approximated by the drag force obtained from studies of uniform flows. The diffusion stress is positive along the flow direction and negative in the directions perpendicular to the flow. In the case of moving particles, this stress could potentially cause particle aggregation in the flow direction and dispersion in the directions perpendicular to the flow. Finally, the diffusion stress is only important when there is a volume fraction gradient, while the PFP stress can be important in inhomogeneous flows with either nonuniform particle concentrations or nonuniform average relative velocities between the phases.

42 ENGINEERING↗

DOE OSTI2024

A Gaussian Process-Based extended Goldak heat source model for finite element simulation of laser powder bed fusion additive manufacturing process

In this study, laser powder bed fusion (L-PBF) additive manufacturing (AM) is a key enabling technology to manufacture highly complex and integrated metallic structures. In L-PBF AM process, the melting of the metal powders and the layers underneath can be governed by either “conduction mode” or “keyhole mode”, with the keyhole mode reportedly leading to porosity and decreased strength and ductility by many studies. In part scale simulations, finite element (FE) model is often used to study the temperature distribution during printing and to predict the residual stress, where a volumetric heat flux with a Gaussian or a double ellipsoidal (Goldak) distribution is often applied as the laser heat source. However, the above heat source models can only capture the melt pool shape in the conduction mode, and fail to capture the transition to keyhole melting mode when the process parameters change. To overcome this inaccuracy, an extended Goldak heat source model is proposed by introducing a laser penetration term as a function of laser parameters obtained from a Gaussian-Process (GP) model. The model is validated by “2D pad” AlSi10Mg L-PBF experiments under a wide range of laser power, scan speed, and laser focus offset, and the results show the model successfully captures the measured melt pool shape in all conditions.

36 MATERIALS SCIENCE↗

DOE OSTI2016

Assessing pellet-clad interaction with VERA: WBN1, cycles 6-7

The Consortium for Advanced Simulation of Light Water Reactors (CASL) aims to provide high fidelity simulations of nuclear reactors through the development of the Virtual Environment for Reactor Applications (VERA). One of the primary goals of CASL is to estimate local cladding failure through pellet-clad interaction (PCI). While this capability is being pursued through several different approaches, the procedure presented here focuses on running independent fuel performance calculations with the BISON fuel performance code being developed by Idaho National Laboratory using the multicycle output data from coupled neutron transport/thermal-hydraulics simulations generated by MPACT, CTF, and MAMBA. This type of approach is consistent with traditional fuel performance analysis methods, where the fuel performance analyses are typically separate from core simulation analyses. Recent work simulating 12 cycles of Watts Bar Unit 1 (WBN1) with VERA are capitalized upon, and stand-alone quarter-core BISON results can be obtained for several cycles, assessing parameters of interest to PCI (e.g., maximum centerline fuel temperature, maximum clad hoop stress, minimum gap thickness, and cumulative damage index). This is accomplished by using the coupled neutronics/thermal-hydraulics output power and temperature distributions as inputs for each BISON fuel rod case. Building on previous work which demonstrated this capability for Cycles 1 through 3, this paper presents the results for Cycles 6 and 7. Cycle 7 is of particular interest because significant CRUD-induced power shift was observed. While the current results do not immediately indicate that any rods are at an extreme risk of failure, further analysis on rods highlighted by this work may shed additional light on the likelihood of failure. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2019

Suite of Comprehensive Rapid Analysis Tools for Estimation (SOCRATES) - 19076

Long-term monitoring and remedy implementation require ongoing assessment and analytics, much of which relies on intensive and time-consuming data integration and engineering, data manipulation, and 'big-data' analytics. Many of these analyses, such as groundwater flow direction, contaminant rate attenuation, and non-parametric statistics, are spatial-temporal in nature and involve coalescing disparate data sets from a variety of sources with measurement frequencies from milliseconds to years, as well as total data volumes often exceeding the petabyte scale. Web-based analytical tools are advantageous for ongoing real-time or near-real-time ingestion of data from various sources for general data access and visualization, and also allow for complex, data rich calculations by leveraging expansive networks and virtually unlimited compute power. A web-based framework offers a natural architecture for such endeavors, as access to compute nodes through auto-scaling and dynamic frameworks enables speedy data manipulation and unlimited access to any data set that is made available via an http request. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

Real-time Distribution Simulation and Application Development for Power Systems Education

To help bridge the gap between traditional power system engineering instruction and emerging industry needs, a new set of classroom and research tools are needed. Real-time simulation tools emulating power system control room software present an opportunity to introduce students to the array of operational considerations, technical challenges, and decision-making associated with power system operations. The GridAPPS-D platform is proposed to support coursework and academic research as it provides an open-source platform for simulation, application development, and software integration. The GridAPPS-D platform, simulation capabilities, development environment, and interactive training are discussed in the context of lessons-learned from implementation for undergraduate and graduate students in the US and India. A series of potential GridAPPS-D supported academic capabilities are introduced.

Active distribution networks, open educational res↗

DOE OSTI2023

A computational parametric study of ducted fuel injection implementation in a heavy-duty diesel engine

Experiments have shown that ducted fuel injection (DFI) effectively reduces soot emissions from direct-injection diesel engines. Although many computational studies have evaluated DFI’s spray development and soot reduction mechanisms in constant volume chambers, only limited computational work on internal combustion engines exists. The DFI duct assembly changes the engine’s in-cylinder flow, spray, and combustion development. Therefore, current production engine designs might not be optimal for achieving the best engine performance with DFI. Here, this work conducted an extensive numerical study to evaluate how parameter changes affect DFI performance. The parameters include swirl ratio, piston geometry, compression ratio (CR), number of injector orifices, split injection strategy, and exhaust gas recirculation (EGR) in a heavy-duty diesel engine utilizing DFI. The combustion and soot emission data from the Sandia compression ignition optical research engine were used for model validation. Simulations showed that an increased swirl ratio resulted in more intense jet flame-piston interaction, slowing down the combustion heat release during the late combustion stage and leading to lower indicated thermal efficiency (ITE) due to higher exhaust losses. A piston-bowl design with a reentrant inner piston edge yielded the highest thermal efficiency, due to the reduced cylinder head heat transfer loss. Additional injector orifices led to higher efficiency owing to a more advanced combustion phasing. Nevertheless, the maximum pressure rise rate (MPRR) and oxides of nitrogen (NO x ) emissions also increased with the number of injector orifices due to more rapid heat release and higher combustion temperature. Implementation of a split injection strategy combined with a higher EGR rate effectively inhibited the excessive MPRR and NO x formation. In general, the study concluded that DFI is not sensitive to most parameter changes but will benefit from future parameter optimization.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2021

Unexpected Scaling of Peak Acceleration for a Yielding Mass-Spring System Subjected to a Triangular Base Acceleration Pulse

The use of bounding scenarios is a common practice that greatly simplifies the design and qualification of structures. However, this approach implicitly assumes that the quantities of interest increase monotonically with the input to the structure, which is not necessarily true for nonlinear structures. This paper surveys the literature for observations of nonmonotonic behavior of nonlinear systems and finds such observations in both the earthquake engineering and applied mechanics literature. Numerical simulations of a single degree-of-freedom mass-spring system with an elastic–plastic spring subjected to a triangular base acceleration pulse are then presented, and it is shown that the relative acceleration of this system scales nonmonotonically with the input magnitude in some cases. The equation of motion for this system is solved symbolically and an approximate expression for the relative acceleration is developed, which qualitatively agrees with the nonmonotonic behavior seen in the numerical results. The nonmonotonicity is investigated and found to be a result of dynamics excited by the discontinuous derivative of the base acceleration pulse, the magnitude of which scales nonmonotonically with the input magnitude due to the fact that the first yield of the spring occurs earlier as the input magnitude is increased. The relevance of this finding within the context of defining bounding scenarios is discussed, and it is recommended that modeling be used to perform a survey of the full range of possible inputs prior to defining bounding scenarios.

42 ENGINEERING↗

DOE OSTI2025

Cross-scale covariance for material property prediction

A simulation can stand its ground against an experiment only if its prediction uncertainty is known. The unknown accuracy of interatomic potentials (IPs) is a major source of prediction uncertainty, severely limiting the use of large-scale classical atomistic simulations in a wide range of scientific and engineering applications. Here we explore covariance between predictions of metal plasticity, from 178 large-scale (~10 8 atoms) molecular dynamics (MD) simulations, and a variety of indicator properties computed at small-scales (≤10 2 atoms). All simulations use the same 178 IPs. In a manner similar to statistical studies in public health, we analyze correlations of strength with indicators, identify the best predictor properties, and build a cross-scale “strength-on-predictors” regression model. This model is then used to estimate regression error over the statistical pool of IPs. Small-scale predictors found to be highly covariant with strength are computed using expensive quantum-accurate calculations and used to predict flow strength, within the statistical error bounds established in our study.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Assessment of Existing Capabilities and Future Needs for Designing Networked Microgrids

This is a review of existing microgrid design tool capabilities, such as the Microgrid Design Tool (MDT), LANL PNNL NRECA Optimal Resilience Model (LPNORM), Distributed Energy Resource-Customer Adoption Model (DER-CAM), Renewable Energy Optimization (REopt), and the Hybrid Optimization Model for Multiple Energy Resources (HOMER). Additionally, other simulation and analysis tools which may provide fundamental support will be examined. These will include GridLAB-DTM, OpenDSS, and the hierarchical Engine for Large-scale Infrastructure Co-Simulation (HELICS). Their applicability to networked microgrid operations will be evaluated, and strengths and gaps of existing tools will be identified. This review will help to determine which elements of the proposed optimal design and operations (OD&D) tool should be formulated from first principles, and which elements should be integrated from past DOE investments.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2017

Verification suite for the application of the limiting surface density method to arrays of 9975 packages

Criticality Safety Engineering at the Savannah River Site's K-Area Complex seeks to reduce the computational burden of analysis for the K-Area Complex. By applying an adapted Limiting Surface Density hand calculation method, the need for lengthy detailed Monte Carlo models of the storage arrays may be reduced or eliminated. After promising initial results using models of idealized but heterogeneous packages, application of the method to large arrays of highly heterogeneous 9975 shipping packages is attempted. A SCALE KENO-VI model of the 9975 package, with true dimensions and radial and axial heterogeneity is developed. A suite of array arrangements representing both arrays with equal number of packages per side (package-cubic) and arrays with unequal number of packages per side (non-cubic) are evaluated in SCALE. The cubic arrangements are used to derive the 9975-specific material and geometry constants for the LSD method. These constants are then used to implement the LSD method and compare the LSD- predicted critical mass per package to SCALE-predicted critical mass per package. Sensitivities to certain material and geometric nuances are evaluated for inclusion in the verification suite. Finally, the ability of the method to predict neutron multiplication changes induced by fissile mass changes is examined as a segue for evaluating credible abnormal conditions with the LSD method. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2016

Simulation methodology for fuel assembly drop accident during handling

Fresh and irradiated fuel assembly handling is a routine part of Nuclear Power Plant activities. Fuel Assembly handling occurs during fresh fuel receipt, core refueling, spent fuel assemblies transfer from spent fuel pool to the storage containers, etc. Westinghouse, ENUSA and KEPCO-NF have collaborated in the technical area of Fuel Assembly Drop Accident (FADA) Simulation to update the approach for estimating the number of failed fuel rods as a consequence of a fresh or irradiated fuel assembly drop, which might occur during a handling accident in a PWR-type reactor. The proposed FADA simulation methodology starts with an assessment of fuel assembly impact velocity at the final stage of the drop in water or air. The vertical, horizontal or inclined drop directions are considered. This velocity is an input to determine a fuel assembly structural response during impact which defines the loading function for the fuel assembly components including spacer grids. The dynamic load on the fuel rod is calculated from the detailed spacer grid model and the corresponding loading function. Finally, the load from the spacer grid spring and dimples are applied to the fuel rod model to determine the cladding parameters for use in compliance to fuel rod cladding failure criteria. A combination of first principal, finite element analysis (FEA) and computational fluid dynamic (CFD) methods has been utilized to create a technically sound methodology. A number of specific and traditional tests for the fuel assembly and its components has been used to confirm that the proposed approaches which describe the physical phenomena related to the fuel assembly drop accident are reasonably conservative.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

Variable-density buoyancy-driven turbulence with asymmetric initial density distribution

The effects of different initial density distributions on the evolution of buoyancy-driven homogeneous variable-density turbulence (HVDT) at low (0.05) and high (0.75) Atwood numbers are studied by using high-resolution direct numerical simulations. HVDT aims to mimic the acceleration-driven Rayleigh–Taylor and shock-driven Richtmyer–Meshkov instabilities and reveals new physics that arise from variable-density effects on the turbulent mixing. Here, the initial amounts of pure light and pure heavy flows are altered primarily to mimic the variable-density turbulence at the different locations of the Rayleigh–Taylor and Richtmyer–Meshkov instabilities’ mixing layers where the amounts of the mixing fluids are not equal. It is found that for the low Atwood number cases, the asymmetric initial density distribution has limited effects on both global and local flow evolution for HVDT. However, at high Atwood number, both global flow evolution and the local flow structures are strongly affected by the initial composition ratio. The flow composed of more light fluid reaches higher turbulent levels and the local statistics reach their fully-developed behavior earlier in the time evolution. In contrast, during the late time decay, where most of the flow is well-mixed, all parameters become independent of the initial composition ratio for both low and high Atwood number cases.

42 ENGINEERING↗

DOE OSTI2021

Investigation of high-pressure turbulent jets using direct numerical simulation

Direct numerical simulations of free round jets at a Reynolds number (Re D ) of 5000, based on jet diameter (D) and jet-exit bulk velocity (U e ), are performed to study jet turbulence characteristics at supercritical pressures. The jet consists of nitrogen (N 2 ) that is injected into N 2 at the same temperature. To understand turbulent mixing, a passive scalar is transported with the flow at unity Schmidt number. Two sets of inflow conditions that model jets issuing from either a smooth contraction nozzle (laminar inflow) or a long pipe nozzle (turbulent inflow) are considered. By changing one parameter at a time, the simulations examine the jet-flow sensitivity to the thermodynamic condition (characterized in terms of the compressibility factor ($Z$) and the normalized isothermal compressibility), inflow condition and ambient pressure (p ∞ ) spanning perfect- to real-gas conditions. The inflow affects flow statistics in the near field (containing the potential core closure and the transition region) as well as further downstream (containing fully developed flow with self-similar statistics) at both atmospheric and supercritical p ∞ . The sensitivity to inflow is larger in the transition region, where the laminar-inflow jets exhibit dominant coherent structures that produce higher mean strain rates and higher turbulent kinetic energy than in turbulent-inflow jets. Decreasing $Z$ at a fixed supercritical p ∞ enhances pressure and density fluctuations (non-dimensionalized by local mean pressure and density, respectively), but the effect on velocity fluctuations depends also on the local flow dynamics. When $Z$ is reduced, large mean strain rates in the transition region of laminar-inflow jets significantly enhance velocity fluctuations (non-dimensionalized by local mean velocity) and scalar mixing, whereas the effects are minimal in jets from turbulent inflow.

42 ENGINEERING↗

DOE OSTI2022

Competing Easy-Axis Anisotropies Impacting Magnetic Tunnel Junction-Based Molecular Spintronics Devices (MTJMSDs)

Molecular spintronics devices (MSDs) attempt to harness molecules’ quantum state, size, and configurable attributes for application in computer devices—a quest that began more than 70 years ago. In the vast number of theoretical studies and limited experimental attempts, MSDs have been found to be suitable for application in memory devices and futuristic quantum computers. MSDs have recently also exhibited intriguing spin photovoltaic-like phenomena, signaling their potential application in cost-effective and novel solar cell technologies. The molecular spintronics field’s major challenge is the lack of mass-fabrication methods producing robust magnetic molecule connections with magnetic electrodes of different anisotropies. Another main challenge is the limitations of conventional theoretical methods for understanding experimental results and designing new devices. Magnetic tunnel junction-based molecular spintronics devices (MTJMSDs) are designed by covalently connecting paramagnetic molecules across an insulating tunneling barrier. The insulating tunneling barrier serves as a mechanical spacer between two ferromagnetic (FM) electrodes of tailorable magnetic anisotropies to allow molecules to undergo many intriguing phenomena. Our experimental studies showed that the paramagnetic molecules could produce strong antiferromagnetic coupling between two FM electrodes, leading to a dramatic large-scale impact on the magnetic electrode itself. Recently, we showed that the Monte Carlo Simulation (MCS) was effective in providing plausible insights into the observation of unusual magnetic domains based on the role of single easy-axis magnetic anisotropy. Here, we experimentally show that the response of a paramagnetic molecule is dramatically different when connected to FM electrodes of different easy-axis anisotropies. Motivated by our experimental studies, here, we report on an MCS study investigating the impact of the simultaneous presence of two easy-axis anisotropies on MTJMSD equilibrium properties. In-plane easy-axis anisotropy produced multiple magnetic phases of opposite spins. The multiple magnetic phases vanished at higher thermal energy, but the MTJMSD still maintained a higher magnetic moment because of anisotropy. The out-of-plane easy-axis anisotropy caused a dominant magnetic phase in the FM electrode rather than multiple magnetic phases. The simultaneous application of equal-magnitude in-plane and out-of-plane easy-axis anisotropies on the same electrode negated the anisotropy effect. Our experimental and MCS study provides insights for designing and understanding new spintronics-based devices.

42 ENGINEERING↗

DOE OSTI2016

Experimental study of the downstream effect of LOCA-generated debris on the pressure drop of full scale fuel assembly

The emergency core cooling system (ECCS) works to cool the reactor core in the event of a postulated loss-of-coolant accident (LOCA). The long-term reactor cooling water of ECCS comes from containment sump. Usually, the containment sump contains one or more screens to prevent the components of the ECCS from debris which could be washed into the sump. Large debris is blocked by the sump screens. However, some of the debris, such as fibrous, chemical precipitates and small particles may be ingested into the ECCS and subsequently into the reactor coolant systems (RCS), which would challenge the long-term core cooling (LTCC) and threaten the reactor core safety. NRC address this concern Generic Letter 2004-02. Westinghouse has evaluated the AP1000 1/4 height fuel assembly head-loss test in response to the issue. AVERA also perform the scaled 1/3 height mock-up fuel assembly test. The debris effect on the pressure drop of half-length mock-up PLUS7 fuel assembly of APR1400 were performed in KHNP. Numerical simulation by RELAP5-3D were performed for the debris-generated core blockage and full core as well as core bypass blockage condition were considered. However, the test in the works mentioned above use the scaled down FA (1/4, 1/3 or 1/2 length of a full FA), and didn't give the relationship between pressure drop and debris batches. The purpose of this study was to extend the debris effect to the full scale PWR fuel assembly, and identified batch by batch effect of debris to head loss of FA. The specific experimental procedure is concluded. The pressure drop test without adding debris is performed at the beginning and then three kinds of debris are added to the test loop step by step to test their effect on the pressure drop of fuel assembly. Two kinds of tests have been carried out, and the quantitative pressure data are obtained and analysis is provided. This work provides reliable and quantitative data for nuclear power plant (NPP) downstream effect analysis and safety analysis program, which is important to for licensing of NPP and further benchmark exercise. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2016

Verification and Validation of CFD Simulations involving Twin Rectangular Jets using Steady RANS in Star-CCM+

Analyses such as component design and safety analysis can be conducted through careful application of computational fluid dynamics (CFD) codes. Before these codes can be considered ready for specific analysis, benchmarking (validation) against applicable datasets is needed. The twin jet water facility (TJWF) developed originally at the University of Tennessee, Knoxville has been purposed for aiding in this type of work. The TJWF allows for large data sets to be created using state of the art non-invasive measurement techniques within a mostly transparent testing volume. Utilizing velocity data collected in the TJWF by Wang et. al., the Star-CCM+ code was benchmarked a using steady Reynolds averaged Navier Stokes (RANS) framework for a combining two jet case. Careful efforts to properly conduct verification of the numerical work was undertaken with a sample reported in this summary. The computational results were validated against the experimental work and found to be within acceptable limits. The parameters of interest were the merging and combining points (MP and CP) of the twin jets, and the velocity profiles at increasing heights after injection. Fig.1 provides a graphical depiction of the twin jets combining within a free shear volume with key geometric factors included. The MP and velocity profiles were predicted with reasonable agreement. Whereas, the combining point of the twin jets was found to be over-predicted by {approx}24% difference. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2017

Methodology of Experiment Design for Improving Transient Critical Heat Flux Model Development

One phenomenon that contributes to nuclear fuel and clad failure mechanisms is the departure from nucleate boiling (DNB). Accident scenarios such as reactivity-initiated accidents (RIAs) provide conditions where the frequency of boiling crisis is much more likely to occur and on much shorter time scales than loss of coolant (LOCA) or loss of flow (LOFA) accidents. Here, boiling crisis in the form of DNB is attributed to transient-induced heat flux values that exceed the critical heat flux (CHF). The resulting effects, such as bubble crowding and eventual dryout, can affect flow boiling heat transfer and compromise the structural integrity of nuclear reactor materials. The focus of the work proposed here lies in developing a rigorous process using out-of-pile testing, under simulated RIA conditions, to determine the transient critical heat flux of cladding materials, without contributions from aforementioned radiation-induced effects that have been observed during in-pile experiments. This study provides a road-map for developing and performing experiments that simulate the power profile and energy deposition aspects of transient events, specific to RIAs, on the laboratory scale. The goal of our project is to observe and measure DNB in an annular flow channel around a single heater rod that is operating and nominal thermal hydraulic conditions of a PWR. It is also an aim of these experiments to determine how to quantify the gap between SS and transient CHF cases, based on state-dependent properties of the cladding material type that is being interrogated. Traditionally, CHF is determined as a combination of four factors: hydraulic diameter, system pressure, mass flow, fluid quality. The authors propose that an additional measureable quantity, which is dependent on the structural properties of the cladding material itself, is partly reasonable for the difference between measured and model-developed CHF values. Using a carefully metered approach will allow for several key concepts to be realized: compensation for transient-based effects to current models for CHF; adaptation of current models for different cladding material types, including Gen-IV concept materials; and implementation of a correction factor to account for both power transient-based and radiation-induced effects on the CHF that have not yet been accounted for.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Turbulence/flame/wall interactions in non-premixed inclined slot-jet flames impinging at a wall using direct numerical simulation

In this paper, three-dimensional turbulent non-premixed oblique slot-jet flames impinging at a wall were investigated using direct numerical simulation (DNS). Two cases are considered with the Damköhler number (Da) of case A being twice that of case B. A 17 species and 73-step mechanism for methane combustion was employed in the simulations. It was found that flame extinction in case B is more prominent compared to case A. Reignition in the lower branch of combustion for case A occurs when the scalar dissipation rate relaxes, while no reignition occurs in the lower branch for case B due to excessive scalar dissipation rate. A method was proposed to identify the flame quenching edges of turbulent non-premixed flames in wall-bounded flows based on the intersections of mixture fraction and OH mass fraction iso-surfaces. The flame/wall interactions were examined in terms of the quenching distance and the wall heat flux along the quenching edges. There is essentially no flame/wall interaction in case B due to the extinction caused by excessive turbulent mixing. In contrast, significant interactions between flames and the wall are observed in case A. The quenching distance is found to be negatively correlated with wall heat flux as previously reported in turbulent premixed flames. The influence of chemical reactions and wall on flow topologies was identified. The FS/U and FC/U topologies are found near flame edges, and the NNN/U topology appears when reignition occurs. The vortex-dominant topologies, FC/U and FS/S, play an increasingly important role as the jet turbulence develops.

42 ENGINEERING↗