Flow Boundary Research Articles

Optimization of Well Placement in Carbon Capture and Storage (CCS): Bayesian Optimization Framework under Permutation Invariance

Carbon Capture and Storage (CCS) stands as a pivotal technological stride toward a sustainable future, with the practice of injecting supercritical CO2 into subsurface formations being already an established practice for enhanced oil recovery operations. The overarching objective of CCS is to protract the operational viability and sustainability of platforms and oilfields, thereby facilitating a seamless transition towards sustainable practices. This study introduces a comprehensive framework for optimizing well placement in CCS operations, employing a derivative-free method known as Bayesian Optimization. The development plan is tailored for scenarios featuring aquifers devoid of flow boundaries, incorporating production wells tasked with controlling pressure buildup and injection wells dedicated to CO2 sequestration. Notably, the wells operate under group control, signifying predefined injection and production targets and constraints that must be adhered to throughout the project’s lifespan. As a result, the objective function remains invariant under specific permutations of the well locations. Our investigation delves into the efficacy of Bayesian Optimization under the introduced permutation invariance. The results reveal that it demonstrates critical efficiency in handling the optimization task extremely fast. In essence, this study advocates for the efficacy of Bayesian Optimization in the context of optimizing well placement for CCS operations, emphasizing its potential as a preferred methodology for enhancing sustainability in the energy sector.

The comprehensive research on the transcritical CO2 thermal management system used in thermostatic transport vehicles

Natural fluid CO2 has become the preferred target of refrigerant alternative in the electrified transportation field in recent years to replace positive temperature coefficient (PTC) heaters, reduce energy consumption, and respect the international ban on Hydrofluorocarbons (HFCs) refrigerants. Thermostatic transport vehicles (TTVs) are well suited to exploit the benefits of CO2 due to their working conditions with a broad range of ambient temperatures, extensive transit distances, and substantial heating capacity requirements in cold climates. In this paper, a transcritical CO2 heat pump air-conditioning system used in TTV was constructed initially with proper control logic. Then following the assessment of the system’s steady-state performance, the dynamic operation process under high thermal inertia was thoroughly examined and analyzed in detail. Ultimately, the optimization of the dynamic response process was realized under varying heat flow boundary conditions. The primary results demonstrated that the rationally designed CO2 thermal management system performed admirably in various operating conditions, particularly in the heating mode where the coefficient of performance (COP) exceeded 3.2. Addressing the high-frequency fluctuation caused by the cargo thermal inertia, the temperature could be maintained with only two starts and stops within 3×105 seconds when the system control logic was refined. The optimized thermal management system of TTV efficiently allowed cargo to achieve the proper temperature as quickly as possible and sustain stability. This study provides an effective theoretical basis for the promotion of the transcritical CO2 thermal management system in the field of TTVs.

The dynamic evolution of powder flow and wall normal stress in different flow pattern silos

Major safety accidents involving silos in industrial processes are closely associated with abnormal wall stresses. The effects of different flow patterns and particle velocity distribution on wall stresses are investigated during silo discharge. The existence of particle velocity retardation layer near the wall boundary in the mass flow was observed by tracer particles and high-speed camera. The silo discharge undergoes a transformation from funnel flow to mass flow at a hopper half angle of 30°. As the outlet diameter increases, the time point of the flow pattern transformation becomes more and more later. The evolution of particle velocity in the central of the funnel flow is related to the outlet velocity wave propagation and the V-shaped surface expansion. The relationship between stress fluctuations and velocity characteristics is established. The flow channel simultaneously expands upwards as the velocity wave propagates and generates a periodic fan-shaped velocity wave at the top. The formation of stress concentration zone at the bin/hopper transition was observed.

Freestream turbulence effects on low Reynolds number NACA 0012 airfoil laminar separation bubble and lift generation

Laminar separation bubbles (LSB's) over the suction surface of a wing at low Reynolds number (O(104)−O(106) based on the airfoil chord length) can significantly affect the aerodynamic performance of the wing, and pose a unique challenge for the predictive capabilities of simulation tools due to their high sensitivity to flow environments and wing surface conditions. In this work a series of two-dimensional (2D) and three-dimensional (3D) low-order, and high-order accurate unstructured-grid-based numerical methods with varying model fidelity levels were used to study LSB physics over a NACA 0012 airfoil both in a clean freestream and in a turbulent freestream at a chord-based Reynolds number of 12,000. Lift production and time-averaged flow fields were compared with available experimental results. A major discovery is that in clean freestream flow a 3D high-order numerical scheme is necessary to capture LSB physics. This is due to the sensitivity of LSB-induced laminar-turbulent transition to flow conditions and boundary geometry at low Reynolds number. In freestream flows with moderate background turbulence (∼5%), 2D simulations failed to capture subtle 3D flow physics due to their intrinsic limitation, but can reasonably predict time-averaged airfoil performance. Similarity and distinction between freestream vortex-LSB interaction in 2D and eddy-LSB interaction in 3D were explained. The role of the Kelvin-Helmholtz instability and Klebanoff modes in the transition of 3D airfoils were shown to be critical for understanding laminar-turbulent transition and LSB formation on airfoils in clean and turbulent freestreams.

Experimental study on the influence of inlet pressure on the flow parameters at the reverse flow boundary of vortex tube and its relation with energy separation performance

The performance of the vortex tube (also known as the Ranque-Hilsch tube) is significantly influenced by the inlet pressure and the internal flow processes, which can be categorized into two distinct flows: the center and the outer flow. The reverse flow boundary, acting as the interface between the two flows, plays a crucial role in energy separation. However, the research on the location of the reverse flow boundary and the associated flow parameters and energy transfer at this boundary are still insufficient. Therefore, this research employs the measurement method of interpolation probe to analyze the distribution of the reverse flow boundary and the specific flow parameters at this boundary. These parameters include static pressure gradient, angular velocity, and static temperature gradient, which substantially impact energy separation. A qualitative analysis of the energy transfer process in the vortex tube by examining mass flow in the reverse area and tangential and axial velocities are further delved. This analysis covers momentum transfer, heat transfer, and turbulent heat transfer processes resulting from the compression and expansion processes. The findings provide a research direction for exploring the energy separation performance of vortex tubes.

Effects of lubricant rheology and flow turbulence on the dynamic properties of hydrodynamic bearing

The present numerical investigation, which takes into account the coexisted consequence of the turbulence domain and non-Newtonian flow rheology, predicts the dynamic properties of a finite hydrodynamic bearing. Under the proper flow boundary conditions, the turbulence (linear) and non-Newtonian hypotheses have revised the Continuity and Navier–Stokes equation. The finite element approach is acclimated to calculate the clearance extent of the finite hydrodynamic bearing using the Galerkin approach and a sturdy iteration strategy. The cylindrical coordinate's versions of continuity and momentum equations are bestowed for the lubricant field of the finite hydrodynamic bearing, postulating the turbulence domain and non-Newtonian flow rheology. Dynamic properties of a finite hydrodynamic bearing have been enumerated using cross-coupled and direct lubricant-film bearing coefficients, whirl ratio and critical mass at various eccentricity ratios for different ranges of Reynolds numbers and different values of flow behavior index of the non-Newtonian rheological model. In comparison to perfect turbulence and non-Newtonian rheological status, the mixed regime demonstrated improved properties. The direct and cross-coupled lubricant-film bearing coefficients are improved by 12.28% and 20.85% respectively while critical mass is enhanced by 20.55% at selected values of Reynolds number and power-law index. The proposed severe fluid film domain improves the stability of finite hydrodynamic bearings significantly.

Dynamical low-rank approximation of the Vlasov–Poisson equation with piecewise linear spatial boundary

Dynamical low-rank approximation (DLRA) for the numerical simulation of Vlasov–Poisson equations is based on separation of space and velocity variables, as proposed in several recent works. The standard approach for the time integration in the DLRA model uses a splitting of the tangent space projector for the low-rank manifold according to the separated variables. It can also be modified to allow for rank-adaptivity. A less studied aspect is the incorporation of boundary conditions in the DLRA model. In this work, a variational formulation of the projector splitting is proposed which allows to handle inflow boundary conditions on spatial domains with piecewise linear boundary. Numerical experiments demonstrate the principle feasibility of this approach.

Dynamical response of the southwestern Laurentide Ice Sheet to rapid Bølling–Allerød warming

Abstract. The shift in climate that occurred between the Last Glacial Maximum (LGM) and the Early Holocene (ca. 18–12 kyr BP) displayed rates of temperature increase similar to present-day warming trends. The most rapid recorded changes in temperature occurred during the abrupt climate oscillations known as the Bølling–Allerød interstadial (14.7–12.9 kyr BP) and the Younger Dryas stadial (12.9–11.7 kyr BP). Reconstructing ice sheet dynamics during these climate oscillations provides the opportunity to assess long-term ice sheet evolution in reaction to a rapidly changing climate. Here, we use glacial geomorphological inversion methods (flowsets) to reconstruct the ice flow dynamics and the marginal retreat pattern of the southwestern sector of the Laurentide Ice Sheet (SWLIS). We combine our reconstruction with a recently compiled regional deglaciation chronology to depict ice flow dynamics that encompass the time period from pre-LGM to the Early Holocene. Our reconstruction portrays three macroscale reorganizations in the orientation and dynamics of ice streaming followed by regional deglaciation associated with rapid warming during the Bølling–Allerød interstadial. Initial westward flow is documented, likely associated with an early set of ice streams that formed during the advance to the LGM. During the LGM ice streaming displays a dominant north to south orientation. Ice sheet thinning at ∼15 ka is associated with a macroscale reorganization in ice stream flow, with a complex of ice streams recording south-eastward flow. A second macroscale reorganization in ice flow is then observed at ∼14 ka, in which southwestern ice flow is restricted to the Hay, Peace, Athabasca, and Churchill river lowlands. Rates of ice sheet retreat then slowed considerably during the Younger Dryas stadial; at this time, the ice margin was situated north of the Canadian Shield boundary and ice flow continued to be sourced from the northeast. Resulting from these changes in ice sheet dynamics, we recognize a three-part pattern of deglacial landform zonation within the SWLIS characterized by active ice margin recession, stagnation and downwasting punctuated by local surging (terrestrial ice sheet collapse): the outer deglacial zone contains large recessional moraines aligned with the direction of active ice margin retreat; the intermediate deglacial zone contains large regions of hummocky and stagnation terrain, in some areas crosscut by the signature of local surges, reflecting punctuated stagnation and downwasting; and the inner deglacial zone contains inset recessional moraines demarcating progressive regional ice margin retreat. We attribute these macroscale changes in ice flow geometry and associated deglacial behaviour to external climatic controls during the Bølling–Allerød and Younger Dryas but also recognize the role of internal (glaciological, lithological and topographic) controls in SWLIS dynamics.

Debris-flow avulsion tendency estimated from boreholes or channel cuts

Debris fans are an important landform for human habitation and development in mountainous regions. While much research in debris flow runout has focused on quantifying recurrence intervals and volume predictions, less research has focused on predicting their paths, which is important to reduce hazards to life and property. This project uses a quantitative stratigraphy tool, the compensation index, to evaluate the likelihood of debris path alteration by avulsion. The compensation index is a quantitative measure of the strength of compensational stacking, which measures the tendency of flow events to preferentially fill topographic lows. A significant constraint on the use of the compensation index for debris-flow hazard analysis is the need for large natural exposures of fan stratigraphy and the difficulty of mapping flow boundaries. Therefore, a more readily available proxy for estimating the compensation index of debris fans is needed, such as vertical stratigraphic sections from borehole data or natural channel incisions. Data from four previous studies was used to correlate trends in unit thickness to the compensation index calculated using data generated from simulated boreholes along the measured cross-sections. It was found that the central tendency and coefficient of variation of unit thickness had moderate correlation to the compensation index, and therefore predictive capacity, in both submarine and subaerial debris fans. By contrast, the same data for fluvial channelized floodplains revealed no correlation to compensation index. Therefore, a limited number of vertical stratigraphic sections could provide a reasonable estimate of the compensation index, and therefore avulsion tendency, of a debris-flow deposit.

Numerical study on influence of surface vegetation on aerodynamics of high-rise buildings

Vegetation on the surface of buildings is becoming popular in modern society. Besides their environmental and visual effect, vegetation can also affect the flow around them. In this study, the impact of surface planted trees on the aerodynamics of high-rise buildings under atmospheric boundary flow is investigated. Large eddy simulation is employed for the quantitative study. The results indicate that fluctuating lift force of building with both continuous ribs and fully planted trees can be 15.3% lower than that with ribs only. A closer examination on the model's surface pressure indicates that trees primarily decrease the mean and fluctuating wind pressures on the model's sidewalls by up to 20% and 26%, respectively. The mean flow field shows the separation bubble of the models with trees are significantly compressed and pushed away from the side face. Such flow structure contributes to the reduction of mean pressure coefficient at the vortex core of the separation bubble and side faces comparing to the reference model. Additionally, results also show that trees can significantly reduce flow fluctuations near the side faces of model, which is responsible of reduction of the fluctuating wind pressure coefficients on the model surfaces.

Rock Breaking Mechanism and Process Optimisation of Jet Cutting Basic Roof Rock under Submerged Jet Condition

The destruction of rock under the condition of a close submerged jet has become a hot topic of scientific research and engineering application in the past decade. With the unremitting efforts of a large number of experts and scholars around the world, gratifying progress has been made in the research of computational fluid dynamics (CFD) on the internal and external flow fields of the jet nozzle, the theoretical derivation of rock mechanics on the fracture initiation and propagation criteria of hydraulic fracturing, and the numerical simulation of jet erosion mechanism under the coupling of fluid and solid fields, however, for the rock mechanics hydraulic fracturing cutting engineering scale of non-oil drilling fracturing technology, the research on the fluid-solid coupling boundary conditions of fracturing fluid and hard dense rock under the flow state conditions of the submerged field inside and outside the borehole is not sufficient. In the calculation of the fluid-solid coupling boundary flow field under the non-submerged jet state, the control equation with Reynolds number between 2300-4000 shall be selected, while it belongs to the laminar flow state in the stage of hole sealing and pressurised fracturing. Therefore, Von-Mises equivalent plastic stress is selected in the mechanical model to calibrate the failure state of the rock-solid boundary, and the control equations of laminar flow and turbulent flow are selected to calibrate the fluid boundary. The mechanism of different stages of rock breaking by hydraulic fracturing jet can be further analysed in detail, and Comsol 6.0 multi-physical field simulation software is selected for verification. The research results will help deepen the understanding of rock breaking mechanism by jet and optimise the selection of parameters for field construction.

Poromechanics Modeling Forecasts Production, Analyzes Productivity Decline

_ This article, written by JPT Technology Editor Chris Carpenter, contains highlights of paper SPE 212200, “Accurate Production Forecasting and Productivity Decline Analysis Using Coupled Full-Field and Near-Wellbore Poromechanics Modeling,” by Yan Li, Bin Wang, and Jiehao Wang, Chevron, et al. The paper has not been peer reviewed. _ Productivity index (PI) decline is caused by different mechanisms in both the wellbore region and the far field. Damage in the wellbore region can be simulated by detailed wellbore modeling. A newly developed full-field and near-wellbore poromechanics coupling scheme is used in the complete paper to model PI degradation against time. Near-wellbore damage and field and well interactions are identified when applying the coupling scheme for a deepwater well. History matching, production forecasting, and safe drawdown limits are derived for operational decisions. Full-Field Reservoir Model and Near-Wellbore Poromechanics Model Coupling Full-field and near-wellbore modeling involves coupling the simulation of a full-field reservoir model with one or more near-wellbore poromechanics models. In the coupled simulation, the full-field reservoir model dictates the changing flow or thermal boundary conditions for the embedded near-wellbore models. Meanwhile, near-wellbore phenomena affect well productivity in the full-field reservoir model, altering the flow and thermal boundary conditions on all near-wellbore models in the same field. While capturing the dynamic interactions between the full-field model and all embedded near-wellbore models is of vital importance, the traditional near-wellbore modeling work flow considers only one-way coupling. This work flow is labor-intensive and lacks the dynamic interplay between the reservoir and near-wellbore models. A novel near-wellbore coupling framework is developed to automate data exchange and capture dynamic interactions between the full-field model and the near-wellbore model. In the full-field reservoir and geomechanics coupling applications, only a reservoir simulator solves reservoir flow and thermal equations and a geomechanics simulator solves solid mechanics equations. The reservoir and geomechanics simulators exchange 3D field data. Because solid mechanics equations are quasistatic, the geomechanics simulator does not take timesteps in full-field coupling. In near-wellbore coupling, however, both reservoir and geomechanics simulators solve reservoir flow and thermal equations. The near-wellbore geomechanics model also may solve solid mechanics or other equations coupled to flow or thermal equations in the near-wellbore model. All physics are included in the near-wellbore model for the coupling. Both field properties (3D) and transient boundary conditions (2D) must be mapped onto the near-wellbore model. It is important to ensure that the full-field and near-wellbore models are consistent. A 3D data-mapping module is developed to map flow and rock properties and initial conditions from the full-field model to the near-wellbore model. During the simulation, the transient pressure/temperature boundary conditions (2D) at the external boundary of the near-wellbore model must be mapped from the full-field model to the near-wellbore model to update transient boundary conditions. The automated data mapping in the coupling scheme eliminates the need for manual mapping of field properties and transient boundary conditions.

A semi-analytical pressure and rate transient analysis model for inner boundary and propped fractures exhibiting dynamic behavior under long-term production conditions

The loss of hydrocarbon production caused by the dynamic behavior of the inner boundary and propped fractures under long-term production conditions has been widely reported in recent studies. However, the quantitative relationships for the variations of the inner boundary and propped fractures have not been determined and incorporated in the semi-analytical models for the pressure and rate transient analysis. This work focuses on describing the variations of the inner boundary and propped fractures and capturing the typical characteristics from the pressure transient curves.A generalized semi-analytical model was developed to characterize the dynamic behavior of the inner boundary and propped fractures under long-term production conditions. The pressure-dependent length shrinkage coefficients, which quantify the length changes of the inner zone and propped fractures, are modified and incorporated into this multi-zone semi-analytical model. With simultaneous numerical iterations and numerical inversions in Laplace and real-time space, the transient solutions to pressure and rate behavior are determined in just a few seconds. The dynamic behavior of the inner boundary and propped fractures on transient pressure curves is divided into five periods: fracture bilinear flow (FR1), dynamic PFs flow (FR2), inner-area linear flow (FR3), dynamic inner boundary flow (FR4), and outer-area dominated linear flow (FR5). The early hump during FR2 period and a positive upward shift during FR4 period are captured on the log-log pressure transient curves, reflecting the dynamic behavior of the inner boundary and propped fractures during the long-term production period.The transient pressure behavior will exhibit greater positive upward trend and the flow rate will be lower with the shrinkage of the inner boundary. The pressure derivative curve will be upward earlier as the inner boundary shrinks more rapidly. The lower permeability caused by the closure of un-propped fractures in the inner zone results in greater upward in pressure derivative curves. If the permeability loss for the dynamic behavior of the inner boundary caused by the closure of un-propped fractures is neglected, the flow rate will be overestimated in the later production period.

Wake Formation Flow Physics and Boundary Layer Analysis on the Sides of the Isosceles Triangular Cylinder with Apex Pointing Downstream

Boundary layer interaction with downstream flow structures was numerically studied to find the region of inactivity behind an 75° isosceles triangular cylinder with apex pointing downstream at intermediate Reynolds numbers (Re = 520, 640, 840 and 1040). The Standard k-ε model in OpenFOAM was used in the study. Numerical results were validated against Particle Image Velocimetry data. Results revealed the stable region of inactivity characterized by low turbulent kinetic energy and vorticity. The onset of secondary vortex and separation point, independent of Reynolds number, was identified. The onset of the secondary vortex was located at (x = 2 mm) from the base and (y = 1.5 mm) from the apex on either side of the cylinder. The ratio of modulus of absolute primary z-component of vorticity, |ω_z^1 |, to the modulus of secondary z-component of vorticity, |ω_z^2 |, was found to be approximately equal to 1.2. This ratio is invariant of the Reynolds number of the study. These findings have practical implications. The unique properties of the inactivity region forms an ideal location that can be used for injecting fluid, placing measurement probe, active flow control and drag reduction. The research problem is formulated in the introduction. Literature is reviewed next providing the background. Details about the range of parameters, governing equations, numerical study details and software used are given in the methodology section. The results section gives the numerical results, verified by mesh refinement test and validated against experimental results. The results are finally discussed in the next section.

Energy transfer in the formation and evolution of energy-containing motions in turbulent boundary layers

This paper examines using direct numerical simulation (DNS) the recovery of the energy-containing scales of motion (energy-containing motions) when these are removed at the inflow boundary condition of a flat plate zero-pressure gradient turbulent boundary layer (ZPG-TBL). The energy-containing motions are identified from the spanwise spectral energy density of the interscale turbulent transport. Analysis of the resulting spanwise spectral energy densitiy of velocity and vorticity suggests that near-wall streaks are the initial energy-containing motions formed with quasi-streamwise vortices only developing further downstream. This result is in agreement with the previous study on turbulent channel flows by Kannadasan et al. [1].

Investigating the use of 3-component-2-dimensional particle image velocimetry fields as inflow boundary condition for the direct numerical simulation of turbulent channel flow

Direct numerical simulation (DNS) of turbulent wall-bounded flows requires long streamwise computational domains to establish the correct spatial evolution of large-scale structures with high fidelity. In contrast, experimental measurements can relatively easily capture large-scale structures but struggle to resolve the dissipative flow scales with high fidelity. One methodology to overcome the shortcomings of each approach is by incorporating experimental velocity field measurements into DNS as an inflow boundary condition. This hybrid approach combines the strengths of DNS and experimental measurements, allowing for a reduction in the streamwise computational domain and accelerated development of large-scale structures in turbulent wall-bounded flows. To this end, this paper reports the results of an investigation to establish the impact of limited spatial resolution and limited near-wall experimental inflow data on the DNS of a wall-bounded turbulent shear flow. Specifically, this study investigates the spatial extent required for the DNS of a turbulent channel flow to recover the turbulent velocity fluctuations and energy when experimental inflow data is typically unable to capture fluctuations down to the viscous sub-layer or the smallest viscous scales (i.e. the Kolmogorov scale or their surrogate viscous scale in wall-bounded turbulent shear slows) is used as the inflow to a DNS. A time-resolved numerically generated experimental field is constructed from a periodic channel flow DNS (PCH-DNS) at Reτ=\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Re_{\ au } =$$\\end{document} 550 and 2300, which is subsequently used as the inflow velocity field for an inflow–outflow boundary conditions DNS. The time-resolved experimental inflow field is generated by appropriately filtering the small scales from the PCH-DNS velocity by integrating over a spatial domain that is representative of a particle image velocimetry interrogation window. This study shows that the recovery of small scales requires a longer domain as the spatial resolution at the inflow decreases with all flow scales recovered and their correct scale-dependent energy is re-established once the flow has developed for 3 channel heights.

A hybrid data-driven framework for loss prediction of MCA airfoils

The Multi Circular Arc airfoils family has gained significant popularity in axial compressor design due to its capability to achieve higher pressure ratios with lower losses compared to conventional NACA airfoils. However, modeling their aerodynamic performance for transonic axial compressors remains a challenging task. This research addresses this issue by developing a data driven model for Multi Circular Arc airfoils loss correlation utilizing an automated blade-to-blade flow solver, StarCCM, under specific flow boundary conditions. Since the data built in this research includes high dimensionality, it is subjected to dimensionality reduction using Principal Component Analysis while retaining as much of its important information as possible. What sets our research apart is the comparative analysis of hybrid machine learning approaches. We specifically combine classification techniques based on shock wave formation on the airfoil’s suction side with Multi-Layer Perceptron neural network methods. This unique combination allows for a comprehensive prediction of MCA loss as the main objective of the research. The results demonstrate that this surrogate-based approach provides accurate loss model with 97% of precision error for the classification which results in 0.99 and 0.94 of R2 for loss model related to the transonic and subsonic airfoils loss prediction.

Measurement and comparison of transient thermal boundary layer of a Local-Heated Mini-Channel under sub- and super-critical conditions

This study is focused on experimentally measuring the transient thermal boundary layer following a sudden heat addition to turbulent flow boundary in supercritical mini-channels. Interferometric diagnostics were performed using a visualization channel to capture the time-evolution boundary layer. The measured pressure drop shows good agreement with numerical iterations according to the friction laws. Transient scenarios have been extracted from the raw interferograms to facilitate quantitative comparison of boundary conditions, such as bulk temperature and heat flux: (1) Localized heat accumulation and a slowly thermal response have been observed in subcritical state; (2) Relatively high density of 1.2 kg/m3 and low temperature changes of 0.008 K are observed in cases with relatively shorter critical distances to the critical point, attributed to the sharp increase in heat capacity and thermal expansion coefficients. (3) A relatively fast thermal convective field is established at larger heat fluxes due to the rapid density shifts of 1.2 kg/m3 and associated flow acceleration in the boundary layer.

Detonation structure in the presence of mixture stratification using reaction-resolved simulations

Many investigations of detonation-based combustors have identified reactant mixture inhomogeneity as having a leading-order impact on wave dynamics and combustion efficiency. To examine this phenomenon in a simplified context, an array of two- and three-dimensional channel detonation simulations are conducted in the present work. The reactant mixture consists of stratified fuel and air, wherein the randomly distributed equivalence ratio field features a characteristic stratification length scale. Detailed chemical kinetics are implemented in an adaptive mesh refinement solution framework where the region near the shock front is resolved with O(100) cells per representative ZND induction length. The results show that in comparison to baseline cases with uniform reactant mixtures, reactant stratification has a marked impact on the detonation structure. Increasing the stratification length scale increases the size and irregularity of the detonation cells, yielding larger variations in wave speed. Triple point collisions in fuel-rich regions lead to local wave speeds above the notional mean CJ speed, but wave passage through inert regions causes the local wave speed and strength to diminish. Conditional statistics show that increasing the stratification length scale increases the variance in pressure and temperature in the primary reaction zone, as well as the variance in heat release over a range of mixture conditions. In addition to the reactant mixture, the impact of the boundary condition behind the detonation is also investigated. The results show that an inflow boundary condition acts to over-drive the wave, leading to higher peak pressures, smaller detonation cells, and increased reactant consumption. On the other hand, cases with a wall behind the wave exhibit weaker waves with lower peak pressures and heat release rates, as well as greater variance in conditional quantities. Comparisons between complementary two- and three-dimensional simulations show reasonable qualitative agreement in wave structure, speed, and conditional statistics.Novelty and significance statement: The present work considers two- and three-dimensional simulations of detonations propagating through stratified reactant mixtures using detailed chemical kinetics with O(100) cells per representative ZND induction length. To the authors’ knowledge, it is the first study to analyze this type of configuration at this resolution and detail, as well as the first to compare two- and three-dimensional simulations of this configuration. The results have implications for the detonation wave behavior observed in practical propulsion and power-generating applications, such as rotating detonation engines.

Long-wave instabilities of sloping stratified exchange flows

We investigate the linear instability of two-layer stratified shear flows in a sloping two-dimensional channel, subject to non-zero longitudinal gravitational forces. We reveal three previously unknown instabilities, distinct from the well-known Kelvin–Helmholtz instability and Holmboe wave instability, in that they have longer wavelengths (of the order of 10 to $10^3$ shear-layer depths) and often slower growth rates. Importantly, they can grow in background flows with gradient Richardson number $\gg 1$ , which offers a new mechanism to sustain turbulence and mixing in strongly stratified flows. These instabilities are shown to be generic and relatively insensitive to Reynolds number, Prandtl number, base flow profile and boundary conditions. The nonlinear evolution of these instabilities is investigated through a forced direct numerical simulation, in which the background momentum and density are sustained. The growth of long unstable waves in background flows initially stable to short wave causes a decrease in the local gradient Richardson number. This leads to local nonlinear processes that result in small-scale overturns resembling Kelvin–Helmholtz billows. Our results establish a new energy exchange pathway, where the mean kinetic energy of a strongly stratified flow is extracted by primary unstable long waves and secondary short waves, and subsequently dissipated into internal energy.