Recent progresses in the development of tubular segmented-in-series solid oxide fuel cells:Experimental and numerical study

Solid oxide fuel cells(SOFCs)have attracted a great deal of interest because they have the highest efficiency without using any noble metal as catalysts among all the fuel cell technologies.However,traditional SOFCs suffer from having a higher volume,current leakage,complex connections,and difficulty in gas sealing.To solve these problems,Rolls-Royce has fabricated a simple design by stacking cells in series on an insulating porous support,resulting in the tubular segmented-in-series solid oxide fuel cells(SIS-SOFCs),which achieved higher output voltage.This work systematically reviews recent advances in the structures,preparation methods,perform-ances,and stability of tubular SIS-SOFCs in experimental and numerical studies.Finally,the challenges and future development of tubular SIS-SOFCs are also discussed.The findings of this work can help guide the direction and inspire innovation of future development in this field.

Numerical study on the cavity dynamics for vertical water entries of twin spheres

In this study, a three dimensional(3D) numerical model of six-degrees-of-freedom(6DOF) is applied to simulate the water entries of twin spheres side-by-side at different lateral distances and time intervals.The turbulence structure is described using the shear-stress transport k-ω(SST k-ω) model, and the volume of fluid(VOF) method is used to track the complex air-liquid interface. The motion of spheres during water entry is simulated using an independent overset grid. The numerical model is verified by comparing the cavity evolution results from simulations and experiments. Numerical results reveal that the time interval between the twin water entries evidently affects cavity expansion and contraction behaviors in the radial direction. However, this influence is significantly weakened by increasing the lateral distance between the two spheres. In synchronous water entries, pressure is reduced on the midline of two cavities during surface closure, which is directly related to the cavity volume. The evolution of vortexes inside the two cavities is analyzed using a velocity vector field, which is affected by the lateral distance and time interval of water entries.

Numerical Study on the Effect of Vortex Generators on the Aerodynamic Drag of a High-Speed Train

A relatively high aerodynamic drag is an important factor that hinders the further acceleration of high-speed trains.Using the shear stress transport(SST)k-ωturbulence model,the effect of various vortex generator types on the aerodynamic characteristics of an ICE2(Inter-city Electricity)train has been investigated.The results indi-cate that the vortex generators with wider triangle,trapezoid,and micro-ramp arranged on the surface of the tail car can significantly change the distribution of surface pressure and affect the vorticity intensity in the wake.This alteration effectively reduces the resistance of the tail car.Meanwhile,the micro-ramp vortex generator with its convergent structure at the rear exhibits enhancedflow-guiding capabilities,resulting in a 15.4%reduction in the drag of the tail car.

Numerical and theoretical study of load transfer behavior during cascading pillar failure

To further study the load transfer mechanism of roofemulti-pillarefloor system during cascading pillar failure(CPF),numerical simulation and theoretical analysis were carried out to study the three CPF modes according to the previous experimental study on treble-pillar specimens,e.g.successive failure mode(SFM),domino failure mode(DFM)and compound failure mode(CFM).Based on the finite element code rock failure process analysis(RFPA^(2D)),numerical models of treble-pillar specimen with different mechanical properties were established to reproduce and verify the experimental results of the three CPF modes.Numerical results show that the elastic rebound of roofefloor system induced by pillar instability causes dynamic disturbance to adjacent pillars,resulting in sudden load increases and sudden jump displacement of adjacent pillars.The phenomena of load transfer in the roofemulti-pillarefloor system,as well as the induced accelerated damage behavior in adjacent pillars,were discovered and studied.In addition,based on the catastrophe theory and the proposed mechanical model of treble-pillar specimen edisc spring group system,a potential function that characterizes the evolution characteristics of roof emulti-pillarefloor system was established.The analytical expressions of sudden jump and energy release of treble-pillar specimenedisc spring group system of the three CPF modes were derived according to the potential function.The numerical and theoretical results show good agreement with the experimental results.This study further reveals the physical essence of load transfer during CPF of roof emulti-pillarefloor system,which provides references for mine design,construction and disaster prevention.

A Numerical Study on the Effect of the Backflow Hole Position on the Performances of a Self-Priming Pump

A self-priming pump is a centrifugal pump that has the ability to prime itself. Typically, its performance dependson the configuration of its reflux hole. In this study, the ANSYS FLUENT software is used to investigate the effectsof three different radial positions of the reflux hole on gas-liquid two-phase distribution, pressure pulsation, andimp during self-priming. The research results indicate that: (1) The effective channel size for the reflux liquid toenter the volute varies depending on the location of the reflux hole. The effect of the impeller rotation on thereflux liquid becomes more obvious as the setting distance of the reflux aperture decreases. (2) The position ofthe reflux hole significantly affects the gas phase mass fraction inside the impeller, resulting in a significant reductionin the time it takes for the mass fraction to exceed 80%. (3) The position of the reflux hole significantly affectsthe average pressure on each monitoring surface. (4) Placing the reflux hole at a excessively distant radial distancecan result in an excessive vertical component. (5) The self-priming performance of the pump can be improved tosome extent by placing the return hole at a small radial distance.

Operation optimization of prefabricated light modular radiant heating system:Thermal resistance analysis and numerical study

The utilization of prefabricated light modular radiant heating system has demonstrated significant increases in heat transfer efficiency and energy conservation capabilities.Within prefabricated building construction,this new heating method presents an opportunity for the development of comprehensive facilities.The parameters for evaluating the effectiveness of such a system are the upper surface layer’s heat flux and temperature.In this paper,thermal resistance analysis calculation based on a simplified model for this unique radiant heating system analysis is presented with the heat transfer mechanism’s evaluation.The results obtained from thermal resistance analysis calculation and numerical simulation indicate that the thermal resistance analysis method is highly accurate with temperature discrepancies ranging from 0.44℃ to−0.44℃ and a heat flux discrepancy of less than 7.54%,which can meet the requirements of practical engineering applications,suggesting a foundation for the prefabricated radiant heating system.

A Hybrid SIR-Fuzzy Model for Epidemic Dynamics:A Numerical Study

This study focuses on the urgent requirement for improved accuracy in diseasemodeling by introducing a newcomputational framework called the Hybrid SIR-Fuzzy Model.By integrating the traditional Susceptible-Infectious-Recovered(SIR)modelwith fuzzy logic,ourmethod effectively addresses the complex nature of epidemic dynamics by accurately accounting for uncertainties and imprecisions in both data and model parameters.The main aim of this research is to provide a model for disease transmission using fuzzy theory,which can successfully address uncertainty in mathematical modeling.Our main emphasis is on the imprecise transmission rate parameter,utilizing a three-part description of its membership level.This enhances the representation of disease processes with greater complexity and tackles the difficulties related to quantifying uncertainty in mathematical models.We investigate equilibrium points for three separate scenarios and perform a comprehensive sensitivity analysis,providing insight into the complex correlation betweenmodel parameters and epidemic results.In order to facilitate a quantitative analysis of the fuzzy model,we propose the implementation of a resilient numerical scheme.The convergence study of the scheme demonstrates its trustworthiness,providing a conditionally positive solution,which represents a significant improvement compared to current forward Euler schemes.The numerical findings demonstrate themodel’s effectiveness in accurately representing the dynamics of disease transmission.Significantly,when the mortality coefficient rises,both the susceptible and infected populations decrease,highlighting the model’s sensitivity to important epidemiological factors.Moreover,there is a direct relationship between higher Holling type rate values and a decrease in the number of individuals who are infected,as well as an increase in the number of susceptible individuals.This correlation offers a significant understanding of how many elements affect the consequences of an epidemic.Our objective is to enhance decision-making in public health by providing a thorough quantitative analysis of the Hybrid SIR-Fuzzy Model.Our approach not only tackles the existing constraints in disease modeling,but also paves the way for additional investigation,providing a vital instrument for researchers and policymakers alike.

Numerical study on local scour characteristics around submarine pipelines in the Yellow River Delta silty sandy soil under waves and currents

Due to their high reliability and cost-efficiency,submarine pipelines are widely used in offshore oil and gas resource engineering.Due to the interaction of waves,currents,seabed,and pipeline structures,the soil around submarine pipelines is prone to local scour,severely affecting their operational safety.With the Yellow River Delta as the research area and based on the renormalized group(RNG)k-εturbulence model and Stokes fifth-order wave theory,this study solves the Navier-Stokes(N-S)equation using the finite difference method.The volume of fluid(VOF)method is used to describe the fluid-free surface,and a threedimensional numerical model of currents and waves-submarine pipeline-silty sandy seabed is established.The rationality of the numerical model is verified using a self-built waveflow flume.On this basis,in this study,the local scour development and characteristics of submarine pipelines in the Yellow River Delta silty sandy seabed in the prototype environment are explored and the influence of the presence of pipelines on hydrodynamic features such as surrounding flow field,shear stress,and turbulence intensity is analyzed.The results indicate that(1)local scour around submarine pipelines can be divided into three stages:rapid scour,slow scour,and stable scour.The maximum scour depth occurs directly below the pipeline,and the shape of the scour pits is asymmetric.(2)As the water depth decreases and the pipeline suspension height increases,the scour becomes more intense.(3)When currents go through a pipeline,a clear stagnation point is formed in front of the pipeline,and the flow velocity is positively correlated with the depth of scour.This study can provide a valuable reference for the protection of submarine pipelines in this area.

Flow field, sedimentation, and erosion characteristics around folded linear HDPE sheet sand fence: Numerical simulation study

Wind and sand hazards are serious in the Milan Gobi area of the Xinjiang section of the Korla Railway. In order to ensure the safe operation of railroads, there is a need for wind and sand protection in heavily sandy areas. The wind and sand flow in the region is notably bi-directional. To shield railroads from sand, a unique sand fence made of folded linear high-density polyethylene(HDPE) is used, aligning with the principle that the dominant wind direction is perpendicular to the fence. This study employed field observations and numerical simulations to investigate the effectiveness of these HDPE sand fences in altering flow field distribution and offering protection. It also explored how these fences affect the deposition and erosion of sand particles. Findings revealed a significant reduction in wind speed near the fence corner;the minimum horizontal wind speed on the leeward side of the first sand fence(LSF) decreased dramatically from 3 m/s to 0.64 m/s. The vortex area on the LSF markedly impacted horizontal wind speeds. Within the LSF, sand deposition was a primary occurrence. As wind speeds increased, the deposition zone shrank, whereas the positive erosion zone expanded. Close to the folded corners of the HDPE sand fence, there was a notable shift from the positive erosion zone to a deposition zone. Field tests and numerical simulations confirmed the high windproof efficiency(WE) and sand resistance efficiency(SE) in the HDPE sand fence. Folded linear HDPE sheet sand fence can effectively slow down the incoming flow and reduce the sand content, thus achieving good wind and sand protection. This study provides essential theoretical guidance for the design and improvement of wind and sand protection systems in railroad engineering.

The Effects of Thickness and Location of PCMon the Building’s Passive Temperature-Control–A Numerical Study

Building energy consumption and building carbon emissions both account for more than 20%of their total national values in China.Building employing phase change material(PCM)for passive temperature control shows a promising prospect in meeting the comfort demand and reducing energy consumption simultaneously.However,there is a lack of more detailed research on the interaction between the location and thickness of PCM and indoor natural convection,as well as indoor temperature distribution.In this study,the numerical model of a passive temperature-controlled building integrating the developed PCM module is established with the help of ANSYS.In which,the actual weather condition of Beijing city is set as the boundary conditions and the indoor natural convection is simulated with the consideration of radiation model.The effects of PCM’s thickness and location on the internal temperature field are analyzed and discussed.The results show that the room could maintain within the human comfort temperature range with the longest ratio of 94.10%and the shortest ratio of 51.04%as integrating PCM.In comparison,the value is only 26.70%without PCM.The room’s maximum temperature fluctuation can also be improved;it could be lowered by 64.4%compared to the normal condition.When the quantity of PCM is sufficient,further increasing the PCM amount results in a temperature fluctuation reduction of less than 0.1°C and does not increase the comfort time.Placing PCM on the wall induces an apparent variation in indoor temperature along the vertical direction.Conversely,placing PCM on the roof can lead to a heat transfer rate difference of up to seven times.The optimal placement of PCM depends on the difference between the environmental and phase change temperatures.If the difference is positive,placing PCM on the roof is more effective;conversely,the opposite holds.According to the results over the entire cycle,PCM application on vertical walls yields better performance.The significant difference in natural convection caused by the same thickness of PCM but different application positions,coupled with the influence of air movement on the melting and solidification of PCM,further impacts indoor temperature fluctuations and comfort.This study can provide guidance for the application location and thickness of PCM,especially for scenarios where temperature regulation is required at a specific time.

Numerical study on gas production via a horizontal well from hydrate reservoirs with different slope angles in the South China Sea

It is important to study the effect of hydrate production on the physical and mechanical properties of low-permeability clayey–silty reservoirs for the largescale exploitation of hydrate reservoirs in the South China Sea.In this study,a multiphysical-field coupling model,combined with actual exploration drilling data and the mechanical experimental data of hydrate cores in the laboratory,was established to investigate the physical and mechanical properties of low-permeability reservoirs with different slope angles during 5-year hydrate production by the depressurization method via a horizontal well.The result shows that the permeability of reservoirs severely affects gas production rate,and the maximum gas production amount of a 20-m-long horizontal well can reach186.8 m3/day during the 5-year hydrate production.Reservoirs with smaller slope angles show higher gas production rates.The depressurization propagation and hydrate dissociation mainly develop along the direction parallel to the slope.Besides,the mean effective stress of reservoirs is concentrated in the near-wellbore area with the on-going hydrate production,and gradually decreases with the increase of the slope angle.Different from the effective stress distribution law,the total reservoir settlement amount first decreases and then increases with the increase of the slope angle.The maximum settlement of reservoirs with a 0°slope angle is up to 3.4 m,and the displacement in the near-wellbore area is as high as2.2 m after 5 years of hydrate production.It is concluded that the pore pressure drop region of low-permeability reservoirs in the South China Sea is limited,and various slope angles further lead to differences in effective stress and strain of reservoirs during hydrate production,resulting in severe uneven settlement of reservoirs.

Numerical and theoretical study of large-scale failure of strata overlying sublevel caving mines with steeply dipping discontinuities

The deformation and fracture evolution mechanisms of the strata overlying mines mined using sublevel caving were studied via numerical simulations.Moreover,an expression for the normal force acting on the side face of a steeply dipping superimposed cantilever beam in the surrounding rock was deduced based on limit equilibrium theory.The results show the following:(1)surface displacement above metal mines with steeply dipping discontinuities shows significant step characteristics,and(2)the behavior of the strata as they fail exhibits superimposition characteristics.Generally,failure first occurs in certain superimposed strata slightly far from the goaf.Subsequently,with the constant downward excavation of the orebody,the superimposed strata become damaged both upwards away from and downwards toward the goaf.This process continues until the deep part of the steeply dipping superimposed strata forms a large-scale deep fracture plane that connects with the goaf.The deep fracture plane generally makes an angle of 12°-20°with the normal to the steeply dipping discontinuities.The effect of the constant outward transfer of strata movement due to the constant outward failure of the superimposed strata in the metal mines with steeply dipping discontinuities causes the scope of the strata movement in these mines to be larger than expected.The strata in the metal mines with steeply dipping discontinuities mainly show flexural toppling failure.However,the steeply dipping structural strata near the goaf mainly exhibit shear slipping failure,in which case the mechanical model used to describe them can be simplified by treating them as steeply dipping superimposed cantilever beams.By taking the steeply dipping superimposed cantilever beam that first experiences failure as the key stratum,the failure scope of the strata(and criteria for the stability of metal mines with steeply dipping discontinuities mined using sublevel caving)can be obtained via iterative computations from the key stratum,moving downward toward and upwards away from the goaf.

Erosion wear at the bend of pipe during tailings slurry transportation:Numerical study considering inlet velocity,particle size and bend angle

Pipeline hydraulic transport is a highly efficient and low energy-consumption method for transporting solids and is commonly used for tailing slurry transport in the mining industry.Erosion wear(EW)remains the main cause of failure in tailings slurry pipeline systems,particularly at bends.EW is a complex phenomenon influenced by numerous factors,but research in this area has been limited.This study performs numerical simulations of slurry transport at the bend by combining computational fluid dynamics and fluid particle tracking using a wear model.Based on the validation of the feasibility of the model,this work focuses on the effects of coupled inlet velocity(IV)ranging from 1.5 to 3.0 m·s^(-1),particle size(PS)ranging from 50 to 650μm,and bend angle(BA)ranging from 45°to 90°on EW at the bend in terms of particle kinetic energy and incidence angle.The results show that the maximum EW rate of the slurry at the bend increases exponentially with IV and PS and first increases and then decreases with the increase in BA with the inflection point at 60°within these parameter ranges.Further comprehensive analysis reveals that the sensitivity level of the three factors to the maximum EW rate is PS>IV>BA,and when IV is 3.0 m/s,PS is 650μm,and BA is 60°,the bend EW is the most severe,and the maximum EW rate is 5.68×10^(-6)kg·m^(-2)·s^(-1).In addition,When PS is below or equal to 450μm,the maximum EW position is mainly at the outlet of the bend.When PS is greater than 450μm,the maximum EW position shifts toward the center of the bend with the increase in BA.Therefore,EW at the bend can be reduced in practice by reducing IV as much as possible and using small particles.

A three-dimensional numerical study on the effect of geometric asymmetry on arcjet thruster performance

In an arcjet thruster,the cathode and constrictor degrade with time,and the electrical arc discharge may become unsymmetrical.In this work,a three-dimensional numerical model of a hydrogen plasma arcjet is developed and validated to study the effect of unsymmetrical electric arc discharge on thruster performance.The unsymmetrical arc discharge is realized by introducing a radial shift of the cathode so that the cathode tip offset is 80μm(25%of the constrictor radius).Simulations are conducted for both axially centered cathode(coaxial)and off-centered cathode(non-coaxial)configurations with identical propellant flow rates and input current.Simulations show asymmetrical arc discharge in the non-coaxial cathode configuration,resulting in azimuthally asymmetric Joule heating,species concentrations,and velocity field.This asymmetry continues as the plasma expands in the divergent section of the nozzle.Temperature,species concentrations,and axial velocity exhibit asymmetric radial distribution at the nozzle exit.The computed Joule heating was found to reduce with cathode shift,and consequently,the thrust and specific impulse of the thruster was decreased by about 6.6%.In the case of the non-coaxial cathode,geometric asymmetry also induces a small side thrust.

Risk assessment and landslide prevention design using numerical modelling——A case study in Qingliu,China

Numerous Quaternary deposits are existed in the mountainous areas of Southwest China,especially in the transition zone between the QinghaiTibet Plateau and the Sichuan Basin,where strong tectonic movements and frequent climatic changes increase the potential landslides.The possible deformation and failure process of potential landslides and their impacts on the surrounding environment are important research topics.Field investigation and monitoring indicate that the Qingliu landslide in Xiameng town,Li County,Sichuan Province,China has been continuously deforming since August 2020.The deformation zone has a maximum deformation depth of approximately 18.9m,a total area of 54,628 m2,and a volume of 34.0×104 m3,which seriously threatens infrastructure projects and dwellings.As a result,understanding the Qingliu landslide evolution process,assessing the hazard risk,and planning disaster prevention measures are of great significance for reducing disaster loss.In this study,the mass movement process and hazard risk of the Qingliu landslide are evaluated,and the effects of different prevention measures are compared and discussed.By using the depth-integrated method,the mass movement of the Qingliu landslide is analyzed.The numerical simulation results indicate that the maximum velocity of the Qingliu landslide is approximately 37.5 m/s,and the duration of the landslide is approximately 90s.The simulated landslide can eventually form a deposited mass with a maximum deposit thickness of 19.4 m and an area of approximately 60,168.3 m2,thereby blocking the river and burying dwellings.Furthermore,a risk assessment of the Qingliu landslide under different forms of protection measures is also produced and discussed by considering the hazard level and economic vulnerability level of the affected area.Setting three layers of anti-slide piles on the deformation zone to reduce the hazard risk of the Qingliu landslide is a better choice.Our results may be useful for planning prevention measures and improving disaster emergency response systems.

A numerical case study on the long-term seismic assessment of reinforced concrete tunnels in corrosive environments

The paper investigates the long-term seismic behaviour of an underground reinforced concrete(RC)metro tunnel in Santiago,Chile,considering the combined effects of chloride-induced corrosion and cumulative,low-amplitude seismic shaking on the structure’s performance.The soil-tunnel response is evaluated with the aid of transient,nonlinear finite element analysis using a two-dimensional(2D)plane strain numerical model that adopts advanced nonlinear models for the simulation of soil and concrete plasticity and the dynamic stiffness behaviour.The effects of corrosion deterioration are demonstrated in terms of time-dependent loss of rebar area and cover concrete stiffness and strength.The study illustrates the influence of ageing and repeated seismic shaking on lining deformation,crack development,and the modal characteristics of the intact and degrading systems.The results indicate that multiple lowamplitude events drive the non-degrading RC tunnel beyond its elastic regime without significant structural response consequences.A noticeable impact of corrosion deterioration on the structure’s seismic performance is revealed,increasing with the number and intensity of earthquake events.Two different tunnel embedment depths are comparatively assessed.The analyses demonstrate larger coseismic section convergence in the case of the deeper tunnel,yet a less pronounced effect of ageing and successive seismic loading compared to the shallow section,which is evident in the RC lining cracks at the end of shaking.

Numerical simulation study on shear resistance of anchorage joints considering tensile-shear fracture criterion of 2G-NPR bolt

2G-NPR bolt (the 2nd generation Negative Poisson’s Ratio bolt) is a new type of bolt with high strength, high toughness and no yield platform. It has signifcant efects on improving the shear strength of jointed rock mass and controlling the stability of surrounding rock. To achieve an accurate simulation of bolted joint shear tests, we have studied a numerical simulation method that takes into account the 2G-NPR bolt's tensile–shear fracture criterion. Firstly, the indoor experimental study on the tensile–shear mechanical properties of 2G-NPR bolt is carried out to explore its mechanical properties under diferent tensile–shear angles, and the fracture criterion of 2G-NPR bolt considering the tensile–shear angle is established. Then, a three-dimensional numerical simulation method considering the tensile–shear mechanical constitutive and fracture criterion of 2G-NPR bolt, the elastoplastic mechanical behavior of surrounding rock and the damage and deterioration of grouting body is proposed. The feasibility and accuracy of the method are verifed by comparing with the indoor shear test results of 2G-NPR bolt anchorage joints. Finally, based on the numerical simulation results, the deformation and stress of the bolt, the distribution of the plastic zone of the rock mass, the stress distribution and the damage of the grouting body are analyzed in detail. The research results can provide a good reference value for the practical engineering application and shear mechanical performance analysis of 2G-NPR bolt.

Numerical Simulation of Thermo-Hydro-Mechanical-Chemical Response Caused by CO_(2)Injection into Saline Geological Formations:A Case Study from the Ordos Project,China

Thermo-hydro-mechanical-chemical(THMC)interactions are prevalent during CO_(2)geological sequestration(CGS).In this study,a sequential coupling THMC numerical simulation program was constructed,which can be used to explore the following issues of CGS:fluid and heat flow,solute transport;stresses,displacements and rock failures related to geo-mechanical effects;equilibrium and kinetic chemical reactions;chemical damage to mechanical properties of the rock.Then,the coupling program was applied to the Ordos CGS Project to study the formation response under the multifield interaction caused by CO_(2)injection.The simulation results show that the mechanical process dominates the short CO_(2)injection period.Specifically,the formation’s permeability near the injection well increases by 43%,due to the reduction of effective stress,which significantly promotes the lateral migration of CO_(2).When the injection rate exceeds 0.15 million tons per year,the cohesion of the reservoir rock is not enough to resist the shear force inside the rock and rock failure may occur.During the subsequent long-term sequestration period(200 years),the influence of mineral reactions gradually increases.Due to calcite dissolution,the shear modulus of caprock is predicted to decrease by 7.6%,which will to some extent increase the risk of rock failure.

A Numerical Simulation Study of Temperature and Pressure Effects on the Breakthrough Pressure of CO_(2)in Unsaturated Low-permeability Rock Core

Gas breakthrough pressure is a key parameter to evaluate the sealing capacity of caprock,and it also plays important roles in safety and capacity of CO_(2)geological storage.Based on the published experimental results,we present numerical simulations on CO_(2)breakthrough pressure in unsaturated low-permeability rock under 9 multiple P-T conditions(which can keep CO_(2)in gaseous,liquid and supercritical states)and thus,a numerical method which can be used to accurately predict CO_(2)breakthrough pressure on rock-core scale is proposed.The simulation results show that CO_(2)breakthrough pressure and breakthrough time are exponential correlated with P-T conditions.Meanwhile,pressure has stronger effects on experimental results than that of temperature.Moreover,we performed sensitivity studies on the pore distribution indexλ(0.6,0.7,0.8,and 0.9)in van Genuchten-Muale model.Results show that with the increase ofλ,CO_(2)breakthrough pressure and breakthrough time both show decreasing trends.In other words,the larger the value ofλis,the better the permeability of the caprock is,and the worse the CO_(2)sealing capacity is.The numerical method established in this study can provide an important reference for the prediction of gas breakthrough pressure on rock-core scale and for related numerical studies.

Flow-Induced Clogging in Microfiltration Membranes: Numerical Modeling and Parametric Study

Microfiltration membrane technology has been widely used in various industries for solid-liquid separation. However, pore clogging remains a persistent challenge. This study employs (CFD) and discrete element method (DEM) models to enhance our understanding of microfiltration membrane clogging. The models were validated by comparing them to experimental data, demonstrating reasonable consistency. Subsequently, a parametric study was conducted on a cross-flow model, exploring the influence of key parameters on clogging. Findings show that clogging is a complex phenomenon affected by various factors. The mean inlet velocity and transmembrane flux were found to directly impact clogging, while the confinement ratio and cosine of the membrane pore entrance angle had an inverse relationship with it. Two clog types were identified: internal (inside the pore) and external (arching at the pore entrance), with the confinement ratio determining the type. This study introduced a dimensionless number as a quantitative clogging indicator based on transmembrane flux, Reynolds number, filtration time, entrance angle cosine, and confinement ratio. While this hypothesis held true in simulations, future studies should explore variations in clogging indicators, and improved modeling of clogging characteristics. Calibration between numerical and physical times and consideration of particle volume fraction will enhance understanding.