To fundamentally alleviate the excavation chamber clogging during slurry tunnel boring machine(TBM)advancing in hard rock,large-diameter short screw conveyor was adopted to slurry TBM of Qingdao Jiaozhou Bay Second Un...To fundamentally alleviate the excavation chamber clogging during slurry tunnel boring machine(TBM)advancing in hard rock,large-diameter short screw conveyor was adopted to slurry TBM of Qingdao Jiaozhou Bay Second Undersea Tunnel.To evaluate the discharging performance of short screw conveyor in different cases,the full-scale transient slurry-rock two-phase model for a short screw conveyor actively discharging rocks was established using computational fluid dynamics-discrete element method(CFD-DEM)coupling approach.In the fluid domain of coupling model,the sliding mesh technology was utilized to describe the rotations of the atmospheric composite cutterhead and the short screw conveyor.In the particle domain of coupling model,the dynamic particle factories were established to produce rock particles with the rotation of the cutterhead.And the accuracy and reliability of the CFD-DEM simulation results were validated via the field test and model test.Furthermore,a comprehensive parameter analysis was conducted to examine the effects of TBM operating parameters,the geometric design of screw conveyor and the size of rocks on the discharging performance of short screw conveyor.Accordingly,a reasonable rotational speed of screw conveyor was suggested and applied to Jiaozhou Bay Second Undersea Tunnel project.The findings in this paper could provide valuable references for addressing the excavation chamber clogging during ultra-large-diameter slurry TBM tunneling in hard rock for similar future.展开更多
A two-phase wedge-sliding model is developed based on the micro-cellular structure and minimum entropy theory of a stable system, and it is used to describe the ingredient distribution of a mixed fluid in a non-unifor...A two-phase wedge-sliding model is developed based on the micro-cellular structure and minimum entropy theory of a stable system, and it is used to describe the ingredient distribution of a mixed fluid in a non-uniform stress field and to analyse its phase drift phenomenon. In the model, the drift-inhibition angle and the expansion-inhibition angle are also deduced and used as evaluating indexes to describe the drifting trend of different ingredients among the mixed fluids. For solving above two indexes of the model, a new calculation method is developed and used to compute the phase distributions of multiphase fluid at peak stress and gradient area stress, respectively. As an example, the flow process of grease in a pipe is analysed by simulation method and used to verify the validity of the model.展开更多
This work aims at comparing surface tension models in VOF(Volume of Fluid) modeling and investigating the effects of gas distributor and gas velocity. Hydrodynamics of a continuous chain of bubbles inside a bubble col...This work aims at comparing surface tension models in VOF(Volume of Fluid) modeling and investigating the effects of gas distributor and gas velocity. Hydrodynamics of a continuous chain of bubbles inside a bubble column reactor was simulated. The grid independence study was first conducted and a grid size of 1.0 mm was adopted in order to minimize the computing time without compromising the accuracy of the results. The predictions were validated by comparing the experimental studies reported in the literature. It was found that all surface tension models can describe the bubble rise and bubble plume in a column with slight deviations.展开更多
Ice particles could form under the continuous impingement of incoming supercooled droplets in icing conditions,which will change the surface roughness to enhance the further heat and mass transfer during icing process...Ice particles could form under the continuous impingement of incoming supercooled droplets in icing conditions,which will change the surface roughness to enhance the further heat and mass transfer during icing process.A fixed-grid porous enthalpy method based on the improved Discrete Phase Model(DPM)and Volume of Fluid(VOF)integrated algorithm is developed to solve the multiphase heat transfer problem to give more detailed demonstration of the formation of initial ice roughness.The algorithms to determine the criterion of transformation from DPM to VOF and the allocation of source items during transformation are improved to the general DPM-VOF algorithm.Two verification cases,namely two glycerine-solution droplets impact and single droplet freeze,are conducted to verify the accuracy and reliability of the enthalpy-DPMVOF method,where the simulation results match well with experiment phenomena.Ice roughness on a NACA0012 airfoil is precisely captured and the effects on convective heat transfer characteristics are preliminarily revealed.The results illustrate that the enthalpy-DPM-VOF method could successfully capture the characteristics of motion and the phase change process of droplet,as well as balance the calculation accuracy and efficiency.展开更多
Liquid permeability of the mushy zone is important for porosity formation during the solidification process. In order to investigate the permeability of the mushy zone, an integrated model was developed by combining t...Liquid permeability of the mushy zone is important for porosity formation during the solidification process. In order to investigate the permeability of the mushy zone, an integrated model was developed by combining the phase field model and computational fluid dynamics (CFD) model. The three-dimensional multigrain dendrite morphology was obtained by using the phase field model. Subsequently, the computer-aided design (CAD) geometry and mesh were generated based on calculated dendrite morphologies. Finally, the permeability of the dendritic mushy zone was obtained by solving the Navier-Stokes and continuity equations in ANSYS Fluent software. As an example, the dendritic mushy zone permeability of Al-4.5wt%Cu alloy and its relationship with the solid fractions were studied in detail. The predicted permeability data can be input to the solidification model on a greater length scale for macro segregation and porosity simulations.展开更多
In the industrial process of producing the strong phosphoric acid(SPA),clarification of the solution is essential to the ultimate product.However,the large viscosity of sediment and the induced interface interaction r...In the industrial process of producing the strong phosphoric acid(SPA),clarification of the solution is essential to the ultimate product.However,the large viscosity of sediment and the induced interface interaction result in difficulties when the SPA is clarified.CFD numerical methodology was applied to simulate internal flow field and performance of the low speed scraper based on Mixture solidliquid two-phase flow model.Sediment deposition was generated by loading solid particles at the bottom of clarifying vessel.The moving mesh and RNG k-εmodel were used to simulate the rotational turbulent flow in clarifying tank.Variables studied,amongst others,were the scraper rotation speed and the mounting height,which could affect the solid suspension height.Features of flow field and solid volume fraction distribution in computational domain were presented and analyzed.The numerical reports of the scraper torque and velocities of inlet and outlet filed were obtained.It seems the torque value of rotatio-nal axis and particle suspending height augment with an increasing rotating speed.Meanwhile,a high revolving speed is good for the deposition discharge.The particle fraction distribution in meridional surface and horizontal surface at fixed rotation speed were analyzed to determine the corresponding optimal installation height.The simulating results reflect the flow field is marginally stirred by the scraper and proper working parameters are obtained,in which case the comprehensive properties of the scraper and the clarifying tank are superior.展开更多
For solid-fluid interaction, one of the phase-density equations in diffuse interface models is degenerated to a "0 = 0" equation when the volume fraction of a certain phase takes the value of zero or unity. ...For solid-fluid interaction, one of the phase-density equations in diffuse interface models is degenerated to a "0 = 0" equation when the volume fraction of a certain phase takes the value of zero or unity. This is because the conservative variables in phasedensity equations include volume fractions. The degeneracy can be avoided by adding an artificial quantity of another material into the pure phase. However, nonphysical waves,such as shear waves in fluids, are introduced by the artificial treatment. In this paper,a transport diffuse interface model, which is able to treat zero/unity volume fractions, is presented for solid-fluid interaction. In the proposed model, a new formulation for phase densities is derived, which is unrelated to volume fractions. Consequently, the new model is able to handle zero/unity volume fractions, and nonphysical waves caused by artificial volume fractions are prevented. One-dimensional and two-dimensional numerical tests demonstrate that more accurate results can be obtained by the proposed model.展开更多
Understanding the dynamics of phase boundaries in fluids requires quantitative knowledge about the microscale processes at the interface.We consider the sharp-interface motion of the compressible two-component flow an...Understanding the dynamics of phase boundaries in fluids requires quantitative knowledge about the microscale processes at the interface.We consider the sharp-interface motion of the compressible two-component flow and propose a heterogeneous multiscale method(HMM)to describe the flow fields accurately.The multiscale approach combines a hyperbolic system of balance laws on the continuum scale with molecular-dynamics(MD)simulations on the microscale level.Notably,the multiscale approach is necessary to compute the interface dynamics because there is—at present—no closed continuum-scale model.The basic HMM relies on a moving-mesh finite-volume method and has been introduced recently for the compressible one-component flow with phase transitions by Magiera and Rohde in(J Comput Phys 469:111551,2022).To overcome the numerical complexity of the MD microscale model,a deep neural network is employed as an efficient surrogate model.The entire approach is finally applied to simulate droplet dynamics for argon-methane mixtures in several space dimensions.To our knowledge,such compressible two-phase dynamics accounting for microscale phase-change transfer rates have not yet been computed.展开更多
Climate change is a reality. The burning of fossil fuels from oil, natural gas and coal is responsible for much of the pollution and the increase in the planet’s average temperature, which has raised discussions on t...Climate change is a reality. The burning of fossil fuels from oil, natural gas and coal is responsible for much of the pollution and the increase in the planet’s average temperature, which has raised discussions on the subject, given the emergencies related to climate. An energy transition to clean and renewable sources is necessary and urgent, but it will not be quick. In this sense, increasing the efficiency of oil extraction from existing sources is crucial, to avoid waste and the drilling of new wells. The purpose of this work was to add diffusive and dispersive terms to the Buckley-Leverett equation in order to incorporate extra phenomena in the temporal evolution between the water-oil and oil-water transitions in the pipeline. For this, the modified Buckley-Leverett equation was discretized via essentially weighted non-oscillatory schemes, coupled with a three-stage Runge-Kutta and a fourth-order centered finite difference methods. Then, computational simulations were performed and the results showed that new features emerge in the transitions, when compared to classical simulations. For instance, the dispersive term inhibits the diffusive term, adding oscillations, which indicates that the absorption of the fluid by the porous medium occurs in a non-homogeneous manner. Therefore, based on research such as this, decisions can be made regarding the replacement of the porous medium or the insertion of new components to delay the replacement.展开更多
The value of form factor k at different drafts is important in predicting full-scale total resistance and speed for different types of ships. In the ITTC community, most organizations predict form factor k using a low...The value of form factor k at different drafts is important in predicting full-scale total resistance and speed for different types of ships. In the ITTC community, most organizations predict form factor k using a low-speed model test. However, this method is problematic for ships with bulbous bows and transom. In this article, a Computational Fluid Dynamics(CFD)-based method is introduced to obtain k for different type of ships at different drafts, and a comparison is made between the CFD method and the model test. The results show that the CFD method produces reasonable k values. A grid generating method and turbulence model are briefly discussed in the context of obtaining a consistent k using CFD.展开更多
The distribution of proppant injected in hydraulic fractures significantly affects the fracture conductivity and well performance.The proppant transport in thin fracturing fluid used during hydraulic fracturing in the...The distribution of proppant injected in hydraulic fractures significantly affects the fracture conductivity and well performance.The proppant transport in thin fracturing fluid used during hydraulic fracturing in the unconventional reservoirs is considerably different from fracturing fluids in the conventional reservoir due to the very low viscosity and quick deposition of the proppants.This paper presents the development of a three-dimensional Computational Fluid Dynamics(CFD)modelling technique for the prediction of proppant-fluid multiphase flow in hydraulic fractures.The proposed model also simulates the fluid leak-off behaviour from the fracture wall.The Euler-Granular and CFD-Discrete Element Method(CFD-DEM)multiphase modelling approach has been applied,and the equations defining the fluid-proppant and inter-proppant interaction have been solved using the finite volume technique.The proppant transport in hydraulic fractures has been studied comprehensively,and the computational modelling results of proppant distribution and other flow properties are in good agreement with the published experimental study.The parametric study is performed to investigate the effect of variation in proppant size,fluid viscosity and fracture width on the proppant transport.Smaller proppants can be injected early,followed by larger proppants to maintain high propping efficiency.This study has enhanced the understanding of the complex flow phenomenon between proppant and fracturing fluid and can play a vital role in hydraulic fracturing design.展开更多
In the present study, considering the transport and transformation processes of variables, a threedimensional water quality model for the river system was established, which coupled the volume of fluid(VOF) method wit...In the present study, considering the transport and transformation processes of variables, a threedimensional water quality model for the river system was established, which coupled the volume of fluid(VOF) method with the k-ε turbulence mathematical model. Then, the water hydrodynamic characteristics and transport processes for BOD_5, NH_(3^-)N and TP were analyzed. The results showed that the water surface of convex bank was a little lower than that of concave bank due to the centrifugal force near the bend, and most concentrations were inferior to the type Ⅴ standard indexes of surface water environmental quality. The model validation indicated that the errors between the simulated and monitored values were comparatively small, satisfying the application demands and providing scientific basis and decision support for the restoration and protection of water quality.展开更多
In this paper, the lattice-Bohzmann method is used to investigate the droplet dynamics after impact on horizontal and inclined solid surface. The two-phase interparticle potential model is employed. The model is found...In this paper, the lattice-Bohzmann method is used to investigate the droplet dynamics after impact on horizontal and inclined solid surface. The two-phase interparticle potential model is employed. The model is found to possess a linear relation between the macroscopic properties （ surface tension σ and contact angle α） and microscopic parameters （ G, G, ）. The flow state of the droplet on the surface is analyzed in detail, and the effects of surface characteristic, impact velocity, impact angle, the viscosity and surface tension of the liquid are investigated, respectively. It is shown that the lattice-Bohzmann method can not only track exactly and automatically the interface, but also the simulation results have a good qualitative agreement with ones of the previous experimental and numerical studies.展开更多
Nutrients supply especially like nutrients and oxygen play vital role in tissue engineering process.It is found that tissue could not grow very well in the middle of the scaffold because few nutrients could transport ...Nutrients supply especially like nutrients and oxygen play vital role in tissue engineering process.It is found that tissue could not grow very well in the middle of the scaffold because few nutrients could transport to the middle.Nutrient limitations would reduce cell proliferation and differentiation.In that case,there is urgent need to understand the nutrient distribution for both in vitro and in vivo study,as no technology is able for researchers to observe the nutrients transport during those process.In this paper,a numerical model coupling with VOF(volume of fluid)model and species transport model together for predicting the distribution of oxygen and glucose in the scaffold after implantation in to the site is developed.Comparing with our previous in vivo tests,the regenerated tissue distribution has a similar trend as oxygen distribution rather than glucose.The reported scaffold manufactured by additive manufacturing provided a good interconnected structure which facilitated the nutrient transportation in the scaffold.Considering nutrient transportation,this numerical model could be used in better understanding the nutrients transportation in the scaffold,and leading to a better understanding of tissue formation in the scaffold.展开更多
Iterative coupled methods are widely used in multi-fidelity simulation of rotating components due to the simple implementation,which iteratively eliminates the errors between the computational fluid dynamics models an...Iterative coupled methods are widely used in multi-fidelity simulation of rotating components due to the simple implementation,which iteratively eliminates the errors between the computational fluid dynamics models and approximate characteristic maps.However,the convergence and accuracy of the iterative coupled method are trapped in characteristic maps.In particular,iterative steps increase sharply as the operation point moves away from the design point.To address these problems,this paper developed an auxiliary iterative coupled method that introduces the static-pressure-auxiliary characteristic maps and modification factor of mass flow into the component-level model.The developed auxiliary method realized the direct transfer of static pressure between the high-fidelity models and the component-level model.Multi-fidelity simulations of the throttle characteristics were carried out using both the auxiliary and traditional iterative coupled methods,and the simulation results were verified using the experimental data.Additionally,the consistency between the auxiliary and traditional iterative coupled methods was confirmed.Subsequently,multi-fidelity simulations of the speed and altitude characteristics were also conducted.The auxiliary and traditional iterative coupled methods were evaluated in terms of convergence speed and accuracy.The evaluation indicated that the auxiliary iterative coupled method significantly reduces iterative steps by approximately 50%at the near-choked state.In general,the auxiliary iterative coupled method is preferred as a development of the traditional iterative coupled method in the near-choked state,and the combined auxiliary-traditional iterative coupled method provides support for successful multi-fidelity simulation in far-off-design conditions.展开更多
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 (...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.展开更多
【目的】棒销式砂磨机的工作过程存在复杂的气体、液体和固体的三相耦合现象,为了使仿真工况更接近实际工况,构建气-液-固三相耦合仿真模型,提高棒销式砂磨机仿真设计的准确性。【方法】采用离散单元法(discrete element method,DEM)和...【目的】棒销式砂磨机的工作过程存在复杂的气体、液体和固体的三相耦合现象,为了使仿真工况更接近实际工况,构建气-液-固三相耦合仿真模型,提高棒销式砂磨机仿真设计的准确性。【方法】采用离散单元法(discrete element method,DEM)和计算流体动力学(computational fluid dynamics,CFD)分别研究固体相和流体相,并引入流体体积模型(volume of fluid model,VOF)区分流体相所包含的液体相和气体相,分析固体相运动方程、流体相控制方程,确定气-液界面的识别方法和耦合计算方法,制定仿真流程;通过单球落水仿真、颗粒群落水仿真试验分析CFD-DEM-VOF三相耦合模型仿真计算的精度,并进行准确性验证;在设置仿真参数、进行网格划分及其无关性分析基础上,针对棒销式砂磨机的CFD-DEM-VOF三相耦合模型进行仿真试验;对流体速度、颗粒总能量和颗粒的速度的仿真结果进行分析,并通过实验验证仿真结果。【结果】在单球落水仿真试验中,根据CFD-DEM-VOF三相耦合模型的仿真结果与根据Stokes定律的理论计算结果基本吻合;在颗粒群落水仿真过程中,液面上升高度的仿真值与理论值之间的相对误差为1.37%,VOF模型的体积守恒性较好;棒销四面体网格边长小于2 mm、研磨桶四面体网格边长小于2.5 mm时,满足网格独立性的精度要求,同时计算量也较少;随着棒销转速的增大,流体速度、颗粒总能量、颗粒平均速度也逐渐增大;当棒销转速为1400~2000 r/min时,CFD-DEM-VOF三相耦合模型流体速度的仿真与实验结果最为接近;当棒销转速为1400~2200 r/min时,CFD-DEM-VOF三相耦合模型的颗粒总能量仿真值与实验值的最大相对误差为1%。【结论】与仅仅采用流体相、固体相单相模型或固-液两相模型相比,采用CFD-DEM-VOF三相耦合模型设计棒销式砂磨机的计算精度和准确性较高,仿真性能好。展开更多
基金supported by the Fundamental Research Funds for the Central Universities(Grant No.2023YJS053)the National Natural Science Foundation of China(Grant No.52278386).
文摘To fundamentally alleviate the excavation chamber clogging during slurry tunnel boring machine(TBM)advancing in hard rock,large-diameter short screw conveyor was adopted to slurry TBM of Qingdao Jiaozhou Bay Second Undersea Tunnel.To evaluate the discharging performance of short screw conveyor in different cases,the full-scale transient slurry-rock two-phase model for a short screw conveyor actively discharging rocks was established using computational fluid dynamics-discrete element method(CFD-DEM)coupling approach.In the fluid domain of coupling model,the sliding mesh technology was utilized to describe the rotations of the atmospheric composite cutterhead and the short screw conveyor.In the particle domain of coupling model,the dynamic particle factories were established to produce rock particles with the rotation of the cutterhead.And the accuracy and reliability of the CFD-DEM simulation results were validated via the field test and model test.Furthermore,a comprehensive parameter analysis was conducted to examine the effects of TBM operating parameters,the geometric design of screw conveyor and the size of rocks on the discharging performance of short screw conveyor.Accordingly,a reasonable rotational speed of screw conveyor was suggested and applied to Jiaozhou Bay Second Undersea Tunnel project.The findings in this paper could provide valuable references for addressing the excavation chamber clogging during ultra-large-diameter slurry TBM tunneling in hard rock for similar future.
基金Project supported by the National Natural Science Foundation of China (Grant No. 51075311)
文摘A two-phase wedge-sliding model is developed based on the micro-cellular structure and minimum entropy theory of a stable system, and it is used to describe the ingredient distribution of a mixed fluid in a non-uniform stress field and to analyse its phase drift phenomenon. In the model, the drift-inhibition angle and the expansion-inhibition angle are also deduced and used as evaluating indexes to describe the drifting trend of different ingredients among the mixed fluids. For solving above two indexes of the model, a new calculation method is developed and used to compute the phase distributions of multiphase fluid at peak stress and gradient area stress, respectively. As an example, the flow process of grease in a pipe is analysed by simulation method and used to verify the validity of the model.
基金Supported by the National Ministry of Science and Technology of China(2017YFB0602401)the National Natural Science Foundation of China(21776173,91834303,U1862201,21625603)the Program of Shanghai Subject Chief Scientists(18XD1402000).
文摘This work aims at comparing surface tension models in VOF(Volume of Fluid) modeling and investigating the effects of gas distributor and gas velocity. Hydrodynamics of a continuous chain of bubbles inside a bubble column reactor was simulated. The grid independence study was first conducted and a grid size of 1.0 mm was adopted in order to minimize the computing time without compromising the accuracy of the results. The predictions were validated by comparing the experimental studies reported in the literature. It was found that all surface tension models can describe the bubble rise and bubble plume in a column with slight deviations.
基金supported by the National Natural Science Foundation of China(No.51706244)National Science and Technology Major Projects of China(No.2017-VIII-0003-0114)。
文摘Ice particles could form under the continuous impingement of incoming supercooled droplets in icing conditions,which will change the surface roughness to enhance the further heat and mass transfer during icing process.A fixed-grid porous enthalpy method based on the improved Discrete Phase Model(DPM)and Volume of Fluid(VOF)integrated algorithm is developed to solve the multiphase heat transfer problem to give more detailed demonstration of the formation of initial ice roughness.The algorithms to determine the criterion of transformation from DPM to VOF and the allocation of source items during transformation are improved to the general DPM-VOF algorithm.Two verification cases,namely two glycerine-solution droplets impact and single droplet freeze,are conducted to verify the accuracy and reliability of the enthalpy-DPMVOF method,where the simulation results match well with experiment phenomena.Ice roughness on a NACA0012 airfoil is precisely captured and the effects on convective heat transfer characteristics are preliminarily revealed.The results illustrate that the enthalpy-DPM-VOF method could successfully capture the characteristics of motion and the phase change process of droplet,as well as balance the calculation accuracy and efficiency.
基金financially supported by the National Key Research and Development Program of China(No.2016YFB0700503)National Natural Science Foundation of China(No.51701013)Beijing Laboratory of Metallic Materials and Processing for Modern Transportation
文摘Liquid permeability of the mushy zone is important for porosity formation during the solidification process. In order to investigate the permeability of the mushy zone, an integrated model was developed by combining the phase field model and computational fluid dynamics (CFD) model. The three-dimensional multigrain dendrite morphology was obtained by using the phase field model. Subsequently, the computer-aided design (CAD) geometry and mesh were generated based on calculated dendrite morphologies. Finally, the permeability of the dendritic mushy zone was obtained by solving the Navier-Stokes and continuity equations in ANSYS Fluent software. As an example, the dendritic mushy zone permeability of Al-4.5wt%Cu alloy and its relationship with the solid fractions were studied in detail. The predicted permeability data can be input to the solidification model on a greater length scale for macro segregation and porosity simulations.
基金Graduate Research and Innovation Program in Jiangsu Province(KYZZ16_0286)
文摘In the industrial process of producing the strong phosphoric acid(SPA),clarification of the solution is essential to the ultimate product.However,the large viscosity of sediment and the induced interface interaction result in difficulties when the SPA is clarified.CFD numerical methodology was applied to simulate internal flow field and performance of the low speed scraper based on Mixture solidliquid two-phase flow model.Sediment deposition was generated by loading solid particles at the bottom of clarifying vessel.The moving mesh and RNG k-εmodel were used to simulate the rotational turbulent flow in clarifying tank.Variables studied,amongst others,were the scraper rotation speed and the mounting height,which could affect the solid suspension height.Features of flow field and solid volume fraction distribution in computational domain were presented and analyzed.The numerical reports of the scraper torque and velocities of inlet and outlet filed were obtained.It seems the torque value of rotatio-nal axis and particle suspending height augment with an increasing rotating speed.Meanwhile,a high revolving speed is good for the deposition discharge.The particle fraction distribution in meridional surface and horizontal surface at fixed rotation speed were analyzed to determine the corresponding optimal installation height.The simulating results reflect the flow field is marginally stirred by the scraper and proper working parameters are obtained,in which case the comprehensive properties of the scraper and the clarifying tank are superior.
基金Project supported by the National Natural Science Foundation of China(Nos.11702029,11771054,U1730118,91852207,and 11801036)the China Postdoctoral Science Foundation(No.2016M600967)
文摘For solid-fluid interaction, one of the phase-density equations in diffuse interface models is degenerated to a "0 = 0" equation when the volume fraction of a certain phase takes the value of zero or unity. This is because the conservative variables in phasedensity equations include volume fractions. The degeneracy can be avoided by adding an artificial quantity of another material into the pure phase. However, nonphysical waves,such as shear waves in fluids, are introduced by the artificial treatment. In this paper,a transport diffuse interface model, which is able to treat zero/unity volume fractions, is presented for solid-fluid interaction. In the proposed model, a new formulation for phase densities is derived, which is unrelated to volume fractions. Consequently, the new model is able to handle zero/unity volume fractions, and nonphysical waves caused by artificial volume fractions are prevented. One-dimensional and two-dimensional numerical tests demonstrate that more accurate results can be obtained by the proposed model.
基金Funding Open Access funding enabled and organized by Projekt DEAL.When preparing this manuscript,the authors have kept the COPE guidelines on how to deal with potential acts of misconduct.The research leading to these results received funding from Deutsche Forschungsgemeinschaft(DFG,German Research Foundation)through the project SFB-TRR 75 with the project number 84292822the DFG under Germanys Excellence Strategy-EXC2075with the project number390740016.
文摘Understanding the dynamics of phase boundaries in fluids requires quantitative knowledge about the microscale processes at the interface.We consider the sharp-interface motion of the compressible two-component flow and propose a heterogeneous multiscale method(HMM)to describe the flow fields accurately.The multiscale approach combines a hyperbolic system of balance laws on the continuum scale with molecular-dynamics(MD)simulations on the microscale level.Notably,the multiscale approach is necessary to compute the interface dynamics because there is—at present—no closed continuum-scale model.The basic HMM relies on a moving-mesh finite-volume method and has been introduced recently for the compressible one-component flow with phase transitions by Magiera and Rohde in(J Comput Phys 469:111551,2022).To overcome the numerical complexity of the MD microscale model,a deep neural network is employed as an efficient surrogate model.The entire approach is finally applied to simulate droplet dynamics for argon-methane mixtures in several space dimensions.To our knowledge,such compressible two-phase dynamics accounting for microscale phase-change transfer rates have not yet been computed.
文摘Climate change is a reality. The burning of fossil fuels from oil, natural gas and coal is responsible for much of the pollution and the increase in the planet’s average temperature, which has raised discussions on the subject, given the emergencies related to climate. An energy transition to clean and renewable sources is necessary and urgent, but it will not be quick. In this sense, increasing the efficiency of oil extraction from existing sources is crucial, to avoid waste and the drilling of new wells. The purpose of this work was to add diffusive and dispersive terms to the Buckley-Leverett equation in order to incorporate extra phenomena in the temporal evolution between the water-oil and oil-water transitions in the pipeline. For this, the modified Buckley-Leverett equation was discretized via essentially weighted non-oscillatory schemes, coupled with a three-stage Runge-Kutta and a fourth-order centered finite difference methods. Then, computational simulations were performed and the results showed that new features emerge in the transitions, when compared to classical simulations. For instance, the dispersive term inhibits the diffusive term, adding oscillations, which indicates that the absorption of the fluid by the porous medium occurs in a non-homogeneous manner. Therefore, based on research such as this, decisions can be made regarding the replacement of the porous medium or the insertion of new components to delay the replacement.
基金Supported by Ministry of Industry and Information(No.K24097)
文摘The value of form factor k at different drafts is important in predicting full-scale total resistance and speed for different types of ships. In the ITTC community, most organizations predict form factor k using a low-speed model test. However, this method is problematic for ships with bulbous bows and transom. In this article, a Computational Fluid Dynamics(CFD)-based method is introduced to obtain k for different type of ships at different drafts, and a comparison is made between the CFD method and the model test. The results show that the CFD method produces reasonable k values. A grid generating method and turbulence model are briefly discussed in the context of obtaining a consistent k using CFD.
文摘The distribution of proppant injected in hydraulic fractures significantly affects the fracture conductivity and well performance.The proppant transport in thin fracturing fluid used during hydraulic fracturing in the unconventional reservoirs is considerably different from fracturing fluids in the conventional reservoir due to the very low viscosity and quick deposition of the proppants.This paper presents the development of a three-dimensional Computational Fluid Dynamics(CFD)modelling technique for the prediction of proppant-fluid multiphase flow in hydraulic fractures.The proposed model also simulates the fluid leak-off behaviour from the fracture wall.The Euler-Granular and CFD-Discrete Element Method(CFD-DEM)multiphase modelling approach has been applied,and the equations defining the fluid-proppant and inter-proppant interaction have been solved using the finite volume technique.The proppant transport in hydraulic fractures has been studied comprehensively,and the computational modelling results of proppant distribution and other flow properties are in good agreement with the published experimental study.The parametric study is performed to investigate the effect of variation in proppant size,fluid viscosity and fracture width on the proppant transport.Smaller proppants can be injected early,followed by larger proppants to maintain high propping efficiency.This study has enhanced the understanding of the complex flow phenomenon between proppant and fracturing fluid and can play a vital role in hydraulic fracturing design.
基金Supported by the Innovative Research Groups of National Natural Science Foundation of China(No.51321065)the Major Science and Technology Program for Water Pollution Control and Treatment(2012ZX07101-008)the National Natural Science Foundation of China(No.51439005)
文摘In the present study, considering the transport and transformation processes of variables, a threedimensional water quality model for the river system was established, which coupled the volume of fluid(VOF) method with the k-ε turbulence mathematical model. Then, the water hydrodynamic characteristics and transport processes for BOD_5, NH_(3^-)N and TP were analyzed. The results showed that the water surface of convex bank was a little lower than that of concave bank due to the centrifugal force near the bend, and most concentrations were inferior to the type Ⅴ standard indexes of surface water environmental quality. The model validation indicated that the errors between the simulated and monitored values were comparatively small, satisfying the application demands and providing scientific basis and decision support for the restoration and protection of water quality.
基金Sponsored by the National Nature Science Foundation of China(Grant No.51276030,51176017)
文摘In this paper, the lattice-Bohzmann method is used to investigate the droplet dynamics after impact on horizontal and inclined solid surface. The two-phase interparticle potential model is employed. The model is found to possess a linear relation between the macroscopic properties （ surface tension σ and contact angle α） and microscopic parameters （ G, G, ）. The flow state of the droplet on the surface is analyzed in detail, and the effects of surface characteristic, impact velocity, impact angle, the viscosity and surface tension of the liquid are investigated, respectively. It is shown that the lattice-Bohzmann method can not only track exactly and automatically the interface, but also the simulation results have a good qualitative agreement with ones of the previous experimental and numerical studies.
基金supported by Versus Arthritis UK(Grant no:21977)European Commission via a H2020-MSCA-RISE programme(BAMOS,Grant no:734156)+1 种基金Innovative UK via Newton Fund(Grant no:102872)Engineering and Physical Science Research Council(EPSRC)via DTP CASE programme(Grant no:EP/T517793/1).
文摘Nutrients supply especially like nutrients and oxygen play vital role in tissue engineering process.It is found that tissue could not grow very well in the middle of the scaffold because few nutrients could transport to the middle.Nutrient limitations would reduce cell proliferation and differentiation.In that case,there is urgent need to understand the nutrient distribution for both in vitro and in vivo study,as no technology is able for researchers to observe the nutrients transport during those process.In this paper,a numerical model coupling with VOF(volume of fluid)model and species transport model together for predicting the distribution of oxygen and glucose in the scaffold after implantation in to the site is developed.Comparing with our previous in vivo tests,the regenerated tissue distribution has a similar trend as oxygen distribution rather than glucose.The reported scaffold manufactured by additive manufacturing provided a good interconnected structure which facilitated the nutrient transportation in the scaffold.Considering nutrient transportation,this numerical model could be used in better understanding the nutrients transportation in the scaffold,and leading to a better understanding of tissue formation in the scaffold.
基金funded by the Science and Technology Innovation Committee Foundation of Shenzhen,China(Nos.JCYJ20200109141403840 and ZDSYS20220527171405012)the National Natural Science Foundation of China(No.52106045)the Pearl River Talent Recruitment Program,China(No.2019CX01Z084)。
文摘Iterative coupled methods are widely used in multi-fidelity simulation of rotating components due to the simple implementation,which iteratively eliminates the errors between the computational fluid dynamics models and approximate characteristic maps.However,the convergence and accuracy of the iterative coupled method are trapped in characteristic maps.In particular,iterative steps increase sharply as the operation point moves away from the design point.To address these problems,this paper developed an auxiliary iterative coupled method that introduces the static-pressure-auxiliary characteristic maps and modification factor of mass flow into the component-level model.The developed auxiliary method realized the direct transfer of static pressure between the high-fidelity models and the component-level model.Multi-fidelity simulations of the throttle characteristics were carried out using both the auxiliary and traditional iterative coupled methods,and the simulation results were verified using the experimental data.Additionally,the consistency between the auxiliary and traditional iterative coupled methods was confirmed.Subsequently,multi-fidelity simulations of the speed and altitude characteristics were also conducted.The auxiliary and traditional iterative coupled methods were evaluated in terms of convergence speed and accuracy.The evaluation indicated that the auxiliary iterative coupled method significantly reduces iterative steps by approximately 50%at the near-choked state.In general,the auxiliary iterative coupled method is preferred as a development of the traditional iterative coupled method in the near-choked state,and the combined auxiliary-traditional iterative coupled method provides support for successful multi-fidelity simulation in far-off-design conditions.
文摘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.
