Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic model reduction”

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,315 records · Page 73

Integrated hydrogeophysical modelling and data assimilation for geoelectrical leak detection

Time-lapse electrical resistivity tomography (ERT) measurements provide indirect observations of hydrological processes in the Earth's shallow subsurface at high spatial and temporal resolution. ERT has been used in the past decades to detect leaks and monitor the evolution of associated contaminant plumes. Specifically, inverted resistivity images allow visualization of the dynamic changes in the structure of the plume. However, existing methods do not allow the direct estimation of leak parameters (e.g. leak rate, location, etc.) and their uncertainties. We propose an ensemble-based data assimilation framework that evaluates proposed hydrological models against observed time-lapse ERT measurements without directly inverting for the resistivities. Each proposed hydrological model is run through the parallel coupled hydro-geophysical simulation code PFLOTRAN-E4D to obtain simulated ERT measurements. The ensemble of model proposals is then updated using an iterative ensemble smoother. In this paper, we demonstrate the proposed framework on synthetic and field ERT data from controlled tracer injection experiments. Our results show that the approach allows joint identification of contaminant source location, initial release time, and solute loading from the cross-borehole time-lapse ERT data, alongside with an assessment of uncertainties in these estimates. We demonstrate a reduction in site-wide uncertainty by comparing the prior and posterior plume mass discharges at a selected image plane. This framework is particularly attractive to sites that have previously undergone extensive geological investigation (e.g., nuclear sites). It is well suited to complement ERT imaging and we discuss practical issues in its application to field problems.

58 GEOSCIENCES↗

Many-Exciton Quantum Dynamics in a Ruddlesden–Popper Tin Iodide

We present a study on the many-body exciton interactions in a Ruddlesden–Popper tin halide, namely, (PEA) 2 SnI 4 (PEA = phenylethylammonium), using coherent two-dimensional electronic spectroscopy. The optical dephasing times of the third-order polarization observed in these systems are determined by exciton many-body interactions and lattice fluctuations. We investigate the excitation-induced dephasing (EID) and observe a significant reduction of the dephasing time with increasing excitation density as compared to its lead counterpart (PEA) 2 PbI 4 , which we have previously reported in a separate publication [J. Chem. Phys. 2020, 153, 164706]. Surprisingly, we find that the EID interaction parameter is four orders of magnitude higher in (PEA) 2 SnI 4 than that in (PEA) 2 PbI 4 . This increase in the EID rate may be due to exciton localization arising from a more statically disordered lattice in the tin derivative. This is supported by the observation of multiple closely spaced exciton states and the broadening of the linewidth with increasing population time (spectral diffusion), which suggests a static disordered structure relative to the highly dynamic lead-halide. Additionally, we find that the exciton nonlinear coherent lineshape shows evidence of a biexcitonic state with low binding energy (<10 meV) not observed in the lead system. We model the lineshapes based on a stochastic scattering theory that accounts for the interaction with a nonstationary population of dark background excitations. Our study provides evidence of differences in the exciton quantum dynamics between tin- and lead-based Ruddlesden–Popper metal halides (RPMHs) and links them to the exciton–exciton interaction strength and the static disorder aspect of the crystalline structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Arctic-associated increased fluctuations of midlatitude winter temperature in the 1.5° and 2.0° warmer world

In recent decades, the interior regions of Eurasia and North America have experienced several unprecedentedly cold winters despite the global surface air temperature increases. One possible explanation of these increasing extreme cold winters comes from the so-called Warm Arctic Cold Continent (WACC) pattern, reflecting the effects of the amplified Arctic warming in driving the circulation change over surrounding continents. This study analyzed reanalysis data and model experiments forced by different levels of anthropogenic forcing. It is found that WACC exists on synoptic scales in observations, model’s historical and even future runs. In the future, the analysis suggests a continued presence of WACC but with a slightly weakened cold extreme due to the overall warming. Warm Arctic events under the warmer climate will be associated with not only a colder continent in East Asia but also a warmer continent, depending on the teleconnection process that is also complicated by the warmer Arctic. Such an increasingly association suggests a reduction in potential predictability of the midlatitude winter anomalies.

54 ENVIRONMENTAL SCIENCES↗

The kinematic properties of Milky Way stellar halo populations

ABSTRACT In the Gaia era, stellar kinematics are extensively used to study Galactic halo stellar populations, to search for halo structures, and to characterize the interface between the halo and hot disc populations. We use distribution function-based models of modern data sets with 6D phase space data to qualitatively describe a variety of kinematic spaces commonly used in the study of the Galactic halo. Furthermore, we quantitatively assess how well each kinematic space can separate radially anisotropic from isotropic halo populations. We find that scaled action space (the ‘action diamond’) is superior to other commonly used kinematic spaces at this task. We present a new, easy to implement selection criterion for members of the radially anisotropic Gaia–Enceladus merger remnant. Assuming a 1:1 ratio of Gaia–Enceladus stars to more isotropic halo, we find our selection achieves a sample purity of 86 per cent in our models with respect to contamination from the more isotropic halo. We compare this criterion to literature criteria, finding that it produces the highest purity in the resulting samples, at the expense of a modest reduction in completeness. We also show that selection biases that underlie nearly all contemporary spectroscopic data sets can noticeably impact the E−Lz distribution of samples in a manner that may be confused for real substructure. We conclude by providing recommendations for how authors should use stellar kinematics in the future to study the Galactic stellar halo.

79 ASTRONOMY AND ASTROPHYSICS↗

Applying Graph Theory to Problems in Air Traffic Management

Graph theory is used to investigate three different problems arising in air traffic management. First, using a polynomial reduction from a graph partitioning problem, it is shown that both the airspace sectorization problem and its incremental counterpart, the sector combination problem are NP-hard, in general, under several simple workload models. Second, using a polynomial time reduction from maximum independent set in graphs, it is shown that for any fixed e, the problem of finding a solution to the minimum delay scheduling problem in traffic flow management that is guaranteed to be within n1-e of the optimal, where n is the number of aircraft in the problem instance, is NP-hard. Finally, a problem arising in precision arrival scheduling is formulated and solved using graph reachability. These results demonstrate that graph theory provides a powerful framework for modeling, reasoning about, and devising algorithmic solutions to diverse problems arising in air traffic management.

precision arrival scheduling↗

Clearance Sensitivity Mitigation in Small Core Compressors

As future aircraft become lighter and more aerodynamically efficient, thrust requirements will decrease, reducing the core size of the engine. Furthermore, in the pursuit of improved fuel burn, engine overall pressure ratio and bypass ratio will increase, further driving down engine core size. These drivers together mean that the core size, or corrected mass flow rate at the compressor exit, for future single-aisle aircraft applications will shrink below 3.0 lbm/s. Traditionally, this small core compressor size is in the domain of axi-centrifugal designs, machines that are typically less efficient and limited to pressure ratios of ~25 due to stress and thermomechanical fatigue in the centrifugal impeller. In this light, NASA, Pratt & Whitney and the Raytheon Technologies Research Center embarked upon a program to develop technologies to enable an all-axial high-pressure compressor with a core size below 3.0 lbm/s and an overall pressure ratio greater than 50. One of the challenges with an all-axial high-pressure compressor at this core size is the small span at the rear of the compressor. As core size is scaled down, the rotor tip clearances do not scale with span, leading to significant efficiency penalties. This paper documents a numerical and experimental effort to mitigate this penalty through design optimization and technology insertion, enabling an estimated 5 to 10% fuel burn reduction relative to 2020 best-in-class. The experimental program consisted of two builds of a high-speed rig: a baseline build and a second build to insert technology to manage large rotor tip gap/span ratios. The results demonstrated a reduction in the sensitivity of the compressor to large clearance ratios, which would be likely at the end-of-life of a small core application. The test campaign demonstrated that, with the insertion of technology, small core compressors can maintain a similar efficiency to current best-in-class large core size compressors. In addition to advancing the state-of-the-art of technology, the program also advanced the modeling standards for multistage compressors with large clearance-to-span ratios. The validation of this approach is described in the paper.

Axial Compressors↗

Wind Tunnel Test of a Risk-Reduction Wing/Fuselage Model to Examine Juncture-Flow Phenomena

A wing/fuselage wind-tunnel model was tested in the Langley 14- by 22-foot Subsonic Wind Tunnel in preparation for a highly-instrumented Juncture Flow Experiment to be conducted in the same facility. This test, which was sponsored by the NASA Transformational Tool and Technologies Project, is part of a comprehensive set of experimental and computational research activities to develop revolutionary, physics-based aeronautics analysis and design capability. The objectives of this particular test were to examine the surface and off-body flow on a generic wing/body combination to: 1) choose a final wing for a future, highly instrumented model, 2) use the results to facilitate unsteady pressure sensor placement on the model, 3) determine the area to be surveyed with an embedded laser-doppler velocimetry (LDV) system, 4) investigate the primary juncture corner- flow separation region using particle image velocimetry (PIV) to see if the particle seeding is adequately entrained and to examine the structure in the separated region, and 5) to determine the similarity of observed flow features with those predicted by computational fluid dynamics (CFD). This report documents the results of the above experiment that specifically address the first three goals. Multiple wing configurations were tested at a chord Reynolds number of 2.4 million. Flow patterns on the surface of the wings and in the region of the wing/fuselage juncture were examined using oil- flow visualization and infrared thermography. A limited number of unsteady pressure sensors on the fuselage around the wing leading and trailing edges were used to identify any dynamic effects of the horseshoe vortex on the flow field. The area of separated flow in the wing/fuselage juncture near the wing trailing edge was observed for all wing configurations at various angles of attack. All of the test objectives were met. The staff of the 14- by 22-foot Subsonic Wind Tunnel provided outstanding support and delivered exceptional value to the experiment, which exceeded expectations. The results of this test will directly inform the planning for the first of a series of instrumented-model tests at the same Reynolds number. These tests will be performed on a slightly larger-scale model with the selected wing, and will include off-body measurements with LDV and PIV, steady and unsteady pressure measurements, and the flow-visualization techniques that are discussed in this report.

Kegerise, Michael A.↗

DCFC + Hydrogen Station Design Optimization [Slides]

Both direct current fast charging (DCFC) and Hydrogen stations are working to create successful long-term business models; however, analyses for DCFC and hydrogen fueling infrastructure are almost always performed separately. This work provides a detailed exploration of the benefit of integrating DCFC and hydrogen stations to lower the total system cost from load balancing and equipment cost sharing. To achieve this we have adapted the REopt optimization framework to simultaneously optimize the design and operation of integrated DCFC and H2 fueling station. Results indicate that 1) combining hydrogen fueling and DCFC stations can significantly reduce lifetime costs compared to separated stations. 2) Co-location with additional site load reduces DCFC costs, however, integration of DCFC with hydrogen provides an even greater cost reduction. 3) Adding PV to combined stations further reduces the lifetime station cost. 4) Capital investments in station combination today can help reduce the cost of operating DCFC tomorrow and 5) product diversification acts as a hedge against variability and enables a more dynamic response to market changes.

08 HYDROGEN↗

Building a machine learning surrogate model for wildfire activities within a global Earth system model

Abstract. Wildfire is an important ecosystem process, influencing land biogeophysical and biogeochemical dynamics and atmospheric composition. Fire-driven loss of vegetation cover, for example, directly modifies the surface energy budget as a consequence of changing albedo, surface roughness, and partitioning of sensible and latent heat fluxes. Carbon dioxide and methane emitted by fires contribute to a positive atmospheric forcing, whereas emissions of carbonaceous aerosols may contribute to surface cooling. Process-based modeling of wildfires in Earth system land models is challenging due to limited understanding of human, climate, and ecosystem controls on fire counts, fire size, and burned area. Integration of mechanistic wildfire models within Earth system models requires careful parameter calibration, which is computationally expensive and subject to equifinality. To explore alternative approaches, we present a deep neural network (DNN) scheme that surrogates the process-based wildfire model with the Energy Exascale Earth System Model (E3SM) interface. The DNN wildfire model accurately simulates observed burned area with over 90 % higher accuracy with a large reduction in parameterization time compared with the current process-based wildfire model. The surrogate wildfire model successfully captured the observed monthly regional burned area during validation period 2011 to 2015 (coefficient of determination, R2=0.93). Since the DNN wildfire model has the same input and output requirements as the E3SM process-based wildfire model, our results demonstrate the applicability of machine learning for high accuracy and efficient large-scale land model development and predictions.

58 GEOSCIENCES↗

Changes in the Magnetic Field Topology and the Dayside/Nightside Reconnection Rates in Response to a Solar Wind Dynamic Pressure Front A Case Study

One of the most significant observations associated with a sharp enhancement in solar wind dynamic pressure, SWP, is the poleward expansion of the auroral oval and the closing of the polar cap. The polar cap shrinking over a wide range of magnetic local times (MLTs), in connection with an observed increase in ionospheric convection and the transpolar potential, led to the conclusion that the nightside reconnection rate is significantly enhanced after a pressure front impact. However, this enhanced tail reconnection has never been directly measured. We demonstrate the effect of a solar wind dynamic pressure front on the polar cap closure, and for the first time, measure the enhanced reconnection rate in the magnetotail, for a case occurring during southward background Interplanetary Magnetic Field (IMF) conditions. We use Polar Ultra-Violet Imager (UVI) measurements to detect the location of the open-closed field line boundary, and combine them with Assimilative Mapping of Ionospheric Electrodynamics (AMIE) potentials to calculate the ionospheric electric field along the polar cap boundary, and thus evaluate the variation of the dayside/nightside reconnection rates. We find a strong response of the polar cap boundary at all available MLTs, exhibiting a significant reduction of the open flux content. We also observe an immediate response of the dayside reconnection rate, plus a phased response, delayed by 15–20min, of the nightside reconnection rate. Finally, we provide comparison of the observations with the results of the Open Geospace General Circulation Model (OpenGGCM), elucidating significant agreements and disagreements

A Boudouridis↗

Modeling ash deposition and shedding during oxy-combustion of coal/rice husk blends at 70% inlet O2

Abstract Co-firing rice husk (RH) and coal with carbon capture using oxy-combustion presents a net carbon negative energy production opportunity. In addition, the high fusion temperature of the non-sticky, silica rich, RH can mitigate ash deposition as well as promote shedding of deposits. To identify the optimum operating conditions, fuel particle sizes, and blend ratios that minimize ash deposition, a Computational Fluid Dynamic methodology with add-on ash deposition and shedding models were employed to predict outer ash deposition and shedding rates during co-combustion of coal/RH in AIR and O 2 /CO 2 (70/30 vol%, OXY70) oxidizer compositions. After ensuring that the fly-ash particle size distributions and particle Stokes numbers near the deposition surface were accurately represented (to model impaction), appropriate models for coal ash and RH ash viscosities that were accurate in the temperature region (1200–1300 K) of interest in this study were identified. A particle viscosity and kinetic energy (PKE) based capture criterion was enforced to model the ash capture. An erosion/shedding criterion that takes the deposit melt fraction and the energy consumed during particle impact into account was also implemented. Deposition rate predictions as well as the deposition rate enhancement (OXY70/AIR) were in good agreement with measured values. While the OXY70 scenario was associated with a significant reduction (60%–70%) in flue gas velocities, it also resulted in larger fly-ash particles. As a result, the PKE distributions of the erosive RH ash were similar in both scenarios and resulted in similar shedding rates.

Energy & Fuels↗

Scalable edge clustering of dynamic graphs via weighted line graphs

Timestamped relational datasets consisting of records (or connections) between pairs of entities are ubiquitous in network science. For applications like peer-to-peer communication, email, various social network interactions, and computer network security, it is useful to organize these records into groups based on how and when they are occurring. Weighted line graphs offer a natural way to model how records are related in such datasets but for large real-world graph topologies, building and utilizing the line graph is prohibitively expensive. Here, we present the framework to cluster the edges of a dynamic graph via the associated line graph that contains two major contributions. The first is a method to work with the line graph implicitly and the second is a distributed scale implementation of an agglomerative hierarchical graph clustering algorithm. We outline a novel hierarchical dynamic graph edge clustering approach that efficiently breaks massive relational datasets into small sets of edges containing events at various timescales. This is in stark contrast to traditional graph clustering algorithms that prioritize highly connected (clique-like) community structures. Our approach relies on constructing a sufficient subgraph of a weighted line graph and applying a hierarchical agglomerative clustering. This approach is related to scalable techniques from spatial clustering, nonlinear-dimension reduction, topological data analysis, and draws particular inspiration from HDBSCAN. As an edge clustering, this method yields an overlapping node clustering. Our algorithm is parallelizable and we demonstrate efficient clustering of a billion-scale, real-world dynamic graph into small edge sets that correlate in topology and time. The entire clustering process for a graph with tens of billions of edges takes just a few minutes of run time on 256 nodes of a distributed compute environment. We argue how the output of the edge clustering is useful for a multitude of data visualization and powerful machine learning tasks, both involving the original massive dynamic graph data and metadata associated with the nodes and edges. Finally, we describe how this approach can be extended to dynamic hypergraphs and dynamic graphs/hypergraphs with unstructured data living on vertices and edges.

Data Analysis↗

Wildfire management decisions outweigh mechanical treatment as the keystone to forest landscape adaptation

Modern land management faces unprecedented uncertainty regarding future climates, novel disturbance regimes, and unanticipated ecological feedbacks. Mitigating this uncertainty requires a cohesive landscape management strategy that utilizes multiple methods to optimize benefits while hedging risks amidst uncertain futures. We used a process-based landscape simulation model (LANDIS-II) to forecast forest management, growth, climate effects, and future wildfire dynamics, and we distilled results using a decision support tool allowing us to examine tradeoffs between alternative management strategies. We developed plausible future management scenarios based on factorial combinations of restoration-oriented thinning prescriptions, prescribed fire, and wildland fire use. Results were assessed continuously for a 100-year simulation period, which provided a unique assessment of tradeoffs and benefits among seven primary topics representing social, ecological, and economic aspects of resilience. Projected climatic changes had a substantial impact on modeled wildfire activity. In the Wildfire Only scenario (no treatments, but including active wildfire and climate change), we observed an upwards inflection point in area burned around mid-century (2060) that had detrimental impacts on total landscape carbon storage. While simulated mechanical treatments (~ 3% area per year) reduced the incidence of high-severity fire, it did not eliminate this inflection completely. Scenarios involving wildland fire use resulted in greater reductions in high-severity fire and a more linear trend in cumulative area burned. Mechanical treatments were beneficial for subtopics under the economic topic given their positive financial return on investment, while wildland fire use scenarios were better for ecological subtopics, primarily due to a greater reduction in high-severity fire. Benefits among the social subtopics were mixed, reflecting the inevitability of tradeoffs in landscapes that we rely on for diverse and countervailing ecosystem services. This study provides evidence that optimal future scenarios will involve a mix of active and passive management strategies, allowing different management tactics to coexist within and among ownerships classes. Our results also emphasize the importance of wildfire management decisions as central to building more robust and resilient future landscapes.

54 ENVIRONMENTAL SCIENCES↗

Process Heating Assessments Using DOE’s Manufacturing Energy Assessment Software for Utility Reduction (MEASUR) Tool Suite

Process heating is the most energy-intensive manufacturing process for most sectors of industry. To quantify energy savings from various energy conservation measures, the Department of Energy (DOE) sponsored the development of the Process Heating Assessment and Survey Tool (PHAST) and similar tools for other industrial systems in the early 2000s. It has been used extensively in the Save Energy Now Program’s Energy Savings Assessments and the Better Plants Program’s In-Plant Trainings. Since the initial development of the legacy tools, both computer operating systems and software development have evolved significantly. Thus, DOE has invested in the modernization of PHAST and other legacy software tools to create the Manufacturing Energy Assessment Software for Utility Reduction (MEASUR) tool suite. MEASUR offers a collection of software tools that can aid manufacturing facilities in improving the efficiency of energy systems and equipment (specifically pumps, fans, steam, and process heating) and in conducting “Energy Treasure Hunts”. Eventually, the tool will also add compressed air and process cooling systems. The Process Heating Assessment (PHA) module of MEASUR is an upgrade of the PHAST tool. PHA provides the means to model fuel-fired, steam-based, and electric process heating systems, covering process heating for most industrial plants in manufacturing sector. It also includes several key upgrades, including the ability to consider multi-component charge loads and account for several different areas of energy losses. The new tool includes a comprehensive flue gas calculator to quantify available heat and heat loss for various gaseous, liquid, and solid fuels and new heat loss calculators. It generates a report and a dynamic Sankey diagram to show the energy consumption in various areas of energy use. MEASUR has significantly improved the user experience by adopting a modern software design. This paper details the structure and workflow of PHA and presents a real-world case study to demonstrate energy savings quantification and MEASUR’s outstanding reporting capabilities.

Nimbalkar, Sachin U.↗

Modeling diurnal and annual ethylene generation from solar-driven electrochemical CO 2 reduction devices

Integrated solar fuels devices for CO 2 reduction (CO 2 R) are a promising technology class towards reducing carbon emissions. Designing integrated CO 2 R solar fuels devices requires careful co-design of electrochemical and photovoltaic components as well as consideration of the diurnal and seasonal effects of solar irradiance, temperature, and other meteorological factors expected for ‘on-sun’ deployment. Using a photovoltaic-electrochemical (PV-EC) platform, we developed a temperature and potential-dependent diurnal and annual model using experimentally-determined CO 2 R performance of Cu-based electrocatalysts, local meteorological data from the National Solar Radiation Database (NSRD), and modeled performance of commercial c-Si PVs. Here we simulated gaseous diurnal product outputs with and without the effects of ambient temperature. From these outputs, we observed seasonal variation in gaseous product generation, with up to two-fold increases in ethylene productivity between the Winter and Summer, analyzed the consequences of dynamic cloud coverage, and identified periods where device cooling/heating mechanisms could be implemented to maximize ethylene generation. Finally, we modeled the annual ethylene generation for a scaled 1 MW solar farm at three different locations (Beijing, CN; Sydney, AUS; Barstow, CA) to determine the consequences of local meteorological climates on PV-EC CO 2 R product output, recording a maximum ethylene output of 18.5 tonne per year at Barstow. Overall, this model presents a critical tool for streamlining the translation of experimental solar-driven electrochemical research to real-world implementation.

14 SOLAR ENERGY↗

Rapid iterative reanalysis for automated design

A method for iterative reanalysis in automated structural design is presented for a finite-element analysis using the direct stiffness approach. A basic feature of the method is that the generalized stiffness and inertia matrices are expressed as functions of structural design parameters, and these generalized matrices are expanded in Taylor series about the initial design. Only the linear terms are retained in the expansions. The method is approximate because it uses static condensation, modal reduction, and the linear Taylor series expansions. The exact linear representation of the expansions of the generalized matrices is also described and a basis for the present method is established. Results of applications of the present method to the recalculation of the natural frequencies of two simple platelike structural models are presented and compared with results obtained by using a commonly applied analysis procedure used as a reference. In general, the results are in good agreement. A comparison of the computer times required for the use of the present method and the reference method indicated that the present method required substantially less time for reanalysis. Although the results presented are for relatively small-order problems, the present method will become more efficient relative to the reference method as the problem size increases. An extension of the present method to static reanalysis is described, ana a basis for unifying the static and dynamic reanalysis procedures is presented.

Bhatia, K. G.↗

Aerodynamic Testing of the Orion Launch Abort Tower Separation with Jettison Motor Jet Interactions

The aerodynamic database for the Orion Launch Abort System (LAS) was developed largely from wind tunnel tests involving powered jet simulations of the rocket exhaust plumes, supported by computational fluid dynamics (CFD) simulations. The LAS contains three solid rocket motors used in various phases of an abort to provide propulsion, steering, and Launch Abort Tower (LAT) jettison from the Crew Module (CM). This paper describes a pair of wind tunnel experiments performed at transonic and supersonic speeds to determine the aerodynamic effects due to proximity and jet interactions during LAT jettison from the CM at the end of an abort. The tests were run using two different scale models at angles of attack from 150deg to 200deg , sideslip angles from -10deg to +10deg , and a range of powered thrust levels from the jettison motors to match various jet simulation parameters with flight values. Separation movements between the CM and LAT included axial and vertical translations as well as relative pitch angle between the two bodies. The paper details aspects of the model design, nozzle scaling methodology, instrumentation, testing procedures, and data reduction. Sample data are shown to highlight trends seen in the results.

Rhode, Matthew N.↗

Large Eddy Simulation of High-Speed, Premixed Ethylene Combustion

A large-eddy simulation / Reynolds-averaged Navier-Stokes (LES/RANS) methodology is used to simulate premixed ethylene-air combustion in a model scramjet designed for dual mode operation and equipped with a cavity for flameholding. A 22-species reduced mechanism for ethylene-air combustion is employed, and the calculations are performed on a mesh containing 93 million cells. Fuel plumes injected at the isolator entrance are processed by the isolator shock train, yielding a premixed fuel-air mixture at an equivalence ratio of 0.42 at the cavity entrance plane. A premixed flame is anchored within the cavity and propagates toward the opposite wall. Near complete combustion of ethylene is obtained. The combustor is highly dynamic, exhibiting a large-scale oscillation in global heat release and mass flow rate with a period of about 2.8 ms. Maximum heat release occurs when the flame front reaches its most downstream extent, as the flame surface area is larger. Minimum heat release is associated with flame propagation toward the cavity and occurs through a reduction in core flow velocity that is correlated with an upstream movement of the shock train. Reasonable agreement between simulation results and available wall pressure, particle image velocimetry, and OH-PLIF data is obtained, but it is not yet clear whether the system-level oscillations seen in the calculations are actually present in the experiment.

Ramesh, Kiran↗