The Bcklund transformation and the generalized Miura transformation for the Volterra lattice equation are constructed by using point symmetry method.As an application,the explicit solution to the lattice equation is...The Bcklund transformation and the generalized Miura transformation for the Volterra lattice equation are constructed by using point symmetry method.As an application,the explicit solution to the lattice equation is obtained.展开更多
A discrete isospectral problem and the associated hierarchy of Lax integrable lattice equations were investigated. A Darboux transformation for the discrete spectral problem was found. Finally, an infinite number of c...A discrete isospectral problem and the associated hierarchy of Lax integrable lattice equations were investigated. A Darboux transformation for the discrete spectral problem was found. Finally, an infinite number of conservation laws were given for the corresponding hierarchy.展开更多
A direct way to construct integrable couplings for discrete systems is presented by use of two semi-direct sum Lie algebras. As their applications, the discrete integrable couplings associated with modified Korteweg-d...A direct way to construct integrable couplings for discrete systems is presented by use of two semi-direct sum Lie algebras. As their applications, the discrete integrable couplings associated with modified Korteweg-de Vries (m-KdV) lattice and two hierarchies of discrete soliton equations are developed. It is also indicated that the study of integrable couplings using semi-direct sums of Lie algebras is an important step towards the complete classification of integrable couplings.展开更多
This paper focuses on studying the symmetry of a practical wave equation on new lattices. It is a new step in that the new lattice equation is applied to reduce the discrete problem of motion of an elastic thin homoge...This paper focuses on studying the symmetry of a practical wave equation on new lattices. It is a new step in that the new lattice equation is applied to reduce the discrete problem of motion of an elastic thin homogeneous bar. The equation of motion of the bar can be changed into a discrete wave equation. With the new lattice equation, the translational and scaling invariant, not only is the infinitesimal transformation given, but the symmetry and Lie algebras are also calculated. We also give a new form of invariant called the ratio invariant, which can reduce the process of the computing invariant with the characteristic equation.展开更多
In this paper, we propose a novel incompressible finite-difference lattice Boltzmann Equation (FDLBE). Because source terms that reflect the interaction between phases can be accurately described, the new model is s...In this paper, we propose a novel incompressible finite-difference lattice Boltzmann Equation (FDLBE). Because source terms that reflect the interaction between phases can be accurately described, the new model is suitable for simulating two-way coupling incompressible multiphase flow The 2-D particle-laden flow over a backward-facing step is chosen as a test case to validate the present method. Favorable results are obtained and the present scheme is shown to have good prospects in practical applications.展开更多
A definition is introduced about traveling waves of 2-1 dimension lattice difference equations. Discrete heat equation is introduced and a discussion is given for the existence of traveling waves. The theory of travel...A definition is introduced about traveling waves of 2-1 dimension lattice difference equations. Discrete heat equation is introduced and a discussion is given for the existence of traveling waves. The theory of traveling waves is extended on 2-1 dimension lattice difference equations. As an application, an example is presented to illustrate the main results.展开更多
A new discrete isospectral problem is introduced,from which a hierarchy of Lax i ntegrable lattice equation is deduced. By using the trace identity,the correspon ding Hamiltonian structure is given and its Liouville i...A new discrete isospectral problem is introduced,from which a hierarchy of Lax i ntegrable lattice equation is deduced. By using the trace identity,the correspon ding Hamiltonian structure is given and its Liouville integrability is proved.展开更多
A lattice Boltzmann flux solver(LBFS)is presented for simulation of fluid flows.Like the conventional computational fluid dynamics(CFD)solvers,the new solver also applies the finite volume method to discretize the gov...A lattice Boltzmann flux solver(LBFS)is presented for simulation of fluid flows.Like the conventional computational fluid dynamics(CFD)solvers,the new solver also applies the finite volume method to discretize the governing differential equations,but the numerical flux at the cell interface is not evaluated by the smooth function approximation or Riemann solvers.Instead,it is evaluated from local solution of lattice Boltzmann equation(LBE)at cell interface.Two versions of LBFS are presented in this paper.One is to locally apply one-dimensional compressible lattice Boltzmann(LB)model along the normal direction to the cell interface for simulation of compressible inviscid flows with shock waves.The other is to locally apply multi-dimensional LB model at cell interface for simulation of incompressible viscous and inviscid flows.The present solver removes the drawbacks of conventional lattice Boltzmann method(LBM)such as limitation to uniform mesh,tie-up of mesh spacing and time interval,limitation to viscous flows.Numerical examples show that the present solver can be well applied to simulate fluid flows with non-uniform mesh and curved boundary.展开更多
A method combining the immersed boundary technique and a multi- relaxation-time (MRT) lattice Boltzmann flux solver (LBFS) is presented for numerical simulation of incompressible flows over circular and elliptic c...A method combining the immersed boundary technique and a multi- relaxation-time (MRT) lattice Boltzmann flux solver (LBFS) is presented for numerical simulation of incompressible flows over circular and elliptic cylinders and NACA 0012 Airfoil. The method uses a simple Cartesian mesh to simulate flows past immersed complicated bodies. With the Chapman-Enskog expansion analysis, a transform is performed between the Navier-Stokes and lattice Boltzmann equations (LBEs). The LBFS is used to discretize the macroscopic differential equations with a finite volume method and evaluate the interface fluxes through local reconstruction of the lattice Boltzmann solution. The immersed boundary technique is used to correct the intermediate velocity around the solid boundary to satisfy the no-slip boundary condition. Agreement of simulation results with the data found in the literature shows reliability of the proposed method in simulating laminar flows on a Cartesian mesh.展开更多
Interaction between turbulence and particles is investigated in a channel flow. The fluid motion is calculated using direct numerical simulation(DNS) with a lattice Boltzmann(LB) method, and particles are tracked in a...Interaction between turbulence and particles is investigated in a channel flow. The fluid motion is calculated using direct numerical simulation(DNS) with a lattice Boltzmann(LB) method, and particles are tracked in a Lagrangian framework through the action of force imposed by the fluid. The particle diameter is smaller than the Kolmogorov length scale, and the point force is used to represent the feedback force of particles on the turbulence. The effects of particles on the turbulence and skin friction coefficient are examined with different particle inertias and mass loadings. Inertial particles suppress intensities of the spanwise and wall-normal components of velocity, and the Reynolds shear stress. It is also found that, relative to the reference particle-free flow,the overall mean skin-friction coefficient is reduced by particles. Changes of near wall turbulent structures such as longer and more regular streamwise low-speed streaks and less ejections and sweeps are the manifestation of drag reduction.展开更多
This paper focuses on a two-dimensional bidirectional pedestrian flow model which involves the next-nearest-neighbor effect. The stability condition and the Korteweg-de Vries (KdV) equation are derived to describe t...This paper focuses on a two-dimensional bidirectional pedestrian flow model which involves the next-nearest-neighbor effect. The stability condition and the Korteweg-de Vries (KdV) equation are derived to describe the density wave of pedestrian congestion by linear stability and nonlinear analysis. Through theoretical analysis, the soliton solution is obtained.展开更多
Large eddy simulation (LES) using the Smagorinsky eddy viscosity model is added to the two-dimensional nine velocity components (D2Q9) lattice Boltzmann equation (LBE) with multi-relaxation-time (MRT) to simul...Large eddy simulation (LES) using the Smagorinsky eddy viscosity model is added to the two-dimensional nine velocity components (D2Q9) lattice Boltzmann equation (LBE) with multi-relaxation-time (MRT) to simulate incompressible turbulent cavity flows with the Reynolds numbers up to 1 × 10^7. To improve the computation efficiency of LBM on the numerical simulations of turbulent flows, the massively parallel computing power from a graphic processing unit (GPU) with a computing unified device architecture (CUDA) is introduced into the MRT-LBE-LES model. The model performs well, compared with the results from others, with an increase of 76 times in computation efficiency. It appears that the higher the Reynolds numbers is, the smaller the Smagorinsky constant should be, if the lattice number is fixed. Also, for a selected high Reynolds number and a selected proper Smagorinsky constant, there is a minimum requirement for the lattice number so that the Smagorinsky eddy viscosity will not be excessively large.展开更多
We present a new reliable analytical study for solving the discontinued problems arising in nanotechnology. Such problems are presented as nonlinear differential-difference equations. The proposed method is based on t...We present a new reliable analytical study for solving the discontinued problems arising in nanotechnology. Such problems are presented as nonlinear differential-difference equations. The proposed method is based on the Laplace trans- form with the homotopy analysis method (HAM). This method is a powerful tool for solving a large amount of problems. This technique provides a series of functions which may converge to the exact solution of the problem. A good agreement between the obtained solution and some well-known results is obtained.展开更多
Many multi-dimensional consistent discrete systems have soliton solutions with nonzero backgrounds, which brings difficulty in the investigation of integrable characteristics. In this paper, we derive infinitely many ...Many multi-dimensional consistent discrete systems have soliton solutions with nonzero backgrounds, which brings difficulty in the investigation of integrable characteristics. In this paper, we derive infinitely many conserved quantities for the lattice potential Korteweg-de Vries equation whose solutions have nonzero backgrounds. The derivation is based on the fact that the scattering data a(z) is independent of discrete space and time and the analytic property of Jost solutions of the discrete Schr5dinger spectral problem. The obtained conserved densities are asymptotic to zero when |n| (or |m|) tends to infinity. To obtain these results, we reconstruct a discrete Riccati equation by using a conformal map which transforms the upper complex plane to the inside of unit circle. Series solution to the Riccati equation is constructed based on the analytic and asymptotic properties of Jost solutions.展开更多
The optical flow analysis of the image sequence based on the formal lattice Boltzmann equation, with different DdQm models, is discussed in this paper. The Mgorithm is based on the lattice Boltzmann method (LBM), wh...The optical flow analysis of the image sequence based on the formal lattice Boltzmann equation, with different DdQm models, is discussed in this paper. The Mgorithm is based on the lattice Boltzmann method (LBM), which is used in computational fluid dynamics theory for the simulation of fluid dynamics. At first, a generalized approximation to the formal lattice Boltzmann equation is discussed. Then the effects of different DdQm models on the results of the optical flow estimation are compared with each other, while calculating the movement vectors of pixels in the image sequence. The experimental results show that the higher dimension DdQm models, e.g., D3Q15 are more effective than those lower dimension ones.展开更多
This paper presents a new finite-volume discretization of a generalised LatticeBoltzmann equation (LBE) on unstructured grids. This equation is the continuumLBE, with the addition of a second order time derivative ter...This paper presents a new finite-volume discretization of a generalised LatticeBoltzmann equation (LBE) on unstructured grids. This equation is the continuumLBE, with the addition of a second order time derivative term (memory), and is derivedfrom a second-order differential form of the semi-discrete Boltzmann equationin its implicit form. The new scheme, named unstructured lattice Boltzmann equationwith memory (ULBEM), can be advanced in time with a larger time-step than the previousunstructured LB formulations, and a theoretical demonstration of the improvedstability is provided. Taylor vortex simulations show that the viscosity is the same aswith standard ULBE and demonstrates that the new scheme improves both stabilityand accuracy. Model validation is also demonstrated by simulating backward-facingstep flow at low and moderate Reynolds numbers, as well as by comparing the reattachmentlength of the recirculating eddy behind the step against experimental andnumerical data available in literature.展开更多
A hierarchy of integrable lattice soliton equations and its Hamiitonian struc ture associated a 3×3 matrix spectral problem are got. An integrable symplectic map is obtained by nonlinearization of Lax pairs and a...A hierarchy of integrable lattice soliton equations and its Hamiitonian struc ture associated a 3×3 matrix spectral problem are got. An integrable symplectic map is obtained by nonlinearization of Lax pairs and ad joint Lax pairs of the hierarchy. Moreover, the solutions to the prototype system of lattice equations in the hierarchy are reduced to the solutions of a system of ordinary differential equations and a simple iterative process of the symplectic map.展开更多
A class of lattice with time delay and nonlocal response is considered.By transforming the lattice delay differential system into an integral equations in a Banach space,we reduces a singular perturbation problem to a...A class of lattice with time delay and nonlocal response is considered.By transforming the lattice delay differential system into an integral equations in a Banach space,we reduces a singular perturbation problem to a regular perturbation problem.Traveling wave solution therefore is obtained by applying Liapunov-Schmidt method and the implicit function theorem.展开更多
A lattice Boltzmann flux solver(LBFS)is presented in this work for simulation of incompressible viscous and inviscid flows.The new solver is based on Chapman-Enskog expansion analysis,which is the bridge to link Navie...A lattice Boltzmann flux solver(LBFS)is presented in this work for simulation of incompressible viscous and inviscid flows.The new solver is based on Chapman-Enskog expansion analysis,which is the bridge to link Navier-Stokes(N-S)equations and lattice Boltzmann equation(LBE).The macroscopic differential equations are discretized by the finite volume method,where the flux at the cell interface is evaluated by local reconstruction of lattice Boltzmann solution from macroscopic flow variables at cell centers.The new solver removes the drawbacks of conventional lattice Boltzmann method such as limitation to uniform mesh,tie-up of mesh spacing and time interval,limitation to viscous flows.LBFS is validated by its application to simulate the viscous decaying vortex flow,the driven cavity flow,the viscous flow past a circular cylinder,and the inviscid flow past a circular cylinder.The obtained numerical results compare very well with available data in the literature,which show that LBFS has the second order of accuracy in space,and can be well applied to viscous and inviscid flow problems with non-uniform mesh and curved boundary.展开更多
We develop a two-relaxation-time (TRT) Lattice Boltzmann model for hydrodynamicequations with variable source terms based on equivalent equilibriumfunctions. A special parametrization of the free relaxation parameter ...We develop a two-relaxation-time (TRT) Lattice Boltzmann model for hydrodynamicequations with variable source terms based on equivalent equilibriumfunctions. A special parametrization of the free relaxation parameter is derived. Itcontrols, in addition to the non-dimensional hydrodynamic numbers, any TRT macroscopicsteady solution and governs the spatial discretization of transient flows. Inthis framework, the multi-reflection approach [16, 18] is generalized and extended forDirichlet velocity, pressure and mixed (pressure/tangential velocity) boundary conditions.We propose second and third-order accurate boundary schemes and adapt themfor corners. The boundary schemes are analyzed for exactness of the parametrization,uniqueness of their steady solutions, support of staggered invariants and for the effectiveaccuracy in case of time dependent boundary conditions and transient flow.When the boundary scheme obeys the parametrization properly, the derived permeabilityvalues become independent of the selected viscosity for any porous structureand can be computed efficiently. The linear interpolations [5, 46] are improved withrespect to this property.展开更多
基金Supported by the Science Research Foundation of Zhanjiang Normal University(L0803)
文摘The Bcklund transformation and the generalized Miura transformation for the Volterra lattice equation are constructed by using point symmetry method.As an application,the explicit solution to the lattice equation is obtained.
基金Project supported by National Natural Science Fundation of China(Grant No .10371070)
文摘A discrete isospectral problem and the associated hierarchy of Lax integrable lattice equations were investigated. A Darboux transformation for the discrete spectral problem was found. Finally, an infinite number of conservation laws were given for the corresponding hierarchy.
文摘A direct way to construct integrable couplings for discrete systems is presented by use of two semi-direct sum Lie algebras. As their applications, the discrete integrable couplings associated with modified Korteweg-de Vries (m-KdV) lattice and two hierarchies of discrete soliton equations are developed. It is also indicated that the study of integrable couplings using semi-direct sums of Lie algebras is an important step towards the complete classification of integrable couplings.
基金Project supported by the National Natural Science Foundation of China (Grant No.10672143)
文摘This paper focuses on studying the symmetry of a practical wave equation on new lattices. It is a new step in that the new lattice equation is applied to reduce the discrete problem of motion of an elastic thin homogeneous bar. The equation of motion of the bar can be changed into a discrete wave equation. With the new lattice equation, the translational and scaling invariant, not only is the infinitesimal transformation given, but the symmetry and Lie algebras are also calculated. We also give a new form of invariant called the ratio invariant, which can reduce the process of the computing invariant with the characteristic equation.
基金The project supported by the National Natural Science Foundation of China(60073044)the State Key Development Programme for Basic Research of China(G1990022207).
文摘In this paper, we propose a novel incompressible finite-difference lattice Boltzmann Equation (FDLBE). Because source terms that reflect the interaction between phases can be accurately described, the new model is suitable for simulating two-way coupling incompressible multiphase flow The 2-D particle-laden flow over a backward-facing step is chosen as a test case to validate the present method. Favorable results are obtained and the present scheme is shown to have good prospects in practical applications.
基金Supported by the National Natural Science Foundation of China(Ill61049)
文摘A definition is introduced about traveling waves of 2-1 dimension lattice difference equations. Discrete heat equation is introduced and a discussion is given for the existence of traveling waves. The theory of traveling waves is extended on 2-1 dimension lattice difference equations. As an application, an example is presented to illustrate the main results.
文摘A new discrete isospectral problem is introduced,from which a hierarchy of Lax i ntegrable lattice equation is deduced. By using the trace identity,the correspon ding Hamiltonian structure is given and its Liouville integrability is proved.
基金Supported by the National Natural Science Foundation of China(11272153)
文摘A lattice Boltzmann flux solver(LBFS)is presented for simulation of fluid flows.Like the conventional computational fluid dynamics(CFD)solvers,the new solver also applies the finite volume method to discretize the governing differential equations,but the numerical flux at the cell interface is not evaluated by the smooth function approximation or Riemann solvers.Instead,it is evaluated from local solution of lattice Boltzmann equation(LBE)at cell interface.Two versions of LBFS are presented in this paper.One is to locally apply one-dimensional compressible lattice Boltzmann(LB)model along the normal direction to the cell interface for simulation of compressible inviscid flows with shock waves.The other is to locally apply multi-dimensional LB model at cell interface for simulation of incompressible viscous and inviscid flows.The present solver removes the drawbacks of conventional lattice Boltzmann method(LBM)such as limitation to uniform mesh,tie-up of mesh spacing and time interval,limitation to viscous flows.Numerical examples show that the present solver can be well applied to simulate fluid flows with non-uniform mesh and curved boundary.
文摘A method combining the immersed boundary technique and a multi- relaxation-time (MRT) lattice Boltzmann flux solver (LBFS) is presented for numerical simulation of incompressible flows over circular and elliptic cylinders and NACA 0012 Airfoil. The method uses a simple Cartesian mesh to simulate flows past immersed complicated bodies. With the Chapman-Enskog expansion analysis, a transform is performed between the Navier-Stokes and lattice Boltzmann equations (LBEs). The LBFS is used to discretize the macroscopic differential equations with a finite volume method and evaluate the interface fluxes through local reconstruction of the lattice Boltzmann solution. The immersed boundary technique is used to correct the intermediate velocity around the solid boundary to satisfy the no-slip boundary condition. Agreement of simulation results with the data found in the literature shows reliability of the proposed method in simulating laminar flows on a Cartesian mesh.
基金Project supported by the National Natural Science Foundation of China(Nos.11572183 and 11272198)
文摘Interaction between turbulence and particles is investigated in a channel flow. The fluid motion is calculated using direct numerical simulation(DNS) with a lattice Boltzmann(LB) method, and particles are tracked in a Lagrangian framework through the action of force imposed by the fluid. The particle diameter is smaller than the Kolmogorov length scale, and the point force is used to represent the feedback force of particles on the turbulence. The effects of particles on the turbulence and skin friction coefficient are examined with different particle inertias and mass loadings. Inertial particles suppress intensities of the spanwise and wall-normal components of velocity, and the Reynolds shear stress. It is also found that, relative to the reference particle-free flow,the overall mean skin-friction coefficient is reduced by particles. Changes of near wall turbulent structures such as longer and more regular streamwise low-speed streaks and less ejections and sweeps are the manifestation of drag reduction.
基金Project supported by the National Natural Science Foundation of China(Grant No.11072117)the Scientific Research Fund of Zhejiang Province,China(Grant No.LY13A010005)+4 种基金the Disciplinary Project of Ningbo City,China(Grant No.SZXL1067)the Scientific Research Fund of Education Department of Zhejiang Province,China(Grant No.Z201119278)the Natural Science Foundation of Ningbo City,China(Grant Nos.2012A610152 and 2012A610038)the K.C.Wong Magna Fund in Ningbo University,Chinathe Research Grant Council,Government of the Hong Kong Administrative Region,China(Grant No.CityU119011)
文摘This paper focuses on a two-dimensional bidirectional pedestrian flow model which involves the next-nearest-neighbor effect. The stability condition and the Korteweg-de Vries (KdV) equation are derived to describe the density wave of pedestrian congestion by linear stability and nonlinear analysis. Through theoretical analysis, the soliton solution is obtained.
基金supported by College of William and Mary,Virginia Institute of Marine Science for the study environment
文摘Large eddy simulation (LES) using the Smagorinsky eddy viscosity model is added to the two-dimensional nine velocity components (D2Q9) lattice Boltzmann equation (LBE) with multi-relaxation-time (MRT) to simulate incompressible turbulent cavity flows with the Reynolds numbers up to 1 × 10^7. To improve the computation efficiency of LBM on the numerical simulations of turbulent flows, the massively parallel computing power from a graphic processing unit (GPU) with a computing unified device architecture (CUDA) is introduced into the MRT-LBE-LES model. The model performs well, compared with the results from others, with an increase of 76 times in computation efficiency. It appears that the higher the Reynolds numbers is, the smaller the Smagorinsky constant should be, if the lattice number is fixed. Also, for a selected high Reynolds number and a selected proper Smagorinsky constant, there is a minimum requirement for the lattice number so that the Smagorinsky eddy viscosity will not be excessively large.
文摘We present a new reliable analytical study for solving the discontinued problems arising in nanotechnology. Such problems are presented as nonlinear differential-difference equations. The proposed method is based on the Laplace trans- form with the homotopy analysis method (HAM). This method is a powerful tool for solving a large amount of problems. This technique provides a series of functions which may converge to the exact solution of the problem. A good agreement between the obtained solution and some well-known results is obtained.
文摘Many multi-dimensional consistent discrete systems have soliton solutions with nonzero backgrounds, which brings difficulty in the investigation of integrable characteristics. In this paper, we derive infinitely many conserved quantities for the lattice potential Korteweg-de Vries equation whose solutions have nonzero backgrounds. The derivation is based on the fact that the scattering data a(z) is independent of discrete space and time and the analytic property of Jost solutions of the discrete Schr5dinger spectral problem. The obtained conserved densities are asymptotic to zero when |n| (or |m|) tends to infinity. To obtain these results, we reconstruct a discrete Riccati equation by using a conformal map which transforms the upper complex plane to the inside of unit circle. Series solution to the Riccati equation is constructed based on the analytic and asymptotic properties of Jost solutions.
基金Project supported by the National Natural Science Foundation of China(Grant No.40976108)the Shanghai Leading Academic Discipline Project(Grant No.J50103)the Innovation Program of Municipal Education Commission of Shanghai Municipality(Grant No.11YZ03)
文摘The optical flow analysis of the image sequence based on the formal lattice Boltzmann equation, with different DdQm models, is discussed in this paper. The Mgorithm is based on the lattice Boltzmann method (LBM), which is used in computational fluid dynamics theory for the simulation of fluid dynamics. At first, a generalized approximation to the formal lattice Boltzmann equation is discussed. Then the effects of different DdQm models on the results of the optical flow estimation are compared with each other, while calculating the movement vectors of pixels in the image sequence. The experimental results show that the higher dimension DdQm models, e.g., D3Q15 are more effective than those lower dimension ones.
文摘This paper presents a new finite-volume discretization of a generalised LatticeBoltzmann equation (LBE) on unstructured grids. This equation is the continuumLBE, with the addition of a second order time derivative term (memory), and is derivedfrom a second-order differential form of the semi-discrete Boltzmann equationin its implicit form. The new scheme, named unstructured lattice Boltzmann equationwith memory (ULBEM), can be advanced in time with a larger time-step than the previousunstructured LB formulations, and a theoretical demonstration of the improvedstability is provided. Taylor vortex simulations show that the viscosity is the same aswith standard ULBE and demonstrates that the new scheme improves both stabilityand accuracy. Model validation is also demonstrated by simulating backward-facingstep flow at low and moderate Reynolds numbers, as well as by comparing the reattachmentlength of the recirculating eddy behind the step against experimental andnumerical data available in literature.
文摘A hierarchy of integrable lattice soliton equations and its Hamiitonian struc ture associated a 3×3 matrix spectral problem are got. An integrable symplectic map is obtained by nonlinearization of Lax pairs and ad joint Lax pairs of the hierarchy. Moreover, the solutions to the prototype system of lattice equations in the hierarchy are reduced to the solutions of a system of ordinary differential equations and a simple iterative process of the symplectic map.
文摘A class of lattice with time delay and nonlocal response is considered.By transforming the lattice delay differential system into an integral equations in a Banach space,we reduces a singular perturbation problem to a regular perturbation problem.Traveling wave solution therefore is obtained by applying Liapunov-Schmidt method and the implicit function theorem.
文摘A lattice Boltzmann flux solver(LBFS)is presented in this work for simulation of incompressible viscous and inviscid flows.The new solver is based on Chapman-Enskog expansion analysis,which is the bridge to link Navier-Stokes(N-S)equations and lattice Boltzmann equation(LBE).The macroscopic differential equations are discretized by the finite volume method,where the flux at the cell interface is evaluated by local reconstruction of lattice Boltzmann solution from macroscopic flow variables at cell centers.The new solver removes the drawbacks of conventional lattice Boltzmann method such as limitation to uniform mesh,tie-up of mesh spacing and time interval,limitation to viscous flows.LBFS is validated by its application to simulate the viscous decaying vortex flow,the driven cavity flow,the viscous flow past a circular cylinder,and the inviscid flow past a circular cylinder.The obtained numerical results compare very well with available data in the literature,which show that LBFS has the second order of accuracy in space,and can be well applied to viscous and inviscid flow problems with non-uniform mesh and curved boundary.
文摘We develop a two-relaxation-time (TRT) Lattice Boltzmann model for hydrodynamicequations with variable source terms based on equivalent equilibriumfunctions. A special parametrization of the free relaxation parameter is derived. Itcontrols, in addition to the non-dimensional hydrodynamic numbers, any TRT macroscopicsteady solution and governs the spatial discretization of transient flows. Inthis framework, the multi-reflection approach [16, 18] is generalized and extended forDirichlet velocity, pressure and mixed (pressure/tangential velocity) boundary conditions.We propose second and third-order accurate boundary schemes and adapt themfor corners. The boundary schemes are analyzed for exactness of the parametrization,uniqueness of their steady solutions, support of staggered invariants and for the effectiveaccuracy in case of time dependent boundary conditions and transient flow.When the boundary scheme obeys the parametrization properly, the derived permeabilityvalues become independent of the selected viscosity for any porous structureand can be computed efficiently. The linear interpolations [5, 46] are improved withrespect to this property.