To solve seismic wave scattering by a large-scale three-dimensional(3-D) hill topography, a fast parallel indirect boundary element method(IBEM) is developed by proposing a new construction method for the wave field, ...To solve seismic wave scattering by a large-scale three-dimensional(3-D) hill topography, a fast parallel indirect boundary element method(IBEM) is developed by proposing a new construction method for the wave field, modifying the generalized minimum residual(GMRES) algorithm and constructing an Open MP plus MPI parallel model. The validations of accuracy and efficiency show that this method can solve 3-D seismic response of a large-scale hill topography for broadband waves, and overcome the weakness of large storage and low efficiency of the traditional IBEM. Based on this new algorithm architecture, taking the broadband scattering of plane SV waves by a large-scale Gaussian-shaped hill of thousands-meters height as an example, the influence of several important parameters is investigated, including the incident frequency, the incident angle and the height-width and length-width ratio of the hill. The numerical results illustrate that the amplification effect on the ground motion by a near-hemispherical hill is more significant than the narrow hill. For low-frequency waves, the scattering effect of the higher hill is more pronounced, and there is only a single peak near the top of the hill. However, for high-frequency waves, rapid spatial variation of displacement amplitude appears on the hill surface.展开更多
The development of two boundary element algorithms for solving 3-D, frictional, and linear elastostatic contact problems is reported in this paper. The algorithms employ nonconforming discreti- zations for solving 3...The development of two boundary element algorithms for solving 3-D, frictional, and linear elastostatic contact problems is reported in this paper. The algorithms employ nonconforming discreti- zations for solving 3-D boundary element models, which provide much needed flexibility in the bound- ary element modeling for 3-D contact problems. These algorithms are implemented in a new 3-D boundary element code and verified using several examples. For the numerical examples studied, the results using the new boundary element algorithms match very well with the results using a commercial finite element code, and clearly demonstrate the feasibility of the new boundary element approach for 3-D contact analysis.展开更多
THE modeling of 2-D terrain effect on magnetotelluric(MT)has been solved using finite ele-ment method(FEM)and boundary element method(BEM).This note uses BEM tosolve the problem.The solution of H_x- and E_x-polarizati...THE modeling of 2-D terrain effect on magnetotelluric(MT)has been solved using finite ele-ment method(FEM)and boundary element method(BEM).This note uses BEM tosolve the problem.The solution of H_x- and E_x-polarization wave can be obtained by the samemethod.Thus,this note only discusses the solution of H_x-polarization wave.展开更多
This paper describes an efficient improvement of the multipole accelerated boundary element method for 3-D capacitance extraction. The overall relations between the positions of 2-D boundary elements are considered in...This paper describes an efficient improvement of the multipole accelerated boundary element method for 3-D capacitance extraction. The overall relations between the positions of 2-D boundary elements are considered instead of only the relations between the center-points of the elements, and a new method of cube partitioning is introduced. Numerical results are presented to demonstrate that the method is accurate and has nearly linear computational growth as O(n), where n is the number of panels/boundary elements. The proposed method is more accurate and much faster than Fastcap.展开更多
Using the boundary element method, the numerical modeling problem of three-dimensional terrain effect on magnetotelluric (MT) field is solved. This modeling technique can be run on PC in the case of adopting special n...Using the boundary element method, the numerical modeling problem of three-dimensional terrain effect on magnetotelluric (MT) field is solved. This modeling technique can be run on PC in the case of adopting special net division. The result of modeling test for 2-D terrain by this modeling technique is basically coincident with that by 2-D modeling technique, but there is a great difference between the results of 3-D and 2-D modeling for 3-D terrain.展开更多
Flapping plates of typical fishlike tail shapes are simulated to investigate their locomotion performance using the multi-block Lattice Boltzmann Method （LBM） and Immersed Boundary （IB） method. Numerical results s...Flapping plates of typical fishlike tail shapes are simulated to investigate their locomotion performance using the multi-block Lattice Boltzmann Method （LBM） and Immersed Boundary （IB） method. Numerical results show that fishlike forked configurations have better locomotion performance compared with unforked plates. Based on our results, the caudal fin in carangi- form mode has greater thrust, and the lunate tail fin in thtmniform mode has higher efficiency. These findings are qualitatively con- sistent with biological observations of fish swimming. Analysis of wake topology shows that the wake of the forked plate consists of a chain of alternating reverse horseshoe-like vortical structures. These structures induce a backward jet and generate a positive thrust. Moreover, this backward jet has a more favorable direction compared with that behind an unforked plate.展开更多
基金National Natural Science Foundation of China under Grant No. 51678390National Natural Science Foundation of China under Grant No. 51708391the Major Science and Technology Projects in Tianjin under Grant No. 18ZXAQSF00110。
文摘To solve seismic wave scattering by a large-scale three-dimensional(3-D) hill topography, a fast parallel indirect boundary element method(IBEM) is developed by proposing a new construction method for the wave field, modifying the generalized minimum residual(GMRES) algorithm and constructing an Open MP plus MPI parallel model. The validations of accuracy and efficiency show that this method can solve 3-D seismic response of a large-scale hill topography for broadband waves, and overcome the weakness of large storage and low efficiency of the traditional IBEM. Based on this new algorithm architecture, taking the broadband scattering of plane SV waves by a large-scale Gaussian-shaped hill of thousands-meters height as an example, the influence of several important parameters is investigated, including the incident frequency, the incident angle and the height-width and length-width ratio of the hill. The numerical results illustrate that the amplification effect on the ground motion by a near-hemispherical hill is more significant than the narrow hill. For low-frequency waves, the scattering effect of the higher hill is more pronounced, and there is only a single peak near the top of the hill. However, for high-frequency waves, rapid spatial variation of displacement amplitude appears on the hill surface.
文摘The development of two boundary element algorithms for solving 3-D, frictional, and linear elastostatic contact problems is reported in this paper. The algorithms employ nonconforming discreti- zations for solving 3-D boundary element models, which provide much needed flexibility in the bound- ary element modeling for 3-D contact problems. These algorithms are implemented in a new 3-D boundary element code and verified using several examples. For the numerical examples studied, the results using the new boundary element algorithms match very well with the results using a commercial finite element code, and clearly demonstrate the feasibility of the new boundary element approach for 3-D contact analysis.
文摘THE modeling of 2-D terrain effect on magnetotelluric(MT)has been solved using finite ele-ment method(FEM)and boundary element method(BEM).This note uses BEM tosolve the problem.The solution of H_x- and E_x-polarization wave can be obtained by the samemethod.Thus,this note only discusses the solution of H_x-polarization wave.
文摘This paper describes an efficient improvement of the multipole accelerated boundary element method for 3-D capacitance extraction. The overall relations between the positions of 2-D boundary elements are considered instead of only the relations between the center-points of the elements, and a new method of cube partitioning is introduced. Numerical results are presented to demonstrate that the method is accurate and has nearly linear computational growth as O(n), where n is the number of panels/boundary elements. The proposed method is more accurate and much faster than Fastcap.
基金Project supported by the National Natural Science Foundation of China
文摘Using the boundary element method, the numerical modeling problem of three-dimensional terrain effect on magnetotelluric (MT) field is solved. This modeling technique can be run on PC in the case of adopting special net division. The result of modeling test for 2-D terrain by this modeling technique is basically coincident with that by 2-D modeling technique, but there is a great difference between the results of 3-D and 2-D modeling for 3-D terrain.
基金the National Natural Science Foundation of China (Grant No. 10832010)the Innovation Project of the Chinese Academy of Sciences (Grant No. KJCX2-YW-L05)the 111 Project (Grant No. B07033)
文摘Flapping plates of typical fishlike tail shapes are simulated to investigate their locomotion performance using the multi-block Lattice Boltzmann Method （LBM） and Immersed Boundary （IB） method. Numerical results show that fishlike forked configurations have better locomotion performance compared with unforked plates. Based on our results, the caudal fin in carangi- form mode has greater thrust, and the lunate tail fin in thtmniform mode has higher efficiency. These findings are qualitatively con- sistent with biological observations of fish swimming. Analysis of wake topology shows that the wake of the forked plate consists of a chain of alternating reverse horseshoe-like vortical structures. These structures induce a backward jet and generate a positive thrust. Moreover, this backward jet has a more favorable direction compared with that behind an unforked plate.
