| 1 | RHEOPLAST 2-D |
|---|
| 2 | |
|---|
| 3 | Copyright (C) 2002 Adam Powell, David Dussault |
|---|
| 4 | |
|---|
| 5 | No, this is not a nose job, that's rhinoplasty. RheoPlast 2-D is for modeling |
|---|
| 6 | rheology of semi-solid metals, which are sort of plastic (at least at high |
|---|
| 7 | solid fraction they are, sort of), and is meant to sound like the semi-solid |
|---|
| 8 | process "rheocast". This package consists of three programs which use the |
|---|
| 9 | PETSc paralles solvers and data management tools with my special adaptations in |
|---|
| 10 | explicit.h, all using periodic boundary conditions. |
|---|
| 11 | |
|---|
| 12 | diffuse is simply a testbed for the fast explicit stuff, and uses the BLAS for |
|---|
| 13 | extreme speed. It can run in 2-D or 3-D, in 3-D it uses the Illuminator's |
|---|
| 14 | Geomview front-end to view the concentration field. |
|---|
| 15 | |
|---|
| 16 | cahnhill is the program used for the Cahn-Hilliard and fluid flow simulations |
|---|
| 17 | shown in the TMS Annual Meeting 2002, and the Physical Metallurgy Gordon |
|---|
| 18 | Research Conference 2002. It has some serious trouble with the pressure, and |
|---|
| 19 | so was not publication quality, and is included here only for "historical" |
|---|
| 20 | purposes. This is different from the others in that it uses symmetry boundary |
|---|
| 21 | conditions. See the comments therein for instructions on duplicating the |
|---|
| 22 | presented results. It is based on David Dussault's master's thesis code, and |
|---|
| 23 | includes the dussolve.c (Dussault's solver) and outputs.c (text viewer) files, |
|---|
| 24 | as well as the display.c program to visualize the whole thing with reflections |
|---|
| 25 | through the symmetry planes. |
|---|
| 26 | |
|---|
| 27 | rheoplast is the big deal. It does Cahn-Hilliard or vector-valued crystal |
|---|
| 28 | solidification/grain boundary phase field, only in 2-D for now but 3-D is |
|---|
| 29 | coming soon. It also does flow using the streamfunction-vorticity formulation, |
|---|
| 30 | thus avoiding the earlier pressure problem and giving the vorticity directly |
|---|
| 31 | for crystal rotation. It is quite cool. Like cahnhill, see comments in the |
|---|
| 32 | code to duplicate the published results. Pretty soon there will be a front end |
|---|
| 33 | based on glade. |
|---|
| 34 | |
|---|
| 35 | So there you have it. It's licensed under the GPL. Share and enjoy! |
|---|
