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.展开更多
In this paper, an overlapping lattice Boltzmann model is introduced and its domain decomposition method, a distributed lattice Boltzmann method is presented. Parallel effectiveness of some programs based on the dist...In this paper, an overlapping lattice Boltzmann model is introduced and its domain decomposition method, a distributed lattice Boltzmann method is presented. Parallel effectiveness of some programs based on the distributed lattice Boltzmann method are analyzed.展开更多
For reducing the core loss of grain oriented silicon steel and improving its aging property, a new method, the LLSA by using Sb as the laser surface alloying element, was investigated, and at proper technique conditio...For reducing the core loss of grain oriented silicon steel and improving its aging property, a new method, the LLSA by using Sb as the laser surface alloying element, was investigated, and at proper technique conditions rather good result was obtained.展开更多
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.展开更多
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.展开更多
Based on the properties of wideband source’s spatial spread,an algorithm for reverberation suppression in beams domain is proposed.In an effective observation range,a compact bank of receiving spatial beams is operat...Based on the properties of wideband source’s spatial spread,an algorithm for reverberation suppression in beams domain is proposed.In an effective observation range,a compact bank of receiving spatial beams is operated to transform the data from element domain to beam domain.Angular Distribution Index(ADI) is obtained by comparing to a reference value which is computed by Monte Carlo Integration.ADI decreases with a stronger target echo.Hence,choosing an appropriate threshold of ADI,the target echo can be separated from reverberation by cutting off the high AD components.The filter is compatible with the procedures of beamforming and detection for the current sonar.The computation complexity of the filter is low.The relationships between ADI and frequency,bandwidth,signal-reverberation ratio,are illustrated by a simulation example.The processing results based on the lake-trial and sea-trial data show that reverberation can be significantly suppressed,and threshold choosing for wideband ADI filter is fault tolerant.展开更多
.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.展开更多
Distributed Lagrange Multiplier/Fictitious Domain (DLM/FD) method was employed to simulate the channel flow with polymer suspension. The polymer molecules were modeled as Finitely Extensible Nonlinear Elastic (FENE...Distributed Lagrange Multiplier/Fictitious Domain (DLM/FD) method was employed to simulate the channel flow with polymer suspension. The polymer molecules were modeled as Finitely Extensible Nonlinear Elastic (FENE) chains and FENE chain cluster. The coiled stretched transition of FENE chains was examined and the change in configuration of FENE chains was presented. The average velocity profile of the fully developed channel flow with 64 FENE chains was given and fitted well with the power-law curve. The change of chain cluster configuration was also simulated. These simulations show that DLM/FD method is capable of simulating the motion of not only FENE chain, but also FENE chain cluster.展开更多
基金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.
文摘In this paper, an overlapping lattice Boltzmann model is introduced and its domain decomposition method, a distributed lattice Boltzmann method is presented. Parallel effectiveness of some programs based on the distributed lattice Boltzmann method are analyzed.
基金National Natural Science FOundation of China! (No. 59974010).
文摘For reducing the core loss of grain oriented silicon steel and improving its aging property, a new method, the LLSA by using Sb as the laser surface alloying element, was investigated, and at proper technique conditions rather good result was obtained.
文摘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.
基金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.
基金supported by the National Natural Science Foundation of China(61302169)
文摘Based on the properties of wideband source’s spatial spread,an algorithm for reverberation suppression in beams domain is proposed.In an effective observation range,a compact bank of receiving spatial beams is operated to transform the data from element domain to beam domain.Angular Distribution Index(ADI) is obtained by comparing to a reference value which is computed by Monte Carlo Integration.ADI decreases with a stronger target echo.Hence,choosing an appropriate threshold of ADI,the target echo can be separated from reverberation by cutting off the high AD components.The filter is compatible with the procedures of beamforming and detection for the current sonar.The computation complexity of the filter is low.The relationships between ADI and frequency,bandwidth,signal-reverberation ratio,are illustrated by a simulation example.The processing results based on the lake-trial and sea-trial data show that reverberation can be significantly suppressed,and threshold choosing for wideband ADI filter is fault tolerant.
基金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.
文摘Distributed Lagrange Multiplier/Fictitious Domain (DLM/FD) method was employed to simulate the channel flow with polymer suspension. The polymer molecules were modeled as Finitely Extensible Nonlinear Elastic (FENE) chains and FENE chain cluster. The coiled stretched transition of FENE chains was examined and the change in configuration of FENE chains was presented. The average velocity profile of the fully developed channel flow with 64 FENE chains was given and fitted well with the power-law curve. The change of chain cluster configuration was also simulated. These simulations show that DLM/FD method is capable of simulating the motion of not only FENE chain, but also FENE chain cluster.