Two-phase anaerobic digestion process is influenced by acid control for hydrogen production, reaction temperature, substrate detention time, sludge activity, and granular formation. Al of these technological parameter...Two-phase anaerobic digestion process is influenced by acid control for hydrogen production, reaction temperature, substrate detention time, sludge activity, and granular formation. Al of these technological parameters are directly related to success or failure of the system operation and treatment effect.展开更多
In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering comp...In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering complex boundary shapes.Utilizing radial basis function point interpolation,the method approximates shape functions for unknown functions within the nodal influence domain.The shape functions constructed by the aforementioned meshless interpolation method haveδ-function properties,which facilitate the handling of essential aspects like the controlled bottom-hole flow pressure in horizontal wells.Moreover,the meshless method offers greater flexibility and freedom compared to grid cell discretization,making it simpler to discretize complex geometries.A variational principle for the flow control equation group is introduced using a weighted least squares meshless method,and the pressure distribution is solved implicitly.Example results demonstrate that the computational outcomes of the meshless point cloud model,which has a relatively small degree of freedom,are in close agreement with those of the Discrete Fracture Model(DFM)employing refined grid partitioning,with pressure calculation accuracy exceeding 98.2%.Compared to high-resolution grid-based computational methods,the meshless method can achieve a better balance between computational efficiency and accuracy.Additionally,the impact of fracture half-length on the productivity of horizontal wells is discussed.The results indicate that increasing the fracture half-length is an effective strategy for enhancing production from the perspective of cumulative oil production.展开更多
A new anaerobic reactor, Jet-loop anaerobic fluidized bed (JLAFB), was designed for treating high-sulfate wastewater. The treatment characteristics, including the effect of influent COD/SO42 ratio and alkalinity and...A new anaerobic reactor, Jet-loop anaerobic fluidized bed (JLAFB), was designed for treating high-sulfate wastewater. The treatment characteristics, including the effect of influent COD/SO42 ratio and alkalinity and sulfide inhibition in reactors, were discussed for a JLAFB and a general anaerobic fiuidized bed (AFB) reactor used as sulfate-reducing phase and methane-producing phase, respectively, in two-phase anaerobic digestion process. The formation of granules in the two reactors was also examined. The results indicated that COD and sulfate removal had different demand of influent COD/SO4^2- ratios. When total COD removal was up to 85%, the ratio was only required up to 1.2, whereas, total sulfate removal up to 95% required it exceeding 3.0. The alkalinity in the two reactors increased linearly with the growth of influent alkalinity. Moreover, the change of influent alkalinity had no significant effect on pH and volatile fatty acids (VFA) in the two reactors. Influent alkalinity kept at 400-500 mg/L could meet the requirement of the treating process. The JLAFB reactor had great advantage in avoiding sulfide and free-H2S accumulation and toxicity inhibition on microorganisms. When sulfate loading rate was up to 8. 1 kg/(m^3.d), the sulfide and free-H2S concentrations in JLAFB reactor were 58.6 and 49.7 mg/L, respectively. Furthermore, the granules, with offwhite color, ellipse shape and diameters of 1.0-3.0 mm, could be developed in JLAFB reactor. In granules, different groups of bacteria were distributed in different layers, and some inorganic metal compounds such as Fe, Ca, Mg etc. were found.展开更多
The two-stage and two-phase anaerobic process (TSTP) composed of hydrolytic acidification reactor,first-order and second-order external circulation anaerobic reactors (EC) was taken to treat methanol wastewater. Test ...The two-stage and two-phase anaerobic process (TSTP) composed of hydrolytic acidification reactor,first-order and second-order external circulation anaerobic reactors (EC) was taken to treat methanol wastewater. Test results show that TSTP process is quick start-up in 51 d, and the maximum VFA of hydrolytic acidification reactor effluent reaches 876 mg/L. Under the condition of volume loading of 6.56 kgCOD/m3·d, COD removal rate of the first-order EC reactor is about 85%, and under the condition of volume loading of 1.02 kgCOD/m3·d, COD removal rate of the second-order EC reactor is about 50%. When the inflow COD of TSTP process is between 7000-11000 mg/L, its effluent COD is lower than 600 mg/L. In the biological conversion process of methanol into methane,the production of acetic acids as an intermediate product can be ignored and the direct production of methane from methanol is predominant.展开更多
The method of two-dimensional viscous space-time conservation element and solution element (CE/SE) can be used to calculate the gas-liquid two-phase interior flow field in pulse detonation engine (PDE). In this paper,...The method of two-dimensional viscous space-time conservation element and solution element (CE/SE) can be used to calculate the gas-liquid two-phase interior flow field in pulse detonation engine (PDE). In this paper, the evolution of the detonation wave and the distribution of its physical parameters were analyzed. The numerical results show that the change of axial velocity of gas is the same as that of detonation pressure. The larger the liquid droplet radius is, the longer the time to get stable detonation wave is. The calculated results coincide with the experimented results better.展开更多
Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on e...Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on element faces.Discontinuity of velocity field leads this method not to conserve mass locally.Moreover,the accuracy and stability of a solution is highly affected by a non-conservative method.In this paper,a three dimensional control volume finite element method is developed for twophase fluid flow simulation which overcomes the deficiency of the standard finite element method,and attains high-orders of accuracy at a reasonable computational cost.Moreover,this method is capable of handling heterogeneity in a very rational way.A fully implicit scheme is applied to temporal discretization of the governing equations to achieve an unconditionally stable solution.The accuracy and efficiency of the method are verified by simulating some waterflooding experiments.Some representative examples are presented to illustrate the capability of the method to simulate two-phase fluid flow in heterogeneous porous media.展开更多
Electrical capacitance tomography(ECT)has been applied to two-phase flow measurement in recent years.Image reconstruction algorithms play an important role in the successful applications of ECT.To solve the ill-posed ...Electrical capacitance tomography(ECT)has been applied to two-phase flow measurement in recent years.Image reconstruction algorithms play an important role in the successful applications of ECT.To solve the ill-posed and nonlinear inverse problem of ECT image reconstruction,a new ECT image reconstruction method based on fast linearized alternating direction method of multipliers(FLADMM)is proposed in this paper.On the basis of theoretical analysis of compressed sensing(CS),the data acquisition of ECT is regarded as a linear measurement process of permittivity distribution signal of pipe section.A new measurement matrix is designed and L1 regularization method is used to convert ECT inverse problem to a convex relaxation problem which contains prior knowledge.A new fast alternating direction method of multipliers which contained linearized idea is employed to minimize the objective function.Simulation data and experimental results indicate that compared with other methods,the quality and speed of reconstructed images are markedly improved.Also,the dynamic experimental results indicate that the proposed algorithm can ful fill the real-time requirement of ECT systems in the application.展开更多
Due to the coupling between the hydrodynamic equation and the phase-field equation in two-phase incompressible flows,it is desirable to develop efficient and high-order accurate numerical schemes that can decouple the...Due to the coupling between the hydrodynamic equation and the phase-field equation in two-phase incompressible flows,it is desirable to develop efficient and high-order accurate numerical schemes that can decouple these two equations.One popular and efficient strategy is to add an explicit stabilizing term to the convective velocity in the phase-field equation to decouple them.The resulting schemes are only first-order accurate in time,and it seems extremely difficult to generalize the idea of stabilization to the second-order or higher version.In this paper,we employ the spectral deferred correction method to improve the temporal accuracy,based on the first-order decoupled and energy-stable scheme constructed by the stabilization idea.The novelty lies in how the decoupling and linear implicit properties are maintained to improve the efficiency.Within the framework of the spatially discretized local discontinuous Galerkin method,the resulting numerical schemes are fully decoupled,efficient,and high-order accurate in both time and space.Numerical experiments are performed to validate the high-order accuracy and efficiency of the methods for solving phase-field models of two-phase incompressible flows.展开更多
The relationship between capillary pressure and saturation plays a critical role in the characterization of two-phase flow and transport in aquifers and oil reservoirs. This relationship is usually determined under th...The relationship between capillary pressure and saturation plays a critical role in the characterization of two-phase flow and transport in aquifers and oil reservoirs. This relationship is usually determined under the static condition, where capillary pressure is the only function of saturation. However,considerable experiments have suggested that the dependence of capillary pressure on desaturation rate is under the dynamic condition. Thus, a more general description of capillary pressure that includes dynamic capillary effect has been approved widely. A comprehensive understanding of the dynamic capillary effect is needed for the investigation of the two-phase flow in porous media by various methods. In general, dynamic capillary effect in porous media can be studied through the laboratory experiment, pore-to macro-scale modeling, and artificial neural network. Here, main principle and research procedures of each method are reviewed in detail. Then, research progress, disadvantages and advantages are discussed, respectively. In addition, upscaling study from pore-to macro-scale are introduced, which explains the difference between laboratory experiment and pore-scale modeling. At last, several future perspectives and recommendations for optimal solution of dynamic capillary effect are presented.展开更多
The knowledge of bubble profiles in gas-liquid two-phase flows is crucial for analyzing the kinetic processes such as heat and mass transfer, and this knowledge is contained in field data obtained by surface-resolved ...The knowledge of bubble profiles in gas-liquid two-phase flows is crucial for analyzing the kinetic processes such as heat and mass transfer, and this knowledge is contained in field data obtained by surface-resolved computational fluid dynamics (CFD) simulations. To obtain this information, an efficient bubble profile reconstruction method based on an improved agglomerative hierarchical clustering (AHC) algorithm is proposed in this paper. The reconstruction method is featured by the implementations of a binary space division preprocessing, which aims to reduce the computational complexity, an adaptive linkage criterion, which guarantees the applicability of the AHC algorithm when dealing with datasets involving either non-uniform or distorted grids, and a stepwise execution strategy, which enables the separation of attached bubbles. To illustrate and verify this method, it was applied to dealing with 3 datasets, 2 of them with pre-specified spherical bubbles and the other obtained by a surface-resolved CFD simulation. Application results indicate that the proposed method is effective even when the data include some non-uniform and distortion.展开更多
By using the interaction of particles, such as the physical principle of the same attract each other and the different repulse each other, a new model of Lattice Boltzmann to simulate the two-phase driven in porous me...By using the interaction of particles, such as the physical principle of the same attract each other and the different repulse each other, a new model of Lattice Boltzmann to simulate the two-phase driven in porous media was discussed. The result shows effectively for the problem of two-phase driven in porous media. Furthermore, the method economizes on computer time, has less fluctuation on boundary surface and takes no average measure.展开更多
A kinetic study of biogas production from Urban Solid Waste (USW) generated in Dar es Salaam city (Tanzania) is presented. An experimental bioreactor simulating mesophilic conditions of most USW landfills was develope...A kinetic study of biogas production from Urban Solid Waste (USW) generated in Dar es Salaam city (Tanzania) is presented. An experimental bioreactor simulating mesophilic conditions of most USW landfills was developed. The goal of the study was to generate the kinetic order of reaction with respect to biodegradable organic waste and use it to model biogas production from food residues mixed with fruit waste. Anaerobic biodegradation was employed under temperature range of 28℃ - 38℃. The main controls were leachate recirculation and pH adjustments to minimize acid inhibitory effects and accelerate waste biodegradation. The experimental setup comprised of three sets of bioreactors. A biodegradation rate law in differential form was proposed and the numerical values of kinetic order and rate constant were determined using initial rate method as 0.994 and 0.3093 mol0.006·day-1, respectively. Results obtained were consistent with that found in literature and model predictions were in reasonable agreement with experimental data.展开更多
A method is provided to achieve an initial basic feasible solution of a linear programming in this paper. This method dose not need introducing any artificial variable, but needs only solving an auxiliary linear progr...A method is provided to achieve an initial basic feasible solution of a linear programming in this paper. This method dose not need introducing any artificial variable, but needs only solving an auxiliary linear programming. Compared with the traditional two-phase method, it has advantages of saving the memories and reducing the computational efforts.展开更多
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.展开更多
In this paper we shall give the characteristic difference methods for two phase displace meat problem in naturally fractured reservoirs.We shall prove the existence,uniqueness of the ap proximate solution and a priori...In this paper we shall give the characteristic difference methods for two phase displace meat problem in naturally fractured reservoirs.We shall prove the existence,uniqueness of the ap proximate solution and a priori discrete L2-error estimates.展开更多
We developed a three-dimensional multi-relaxation-time lattice Boltzmann method for incompressible and immiscible two-phase flow by coupling with a front-tracking technique. The flow field was simulated by using an Eu...We developed a three-dimensional multi-relaxation-time lattice Boltzmann method for incompressible and immiscible two-phase flow by coupling with a front-tracking technique. The flow field was simulated by using an Eulerian grid, an adaptive unstructured triangular Lagrangian grid was applied to track explicitly the motion of the two-fluid interface, and an indicator function was introduced to update accurately the fluid properties. The surface tension was computed directly on a triangular Lagrangian grid, and then the surface tension was distributed to the background Eulerian grid. Three benchmarks of two-phase flow, including the Laplace law for a stationary drop, the oscillation of a three-dimensional ellipsoidal drop, and the drop deformation in a shear flow, were simulated to validate the present model.展开更多
In this paper, an improved incompressible multi-relaxation-time lattice Boltzmann-front tracking approach is proposed to simulate two-phase flow with a sharp interface, where the surface tension is implemented. The la...In this paper, an improved incompressible multi-relaxation-time lattice Boltzmann-front tracking approach is proposed to simulate two-phase flow with a sharp interface, where the surface tension is implemented. The lattice Boltzmann method is used to simulate the incompressible flow with a stationary Eulerian grid, an additional moving Lagrangian grid is adopted to track explicitly the motion of the interface, and an indicator function is introduced to update the fluid properties accurately. The interface is represented by using a four-order Lagrange polynomial through fitting a set of discrete marker points, and then the surface tension is directly computed by using the normal vector and curvature of the interface. Two benchmark problems, including Laplace's law for a stationary bubble and the dispersion relation of the capillary wave between two fluids are conducted for validation. Excellent agreement is obtained between the numerical simulations and the theoretical results in the two cases.展开更多
Motivated by inconveniences of present hybrid methods,a gradient-augmented hybrid interface capturing method(GAHM) is presented for incompressible two-phase flow.A front tracking method(FTM) is used as the skeleto...Motivated by inconveniences of present hybrid methods,a gradient-augmented hybrid interface capturing method(GAHM) is presented for incompressible two-phase flow.A front tracking method(FTM) is used as the skeleton of the GAHM for low mass loss and resources.Smooth eulerian level set values are calculated from the FTM interface,and are used for a local interface reconstruction.The reconstruction avoids marker particle redistribution and enables an automatic treatment of interfacial topology change.The cubic Hermit interpolation is employed in all steps of the GAHM to capture subgrid structures within a single spacial cell.The performance of the GAHM is carefully evaluated in a benchmark test.Results show significant improvements of mass loss,clear subgrid structures,highly accurate derivatives(normals and curvatures) and low cost.The GAHM is further coupled with an incompressible multiphase flow solver,Super CE/SE,for more complex and practical applications.The updated solver is evaluated through comparison with an early droplet research.展开更多
A lattice Boltzmann method for gas–liquid two-phase flow involving non-Newtonian fluids is developed. Bubble formation in a flow-focusing microchannel is simulated by the method. The influences of flow rate ratio, su...A lattice Boltzmann method for gas–liquid two-phase flow involving non-Newtonian fluids is developed. Bubble formation in a flow-focusing microchannel is simulated by the method. The influences of flow rate ratio, surface tension,wetting properties, and rheological characteristics of the fluid on the two-phase flow are analyzed. The results indicate that the flow pattern transfers from slug flow to dry-plug flow with a sufficiently small capillary number. Due to the presence of three-phase contact lines, the contact angle has a more significant effect on the dry-plug flow pattern than on the slug flow pattern. The deformation of the front and rear meniscus of a bubble in the shear-thinning fluid can be explained by the variation of the capillary number. The reduced viscosity and increased contact angle are beneficial for the drag reduction in a microchannel. It also demonstrates the effectiveness of the current method to simulate the gas–liquid two-phase flow in a microchannel.展开更多
In this paper the transient two-phase flow equations and their eigenvalues are first introduced. The flux vector is then split into subvectors which just contain a specially signed eigenvalue. Using one-sided spatial ...In this paper the transient two-phase flow equations and their eigenvalues are first introduced. The flux vector is then split into subvectors which just contain a specially signed eigenvalue. Using one-sided spatial difference operators finite difference equations and their solutions are obtained. Finally comparison with experiment shows the predicted results produce good agreement with experimental data.展开更多
基金Supported by Yunnan S&T Innovation Platform Construction Project(2013DH041)National Natural Science Foundation of China(51366015)Specialized Research Fund for the Doctoral Program of Higher Education of China(20135303110001)~~
文摘Two-phase anaerobic digestion process is influenced by acid control for hydrogen production, reaction temperature, substrate detention time, sludge activity, and granular formation. Al of these technological parameters are directly related to success or failure of the system operation and treatment effect.
文摘In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering complex boundary shapes.Utilizing radial basis function point interpolation,the method approximates shape functions for unknown functions within the nodal influence domain.The shape functions constructed by the aforementioned meshless interpolation method haveδ-function properties,which facilitate the handling of essential aspects like the controlled bottom-hole flow pressure in horizontal wells.Moreover,the meshless method offers greater flexibility and freedom compared to grid cell discretization,making it simpler to discretize complex geometries.A variational principle for the flow control equation group is introduced using a weighted least squares meshless method,and the pressure distribution is solved implicitly.Example results demonstrate that the computational outcomes of the meshless point cloud model,which has a relatively small degree of freedom,are in close agreement with those of the Discrete Fracture Model(DFM)employing refined grid partitioning,with pressure calculation accuracy exceeding 98.2%.Compared to high-resolution grid-based computational methods,the meshless method can achieve a better balance between computational efficiency and accuracy.Additionally,the impact of fracture half-length on the productivity of horizontal wells is discussed.The results indicate that increasing the fracture half-length is an effective strategy for enhancing production from the perspective of cumulative oil production.
基金Project supported by the National Natural Science Foundation of China(No. 50278036)the Natural Science Foundation of Guangdong Province (No. 04105951)
文摘A new anaerobic reactor, Jet-loop anaerobic fluidized bed (JLAFB), was designed for treating high-sulfate wastewater. The treatment characteristics, including the effect of influent COD/SO42 ratio and alkalinity and sulfide inhibition in reactors, were discussed for a JLAFB and a general anaerobic fiuidized bed (AFB) reactor used as sulfate-reducing phase and methane-producing phase, respectively, in two-phase anaerobic digestion process. The formation of granules in the two reactors was also examined. The results indicated that COD and sulfate removal had different demand of influent COD/SO4^2- ratios. When total COD removal was up to 85%, the ratio was only required up to 1.2, whereas, total sulfate removal up to 95% required it exceeding 3.0. The alkalinity in the two reactors increased linearly with the growth of influent alkalinity. Moreover, the change of influent alkalinity had no significant effect on pH and volatile fatty acids (VFA) in the two reactors. Influent alkalinity kept at 400-500 mg/L could meet the requirement of the treating process. The JLAFB reactor had great advantage in avoiding sulfide and free-H2S accumulation and toxicity inhibition on microorganisms. When sulfate loading rate was up to 8. 1 kg/(m^3.d), the sulfide and free-H2S concentrations in JLAFB reactor were 58.6 and 49.7 mg/L, respectively. Furthermore, the granules, with offwhite color, ellipse shape and diameters of 1.0-3.0 mm, could be developed in JLAFB reactor. In granules, different groups of bacteria were distributed in different layers, and some inorganic metal compounds such as Fe, Ca, Mg etc. were found.
基金Sponsored by the National Hi-Tech Research and Development Program of China (Grant No.2003AA601090)Projects of Development Plan of the State Key Fundamental Research of China (Grant No.2004CB4185)
文摘The two-stage and two-phase anaerobic process (TSTP) composed of hydrolytic acidification reactor,first-order and second-order external circulation anaerobic reactors (EC) was taken to treat methanol wastewater. Test results show that TSTP process is quick start-up in 51 d, and the maximum VFA of hydrolytic acidification reactor effluent reaches 876 mg/L. Under the condition of volume loading of 6.56 kgCOD/m3·d, COD removal rate of the first-order EC reactor is about 85%, and under the condition of volume loading of 1.02 kgCOD/m3·d, COD removal rate of the second-order EC reactor is about 50%. When the inflow COD of TSTP process is between 7000-11000 mg/L, its effluent COD is lower than 600 mg/L. In the biological conversion process of methanol into methane,the production of acetic acids as an intermediate product can be ignored and the direct production of methane from methanol is predominant.
基金Sponsored by the National Natural Science Foundation of China (10672080)
文摘The method of two-dimensional viscous space-time conservation element and solution element (CE/SE) can be used to calculate the gas-liquid two-phase interior flow field in pulse detonation engine (PDE). In this paper, the evolution of the detonation wave and the distribution of its physical parameters were analyzed. The numerical results show that the change of axial velocity of gas is the same as that of detonation pressure. The larger the liquid droplet radius is, the longer the time to get stable detonation wave is. The calculated results coincide with the experimented results better.
基金Iranian Offshore Oil Company (IOOC) for financial support of this work
文摘Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on element faces.Discontinuity of velocity field leads this method not to conserve mass locally.Moreover,the accuracy and stability of a solution is highly affected by a non-conservative method.In this paper,a three dimensional control volume finite element method is developed for twophase fluid flow simulation which overcomes the deficiency of the standard finite element method,and attains high-orders of accuracy at a reasonable computational cost.Moreover,this method is capable of handling heterogeneity in a very rational way.A fully implicit scheme is applied to temporal discretization of the governing equations to achieve an unconditionally stable solution.The accuracy and efficiency of the method are verified by simulating some waterflooding experiments.Some representative examples are presented to illustrate the capability of the method to simulate two-phase fluid flow in heterogeneous porous media.
基金Supported by the National Natural Science Foundation of China(61203021)the Key Science and Technology Program of Liaoning Province(2011216011)+1 种基金the Natural Science Foundation of Liaoning Province(2013020024)the Program for Liaoning Excellent Talents in Universities(LJQ2015061)
文摘Electrical capacitance tomography(ECT)has been applied to two-phase flow measurement in recent years.Image reconstruction algorithms play an important role in the successful applications of ECT.To solve the ill-posed and nonlinear inverse problem of ECT image reconstruction,a new ECT image reconstruction method based on fast linearized alternating direction method of multipliers(FLADMM)is proposed in this paper.On the basis of theoretical analysis of compressed sensing(CS),the data acquisition of ECT is regarded as a linear measurement process of permittivity distribution signal of pipe section.A new measurement matrix is designed and L1 regularization method is used to convert ECT inverse problem to a convex relaxation problem which contains prior knowledge.A new fast alternating direction method of multipliers which contained linearized idea is employed to minimize the objective function.Simulation data and experimental results indicate that compared with other methods,the quality and speed of reconstructed images are markedly improved.Also,the dynamic experimental results indicate that the proposed algorithm can ful fill the real-time requirement of ECT systems in the application.
基金supported by the NSFC Grant no.12271492the Natural Science Foundation of Henan Province of China Grant no.222300420550+1 种基金supported by the NSFC Grant no.12271498the National Key R&D Program of China Grant no.2022YFA1005202/2022YFA1005200.
文摘Due to the coupling between the hydrodynamic equation and the phase-field equation in two-phase incompressible flows,it is desirable to develop efficient and high-order accurate numerical schemes that can decouple these two equations.One popular and efficient strategy is to add an explicit stabilizing term to the convective velocity in the phase-field equation to decouple them.The resulting schemes are only first-order accurate in time,and it seems extremely difficult to generalize the idea of stabilization to the second-order or higher version.In this paper,we employ the spectral deferred correction method to improve the temporal accuracy,based on the first-order decoupled and energy-stable scheme constructed by the stabilization idea.The novelty lies in how the decoupling and linear implicit properties are maintained to improve the efficiency.Within the framework of the spatially discretized local discontinuous Galerkin method,the resulting numerical schemes are fully decoupled,efficient,and high-order accurate in both time and space.Numerical experiments are performed to validate the high-order accuracy and efficiency of the methods for solving phase-field models of two-phase incompressible flows.
基金financially supported by the National Natural Science Foundation of China (No. 42102149)the Fundamental Research Funds for the Central Universities (No. 2462021YXZZ005)。
文摘The relationship between capillary pressure and saturation plays a critical role in the characterization of two-phase flow and transport in aquifers and oil reservoirs. This relationship is usually determined under the static condition, where capillary pressure is the only function of saturation. However,considerable experiments have suggested that the dependence of capillary pressure on desaturation rate is under the dynamic condition. Thus, a more general description of capillary pressure that includes dynamic capillary effect has been approved widely. A comprehensive understanding of the dynamic capillary effect is needed for the investigation of the two-phase flow in porous media by various methods. In general, dynamic capillary effect in porous media can be studied through the laboratory experiment, pore-to macro-scale modeling, and artificial neural network. Here, main principle and research procedures of each method are reviewed in detail. Then, research progress, disadvantages and advantages are discussed, respectively. In addition, upscaling study from pore-to macro-scale are introduced, which explains the difference between laboratory experiment and pore-scale modeling. At last, several future perspectives and recommendations for optimal solution of dynamic capillary effect are presented.
基金Projects(51634010,51676211) supported by the National Natural Science Foundation of ChinaProject(2017SK2253) supported by the Key Research and Development Program of Hunan Province,China
文摘The knowledge of bubble profiles in gas-liquid two-phase flows is crucial for analyzing the kinetic processes such as heat and mass transfer, and this knowledge is contained in field data obtained by surface-resolved computational fluid dynamics (CFD) simulations. To obtain this information, an efficient bubble profile reconstruction method based on an improved agglomerative hierarchical clustering (AHC) algorithm is proposed in this paper. The reconstruction method is featured by the implementations of a binary space division preprocessing, which aims to reduce the computational complexity, an adaptive linkage criterion, which guarantees the applicability of the AHC algorithm when dealing with datasets involving either non-uniform or distorted grids, and a stepwise execution strategy, which enables the separation of attached bubbles. To illustrate and verify this method, it was applied to dealing with 3 datasets, 2 of them with pre-specified spherical bubbles and the other obtained by a surface-resolved CFD simulation. Application results indicate that the proposed method is effective even when the data include some non-uniform and distortion.
文摘By using the interaction of particles, such as the physical principle of the same attract each other and the different repulse each other, a new model of Lattice Boltzmann to simulate the two-phase driven in porous media was discussed. The result shows effectively for the problem of two-phase driven in porous media. Furthermore, the method economizes on computer time, has less fluctuation on boundary surface and takes no average measure.
文摘A kinetic study of biogas production from Urban Solid Waste (USW) generated in Dar es Salaam city (Tanzania) is presented. An experimental bioreactor simulating mesophilic conditions of most USW landfills was developed. The goal of the study was to generate the kinetic order of reaction with respect to biodegradable organic waste and use it to model biogas production from food residues mixed with fruit waste. Anaerobic biodegradation was employed under temperature range of 28℃ - 38℃. The main controls were leachate recirculation and pH adjustments to minimize acid inhibitory effects and accelerate waste biodegradation. The experimental setup comprised of three sets of bioreactors. A biodegradation rate law in differential form was proposed and the numerical values of kinetic order and rate constant were determined using initial rate method as 0.994 and 0.3093 mol0.006·day-1, respectively. Results obtained were consistent with that found in literature and model predictions were in reasonable agreement with experimental data.
文摘A method is provided to achieve an initial basic feasible solution of a linear programming in this paper. This method dose not need introducing any artificial variable, but needs only solving an auxiliary linear programming. Compared with the traditional two-phase method, it has advantages of saving the memories and reducing the computational efforts.
文摘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 National Natural Science Foundation of ChinaChina Key Project for Basic Researches
文摘In this paper we shall give the characteristic difference methods for two phase displace meat problem in naturally fractured reservoirs.We shall prove the existence,uniqueness of the ap proximate solution and a priori discrete L2-error estimates.
基金supported by the National Natural Science Foundation of China(Grant No.11572062)the Fundamental Research Funds for the Central Universities,China(Grant No.CDJZR13248801)+2 种基金the Program for Changjiang Scholars and Innovative Research Team in University,China(Grant No.IRT13043)Key Laboratory of Functional Crystals and Laser Technology,TIPCChinese Academy of Sciences
文摘We developed a three-dimensional multi-relaxation-time lattice Boltzmann method for incompressible and immiscible two-phase flow by coupling with a front-tracking technique. The flow field was simulated by using an Eulerian grid, an adaptive unstructured triangular Lagrangian grid was applied to track explicitly the motion of the two-fluid interface, and an indicator function was introduced to update accurately the fluid properties. The surface tension was computed directly on a triangular Lagrangian grid, and then the surface tension was distributed to the background Eulerian grid. Three benchmarks of two-phase flow, including the Laplace law for a stationary drop, the oscillation of a three-dimensional ellipsoidal drop, and the drop deformation in a shear flow, were simulated to validate the present model.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.10872222 and 50921063)the Specialized Research Fund for the Doctoral Program of Higher Education of China(Grant No.20110191110037)the Fundamental Research Funds for the Central Universities,China(Grant Nos.CDJXS11240011 and CDJXS10241103)
文摘In this paper, an improved incompressible multi-relaxation-time lattice Boltzmann-front tracking approach is proposed to simulate two-phase flow with a sharp interface, where the surface tension is implemented. The lattice Boltzmann method is used to simulate the incompressible flow with a stationary Eulerian grid, an additional moving Lagrangian grid is adopted to track explicitly the motion of the interface, and an indicator function is introduced to update the fluid properties accurately. The interface is represented by using a four-order Lagrange polynomial through fitting a set of discrete marker points, and then the surface tension is directly computed by using the normal vector and curvature of the interface. Two benchmark problems, including Laplace's law for a stationary bubble and the dispersion relation of the capillary wave between two fluids are conducted for validation. Excellent agreement is obtained between the numerical simulations and the theoretical results in the two cases.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.10972010,11028206,11371069,11372052,11402029,and 11472060)the Science and Technology Development Foundation of China Academy of Engineering Physics(CAEP),China(Grant No.2014B0201030)the Defense Industrial Technology Development Program of China(Grant No.B1520132012)
文摘Motivated by inconveniences of present hybrid methods,a gradient-augmented hybrid interface capturing method(GAHM) is presented for incompressible two-phase flow.A front tracking method(FTM) is used as the skeleton of the GAHM for low mass loss and resources.Smooth eulerian level set values are calculated from the FTM interface,and are used for a local interface reconstruction.The reconstruction avoids marker particle redistribution and enables an automatic treatment of interfacial topology change.The cubic Hermit interpolation is employed in all steps of the GAHM to capture subgrid structures within a single spacial cell.The performance of the GAHM is carefully evaluated in a benchmark test.Results show significant improvements of mass loss,clear subgrid structures,highly accurate derivatives(normals and curvatures) and low cost.The GAHM is further coupled with an incompressible multiphase flow solver,Super CE/SE,for more complex and practical applications.The updated solver is evaluated through comparison with an early droplet research.
基金Project supported by the National Natural Science Foundation of China (Grant No. 51775077)。
文摘A lattice Boltzmann method for gas–liquid two-phase flow involving non-Newtonian fluids is developed. Bubble formation in a flow-focusing microchannel is simulated by the method. The influences of flow rate ratio, surface tension,wetting properties, and rheological characteristics of the fluid on the two-phase flow are analyzed. The results indicate that the flow pattern transfers from slug flow to dry-plug flow with a sufficiently small capillary number. Due to the presence of three-phase contact lines, the contact angle has a more significant effect on the dry-plug flow pattern than on the slug flow pattern. The deformation of the front and rear meniscus of a bubble in the shear-thinning fluid can be explained by the variation of the capillary number. The reduced viscosity and increased contact angle are beneficial for the drag reduction in a microchannel. It also demonstrates the effectiveness of the current method to simulate the gas–liquid two-phase flow in a microchannel.
文摘In this paper the transient two-phase flow equations and their eigenvalues are first introduced. The flux vector is then split into subvectors which just contain a specially signed eigenvalue. Using one-sided spatial difference operators finite difference equations and their solutions are obtained. Finally comparison with experiment shows the predicted results produce good agreement with experimental data.