This article reports recent developments and advances in the simulation of the CO2-formation fluid displacement behaviour at the pore scale of subsurface porous media. Roughly, there are three effective visualization ...This article reports recent developments and advances in the simulation of the CO2-formation fluid displacement behaviour at the pore scale of subsurface porous media. Roughly, there are three effective visualization approaches to detect and observe the CO2-formation fluid displacement mechanism at the micro-scale, namely, magnetic resonance imaging, X-ray computed tomography and fabricated micromodels, but they are not capable of investigating the dis- placement process at the nano-scale. Though a lab-on-chip approach for the direct visualization of the fluid flow behaviour in nanoscale channels has been developed using an advanced epi-fluorescence microscopy method combined with a nanofluidic chip, it is still a qualitative analysis method. The lattice Boltzmann method (LBM) can simulate the CO2 displacement processes in a two-dimensional or three-dimensional (3D) pore structure, but until now, the CO2 displace- ment mechanisms had not been thoroughly investigated and the 3D pore structure of real rock had not been directly taken into account in the simulation of the CO2 displacement process. The status of research on the applications of CO2 displacement to enhance shale gas recovery is also analyzed in this paper. The coupling of molecular dynamics and LBM in tandem is proposed to simulate the CO2-shale gas displacement process based on the 3D digital model of shale obtained from focused ion beams and scanning electron microscopy.展开更多
To promote/inhibit ice formation in the natural environment and industrial systems,the growth and evolution process of ice single-crystal nuclei were simulated using the phase field-lattice Boltzmann method(PF-LBM),an...To promote/inhibit ice formation in the natural environment and industrial systems,the growth and evolution process of ice single-crystal nuclei were simulated using the phase field-lattice Boltzmann method(PF-LBM),and the influence of a background flow field on the growth of single-crystal nucleus dendrites was also analyzed.The results show that the flow field makes dendrite growth asymmetric.The growth of dendrites is more developed on the upstream side than on the downstream side.The dendrite tip growth rate and tip radius are greater on the upstream side than on the downstream side.The solid phase ratio is greater with a background flow field than without one.The higher the flow velocity is,the more developed the dendrites on the upstream side,the faster the dendrites grow,and the higher the dendrite tip growth rate.The dendrites on the backflow side have a lower flow rate and a lower degree of supercooling than those on the upstream side,which inhibits the solidification process,the growth rate is slow,and the dendrites are underdeveloped.展开更多
This paper presents a KIVA-3 code based numerical model for three-dimensional transient intake flow in the intake port-valve-cylinder system of internal combustion engine using body-fitted technique, which can be used...This paper presents a KIVA-3 code based numerical model for three-dimensional transient intake flow in the intake port-valve-cylinder system of internal combustion engine using body-fitted technique, which can be used in numerical study on internal combustion engine with vertical and inclined valves, and has higher calculation precision. A numerical simulation (on the intake process of a two-valve engine with a semi-sphere combustion chamber and a radial intake port) is provided for analysis of the velocity field and pressure field of different plane at different crank angles. The results revealed the formation of the tumble motion, the evolution of flow field parameters and the variation of tumble ratios as important information for the design of engine in-take system.展开更多
The lattice Boltzmann method (LBM) is employed to simulate the uniform flow past a circular cylinder. The performance of the two-dimensional LBM model on the prediction of force coefficients and vortex shedding freque...The lattice Boltzmann method (LBM) is employed to simulate the uniform flow past a circular cylinder. The performance of the two-dimensional LBM model on the prediction of force coefficients and vortex shedding frequency is investigated. The local grid refinement technique and second-order boundary condition for curved walls are applied in the calculations. It is found that the calculated vortex shedding frequency, drag coefficient and lift coefficient are consistent with experimental results at Reynolds nu...展开更多
The numerical simulation using the multiple relaxation time lattice Boltzmann method (MRT-LBM) is carried out for the purpose of investigating the two-dimensional flow around three circular cylinders. Among these th...The numerical simulation using the multiple relaxation time lattice Boltzmann method (MRT-LBM) is carried out for the purpose of investigating the two-dimensional flow around three circular cylinders. Among these three circular cylinders, one of the three cylinders on which a forced in-line vibrating is used to do this research and attempt to find out the effects of the moving cylinder and the other two rigid cylinders on the wake characteristics and vortex formation. As a benchmark problem to discuss the problem of lift coefficient r.m.s for these cylinders with spacing ratios T/ D between other rigid side-by-side cylinders, and the calculation is carried out with two compared cases at Reynolds number of 100, two of the cylinders are rigid and the other one is an in-line vibrated cylinder lying downstream, in addition, forced vibrating amplitude and frequency are A/D = 0.5 and fv= 0.4 (where A is the forced amplitude, D is the cylinder diameter, and fv stands for the vibrating frequency, respectively). The calculated results not only indicate that the spacing ratios T/D (T is the center-to-center spacing between the two upstream cylinders) have influence on the wake patterns and the formation of vortex shedding, but also analyze the lift coefficient r.m.s for the three cylinders with the spacing ratios S/D (where S is the center-to-center spacing between the center of upstream two side-by-side cylinders and downstream cylinder).展开更多
Based on the flow simulation in the condensing zone of compact spinning with lattice apron and a bead-elastic rod dynamic model of the flexible fiber,trajectories of fibers with different negative pressure are simulat...Based on the flow simulation in the condensing zone of compact spinning with lattice apron and a bead-elastic rod dynamic model of the flexible fiber,trajectories of fibers with different negative pressure are simulated by specially designed Matlab procedure.Then displacement components of fibers at YZ profile under different negative pressure conditions are extracted and compared.The results show that the fibers of different initial positions gradually converge,and are interlaced for position change in yarn cross-section,caused by the airflow in the condensing zone.Finally,compact-spun yarn with different negative pressure and conventional ring spun yarn are produced and their twists are tested.Both the results of simulation and experiments illustrate the existence of additional twists.Also the relationship between additional twists and negative pressure is verified.展开更多
The effects of temperature on atomic anti-site behaviors in L12-Ni3(AlFe) phases were studied using microscopic phase-field dynamic model in precipitation progress of Ni75Al20Fe5 alloy.The results show that with the i...The effects of temperature on atomic anti-site behaviors in L12-Ni3(AlFe) phases were studied using microscopic phase-field dynamic model in precipitation progress of Ni75Al20Fe5 alloy.The results show that with the increase of temperature,the formation of NiAl and AlNi anti-sites is much easier in Ni3(AlFe),and Ni and Al anti-site atoms show clearly stronger temperature-dependent than Fe anti-site atoms.The evolution progress of anti-site atoms is completed at the initial growth stage of L12-Ni3(AlFe) phases.The site occupation probabilities of Ni atoms on the sublattice A(NiNi,face centers sites of FCC),and Al and Fe atoms on the sublattice B(AlAl and FeAl,corners sites of FCC) all present the degressive tendency with the temperature increasing.Fe atoms mainly prefer to occupy the Al sublattice at the whole temperature range.展开更多
Numerical simulation of multiphase flows in processing equipment in industry with two-fluid models and Eulerian-Lagrangian approaches requires the constitutive equations describing the interactions between the dispers...Numerical simulation of multiphase flows in processing equipment in industry with two-fluid models and Eulerian-Lagrangian approaches requires the constitutive equations describing the interactions between the dispersed phase of high concentration and the continuous phase. The status of research on the forces on dispersed solid and fluid particles is reviewed in this article. As compared with the knowledge on drag of single solid particles study on panicle swarms and on other forces is not sufficient to meet the demand of reliable and efficient numerical simulation of multiphase flows. Thus, thorough study on the panicle swarms becomes the key to accurate multi-scale simulation of multiphase flows. Besides, the development of efficient algorithm dealing with the non-uniformity on both equipment and mesoscopic scales is recognized as an important issue to be resolved. The research topics in the near future are suggested.展开更多
The adsorption behavior of symmetric triblock copolymers, Am/2BnAm/2, from a nonselective solvent at solid-liquid interface has been studied by Monte Carlo simulations on a simple lattice model. Either segment A or se...The adsorption behavior of symmetric triblock copolymers, Am/2BnAm/2, from a nonselective solvent at solid-liquid interface has been studied by Monte Carlo simulations on a simple lattice model. Either segment A or segment B is attractive, while the other is non-attractive to the surface. Influences of the adsorption energy, bulk concentration, chain composition and chain length on the microstructure of adsorbed layers are presented. The results show that the total surface coverage and the adsorption amount increases monotonically as the bulk concentration increases. The larger the adsorption energy and the higher the fraction of adsorbing segments, the higher the total surface coverage is exhibited. The product of surface coverage and the proportion of non-attractive segments are nearly independent of the chain length, and the logarithm of the adsorption amount is a linear function of the reciprocal of the reduced temperature. When the adsorption energy is larger, the adsorption amount exhibits a maximum as the fraction of adsorbing segment increases. The adsorption isotherms of copolymers with different length of non-attractive segments can be mapped onto a single curve under given adsorption energy. The adsorption layer thickness decreases as the adsorption energy and the fraction of adsorbing segments increases, but it increases as the length of non-attractive segments increases. The tails mainly govern the adsorption layer thickness.展开更多
Phase transition between nematic and isotropic liquid crystal is a very weak first order phase transition.We avoid to use the normal Landau-de Gennes's free energy that reduces a strong first order transition, and...Phase transition between nematic and isotropic liquid crystal is a very weak first order phase transition.We avoid to use the normal Landau-de Gennes's free energy that reduces a strong first order transition, and set up adata base of free energy calculated by means of Tao-Sheng Lin's extended molecular field theory that can explain theexperiments of the equilibrium properties of nematic liquid crystal very well. Then we use the free energy method oflattice Boltzmann developed by Oxford group to study the phase decomposition, pattern formation in the flow of theliquid crystal near transition temperature.展开更多
Hydroelasticity caused by water impact is of concem in many applications of ocean engineering/naval architect and is a complicated physical phenomenon. The authors have developed a coupled Eulerian scheme with Lagrang...Hydroelasticity caused by water impact is of concem in many applications of ocean engineering/naval architect and is a complicated physical phenomenon. The authors have developed a coupled Eulerian scheme with Lagrangian particles to combine advantages and to compensate disadvantages in both grid based method and particle based method. In this study, the developed numerical model was applied to hydroelastic problems due to impact pressure such as water entry of an elastic cylinder and elastic tanker motion in wave. The authors showed the numerical results which is overall agreement with experimental results. The proposed numerical scheme can be useful and effectiveness to evaluate hydroelasticity and ship-wave interaction in nonlinear wave motion with breaking.展开更多
This paper presents the simulation results of the wind environment around a single high-rise building and that around two tall buildings in tandem arrangement by using the lattice Boltzmann method with an aim to under...This paper presents the simulation results of the wind environment around a single high-rise building and that around two tall buildings in tandem arrangement by using the lattice Boltzmann method with an aim to understand the ventilation issues around high-rise buildings in an urban environment.We analyzed the velocity distribution around the buildings and performed numericl simulations to reveal the formation and evolution law of the complex vortex system around the high-rise buildings.Numerical simulation results manifest a periodicity phenamenon in the process of the vortex evolution.For the case of two high-rise buildings,wind velocity in the space between the two buildings is very small,which is nearly a silent regime.Wind velocity above the front building is relatively larger and the maximum wind velocity is approximately 2.5 times the incoming wind velocity.The numerical results can be used in layout planning of high-rise residential buildings to create better environment for ventilation purpose in an urban area.展开更多
A lattice Boltzmann model is presented to simulate the deformation and motions of a red blood cell (RBC) in a shear flow. The curvatures of the membrane of a static RBC with different chemical potentiM drops calcula...A lattice Boltzmann model is presented to simulate the deformation and motions of a red blood cell (RBC) in a shear flow. The curvatures of the membrane of a static RBC with different chemical potentiM drops calculated by our model agree with those computed by a shooting method very well. Our simulation results show that in a shear flow, biconcave RBC becomes highly flattened and undergoes tank-treading motion. With intrinsically parallel dynamics, this lattice Boltzmann method is expected to find wide applications to both single and multi-vesicles suspension as well as complex open membranes in various fluid flows for a wide range of Reynolds numbers.展开更多
Double diffusion convection in a cavity with a hot square obstacle inside is simulated using the lattice Boltzmann method. The results are presented for the Rayleigh numbers 104,105 and 106, the Lewis numbers 0.1, 2 a...Double diffusion convection in a cavity with a hot square obstacle inside is simulated using the lattice Boltzmann method. The results are presented for the Rayleigh numbers 104,105 and 106, the Lewis numbers 0.1, 2 and 10 and aspect ratio A(obstacle height/cavity height) of 0.2, 0.4 and 0.6 for a range of buoyancy number N = 0 to- 4 with the effect of opposing flow. The results indicate that for |N| b 1, the Nusselt and Sherwood numbers decrease as buoyancy ratio increases, while for |N| N 1, they increase with |N|. As the Lewis number increases, higher buoyancy ratio is required to overcome the thermal effects and the minimum value of the Nusselt and Sherwood numbers occur at higher buoyancy ratios. The increase in the Rayleigh or Lewis number results in the formation of the multi-cell flow in the enclosure and the vortices will vanish as |N| increases.展开更多
By using a hybrid lattice-Boltzmann–finite-difference method(hybrid LBM–FDM method),three-dimensional simulations of solutal interfacial convection were conducted for the process of CO2absorption into ethanol.A self...By using a hybrid lattice-Boltzmann–finite-difference method(hybrid LBM–FDM method),three-dimensional simulations of solutal interfacial convection were conducted for the process of CO2absorption into ethanol.A self-renewal interface model is adopted as an interfacial perturbation model.The simulation results revealed some three-dimensional features of the induced interfacial convection,such as the development of diverging cellular flow and Rayleigh plume-like convection in liquid phase.The concentration distribution of the simulation result is validated and found to be in well agreement with the Schlieren visualization results qualitatively.Additionally,the mass transfer enhancements by interfacial convection were investigated via both simulation and experiment for the absorption process,and the mass transfer is shown to be enhanced by the interfacial convection by about two-fold comparing with that by diffusion.展开更多
A family of coupled map lattice (CML) models has been developed to simulate the evolutional mechanism of interactions of convection, diffusion, and dispersion in both weakly and strongly coupled cases. Not only cohe...A family of coupled map lattice (CML) models has been developed to simulate the evolutional mechanism of interactions of convection, diffusion, and dispersion in both weakly and strongly coupled cases. Not only coherent and turbulent properties as well as their relations, but also the transitional states between convection dominating, diffusion dominating and dispersion dominating are analyzed to demonstrate the essential characteristics of any state. Numerical results show that the models are capable of simulating both layered coupling and stochastic mechanism, and lead us to understand whether or not turbulence coherent structure is formed by modulation of wave packet. The duality of wave and particle characters of turbulence is illustrated in the numerical simulation; a sketch picture is given to explain the questions associated with the turbulent inverse cascade, which is the result of the mutual interactions among the physical factors of nonlinearity, dissipation and dispersion.展开更多
基金The authors gratefully acknowledge the financial support of the National Natural Science Foundation of China (Grant Nos. 51374213 and 51674251), the State Key Research Development Program of China (Grant No. 2016YFC0600705), the National Natural Science Fund for Distinguished Young Scholars of China (Grant No. 51125017), the Fund for Innovative Research and Development Group Program of Jiangsu Province (Grant No. 2014- 27), the Science Fund for Creative Research Groups of the National Natural Science Foundation of China (Grant No. 51421003), and the Priority Academic Program Development of the Jiangsu Higher Education Institutions (PAPD 2014).
文摘This article reports recent developments and advances in the simulation of the CO2-formation fluid displacement behaviour at the pore scale of subsurface porous media. Roughly, there are three effective visualization approaches to detect and observe the CO2-formation fluid displacement mechanism at the micro-scale, namely, magnetic resonance imaging, X-ray computed tomography and fabricated micromodels, but they are not capable of investigating the dis- placement process at the nano-scale. Though a lab-on-chip approach for the direct visualization of the fluid flow behaviour in nanoscale channels has been developed using an advanced epi-fluorescence microscopy method combined with a nanofluidic chip, it is still a qualitative analysis method. The lattice Boltzmann method (LBM) can simulate the CO2 displacement processes in a two-dimensional or three-dimensional (3D) pore structure, but until now, the CO2 displace- ment mechanisms had not been thoroughly investigated and the 3D pore structure of real rock had not been directly taken into account in the simulation of the CO2 displacement process. The status of research on the applications of CO2 displacement to enhance shale gas recovery is also analyzed in this paper. The coupling of molecular dynamics and LBM in tandem is proposed to simulate the CO2-shale gas displacement process based on the 3D digital model of shale obtained from focused ion beams and scanning electron microscopy.
基金The National Key Research and Development Program of China(No.2020YFB1506203)。
文摘To promote/inhibit ice formation in the natural environment and industrial systems,the growth and evolution process of ice single-crystal nuclei were simulated using the phase field-lattice Boltzmann method(PF-LBM),and the influence of a background flow field on the growth of single-crystal nucleus dendrites was also analyzed.The results show that the flow field makes dendrite growth asymmetric.The growth of dendrites is more developed on the upstream side than on the downstream side.The dendrite tip growth rate and tip radius are greater on the upstream side than on the downstream side.The solid phase ratio is greater with a background flow field than without one.The higher the flow velocity is,the more developed the dendrites on the upstream side,the faster the dendrites grow,and the higher the dendrite tip growth rate.The dendrites on the backflow side have a lower flow rate and a lower degree of supercooling than those on the upstream side,which inhibits the solidification process,the growth rate is slow,and the dendrites are underdeveloped.
文摘This paper presents a KIVA-3 code based numerical model for three-dimensional transient intake flow in the intake port-valve-cylinder system of internal combustion engine using body-fitted technique, which can be used in numerical study on internal combustion engine with vertical and inclined valves, and has higher calculation precision. A numerical simulation (on the intake process of a two-valve engine with a semi-sphere combustion chamber and a radial intake port) is provided for analysis of the velocity field and pressure field of different plane at different crank angles. The results revealed the formation of the tumble motion, the evolution of flow field parameters and the variation of tumble ratios as important information for the design of engine in-take system.
基金Supported by the Specialized Research Fund for the Doctoral Program of Higher Education (No. 20060056036)
文摘The lattice Boltzmann method (LBM) is employed to simulate the uniform flow past a circular cylinder. The performance of the two-dimensional LBM model on the prediction of force coefficients and vortex shedding frequency is investigated. The local grid refinement technique and second-order boundary condition for curved walls are applied in the calculations. It is found that the calculated vortex shedding frequency, drag coefficient and lift coefficient are consistent with experimental results at Reynolds nu...
基金Support by the National Natural Science Foundation of China under Grant Nos.10932010 and 11072220the Natural Science Foundation of Zhejiang Province under Grant Nos.Y607425,Z6090556the Foundation Project for Youths of Zhijiang Normal University under Grant No.KJ20090102
文摘The numerical simulation using the multiple relaxation time lattice Boltzmann method (MRT-LBM) is carried out for the purpose of investigating the two-dimensional flow around three circular cylinders. Among these three circular cylinders, one of the three cylinders on which a forced in-line vibrating is used to do this research and attempt to find out the effects of the moving cylinder and the other two rigid cylinders on the wake characteristics and vortex formation. As a benchmark problem to discuss the problem of lift coefficient r.m.s for these cylinders with spacing ratios T/ D between other rigid side-by-side cylinders, and the calculation is carried out with two compared cases at Reynolds number of 100, two of the cylinders are rigid and the other one is an in-line vibrated cylinder lying downstream, in addition, forced vibrating amplitude and frequency are A/D = 0.5 and fv= 0.4 (where A is the forced amplitude, D is the cylinder diameter, and fv stands for the vibrating frequency, respectively). The calculated results not only indicate that the spacing ratios T/D (T is the center-to-center spacing between the two upstream cylinders) have influence on the wake patterns and the formation of vortex shedding, but also analyze the lift coefficient r.m.s for the three cylinders with the spacing ratios S/D (where S is the center-to-center spacing between the center of upstream two side-by-side cylinders and downstream cylinder).
文摘Based on the flow simulation in the condensing zone of compact spinning with lattice apron and a bead-elastic rod dynamic model of the flexible fiber,trajectories of fibers with different negative pressure are simulated by specially designed Matlab procedure.Then displacement components of fibers at YZ profile under different negative pressure conditions are extracted and compared.The results show that the fibers of different initial positions gradually converge,and are interlaced for position change in yarn cross-section,caused by the airflow in the condensing zone.Finally,compact-spun yarn with different negative pressure and conventional ring spun yarn are produced and their twists are tested.Both the results of simulation and experiments illustrate the existence of additional twists.Also the relationship between additional twists and negative pressure is verified.
基金Project(50671084) supported by the National Natural Science Foundation of ChinaProject(2009021028) supported by Science and Technique Foundation for Young Scholars of Shanxi Province, ChinaProject(20100470125) supported by National Science Foundation for Post-doctoral Scientists of China
文摘The effects of temperature on atomic anti-site behaviors in L12-Ni3(AlFe) phases were studied using microscopic phase-field dynamic model in precipitation progress of Ni75Al20Fe5 alloy.The results show that with the increase of temperature,the formation of NiAl and AlNi anti-sites is much easier in Ni3(AlFe),and Ni and Al anti-site atoms show clearly stronger temperature-dependent than Fe anti-site atoms.The evolution progress of anti-site atoms is completed at the initial growth stage of L12-Ni3(AlFe) phases.The site occupation probabilities of Ni atoms on the sublattice A(NiNi,face centers sites of FCC),and Al and Fe atoms on the sublattice B(AlAl and FeAl,corners sites of FCC) all present the degressive tendency with the temperature increasing.Fe atoms mainly prefer to occupy the Al sublattice at the whole temperature range.
基金Supported by the National Natural Science Foundation of China (20490206, 20676134), the National Basic Research Program of China (2009CB623406, 2004CB217604), the National Project of Scientific and Technical Supporting Program (2008BAF33B03) and the National High Technology Research and Development Program of China (2007AA060904)..
文摘Numerical simulation of multiphase flows in processing equipment in industry with two-fluid models and Eulerian-Lagrangian approaches requires the constitutive equations describing the interactions between the dispersed phase of high concentration and the continuous phase. The status of research on the forces on dispersed solid and fluid particles is reviewed in this article. As compared with the knowledge on drag of single solid particles study on panicle swarms and on other forces is not sufficient to meet the demand of reliable and efficient numerical simulation of multiphase flows. Thus, thorough study on the panicle swarms becomes the key to accurate multi-scale simulation of multiphase flows. Besides, the development of efficient algorithm dealing with the non-uniformity on both equipment and mesoscopic scales is recognized as an important issue to be resolved. The research topics in the near future are suggested.
基金the National Natural Science Foundation of China (No. 20025618, No. 20236010) Shanghai Municipal Education Commission of China.
文摘The adsorption behavior of symmetric triblock copolymers, Am/2BnAm/2, from a nonselective solvent at solid-liquid interface has been studied by Monte Carlo simulations on a simple lattice model. Either segment A or segment B is attractive, while the other is non-attractive to the surface. Influences of the adsorption energy, bulk concentration, chain composition and chain length on the microstructure of adsorbed layers are presented. The results show that the total surface coverage and the adsorption amount increases monotonically as the bulk concentration increases. The larger the adsorption energy and the higher the fraction of adsorbing segments, the higher the total surface coverage is exhibited. The product of surface coverage and the proportion of non-attractive segments are nearly independent of the chain length, and the logarithm of the adsorption amount is a linear function of the reciprocal of the reduced temperature. When the adsorption energy is larger, the adsorption amount exhibits a maximum as the fraction of adsorbing segment increases. The adsorption isotherms of copolymers with different length of non-attractive segments can be mapped onto a single curve under given adsorption energy. The adsorption layer thickness decreases as the adsorption energy and the fraction of adsorbing segments increases, but it increases as the length of non-attractive segments increases. The tails mainly govern the adsorption layer thickness.
文摘Phase transition between nematic and isotropic liquid crystal is a very weak first order phase transition.We avoid to use the normal Landau-de Gennes's free energy that reduces a strong first order transition, and set up adata base of free energy calculated by means of Tao-Sheng Lin's extended molecular field theory that can explain theexperiments of the equilibrium properties of nematic liquid crystal very well. Then we use the free energy method oflattice Boltzmann developed by Oxford group to study the phase decomposition, pattern formation in the flow of theliquid crystal near transition temperature.
文摘Hydroelasticity caused by water impact is of concem in many applications of ocean engineering/naval architect and is a complicated physical phenomenon. The authors have developed a coupled Eulerian scheme with Lagrangian particles to combine advantages and to compensate disadvantages in both grid based method and particle based method. In this study, the developed numerical model was applied to hydroelastic problems due to impact pressure such as water entry of an elastic cylinder and elastic tanker motion in wave. The authors showed the numerical results which is overall agreement with experimental results. The proposed numerical scheme can be useful and effectiveness to evaluate hydroelasticity and ship-wave interaction in nonlinear wave motion with breaking.
基金Funded by the Natural Science Foundation of Tianjin Municipality (No. 06YFJMJC05300)the Science and Technology Development Foundation for Universities of Tianjin Municipality (No. 20060823)
文摘This paper presents the simulation results of the wind environment around a single high-rise building and that around two tall buildings in tandem arrangement by using the lattice Boltzmann method with an aim to understand the ventilation issues around high-rise buildings in an urban environment.We analyzed the velocity distribution around the buildings and performed numericl simulations to reveal the formation and evolution law of the complex vortex system around the high-rise buildings.Numerical simulation results manifest a periodicity phenamenon in the process of the vortex evolution.For the case of two high-rise buildings,wind velocity in the space between the two buildings is very small,which is nearly a silent regime.Wind velocity above the front building is relatively larger and the maximum wind velocity is approximately 2.5 times the incoming wind velocity.The numerical results can be used in layout planning of high-rise residential buildings to create better environment for ventilation purpose in an urban area.
基金supported by National Natural Science Foundation of China under Grant No. 10747004the Guangxi Science Foundation under Grant Nos. 0640064 and 0542045
文摘A lattice Boltzmann model is presented to simulate the deformation and motions of a red blood cell (RBC) in a shear flow. The curvatures of the membrane of a static RBC with different chemical potentiM drops calculated by our model agree with those computed by a shooting method very well. Our simulation results show that in a shear flow, biconcave RBC becomes highly flattened and undergoes tank-treading motion. With intrinsically parallel dynamics, this lattice Boltzmann method is expected to find wide applications to both single and multi-vesicles suspension as well as complex open membranes in various fluid flows for a wide range of Reynolds numbers.
文摘Double diffusion convection in a cavity with a hot square obstacle inside is simulated using the lattice Boltzmann method. The results are presented for the Rayleigh numbers 104,105 and 106, the Lewis numbers 0.1, 2 and 10 and aspect ratio A(obstacle height/cavity height) of 0.2, 0.4 and 0.6 for a range of buoyancy number N = 0 to- 4 with the effect of opposing flow. The results indicate that for |N| b 1, the Nusselt and Sherwood numbers decrease as buoyancy ratio increases, while for |N| N 1, they increase with |N|. As the Lewis number increases, higher buoyancy ratio is required to overcome the thermal effects and the minimum value of the Nusselt and Sherwood numbers occur at higher buoyancy ratios. The increase in the Rayleigh or Lewis number results in the formation of the multi-cell flow in the enclosure and the vortices will vanish as |N| increases.
基金Supported by the National Natural Science Foundation of China(20736005)
文摘By using a hybrid lattice-Boltzmann–finite-difference method(hybrid LBM–FDM method),three-dimensional simulations of solutal interfacial convection were conducted for the process of CO2absorption into ethanol.A self-renewal interface model is adopted as an interfacial perturbation model.The simulation results revealed some three-dimensional features of the induced interfacial convection,such as the development of diverging cellular flow and Rayleigh plume-like convection in liquid phase.The concentration distribution of the simulation result is validated and found to be in well agreement with the Schlieren visualization results qualitatively.Additionally,the mass transfer enhancements by interfacial convection were investigated via both simulation and experiment for the absorption process,and the mass transfer is shown to be enhanced by the interfacial convection by about two-fold comparing with that by diffusion.
基金supported by National Natural Science Foundation of China under Grant No.40535025
文摘A family of coupled map lattice (CML) models has been developed to simulate the evolutional mechanism of interactions of convection, diffusion, and dispersion in both weakly and strongly coupled cases. Not only coherent and turbulent properties as well as their relations, but also the transitional states between convection dominating, diffusion dominating and dispersion dominating are analyzed to demonstrate the essential characteristics of any state. Numerical results show that the models are capable of simulating both layered coupling and stochastic mechanism, and lead us to understand whether or not turbulence coherent structure is formed by modulation of wave packet. The duality of wave and particle characters of turbulence is illustrated in the numerical simulation; a sketch picture is given to explain the questions associated with the turbulent inverse cascade, which is the result of the mutual interactions among the physical factors of nonlinearity, dissipation and dispersion.