| 8 | |
| 9 | The hard pulsating sphere has velocity variation. Its amplitude is 5km/s and its period is 1 year. |
| 10 | |
| 11 | The radius of the photosphere is also fixed. In reality, it must change but it is difficult to calculate. Within the photosphere, reduced gravity (80% of the original gravity) is being used. The reduced gravity has its physical stand - higher opacity inside the star but it is also because the resolution is poor. The simulation is resolving 2~3 order of magnitude density drop within 3 cells. |
| 12 | |
| 13 | Between the photosphere and the dust formation zone, it is the gas opacity ($2\times10^{-4}cm^{2}/g$). The dust formation radii is calculated by |
| 14 | |
| 15 | $R_d=\sqrt{\frac{\alpha L \exp{-\tau}}{4 \pi \sigma T_{sub}^{4}}}$ |