boundary element method (bem)
TRANSCRIPT
Boundary Element Method (BEM)
Zoran Ilievski
Wednesday 28th June, 2006HG 6.96 (TU/e)
2
Talk Overview• The idea of BEM and its advantages
• The 2D potential problem
• Numerical implementation
3
Idea of BEM
4
Idea of BEM
5
Advantages of BEM1) Reduction of problem dimension by 1
• Less data preparation time.• Easier to change the applied mesh.• Useful for problems that require re-meshing.
6
Advantages of BEM2) High Accuracy • Stresses are accurate as there are
no approximations imposed on the solution in interior domain points.
• Suitable for modeling problems of rapidly changing stresses.
7
Advantages of BEM3) Less computer time and storage
• For the same level of accuracy as other methods BEM uses less number of nodes and elements.
8
Advantages of BEM4) Filter out unwanted information.• Internal points of the domain are
optional.
• Focus on particularinternal region.
• Further reduces computer time.
9
Advantages of BEM1. Reduction of problem dimension by 1.2. High Accuracy.3. Less computer time and storage.4. Filter out unwanted information and so
focus on section of the domain you are interested in.
BEM is an attractive option.
10
The 2D potential problem• Where can BEM be applied?• Two important functions.• Description of the domain.• Mapping of higher to lower dimensions.• Satisfaction of the Laplace equations
and how to deal with a singularity.• The boundary integral equation (BIE)
11
The 2D potential problemWhere can BEM be applied?
Where any potential problem is governed by a differential equationthat satisfies the Laplace equation.
(or any other behavior that has a related fundamental solution)
e.g. The following can be analyzed with the Laplace equation: fluid flow, torsion of bars,
diffusion and steady state heat conduction.
12
The 2D potential problem
02
2
2
22
yx
.2
The Laplace equation for 2D
yx,
Laplacian operatorPotential functionCartesian coordinate
axis
13
The 2D potential problemTwo important functions.
)(2 PQ
The function describing the property under analysis. e.g. heat. (Unknown)
The fundamental solution of the Laplace equation. (These are well known)
14
The 2D potential problem
),(1ln
21),(
QprQp
),(1ln
21),(
QprQp
Fundamental solution of the 2DLaplace equation for a concentrated source point at p is
22),( QpQp yYxXQpr
Where
Description of the domain
15
The 2D potential problemMapping of higher to lower dimensions
dnn
dAA
22
• Boundary of any domain is of a dimension 1 less than of the domain.
• In BEM the problem is moved from within the domain to its boundary.
• This means you must, in this case, map Area to Line.• The well known ‘Greens Second Identity’ is used to do this.
,
nn
have continuous 1st and 2nd derivatives.unknown potential at any point.
known fundamental solution at any point.unit outward
normal.derivative in the direction of
normal.
16
The 2D potential problemSatisfying the Laplace equation
The unknown
02 will satisfy
everywhere in the solution domain.
The known fundamental solution
satisfies
02 everywhere except the point p where it is singular.
),(1ln
21),(
QprQp
22),( QpQp yYxXQpr
17
The 2D potential problemHow to deal with the singularity
d
nndA
AA
22
• Surround p with a small circle of radius ε , thenexamine solution as ε 0
• New area is (A – Aε )• New boundary is (Γ + Γε )
Within area (A – Aε)
02 02 &
The left hand side of the equation is now 0 and the right is now …
18
The 2D potential problemHow to deal with the singularity
d
nnd
nn
0
The second term must be evaluated and to do this let
dd
rnr
rn
21.
And use the fact
that
19
The 2D potential problemHow to deal with the singularity
.1)2(21
1ln121 2
0
dn
dnn
02/11
)(PCEvaluated with p in the
domain,on the boundary (Smooth surface),and outside the boundary.
2
)( PC For coarse surfaces
20
The 2D potential problemThe boundary integral equation
)()(),()(),()()( 12 QdQQPKQdn
QPKPPC
),(),(
),(),(
2
1
QPQPKnQPQPK
Where K1 and K2 are the known fundamental solutions and are equal to
2
)( PC
21
The 2D potential problem• BEM can be applied where any potential problem is governed by
a differential equation that satisfies the Laplace equation. In this case the 2D form.
• A potential problem can be mapped from higher to lower
dimension using Green’s second identity.
• Shown how to deal with the case of the singularity point.
• Derived the boundary integral equation (BIE)
)()(),()(),()()( 12 QdQQPKQdn
QPKPPC
22
Numerical Implementation• Dirichlet, Neumann and mixed case.• Discretisation• Reduction to a form Ax=B
23
Numerical ImplementationDirichlet, Neumann and mixed case.
)()(),()(),()()( 12 QdQQPKQdn
QPKPPC
The unknowns of the above are values on the boundary and are n
,
Dirichlet Problem
Neumann Problem
n
is given every point Q on the boundary.
is given every point Q on the boundary.
Mixed case – Either are given at point Q
24
Numerical ImplementationDiscretisation
N
jjjij
N
jjji
ji
jj
dQPKQdQPKnQ
P1
11
221 ),()(),(
)()(
Unknowns
25
Numerical ImplementationDiscretisation
N
jijj
N
jij
ji KQK
nQ
P1
11
221 )(
)()(
jjiij dQPKKj
),(22 j
jjiij dQPKK ),(11
Unknowns
Let
26
Numerical Implementation
N
j
jijj
N
jijij n
QKQK
12
121
1
)()(
jiwhenQP ji )()(
BzAx
27
Numerical Implementation
BzAx
cAx
Dirichlet ProblemBzc
Neumann Problem
Mixed case
Matrix A and vector C are known
Matrix B and vector C are known
Unknowns and knowns can be separated in to same form as above
28
Numerical Implementation As each point p in the domain is expressed in terms of
the boundary values, once all boundary values are known ANYpotential value within the domain can now be found.
)()(),()(),()()( 12 QdQQPKQdn
QPKPPC
29
Boundary Element Method (BEM)
THE END
Book: The Boundary Element Method in Engineering A.A.BECKER