SAMR Phase Field Modeling (2007/08/14)
Results for uniform grid case: theta = 90
- Domain size: 200x200
- Node spacing: 0.025
- Time step size (dt): 0.000195313 (computed in terms of dx, etc.)
- Time integration scheme: explicit 4th-order Runge-Kutta
- Spatial discretization: all 2nd-order, effective Laplacian is 9pt
- Mobility: 3000
- Strength of anisotropy: delta = 0.04
- Number of symmetry planes: j = 6
- Dimensionless latent heat: K = 2
- Dimensionless thermal conductivity: kappa = 1
- Initial undercooling: T0 = 272
- Equilibrium temperature: Teq = 273
- Rotation of system: theta = 90 degrees
- Initial solid: circle of radius 0.05
- No heat flux at boundaries. Homogeneous Neumann boundary conditions for phase field.
NOTES
- I didn't put the results after the simulation hit the wall because my computational cell is larger than the one Chin Yi used.
- There seem to be grid anisotropy effects here. Growth is faster along the y-axis.
- Click on image to see enlarged figure.
|
Time = 0.0000 Area = 0.0075 |
|
|
Time = 0.00976562 Area = 0.0784663 |
|
|
Time = 0.0195312 Area = 0.185545 |
|
|
Time = 0.0390625 Area = 0.422308 |
|
|
Time = 0.078125 Area = 1.01125 |
|
|
Time = 0.107422 Area = 1.55483 |
|
|
Time = 0.15625 Area = 2.65573 |
|
|
Time = 0.205078 Area = 3.99858 |
|
|
Time = 0.302734 Area = 7.0818 |
|
Results for uniform grid case: theta = 20
- Domain size: 200x200
- Node spacing: 0.025
- Time step size (dt): 0.000195313 (computed in terms of dx, etc.)
- Time integration scheme: explicit 4th-order Runge-Kutta
- Spatial discretization: all 2nd-order, effective Laplacian is 9pt
- Mobility: 3000
- Strength of anisotropy: delta = 0.04
- Number of symmetry planes: j = 6
- Dimensionless latent heat: K = 2
- Dimensionless thermal conductivity: kappa = 1
- Initial undercooling: T0 = 272
- Equilibrium temperature: Teq = 273
- Rotation of system: theta = 20 degrees
- Initial solid: circle of radius 0.05
- No heat flux at boundaries. Homogeneous Neumann boundary conditions for phase field.
NOTES
- I didn't put the results after the simulation hit the wall because my computational cell is larger than the one Chin Yi used.
- Grid anisotropy effects seem to be less important for this case.
- Click on image to see enlarged figure.
|
Time = 0.0000 Area = 0.0075 |
|
|
Time = 0.00976562 Area = 0.076994 |
|
|
Time = 0.0195312 Area = 0.182464 |
|
|
Time = 0.0325625 Area = 0.410688 |
|
|
Time = 0.078125 Area = 0.961494 |
|
|
Time = 0.107422 Area = 1.46193 |
|
|
Time = 0.15625 Area = 2.46931 |
|
|
Time = 0.205078 Area = 3.68651 |
|
|
Time = 0.302734 Area = 6.74154 |
|
