Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “maximum”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Beyond the Holographic Entropy Cone via Cycle Flows

Motivated by bit threads, we introduce a new prescription for computing entropy vectors outside the holographic entropy cone. By utilizing cycle flows on directed graphs, we show that the maximum cycle flow associated to any subset of vertices, which corresponds to a subsystem, manifestly obeys purification symmetry. Furthermore, by restricting ourselves to a subclass of directed graphs, we prove that the maximum cycle flow obeys both subadditivity and strong subadditivity, thereby establishing it as a viable candidate for the entropy associated to the subsystem. Lastly, we demonstrate how our model generalizes the entropy vectors obtainable via conventional flows in undirected graphs, as well as conjecture that our model similarly generalizes the entropy vectors arising from hypergraphs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

The substantial role of May soil temperature over Central Asia for summer surface air temperature variation and prediction over Northeastern China

The slowly varying soil temperature can exert local and nonlocal influences on regional climate system, and may thus provide a critical source of subseasonal-to-seasonal climate prediction. In this study, we identify that soil temperature in May over the key region of Central Asia (42 °N–50 °N, 62 °E–80 °E, KRCA) from Noah, Mosaic, CLM and ERA-interim datasets is closely linked to variations of the surface air temperature, daily maximum temperature and hot days over Northeastern China in summer (June-July-August), with correlation coefficients of regional average detrended time series ranging from 0.42 to 0.54, and all significant at the 99% confidence level for the period of 1979-2018. The possible physical mechanism behind the substantial downstream impacts of soil temperature over Central Asia are explored via diagnostical analysis combined with regional climate model experiments. Warmer soil temperature in May over the KRCA tends to cause positive anomalies of geopotential height in summer over Northeastern China through the Rossby wave propagation, and associated stronger subsidence warming, less cloud cover, more solar radiation reaching the surface, higher planetary boundary layer, and stronger thermal advection at 850 hPa, which provide favorable conditions for warmer surface air temperature particularly in the daytime as well as more hot days. Here this study further reveals that soil temperature over Central Asia in May makes important contribution to prediction of summer surface air temperature, daily maximum temperature and hot days over Northeastern China in terms of regional average time series and spatial patterns. Our findings highlight the previously-unknown substantial role of antecedent soil temperature condition over Central Asia for summer surface air temperature variation and prediction over Northeastern China.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Recalibration of missing low-frequency variability and trends in the North Atlantic Oscillation

Abstract Multi-decadal trends in the wintertime North Atlantic Oscillation (NAO) are under-represented by coupled general circulation models (CGCMs), consistent with a lack of autocorrelation in their NAO index series. This study proposes and tests two simple “reddening” approaches for correcting this problem in simulated indices based on simple one parameter short-term (AR; Auto-Regressive order 1) and long-term (FD; Fractional-Difference) time series filters. Using CGCMs from the Coupled Model Intercomparison Project Phase 6 (CMIP6), the FD filter successfully improves the autocorrelation structure of the NAO, and in turn the simulation of extreme trends, while the AR filter is less successful. The 1963–1993 NAO trend is the maximum 31-year trend in the historical period. Raw CGCMs underestimate the likelihood of this trend by a factor of ten but this discrepancy is corrected after reddening. CMIP6 future projections show that long-term (2024–2094) NAO ensemble mean trends systematically increase with the magnitude of radiative forcing: -2.4 to 3.5 hPa/century for low-to-high forcing after reddening (more than double the range using raw output). The related likelihood of future maximum 31year trends comparable to 1963–1993 ranges from 3 to 7% whereas none of these CMIP6 projections simulate this without reddening. Near-term projections of the next 31 years (2024–2054) are less sensitive than long term trends to the future scenario, showing weak-to-no forced trend. However, reddening increases the ensemble range by 74% (to +/-1 standard deviation/decade), which could increase/decrease regional climate change signals in the Northern Hemisphere by magnitudes that are underestimated when using raw CGCM output.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

Developing a Control Strategy for Minimum Airflow Setting Considering CO2 Level and Energy Consumption in a Variable Air Volume System

In an office building equipped with a Variable Air Volume (VAV) system, this paper introduces a novel method for controlling the minimum supply airflow fraction in each zone’s VAV box, having a capability to consider indoor CO2 level and energy consumption. The EnergyPlus simulation using the medium office prototype model was employed, which evaluated the performance of the energy and CO2 concentration for five VAV box airflow control strategies. The paper focuses on CO2 concentration-based airflow control method and compares it with other four methods including conventional single-max, reduced minimum single-max, demand-controlled ventilation(DCV), and dualmax control methods according to guidelines and common practices. The newly proposed control strategy directly correlates the minimum airflow fraction to CO2 concentration. A general trend emerged when comparing CO2 concentrations—lower minimum airflow fractions were associated with higher concentrations. The proposed control method effectively maintained low CO2 concentrations and enabled a lower airflow fraction contributing to energy consumption reduction. It was confirmed that heating energy consumption in climate zone 4A, 5B, and 6A showed a maximum saving of approximately 30% compared to the conventional single-max and dual max control strategies. It was found that cooling energy consumption in climate zone 4A and 6A can achieve a maximum saving of approximately 10% compared to the conventional control strategies. The proposed CO2 concentration-based control logic is promising as it not only improves the indoor air quality lowering the CO2 concentration in the occupied spaces, but also contributes to HVAC energy savings.

Lee, Jong Man↗
DOE OSTI2021

Extraction of gold, mercury and lead on CL resin in HCl

The uptake behavior of 197 mHg, 194 Au and 212 Pb in HCl on CL resin (Triskem) was studied. The extraction of mercury and gold is high (D w > 5,000) at all concentrations of HCl, with a maximum extraction of gold at 4 M HCl (D w ~ 53,000) and maximum extraction of mercury at 0.01 M HCl (D w ~ 44,000). Lead has moderate uptake on CL resin (D w ~ 100) from lower concentrations of HCl (0.001 to 0.05 M) and no uptake on the resin in concentrations ≥ 0.4 M. The uptake kinetics are similar for all three isotopes with high weight distribution ratios (>1,000 for 197 mHg, 194 Au and > 100 for 212 Pb) in ~ 30 s and equilibrium established at ~ 5 min. A column study was performed to demonstrate that 197m Hg, 194 Au and 212 Pb could be retained on CL resin in dilute HCl and 212 Pb eluted in 1 M HCl with no detectable 197m Hg or 194 Au. Elution of 197 mHg and 194 Au was performed with concentrated HNO 3 , but with significant losses to the column, particularly for 194 Au.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Extraction of radium and actinium with dibenzo-21-crown-7 resin and Rose Bengal

The extraction of 224 Ra and 228 Ac with dibenzo-21-crown-7 resin was studied in solutions with Rose Bengal and identical solutions without Rose Bengal (blanks) from pH 3 to 12. Batch studies demonstrated that actinium is extracted at high pH (> 10) from both solutions but with a significantly higher extraction from the Rose Bengal solutions (maximum extraction ~ 60,000) than the blank solutions (maximum extraction ~ 2,000). In the Rose Bengal solutions, there was a weak extraction of actinium at lower pH (4–10) as well. There was no uptake of radium from the blank solutions, and a weak extraction at all pH levels from the Rose Bengal solutions. Kinetics studies were performed for actinium only, as radium had no significant extraction. Two column studies are demonstrated for the separation of radium from actinium. These separations are moderately effective with ~ 78% of the actinium recovered with no detectable radium.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Processing Compressed Expanded Natural Graphite for Phase Change Material Composites

Phase change materials (PCMs) are used in various thermal energy storage applications but are limited by their low thermal conductivity. One method to increase conductivity involves impregnating organic PCMs into highly porous conductive matrix materials. Of these materials, compressed expanded natural graphite (CENG) matrices have received the most attention. Despite this attention, the effect that CENG processing has on PCM saturation and overall matrix thermal conductivity has not been fully investigated. Therefore, the effect of the heat treatment process used to expand intercalated graphite flakes is evaluated here. Higher heat treatment temperatures yielded higher saturation rates and overall saturation at similar matrix porosities. For example, increasing temperature from 300 °C to 700 °C resulted in approximately 60%-70% increase in pore saturation after 100 minutes of soaking. The exposure time to heat treatment had less of an effect on PCM saturation. The exposure time had negligible effect above 30 min and above 500 °C heating temperatures. However, because the expanded graphite was found to oxidize around 700 °C, the use of longer exposure time in manufacturing applications can be beneficial if a shortened impregnation time is needed. Heat treatment conditions did not impact thermal conductivity. The composite latent heat of fusion was also reduced approximately proportionally to the PCM mass fraction. A local maximum in axial thermal conductivity was observed at around 83% porosity, which is similar to previous studies. The observed conductivity at this maximum was a factor of 81 times greater than the conductivity of the PCM.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Fluid Dynamics Effects on Microstructure Prediction in Single-Laser Tracks for Additive Manufacturing of IN625

Single-track laser fusion were simulated using a heat-transfer-solidification-only (HTS) model and its extension with fluid dynamics (HTS_FD) model using a parallel open-source code, which included laminar fluid dynamics, flat-free surface of the molten alloy, heat transfer, phase-change, evaporation, and surface tension phenomena. The results illustrate that the fluid dynamics affects the solidification and ensuing microstructure. For the HTS_FD simulations, thermal gradient, G was found to exhibit a maximum at the extremity of the solidified pool ( i.e. , at the free surface), while for HTS simulations, G exhibited a maximum around the entire edge of the solidified pool. HTS_FD simulations predicted a wider range of cooling rates than the HTS simulations, exhibited an increased spread in the solidification speed, V variation within the melt-pool with respect to the HTS model results. Primary dendrite arm spacing (PDAS) were evaluated based on power law correlations and marginal stability theory models using the ( G , V ) from HTS and HTS_FD simulations to quantify the effect of the fluid dynamics on the microstructure. At low-laser powers and low-scan speeds, the PDAS obtained with the fluid dynamics model (HTS_FD) was larger by more than 30 pct with respect to the PDAS calculated with the simple HTS model. A new PDAS correlation, i.e. , \( \lambda_{1} \left[ {\mu {\text{m}}} \right] = 832\;G\left[ {\text{K/m}} \right]^{ - 0.5} V\left[ {\text{m/s}} \right]^{ - 0.25} \) , which uses the ( G , V ) results from the HTS_FD model was developed and validated against experimental results.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Residual Stress in Cold Spray SS304L Measured Via Neutron Diffraction and Comparison of Analytical Models to Predict the Residual Stress

Here, this study employs neutron diffraction to investigate the relationship between residual stress and coating thickness in cold sprayed 304L austenitic stainless steel. Results show that shot peening predominantly impacts the residual stress profile, leading to substantial in-plane compressive force. The impact of laser heating, a widely used method to alter cold spray's microstructural properties, on the coating's residual stress is also analyzed. The findings indicate that the maximum compressive residual stress in the in-plane component is mainly independent of coating thickness, which suggests that the material properties determine the maximum residual stress. The cold sprayed deposits possessed compressive, nearly biaxial strain and stresses. After laser heating, these stresses were replaced by tensile residual stresses. Two analytical models, the Tsui and Clyne and the Boruah models, for predicting residual stresses are also evaluated, and both models provide reasonable fits to the experimental data. At this point, the deviations between the experimental results and the models are principally caused by the inability of the current models to address plastic deformation and relaxation, and the residual stresses generated by thermal gradients.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Leveraging Extremal Dependence to Better Characterize the 2021 Pacific Northwest Heatwave

Abstract In late June, 2021, a devastating heatwave affected the US Pacific Northwest and western Canada, breaking numerous all-time temperature records by large margins and directly causing hundreds of fatalities. The observed 2021 daily maximum temperature across much of the U.S. Pacific Northwest exceeded upper bound estimates obtained from single-station temperature records even after accounting for anthropogenic climate change, meaning that the event could not have been predicted under standard univariate extreme value analysis assumptions. In this work, we utilize a flexible spatial extremes model that considers all stations across the Pacific Northwest domain and accounts for the fact that many stations simultaneously experience extreme temperatures. Our analysis incorporates the effects of anthropogenic forcing and natural climate variability in order to better characterize time-varying changes in the distribution of daily temperature extremes. We show that greenhouse gas forcing, drought conditions and large-scale atmospheric modes of variability all have significant impact on summertime maximum temperatures in this region. Our model represents a significant improvement over corresponding single-station analysis, and our posterior medians of the upper bounds are able to anticipate more than 96% of the observed 2021 high station temperatures after properly accounting for extremal dependence. Supplementary materials accompanying this paper appear online.

Zhang, Likun (ORCID:0000000154584556)↗
DOE OSTI2024

Reductive catalytic fractionation of cotton stalks: catalytic strategy for tuning the selectivity of phenolic monomers

Pretreatment of lignocellulosic biomass is a primary step to delink lignin from the lignin-carbohydrate complex for bioethanol production and other value-added products. Presently, these processes yield lignin (technical) with significant structural changes compared to native lignin and make it difficult to valorize. Reductive catalytic fractionation (RCF) is an advanced pretreatment process to valorize native lignin to selective phenolic monomers before capitalizing on carbohydrates. Herein, cotton stalks (CS) were pretreated with ethanol: water mixture catalyzed by bimetallic, Ru-Ni/HY (RNY) and Ru-Fe/HY (RFY), and trimetallic, Ru-Fe-Ni/HY (RFNY), catalysts using the endogenously produced H 2 . The favorable catalytic properties of the RFNY catalyst, including high reducibility, acidity, metal dispersion, the synergistic effect of all the metals, and appropriate pore and particle size resulted in the highest catalytic activity. Maximum 92% delignification efficiency was achieved with 91% holocellulose retention. The extracted lignin after depolymerization in the same step produced a 19 wt.% yield of phenolic monomers. Increasing the duration of the reaction from 2 to 6 h increased the selectivity to saturated alkyl chain phenolics (1 to 16 wt.%) at the expense of hemicellulose (50 to 9%). The fate of carbohydrates left behind after the RCF (CS-RCF) process was analyzed through enzymatic saccharification and compared with raw CS and acid-pretreated CS. In conclusion, CS-RCF yielded maximum glucose yield (45 g/L) within an incubation period of 72 h.

09 BIOMASS FUELS↗
DOE OSTI2024

Newton-Okounkov bodies of chemical reaction systems

Despite their noted potential in polynomial-system solving, there are few concrete examples of Newton-Okounkov bodies arising from applications. Accordingly, in this paper, we introduce a new application of Newton-Okounkov body theory to the study of chemical reaction networks and compute several examples. An important invariant of a chemical reaction network is its maximum number of positive steady states Here, we introduce a new upper bound on this number, namely the ‘Newton-Okounkov body bound’ of a chemical reaction network. Further, through explicit examples, we show that the Newton-Okounkov body bound of a network gives a good upper bound on its maximum number of positive steady states. We also compare this Newton-Okounkov body bound to a related upper bound, namely the mixed volume of a chemical reaction network, and find that it often achieves better bounds.

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

An energetic link between order and strength in metals: A nanocrystalline strength limit in high-entropy alloys and intermetallic compounds

The metallurgy and materials communities have long understood and exploited fundamental links between chemical and structural ordering in metallic solids to tailor their mechanical properties. We extend these ideas to include prediction of the nanocrystalline strength limit in high-entropy alloys and intermetallic compounds, where a breakdown occurs in the classical Hall-Petch strengthening behavior. The highest reported strength achievable through alloying has rapidly climbed and given rise to new classifications of materials with extraordinary properties, with a notable case being nanocrystalline metals. High-entropy alloys (chemically disordered, concentrated solid solutions) and intermetallic compounds are two boundary cases of how tailored order can be used to manipulate mechanical behavior. Here, we show that the complex electronic-structure mechanisms governing the peak strength of alloys and pure metals can be reduced to a few physically meaningful parameters based on their atomic arrangements and used – with no fitting parameters – to predict the maximum strength of these materials. This includes a generalized energy-based accounting for the degree of structural and chemical ordering that allows for rapid and reasonably accurate prediction of peak strength (validated in the nanocrystalline limit) as a function of temperature. Predictions of maximum strength based on the activation energy (with all materials properties derived from DFT calculations or experiments) for a stress-driven transition to an amorphous state is shown to accurately describe the breakdown in Hall-Petch behavior at the smallest crystallite sizes for pure metals, intermetallic compounds, high-entropy alloys, and metallic glasses. Further, this activation energy is also shown to be directly proportional to interstitial electronic charge density, which is a good predictor of ductility, stiffness (moduli), and phase stability in high-entropy alloys and solid metals generally. The proposed framework suggests the possibility of coupling ordering and intrinsic strength to mechanisms like dislocation nucleation, hydrogen embrittlement, and transport properties, such as through correlations between the activation energies for amorphization with stacking-fault and grain boundary energies. It additionally opens the prospect for greatly accelerated structural materials design and development to address materials challenges limiting more sustainable and efficient use of energy.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Which way does the dendrite grow? Competition among epitaxy, preferred growth direction, and thermal gradients in powder bed fusion additive manufacturing

The as-processed microstructure of metal alloy parts manufactured through laser powder bed fusion (LPBF) is heavily derived from the cellular dendritic solidification. The growth direction of dendrites within the melt pool is determined through competition among epitaxial growth, preferred growth directions, and maximum thermal gradients. However, the dominant factor and the specific role of each in developing melt pool microstructures remain unknown. Here, in this study, we performed single laser track scans on an SS316L single crystal substrate and combined experimental characterization of microstructure and crystal orientations with Computational Fluid Dynamics simulations of thermal gradients to evaluate the role of each factor in determining dendritic growth direction and evolution. Our results reveal that epitaxial growth dominates microstructure development by preferentially growing along a single 〈100〉 variant of the single crystal substrate adjacent to the melt pool boundary. Under LPBF’s highly curved and rapidly evolving thermal field, this preferential dendrite variant selection and its continued growth from the melt pool boundary to the centerline are governed by the local temperature gradient magnitude at the solid-liquid interface, rather than by the instantaneous maximum temperature gradient direction alone. Using these findings, we successfully predict changes in the dendrite growth direction with changing laser scan direction on a single crystal substrate, and show that the geometric melt pool centerline can deviate from the microstructural centerline because asymmetric local temperature gradient magnitudes transiently limit growth, resulting in different dendrite travel distances on each side of the melt pool.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Geometrical effects on the measured neutron residual stress in additively manufactured cold spray AA6061

Aluminum cold spray additive builds were made with varying wall widths and geometries to explore the effects of build geometry on residual stress development in the additively sprayed material. Deposition rates of over 750 g/hr were achieved without any cracking or delamination of the additively manufactured components. Neutron residual stress measurements revealed a maximum tensile residual stress of 41 MPa at the interface of the substrate and a maximum compressive residual stress of −35 MPa in the cold sprayed material. While there was not found to be significant residual stress changes between the cold sprayed walls, the cold sprayed cylinder had varying residual stress throughout the diameter with the outer diameter having a higher tensile residual stress than the interior. The low residual stress and minimal geometrical dependency demonstrate further part size is possible without concerns for component failure in AA6061. The application of common coating residual stress models to full scale additive builds is evaluated and discussed. The lack of residual stress evolution was studied by in-situ coating property measurement that revealed that the deposition stress of the AA6061 sprays reduced with increased layer thickness. In conclusion, this study provides insight into the evolution of residual stress in cold spray additive components and how both the mechanics and material responses differ from coatings.

Additive manufacturing↗
DOE OSTI2025

Karhunen–Loève deep learning method for surrogate modeling and approximate Bayesian parameter estimation

We evaluate the performance of the Karhunen-Loève Deep Neural Network (KL-DNN) framework for surrogate modeling and approximate Bayesian parameter estimation in partial differential equation models. In the surrogate model, the Karhunen-Loève (KL) expansions are used for the dimensionality reduction of the number of unknown parameters and variables, and a deep neural network is employed to relate the reduced space of parameters to that of the state variables. The KL-DNN surrogate model is used to formulate a maximum-a-posteriori-like least-squares problem, which is randomized to draw samples of the posterior distribution of the parameters. We test the proposed framework for a hypothetical unconfined aquifer via comparison with the forward MODFLOW and inverse PEST++ iterative ensemble smoother (IES) solutions as well as the state-of-the-art Fourier neural operator (FNO) and deep operator networks (DeepONets) operator learning surrogate models. Our results show that the KL-DNN surrogate model outperforms FNO and DeepONet for forward predictions. For solving inverse problems, the randomized algorithm provides the same or more accurate Bayesian predictions of the parameters than IES as evidenced by the higher log-predictive probability of both the estimated parameter field and the forecast hydraulic head. The posterior mean obtained from the randomized algorithm is closer to the reference parameter field than that obtained with FNO as the maximum a posteriori estimate.

Approximate Bayesian inference↗
DOE OSTI2021

A validated thermal and biological model for predicting algal productivity in large scale outdoor cultivation systems

Proper assessment of the sustainability of algal products requires an understanding of algal growth rate with regional and temporal resolution. This study developed a bulk growth model that utilizes geospatial data, reactor geometry inputs, and a maximum of six strain-specific parameters to predict algal productivity with temporal and regional resolution. The model requires local weather data to calculate the time-resolved culture temperature which is combined with a biological model based on maximum theoretical carbon fixation to predict the time-resolved algal concentration. Validation of the model against experimental data illustrates an average ac- curacy of 0.27% ± 5.32% when comparing algal concentration at harvest over 45 days of cultivation data for Galdieria sulphuraria 5587.1, and Galdieria sulphuraria Soos grown in photobioreactors (PBRs) in Mesa, AZ (33.4152° N, 111.8315° W). For open raceway ponds (ORPs), an accuracy of 0.9% ± 2.35% was achieved when comparing algal concentration at harvest over 50 days of cultivation data for Chlorella vulgaris, Desmodesmus intermedius, Nannochloropsis oceanica grown in Atlanta, GA (33.7490° N, 84.3880° W) and Mesa, AZ (33.4152° N, 111.8315° W). The validation trials spanned three seasons (winter, spring, and summer) to ensure model ac- curacy throughout the year, and concentration between harvests was modeled with hourly resolution. The validated model is used to assess the productivity at several case-study locations in the United States with Chlorella vulgaris, grown in open raceway ponds and Galdieria sulphuraria grown in vertical flat panel photo- bioreactors. In conclusion, the model has been published as an open-source tool, in an effort to increase the fidelity of future studies that rely on outdoor algae cultivation for feedstock production and allow for time resolved results for future techno-economic and life cycle assessments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Large-eddy simulation study on cycle-to-cycle variation of knocking combustion in a spark-ignition engine

The cycle-to-cycle variation in the knock intensity is commonly encountered under abnormal combustion conditions. The severity of these abnormal combustion events can vary significantly, and the efficiency of engines at high loads is limited in practice by heavy knocking phenomena. Since, a thorough analysis of such recurrent but non-cyclic phenomena via experiments alone becomes highly cumbersome, in the present work, a multi-cycle large-eddy simulation study was performed to quantitatively predict cyclic variability in the combustion process and cyclic knock intensity variability in a direct injection spark-ignition engine. To account for the turbulence-chemistry interaction effects on flame propagation, the G-equation combustion model was used. Detailed chemistry was solved outside the flame front with a toluene primary reference fuel skeletal kinetic mechanism. For both the mild knock and heavy knock conditions, the numerical results were validated against experimental measurements. Based on the simulation results, a correlation analysis was performed considering combustion phasing, peak cylinder pressure and maximum amplitude of pressure oscillation. Furthermore, a detailed three-dimensional spatial analysis illustrated the evolution of auto-ignition kernel development and propagation of pressure waves during knocking combustion for three typical cycles with different knock intensities. In this process, it was found that an early occurrence of auto-ignition in the end gas was prone to high knock intensity. Although multiple auto-ignition kernels were observed in different cycles, the degree of coupling between chemical heat release and pressure waves varied, thereby leading to different maximum amplitude of pressure oscillation values.

42 ENGINEERING↗