An inverse analysis of fluid flow through granular media using differentiable lattice Boltzmann method

This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.

Computational fluid dynamics modeling of rapid pyrolysis of solid waste magnesium nitrate hydrate under different injection methods

This study developed a numerical model to efficiently treat solid waste magnesium nitrate hydrate through multi-step chemical reactions.The model simulates two-phase flow,heat,and mass transfer processes in a pyrolysis furnace to improve the decomposition rate of magnesium nitrate.The performance of multi-nozzle and single-nozzle injection methods was evaluated,and the effects of primary and secondary nozzle flow ratios,velocity ratios,and secondary nozzle inclination angles on the decomposition rate were investigated.Results indicate that multi-nozzle injection has a higher conversion efficiency and decomposition rate than single-nozzle injection,with a 10.3%higher conversion rate under the design parameters.The decomposition rate is primarily dependent on the average residence time of particles,which can be increased by decreasing flow rate and velocity ratios and increasing the inclination angle of secondary nozzles.The optimal parameters are injection flow ratio of 40%,injection velocity ratio of 0.6,and secondary nozzle inclination of 30°,corresponding to a maximum decomposition rate of 99.33%.

Optimization method of fracturing fluid volume intensity for SRV fracturing technique in shale oil reservoir based on forced imbibition:A case study of well X-1 in Biyang Sag of Nanxiang Basin,China

An optimization method of fracturing fluid volume strength was introduced taking well X-1 in Biyang Sag of Nanxiang Basin as an example.The characteristic curves of capillary pressure and relative permeability were obtained from history matching between forced imbibition experimental data and core-scale reservoir simulation results and taken into a large scale reservoir model to mimic the forced imbibition behavior during the well shut-in period after fracturing.The optimization of the stimulated reservoir volume(SRV)fracturing fluid volume strength should meet the requirements of estimated ultimate recovery(EUR),increased oil recovery by forced imbibition and enhancement of formation pressure and the fluid volume strength of fracturing fluid should be controlled around a critical value to avoid either insufficiency of imbibition displacement caused by insufficient fluid amount or increase of costs and potential formation damage caused by excessive fluid amount.Reservoir simulation results showed that SRV fracturing fluid volume strength positively correlated with single-well EUR and an optimal fluid volume strength existed,above which the single-well EUR increase rate kept decreasing.An optimized increase of SRV fracturing fluid volume and shut-in time would effectively increase the formation pressure and enhance well production.Field test results of well X-1 proved the practicality of established optimization method of SRV fracturing fluid volume strength on significant enhancement of shale oil well production.

Experimental study on secondary air mixing along the bed height in a circulating fluidized bed with a multitracer-gas method

A multitracer-gas method was proposed to study the secondary air(SA)mixing along the bed height in a circulating fluidized bed(CFB)using carbon monoxide(CO),oxygen(O_(2)),and carbon dioxide(CO_(2))as tracer gases.Experiments were carried out on a cold CFB test rig with a cross-section of 0.42 m×0.73 m and a height of 5.50 m.The effects of superficial velocity,SA ratio,bed inventory,and particle diameter on the SA mixing were investigated.The results indicate that there are some differences in the measurement results obtained using different tracer gases,wherein the deviation between CO and CO_(2) ranges from 42%to 66%and that between O_(2) and CO_(2) ranges from 45%to 71%in the lower part of the fluidized bed.However,these differences became less pronounced as the bed height increased.Besides,the high solid concentration and fine particle diameter in the CFB may weaken the difference.The measurement results of different tracer gases show the same trends under the variation of operating parameters.Increasing superficial velocity and SA ratio and decreasing particle diameter result in better mixing of the SA.The effect of bed inventory on SA mixing is not monotonic.

A two-phase type-curve method with multiscale fluid transport mechanisms in hydraulically fractured shale reservoirs

The quantitative understanding of hydraulic fracture(HF)properties guides accurate production forecasts and reserve estimation.Type curve is a powerful technique to characterize HF and reservoir properties from flowback and long-term production data.However,two-phase flow of water and hydrocarbon after an HF stimulation together with the complex transport mechanisms in shale nanopores exacerbate the nonlinearity of the transport equation,causing errors in type-curve analysis.Accordingly,we propose a new two-phase type-curve method to estimate HF properties,such as HF volume and permeability of fracture,through the analysis of flowback data of multi-fractured shale wells.The proposed type curve is based on a semianalytical solution that couples the two-phase flow from the matrix with the flow in HF by incorporating matrix influx,slippage effect,stress dependence,and the spatial variation of fluid properties in inorganic and organic pores.For the first time,multiple fluid transport mechanisms are considered into two-phase type-curve analysis for shale reservoirs.We analyze the flowback data from a multi-fractured horizontal well in a shale gas reservoir to verify the field application of the proposed method.The results show that the fracture properties calculated by the type-curve method are in good agreement with the long-time production data.

Effect of fracture fluid flowback on shale microfractures using CT scanning

The field data of shale fracturing demonstrate that the flowback performance of fracturing fluid is different from that of conventional reservoirs,where the flowback rate of shale fracturing fluid is lower than that of conventional reservoirs.At the early stage of flowback,there is no single-phase flow of the liquid phase in shale,but rather a gas-water two-phase flow,such that the single-phase flow model for tight oil and gas reservoirs is not applicable.In this study,pores and microfractures are extracted based on the experimental results of computed tomography(CT)scanning,and a spatial model of microfractures is established.Then,the influence of rough microfracture surfaces on the flow is corrected using the modified cubic law,which was modified by introducing the average deviation of the microfracture height as a roughness factor to consider the influence of microfracture surface roughness.The flow in the fracture network is simulated using the modified cubic law and the lattice Boltzmann method(LBM).The results obtained demonstrate that most of the fracturing fluid is retained in the shale microfractures,which explains the low fracturing fluid flowback rate in shale hydraulic fracturing.

Numerical Assessment of Nanofluid Natural Convection Using Local RBF Method Coupled with an Artificial Compressibility Model

In this paper,natural heat convection inside square and equilateral triangular cavities was studied using a meshless method based on collocation local radial basis function(RBF).The nanofluids used were Cu-water or Al_(2)O_(3)-water mixture with nanoparticle volume fractions range of 0≤φ≤0.2.A system of continuity,momentum,and energy partial differential equations was used in modeling the flow and temperature behavior of the fluids.Partial derivatives in the governing equations were approximated using the RBF method.The artificial compressibility model was implemented to overcome the pressure velocity coupling problem that occurs in such equations.Themain goal of this work was to present a simple and efficient method to deal with complex geometries for a variety of problem conditions.To assess the accuracy of the proposed method,several test cases of natural convection in square and triangular cavities were selected.For Rayleigh numbers ranging from 103 to 105,a validation test of natural convection of Cu-water in a square cavity was used.The numerical investigation was then extended to Rayleigh number 106,as well as Al_(2)O_(3)-water nanofluid with a volume fraction range of 0≤φ≤0.2.In a second investigation,the same nanofluids were used in a triangular cavitywith varying volume fractions to test the proposed meshless approach on non-rectangular geometries.The numerical results appear to be in agreement with those from earlier investigations.Furthermore,the suggested meshless method was found to be stable and accurate,demonstrating that it may be a viable alternative for solving natural heat transfer equations of nanofluids in enclosures with irregular geometries.

Computational fluid dynamics-discrete element method simulation of stirred tank reactor for graphene production

Liquid phase exfoliation(LPE)process for graphene production is usually carried out in stirred tank reactor and the interactions between the solvent and the graphite particles are important as to improve the production efficiency.In this paper,these interactions were revealed by computational fluid dynamics–discrete element method(CFD-DEM)method.Based on simulation results,both liquid phase flow hydrodynamics and particle motion behavior have been analyzed,which gave the general information of the multiphase flow behavior inside the stirred tank reactor as to graphene production.By calculating the threshold at the beginning of graphite exfoliation process,the shear force from the slip velocity was determined as the active force.These results can support the optimization of the graphene production process.

A Low Mach Number IMEX Flux Splitting for the Level Set Ghost Fluid Method

Considering droplet phenomena at low Mach numbers,large differences in the magnitude of the occurring characteristic waves are presented.As acoustic phenomena often play a minor role in such applications,classical explicit schemes which resolve these waves suffer from a very restrictive timestep restriction.In this work,a novel scheme based on a specific level set ghost fluid method and an implicit-explicit(IMEX)flux splitting is proposed to overcome this timestep restriction.A fully implicit narrow band around the sharp phase interface is combined with a splitting of the convective and acoustic phenomena away from the interface.In this part of the domain,the IMEX Runge-Kutta time discretization and the high order discontinuous Galerkin spectral element method are applied to achieve high accuracies in the bulk phases.It is shown that for low Mach numbers a significant gain in computational time can be achieved compared to a fully explicit method.Applica-tions to typical droplet dynamic phenomena validate the proposed method and illustrate its capabilities.

A Review on Sources,Extractions and Analysis Methods of a Sustainable Biomaterial:Tannins

Condensed and hydrolysable tannins are non-toxic natural polyphenols that are a commercial commodity industrialized for tanning hides to obtain leather and for a growing number of other industrial applications mainly to substitute petroleum-based products.They are a definite class of sustainable materials of the forestry industry.They have been in operation for hundreds of years to manufacture leather and now for a growing number of applications in a variety of other industries,such as wood adhesives,metal coating,pharmaceutical/medical applications and several others.This review presents the main sources,either already or potentially commercial of this forestry by-materials,their industrial and laboratory extraction systems,their systems of analysis with their advantages and drawbacks,be these methods so simple to even appear primitive but nonetheless of proven effectiveness,or very modern and instrumental.It constitutes a basic but essential summary of what is necessary to know of these sustainable materials.In doing so,the review highlights some of the main challenges that remain to be addressed to deliver the quality and economics of tannin supply necessary to fulfill the industrial production requirements for some materials-based uses.

Drilling-based measuring method for the c-φ parameter of rock and its field application

The technology of drilling tests makes it possible to obtain the strength parameter of rock accurately in situ. In this paper, a new rock cutting analysis model that considers the influence of the rock crushing zone(RCZ) is built. The formula for an ultimate cutting force is established based on the limit equilibrium principle. The relationship between digital drilling parameters(DDP) and the c-φ parameter(DDP-cφ formula, where c refers to the cohesion and φ refers to the internal friction angle) is derived, and the response of drilling parameters and cutting ratio to the strength parameters is analyzed. The drillingbased measuring method for the c-φ parameter of rock is constructed. The laboratory verification test is then completed, and the difference in results between the drilling test and the compression test is less than 6%. On this basis, in-situ rock drilling tests in a traffic tunnel and a coal mine roadway are carried out, and the strength parameters of the surrounding rock are effectively tested. The average difference ratio of the results is less than 11%, which verifies the effectiveness of the proposed method for obtaining the strength parameters based on digital drilling. This study provides methodological support for field testing of rock strength parameters.

Reservoir quality evaluation of the Narimba Formation in Bass Basin,Australia:Implications from petrophysical analysis,sedimentological features,capillary pressure and wetting fluid saturation

The evaluation of reservoir quality was accomplished on the Late Paleocene to Early Eocene Narimba Formation in Bass Basin,Australia.This study involved combination methods such as petrophysical analysis,petrography and sedimentological studies,reservoir quality and fluid flow units from derivative parameters,and capillary pressure and wetting fluid saturation relationship.Textural and diagenetic features are affecting the reservoir quality.Cementation,compaction,and presence of clay minerals such as kaolinite are found to reduce the quality while dissolution and secondary porosity are noticed to improve it.It is believed that the Narimba Formation is a potential reservoir with a wide range of porosity and permeability.Porosity ranges from 3.1%to 25.4%with a mean of 15.84%,while permeability ranges between 0.01 mD and 510 mD,with a mean of 31.05 mD.Based on the heterogenous lithology,the formation has been categorized into five groups based on permeability variations.Group I showed an excellent to good quality reservoir with coarse grains.The impacts of both textural and diagenetic features improve the reservoir and producing higher reservoir quality index(RQI)and flow zone indicators(FZI)as well as mostly mega pores.The non-wetting fluid migration has the higher possibility to flow in the formation while displacement pressure recorded as zero.Group II showed a fair quality reservoir with lower petrophysical properties in macro pores.The irreducible water saturation is increasing while the textural and digenetic properties are still enhancing the reservoir quality.Group III reflects lower quality reservoir with mostly macro pores and higher displacement pressure.It may indicate smaller grain size and increasing amount of cement and clay minerals.Group IV,and V are interpreted as a poor-quality reservoir that has lower RQI and FZI.The textural and digenetic features are negatively affecting the reservoir and are leading to smaller pore size and pore throat radii(r35)values to be within the range of macro,meso-,micro-,and nano pores.The capillary displacement pressure curves of the three groups show increases reaching the maximum value of 400 psia in group V.Agreement with the classification of permeability,r35 values,and pore type can be used in identifying the quality of reservoir.

An MPI parallel DEM-IMB-LBM framework for simulating fluid-solid interaction problems

The high-resolution DEM-IMB-LBM model can accurately describe pore-scale fluid-solid interactions,but its potential for use in geotechnical engineering analysis has not been fully unleashed due to its prohibitive computational costs.To overcome this limitation,a message passing interface(MPI)parallel DEM-IMB-LBM framework is proposed aimed at enhancing computation efficiency.This framework utilises a static domain decomposition scheme,with the entire computation domain being decomposed into multiple subdomains according to predefined processors.A detailed parallel strategy is employed for both contact detection and hydrodynamic force calculation.In particular,a particle ID re-numbering scheme is proposed to handle particle transitions across sub-domain interfaces.Two benchmarks are conducted to validate the accuracy and overall performance of the proposed framework.Subsequently,the framework is applied to simulate scenarios involving multi-particle sedimentation and submarine landslides.The numerical examples effectively demonstrate the robustness and applicability of the MPI parallel DEM-IMB-LBM framework.

A comparison study on structure-function relationship of polysaccharides obtained from sea buckthorn berries using different methods:antioxidant and bile acid-binding capacity

In this study,the structural characters,antioxidant activities and bile acid-binding ability of sea buckthorn polysaccharides(HRPs)obtained by the commonly used hot water(HRP-W),pressurized hot water(HRP-H),ultrasonic(HRP-U),acid(HRP-C)and alkali(HRP-A)assisted extraction methods were investigated.The results demonstrated that extraction methods had significant effects on extraction yield,monosaccharide composition,molecular weight,particle size,triple-helical structure,and surface morphology of HRPs except for the major linkage bands.Thermogravimetric analysis showed that HRP-U with filamentous reticular microstructure exhibited better thermal stability.The HRP-A with the lowest molecular weight and highest arabinose content possessed the best antioxidant activities.Moreover,the rheological analysis indicated that HRPs with higher galacturonic acid content and molecular weight showed higher viscosity and stronger crosslinking network(HRP-C,HRP-W and HRP-U),which exhibited stronger bile acid binding capacity.The present findings provide scientific evidence in the preparation technology of sea buckthorn polysaccharides with good antioxidant and bile acid binding capacity which are related to the structure affected by the extraction methods.

Synthetic polymers:A review of applications in drilling fluids

With the growth of deep drilling and the complexity of the well profile,the requirements for a more complete and efficient exploitation of productive formations increase,which increases the risk of various complications.Currently,reagents based on modified natural polymers(which are naturally occurring compounds)and synthetic polymers(SPs)which are polymeric compounds created industrially,are widely used to prevent emerging complications in the drilling process.However,compared to modified natural polymers,SPs form a family of high-molecular-weight compounds that are fully synthesized by undergoing chemical polymerization reactions.SPs provide substantial flexibility in their design.Moreover,their size and chemical composition can be adjusted to provide properties for nearly all the functional objectives of drilling fluids.They can be classified based on chemical ingredients,type of reaction,and their responses to heating.However,some of SPs,due to their structural characteristics,have a high cost,a poor temperature and salt resistance in drilling fluids,and degradation begins when the temperature reaches 130℃.These drawbacks prevent SP use in some medium and deep wells.Thus,this review addresses the historical development,the characteristics,manufacturing methods,classification,and the applications of SPs in drilling fluids.The contributions of SPs as additives to drilling fluids to enhance rheology,filtrate generation,carrying of cuttings,fluid lubricity,and clay/shale stability are explained in detail.The mechanisms,impacts,and advances achieved when SPs are added to drilling fluids are also described.The typical challenges encountered by SPs when deployed in drilling fluids and their advantages and drawbacks are also discussed.Economic issues also impact the applications of SPs in drilling fluids.Consequently,the cost of the most relevant SPs,and the monomers used in their synthesis,are assessed.Environmental impacts of SPs when deployed in drilling fluids,and their manufacturing processes are identified,together with advances in SP-treatment methods aimed at reducing those impacts.Recommendations for required future research addressing SP property and performance gaps are provided.

Effect of bubble morphology and behavior on power consumption in non-Newtonian fluids’aeration process

Due to a prolonged operation time and low mass transfer efficiency, the primary challenge in the aeration process of non-Newtonian fluids is the high energy consumption, which is closely related to the form and rate of impeller, ventilation, rheological properties and bubble morphology in the reactor. In this perspective, through optimal computational fluid dynamics models and experiments, the relationship between power consumption, volumetric mass transfer rate(kLa) and initial bubble size(d0) was constructed to establish an efficient operation mode for the aeration process of non-Newtonian fluids. It was found that reducing the d0could significantly increase the oxygen mass transfer rate, resulting in an obvious decrease in the ventilation volume and impeller speed. When d0was regulated within 2-5 mm,an optimal kLa could be achieved, and 21% of power consumption could be saved, compared to the case of bubbles with a diameter of 10 mm.

A particle-resolved heat-particle-fluid coupling model by DEM-IMB-LBM

Multifield coupling is frequently encountered and also an active area of research in geotechnical engineering.In this work,a particle-resolved direct numerical simulation(PR-DNS)technique is extended to simulate particle-fluid interaction problems involving heat transfer at the grain level.In this extended technique,an immersed moving boundary(IMB)scheme is used to couple the discrete element method(DEM)and lattice Boltzmann method(LBM),while a recently proposed Dirichlet-type thermal boundary condition is also adapted to account for heat transfer between fluid phase and solid particles.The resulting DEM-IBM-LBM model is robust to simulate moving curved boundaries with constant temperature in thermal flows.To facilitate the understanding and implementation of this coupled model for non-isothermal problems,a complete list is given for the conversion of relevant physical variables to lattice units.Then,benchmark tests,including a single-particle sedimentation and a two-particle drafting-kissing-tumbling(DKT)simulation with heat transfer,are carried out to validate the accuracy of our coupled technique.To further investigate the role of heat transfer in particle-laden flows,two multiple-particle problems with heat transfer are performed.Numerical examples demonstrate that the proposed coupling model is a promising high-resolution approach for simulating the heat-particle-fluid coupling at the grain level.

High-temperature corrosion of phosphate completion fluid and corrosion inhibition method by membrane transformation

By analyzing the corrosion of phosphate completion fluid on the P110 steel at 170 °C, the high-temperature corrosion mechanism of phosphate completion fluid was revealed, and a corrosion inhibition method by membrane transformation was proposed and an efficient membrane-forming agent was selected. Scanning electron microscope (SEM) images, X-ray energy spectrum and X-ray diffraction results were used to characterize the microscopic morphology, elemental composition and phase composition of the precipitation membrane on the surface of the test piece. The effect and mechanism of corrosion inhibition by membrane transformation were clarified. The phosphate completion fluid eroded the test piece by high-temperature water vapor and its hydrolyzed products to form a membrane of iron phosphate corrosion product. By changing the corrosion reaction path, the Zn2+ membrane-forming agent could generate KZnPO4 precipitation membrane with high temperature resistance, uniform thickness and tight crystal packing on the surface of the test piece, which could inhibit the corrosion of the test piece, with efficiency up to 69.63%. The Cu2+ membrane-forming agent electrochemically reacted with Fe to precipitate trace elemental Cu on the surface of the test piece, thus forming a protective membrane, which could inhibit metal corrosion, with efficiency up to 96.64%, but the wear resistance was poor. After combining 0.05% Cu2+ and 0.25% Zn2+, a composite protective membrane of KZnPO4 crystal and elemental Cu was formed on the surface of the test piece. The corrosion inhibition efficiency reached 93.03%, which ensured the high corrosion inhibition efficiency and generated a precipitation membrane resistant to temperature and wear.

Application of Fluid Phase Recovery Method in the Early Development of BZ26-6 Volatile Oil Field in Bohai Sea

BZ26-6 Oilfield is a kind of deep metamorphic rock buried-hill volatile oilfield in Bohai Sea, China. Its early development plan is restricted due to the simultaneous production of oil and gas in large sections of reservoirs, unclear understanding of formation fluid properties and uncertainty of gas-oil interface. Through theoretical research on phase recovery and experimental analysis of crude oil phase characteristics in the original formation, characteristic parameters of the equilibrium condensate gas fluid are restored and calculated. Through the superimposed phase diagram of volatile oil and condensate gas, BZ26-6 Oilfield is determined to be a volatile oil reservoir with a condensate gas cap, with formation pressure and saturation pressure of 36.1 MPa, respectively. Based on the research results of oil-gas phase behavior characteristics, the thermodynamic equations and equation of state are jointly used to solve the problem, and the content change curves of each component at different depths are drawn. Combined with the sensitivity analysis of numerical simulation, the gas-oil interface is determined to be -3726 m above sea level. The fluid phase analysis software, Fluidmodeler, is used to simulate volatile oil degassing and condensate gas separation experiments. In combination with oil and gas production data obtained through the production test, the specific oil recovery index and the specific gas recovery index are determined to be 0.408 m3/(MPa·d·m) and 1195 m3/(MPa·d·m), respectively. And the reasonable production capacity prediction is conducted on the early development of BZ26-6 Oilfield. The research results can provide a theoretical basis for the efficient development of similar complex oil and gas reservoirs.

Endoscopic intramural cystogastrostomy for treatment of peripancreatic fluid collection: A viewpoint from a surgeon

Percutaneous or endoscopic drainage is the initial choice for the treatment of peripancreatic fluid collection in symptomatic patients.Endoscopic transgastric fenestration(ETGF)was first reported for the management of pancreatic pseu-docysts of 20 patients in 2008.From a surgeon’s viewpoint,ETGF is a similar procedure to cystogastrostomy in that they both produce a wide outlet orifice for the drainage of fluid and necrotic debris.ETGF can be performed at least 4 wk after the initial onset of acute pancreatitis and it has a high priority over the surgical approach.However,the surgical approach usually has a better success rate because surgical cystogastrostomy has a wider outlet(>6 cm vs 2 cm)than ETGF.However,percutaneous or endoscopic drainage,ETGF,and surgical approach offer various treatment options for peripancreatic fluid collection patients based on their conditions.