mathematical modeling of particle nucleation and growth in...
TRANSCRIPT
![Page 1: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/1.jpg)
1
Delft Institute of Applied MathematicsDelft University of Technology
Mathematical modeling of particlenucleation and growth in metallicalloys
Dennis den Ouden
14-12-2009
![Page 2: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/2.jpg)
14-12-2009 2
Delft Institute of Applied Mathematics
Contents
• Something about metallurgy• A model for nucleation• A model for deformations• Combining the models• Results• Conclusion
![Page 3: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/3.jpg)
14-12-2009 3
Delft Institute of Applied Mathematics
Something about alloys
Binary Ternary Quaternary Complex
Quasi-Binary
Quasi-Binary
Quasi-Binary Quasi-Ternary
Quasi-Ternary
Quasi-Quaternary
![Page 4: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/4.jpg)
14-12-2009 3
Delft Institute of Applied Mathematics
Something about alloys
Binary Ternary Quaternary Complex
Quasi-Binary
Quasi-Binary
Quasi-Binary Quasi-Ternary
Quasi-Ternary
Quasi-Quaternary
![Page 5: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/5.jpg)
14-12-2009 4
Delft Institute of Applied Mathematics
Something about deformations
![Page 6: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/6.jpg)
14-12-2009 4
Delft Institute of Applied Mathematics
Something about deformations
![Page 7: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/7.jpg)
14-12-2009 5
Delft Institute of Applied Mathematics
Models for nucleation
Two models, with small differences• Myhr and Grong (2000)• Robson et al. (2003)
![Page 8: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/8.jpg)
14-12-2009 5
Delft Institute of Applied Mathematics
Models for nucleation
Two models, with small differences• Myhr and Grong (2000)• Robson et al. (2003)
Comparison:• Basic model == Myhr and Grong (2000)• Adapted using Robson et al. (2003)
![Page 9: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/9.jpg)
14-12-2009 6
Delft Institute of Applied Mathematics
Governing DE
Changein time
=Change
by growth+
Productionrate
![Page 10: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/10.jpg)
14-12-2009 6
Delft Institute of Applied Mathematics
Governing DE
Changein time
=Change
by growth+
Productionrate
⇓ ⇓ ⇓
∂N
∂t= −
∂ (Nv)
∂r+ S
![Page 11: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/11.jpg)
14-12-2009 7
Delft Institute of Applied Mathematics
Unknowns
• Growth rate v
• Production term S
![Page 12: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/12.jpg)
14-12-2009 8
Delft Institute of Applied Mathematics
Growth of spherical particles
Directly influence by
• Mean concentration C
• Interface concentrationCi
• Internal concentrationCp
• Diffusion coefficient D
Cp C
Ci
Matrix
Particle
![Page 13: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/13.jpg)
14-12-2009 8
Delft Institute of Applied Mathematics
Growth of spherical particles
v(r, t) =C − Ci
Cp − Ci
D
r Cp C
Ci
Matrix
Particle
![Page 14: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/14.jpg)
14-12-2009 9
Delft Institute of Applied Mathematics
A special particle
A particle with radius r∗ that will neither grow ordissolve:
v(r∗) = 0
![Page 15: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/15.jpg)
14-12-2009 9
Delft Institute of Applied Mathematics
A special particle
A particle with radius r∗ that will neither grow ordissolve:
v(r∗) = 0
Solved for r∗:
r∗ =2γαβVm
RT
(
ln
(
C
Ce
))−1
![Page 16: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/16.jpg)
14-12-2009 9
Delft Institute of Applied Mathematics
A special particle
A particle with radius r∗ that will neither grow ordissolve:
v(r∗) = 0
Solved for r∗:
r∗ =2γαβVm
RT
(
ln
(
C
Ce
))−1
Definition: Critical particle radius
![Page 17: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/17.jpg)
14-12-2009 10
Delft Institute of Applied Mathematics
Production term
• Indicates the number of particles thatnucleate over the whole domain
• Influenced by critical radius r∗
• Influenced by nucleation rate j
![Page 18: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/18.jpg)
14-12-2009 10
Delft Institute of Applied Mathematics
Production term
• Indicates the number of particles thatnucleate over the whole domain
• Influenced by critical radius r∗
• Influenced by nucleation rate j
Kampmann et al. (1987):
S(r, t) =
{
j(t) if r = r∗ + ∆r∗,0 otherwise.
![Page 19: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/19.jpg)
14-12-2009 11
Delft Institute of Applied Mathematics
Nucleation rate
The number of particles that nucleate with radiusr∗ + ∆r∗:
• Influenced by diffusion• Only if some barrier has been overcome
![Page 20: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/20.jpg)
14-12-2009 11
Delft Institute of Applied Mathematics
Nucleation rate
The number of particles that nucleate with radiusr∗ + ∆r∗:
• Influenced by diffusion• Only if some barrier has been overcome
j = j0 exp
(
−∆G∗
het
RT
)
exp
(
−Qd
RT
)
![Page 21: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/21.jpg)
14-12-2009 12
Delft Institute of Applied Mathematics
Nucleation energy barrier
• Chemical composition• Misfit strain energy
![Page 22: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/22.jpg)
14-12-2009 12
Delft Institute of Applied Mathematics
Nucleation energy barrier
• Chemical composition• Misfit strain energy
∆G∗het =
A30
(∆Gv + ∆Gm
s)2
![Page 23: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/23.jpg)
14-12-2009 13
Delft Institute of Applied Mathematics
Model overview
• Governing DE:
∂N
∂t= −
∂ (Nv)
∂r+ S
• Source term:
S(r, t) =
{
j(t) if r = r∗ + ∆r∗,0 otherwise.
![Page 24: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/24.jpg)
14-12-2009 13
Delft Institute of Applied Mathematics
Model overview
• Governing DE:
∂N
∂t= −
∂ (Nv)
∂r+ S
• Growth rate:
v =C − Ci
Cp − Ci
D
r
![Page 25: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/25.jpg)
14-12-2009 14
Delft Institute of Applied Mathematics
Elastic deformations
![Page 26: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/26.jpg)
14-12-2009 14
Delft Institute of Applied Mathematics
Elastic deformations
![Page 27: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/27.jpg)
14-12-2009 15
Delft Institute of Applied Mathematics
Assumptions
Rotation symmetry:• No deformations in tangential direction• No deformation at center axis in radial
direction• All derivatives in tangential direction vanish
![Page 28: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/28.jpg)
14-12-2009 15
Delft Institute of Applied Mathematics
Assumptions
Rotation symmetry:• No deformations in tangential direction• No deformation at center axis in radial
direction• All derivatives in tangential direction vanish
uθ = 0 uη(0, θ, z) = 0∂(.)
∂θ= 0
![Page 29: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/29.jpg)
14-12-2009 16
Delft Institute of Applied Mathematics
Strain and deformation
Chau and Wei (2000):
εηη =∂uη
∂ηεθθ =
uη
η
εzz =∂uz
∂zεηθ = 0
εηz =1
2
(
∂uη
∂z+
∂uz
∂η
)
εθz = 0
![Page 30: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/30.jpg)
14-12-2009 17
Delft Institute of Applied Mathematics
Stress and strain
Hook’s Law:
σαβ = δαβλ (εηη + εθθ + εzz) + 2µεαβ
![Page 31: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/31.jpg)
14-12-2009 17
Delft Institute of Applied Mathematics
Stress and strain
Hook’s Law:
σαβ = δαβ λ (εηη + εθθ + εzz) + 2 µ εαβ
Stiffness matrix Shear modulus
![Page 32: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/32.jpg)
14-12-2009 18
Delft Institute of Applied Mathematics
Force balance
Jaeger et al. (2007):
∂σηη
∂η+
∂σηz
∂z+
σηη − σθθ
η+ bη = 0
∂σηz
∂η+
∂σzz
∂z+
σηz
η+ bz = 0
![Page 33: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/33.jpg)
14-12-2009 19
Delft Institute of Applied Mathematics
Boundary conditions
• Symmetry condition:
uη(0, θ, z) = 0
![Page 34: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/34.jpg)
14-12-2009 19
Delft Institute of Applied Mathematics
Boundary conditions
• Symmetry condition:
uη(0, θ, z) = 0
• Fixed boundaries:
uα(η, θ, z) = 0
![Page 35: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/35.jpg)
14-12-2009 19
Delft Institute of Applied Mathematics
Boundary conditions
• Symmetry condition:
uη(0, θ, z) = 0
• Fixed boundaries:
uα(η, θ, z) = 0
• Moving boundaries:(
σ(η, θ, z))
α· n = fα(η, θ, z)
![Page 36: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/36.jpg)
14-12-2009 20
Delft Institute of Applied Mathematics
Coupling the models
Remember the nucleation energy barrier:
∆G∗het =
A30
(
∆Gv + ∆Gms
)2
Misfit strain energy
![Page 37: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/37.jpg)
14-12-2009 20
Delft Institute of Applied Mathematics
Coupling the models
Remember the nucleation energy barrier:
∆G∗het =
A30
(
∆Gv + ∆Gms
)2
Misfit strain energy
Question:
Is there something like elastic strain energy?
![Page 38: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/38.jpg)
14-12-2009 21
Delft Institute of Applied Mathematics
Coupling the models (2)
Answer:YES!!!
![Page 39: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/39.jpg)
14-12-2009 21
Delft Institute of Applied Mathematics
Coupling the models (2)
Answer:YES!!!
Solution:
∆Gels =
1
2σ : ε
![Page 40: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/40.jpg)
14-12-2009 21
Delft Institute of Applied Mathematics
Coupling the models (2)
Answer:YES!!!
Solution:
∆Gels =
1
2σ : ε
and:
∆G∗het =
A30
(
∆Gv + ∆Gms +∆Gel
s
)2
![Page 41: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/41.jpg)
14-12-2009 22
Delft Institute of Applied Mathematics
Coupling the models (3)
Question:Is there also a reverse coupling?
![Page 42: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/42.jpg)
14-12-2009 22
Delft Institute of Applied Mathematics
Coupling the models (3)
Question:Is there also a reverse coupling?
Answer:YES !!!
![Page 43: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/43.jpg)
14-12-2009 22
Delft Institute of Applied Mathematics
Coupling the models (3)
Question:Is there also a reverse coupling?
Solution by Pal (2005):
µ = µm +
(
15(1 − νm)(µp − µm)
2µp(4 − 5νm) + µm(7 − 5νm)
)
µm f
E = Em + (10β1(1 + νm) + β2(1 − 2νm)) Em f
λ = µE − 2µ
3µ − E
![Page 44: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/44.jpg)
14-12-2009 23
Delft Institute of Applied Mathematics
Recap
Two models:• Nucleation model• Elastic model
Two couplings:• From elastic to nucleation• From nucleation to elastic
![Page 45: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/45.jpg)
14-12-2009 24
Delft Institute of Applied Mathematics
Numerical methods
Nucleation model:• Upwind scheme
• IMEX-θ method with θ = 1
2
![Page 46: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/46.jpg)
14-12-2009 24
Delft Institute of Applied Mathematics
Numerical methods
Nucleation model:• Upwind scheme
• IMEX-θ method with θ = 1
2
(
I −1
2
∆t
∆rAn
)
~Nn+1 =
(
I +1
2
∆t
∆rAn
)
~Nn + ∆t~Sn
![Page 47: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/47.jpg)
14-12-2009 25
Delft Institute of Applied Mathematics
Numerical methods (2)
Elastic model:• Finite Element Method• Linear elements• Use of rotation symmetry
![Page 48: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/48.jpg)
14-12-2009 25
Delft Institute of Applied Mathematics
Numerical methods (2)
0 0.5 1 1.5 2 2.5 3
x 10−3
0
0.005
0.01
0.015
0.02
0.025
0.03z
(m)
η (m)
![Page 49: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/49.jpg)
14-12-2009 25
Delft Institute of Applied Mathematics
Numerical methods (2)
Equation:[
Sηη Sηz
Szη Szz
][
uη
uz
]
=
[
qη
qz
]
![Page 50: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/50.jpg)
14-12-2009 26
Delft Institute of Applied Mathematics
Numerical methods
Algorithm:
1. Set all constants;
2. Set all initial values;
![Page 51: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/51.jpg)
14-12-2009 26
Delft Institute of Applied Mathematics
Numerical methods
Algorithm:
1. Set all constants;
2. Set all initial values;
3. For each time step:(a) Calculate elastic parameters;(b) Build matrices for elastic deformation;(c) Calculate elastic deformations;(d) Calculate elastic strain energy;
![Page 52: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/52.jpg)
14-12-2009 26
Delft Institute of Applied Mathematics
Numerical methods
Algorithm:. . .
3. For each time step:. . .
(e) For each point:i. Calculate nucleation parameters;ii. Calculate matrices for nucleation;iii. Calculate nucleation.
![Page 53: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/53.jpg)
14-12-2009 27
Delft Institute of Applied Mathematics
Simulation
Material:• Aluminum alloy AA 6082• Mg2Si particles
![Page 54: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/54.jpg)
14-12-2009 27
Delft Institute of Applied Mathematics
Simulation
Material:• Aluminum alloy AA 6082• Mg2Si particles
Shape:• Cylindrical• Height 30 millimeter• Radius 3 millimeter
![Page 55: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/55.jpg)
14-12-2009 28
Delft Institute of Applied Mathematics
Simulation (2)
Time• Total of 3000 seconds• Time step of 0.5 seconds
![Page 56: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/56.jpg)
14-12-2009 28
Delft Institute of Applied Mathematics
Simulation (2)
Time• Total of 3000 seconds• Time step of 0.5 seconds
Test:• Tensile test• Bottom axial and radial fixed• Top radial fixed
• Axial force at top of 6 million N/m2
• Sides free
![Page 57: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/57.jpg)
14-12-2009 29
Delft Institute of Applied Mathematics
Typical deformations: Axial
0 0.5 1 1.5 2 2.5 3
x 10−3
0
0.005
0.01
0.015
0.02
0.025
0.03
0
1
2
3
4
5
6
7
8x 10
−6Deformation uz (m)z
(m)
η (m)
![Page 58: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/58.jpg)
14-12-2009 29
Delft Institute of Applied Mathematics
Typical deformations: Radial
0 0.5 1 1.5 2 2.5 3
x 10−3
0
0.005
0.01
0.015
0.02
0.025
0.03
−2.5
−2
−1.5
−1
−0.5
0
0.5
1
1.5
2
2.5
3x 10
−7
z(m
)
η (m)
Deformation uη (m)
![Page 59: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/59.jpg)
14-12-2009 29
Delft Institute of Applied Mathematics
Typical deformations: Energy
0 0.5 1 1.5 2 2.5 3
x 10−3
0
0.005
0.01
0.015
0.02
0.025
0.03
2
4
6
8
10
12
14x 10
4Elastic strain energy (N/m2)z
(m)
η (m)
![Page 60: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/60.jpg)
14-12-2009 30
Delft Institute of Applied Mathematics
Nucleation results: Nucleation rate
0
1
2
3
x 10−3
0
0.01
0.02
0.0310
0
101
102
103
0.5
1
1.5
2
2.5
3
x 10−7Particle nucleation rate
time
(s)
Per
cent
age
(%)
z (m) η (m)
![Page 61: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/61.jpg)
14-12-2009 30
Delft Institute of Applied Mathematics
Nucleation results: Number density
0
1
2
3
x 10−3
0
0.01
0.02
0.0310
0
101
102
103
0.5
1
1.5
2
2.5
3
x 10−7Particle number density
time
(s)
Per
cent
age
(%)
z (m) η (m)
![Page 62: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/62.jpg)
14-12-2009 30
Delft Institute of Applied Mathematics
Nucleation results: Concentration
0
1
2
3
x 10−3
0
0.01
0.02
0.0310
0
101
102
103
−2
−1.8
−1.6
−1.4
−1.2
−1
−0.8
−0.6
−0.4
−0.2
0x 10
−7Mean concentration
time
(s)
Per
cent
age
(%)
z (m) η (m)
![Page 63: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/63.jpg)
14-12-2009 31
Delft Institute of Applied Mathematics
Reflection
Are the results anomalies during simulation?
![Page 64: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/64.jpg)
14-12-2009 31
Delft Institute of Applied Mathematics
Reflection
Are the results anomalies during simulation?
Increase force to test for similar behavior.
F = 6 × 109N
m2
![Page 65: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/65.jpg)
14-12-2009 31
Delft Institute of Applied Mathematics
Reflection
Are the results anomalies during simulation?
Increase force to test for similar behavior.
F = 6 × 109N
m2
Physically no longer elasticity
![Page 66: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/66.jpg)
14-12-2009 32
Delft Institute of Applied Mathematics
New elastic stain energy
0 0.5 1 1.5 2 2.5 3
x 10−3
0
0.005
0.01
0.015
0.02
0.025
0.03
1
2
3
4
5
6
7
8
9
10
11
12x 10
10Elastic strain energy (N/m2)z
(m)
η (m)
![Page 67: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/67.jpg)
14-12-2009 33
Delft Institute of Applied Mathematics
Nucleation results: Number density
0
1
2
3
x 10−3
0
0.01
0.02
0.0310
0
101
102
103
0.005
0.01
0.015
0.02
0.025
0.03
0.035
Particle number density
time
(s)
Per
cent
age
(%)
z (m) η (m)
![Page 68: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/68.jpg)
14-12-2009 34
Delft Institute of Applied Mathematics
Conclusions
• Two separate nucleation models combined
![Page 69: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/69.jpg)
14-12-2009 34
Delft Institute of Applied Mathematics
Conclusions
• Two separate nucleation models combined• Formulated model for elastic deformations
![Page 70: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/70.jpg)
14-12-2009 34
Delft Institute of Applied Mathematics
Conclusions
• Two separate nucleation models combined• Formulated model for elastic deformations• Coupling between nucleation and
deformations
![Page 71: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/71.jpg)
14-12-2009 34
Delft Institute of Applied Mathematics
Conclusions
• Two separate nucleation models combined• Formulated model for elastic deformations• Coupling between nucleation and
deformations• Simulations show influence of deformations
on nucleation
![Page 72: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/72.jpg)
14-12-2009 35
Delft Institute of Applied Mathematics
Recommended future work
• Extension to multiple particle configurations
![Page 73: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/73.jpg)
14-12-2009 35
Delft Institute of Applied Mathematics
Recommended future work
• Extension to multiple particle configurations• Adaption to other alloys
![Page 74: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/74.jpg)
14-12-2009 35
Delft Institute of Applied Mathematics
Recommended future work
• Extension to multiple particle configurations• Adaption to other alloys• Improving numerical techniques
![Page 75: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/75.jpg)
14-12-2009 35
Delft Institute of Applied Mathematics
Recommended future work
• Extension to multiple particle configurations• Adaption to other alloys• Improving numerical techniques• Comparison with experimental data
![Page 76: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/76.jpg)
14-12-2009 35
Delft Institute of Applied Mathematics
Recommended future work
• Extension to multiple particle configurations• Adaption to other alloys• Improving numerical techniques• Comparison with experimental data• Including plastic deformations
![Page 77: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/77.jpg)
14-12-2009 35
Delft Institute of Applied Mathematics
Recommended future work
• Extension to multiple particle configurations• Adaption to other alloys• Improving numerical techniques• Comparison with experimental data• Including plastic deformations• Including homogeneous nucleation
![Page 78: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/78.jpg)
14-12-2009 35
Delft Institute of Applied Mathematics
Recommended future work
• Extension to multiple particle configurations• Adaption to other alloys• Improving numerical techniques• Comparison with experimental data• Including plastic deformations• Including homogeneous nucleation• Including grain prediction models
![Page 79: Mathematical modeling of particle nucleation and growth in ...homepage.tudelft.nl/d2b4e/numanal/ouden_presentation2.pdf · Mathematical modeling of particle nucleation and growth](https://reader034.vdocument.in/reader034/viewer/2022042417/5f32cfd82dadd749ff5f4df0/html5/thumbnails/79.jpg)
References
Chau, K.T. and Wei, X.X. . Finite solid circular cylinders sub-
jected to arbitrary surface load. Part I – analytic solution.
International Journal of Solids and Structures, 37(40):5707–
5732, 2000.
Jaeger, J.C. , Cook, N.G.W. , and Zimmerman, R.W. .
Fundamentals of rock mechanics. Wiley, Blackwell, 2007.
Kampmann, R. , Eckerlebe, H. , and Wagner, R. . In:
Materials Research Society Symposium Proceedings, vol-
ume 57, page 525. MRS, 1987.
Myhr, O. R. and Grong, Ø. . Modelling of non-isothermal trans-
formations in alloys containing a particle distribution. Acta
Materialia, 48(7):1605–1615, 2000.
Pal, R. . New models for effective Young’s modulus of partic-
ulate composites. Composites Part B: Engineering, 36(6-7):
513–523, 2005.
Robson, J.D. , Jones, M.J. , and Prangnell, P.B. . Extension of
the N-model to predict competing homogeneous and hetero-
geneous precipitation in Al-Sc alloys. Acta Materialia, 51(5):
1453–1468, 2003.
35-1