Page 180 - Fundamentals of Computational Geoscience Numerical Methods and Algorithms

170 7 Simulating Thermal and Chemical Effects of Intruded Magma Solidification Problems

computational domain. Under the above temperature and pressure conditions of the
3
intruded magma, the density of the released volatile fluids is 374.7 kg/m (Haar
3
et al. 1984), while the density of the albite melt is 2700 kg/m . In addition, the mole
fraction of the released volatile fluids from the albite melt can be determined as
0.5306 from Eqs. (7.10) and (7.11), which results in a mass source of 1.433 × 10 –8
3
kg/(m × s) for the released volatile fluids in the computation.
Figure 7.7 shows the temperature distributions of the dike-like magma intrusion
problem at four different time instants. It is clear that with increasing time, the total
temperature localization area generated by the intruded magma becomes larger and

o
500 T ( C )
o t 1= . 8765 × 10 11 s
+ t = . 4 6673× 10 11 s
440 11
* t = . 9 3826 × 10 s
x t 1= . 8765 × 10 12 s
380

320

260

200
0 8.4 16.8 25.2 33.6 42.0
x (km)
(Along the cross-section of y = 3000 m)

o
T ( C )
620
o t 1= . 8765 × 10 11 s
500 + t = . 4 6673× 10 11 s
t = . 9 3826 × 10 11 s
*
x t 1= . 8765 × 10 12 s
380
260

140

20
0 2.0 4.0 6.0 8.0 10.0
y (km)
(Along the cross-section of x = 20,100 m)
Fig. 7.8 Temperature distributions of the dike-like magma intrusion problem along two different
cross-sections