lecture 2. how to model: numerical methods...lecture 2. how to model: numerical methods outline •...
TRANSCRIPT
![Page 1: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/1.jpg)
Lecture 2. How to model:
Numerical methods
Outline
• Brief overview and comparison of methods
• FEM LAPEX
• FEM SLIM3D
• Petrophysical modeling
• Supplementary: details for SLIM3D
![Page 2: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/2.jpg)
chemijij
effii
p HAx
T
xDt
DTC
1)(
ij
IIeff
ij
k
kij
i
j
j
id
TPDt
D
Gx
v
x
v
x
vij
),,(2
1
2
1
3
1)(
2
1
Dt
Dvg
xx
P ii
j
ij
i
Full set of equations
mass
momentum
energy
01
i
i
x
v
Dt
DT
Dt
DP
K
![Page 3: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/3.jpg)
Final effective viscosity
IIPNL
eff
/)(2
1
exp
LL L II
HB
RT
exp
n NN N II
HB
RT
))1(exp(P
IIPPP
RT
HB
PcII 0 if
0 PcII if
![Page 4: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/4.jpg)
Boundary conditions
![Page 5: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/5.jpg)
General case
D
S
![Page 6: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/6.jpg)
Boundary conditions
Kinematic boundary conditions
Dynamic boundary conditions:
Free surface
Free slip
![Page 7: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/7.jpg)
Numerical methods
![Page 8: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/8.jpg)
According to the type of parameterization in space:
FDM, FEM, FVM, SM, BEM etc.
According to the type of parameterization in time:
Explicit, Implicit
According to how mesh changes (if) within a
deforming body:
Lagrangian, Eulerian, Arbitrary Lagrangian
Eulerian (ALE)
![Page 9: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/9.jpg)
Brief Comparison of Methods
Finite Difference Method
(FDM) :
FDM approximates an
operator (e.g., the
derivative)
Finite Element Method
(FEM) :
FEM uses exact operators
but approximates the
solution basis functions.
![Page 10: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/10.jpg)
FD Staggered grid
![Page 11: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/11.jpg)
Finite ElementsTetrahedron
Hexahedron
![Page 12: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/12.jpg)
Finite Elements
Non-uniform
meshes
![Page 13: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/13.jpg)
![Page 14: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/14.jpg)
Spectral Methods (SM):
Spectral methods use
global basis functions to
approximate a solution
across the entire
domain.
Finite Element Methods
(FEM):
FEM use compact basis
functions to approximate
a solution on individual
elements.
Brief Comparison of Methods
![Page 15: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/15.jpg)
Explicit vrs. Implicit
),( tXFdt
dX
![Page 16: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/16.jpg)
Explicit vrs. Implicit
),( tXFdt
dX
)2/),2/(()()(
ttttXFt
tXttX
Should be:
![Page 17: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/17.jpg)
Explicit vrs. Implicit
))),(()()(
ttXFt
tXttX
),( tXFdt
dX
Explicit approximation:
![Page 18: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/18.jpg)
Explicit vrs. Implicit
tttXFtXttX ))),(()()(
),( tXFdt
dX
Explicit approximation:
![Page 19: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/19.jpg)
2-D Thermomechanical Modelling
Explicit finite element algorithm
Basic calculational cycle:m ·d /dt = V F F
- solution of full dynamic equation of motion- calculations in Lagrangian coordinates- remeshing when grid is too distorted- no problems with highly non-linear rheology
General model setupComplex visco-elasto-plastic rheologyT=0, =0xz = zzT or , 0 xz zz = , - Archim.force T/ z = const T/ x=00xz = T/ x=00xz = VxVx
Governing equations:
tLAxTxtTC plasticityCoulombMohrordtdG xvKtp gxxptv
ijijiip ijijij ii ijijiiinert
)(2121
tectoi
iner
iineri
j
ij
i
FDt
Dv
Dt
Dvg
xx
p
,
Dynamic relaxation:
Modified FLAC = LAPEX
(Babeyko et al, EPSL2002)
![Page 20: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/20.jpg)
mF
ij ij
ijij
a = F/m
x, ij ij
ij
F
Explicit finite element method
![Page 21: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/21.jpg)
21
Markers track material and history properties
Upper crust
Lower crust
Mantle
![Page 22: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/22.jpg)
Benchmark: Rayleigh-Taylor instabilityvan Keken et al. (1997)
![Page 23: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/23.jpg)
Sand-box benchmark movie
Friction angle 19°
Weaker layer
Movie attached
![Page 24: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/24.jpg)
Sand-box benchmark movie
Movie attached
![Page 25: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/25.jpg)
2D models
2.5D models
3D- models
Fully 3D models
/ x / x / x3 3 1,2
v3 , / x3 ij
Current state.
no restrictions
v3 , / x3
x1
x2
x1
x2
x3
x1
x2
x3
x1
x2
x3
v -i velocity vector component stress tensor component, ij -
![Page 26: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/26.jpg)
Lagrangian 2.5D FE
Eulerian FD
x1
x2
x3
Simplified 3D concept.
![Page 27: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/27.jpg)
Explicit method vs. implicit
• Advantages
– Easy to implement, small computational efforts per
time step.
– No global matrices. Low memory requirements.
– Even highly nonlinear constitutive laws are always
followed in a valid physical way and without additional
iterations.
– Straightforward way to add new effects (melting,
shear heating, . . . . )
– Easy to parallelize.
• Disadvantages
– Small technical time-step (order of a year)
![Page 28: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/28.jpg)
Implicit ALE FEM SLIM3D(Popov and Sobolev, 2008)
![Page 29: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/29.jpg)
Physical background
Balance equations
1ˆ Elastic strain:
2
1 Viscous strain:
2
Plastic strain:
el vs pl
ij ij ij ij
el
ij ij
vs
ij ij
eff
pl
ij
ij
G
Q
Deformation mechanisms
Popov and Sobolev (2008)
Mohr-Coulomb
Momentum: 0
Energy:
ij
i
j
i
i
g zx
qDUr
Dt x
![Page 30: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/30.jpg)
Numerical background
Discretization by
Finite Element Method
Fast implicit time stepping
+ Newton-Raphson solver
Remapping of
entire fields by
Particle-In-Cell
technique
Arbitrary Lagrangian-Eulerian
kinematical formulation
Popov and Sobolev (2008)
1
1
k k k ku u K r
r Residual Vector
rK Tangent Matrixu
![Page 31: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/31.jpg)
Numerical benchmarks
![Page 32: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/32.jpg)
Numerical benchmarksMovie attached
![Page 33: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/33.jpg)
Numerical benchmarksMovie attached
![Page 34: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/34.jpg)
Numerical benchmarksMovie attached
![Page 35: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/35.jpg)
Numerical benchmarksMovie attached
![Page 36: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/36.jpg)
Numerical benchmarksMovie attached
![Page 37: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/37.jpg)
Solving Stokes equations with code Rhea
(adaptive mesh refinement similar to ASPECT)
Burstedde et al.,2008-2010
![Page 38: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/38.jpg)
Mesh refinement: octree discretization
![Page 39: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/39.jpg)
Solving Stokes equations with code Rhea
Stadler et al., 2010
![Page 40: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/40.jpg)
Open codes
![Page 41: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/41.jpg)
CitComCU. A finite element E parallel code capable of modelling thermo-chemicalconvection in a 3-D domain appropriate for convection within the Earth’s mantle.
Developed from CitCom (Moresi and Solomatov, 1995; Moresi et al., 1996).
CitComS.A finite element E code designed to solve thermal convection problems relevant to Earth’s mantle in 3-D spherical geometry, developed from CitCom by
Zhong et al.(2000).
Ellipsis3D. A 3-D particle-in-cell E finite element solid modelling code for viscoelastoplastic materials, as described in O’Neill et al. (2006).
Gale. An Arbitrary Lagrangian Eulerian (ALE) code that solves problems related to
orogenesis, rifting, and subduction with coupling to surface erosion models. This
is an application of the Underworld platform listed below.
PyLith . A finite element code for the solution of viscoelastic/ plastic deformation that
was designed for small-strain lithospheric modeling problems.
SNAC is a L explicit finite difference code for modelling a finitely deforming elasto-
visco-plastic solid in 3D.
------------------------------------------------------------------------------------------------------
Available from http://milamin.org/.
MILAMIN. A finite element method implementation in MATLAB that is capable of
modelling viscous flow with large number of degrees of freedom on a normal
computer Dabrowski et al. (2008).
Available from CIG (http://geodynamics.org )
![Page 42: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/42.jpg)
Open code Aspect
![Page 43: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/43.jpg)
Ax
T
xDt
DTC IIII
ii
p
)(
ij
eff
ij
i
j
j
iij
Dt
D
Gx
v
x
v
2
1
2
1)(
2
1
Dt
DvgTP
xx
P ii
j
ij
i
),(
Full set of equations
mass
momentum
energy
01
i
i
x
v
Dt
DT
Dt
DP
K
![Page 44: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/44.jpg)
Petrophysical modeling
![Page 45: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/45.jpg)
Goals of the petrophysical modeling
To establish link between rock composition and its physical
properties.
Direct problems:prediction of the density and seismic structure (also
anisotropic)
incorporation in the thermomechanical modeling
Inverse problem:interpretation of seismic velocities in terms of
composition
![Page 46: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/46.jpg)
Petrophysical modeling
Internally-consistent dataset of thermodynamic
properties of minerals and solid solutions(Holland and Powell ‘90, Sobolev and Babeyko ’94)
Gibbs free energy minimization algorithmAfter de Capitani and Brown ‘88
Equilibrium mineralogical composition of a rockgiven chemical composition and PT-conditions
Density and elastic properties optionally with cracks and anisotropy
SiO2
Al2O3
Fe2O3
MgO + (P,T)
CaO
FeO
Na2O
K2O
![Page 47: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/47.jpg)
Gibbs energy
![Page 48: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/48.jpg)
Solid solutions model
![Page 49: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/49.jpg)
Density P-T diagram for average gabbro composition
![Page 50: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/50.jpg)
Supplement: details for FEM SLIM3D
(Popov and Sobolev, PEPI, 2008)
![Page 51: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/51.jpg)
Finite element discretization (SLIM3D)
![Page 52: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/52.jpg)
Time discretization
and nonlinear solution (SLIM3D)
![Page 53: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/53.jpg)
Objective stress integration (SLIM3D)
![Page 54: Lecture 2. How to model: Numerical methods...Lecture 2. How to model: Numerical methods Outline • Brief overview and comparison of methods ... FDM, FEM, FVM, SM, BEM etc. According](https://reader033.vdocument.in/reader033/viewer/2022060916/60a9496cf1809f525f2392aa/html5/thumbnails/54.jpg)
Linearization and tangent operator
(SLIM3D)