Front Propagation Research Articles

Staircase resiliency in a fluctuating cellular array.

Inhomogeneous mixing by stationary convective cells set in a fixed array is a particularly simple route to layering. Layered profile structures, or staircases, have been observed in many systems, including drift-wave turbulence in magnetic confinement devices. The simplest type of staircase occurs in passive-scalar advection, due to the existence and interplay of two disparate timescales, the cell turn-over (τ_{H}), and the cell diffusion (τ_{D}) time. In this simple system, we study the resiliency of the staircase structure in the presence of global transverse shear and weak vortex scattering. The fixed cellular array is then generalized to a fluctuating vortex array in a series of numerical experiments. The focus is on regimes of low-modest effective Reynolds numbers, as found in magnetic fusion devices. By systematically perturbing the elements of the vortex array, we learn that staircases form and are resilient (although steps become less regular, due to cell mergers) over a broad range of Reynolds numbers. The criteria for resiliency are (a) τ_{D}≫τ_{H} and (b) a sufficiently high profile curvature (κ≥1.5). We learn that scalar concentration travels along regions of shear, thus staircase barriers form first, and scalar concentration "homogenizes" in vortices later. The scattering of vortices induces a lower effective speed of scalar concentration front propagation. The paths are those of the least time. We observe that if background diffusion is kept fixed, the cell geometric properties can be used to derive an approximation for the effective diffusivity of the scalar. The effective diffusivity of the fluctuating vortex array does not deviate significantly from that of the fixed cellular array.

Propagation of Combustion Front within Fractured Shale and Its Influence on Shale Structure and Crude Oil Properties: An Experimental Study

Summary The in-situ combustion (ISC) technique has emerged as a significant approach for shale oil production. However, currently, there is a lack of experimental evidence supporting the stable propagation of combustion front within fractured shale. This study aimed to investigate the combustion characteristics within fractured shale by using a self-designed combustion tube (CT) and an experimental scheme. Subsequently, an analysis of shale structure and oil properties was conducted. The results demonstrated that while the combustion front could propagate through shale with a single fracture width of approximately 43 μm, the combustion intensity gradually diminished, leading to an inability to sustain stable propagation in the later part of the oil-detritus mixtures. The combustion intensity within the shale was enhanced by preheating the shale at 250°C, resulting in an improved oil recovery from 67.8% to 77.9%. The findings indicated that the complex fractured shale allowed for the stable propagation of the combustion front without a significant decrease in combustion intensity. Moreover, the T2 spectrum analysis of shale revealed a gradual expansion of the pore-fracture structure and improved shale connectivity after combustion. The T1-T2 response illustrated the transformation of solid and heavy components into lighter components. Furthermore, the content of saturates and H in the oil increased after combustion, whereas there was a significant decrease in resins, O, and S. Overall, this study provided technical evidence supporting the feasibility of employing the ISC technique for the development of shale oil reservoirs with additional fractures.

Experimental and 2D fluid simulation of a streamer discharge in air over a water surface

The high reactivity and attractive properties of streamer discharges make them useful in many applications based on plasma-surface interactions. Therefore, understanding the mechanisms governing the propagation of a streamer discharge as well as its properties is an essential task. This paper presents the development and application of a 2D fluid model to the simulation of discharges triggered at the air-water interface by a pulsed nanosecond high voltage. Experimental characterization using 1 ns-time-resolved imaging reveals rapid transitions from a homogeneous disc to a ring and finally to dots during the discharge process. The simulation enables the determination of the spatio-temporal dynamics of the E-field and electron density, highlighting that the discharge reaches the liquid surface in less than 1 ns, triggering a radial surface discharge. As the discharge propagates along/over the water surface, a sheath forms behind its head. Furthermore, the simulation elucidates the transitions from disc to ring and from ring to dots. The former transition arises from the ionization front’s propagation speed, where an initial disc-like feature changes to a ring due to the decreasing E-field strength. The ring-to-dots transition results from the destabilization caused by radial electron avalanches as the discharge head reaches a radius of ∼1.5 mm. The simulation is further utilized to estimate a charge number and a charge content in the discharge head. This work contributes to a better understanding of discharge propagation in air near a dielectric surface, with the agreement between simulation and experiment validating the model in its present version.

Research of Ni—Al metal films using various physical methods

Ni—Al films obtained by resistive thermal evaporation using a "ВУП-5" vacuum station are studied. The main stages of film synthesis are established. The speed of propagation of the synthesis wave front and the radial speed of the islands are determined. Surface grain analysis and Fourier analysis are performed, providing information on the structural properties of the surface. The electrical conductivity of the thin metal film under study is determined. Keywords:Ni—Al thin metal film, resistive method, film synthesis, microstructure, coalescence, island structure, eddy current method

Redox cationic frontal polymerization: a rapid curing approach for carbon fiber-reinforced composites with high fiber content

Conventional frontal polymerization processes for epoxy-based composites rely on cations and radicals generated by a short (and local) light or heat stimulus in the presence of an iodonium salt and a radical thermal initiator. However, due to heat losses, the propagation of the exothermic curing front is often limited by sample geometry and filler concentration. Redox cationic frontal polymerization (RCFP) is a promising approach to radically expand the composition and design options of frontally cured epoxy-based composites. By adding stannous octoate as reducing agent, a higher number of radicals and cations are generated at lower temperature, which yields highly cured composite even at elevated filler content. In the current study, RCFP was used to cure standard unidirectional carbon fiber-reinforced composites based on a commercially available epoxy resin and the properties were compared with its anhydride hardener-cured counterpart. Cure degree and thermal properties of the resins were determined by ATR FT-IR spectroscopy and DMA analysis. Subsequently, unidirectional composites with a fiber volume content of ~ 60% were produced via vacuum infusion and subjected to DMA, tensile, compression, and inter-laminar shear tests. The results showed a remarkable similarity between mechanical properties of RCFP and anhydride hardener-cured composites. The RCFP-cured composites exhibited even a higher damping resistance and compression strength than anhydride hardener-cured composites. The results show that RCFP allows for a significant reduction in the curing time (from several hours to 60 min), while it yields composites with properties comparable to classic anhydride-cured systems.Graphical abstract

Orographic Controls on Extreme Precipitation Associated with a Mei-Yu Front

Abstract Taiwan regularly receives extreme rainfall due to seasonal mei-yu fronts that are modified by Taiwan’s complex topography. One such case occurred between 1 and 3 June 2017 when a mei-yu front contributed to flooding and landslides from over 600 mm of rainfall in 12 h near the Taipei basin, and over 1500 mm of rainfall in 2 days near the Central Mountain Range (CMR). This mei-yu event is simulated using the Weather Research and Forecasting (WRF) Model with halved terrain as a sensitivity test to investigate the orographic mechanisms that modify the intensity, duration, and location of extreme rainfall. The reduction in WRF terrain height produced a decrease in rainfall duration and accumulation in northern Taiwan and a decrease in rainfall duration, intensity, and accumulation over the CMR. The reductions in northern Taiwan are linked to a weaker orographic barrier jet resulting from a lowered terrain height. The reductions in rainfall intensity and duration over the CMR are partially explained by a lack of orographic enhancements to mei-yu frontal convergence near the terrain. A prominent feature missing with the reduced terrain is a redirection of postfrontal westerly winds attributed to orographic deformation, i.e., the redirection of flow due to upstream topography. Orographically deforming winds converge with prefrontal flow to maintain the mei-yu front. In both regions, the decrease in mei-yu front propagation speed is linked to increased rainfall duration. These orographic features will be further explored using observations captured during the 2022 Prediction of Rainfall Extremes Campaign in the Pacific (PRECIP) field campaign. Significance Statement This study examines the impact of terrain on rainfall intensity, duration, and location. A mei-yu front, an East Asian weather front known for producing heavy, long-lasting rainfall, was simulated for an extreme rain event in Taiwan with mountain heights halved as a sensitivity test. Reducing terrain decreased rainfall duration in northern and central Taiwan. Decreases in rainfall duration for both regions is attributed to increased mei-yu front propagation speed. This increase in northern Taiwan is attributed to a weakened barrier jet, a low-level jet induced by flow blocked by the steep mountains of Taiwan. A unique finding of this work is a change in winds north of the front controlling movement of the front near the mountains in central Taiwan.

Propagation of extended fractures by local nucleation and rapid transverse expansion of crack-front distortion.

Fractures are ubiquitous and can lead to the catastrophic material failure of materials. Although fracturing in a two-dimensional plane is well understood, all fractures are extended in and propagate through three-dimensional space. Moreover, their behaviour is complex. Here we show that the forward propagation of a fracture front occurs through an initial rupture, nucleated at some localized position, followed by a very rapid transverse expansion at velocities as high as the Rayleigh-wave speed. We study fracturing in a circular geometry that achieves an uninterrupted extended fracture front and use a fluid to control the loading conditions that determine the amplitude of the forward jump. We find that this amplitude correlates with the transverse velocity. Dynamic rupture simulations capture the observations for only a high transverse velocity. These results highlight the importance of transverse dynamics in the forward propagation of an extended fracture.

OME[formula omitted]/air mixtures: The structure and propagation mechanism of their laminar premixed flames

Oxymethylene dimethyl ethers (CH3O(CH2O)xCH3 or OMEx) are synthetic substitutes of diesel and jet fuel. They can achieve soot-less combustion and have been shown to exhibit favourable ignition characteristics. With the aim to explore their combustion physics, the structure of their flames is analysed using the algorithmic tools of the Computational Singular Perturbation (CSP) method. It is demonstrated that the flame structure and the mechanism for its propagation is encapsulated in the fastest explosive mode, introduced by chemical dynamics in a narrow region of the flame front. This narrow region extends around the point where the maximum heat release rate (MHRR) is recorded. CSP diagnostics identify (i) reactions H + O2→ O + OH and CO + OH → CO2 + H beyond MHRR as initiating an upstream diffusion of heat and H radicals, both of which are generated by the exothermic reaction H2 + OH → H2O + H, (ii) a downstream convective motion that ensues and (iii) chemical activity that is initiated ahead of MHRR and sustained by the convective transport of fuel and oxygen. CSP also identifies differences in the action of the upstream diffusion of heat and H radicals and of the downstream transport of fuel and oxidiser. It is also shown that the fastest explosive mode far from the neighbourhood of MHRR has no influence on the structure and propagation of the flame front. Although OME2−4/air mixtures exhibit drastically different dynamics in the context of homogeneous autoignition, very small differences are recorded in terms of flame structure and propagation. The findings suggests that, since practically all major processes driving the flame front propagation are similar to those of hydrocarbons, OMEx fuels can be used as “drop-in” fuels in currently employed burners.

Flame evolution and pressure dynamics of premixed stoichiometric ammonia/hydrogen/air in a closed duct

Ammonia is regarded as a promising energy carrier due to its zero-carbon emissions and its suitability for long-distance, large-scale storage, and transportation. Ammonia/hydrogen mixed combustion is an important way to solve the problem of high ignition temperature and low flame speed in the process of ammonia combustion. This study systematically investigates the evolution of flames and pressure dynamics in a closed duct for ammonia/hydrogen/air mixtures with varying hydrogen blending ratio (0–100 %). Simultaneously, the study assesses the impact of flame instability and heat release on the evolution of ammonia/hydrogen/air flames and pressure dynamics. The results show that the flame tip speed and overpressure increase with the increase of hydrogen blending ratio and initial pressure, and when the hydrogen blending ratio is more than 25 %, the growth range of flame tip speed peak and overpressure peak value increases significantly. At the same time, the flame thickness decreases after hydrogen blending, the hydrodynamic instability increases, and the degree of wrinkles on the flame surface increases, which promotes flame acceleration. The Froude number increases significantly, the influence of buoyancy instability gradually weakens, and the flame core does not undergo significant displacement. In the later stage of flame propagation, both the flame front propagation velocity and pressure show significant pulsation phenomena. The heat release calculation results show that hydrogen blending significantly increases the adiabatic flame temperature. Consistent with the evolution law of overpressure, when the hydrogen blending ratio is greater than 25 %, the net heat release increases significantly. The elementary reaction with the highest heat release rate changes from the chain termination reaction R41 to the chain termination reaction R3. The research results can provide theoretical basis and experimental data support for the formulation of safety standards during the regulation and utilization of the combustion structure of ammonia/hydrogen/air mixtures.

Development and assessment of a three-field subchannel code for thermo-hydraulic transient analysis under reflood conditions

A three-field twelve-equation subchannel code was developed that contains separate energy equations for the continuous liquid and droplet phases. Several physical models were modified and optimized to accommodate the code and improve predictions of reflooding. Eight RBHT tests were used to assess the code. Comparisons show that the predictions of the code are in good agreement with experimental data on the quench front propagation, peak cladding temperature, and temperatures of the cladding, vapor and spacer grid. Moreover, the code can capture trends in droplet size and velocity and obtain the continuous liquid-droplet temperature difference. The impact of new interfacial heat transfer model, droplet impingement heat transfer model and energy equations on the reflooding calculation are analyzed. Overall, the code is applicable to the simulation of the reflooding phenomenon, and the results can provide guidance for model modification. Meanwhile, the pre-assessment of droplet temperature provides new perspectives for the code development.

An efficient deep learning-based workflow for CO2 plume imaging considering model uncertainties with distributed pressure and temperature measurements

Monitoring CO2 plumes throughout the operation of geologic CO2 sequestration projects is essential to environmental safety. The evolution of underground CO2 saturation can be predicted using high-fidelity numerical simulations. However, high-fidelity simulations can be prohibitively expensive to compute. As a result of recent developments in data-driven models, rapid predictions of the CO2 plume can now be made using readily available pressure and temperature measurements. This study presents a novel deep learning-based workflow for efficiently visualizing CO2 plumes in near real-time while considering their uncertainties.In our deep learning workflow, we visualize the CO2 plume images in the reservoir as a propagating saturation front, represented by ‘onset time’, using field measurements, including downhole pressure and temperature. At a given location, the ‘onset time’ is the calendar time when the CO2 saturation exceeds a certain threshold. Therefore, a single image of the CO2 front propagation is captured using the ‘onset time’ rather than storing multiple CO2 saturation images at different time steps. The use of ‘onset time’ significantly reduces memory and computational cost of deep learning-based framework, enabling large-scale field applications.We use a variational autoencoder-decoder (VAE) network to compress high dimensional ‘onset time’ images into low dimensional latent variables while considering uncertainties of the predicted images. The use of VAE and onset time, simplifies the overall neural network architecture and significantly enhances the training efficiency. To estimate the latent variables of the VAE network, we train a feed forward neural network model that incorporates available monitoring data such as downhole pressure and temperature measurements. The estimated latent variables are then fed into a trained decoder network to generate 3D onset time images, visualizing the propagation of CO2 plume in near real time.The proposed workflow is applied to both synthetic and field cases, where the field application is a large-scale geological carbon storage project in a carbonate reef reservoir in the Northern Niagaran Pinnacle Reef Trend in Michigan, USA. The field measurements include distributed temperature sensing (DTS) and distributed pressure responses from down-hole gauges along the monitoring well. The predicted CO2 plume images provided by the proposed workflow are shown to be consistent with the simulation results of traditional history matched model using genetic algorithm. Our deep learning-based framework can predict the CO2 front propagation in terms of onset time in seconds, making it well-suited for real time decision-making and operational optimization.

Analytical modeling of phase change in a composite wall comprising two distinct phase change materials in series

Most of the past literature on solid-liquid phase change heat transfer modeling focuses on a single phase change material (PCM), whereas, a sandwich of two or more PCMs arranged in series may be of interest in several applications. The modeling of heat transfer and phase change in such a system is complicated by the presence of multiple phase change fronts propagating at the same time. This work presents a theoretical analysis of the problem of melting of a two-PCM stack being heated from one end. Depending on whether the lower-melting PCM is located next to or away from the heat source, two distinct cases are considered. The propagation of the phase change fronts in each case is divided into several stages, each of which is characterized by a distinct phase change and/or sensible heating processes. The transient temperature field in each stage is determined in the form of infinite series solutions, and the propagation of phase change front over time is then determined using energy conservation at the phase change interface. In this manner, explicit expressions are derived for the time taken for each layer to completely melt. Results are shown to be in good agreement with finite-element simulations, and with the exact Stefan solution under special conditions. Key non-dimensional parameters that influence the nature of the two-PCM phase change process are identified, and their impact on phase change performance is analyzed. Quantitative predictions are presented for the conditions under which both PCMs finish melting at the same time, which is a favorable outcome for efficient energy storage.

Numerical Simulation of integrated three-dimensional hydraulic fracture propagation and proppant transport in multi-well pad fracturing

Multi-well pad fracturing is an indispensable technology for the sustainable development of unconventional oil and gas reservoirs, but the integrated simulation of fracture propagation and proppant migration in massive fracturing treatments remains challenging, including complicated numerical implementations and computational bottlenecks due to its multi-physics, fracture front propagation, size differences between proppant particles and hydraulic fractures. In this work, an integrated model of fracture propagation and proppant migration in multi-well pad fracturing was established based on the three-dimensional discontinuous displacement method for fracturing deformation computations and the Eulerian-Eulerian framework for capturing proppant concentration. The proposed simulator not only accounts for key physical processes such as stress shadow, dynamic flow rate distribution in the wellbore, fluid leak-off, perforation erosion, fractures evolution, proppant bridging, and gravitational settling, but also couples the simulation of 3D planar propagating fractures and proppant transport at well-pad scale, and presents excellent computational efficiency and promising applications than conventional fracturing models. To explore the fracture geometries and proppant distribution in typically stacked layers, single-well, multi-lateral wells, and multi-well pad fracturing are simulated from simple to complex. The findings showed that fracture propagation is jointly affected by intra-stage, inter-stage, and inter-well stress interference mechanisms, and fracture geometries are characterized by unevenness, asymmetry, and heel bias. Proppant transport is controlled by the fracture width and geometries, and most of the proppant is concentrated near the wellbore area and bridging at the narrow fractures. Fracture geometries and proppant distribution are sensitive to the landing locations of horizontal wells because hydraulic fractures always grow in the path of least resistance. Although close lateral well-spacing is beneficial to improve the stimulation efficiency, there are some challenges such as overlapped SRV, enhanced inter-well interference, and severe fracturing hits. Compared with the layout of stacked wells, staggered wells are not only beneficial to reduce the vertical inter-well interference, but also improve reservoir stimulation efficiency. Overall, the stimulated area of the high in-situ stress layers is relatively low even with intentionally landed wells.

Strontium concentration spikes as evidence for fluid-driven carbonate dissolution during subduction

Carbonate dissolution, the process by which carbonate minerals are stoichiometrically dissolved by an aqueous fluid, likely plays an important role in liberating carbon from subducting slabs and may thus make significant contributions to global carbon fluxes. It is therefore necessary to understand the fluid infiltration and flow geometries associated with this process. Despite its widespread use as a powerful tracer of fluid-rock interaction, strontium (Sr) has yet to be widely applied to carbonate dissolution scenarios. Previous work has identified a partially altered metacarbonate layer from Syros, Greece, that has undergone significant carbonate dissolution leading to ∼90% CO2 loss (Ague and Nicolescu, 2014). A bulk-rock compositional profile along this layer has a distinctive spike-shaped Sr profile. Sr concentrations are uniform in the altered region behind the reaction front, increase to form a spike just ahead of the front, and taper off to a baseline value moving farther into the unaltered metacarbonate rock. The origin of this type of spike shape and its potential utility for better understanding carbonate dissolution processes remain unaddressed. Here, we present a numerical diffusion-reaction model which demonstrates that both carbonate dissolution and multidimensional transport of Sr, in which the primary advective flow direction is perpendicular to the propagation of the diffusion front, are required to produce this shape. This model, as well as newly acquired Sr concentration data for calcite (former aragonite) and epidote, also indicate that the infiltration of a fluid with an elevated Sr concentration is necessary. We suggest that such spike-shaped features may aid in identifying carbonate dissolution in other settings, such as vein-selvage systems in metacarbonate rocks or reaction rinds on carbonate-bearing mélange blocks.

Observations of river-wave interactions at a small-scale river mouth

Small coastal rivers often discharge directly into the surfzone, where the fate of freshwater and river-borne materials (e.g., sediment, nutrients and contaminants) is primarily determined by interactions between buoyant plume processes and the nearshore wave forcing. Understanding river-wave interactions in the surfzone is essential to assess coastal water quality and impacts on ecosystem health along adjacent shorelines. These interactions are modulated by the variability in river discharge, offshore wave climate, tidal forcing, and surfzone dynamics, and therefore occur over a wide range of timescales. However, the different timescales at which river-wave interactions occur have not been fully investigated. Here we use in-situ and remote (drone-based) observations collected at the mouth of the Maipo River, a small-scale river system in central Chile, to investigate interactions between the river outflow and the incident wave forcing that may influence the distribution of freshwater along the coast. We focus on interactions occurring at infragravity (wave groups), tidal, and synoptic (offshore wave climate) timescales. The observations included inlet and surfzone deployments, and documented low river discharge conditions (QR∼ 20–27 m3s−1) with variable wave forcing (Hs∼ 1–3 m). We observed that the salinity signal at the river mouth is strongly modulated by the tide, with a freshwater plume forming solely during the ebb. Large oscillations associated to infragravity (IG) motions are continuously observed on top of the tidal variability. These oscillations are evidenced in surfzone salinity, inlet water levels, and inlet velocity during the early ebb when the freshwater plume starts to develop. As the plume evolves into a strong outflow jet during the late ebb, they are not observed in the inlet and become restricted to the surfzone. The combined analysis of drone imagery and in-situ observations suggests that the observed variability in surfzone salinity at IG frequencies is associated with the arrival of wave groups and the propagation of wave fronts in the plume area. In particular, during the late ebb, the large oscillations in surfzone salinity (amplitude of 15–20 psu) are explained by a contraction and expansion of the outflow plume in response to the variability in the onshore wave forcing at the timescale of wave groups. On synoptic timescales (hours to days), the surfzone salinity responded to the intensity of the offshore wave forcing (∼ 20 m depth), with higher waves promoting lower salinity as freshwater is mostly retained near the coast. A wave-river momentum comparison confirms that the wave forcing had the potential to trap freshwater in the surfzone. This study shows that nearshore waves may influence the dynamics of small river plumes and the fate of freshwater (and associated terrestrial materials) on a wide range of timescales, from infragravity (1–5 min) to hours and days.

Front instabilities in the presence of convection due to thermal and compositional gradients.

Reaction fronts separate fluids of different densities due to thermal and compositional gradients that may lead to convection. The stability of convectionless flat fronts propagating in the vertical direction depends not only on fluid properties but also in the dynamics of a front evolution equation. In this work, we analyze fronts described by the Kuramoto-Sivashinsky (KS) equation coupled to hydrodynamics. Without density gradients, the KS equation has a flat front solution that is unstable to perturbations of long wavelengths. Buoyancy enhances this instability if a fluid of lower density is underneath a denser fluid. In the reverse situation, with the denser fluid underneath, the front can be stabilized with appropriate thermal and compositional gradients. However, in this situation, a different instability develops for large enough thermal gradients. We also solve numerically the nonlinear KS equation coupled to the Navier-Stokes equations to analyze the front propagation in two-dimensional rectangular domains. As convection takes place, the reaction front curves, increasing its velocity.

ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ФРОНТА ПЛАМЕНИ ГАЗОВОЗДУШНОЙ СТРАТИФИЦИРОВАННОЙ СМЕСИ

The study conducts a test and numerical modeling of fl ame front propagation along a stratifi ed gas mixture in a small transverse fl at channel made by two parallel plates. Images of the fl ame front were obtained through direct photography. A shape change was detected during the fl ame oscillatory propagation along a long channel. The distribution fi elds of temperature, gas concentration, pressure, velocity, and streamline were constructed. The simulation fi ndings have been shown to qualitatively coincide with the test fi ndings.

RESEARCH OF THE UNSTEADY MOVEMENT OF LIQUID IN OPEN CHANNELS

Summary. In this work, the transient processes were studied and quantitative calculations were performed on the example of irrigation canals located in the south of Ukraine. Field, laboratory, and mathematical studies of the hydraulic structures' throughput capacity were carried out. Transient processes that arise in the process of re-regulation of the water consumption regime or as a result of emergencies, which are important aspects of the functioning of irrigation canals, were investigated. The basis of such processes is the unsteady fluid motion. To study the peculiarities of unsteady motion in irrigation canals, the general characteristics of the filling wave were determined, the instantaneous profiles of the free surface of the flow and changes in water levels over time in fixed sections along the canal were analyzed, and theoretical calculations were compared with the results of field observations. Changes in the average velocity of the filling wave front propagation along the channel were also determined, and the conditions for the transformation of the maximum flow along the entire length of the channel were studied. The algorithm for solving the problem uses a previously developed model based on the numerical integration of the Saint-Venant equations. Numerical calculations of transient processes were performed using a program that made it possible to calculate optimal control options, evaluate the effectiveness of the chosen method of automated water supply regulation, and calculate typical operating modes of the canal. As an example, to solve the problem of water distribution management, various schemes of automatic regulation of water flows and levels on the North Crimean Canal were considered. The results of the calculations were experimentally verified in real conditions with the same initial data. The experiments were carried out at the facilities located on the Razdolnensky Canal. The results of field observations were modeled and tested in the hydraulic tray of the department's laboratory. The field and laboratory studies confirmed the results of preliminary calculations. Keywords: hydraulic structures, unsteady motion, maximum flow transformation, water resources, irrigation systems.

Relationship between spatially heterogeneous reaction dynamics and photochemical kinetics in single crystals of anthracene derivatives.

Understanding physicochemical property changes based on reaction kinetics is required to design materials exhibiting desired functions at arbitrary timings. In this work, we investigated the photodimerization of anthracene derivatives in single crystals. Single crystals of 9-cyanoanthracene (9CA) and 9-anthraldehyde (9AA) exhibited reaction front propagation on the optical length scale, while 9-methylanthracene and 9-acetylanthracene crystals underwent spatially homogeneous conversion. Moreover, the sigmoidal behavior in the absorbance change associated with the reaction was much pronounced in the case of 9CA and 9AA and correlated with the observation of heterogeneous reaction progress. A kinetic analysis based on the Finke-Watzky model showed that the effective quantum yield of the photochemical reaction changes by more than an order of magnitude during the course of the reaction in 9CA and 9AA. Both the reaction front propagation and nonlinear kinetic behavior could be rationalized in terms of the difference in the cooperativity of the reactions. We propose a plausible mechanism for the heterogeneous reaction progress in single crystals that depends on the magnitude of the conformational change required for reaction. Our results provide useful information to understand the connection between photochemical reaction progress in the crystalline phase and the dynamic changes in the physicochemical properties.

Accelerated ignition-shock coupling and deflagration to detonation transition by ozone kinetic enhancement of dimethyl ether mixture

This study aims to understand the effect of kinetic enhancement by ozone addition on the Deflagration to Detonation Transition (DDT) and ignition-shock coupling in a microchannel using dimethyl ether (DME). Simultaneous high-speed chemiluminescence and shadowgraph imaging are carried out to analyze the flame front propagation, shock wave structures, autoignition, and ignition-shock coupling during the DDT process. The experimental results show that for all cases, with and without ozone kinetic enhancement, autoignition and DDT always occur in a reactivity gradient region with oblique shock cluster between the flame front and the flame acceleration-induced precursor shock. The observed DDT occurs via the Zeldovich hot spot gradient mechanism through the following stages: (1) Oblique shock cluster formation and acoustic compression between the flame front and the precursor shock, (2) autoignition initiation in the oblique shock cluster region, (3) spontaneous ignition wave propagation, and (4) ignition-shock coupling. It is found that autoignition location, oblique shock strength, induction zone length, ignition-shock coupling, and DDT onset flame velocity are very different from ozone kinetic enhancement. As the ozone addition is higher, the autoignition is advanced and the ignition location moves farther ahead of the flame front at a lower flame front propagation speed and weaker oblique shock clusters. Two different ignition-shock coupling modes are identified depending on ozone concentration. Without ozone kinetic enhancement, DDT occurs via ignition wave coupling with the precursor shock. However, ozone addition can accelerate the autoignition and initiate a direct ignition-shock coupling for DDT without interaction with the precursor shock. The present results show that kinetic enhancement in autoignition using an ignition enhancer such as ozone or plasma is critical and effective in accelerating and controlling DDT.