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.展开更多
This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the...This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the rigid structure is taken as ＂fictitious＂ fluid with zero strain rate. Both fluid and structure are described by velocity and pressure. The whole domain, including fluid region and structure region, is modeled by the incompressible Navier-Stokes equations which are discretized with fixed Eulerian mesh. However, to keep the structure＇ s rigid body shape and behavior, a rigid body constraint is enforced on the ＂fictitious＂ fluid domain by use of the Distributed Lagrange Multipher/Fictitious Domain （DLM/ FD） method which is originally introduced to solve particulate flow problems by Glowinski et al. For the verification of the model presented herein, a 2D numerical wave tank is established to simulate small amplitude wave propagations, and then numerical results are compared with analytical solutions. Finally, a 2D example of fluid-structure interaction under wave dynamic forces provides convincing evidences for the method excellent solution quality and fidelity.展开更多
This paper discusses a fictitious domain method for the linear Dirichlet problem and its applications to the generalized Stokes problem. This method treats Dirichlet boundary condit ion via a Lagrange multiplier tec...This paper discusses a fictitious domain method for the linear Dirichlet problem and its applications to the generalized Stokes problem. This method treats Dirichlet boundary condit ion via a Lagrange multiplier technique and is well suited to the no-slip bound ary condition in viscous flow problems. In order to improve the accuracy of solu tions, meshes are refined according to the a posteriori error estimate. The mini -element discretization is applied to solve the generalized Stokes problem. Fin ally, some numerical results to validate this method are presented for partial d ifferential equations with Dirichlet boundary condition.展开更多
Abstract. In this paper which is motivated by computation on parallel machine, we showthat the superconvergence results of the finite element method(FEM) with Lagrange mul-tipliers based on domain decomposition method...Abstract. In this paper which is motivated by computation on parallel machine, we showthat the superconvergence results of the finite element method(FEM) with Lagrange mul-tipliers based on domain decomposition method (DDM) with nonmatching grids can becarried over to parabolic problems. The main idea of this paper is to achieve the combina-tion of parallel computational method with the higher accuracy technique by interpolationfinite element postprocessing.展开更多
针对目前极化敏感面阵空域-极化域联合谱估计运算量大、耗时长的问题,提出一种降维求根MUSIC(Multiple Signal Classification)优化算法。通过对接收信号进行降维处理,提出新的求解模型将传统四维MUSIC转化为两个一维求根MUSIC求解空域...针对目前极化敏感面阵空域-极化域联合谱估计运算量大、耗时长的问题,提出一种降维求根MUSIC(Multiple Signal Classification)优化算法。通过对接收信号进行降维处理,提出新的求解模型将传统四维MUSIC转化为两个一维求根MUSIC求解空域波达方向和引用已求解出的空域信息结合拉格朗日乘子法解决来波信号极化信息估计问题。相比传统的4D-MUSIC和秩亏MUSIC,所提算法在不损失估计精度的前提下提高了运算速度,降低了运算复杂度,无需谱峰搜索过程,消除了因搜索步长而导致的量化误差。对日后大规模阵列计算及MIMO(Multiple Input Multiple Output)雷达引入提供快速求解方法。仿真实验表明,所提算法在低信噪比0 dB下空域误差约为0.85°,速度相比秩亏MUSIC提升了约64.7%,验证了该算法的有效性和高精度性。展开更多
The sedimentation of a single circular particle between two parallel walls was studied by means of direct numerical simulation (DNS) and experiment. The improved implementation of distributed Lagrange multiplier/ficti...The sedimentation of a single circular particle between two parallel walls was studied by means of direct numerical simulation (DNS) and experiment. The improved implementation of distributed Lagrange multiplier/fictitious domain method used in our DNS is a promising new way for simulation of particulate flows. The settling behaviors of the particle are presented ranging in Reynolds number from 0 to about 700, which showed that our results for low Reynolds numbers agreed well with that reported before. Nevertheless, for higher Reynolds numbers our results were different from theirs. The long-term mean equilibrium positions in our results were all on the centerline, but not at off-center position as reported before. In order to validate our simulation, experiments were also conducted. The results showed that the sedimenting behavior simulated in this paper agreed well with our experiment result.展开更多
An improved implementation of Distributed Lagrange multiplier/fictitious domain method was presented and used to simulate the interactions between two circular particles sedimenting in a two_dimensional channel. The s...An improved implementation of Distributed Lagrange multiplier/fictitious domain method was presented and used to simulate the interactions between two circular particles sedimenting in a two_dimensional channel. The simulation results were verified by comparison with experiments. The results show that the interactions between two particles with different sizes can be described as drafting, kissing, tumbling and separating. Only for small diameter ratio, the two particles will interact undergoing repeated DKT (Drafting, Kissing and Tumbling) process. Otherwise, the two particles will separate after their tumbling. The results also show that, during the interaction process, the motion of the small particle is strongly affected while the large particle is affected slightly.展开更多
.In this paper,an augmented Lagrangian Uzawa iterative method is developed and analyzed for solving a class of double saddle-point systems with semidefinite(2,2)block.Convergence of the iterativemethod is proved under....In this paper,an augmented Lagrangian Uzawa iterative method is developed and analyzed for solving a class of double saddle-point systems with semidefinite(2,2)block.Convergence of the iterativemethod is proved under the assumption that the double saddle-point problem exists a unique solution.An application of the iterative method to the double saddle-point systems arising from the distributed Lagrange multiplier/fictitious domain(DLM/FD)finite element method for solving elliptic interface problems is also presented,in which the existence and uniqueness of the double saddle-point system is guaranteed by the analysis of the DLM/FD finite element method.Numerical experiments are conducted to validate the theoretical results and to study the performance of the proposed iterative method.展开更多
An improved implementation of distributed multiplier/fictitious domain method is presented for the direct numerical simulation of particulate flow. The key improvement is to replace a finite element triangulation for...An improved implementation of distributed multiplier/fictitious domain method is presented for the direct numerical simulation of particulate flow. The key improvement is to replace a finite element triangulation for the velocity and a “twice coarser' triangulation for the pressure with a rectangular discretization for the velocity and pressure. For code validation, the sedimentation of a single particle in a two dimensional channel was simulated. The results showed that the simulation is independent of the mesh size as well as the time step. The comparison between experimental data and this simulation showed that our code can give a more accurate simulation on the motion of particles than previous DLM code. The code was then applied to simulate the sedimentation of 600 particles in a rectangular box. The falling course is presented and discussed. At the same time, this simulation also demonstrates that the method presented in this paper can be used for solving the initial problems involving a lager number of particles exactly with computing durations kept at acceptable levels.展开更多
In this paper, a new method-Distributed La-grange Multiplier/fictitiousdomain (DLM) method for partic-ulate flows was improved. A rectangular mesh was introduced todiscretize the domain, and the buoyant force was cons...In this paper, a new method-Distributed La-grange Multiplier/fictitiousdomain (DLM) method for partic-ulate flows was improved. A rectangular mesh was introduced todiscretize the domain, and the buoyant force was considered for predicting the positions ofparticles. In order to validate the presented algorithm, the sedimentation of a single circularparticle was simulated using different mesh sizes and time steps firstly. The results show that thesimulation is independent of the mesh size as well as the time step. And then, the interactionsbetween two falling particles, including drafting, kissing and tumbling, and the sedimentation of 18particles also have been simulated with the code.展开更多
基金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.
基金This study is supported by the National Natural Science Foundation of China (Grant No50579046) the Science Foundation of Tianjin Municipal Commission of Science and Technology (Grant No043114711)
文摘This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the rigid structure is taken as ＂fictitious＂ fluid with zero strain rate. Both fluid and structure are described by velocity and pressure. The whole domain, including fluid region and structure region, is modeled by the incompressible Navier-Stokes equations which are discretized with fixed Eulerian mesh. However, to keep the structure＇ s rigid body shape and behavior, a rigid body constraint is enforced on the ＂fictitious＂ fluid domain by use of the Distributed Lagrange Multipher/Fictitious Domain （DLM/ FD） method which is originally introduced to solve particulate flow problems by Glowinski et al. For the verification of the model presented herein, a 2D numerical wave tank is established to simulate small amplitude wave propagations, and then numerical results are compared with analytical solutions. Finally, a 2D example of fluid-structure interaction under wave dynamic forces provides convincing evidences for the method excellent solution quality and fidelity.
文摘This paper discusses a fictitious domain method for the linear Dirichlet problem and its applications to the generalized Stokes problem. This method treats Dirichlet boundary condit ion via a Lagrange multiplier technique and is well suited to the no-slip bound ary condition in viscous flow problems. In order to improve the accuracy of solu tions, meshes are refined according to the a posteriori error estimate. The mini -element discretization is applied to solve the generalized Stokes problem. Fin ally, some numerical results to validate this method are presented for partial d ifferential equations with Dirichlet boundary condition.
文摘Abstract. In this paper which is motivated by computation on parallel machine, we showthat the superconvergence results of the finite element method(FEM) with Lagrange mul-tipliers based on domain decomposition method (DDM) with nonmatching grids can becarried over to parabolic problems. The main idea of this paper is to achieve the combina-tion of parallel computational method with the higher accuracy technique by interpolationfinite element postprocessing.
文摘The sedimentation of a single circular particle between two parallel walls was studied by means of direct numerical simulation (DNS) and experiment. The improved implementation of distributed Lagrange multiplier/fictitious domain method used in our DNS is a promising new way for simulation of particulate flows. The settling behaviors of the particle are presented ranging in Reynolds number from 0 to about 700, which showed that our results for low Reynolds numbers agreed well with that reported before. Nevertheless, for higher Reynolds numbers our results were different from theirs. The long-term mean equilibrium positions in our results were all on the centerline, but not at off-center position as reported before. In order to validate our simulation, experiments were also conducted. The results showed that the sedimenting behavior simulated in this paper agreed well with our experiment result.
文摘An improved implementation of Distributed Lagrange multiplier/fictitious domain method was presented and used to simulate the interactions between two circular particles sedimenting in a two_dimensional channel. The simulation results were verified by comparison with experiments. The results show that the interactions between two particles with different sizes can be described as drafting, kissing, tumbling and separating. Only for small diameter ratio, the two particles will interact undergoing repeated DKT (Drafting, Kissing and Tumbling) process. Otherwise, the two particles will separate after their tumbling. The results also show that, during the interaction process, the motion of the small particle is strongly affected while the large particle is affected slightly.
基金supported by the 10 plus 10 project of Tongji University(No.4260141304/004/010).
文摘.In this paper,an augmented Lagrangian Uzawa iterative method is developed and analyzed for solving a class of double saddle-point systems with semidefinite(2,2)block.Convergence of the iterativemethod is proved under the assumption that the double saddle-point problem exists a unique solution.An application of the iterative method to the double saddle-point systems arising from the distributed Lagrange multiplier/fictitious domain(DLM/FD)finite element method for solving elliptic interface problems is also presented,in which the existence and uniqueness of the double saddle-point system is guaranteed by the analysis of the DLM/FD finite element method.Numerical experiments are conducted to validate the theoretical results and to study the performance of the proposed iterative method.
基金TheNationalNaturalSciencesFoundationforOutstandingYouthofChina (No .19925210)andZhejiangProvincialNaturalScienceFoundationofChina(No .10 10 4 7)
文摘An improved implementation of distributed multiplier/fictitious domain method is presented for the direct numerical simulation of particulate flow. The key improvement is to replace a finite element triangulation for the velocity and a “twice coarser' triangulation for the pressure with a rectangular discretization for the velocity and pressure. For code validation, the sedimentation of a single particle in a two dimensional channel was simulated. The results showed that the simulation is independent of the mesh size as well as the time step. The comparison between experimental data and this simulation showed that our code can give a more accurate simulation on the motion of particles than previous DLM code. The code was then applied to simulate the sedimentation of 600 particles in a rectangular box. The falling course is presented and discussed. At the same time, this simulation also demonstrates that the method presented in this paper can be used for solving the initial problems involving a lager number of particles exactly with computing durations kept at acceptable levels.
文摘In this paper, a new method-Distributed La-grange Multiplier/fictitiousdomain (DLM) method for partic-ulate flows was improved. A rectangular mesh was introduced todiscretize the domain, and the buoyant force was considered for predicting the positions ofparticles. In order to validate the presented algorithm, the sedimentation of a single circularparticle was simulated using different mesh sizes and time steps firstly. The results show that thesimulation is independent of the mesh size as well as the time step. And then, the interactionsbetween two falling particles, including drafting, kissing and tumbling, and the sedimentation of 18particles also have been simulated with the code.
