Flow Rate Distribution Research Articles

Multi-scale modeling of aerosol transport in a mouth-to-truncated bronchial tree system

Computational fluid particle dynamics (CFPD) is widely employed to predict aerosol transport in a truncated bronchial tree model on account of its capacity to reveal details of flow field and particle movement. However, setting a physiologically consistent boundary condition in the CFPD for the idealized or image-based truncated bronchial tree model is still a challenge. This paper proposes a multi-scale modeling method, which contains an Extend-Bronchial tree-Network (EBN) boundary condition for a mouth-to-truncated bronchi system. The comparison between EBN boundary condition and a commonly used uniform pressure (UP) boundary condition is conducted. Subsequently, EBN method is used to study the nano-micron (100nm-10μm) particles transport in the mouth-to-truncated bronchi model at different inhalation volume rates (15, 60, 90L/min). Results show that EBN method is more physiologically rational and two methods differ in flow distribution in lobes, vortex structure, and particle transport. The maximum difference in flow rate distribution in lobes between two methods is about 20%, while the maximum relative disparity of particle penetration fraction from lobes and deposition fraction in the TLB is about 93% and 30%, respectively. Meanwhile, this paper reveals the variation of deposition fraction and penetration fraction with the changes in particle diameter and inhalation volume. Deposition efficiency, deposition hotspots and deposition mechanism are also analyzed with inlet Stokes number (Stk) and Reynolds number (Re). This research establishes a foundation for the simulation of aerosol transport in a whole respiratory tract and provides references for inhalation drug delivery and air pollutant management.

A method to calculate the two-phase distribution in a microchannel heat exchanger

Maldistribution seriously impacts the heat transfer performance of the microchannel heat exchanger (MCHX). Fully understanding the two-phase distribution in the heat exchanger is important for advancing academic research and engineering applications. This study introduces a novel method for quantifying two-phase distribution in a microchannel heat exchanger. An experimental setup was developed to measure the local vapor mass fraction in the heat exchanger header. Capacitance signals were measured under inlet vapor mass fractions from 0 to 1, and inlet flow rates of 10, 15, and 25 g s−1 corresponding to a mass flux of 17.47, 26.2, and 43.66 kg m−2s−1, respectively. The local vapor mass fraction in the header was estimated using the capacitance measurements. The mass flow rate in the header, the microchannel tube, and the vapor mass fraction in the tube were calculated using the proposed model. The calculation model was validated against literature data, and the results were analyzed. The analysis reveals the characteristics of vapor mass fraction and mass flow rate distribution in the MCHX and further elaborates on the effects of phase separation, entrainment ratio, and pressure drop balance on the distribution. The proposed method can evaluate distribution in the header and tubes of microchannel heat exchangers, and it is also applicable to other types of two-phase flow devices.

A Comprehensive Flow–Mass–Thermal–Electrochemical Coupling Model for a VRFB Stack and Its Application in a Stack Temperature Control Strategy

In this work, a comprehensive multi-physics electrochemical hybrid stack model is developed for a vanadium redox flow battery (VRFB) stack considering electrolyte flow, mass transport, electrochemical reactions, shunt currents, and as heat generation and transfer simultaneously. Compared with other VRFB stack models, this model is more comprehensive in considering the influence of multiple factors. Based on the established model, the electrolyte flow rate distribution across cells in the stack is investigated. The distribution and variation in shunt currents, single-cell current and single-cell voltage are analyzed. The distribution and variation in temperature and heat generation and heat transfer are also researched. It can be found that the VRFB stack temperature will exceed 40 °C when operating at 60 A and 100 mA cm−2 at an ambient temperature of 30 °C, which will lead to electrolyte ion precipitation, affecting the performance and safety of the battery. To control the stack temperature below 40 °C, a new tank cooling control strategy is proposed, and the suitable starting cooling point and the controlled temperature are specified. Compared with the common room cooling strategy, the new tank cooling strategy reduces energy consumption by 27.18% during 20 charge–discharge cycles.

CFD–DEM modeling of particle segregation behavior in a simulated flash smelting furnace

Particle behavior plays a crucial role in flash smelting furnaces, but limited understanding is hindering further development of this modern pyrometallurgical process. This study presents a two-way coupled CFD–DEM approach to investigate particle distribution, collision, and segregation behaviors within a furnace. The simulation results are validated against experimental data, demonstrating good agreement. The effects of the particle mass–flow rate and size distribution are then analyzed. The results indicate that small particles are carried by airflows to the center of the particle-dense region, whereas large particles remain at the outer part, worsening segregation. Increasing the mass–flow rate significantly increases the particle number density and number of particle collisions but barely influences size-induced segregation. In contrast, enlarging the size of fed particles reduces particle number density and broadens dense distribution within a safe range. This ensures the necessary space and high mixing index of the particles, which benefits high-rate smelting processes.

Transient modeling of stratified thermal storage tanks: Comparison of 1D models and the Advanced Flowrate Distribution method

Thermal energy storage (TES) is one of the key technologies for enabling a higher deployment of renewable energy. In this context, the present study analyzes the modeling strategies of one of the most common TES systems: stratified thermal storage tanks. These systems are essential to many solar thermal installations and heat pumps, among other clean energy technologies. Three different one-dimensional tank models are compared by their computing speed and resilience to long time steps. Two of the models analyzed are numerical, one being explicit and the other one implicit, and the other is analytical. The models are validated against data from experiments carried out considering small-scale stratified tanks, showing that their performance can be improved by using the Advanced Flowrate Distribution (AFD) method. The results show that the analytical model maintains its accuracy with longer time steps and is robust against divergence. Conversely, the numerical models show equivalent performance for short time steps, while the computation time is reduced. Although the AFD method shows promising results by achieving an improvement of 43% in terms of Dynamic Time Warping, its parameter optimization must be generalized for different tank designs, flow rates, and temperatures.

Multi-objective design optimization and physics-based sensitivity analysis of field emission electric propulsion for CubeSat platforms

Field-emission electric propulsion is an electrostatic space electric propulsion technology that offers various advantageous features including efficient design, high specific impulse, and versatile thrust capabilities ranging from micro-Newton to milli-Newton levels. These characteristics make this type of propulsion a promising technology for small satellite platforms, enabling precise attitude control, orbit maintenance, and de-orbiting through ionization and acceleration of a liquid metal propellant. The growing demand for small propulsion systems in CubeSat platforms has spurred significant progress in modeling and characterizing field emission electric propulsion thrusters to enhance their overall performance. However, little study has been conducted to investigate the effect of geometric configurations on electric fields or expelled ion trajectories for design optimization. In this study, multi-objective design optimization is performed by incorporating electrostatic simulation coupled with an analytical performance model into evolutionary algorithms based on prediction from surrogate modeling, aiming to optimize the thruster emission design to maximize thruster performance. Physical insights into the key design factors influencing the performance of field emission electric propulsion have been gained by probing into the interaction between ion particles and electric field behavior within the thruster. It has been found that the length of the emitter tip has a significant effect on plume divergence, i.e., a longer emitter tip under the influence of electric field at higher emitter current tends to result in lower initial acceleration of emitted ions and subsequently wider spread or divergence of the ion beam. A shorter emitter tip, on the other hand, generates a sharper E-field gradient, resulting in a more focused and narrower ion beam. Additionally, sensitivity analysis has identified the mass flow rate and potential distributions as the most influential design factors on performance due to the active roles they play in the performance generation process.

Effects of slit geometric parameters on spray characteristics of double-slit pintle injectors

Pintle injectors have garnered significant research attention in recent years, particularly for their applicability in reusable launch vehicles, owing to their wide thrust control range and excellent combustion stability. While research has explored the characteristics of pintle injectors in the context of developing these components for actual engine applications, studies focusing on the effects of design parameters on injector performance have been limited. This study investigated the effects of slit geometric parameters, specifically the blockage factor (BF) and slit area ratio (γ), on the spray characteristics of double-slit pintle injectors. Cold-flow tests were conducted using planar pintle injectors with water and ethanol as simulants. The spray angle and Sauter mean diameter (SMD) were measured using the shadowgraph technique, and the distribution of mass flow rate and mixture ratio was analyzed using a mechanical patternator. The experimental results revealed that two distinct streams were injected at different angles from each row of slits, resulting in a division of spray shape, SMD, and mass flow distribution into three regions based on the two streams. These spray angles, termed primary and secondary spray angles, were quantified as functions of the local momentum ratio, determined by BF and γ. To correlate the spray characteristics with combustion performance, mixing quality and a representative droplet size metric, the integral Sauter mean diameter (ID32), were presented. The study found that higher values of BF and γ corresponded to improved mixing quality.

A novel solar tower simulator for hydrogen regeneration by thermo-catalytic cracking of Ammonia: Energy, Exergy, and CFD investigations

This paper presents a novel concept of a central solar tower-based hydrogen regeneration system by thermo-catalytic cracking of ammonia. This includes a closed volumetric receiver for ammonia heating and cracking, a 2.45 kWe high-concentration radiation source, and an air-cooling system for the radiation source. The thermodynamic analysis shows that the net solar-to-hydrogen generation efficiency, first law efficiency, and exergy efficiency are 31 %, 82 %, and 79.8 %, respectively, for the optical efficiency of 35 %, the conversion efficiency of 99 %, and the waste heat recovery fraction of 20 %. The simulations with an ammonia flow rate of 1 kg/hr and flux concentration of 1200 Suns show the regeneration of 1.27 kg H2/day at Jodhpur and 1.57 kg H2/day at Ladakh. Monte Carlo ray tracing simulations revealed that the developed radiation source provides an average flux of ∼ 2 MW/m2 onto a spot size of 20 mm with an efficiency of about 31 %. Reynolds-averaged Navier Stokes-based analyses are used to design a novel air-cooling system for the radiation source. The computations show that (a) an unequal distribution of air mass flow rates at the front mḟ and back mḃ, should be preferred with mḟ>mḃ, and (b) for a total mass flowrate of 0.1 kg/s with mḟmḃ = 4, the maximum surface temperature of reflectors and Xe-arc lamps is less than 315 K and 610 K, respectively. The investigations revealed the feasibility of the concept and provided valuable insight for the planned experimental assessment and scale-up of the solar tower simulator.

Effect of liquid-vapor two-phase flow maldistribution on the thermal performance of brazed plate heat exchangers

The effect of liquid-vapor two-phase flow maldistribution on the thermal performance of brazed plate heat exchangers (BPHEs) is investigated experimentally and numerically. The two-phase refrigerant flow distributions among different channels of the BPHEs are first experimentally quantified based on the infrared radiation images of the heat exchanger sidewalls. To assess the performance degradation due to two-phase flow maldistribution, the thermal performance of the maldistributed BPHEs is measured and compared with that of the uniform-distributed BPHEs, predicted using an experimentally validated BPHE model. Several factors that affect the performance degradation are analyzed, including refrigerant inlet vapor quality and flow rate, exit superheat, and water-side flow distribution. The results show that, under the tested conditions, the two-phase refrigerant flow maldistribution causes up to a 29.8 % degradation in heat exchanger capacity and a 4 K decrease in evaporation temperature. Furthermore, the vulnerability of the BPHEs to two-phase flow maldistribution is significantly influenced by the operating conditions. This is closely related to the varying size ratios of superheated and two-phase regions within the BPHEs under different operating conditions. Reducing the refrigerant exit superheat is verified to help overcome the adverse effect brought by the flow maldistribution in the evaporators, but brings challenges in the system control. Matching the flow distributions of liquid refrigerant and water also proves beneficial to the BPHE evaporators.

Scaled experiment and numerical study on the effect of a novel makeup air system on smoke control in atrium fires

In the field of fire safety engineering, makeup air systems design for atriums is always a concern for researchers and engineers. A novel makeup air system integrating breathing zone and underfloor air supply (BUMS) was applied in the investigated atrium in this study. However, the flow rates proportional relationship of makeup air required for effective smoke control on each floor remains unclear. Scaled experiments and numerical simulations were carried out to address this problem in this study. The results show that the BUMS effectively controlled temperature, followed by visibility, with minimal impact on carbon monoxide (CO) concentrations. When the flow rate of makeup air was equal on each floor, the sizes of the safe evacuation passageways created by the BUMS on the third and fourth floors were not conducive to the safe evacuation of personnel. Based on the proportion of the CO concentration of each floor under natural ventilation, distribution proportions of makeup air flow rate were obtained. The effective action area was used to help determine the upper limit flow rate ratio for each floor. Calculation results revealed that when the flow rate distribution proportion for the first to fourth floors was 15 %, 27 %, 28 % and 30 %, respectively, the smoke control effectiveness has been improved. These findings are expected are expected to provide a theoretical basis for designing makeup air systems in atrium and large-space buildings.

Non‐isothermal flow dynamics during reverse roll coating phenomena of non‐Newtonian polymer

AbstractRoll coating plays a major role in industrial coating, including wallpapers, plastic and photographic films, sticky tapes, magnetic recordings, wrapping magazines and books, and so on. The current study proposes a mathematical model for a non‐isothermal, incompressible Sutterby fluid flowing through a narrow gap between two heated, counterrotating rollers. Lubrication approximation theory is used to simplify nondimensional expressions. The perturbation technique provides exact results for velocity profile, temperature, flow rate, and pressure gradient, whereas the numerical technique (Simpson rule) is used to compute the pressure profile and flow rate, respectively. The effects of the involved parameters on various physical characteristics like pressure, flow rate, temperature, pressure gradient, force, and power input are depicted in graphs and tabular form. A mechanism for controlling the coating thickness, power input, flow rate, separation force, and pressure distribution is provided by the material properties involved. With variations to a Sutterby fluid parameter and the velocity ratio K, both the pressure gradient and pressure decrease. It is significant to note that the temperature distribution is controlled by the velocities ratio and Brinkman number. Moreover, the separation point is shifted towards the nip area and the coating thickness on the web reduces with increasing velocities ratio K.

UPDATED MATHEMATICAL MODEL OF THE TRUNK PIPELINE SECTION WITH LOOPING IN CASE OF LEAKAGE

Pipeline transport is the main method of moving volumes of oil and oil products, as it has several advantages, for example, quality safety and high throughput. At the same time, it also has some disadvantages, which can include the difficulty of increasing throughput without high capital investments and possible leaks and unauthorized tie-ins into pipelines, which are sometimes difficult to identify. Looping is one of the main ways to increase the throughput of the existing pipeline system, when the need to increase the volume of oil and oil products supplies requires the construction of several parallel pipelines — looping. At the same time, the basic mathematical model for calculating the section of the trunk oil pipeline with looping does not fully cover all the physical processes associated with the distribution of flows in the parallel section, which arise due to turbulences in the connection nodes, as well as head loss along the length of the jumpers. Taking into account the resulting flow redistributions, it becomes difficult to detect a potential leak, since any of the detection methods can not solve the problem with high accuracy with determining the leak coordinate with an unknown distribution of fluid flow rates in a parallel section. As part of the study, an updated mathematical model of the oil export pipeline section with looping was developed, which took into account the presence of a potential leak in the parallel section. Mathematical modeling was carried out on the basis of the assumption of a known coordinate, the effect of leakage on the redistribution of flows in a parallel section was estimated. To assess the impact of the leak, various coordinates and volumetric flow rates were set, based on the calculation results, graphs of the flow rate in looping versus its length and internal diameter were plotted. The study is the initial phase of developing a hydraulic location leak detection system for the oil export pipeline section with looping. The need to refine the mathematical model for determining the flow rate in the parallel section is an urgent and necessary task, the solution of which will allow identifying potential leaks with reliable accuracy.

Design optimisation of a variable flow CO2 pipeline – A statistical approach

Pipeline transport has emerged as the most cost-effective method for transporting CO2 onshore. The CO2 flow rate is a key factor driving transport costs, underscoring the need to understand the impact of flow rate variability in CO2 pipeline networks on their design and economics. This paper presents an optimal pipeline design methodology for CO2 pipelines operating under variable flow that considers the probability distribution of CO2 flow rate and the pipeline length. The results imply that pipelines designed for optimal performance under variable flow rates often demand a higher level of overdesign compared to pipelines intended for steady-state conditions. Decision-makers must balance the trade-offs between pipeline oversizing and installing multiple pressure boosting stations, especially applicable to large transportation distances and projects of extended duration. The examination of different approaches to pipeline design reveals that a variable flow pipeline design based on mean flow rate is not recommended, because it is incapable of handling the maximum flow rate. This drawback is overcome by adopting the variable flow stochastic pipeline design presented in this paper.

The heat transfer enhancement mechanism of W-shaped micro-ribs on impingement cooling

In order to improve the performance of impingement cooling in turbine blades, W-shaped micro-ribs are arranged on the target surface to weaken the influence of cross-flow and generate longitudinal vortexes. The k-ω SST turbulence model is used for numerical simulation, and the mechanism of W-ribs on the enhancement of air impingement cooling performance and the influence of W-ribs geometry are studied. Under the condition of jet flow Reynolds number 4000 to 10000, based on the W-ribs with an angle α = 120°, width W = 0.1 mm and height H = 0.3 mm, the effects of angle (α = 90°, 120°, 150°, 180°), width (W = 0.05, 0.1, 0.2, 0.4 mm) and height (H = 1, 3, 4, 8 mm) on flow and heat transfer are investigated respectively. The results show that the W-ribs can separate the flow to both sides while generating longitudinal vortexes, which can reduce the influence of crossflow and enhance the sweep of the target surface to achieve the effect of enhancing heat transfer. α is the key factor affecting the generation of longitudinal vortex. W has little effect on the flow structure, but mainly affects the vorticity. H mainly affects the distribution of mass flow rate, which has a direct impact on the uniformity of upstream and downstream heat transfer. The W-ribs with α = 120°, W = 0.1 mm and H = 0.4 mm is recommended because of its excellent heat transfer performance and high comprehensive heat transfer performance, which are 39.7–48.9 % and 1.5–9.7 % higher than the smooth target surface, respectively.

Effect of number of channels on the fluctuations and distribution of humid-nitrogen condensing flow in a parallel flow field

Pressure drop and flow rate fluctuations of condensation-induced gas–liquid two-phase flow through a single microchannel have been widely recognized, which are quite different from that through parallel-connected multiple microchannels, as in the flow field of proton exchange membrane fuel cells and heat exchangers. This work investigates the effects of the number of channels on the fluctuations of the humid-nitrogen condensing flow and the uniformity of the flow rate distribution across the channels. A parallel flow field with eight channels is designed; the arrangement of these channels is axisymmetric, ensuring identical flow routes and thermal boundary conditions among individual channels. The number of open channels in the flow field could be adjusted from one to eight. The observed pressure drop fluctuated between 0.20 kPa and 0.68 kPa in the single-channel flow field, compared to 0.22 kPa to 0.43 kPa in the two-channel flow field and 0.28 kPa to 0.48 kPa in the eight-channel flow field. This indicates that the amplitude of pressure drop fluctuations decreases as the number of channels increases. In contrast, both the flow rate fluctuation and maldistribution are most significant in the eight-channel arrangement. It is also found that the effect of channel number on condensing flow dynamics is noticeable when both the velocity and humidity of the inlet gas are low.

Evaluation of plate heat exchangers comprising sections with different chevron angle arrangements

Gasketed Plate Heat Exchangers (GPHEs) are essential equipment for thermal control of oil facilities owing to their compactness and well-organized assembly concept. In order to satisfy heat loads of several megawatts, numerous and large plates are necessary. The application of GPHEs with plentiful plates yields flow maldistribution which affects GPHE thermal–hydraulic performance. With the purpose of matching different heat loads and processing conditions or as an alternative to reduce maldistribution effects, GPHEs with mixed arrangements, comprising sections where channels contain different geometrical features, may be selected. This works addresses the evaluation of flow and pressure distributions in GPHEs with mixed arrangements. The pressure drop in each GPHE section is initially approximated by the pressure drop in parallel pipes, neglecting the pressure drop in the ports. With this assumption, the inlet flow rate is divided into two different parts. The pressure distributions in each section are assumed to maintain flow features of m2-model for U-type arrangement. The friction factor for each GPHE segment was determined with independent experiments including only one section type and with a small number of channels as to avoid flow maldistribution. Friction factor correlations were determined for single section GPHEs where each channel can be formed with plates containing different chevron angles: 30°/30°, 60°/60°, 66°/66°, 30°/60° and 30°/66°, with channel Reynolds number ranging from 260 to 3,080. Subsequently, experiments with GPHEs containing mixed sections with high and low pressure drop channels were performed: 60°/60° and 30°/30°, 30°/30° and 66°/66°; with total number of channels up to 180. The channel pressure drop was measured in regularly spaced locations. The pressure drop in each GPHE segment has been predicted with the fractional RMS deviation equal to 2.2%. Although the mass flow rate per channel has not been measured, it is expected that the predicted flow rate distribution can reasonably stand for the real one and, therefore, assisting in proper thermal performance evaluation.

Analysis of geothermal heat recovery from abandoned coal mine water for clean heating and cooling: A case from Shandong, China

Exploiting heat from abandoned mine water using heat pumps is attracting increasing attention worldwide. China has a large number of abandoned coal mines, but there are few applications of geothermal heat recovery from mine water. Here, we report a new case from Shandong, North China. Based on actual geological and underground space structures, numerical analysis is performed to study the flow and heat transfer in the complex connection network among roadways. For given temperature differences on the source side, the underground heat transfer is greatly influenced by flow rates in roadways. Under the acceptable performances and long-term stability of heat pumps, the maximum heat transfer capacity of roadways reaches 697.3 kW at the flow rate of 100 m3/h, with the power per unit length of 151.3 W/m, indicating a promising heat recovery potential. When switching between the extraction and injection well, both fluid directions and the distribution of flow rates in roadways change obviously, thereby affecting the underground heat transfer. Besides dominant flow paths, there exist weak flow paths, which can be identified by numerical analysis, and necessary blocking can be applied to optimize the flow and heat transfer. Excessively long horizontal straight roadways for heat exchange should be avoided.

Optimization of the first wall cooling circuits to achieve an efficient coolant flow distribution in the DCLL breeding blanket

In the current design proposal of the Dual Coolant Lithium Lead (DCLL) breeding blanket for the EU DEMO reactor, the first wall is cooled down by a succession of He channels arranged along the poloidal length. In this scheme, since He extracts around 30% of the blanket thermal power, a conservative oversizing of the coolant mass flow rate compelled by highly heterogeneous heat loads, can severely jeopardize the blanket thermal efficiency and the consumption of the auxiliary systems. A meticulous design of the hydraulic network can enhance the cooling system efficiency by ensuring that each channel is fed with a mass flow rate according to its specific needs.In this study, a novel approach employing genetic algorithms has been developed to optimize the geometry of the first wall, encompassing the feeding/collecting manifolds and channels. This optimization aims to replicate the mass flow rate distribution obtained from TOMFLOW (Tool for Optimization of Helium Mass Flow Rate Distribution), a code dedicated to determine the minimum mass flow rate required in each channel to maintain the EUROFER steel structure below 550 °C.

Multi-Objective Topology Optimization of Conjugate Heat Transfer Using Level Sets and Anisotropic Mesh Adaptation

This study proposes a new computational framework for the multi-objective topology optimization of conjugate heat transfer systems using a continuous adjoint approach. It relies on a monolithic solver for the coupled steady-state Navier–Stokes and heat equations, which combines finite elements stabilized by the variational multi-scale method, level set representations of the fluid–solid interfaces and immersed modeling of heterogeneous materials (fluid–solid) to ensure that the proper amount of heat is exchanged to the ambient fluid by solid objects in arbitrary geometry. At each optimization iteration, anisotropic mesh adaptation is applied in near-wall regions automatically captured by the level set. This considerably cuts the computational effort associated with calling the finite element solver, in comparison to traditional topology optimization algorithms operating on isotropic grids with a comparable refinement level. Given that we operate within the constraint of a specified number of nodes in the mesh, this allows not only to improve the accuracy of interface representation and motion but also to retain the high fidelity of the numerical solutions at the grid points just adjacent to the interface. Finally, the remeshing and resolution steps both run within a highly parallel environment, which makes it possible for the proposed algorithm to tackle large-scale problems in three dimensions with several tens of millions of state degrees of freedom. The developed solver is validated first by minimizing dissipation in a flow splitter device, for which the method delivers relevant optimal designs over a wide range of volume constraints and flow rate distributions over the multiple outlet orifices but yields better accuracy compared to reference data from literature obtained using uniform meshes (in the sense that the layouts are more smooth, and the solutions are better resolved). The scheme is then applied to a two-dimensional heat transfer problem, using bi-objective cost functionals combining flow resistance and thermal recoverable power. A comprehensive parametric study reveals a complex arrangement of optimal solutions on the Pareto front, with multiple branches of symmetric and asymmetric designs, some of them previously unreported. Finally, the algorithmic developments are substantiated with several three-dimensional numerical examples tackled under fixed weights for heat transfer and flow resistance, for which we show that the optimal layouts computed at low Reynolds number, that are intrinsically relevant to a broad range of microfluidic application, can also serve as smooth solutions to high-Reynolds-number engineering problems of practical interest.

Study on the Influence of Perforating Parameters on the Flow Rate and Stress Distribution of Multi-Fracture Competitive Propagation

It is of great significance to investigate the flow rate and stress distribution of multi-fracture propagation for the optimization of perforation parameters and fracture parameters. Considering the coupling of rock deformation, fracture direction and fluid flow in multi-fracture scenarios, a mathematical model and solution program for the flow and stress distribution of multiple fractures are established, and the analytical model is used for comparison and verification. The effects of perforation cluster number, cluster spacing, perforation diameter on fracture extension trajectory, fracture width, flow rate of each fracture and stress field are studied by the model. The results show that, as the number of perforating clusters increases, the inner fracture is inhibited more severely with less width, length and flow distribution, as well as lower bottom hole pressure. With the increase in cluster spacing, the stress interference between whole fractures is weakened and the flow distribution of the inner fracture is increased with lower bottom hole pressure. With the decrease in perforation diameter, the inhibition effect of inside fractures is weakened, while the inhibition effect of outside fractures, the flow distribution of inside fractures and the bottom hole pressure are increased. The uniform propagation of multiple fractures can be promoted by decreasing the perforation clusters’ number and perforation diameter or increasing fracture spacing.