首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 15 毫秒
1.
    
A D‐BEM approach, based on time‐weighting residuals, is developed for the solution of two‐dimensional scalar wave propagation problems. The modified basic equation of the D‐BEM formulation is generated by weighting, with respect to time, the basic D‐BEM equation, under the assumption of linear and cubic time variation for the potential and for the flux. A constant time‐weighting function is adopted. The time integration reduces the order of the time‐derivative that appears in the domain integral; as a consequence, the initial conditions are directly taken into account. An assessment of the potentialities of the proposed formulation is provided by the examples included at the end of the work. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

2.
    
A modified version of an exact Non‐reflecting Boundary Condition (NRBC) first derived by Grote and Keller is implemented in a finite element formulation for the scalar wave equation. The NRBC annihilate the first N wave harmonics on a spherical truncation boundary, and may be viewed as an extension of the second‐order local boundary condition derived by Bayliss and Turkel. Two alternative finite element formulations are given. In the first, the boundary operator is implemented directly as a ‘natural’ boundary condition in the weak form of the initial–boundary value problem. In the second, the operator is implemented indirectly by introducing auxiliary variables on the truncation boundary. Several versions of implicit and explicit time‐integration schemes are presented for solution of the finite element semidiscrete equations concurrently with the first‐order differential equations associated with the NRBC and an auxiliary variable. Numerical studies are performed to assess the accuracy and convergence properties of the NRBC when implemented in the finite element method. The results demonstrate that the finite element formulation of the (modified) NRBC is remarkably robust, and highly accurate. Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

3.
    
Asymptotic and exact local radiation boundary conditions (RBC) for the scalar time‐dependent wave equation, first derived by Hagstrom and Hariharan, are reformulated as an auxiliary Cauchy problem for each radial harmonic on a spherical boundary. The reformulation is based on the hierarchy of local boundary operators used by Bayliss and Turkel which satisfy truncations of an asymptotic expansion for each radial harmonic. The residuals of the local operators are determined from the solution of parallel systems of linear first‐order temporal equations. A decomposition into orthogonal transverse modes on the spherical boundary is used so that the residual functions may be computed efficiently and concurrently without altering the local character of the finite element equations. Since the auxiliary functions are based on residuals of an asymptotic expansion, the proposed method has the ability to vary separately the radial and transverse modal orders of the RBC. With the number of equations in the auxiliary Cauchy problem equal to the transverse mode number, this reformulation is exact. In this form, the equivalence with the closely related non‐reflecting boundary condition of Grote and Keller is shown. If fewer equations are used, then the boundary conditions form high‐order accurate asymptotic approximations to the exact condition, with corresponding reduction in work and memory. Numerical studies are performed to assess the accuracy and convergence properties of the exact and asymptotic versions of the RBC. The results demonstrate that the asymptotic formulation has dramatically improved accuracy for time domain simulations compared to standard boundary treatments and improved efficiency over the exact condition. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

4.
    
A family of radiation boundary conditions for the wave equation is derived by truncating a rational function approximation of the corresponding plane wave representation, and it is demonstrated how these boundary conditions can be formulated in terms of fictitious surface densities, governed by second‐order wave equations on the radiating surface. Several well‐established radiation boundary conditions appear as special cases, corresponding to different choices of the coefficients in the rational approximation. The relation between these choices is established, and an explicit formulation in terms of selected directions with ideal transmission is presented. A mechanical interpretation of the fictitious surface densities enables identification of suitable conditions at corners and boundaries of the radiating surface. Numerical examples illustrate excellent results with one or two fictitious layers with suitable corner and boundary conditions. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

5.
A time-discontinuous Galerkin space–time finite element method is formulated for the exterior structural acoustics problem in two space dimensions. The problem is posed over a bounded computational domain with local time-dependent radiation (absorbing) boundary conditions applied to the fluid truncation boundary. Absorbing boundary conditions are incorporated as ‘natural’ boundary conditions in the space–time variational equation, i.e. they are enforced weakly in both space and time. Following Bayliss and Turkel, time-dependent radiation boundary conditions for the two-dimensional wave equation are developed from an asymptotic approximation to the exact solution in the frequency domain expressed in negative powers of a non-dimensional wavenumber. In this paper, we undertake a brief development of the time-dependent radiation boundary conditions, establishing their relationship to the exact impedance (Dirichlet-to-Neumann map) for the acoustic fluid, and characterize their accuracy when implemented in our space–time finite element formulation for transient structural acoustics. Stability estimates are reported together with an analysis of the positive form of the matrix problem emanating from the space–time variational equations for the coupled fluid-structure system. Several numerical simulations of transient radiation and scattering in two space dimensions are presented to demonstrate the effectiveness of the space–time method.  相似文献   

6.
    
  相似文献   

7.
    
In this paper, we propose an efficient strategy to compute nearly singular integrals over planar triangles in R 3 arising in boundary element method collocation. The strategy is based on a proper use of various non‐linear transformations, which smooth or move away or quite eliminate all the singularities close to the domain of integration. We will deal with near singularities of the form 1/r, 1/r2 and 1/r3, r=∥ x ? y ∥ being the distance between a fixed near observation point x and a generic point y of a triangular element. Extensive numerical tests and comparisons with some already existing methods show that the approach proposed here is highly efficient and competitive. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

8.
This work is concerned with the development of a D-BEM approach to the solution of 2D scalar wave propagation problems. The time-marching process can be accomplished with the use of the Houbolt method, as usual, or with the use of the Newmark method. Special attention was devoted to the development of a procedure that allows for the computation of the initial conditions contributions. In order to verify the applicability of the Newmark method and also the correctness of the expressions concerned with the computation of the initial conditions contributions, four examples are presented and the D-BEM results are compared with the corresponding analytical solutions.  相似文献   

9.
    
This work is concerned with the development of different domain‐BEM (D‐BEM) approaches to the solution of two‐dimensional diffusion problems. In the first approach, the process of time marching is accomplished with a combination of the finite difference and the Houbolt methods. The second approach starts by weighting, with respect to time, the basic D‐BEM equation, under the assumption of linear and constant time variation for the temperature and for the heat flux, respectively. A constant time weighting function is adopted. The time integration reduces the order of the time derivative that appears in the domain integral; as a consequence, the initial conditions are directly taken into account. Four examples are presented to verify the applicability of the proposed approaches, and the D‐BEM results are compared with the corresponding analytical solutions.Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

10.
    
A novel substructuring method is developed for the coupling of boundary element and finite element subdomains in order to model three‐dimensional multi‐region elastodynamic problems in the time domain. The proposed procedure is based on the interface stiffness matrix approach for static multi‐region problems using variational principles together with the concept of Duhamel integrals. Unit impulses are applied at the boundary of each region in order to evaluate the impulse response matrices of the Duhamel (convolution) integrals. Although the method is not restricted to a special discretization technique, the regions are discretized using the boundary element method combined with the convolution quadrature method. This results in a time‐domain methodology with the advantages of performing computations in the Laplace domain, which produces very accurate and stable results as verified on test examples. In addition, the assembly of the boundary element regions and the coupling to finite elements are greatly simplified and more efficient. Finally, practical applications in the area of soil–structure interaction and tunneling problems are shown. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

11.
    
We present a new solution to accelerate the boundary integral equation method (BIEM). The calculation time of the BIEM is dominated by the evaluation of the layer potential in the boundary integral equation. We performed this task using MDGRAPE‐2, a special‐purpose computer designed for molecular dynamics simulations. MDGRAPE‐2 calculates pairwise interactions among particles (e.g. atoms and ions) using hardwired‐pipeline processors. We combined this hardware with an iterative solver. During the iteration process, MDGRAPE‐2 evaluates the layer potential. The rest of the calculation is performed on a conventional PC connected to MDGRAPE‐2. We applied this solution to the Laplace and Helmholtz equations in three dimensions. Numerical tests showed that BIEM is accelerated by a factor of 10–100. Our rather naive solution has a calculation cost of O(N2 × Niter), where N is the number of unknowns and Niter is the number of iterations. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

12.
    
An accurate absorbing boundary condition (ABC) is developed in frequency domain for finite element analysis of scalar wave propagation in unbounded layered half-space. The proposed ABC is H-shaped line that consists of two parts: a new ABC at horizontal bottom boundary of finite domain to replace semiinfinite strip below horizontal boundary and between two vertical boundaries, and a general consistent ABC at vertical lateral boundary to replace semiinfinite layered half-space outside vertical boundary. The key point for constructing the ABC is that a new continued fraction (CF) is presented to expand dynamic stiffness of underlying half-space, and the CF-based stress-displacement relationship is then transformed into an auxiliary variable system with square of horizontal wavenumber. The ABC has only one undetermined real parameter that is the CF-order independent of frequency and incidence angle of propagating outgoing waves. The parameter can be chosen relatively small value to achieve an accurate ABC. Moreover, the ABC can couple seamlessly with finite element method of finite domain. The finite domain can be chosen very small size due to high accuracy of the ABC. Numerical examples are finally given to demonstrate the effectiveness of the ABC.  相似文献   

13.
    
In this paper, a new effective boundary node method is presented for the solution of acoustic problems, directly in time domain, using exponential basis functions. Unlike many other methods using boundary information, the final coefficient matrix is sparse. The formulation is well suited for domains whose extent is relatively larger than the distance traveled by the acoustic wave in an increment of time. The exponential basis functions used satisfy the time‐space governing equation. This helps to choose a relatively large time increment and a moderate number of boundary points, which leads to reduction of computation time. The computation is performed incrementally using a weighted residual in time. Through a series of numerical examples, it is shown that the method, when combined with a domain decomposition strategy, is effectively capable of solving various 1‐ to 3‐dimensional acoustic problems.  相似文献   

14.
    
In this paper the numerical solution for ideal irrotational incompressible flow around axisymmetric 3D obstacles is discussed with the artificial boundary method. By introducing an artificial boundary, we divide the exterior unbounded domain into a bounded part and an unbounded part. After setting up a proper artificial boundary condition, we get an approximate reduced problem defined on the bounded part. Both non‐local and local artificial boundary conditions are designed. Numerical experiment is also presented, and its result demonstrates the effectiveness of these artificial boundary conditions. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

15.
    
This article advocates a general procedure for the numerical investigation of pseudo‐rigid bodies. The equations of motion for pseudo‐rigid bodies are shown to be mathematically equivalent to those corresponding to certain constant‐strain finite element approximations for general deformable continua. A straightforward algorithmic implementation is achieved in a classical finite element framework. Also, a penalty formulation is suggested for modelling contact between pseudo‐rigid bodies. Representative planar simulations using a non‐linear elastic model demonstrate the predictive capacity of the pseudo‐rigid theory, as well as the robustness of the proposed computational procedure. Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

16.
    
We develop a three-dimensional shape optimization (SO) framework for the wave equation with taking the unsteadiness into account. Resorting to the adjoint variable method, we derive the shape derivative (SD) with respect to a deformation (perturbation) of an arbitrary point on the target surface of acoustic scatterers. Successively, we represent the target surface with non-uniform rational B-spline patches and then discretize the SD in term of the associated control points (CPs), which are useful for manipulating a surface. To solve both the primary and adjoint problems, we apply the time-domain boundary element method (TDBEM) because it is the most appropriate when the analysis domain is the ambient air and thus infinitely large. The issues of the severe computational cost and instability of the TDBEM are resolved by exploiting the fast and stable TDBEM proposed by the present authors. Instead, since the TDBEM is mesh-based and employs the piecewise-constant element for space, we introduce some approximations in evaluating the discretized SD from the two solutions of TDBEM. By regarding the evaluation scheme as the computation of the gradient of the objective functional, given as the summation of the absolute value of the sound pressure over the predefined observation points, we can solve SO problems with a gradient-based non-linear optimization solver. To assess the developed SO system, we performed several numerical experiments from the perspective of verification and application with satisfactory results.  相似文献   

17.
论Helmholtz方程的一类边界积分方程的合理性   总被引：5，自引：0，他引：5  
本文导出了Helmholtz 方程超定边值问题有解的一个充要条件,和用非解析开拓法证明了文[1]中的Helmholtz 方程在外域中的解的边界积分表示式的合理性,并将此类边界积分表示式推广用于带空洞的有限域。这样就比较严密而又浅近地证明了基于该表示式建立起来的间接变量和直接变量边界积分方程的合理性。  相似文献   

18.
    
In this paper we consider Dirichlet or Neumann wave propagation problems reformulated in terms of boundary integral equations with retarded potential. Starting from a natural energy identity, a space–time weak formulation for 1D integral problems is briefly introduced, and continuity and coerciveness properties of the related bilinear form are proved. Then, a theoretical analysis of an extension of the introduced formulation for 2D problems is proposed, pointing out the novelty with respect to existing literature results. At last, various numerical simulations will be presented and discussed, showing unconditional stability of the space–time Galerkin boundary element method applied to the energetic weak problem. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

19.
    
One of the important issues in the implementation of the differential quadrature method is the imposition of the given boundary conditions. There may be multiple boundary conditions involving higher‐order derivatives at the boundary points. The boundary conditions can be imposed by modifying the weighting coefficient matrices directly. However, the existing method is not robust and is known to have many limitations. In this paper, a systematic procedure is proposed to construct the modified weighting coefficient matrices to overcome these limitations. The given boundary conditions are imposed exactly. Furthermore, it is found that the numerical results depend only on those sampling grid points where the differential quadrature analogous equations of the governing differential equations are established. The other sampling grid points with no associated boundary conditions are not essential. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

20.
    
We consider the efficient numerical solution of the three‐dimensional wave equation with Neumann boundary conditions via time‐domain boundary integral equations. A space‐time Galerkin method with C∞‐smooth, compactly supported basis functions in time and piecewise polynomial basis functions in space is employed. We discuss the structure of the system matrix and its efficient parallel assembly. Different preconditioning strategies for the solution of the arising systems with block Hessenberg matrices are proposed and investigated numerically. Furthermore, a C++ implementation parallelized by OpenMP and MPI in shared and distributed memory, respectively, is presented. The code is part of the boundary element library BEM4I. Results of numerical experiments including convergence and scalability tests up to a thousand cores on a cluster are provided. The presented implementation shows good parallel scalability of the system matrix assembly. Moreover, the proposed algebraic preconditioner in combination with the FGMRES solver leads to a significant reduction of the computational time. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号