A new coupled finite element formulation is proposed to calculate a conformation tensor model in two complex flows:a planar contraction flow and a planar flow around a symmetrically placed cylinder.The components of c...A new coupled finite element formulation is proposed to calculate a conformation tensor model in two complex flows:a planar contraction flow and a planar flow around a symmetrically placed cylinder.The components of conformation tensor are first computed together with the velocity and pressure to describe the change of morphology of polymer chain coils in flow fields.Macroscopic quantities of viscoelastic flow are then calculated based on the conformation tensor.Comparisons between the numerical simulations and experiments for stress patterns and velocity profiles are carried out to prove the validity of the method.展开更多
Three types of previously used numerical methods are revisited for computing the streamfunctionψand velocity potentialχfrom the horizontal velocity v in limited domains.The first type,called the SOR-based method,use...Three types of previously used numerical methods are revisited for computing the streamfunctionψand velocity potentialχfrom the horizontal velocity v in limited domains.The first type,called the SOR-based method,uses a classical successive over-relaxation(SOR)scheme to computeψ(orχ)first with an arbitrary boundary condition(BC)and thenχ(orψ)with the BC derived from v.The second type,called the spectral method,uses spectral formulations to construct the inner part of(ψ,χ)-the inversion of(vorticity,divergence)with a homogeneous BC,and then the remaining harmonic part of(ψ,χ)with BCs from v.The third type,called the integral method,uses integral formulas to compute the internally induced(ψ,χ)-the inversion of domain-internal(vorticity,divergence)using the free-space Greenꞌs function without BCs and then the remaining harmonicψ(orχ)with BCs from v minus the internally-induced part.Although these methods have previously been successfully applied to flows in large-scale and synoptic-scale domains,their accuracy is compromised when applied to complex flows over mesoscale domains,as shown in this paper.To resolve this problem,two hybrid approaches,the integral-SOR method and the integral-spectral method,are developed by combining the first step of the integral method with the second step adopted from the SOR-based and spectral methods,respectively.Upon testing these methods on real-case complex flows,the integral-SOR method is significantly more accurate than the integral-spectral method,noting that the latter is still generally more accurate than the three previously-used methods.The integral-SOR method is recommended for future applications and diagnostic studies of complex flows.展开更多
Sediment incipience under flows passing a backward-facing step was studied. A series of experiments were conducted to measure scouring depth, probability of sediment incipience, and instantaneous flow velocity field d...Sediment incipience under flows passing a backward-facing step was studied. A series of experiments were conducted to measure scouring depth, probability of sediment incipience, and instantaneous flow velocity field downstream of a backward-facing step. Instantaneous flow velocity fields were measured by using Particle Image Velocimetry (PIV), and an image processing method for determining probability of sediment incipience was employed to analyze the experimental data. The experimental results showed that the probability of sediment incipience was the highest near the reattachment point, even though the near-wall instantaneous flow velocity and the Reynolds stress were both much higher further downstream of the backward-facing step. The possible me- chanisms are discussed for the sediment incipience near the reattachment point.展开更多
Since proposed,the self-similarity variables based genuinely multidimensional Riemann solver is attracting more attentions due to its high resolution in multidimensional complex flows.However,it needs numerous logical...Since proposed,the self-similarity variables based genuinely multidimensional Riemann solver is attracting more attentions due to its high resolution in multidimensional complex flows.However,it needs numerous logical operations in supersonic cases,which limit the method’s applicability in engineering problems greatly.In order to overcome this defect,a hybrid multidimensional Riemann solver,called HMTHS(Hybrid of MulTv and multidimensional HLL scheme based on Self-similar structures),is proposed.It simulates the strongly interacting zone by adopting the MHLLES(Multidimensional Harten-Lax-van Leer-Eifeldt scheme based on Self-similar structures)scheme at subsonic speeds,which is with a high resolution by considering the second moment in the similarity variables.Also,it adopts the MULTV(Multidimensional Toro and Vasquez)scheme,which is with a high resolution in capturing discontinuities,to simulate the flux at supersonic speeds.Systematic numerical experiments,including both one-dimensional cases and twodimensional cases,are conducted.One-dimensional moving contact discontinuity case and sod shock tube case suggest that HMTHS can accurately capture one-dimensional expansion waves,shock waves,and linear contact discontinuities.Two-dimensional cases,such as the double Mach reflection case,the supersonic shock/boundary layer interaction case,the hypersonic flow over the cylinder case,and the hypersonic viscous flow over the double-ellipsoid case,indicate that the HMTHS scheme is with a high resolution in simulating multidimensional complex flows.Therefore,it is promising to be widely applied in both scholar and engineering areas.展开更多
In this letter,we present a lattice Boltzmann simulation for complex flow in a solar wall system which includes porous media flow and heat transfer,specifically for solar energy utilization through an unglazed transpi...In this letter,we present a lattice Boltzmann simulation for complex flow in a solar wall system which includes porous media flow and heat transfer,specifically for solar energy utilization through an unglazed transpired solar air collector(UTC).Besides the lattice Boltzmann equation(LBE) for time evolution of particle distribution function for fluid field,we introduce an analogy,LBE for time evolution of distribution function for temperature.Both temperature fields of fluid(air) and solid(porous media) are modeled.We study the effects of fan velocity,solar radiation intensity,porosity,etc.on the thermal performance of the UTC.In general,our simulation results are in good agreement with what in literature.With the current system setting,both fan velocity and solar radiation intensity have significant effect on the thermal performance of the UTC.However,it is shown that the porosity has negligible effect on the heat collector indicating the current system setting might not be realistic.Further examinations of thermal performance in different UTC systems are ongoing.The results are expected to present in near future.展开更多
The multi-scale structures of complex flows have been great challenges to both theoretical and engineer-ing researches, and multi-scale modeling is the natural way in response. Particle methods (PMs) are ideal constit...The multi-scale structures of complex flows have been great challenges to both theoretical and engineer-ing researches, and multi-scale modeling is the natural way in response. Particle methods (PMs) are ideal constitutors and powerful probes of multi-scale models, owing to their physical insight and computational simplicity. In this paper, the role of different PMs for multi-scale modeling of complex flows is critically reviewed and possible development of PMs in this background is prospected, with the emphasis on pseudo-particle modeling (PPM). The performances of some different PMs are compared in simulations and new devel-opment in the fundamentals and applications of PPM is also reported, demonstrating PPM as a unique PM for multi-scale modeling.展开更多
Large-eddy simulation(LES) is under its rapid development and is recognized as a possible second generation of CFD methods used in engineering.Large-eddy simulation of two-phase flows and combustion is particularly im...Large-eddy simulation(LES) is under its rapid development and is recognized as a possible second generation of CFD methods used in engineering.Large-eddy simulation of two-phase flows and combustion is particularly important for engineering applications.Some investigators,including the present authors,give their review on LES of spray combustion in gas-turbine combustors and internal combustion engines.However,up to now only a few papers are related to the state-of-the-art on LES of gas-particle flows and combustion.In this paper a review of the advances in LES of complex gas-particle flows and coal combustion is presented.Different sub-grid scale(SGS) stress models and combustion models are described,some of the main results are summarized,and some research needs are discussed.展开更多
Complex flow around floating structures is a highly nonlinear problem,and it is a typical feature in ship and ocean engineering.Traditional experimental methods and potential flow theory have limitations in predicting...Complex flow around floating structures is a highly nonlinear problem,and it is a typical feature in ship and ocean engineering.Traditional experimental methods and potential flow theory have limitations in predicting complex viscous flows.With the improvement of high-performance computing and the development of numerical techniques,computational fluid dynamics(CFD)has become increasingly powerful in predicting the complex viscous flow around floating structures.This paper reviews the recent progress in CFD techniques for numerical solutions of typical complex viscous flows in ship and ocean engineering.Applications to free-surface flows,breaking bow waves of high-speed ship,ship hull-propeller-rudder interaction,vortexinduced vibration of risers,vortex-induced motions of deep-draft platforms,and floating offshore wind turbines are discussed.Typical techniques,including volume of fluid for sharp interface,dynamic overset grid,detached eddy simulation,and fluid-structure coupling,are reviewed along with their applications.Some novel techniques,such as high-efficiency Cartesian grid method and GPU acceleration technique,are discussed in the last part as the future perspective for further enhancement of accuracy and efficiency for CFD simulations of complex flow in ship and ocean engineering.展开更多
Complex networks have established themselves in recent years as being particularly suitable and flexible for representing and modelling many complex natural and artificial systems. Oil-water two-phase flow is one of t...Complex networks have established themselves in recent years as being particularly suitable and flexible for representing and modelling many complex natural and artificial systems. Oil-water two-phase flow is one of the most complex systems. In this paper, we use complex networks to study the inclined oil water two-phase flow. Two different complex network construction methods are proposed to build two types of networks, i.e. the flow pattern complex network （FPCN） and fluid dynamic complex network （FDCN）. Through detecting the community structure of FPCN by the community-detection algorithm based on K-means clustering, useful and interesting results are found which can be used for identifying three inclined oil-water flow patterns. To investigate the dynamic characteristics of the inclined oil-water two-phase flow, we construct 48 FDCNs under different flow conditions, and find that the power-law exponent and the network information entropy, which are sensitive to the flow pattern transition, can both characterize the nonlinear dynamics of the inclined oil-water two-phase flow. In this paper, from a new perspective, we not only introduce a complex network theory into the study of the oil-water two-phase flow but also indicate that the complex network may be a powerful tool for exploring nonlinear time series in practice.展开更多
The incompressible two-phase flows are simulated using combination of an etching multiblock method and a diffuse interface （DI） model, particularly in the com- plex domain that can be decomposed into multiple rectan...The incompressible two-phase flows are simulated using combination of an etching multiblock method and a diffuse interface （DI） model, particularly in the com- plex domain that can be decomposed into multiple rectangular subdomains. The etching multiblock method allows natural communications between the connected subdomains and the efficient parallel computation. The DI model can consider two-phase flows with a large density ratio, and simulate the flows with the moving contact line （MCL） when a geometric formulation of the MCL model is included. Therefore, combination of the etch- ing method and the DI model has potential to deal with a variety of two-phase flows in industrial applications. The performance is examined through a series of numerical exper- iments. The convergence of the etching method is firstly tested by simulating single-phase flows past a square cylinder, and the method for the multiphase flow simulation is vali- dated by investing drops dripping from a pore. The numerical results are compared with either those from other researchers or experimental data. Good agreement is achieved. The method is also used to investigate the impact of a droplet on a grooved substrate and droplet generation in flow focusing devices.展开更多
Viscoelastic flows play an important role in numerous engineering fields,and the multiscale algorithms for simulating viscoelastic flows have received significant attention in order to deepen our understanding of the ...Viscoelastic flows play an important role in numerous engineering fields,and the multiscale algorithms for simulating viscoelastic flows have received significant attention in order to deepen our understanding of the nonlinear dynamic behaviors of viscoelastic fluids.However,traditional grid-based multiscale methods are confined to simple viscoelastic flows with short relaxation time,and there is a lack of uniform multiscale scheme available for coupling different solvers in the simulations of viscoelastic fluids.In this paper,a universal multiscale method coupling an improved smoothed particle hydrodynamics(SPH)and multiscale universal interface(MUI)library is presented for viscoelastic flows.The proposed multiscale method builds on an improved SPH method and leverages the MUI library to facilitate the exchange of information among different solvers in the overlapping domain.We test the capability and flexibility of the presented multiscale method to deal with complex viscoelastic flows by solving different multiscale problems of viscoelastic flows.In the first example,the simulation of a viscoelastic Poiseuille flow is carried out by two coupled improved SPH methods with different spatial resolutions.The effects of exchanging different physical quantities on the numerical results in both the upper and lower domains are also investigated as well as the absolute errors in the overlapping domain.In the second example,the complex Wannier flow with different Weissenberg numbers is further simulated by two improved SPH methods and coupling the improved SPH method and the dissipative particle dynamics(DPD)method.The numerical results show that the physical quantities for viscoelastic flows obtained by the presented multiscale method are in consistence with those obtained by a single solver in the overlapping domain.Moreover,transferring different physical quantities has an important effect on the numerical results.展开更多
The regions with shear stress and mean velocity gradient of opposite sign often exist in complex turbulent shear flows.In these cases,the eddy viscosity hypothesis breaks down.Hinze regards the,departure from eddy vis...The regions with shear stress and mean velocity gradient of opposite sign often exist in complex turbulent shear flows.In these cases,the eddy viscosity hypothesis breaks down.Hinze regards the,departure from eddy viscosity hypothesis as a result from transportation of mean momentum over distance by the large structures and arrives at a shear stress expression including the second order derivatives of the mean velocity.However,his expression greatly overestimates the shear stress.This implies that the flow particles are unlikely to have enough memory of the mean momentum over distance.By assuming the departure from eddy viscosity hypothesis as a result from transportation of the shear stress contained in smaller eddies over distance by the large structures,the present author has arrived at a new shear stress expression.The shear stress estimated so far is in good agreement with the experiments.展开更多
The effects of packing configurations on the phase transition of straight granular chute flow with two bottlenecks axe studied. The granular flow shows a dilute- to-dense flow transition when the channel width is vari...The effects of packing configurations on the phase transition of straight granular chute flow with two bottlenecks axe studied. The granular flow shows a dilute- to-dense flow transition when the channel width is varied, accompanied with a peculiar bistable phenomenon. The bistable phenomenon is induced by the initial packing config- uration of particles. When the packing is dense, the initial flux is small and will induce a dense flow. When the packing is loose, the initial flux is large and will induce a di- lute flow. The fabric network of granulax packing is analyzed from a complex network perspective. The degree distribution shows quantitatively different characteristics for the configurations. A two-dimensional （2D） packing clustering coefficient is defined to better quantify the fabric network.展开更多
Boundary conditions （BCs） play an essential role in lattice Boltzmann （LB） simulations. This paper investigates several most commonly applied BCs by evaluating the relative L2-norm errors of the LB simulations for...Boundary conditions （BCs） play an essential role in lattice Boltzmann （LB） simulations. This paper investigates several most commonly applied BCs by evaluating the relative L2-norm errors of the LB simulations for two-dimensional （2-D） Poiseuille flow. It is found that the relative L2-norm error resulting from FHML＇s BC is smaller than that from other BCs as a whole. Then, based on the FHML＇s BC, it formulates an LB model for simulating fluid flows in 2-D channel with complex geometries. Afterwards, the flows between two inclined plates, in a pulmonary blood vessel and in a blood vessel with local expansion region, are simulated. The numerical results are in good agreement with the analytical predictions and clearly show that the model is effective. It is expected that the model can be extended to simulate some real biologic flows, such as blood flows in arteries, vessels with stenosises, aneurysms and bifurcations,展开更多
A new centrality measure for complex networks, called resource flow centrality, is pro- posed in this paper. This centrality measure is based on the concept of the resource flow in net- works. It not only can be appli...A new centrality measure for complex networks, called resource flow centrality, is pro- posed in this paper. This centrality measure is based on the concept of the resource flow in net- works. It not only can be applied to the connected networks, but also the disconnected networks. Moreover, it overcomes some disadvantages of several common centrality measures. The perform- ance of the proposed measure is compared with some standard centrality measures using a classic dataset and the results indicate the proposed measure performs more reasonably. The statistical dis- tribution of the proposed centrality is investigated by experiments on large scale computer generated graphs and two networks from the real world.展开更多
We extend the complexity entropy causality plane（CECP） to propose a multi-scale complexity entropy causality plane（MS-CECP） and further use the proposed method to discriminate the deterministic characteristics of ...We extend the complexity entropy causality plane（CECP） to propose a multi-scale complexity entropy causality plane（MS-CECP） and further use the proposed method to discriminate the deterministic characteristics of different oil-in-water flows. We first take several typical time series for example to investigate the characteristic of the MS-CECP and find that the MS-CECP not only describes the continuous loss of dynamical structure with the increase of scale, but also reflects the determinacy of the system. Then we calculate the MS-CECP for the conductance fluctuating signals measured from oil–water two-phase flow loop test facility. The results indicate that the MS-CECP could be an intrinsic measure for indicating oil-in-water two-phase flow structures.展开更多
In this paper, we mainly consider the complexity of the k-splittable flow minimizing congestion problem. We give some complexity results. For the k-splittable flow problem, the existence of a feasible solution is stro...In this paper, we mainly consider the complexity of the k-splittable flow minimizing congestion problem. We give some complexity results. For the k-splittable flow problem, the existence of a feasible solution is strongly NP-hard. When the number of the source nodes is an input, for the uniformly exactly k-splittable flow problem, obtaining an approximation algorithm with performance ratio better than (√5+1)/2 is NP-hard. When k is an input, for single commodity k-splittable flow problem, obtaining an algorithm with performance ratio better than is NP-hard. In the last of the paper, we study the relationship of minimizing congestion and minimizing number of rounds in the k-splittable flow problem. The smaller the congestion is, the smaller the number of rounds.展开更多
基金This work was financially supported by the National Natural Science Foundation of China(Nos.20204007 and 50390090)the Doctoral Foundation of National Education Committee of China(No.20030248008)the 863 Project of China(No.2002AA336120).
文摘A new coupled finite element formulation is proposed to calculate a conformation tensor model in two complex flows:a planar contraction flow and a planar flow around a symmetrically placed cylinder.The components of conformation tensor are first computed together with the velocity and pressure to describe the change of morphology of polymer chain coils in flow fields.Macroscopic quantities of viscoelastic flow are then calculated based on the conformation tensor.Comparisons between the numerical simulations and experiments for stress patterns and velocity profiles are carried out to prove the validity of the method.
基金supported by the National Natural Science Foundation of China under Grant Nos. 91937301, 41875074, and 41675060the Second Tibetan Plateau Comprehensive Scientific Expedition 2019QZKK0104+1 种基金the National Key Scientific and Technological Infrastructure Project “EarthLab”provided by NOAA/OAR under NOAA–OU Cooperative Agreement #NA16OAR4320072, U.S. Department of Commerce
文摘Three types of previously used numerical methods are revisited for computing the streamfunctionψand velocity potentialχfrom the horizontal velocity v in limited domains.The first type,called the SOR-based method,uses a classical successive over-relaxation(SOR)scheme to computeψ(orχ)first with an arbitrary boundary condition(BC)and thenχ(orψ)with the BC derived from v.The second type,called the spectral method,uses spectral formulations to construct the inner part of(ψ,χ)-the inversion of(vorticity,divergence)with a homogeneous BC,and then the remaining harmonic part of(ψ,χ)with BCs from v.The third type,called the integral method,uses integral formulas to compute the internally induced(ψ,χ)-the inversion of domain-internal(vorticity,divergence)using the free-space Greenꞌs function without BCs and then the remaining harmonicψ(orχ)with BCs from v minus the internally-induced part.Although these methods have previously been successfully applied to flows in large-scale and synoptic-scale domains,their accuracy is compromised when applied to complex flows over mesoscale domains,as shown in this paper.To resolve this problem,two hybrid approaches,the integral-SOR method and the integral-spectral method,are developed by combining the first step of the integral method with the second step adopted from the SOR-based and spectral methods,respectively.Upon testing these methods on real-case complex flows,the integral-SOR method is significantly more accurate than the integral-spectral method,noting that the latter is still generally more accurate than the three previously-used methods.The integral-SOR method is recommended for future applications and diagnostic studies of complex flows.
基金National Natural Science Foundation of China (No.10602017)Maritime Research Center and DHI-NTU Center of Nanyang Technological University, Singapore
文摘Sediment incipience under flows passing a backward-facing step was studied. A series of experiments were conducted to measure scouring depth, probability of sediment incipience, and instantaneous flow velocity field downstream of a backward-facing step. Instantaneous flow velocity fields were measured by using Particle Image Velocimetry (PIV), and an image processing method for determining probability of sediment incipience was employed to analyze the experimental data. The experimental results showed that the probability of sediment incipience was the highest near the reattachment point, even though the near-wall instantaneous flow velocity and the Reynolds stress were both much higher further downstream of the backward-facing step. The possible me- chanisms are discussed for the sediment incipience near the reattachment point.
基金This study was co-supported by National Natural Science Foundation of China(Nos.11902265 and 11972308)Natural Science Foundation of Shaanxi Province of China(No.2019JQ-376)the Fundamental Research Funds for the Central Universities of China(Nos.G2018KY0304 and G2018KY0308).
文摘Since proposed,the self-similarity variables based genuinely multidimensional Riemann solver is attracting more attentions due to its high resolution in multidimensional complex flows.However,it needs numerous logical operations in supersonic cases,which limit the method’s applicability in engineering problems greatly.In order to overcome this defect,a hybrid multidimensional Riemann solver,called HMTHS(Hybrid of MulTv and multidimensional HLL scheme based on Self-similar structures),is proposed.It simulates the strongly interacting zone by adopting the MHLLES(Multidimensional Harten-Lax-van Leer-Eifeldt scheme based on Self-similar structures)scheme at subsonic speeds,which is with a high resolution by considering the second moment in the similarity variables.Also,it adopts the MULTV(Multidimensional Toro and Vasquez)scheme,which is with a high resolution in capturing discontinuities,to simulate the flux at supersonic speeds.Systematic numerical experiments,including both one-dimensional cases and twodimensional cases,are conducted.One-dimensional moving contact discontinuity case and sod shock tube case suggest that HMTHS can accurately capture one-dimensional expansion waves,shock waves,and linear contact discontinuities.Two-dimensional cases,such as the double Mach reflection case,the supersonic shock/boundary layer interaction case,the hypersonic flow over the cylinder case,and the hypersonic viscous flow over the double-ellipsoid case,indicate that the HMTHS scheme is with a high resolution in simulating multidimensional complex flows.Therefore,it is promising to be widely applied in both scholar and engineering areas.
基金Supported by the National Nature Science Foundation of China under Grant Nos. 10932010,11072220,11072229,U1262109,51176172,and 10972208
文摘In this letter,we present a lattice Boltzmann simulation for complex flow in a solar wall system which includes porous media flow and heat transfer,specifically for solar energy utilization through an unglazed transpired solar air collector(UTC).Besides the lattice Boltzmann equation(LBE) for time evolution of particle distribution function for fluid field,we introduce an analogy,LBE for time evolution of distribution function for temperature.Both temperature fields of fluid(air) and solid(porous media) are modeled.We study the effects of fan velocity,solar radiation intensity,porosity,etc.on the thermal performance of the UTC.In general,our simulation results are in good agreement with what in literature.With the current system setting,both fan velocity and solar radiation intensity have significant effect on the thermal performance of the UTC.However,it is shown that the porosity has negligible effect on the heat collector indicating the current system setting might not be realistic.Further examinations of thermal performance in different UTC systems are ongoing.The results are expected to present in near future.
基金the National Natural Science Foundation of China(Grant Nos.20336040 , 20490201)the Chinese Academy of Sciences(Grant No.INF105-SCE-2-07).
文摘The multi-scale structures of complex flows have been great challenges to both theoretical and engineer-ing researches, and multi-scale modeling is the natural way in response. Particle methods (PMs) are ideal constitutors and powerful probes of multi-scale models, owing to their physical insight and computational simplicity. In this paper, the role of different PMs for multi-scale modeling of complex flows is critically reviewed and possible development of PMs in this background is prospected, with the emphasis on pseudo-particle modeling (PPM). The performances of some different PMs are compared in simulations and new devel-opment in the fundamentals and applications of PPM is also reported, demonstrating PPM as a unique PM for multi-scale modeling.
基金Supported by the National Natural Science Foundation of China (50606026,50736006)the Foundation of State Key Laboratory of Engines,Tianjin University (K-2010-07)
文摘Large-eddy simulation(LES) is under its rapid development and is recognized as a possible second generation of CFD methods used in engineering.Large-eddy simulation of two-phase flows and combustion is particularly important for engineering applications.Some investigators,including the present authors,give their review on LES of spray combustion in gas-turbine combustors and internal combustion engines.However,up to now only a few papers are related to the state-of-the-art on LES of gas-particle flows and combustion.In this paper a review of the advances in LES of complex gas-particle flows and coal combustion is presented.Different sub-grid scale(SGS) stress models and combustion models are described,some of the main results are summarized,and some research needs are discussed.
基金supported by the National Natural Science Foundation of China(51809169,51879159)Chang Jiang Scholars Program(T2014099)+2 种基金Shanghai Excellent Academic Leaders Program(17XD1402300)Innovative Special Project of Numerical Tank of Ministry of Industry and Information Technology of China(2016-23/09)National Key Research and Development Program of China(2019YFB1704203,2019YFC0312400).
文摘Complex flow around floating structures is a highly nonlinear problem,and it is a typical feature in ship and ocean engineering.Traditional experimental methods and potential flow theory have limitations in predicting complex viscous flows.With the improvement of high-performance computing and the development of numerical techniques,computational fluid dynamics(CFD)has become increasingly powerful in predicting the complex viscous flow around floating structures.This paper reviews the recent progress in CFD techniques for numerical solutions of typical complex viscous flows in ship and ocean engineering.Applications to free-surface flows,breaking bow waves of high-speed ship,ship hull-propeller-rudder interaction,vortexinduced vibration of risers,vortex-induced motions of deep-draft platforms,and floating offshore wind turbines are discussed.Typical techniques,including volume of fluid for sharp interface,dynamic overset grid,detached eddy simulation,and fluid-structure coupling,are reviewed along with their applications.Some novel techniques,such as high-efficiency Cartesian grid method and GPU acceleration technique,are discussed in the last part as the future perspective for further enhancement of accuracy and efficiency for CFD simulations of complex flow in ship and ocean engineering.
基金Project supported by the National Natural Science Foundation of China (Grant Nos 50674070 and 60374041)the National High Technology Research and Development Program of China (Grant No 2007AA06Z231)
文摘Complex networks have established themselves in recent years as being particularly suitable and flexible for representing and modelling many complex natural and artificial systems. Oil-water two-phase flow is one of the most complex systems. In this paper, we use complex networks to study the inclined oil water two-phase flow. Two different complex network construction methods are proposed to build two types of networks, i.e. the flow pattern complex network （FPCN） and fluid dynamic complex network （FDCN）. Through detecting the community structure of FPCN by the community-detection algorithm based on K-means clustering, useful and interesting results are found which can be used for identifying three inclined oil-water flow patterns. To investigate the dynamic characteristics of the inclined oil-water two-phase flow, we construct 48 FDCNs under different flow conditions, and find that the power-law exponent and the network information entropy, which are sensitive to the flow pattern transition, can both characterize the nonlinear dynamics of the inclined oil-water two-phase flow. In this paper, from a new perspective, we not only introduce a complex network theory into the study of the oil-water two-phase flow but also indicate that the complex network may be a powerful tool for exploring nonlinear time series in practice.
基金Project supported by the National Natural Science Foundation of China(No.11425210)the Fundamental Research Funds for the Central Universities(No.WK2090050025)
文摘The incompressible two-phase flows are simulated using combination of an etching multiblock method and a diffuse interface （DI） model, particularly in the com- plex domain that can be decomposed into multiple rectangular subdomains. The etching multiblock method allows natural communications between the connected subdomains and the efficient parallel computation. The DI model can consider two-phase flows with a large density ratio, and simulate the flows with the moving contact line （MCL） when a geometric formulation of the MCL model is included. Therefore, combination of the etch- ing method and the DI model has potential to deal with a variety of two-phase flows in industrial applications. The performance is examined through a series of numerical exper- iments. The convergence of the etching method is firstly tested by simulating single-phase flows past a square cylinder, and the method for the multiphase flow simulation is vali- dated by investing drops dripping from a pore. The numerical results are compared with either those from other researchers or experimental data. Good agreement is achieved. The method is also used to investigate the impact of a droplet on a grooved substrate and droplet generation in flow focusing devices.
基金Project supported by the National Natural Science Foundation of China(No.52109068)the Water Conservancy Technology Project of Jiangsu Province of China(No.2022060)。
文摘Viscoelastic flows play an important role in numerous engineering fields,and the multiscale algorithms for simulating viscoelastic flows have received significant attention in order to deepen our understanding of the nonlinear dynamic behaviors of viscoelastic fluids.However,traditional grid-based multiscale methods are confined to simple viscoelastic flows with short relaxation time,and there is a lack of uniform multiscale scheme available for coupling different solvers in the simulations of viscoelastic fluids.In this paper,a universal multiscale method coupling an improved smoothed particle hydrodynamics(SPH)and multiscale universal interface(MUI)library is presented for viscoelastic flows.The proposed multiscale method builds on an improved SPH method and leverages the MUI library to facilitate the exchange of information among different solvers in the overlapping domain.We test the capability and flexibility of the presented multiscale method to deal with complex viscoelastic flows by solving different multiscale problems of viscoelastic flows.In the first example,the simulation of a viscoelastic Poiseuille flow is carried out by two coupled improved SPH methods with different spatial resolutions.The effects of exchanging different physical quantities on the numerical results in both the upper and lower domains are also investigated as well as the absolute errors in the overlapping domain.In the second example,the complex Wannier flow with different Weissenberg numbers is further simulated by two improved SPH methods and coupling the improved SPH method and the dissipative particle dynamics(DPD)method.The numerical results show that the physical quantities for viscoelastic flows obtained by the presented multiscale method are in consistence with those obtained by a single solver in the overlapping domain.Moreover,transferring different physical quantities has an important effect on the numerical results.
文摘The regions with shear stress and mean velocity gradient of opposite sign often exist in complex turbulent shear flows.In these cases,the eddy viscosity hypothesis breaks down.Hinze regards the,departure from eddy viscosity hypothesis as a result from transportation of mean momentum over distance by the large structures and arrives at a shear stress expression including the second order derivatives of the mean velocity.However,his expression greatly overestimates the shear stress.This implies that the flow particles are unlikely to have enough memory of the mean momentum over distance.By assuming the departure from eddy viscosity hypothesis as a result from transportation of the shear stress contained in smaller eddies over distance by the large structures,the present author has arrived at a new shear stress expression.The shear stress estimated so far is in good agreement with the experiments.
基金Project supported by the National Natural Science Foundation of China(Nos.11034010 and 71171185)
文摘The effects of packing configurations on the phase transition of straight granular chute flow with two bottlenecks axe studied. The granular flow shows a dilute- to-dense flow transition when the channel width is varied, accompanied with a peculiar bistable phenomenon. The bistable phenomenon is induced by the initial packing config- uration of particles. When the packing is dense, the initial flux is small and will induce a dense flow. When the packing is loose, the initial flux is large and will induce a di- lute flow. The fabric network of granulax packing is analyzed from a complex network perspective. The degree distribution shows quantitatively different characteristics for the configurations. A two-dimensional （2D） packing clustering coefficient is defined to better quantify the fabric network.
基金Project supported by the National Natural Science Foundation of China (Grant No 10765002)Guangxi Natural Science Foundation (Grant No 0542045)
文摘Boundary conditions （BCs） play an essential role in lattice Boltzmann （LB） simulations. This paper investigates several most commonly applied BCs by evaluating the relative L2-norm errors of the LB simulations for two-dimensional （2-D） Poiseuille flow. It is found that the relative L2-norm error resulting from FHML＇s BC is smaller than that from other BCs as a whole. Then, based on the FHML＇s BC, it formulates an LB model for simulating fluid flows in 2-D channel with complex geometries. Afterwards, the flows between two inclined plates, in a pulmonary blood vessel and in a blood vessel with local expansion region, are simulated. The numerical results are in good agreement with the analytical predictions and clearly show that the model is effective. It is expected that the model can be extended to simulate some real biologic flows, such as blood flows in arteries, vessels with stenosises, aneurysms and bifurcations,
基金Supported by the National Natural Science Foundation of China(61272119,61203372)
文摘A new centrality measure for complex networks, called resource flow centrality, is pro- posed in this paper. This centrality measure is based on the concept of the resource flow in net- works. It not only can be applied to the connected networks, but also the disconnected networks. Moreover, it overcomes some disadvantages of several common centrality measures. The perform- ance of the proposed measure is compared with some standard centrality measures using a classic dataset and the results indicate the proposed measure performs more reasonably. The statistical dis- tribution of the proposed centrality is investigated by experiments on large scale computer generated graphs and two networks from the real world.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.41174109 and 61104148)the National Science and Technology Major Project of China(Grant No.2011ZX05020-006)the Zhejiang Key Discipline of Instrument Science and Technology,China(Grant No.JL130106)
文摘We extend the complexity entropy causality plane（CECP） to propose a multi-scale complexity entropy causality plane（MS-CECP） and further use the proposed method to discriminate the deterministic characteristics of different oil-in-water flows. We first take several typical time series for example to investigate the characteristic of the MS-CECP and find that the MS-CECP not only describes the continuous loss of dynamical structure with the increase of scale, but also reflects the determinacy of the system. Then we calculate the MS-CECP for the conductance fluctuating signals measured from oil–water two-phase flow loop test facility. The results indicate that the MS-CECP could be an intrinsic measure for indicating oil-in-water two-phase flow structures.
文摘In this paper, we mainly consider the complexity of the k-splittable flow minimizing congestion problem. We give some complexity results. For the k-splittable flow problem, the existence of a feasible solution is strongly NP-hard. When the number of the source nodes is an input, for the uniformly exactly k-splittable flow problem, obtaining an approximation algorithm with performance ratio better than (√5+1)/2 is NP-hard. When k is an input, for single commodity k-splittable flow problem, obtaining an algorithm with performance ratio better than is NP-hard. In the last of the paper, we study the relationship of minimizing congestion and minimizing number of rounds in the k-splittable flow problem. The smaller the congestion is, the smaller the number of rounds.
