The interfacial effects on flow and heat transfer on micro/nano scale are discussed in this paper. Dif- ferent from bulk cases where interfaces can be simply treated as a boundary, the interfacial effects are not limi...The interfacial effects on flow and heat transfer on micro/nano scale are discussed in this paper. Dif- ferent from bulk cases where interfaces can be simply treated as a boundary, the interfacial effects are not limited to the interface on a microscale but could extend into a significant, even the whole domain of the flow and heat transfer field when the characteristic size of the domain is close to the mean free path (MFP) of the carriers inside an object. Most of microscale thermal phenomena result from interfa- cial interactions. Any changes in the interactions between the object and boundary particles, such as the force between fluid and solid wall particles, microstructure of interfaces, could affect thermal properties, flow and heat transfer characteristics and hence change thermal conductivity, velocity and temperature profiles, friction coefficient and thermal radiative properties, etc. The properties of nano- structure or flow and heat transfer features of fluid in micro/nanostructures not only depend on them- selves, but also on the interaction with the interface because the interface impact can go deep inside the flow. The same fluid, same channel geometry but different wall materials could have different flow and heat transport characteristics on microscale.展开更多
In this paper,a class of new immersed interface finite element methods (IIFEM) is developed to solve elasticity interface problems with homogeneous and non-homogeneous jump conditions in two dimensions.Simple non-body...In this paper,a class of new immersed interface finite element methods (IIFEM) is developed to solve elasticity interface problems with homogeneous and non-homogeneous jump conditions in two dimensions.Simple non-body-fitted meshes are used.For homogeneous jump conditions,both non-conforming and conforming basis functions are constructed in such a way that they satisfy the natural jump conditions. For non-homogeneous jump conditions,a pair of functions that satisfy the same non-homogeneous jump conditions are constructed using a level-set representation of the interface.With such a pair of functions,the discontinuities across the interface in the solution and flux are removed;and an equivalent elasticity interface problem with homogeneous jump conditions is formulated.Numerical examples are presented to demonstrate that such methods have second order convergence.展开更多
The distributed Lagrange multiplier/fictitious domain(DLM/FD)-mixed finite element method is developed and analyzed in this paper for a transient Stokes interface problem with jump coefficients.The semi-and fully disc...The distributed Lagrange multiplier/fictitious domain(DLM/FD)-mixed finite element method is developed and analyzed in this paper for a transient Stokes interface problem with jump coefficients.The semi-and fully discrete DLM/FD-mixed finite element scheme are developed for the first time for this problem with a moving interface,where the arbitrary Lagrangian-Eulerian(ALE)technique is employed to deal with the moving and immersed subdomain.Stability and optimal convergence properties are obtained for both schemes.Numerical experiments are carried out for different scenarios of jump coefficients,and all theoretical results are validated.展开更多
In this paper,we propose a numerical method for solving parabolic interface problems with nonhomogeneous flux jump condition and nonlinear jump condition.The main idea is to use traditional finite element method on se...In this paper,we propose a numerical method for solving parabolic interface problems with nonhomogeneous flux jump condition and nonlinear jump condition.The main idea is to use traditional finite element method on semi-Cartesian mesh coupled with Newton’s method to handle nonlinearity.It is easy to implement even though variable coefficients are used in the jump condition instead of constant in previous work for elliptic interface problem.Numerical experiments show that our method is about second order accurate in the L1 norm.展开更多
In this paper we present a one dimensional second order accurate method to solve Elliptic equations with discontinuous coefficients on an arbitrary interface.Second order accuracy for the first derivative is obtained ...In this paper we present a one dimensional second order accurate method to solve Elliptic equations with discontinuous coefficients on an arbitrary interface.Second order accuracy for the first derivative is obtained as well.The method is based on the Ghost Fluid Method,making use of ghost points on which the value is defined by suitable interface conditions.The multi-domain formulation is adopted,where the problem is split in two sub-problems and interface conditions will be enforced to close the problem.Interface conditions are relaxed together with the internal equations(following the approach proposed in[10]in the case of smooth coefficients),leading to an iterative method on all the set of grid values(inside points and ghost points).A multigrid approach with a suitable definition of the restriction operator is provided.The restriction of the defect is performed separately for both sub-problems,providing a convergence factor close to the one measured in the case of smooth coefficient and independent on the magnitude of the jump in the coefficient.Numerical tests will confirm the second order accuracy.Although the method is proposed in one dimension,the extension in higher dimension is currently underway[12]and it will be carried out by combining the discretization of[10]with the multigrid approach of[11]for Elliptic problems with non-eliminated boundary conditions in arbitrary domain.展开更多
In framework of the fictitious domain methods with immersed interfaces for the elasticity problem,the present contribution is to study and numerically validate the jump-integrated boundary conditions method with sharp...In framework of the fictitious domain methods with immersed interfaces for the elasticity problem,the present contribution is to study and numerically validate the jump-integrated boundary conditions method with sharp interface for the vector elasticity system discretized by a proposed finite volume method.The main idea of the fictitious domain approach consists in embedding the original domain of study into a geometrically larger and simpler one called the fictitious domain.Here,we present a cell-centered finite volume method to discretize the fictitious domain problem.The proposed method is numerically validated for different test cases.This work can be considered as a first step before more challenging problems such as fluid-structure interactions or moving interface problems.展开更多
We develop the immersed interface method(IIM)to simulate a two-fluid flow of two immiscible fluids with different density and viscosity.Due to the surface tension and the discontinuous fluid properties,the two-fluid f...We develop the immersed interface method(IIM)to simulate a two-fluid flow of two immiscible fluids with different density and viscosity.Due to the surface tension and the discontinuous fluid properties,the two-fluid flow has nonsmooth velocity and discontinuous pressure across the moving sharp interface separating the two fluids.The IIM computes the flow on a fixed Cartesian grid by incorporating into numerical schemes the necessary jump conditions induced by the interface.We present how to compute these necessary jump conditions from the analytical principal jump conditions derived in[Xu,DCDS,Supplement 2009,pp.838-845].We test our method on some canonical two-fluid flows.The results demonstrate that the method can handle large density and viscosity ratios,is second-order accurate in the infinity norm,and conserves mass inside a closed interface.展开更多
Based on an error estimate in terms of element edge vectors on arbitrary unstructured simplex meshes,we propose a new edge-based anisotropic mesh refinement algorithm.As the mesh adaptation indicator,the error estimat...Based on an error estimate in terms of element edge vectors on arbitrary unstructured simplex meshes,we propose a new edge-based anisotropic mesh refinement algorithm.As the mesh adaptation indicator,the error estimate involves only the gradient of error rather than higher order derivatives.The preferred refinement edge is chosen to reduce the maximal term in the error estimate.The algorithm is implemented in both two-and three-dimensional cases,and applied to the singular function interpolation and the elliptic interface problem.The numerical results demonstrate that the convergence order obtained by using the proposed anisotropic mesh refinement algorithm can be higher than that given by the isotropic one.展开更多
An iterative solver based on the immersed interface method is proposed to solve the pressure in a two-fluid flow on a Cartesian grid with second-order accuracy in the infinity norm.The iteration is constructed by intr...An iterative solver based on the immersed interface method is proposed to solve the pressure in a two-fluid flow on a Cartesian grid with second-order accuracy in the infinity norm.The iteration is constructed by introducing an unsteady term in the pressure Poisson equation.In each iteration step,a Helmholtz equation is solved on the Cartesian grid using FFT.The combination of the iteration and the immersed interface method enables the solver to handle various jump conditions across twofluid interfaces.This solver can also be used to solve Poisson equations on irregular domains.展开更多
Based on the lattice Boltzmann method,a lattice Boltzmann(LB) model of the ski-jump jet two-phase flow is developed first and the corresponding boundary conditions are studied.A simple case study of a droplet horizont...Based on the lattice Boltzmann method,a lattice Boltzmann(LB) model of the ski-jump jet two-phase flow is developed first and the corresponding boundary conditions are studied.A simple case study of a droplet horizontal movement calculation is carried out to test and verify the model,where level set method is used to track and reconstruct the moving droplet free surface. Then,we numerically simulate a two dimensional flow field of the ski-jump jet with the LB model,derive the moving surface and velocity vector field of the jet flow.The simulation results are very consistent with the physical mechanisms.The effectiveness and reliability of the model are demonstrated by the numerical examples.展开更多
The aim of the present work is to investigate the numerical modeling of interfacial cracks that may appear at the interface between two isotropic elastic materials. The extended finite element method is employed to an...The aim of the present work is to investigate the numerical modeling of interfacial cracks that may appear at the interface between two isotropic elastic materials. The extended finite element method is employed to analyze brittle and bi-material interfacial fatigue crack growth by computing the mixed mode stress intensity factors(SIF). Three different approaches are introduced to compute the SIFs. In the first one, mixed mode SIF is deduced from the computation of the contour integral as per the classical J-integral method,whereas a displacement method is used to evaluate the SIF by using either one or two displacement jumps located along the crack path in the second and third approaches. The displacement jump method is rather classical for mono-materials,but has to our knowledge not been used up to now for a bimaterial. Hence, use of displacement jump for characterizing bi-material cracks constitutes the main contribution of the present study. Several benchmark tests including parametric studies are performed to show the effectiveness of these computational methodologies for SIF considering static and fatigue problems of bi-material structures. It is found that results based on the displacement jump methods are in a very good agreement with those of exact solutions, such as for the J-integral method, but with a larger domain of applicability and a better numerical efficiency(less time consuming and less spurious boundary effect).展开更多
基金Supported by the National Natural Science Foundation of China (Grant No. 50376025)
文摘The interfacial effects on flow and heat transfer on micro/nano scale are discussed in this paper. Dif- ferent from bulk cases where interfaces can be simply treated as a boundary, the interfacial effects are not limited to the interface on a microscale but could extend into a significant, even the whole domain of the flow and heat transfer field when the characteristic size of the domain is close to the mean free path (MFP) of the carriers inside an object. Most of microscale thermal phenomena result from interfa- cial interactions. Any changes in the interactions between the object and boundary particles, such as the force between fluid and solid wall particles, microstructure of interfaces, could affect thermal properties, flow and heat transfer characteristics and hence change thermal conductivity, velocity and temperature profiles, friction coefficient and thermal radiative properties, etc. The properties of nano- structure or flow and heat transfer features of fluid in micro/nanostructures not only depend on them- selves, but also on the interaction with the interface because the interface impact can go deep inside the flow. The same fluid, same channel geometry but different wall materials could have different flow and heat transport characteristics on microscale.
基金supported by the US ARO grants 49308-MA and 56349-MAthe US AFSOR grant FA9550-06-1-024+1 种基金he US NSF grant DMS-0911434the State Key Laboratory of Scientific and Engineering Computing of Chinese Academy of Sciences during a visit by Z.Li between July-August,2008.
文摘In this paper,a class of new immersed interface finite element methods (IIFEM) is developed to solve elasticity interface problems with homogeneous and non-homogeneous jump conditions in two dimensions.Simple non-body-fitted meshes are used.For homogeneous jump conditions,both non-conforming and conforming basis functions are constructed in such a way that they satisfy the natural jump conditions. For non-homogeneous jump conditions,a pair of functions that satisfy the same non-homogeneous jump conditions are constructed using a level-set representation of the interface.With such a pair of functions,the discontinuities across the interface in the solution and flux are removed;and an equivalent elasticity interface problem with homogeneous jump conditions is formulated.Numerical examples are presented to demonstrate that such methods have second order convergence.
基金P.Sun was supported by NSF Grant DMS-1418806C.S.Zhang was partially supported by the National Key Research and Development Program of China(Grant No.2016YFB0201304)+1 种基金the Major Research Plan of National Natural Science Foundation of China(Grant Nos.91430215,91530323)the Key Research Program of Frontier Sciences of CAS.
文摘The distributed Lagrange multiplier/fictitious domain(DLM/FD)-mixed finite element method is developed and analyzed in this paper for a transient Stokes interface problem with jump coefficients.The semi-and fully discrete DLM/FD-mixed finite element scheme are developed for the first time for this problem with a moving interface,where the arbitrary Lagrangian-Eulerian(ALE)technique is employed to deal with the moving and immersed subdomain.Stability and optimal convergence properties are obtained for both schemes.Numerical experiments are carried out for different scenarios of jump coefficients,and all theoretical results are validated.
基金L.Shi’s research is supported by National Natural Science Foundation of China(No.11701569)S.Hou’s research is supported by Dr.Walter Koss Professorship made available through Louisiana Board of RegentsL.Wang’s research is supported by Science Foundations of China University of Petroleum-Beijing(No.2462015BJB05).
文摘In this paper,we propose a numerical method for solving parabolic interface problems with nonhomogeneous flux jump condition and nonlinear jump condition.The main idea is to use traditional finite element method on semi-Cartesian mesh coupled with Newton’s method to handle nonlinearity.It is easy to implement even though variable coefficients are used in the jump condition instead of constant in previous work for elliptic interface problem.Numerical experiments show that our method is about second order accurate in the L1 norm.
文摘In this paper we present a one dimensional second order accurate method to solve Elliptic equations with discontinuous coefficients on an arbitrary interface.Second order accuracy for the first derivative is obtained as well.The method is based on the Ghost Fluid Method,making use of ghost points on which the value is defined by suitable interface conditions.The multi-domain formulation is adopted,where the problem is split in two sub-problems and interface conditions will be enforced to close the problem.Interface conditions are relaxed together with the internal equations(following the approach proposed in[10]in the case of smooth coefficients),leading to an iterative method on all the set of grid values(inside points and ghost points).A multigrid approach with a suitable definition of the restriction operator is provided.The restriction of the defect is performed separately for both sub-problems,providing a convergence factor close to the one measured in the case of smooth coefficient and independent on the magnitude of the jump in the coefficient.Numerical tests will confirm the second order accuracy.Although the method is proposed in one dimension,the extension in higher dimension is currently underway[12]and it will be carried out by combining the discretization of[10]with the multigrid approach of[11]for Elliptic problems with non-eliminated boundary conditions in arbitrary domain.
文摘In framework of the fictitious domain methods with immersed interfaces for the elasticity problem,the present contribution is to study and numerically validate the jump-integrated boundary conditions method with sharp interface for the vector elasticity system discretized by a proposed finite volume method.The main idea of the fictitious domain approach consists in embedding the original domain of study into a geometrically larger and simpler one called the fictitious domain.Here,we present a cell-centered finite volume method to discretize the fictitious domain problem.The proposed method is numerically validated for different test cases.This work can be considered as a first step before more challenging problems such as fluid-structure interactions or moving interface problems.
基金the support of this work by the NSF grant DMS 0915237.
文摘We develop the immersed interface method(IIM)to simulate a two-fluid flow of two immiscible fluids with different density and viscosity.Due to the surface tension and the discontinuous fluid properties,the two-fluid flow has nonsmooth velocity and discontinuous pressure across the moving sharp interface separating the two fluids.The IIM computes the flow on a fixed Cartesian grid by incorporating into numerical schemes the necessary jump conditions induced by the interface.We present how to compute these necessary jump conditions from the analytical principal jump conditions derived in[Xu,DCDS,Supplement 2009,pp.838-845].We test our method on some canonical two-fluid flows.The results demonstrate that the method can handle large density and viscosity ratios,is second-order accurate in the infinity norm,and conserves mass inside a closed interface.
基金supported by the National Basic Research Program under the Grant 2005CB321701the National Science Foundation of China under the grant 10771008 and 10771211partial supported by A Foundation for the Author of National Excellent Doctoral Dissertation of PRC.
文摘Based on an error estimate in terms of element edge vectors on arbitrary unstructured simplex meshes,we propose a new edge-based anisotropic mesh refinement algorithm.As the mesh adaptation indicator,the error estimate involves only the gradient of error rather than higher order derivatives.The preferred refinement edge is chosen to reduce the maximal term in the error estimate.The algorithm is implemented in both two-and three-dimensional cases,and applied to the singular function interpolation and the elliptic interface problem.The numerical results demonstrate that the convergence order obtained by using the proposed anisotropic mesh refinement algorithm can be higher than that given by the isotropic one.
基金the support of this work by the NSF grant DMS 0915237.
文摘An iterative solver based on the immersed interface method is proposed to solve the pressure in a two-fluid flow on a Cartesian grid with second-order accuracy in the infinity norm.The iteration is constructed by introducing an unsteady term in the pressure Poisson equation.In each iteration step,a Helmholtz equation is solved on the Cartesian grid using FFT.The combination of the iteration and the immersed interface method enables the solver to handle various jump conditions across twofluid interfaces.This solver can also be used to solve Poisson equations on irregular domains.
基金supported by the National Natural Science Foundation of China(Grant No.50579083)
文摘Based on the lattice Boltzmann method,a lattice Boltzmann(LB) model of the ski-jump jet two-phase flow is developed first and the corresponding boundary conditions are studied.A simple case study of a droplet horizontal movement calculation is carried out to test and verify the model,where level set method is used to track and reconstruct the moving droplet free surface. Then,we numerically simulate a two dimensional flow field of the ski-jump jet with the LB model,derive the moving surface and velocity vector field of the jet flow.The simulation results are very consistent with the physical mechanisms.The effectiveness and reliability of the model are demonstrated by the numerical examples.
文摘The aim of the present work is to investigate the numerical modeling of interfacial cracks that may appear at the interface between two isotropic elastic materials. The extended finite element method is employed to analyze brittle and bi-material interfacial fatigue crack growth by computing the mixed mode stress intensity factors(SIF). Three different approaches are introduced to compute the SIFs. In the first one, mixed mode SIF is deduced from the computation of the contour integral as per the classical J-integral method,whereas a displacement method is used to evaluate the SIF by using either one or two displacement jumps located along the crack path in the second and third approaches. The displacement jump method is rather classical for mono-materials,but has to our knowledge not been used up to now for a bimaterial. Hence, use of displacement jump for characterizing bi-material cracks constitutes the main contribution of the present study. Several benchmark tests including parametric studies are performed to show the effectiveness of these computational methodologies for SIF considering static and fatigue problems of bi-material structures. It is found that results based on the displacement jump methods are in a very good agreement with those of exact solutions, such as for the J-integral method, but with a larger domain of applicability and a better numerical efficiency(less time consuming and less spurious boundary effect).
