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.展开更多
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.展开更多
A cold-model vertical multi-tube circulating fluidized bed evaporator was designed and built to conduct a visualization study on the pressure drop of a liquid–solid two-phase flow and the corresponding particle distr...A cold-model vertical multi-tube circulating fluidized bed evaporator was designed and built to conduct a visualization study on the pressure drop of a liquid–solid two-phase flow and the corresponding particle distribution.Water and polyformaldehyde particle(POM)were used as the liquid and solid phases,respectively.The effects of operating parameters such as the amount of added particles,circulating flow rate,and particle size were systematically investigated.The results showed that the addition of the particles increased the pressure drop in the vertical tube bundle.The maximum pressure drop ratios were 18.65%,21.15%,18.00%,and 21.15%within the experimental range of the amount of added particles for POM1,POM2,POM3,and POM4,respectively.The pressure drop ratio basically decreased with the increase in the circulating flow rate but fluctuated with the increase in the amount of added particles and particle size.The difference in pressure drop ratio decreased with the increase in the circulating flow rate.As the amount of added particles increased,the difference in pressure drop ratio fluctuated at low circulating flow rate but basically decreased at high circulating flow rate.The pressure drop in the vertical tube bundle accounted for about 70%of the overall pressure drop in the up-flow heating chamber and was the main component of the overall pressure within the experimental range.Three-dimensional phase diagrams were established to display the variation ranges of the pressure drop and pressure drop ratio in the vertical tube bundle corresponding to the operating parameters.The research results can provide some reference for the application of the fluidized bed heat transfer technology in the industry.展开更多
In order to optimize the design of the submerged combustion vaporizer(SCV), an experimental apparatus was set up to investigate the heat transfer character outside the tube bundle in SCV. Several experiments were cond...In order to optimize the design of the submerged combustion vaporizer(SCV), an experimental apparatus was set up to investigate the heat transfer character outside the tube bundle in SCV. Several experiments were conducted using water and CO_2 as the heat transfer media in the tubes, respectively. The results indicated that hot air flux, the initial liquid level height and the tube pitch ratio had great influence on the heat transfer coefficient outside the tube bundle(ho). Finally, the air flux associated factor β and height associated factor γ were introduced to propose a new hocorrelation. After verified by experiments using cold water, high pressure CO_2 and liquid N_2 as heat transfer media, respectively, it was found that the biggest deviation between the predicted and the experimental values was less than 25%.展开更多
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.展开更多
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.展开更多
A nonlocal study of the vibration responses of functionally graded(FG)beams supported by a viscoelastic Winkler-Pasternak foundation is presented.The damping responses of both the Winkler and Pasternak layers of the f...A nonlocal study of the vibration responses of functionally graded(FG)beams supported by a viscoelastic Winkler-Pasternak foundation is presented.The damping responses of both the Winkler and Pasternak layers of the foundation are considered in the formulation,which were not considered in most literature on this subject,and the bending deformation of the beams and the elastic and damping responses of the foundation as nonlocal by uniting the equivalently differential formulation of well-posed strain-driven(ε-D)and stress-driven(σ-D)two-phase local/nonlocal integral models with constitutive constraints are comprehensively considered,which can address both the stiffness softening and toughing effects due to scale reduction.The generalized differential quadrature method(GDQM)is used to solve the complex eigenvalue problem.After verifying the solution procedure,a series of benchmark results for the vibration frequency of different bounded FG beams supported by the foundation are obtained.Subsequently,the effects of the nonlocality of the foundation on the undamped/damping vibration frequency of the beams are examined.展开更多
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 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.展开更多
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.展开更多
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.展开更多
Air entrapped in liquid metal during the mold filling process seriously affects the casting quality, thus it is important to track its behavior in the mold cavity. A liquid-gas two-phase flow model is developed to des...Air entrapped in liquid metal during the mold filling process seriously affects the casting quality, thus it is important to track its behavior in the mold cavity. A liquid-gas two-phase flow model is developed to describe the mold filling process and predict the air entrapment defect. The model is based on the combination of SOLA and Level Set Method. The pressure and velocity fields are calculated by SOLA,and the interface movement is simulated by Level Set method as the most common interface tracking method in recent years.In order to validate the feasibility of the model,the liquid-gas two-phase simulation results were tested by the broken dam problem and the S-shaped experiment. Comparison between the experiments and simulation results show that Level Set method might be a very promising tool in two-phase flow simulation during the mold filling process.展开更多
Background The interaction of gas and liquid can produce many interesting phenomena,such as bubbles rising from the bottom of the liquid.The simulation of two-phase fluids is a challenging topic in computer graphics.T...Background The interaction of gas and liquid can produce many interesting phenomena,such as bubbles rising from the bottom of the liquid.The simulation of two-phase fluids is a challenging topic in computer graphics.To animate the interaction of a gas and liquid,MultiFLIP samples the two types of particles,and a Euler grid is used to track the interface of the liquid and gas.However,MultiFLIP uses the fluid implicit particle(FLIP)method to interpolate the velocities of particles into the Euler grid,which suffer from additional noise and instability.Methods To solve the problem caused by fluid implicit particles(FLIP),we present a novel velocity transport technique for two individual particles based on the affine particle-in-cell(APIC)method.First,we design a weighed coupling method for interpolating the velocities of liquid and gas particles to the Euler grid such that we can apply the APIC method to the simulation of a two-phase fluid.Second,we introduce a narrowband method to our system because MultiFLIP is a time-consuming approach owing to the large number of particles.Results Experiments show that our method is well integrated with the APIC method and provides a visually credible two-phase fluid animation.Conclusions The proposed method can successfully handle the simulation of a two phase fluid.展开更多
An explicit,time-dependent variable grid finite difference method is introduced and analyzed for approximating the solution of a scalar conservation law in two dimension. The scheme is stable,and the numerical solutio...An explicit,time-dependent variable grid finite difference method is introduced and analyzed for approximating the solution of a scalar conservation law in two dimension. The scheme is stable,and the numerical solution is proved to converges to the relevant physical solution.展开更多
基金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.
基金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.
基金supported by the open foundation of State Key Laboratory of Chemical Engineering (SKL-ChE-18B03)the Municipal Science and Technology Commission of Tianjin (No. 2009ZCKFGX01900)
文摘A cold-model vertical multi-tube circulating fluidized bed evaporator was designed and built to conduct a visualization study on the pressure drop of a liquid–solid two-phase flow and the corresponding particle distribution.Water and polyformaldehyde particle(POM)were used as the liquid and solid phases,respectively.The effects of operating parameters such as the amount of added particles,circulating flow rate,and particle size were systematically investigated.The results showed that the addition of the particles increased the pressure drop in the vertical tube bundle.The maximum pressure drop ratios were 18.65%,21.15%,18.00%,and 21.15%within the experimental range of the amount of added particles for POM1,POM2,POM3,and POM4,respectively.The pressure drop ratio basically decreased with the increase in the circulating flow rate but fluctuated with the increase in the amount of added particles and particle size.The difference in pressure drop ratio decreased with the increase in the circulating flow rate.As the amount of added particles increased,the difference in pressure drop ratio fluctuated at low circulating flow rate but basically decreased at high circulating flow rate.The pressure drop in the vertical tube bundle accounted for about 70%of the overall pressure drop in the up-flow heating chamber and was the main component of the overall pressure within the experimental range.Three-dimensional phase diagrams were established to display the variation ranges of the pressure drop and pressure drop ratio in the vertical tube bundle corresponding to the operating parameters.The research results can provide some reference for the application of the fluidized bed heat transfer technology in the industry.
文摘In order to optimize the design of the submerged combustion vaporizer(SCV), an experimental apparatus was set up to investigate the heat transfer character outside the tube bundle in SCV. Several experiments were conducted using water and CO_2 as the heat transfer media in the tubes, respectively. The results indicated that hot air flux, the initial liquid level height and the tube pitch ratio had great influence on the heat transfer coefficient outside the tube bundle(ho). Finally, the air flux associated factor β and height associated factor γ were introduced to propose a new hocorrelation. After verified by experiments using cold water, high pressure CO_2 and liquid N_2 as heat transfer media, respectively, it was found that the biggest deviation between the predicted and the experimental values was less than 25%.
基金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.
基金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.
基金the National Natural Science Foundation of China(No.12172169)the China Scholarship Council(CSC)(No.202006830038)the Natural Sciences and Engineering Research Council of Canada(No.RGPIN-2017-03716115112)。
文摘A nonlocal study of the vibration responses of functionally graded(FG)beams supported by a viscoelastic Winkler-Pasternak foundation is presented.The damping responses of both the Winkler and Pasternak layers of the foundation are considered in the formulation,which were not considered in most literature on this subject,and the bending deformation of the beams and the elastic and damping responses of the foundation as nonlocal by uniting the equivalently differential formulation of well-posed strain-driven(ε-D)and stress-driven(σ-D)two-phase local/nonlocal integral models with constitutive constraints are comprehensively considered,which can address both the stiffness softening and toughing effects due to scale reduction.The generalized differential quadrature method(GDQM)is used to solve the complex eigenvalue problem.After verifying the solution procedure,a series of benchmark results for the vibration frequency of different bounded FG beams supported by the foundation are obtained.Subsequently,the effects of the nonlocality of the foundation on the undamped/damping vibration frequency of the beams are examined.
文摘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 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.
基金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.
文摘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.
基金National High Technology Research and Development Program of China (863program) (2006AA04Z140)National Natural Science Foundation of China (NSFC) (50605024)
文摘Air entrapped in liquid metal during the mold filling process seriously affects the casting quality, thus it is important to track its behavior in the mold cavity. A liquid-gas two-phase flow model is developed to describe the mold filling process and predict the air entrapment defect. The model is based on the combination of SOLA and Level Set Method. The pressure and velocity fields are calculated by SOLA,and the interface movement is simulated by Level Set method as the most common interface tracking method in recent years.In order to validate the feasibility of the model,the liquid-gas two-phase simulation results were tested by the broken dam problem and the S-shaped experiment. Comparison between the experiments and simulation results show that Level Set method might be a very promising tool in two-phase flow simulation during the mold filling process.
基金National High Technology R&D Program of China(2017YFB1002701,M2019YFB1600702)NSFC(62072449)+1 种基金Science and Technology Development Fund,Macao SAR(0018/2019/AKP,0008/2019/AGJ,SKL-IOTSC-2018-2020)University of Macao Grant(MYRG2019-00006-FST).
文摘Background The interaction of gas and liquid can produce many interesting phenomena,such as bubbles rising from the bottom of the liquid.The simulation of two-phase fluids is a challenging topic in computer graphics.To animate the interaction of a gas and liquid,MultiFLIP samples the two types of particles,and a Euler grid is used to track the interface of the liquid and gas.However,MultiFLIP uses the fluid implicit particle(FLIP)method to interpolate the velocities of particles into the Euler grid,which suffer from additional noise and instability.Methods To solve the problem caused by fluid implicit particles(FLIP),we present a novel velocity transport technique for two individual particles based on the affine particle-in-cell(APIC)method.First,we design a weighed coupling method for interpolating the velocities of liquid and gas particles to the Euler grid such that we can apply the APIC method to the simulation of a two-phase fluid.Second,we introduce a narrowband method to our system because MultiFLIP is a time-consuming approach owing to the large number of particles.Results Experiments show that our method is well integrated with the APIC method and provides a visually credible two-phase fluid animation.Conclusions The proposed method can successfully handle the simulation of a two phase fluid.
文摘An explicit,time-dependent variable grid finite difference method is introduced and analyzed for approximating the solution of a scalar conservation law in two dimension. The scheme is stable,and the numerical solution is proved to converges to the relevant physical solution.