5. Numerical Results
5.1. Vaporization Test Problem
In this section, three test problems are solved to check the validity of the proposed algorithm and the high accuracy expected. The first example is a multi-moving boundary problem  . The thermo-physical properties are listed below in Table 1.
The solid material herein is of a low thermal conductivity and so the vaporization occurs before the moving boundary separating liquid and solid reaches the adiabatic boundary. The result due to the present algorithm is shown in Figure 3. In this figure, both vapor/liquid and liquid/solid moving boundaries are plotted both on the same figure. The liquid/solid moving boundary is in the left side, while the other one is on the right. From this figure, one can clearly determine the melting, vapor and the time at which the process ends. These times are determined and summarized in Table 2, at two different time steps. The absolute errors between the results due to the present method compared with the finite elements method are also presented in the same table. As we see, the absolute errors decrease by decreasing the time step. On the other hand by decreasing the time step, the number of iterations increases. The number of iterations are shown in Table 2, corresponding to the two time steps used. It is clear that the increase in the number of iterations is not so much but the accuracy improved to
Figure 2. Flow chart for THREEPHASE subroutine.
Figure 3. Moving boundaries location for vaporization test problem.
Table 1. Numerical data for vaporization test problem.
Table 2. Comparison between melting, vaporization and end time.
nearly the half.
5.2. Ablation Test Problem-1
A solid medium initially at uniform temperature, , the boundary exposed to two different cases of input heat flux, constant, and linear, respectively. The domain in the present problem is still fixed while appearing two moving interfaces, solid-liquid and vapor (gas)-liquid (Ablation surface). The results due to the present algorithm are shown in Figures 4-6, respectively and the results are compared with the results due to the source and sink method  . Figure 4 shows the movement of solid-liquid due to constant input heat flux, and the resulting ablation surface due to the same input heat flux is shown in Figure 5. The same results due to linear case are plotted on the same plot as shown in Figure 6. From the computations and figures, it is found that the solid-liquid interface has the same behavior in both constant and linear cases of input heat flux that is concave upward. In case of linear heat flux input this concavity becomes more apparent than the constant case. In the contrary, the vapor (gas)/liquid interface behaves concave downward but in linear case this concavity increases.
5.3. Ablation Test Problem-2
This problem is for a long enough solid mold initially at a uniform temperature. The surface exposed to two different high input heat flux, constant and linear and the vapor is removed as soon as it formed- ablation problem-. Two specific heat flux boundary condition are chosen in the present computations, namely, and the following values used in the calculation  (Table 3).
The ablation thickness due to the present method compared with the corresponding from the heat balance integral method is shown in Figure 7. It is clear that the ablation thickness is concave downward in both case of
Figure 4. Solid/liquid due to Constant heat flux-problem-1.
Figure 5. Ablated surface due to constant heat flux-problem-1.
Figure 6. Solid/liquid and Ablated surface due to linear heat flux-problem-1.
Figure 7. Ablation thickness for ablation test problem-2.
Table 3. Thermo-physical parameters.
input heat flux and at the same time, the ablation thickness in case of linear heat flux is higher than that for constant case. There is a good agreement between the two methods in both cases.
The importance of the present paper comes from its dealing with practical applications of wide range of our daily life. The boundary integral equation method is not so new but it is used as a mathematical tool due to its simplicity in use. Based on this method, a generalized numerical algorithm and computer code are developed to solve such applications. It is found from the computations that the developed algorithm and the code are so simple to handle and that an acceptable accuracy is obtained. Also by decreasing the time step, and a little bit increase of iterations, the absolute errors are decreased to the half nearly. Finally, the algorithm and subsequently the code can be easily modified to cover higher dimensional problems with acceptable accuracy, which can be improved by decreasing the time step, but on the other hand, stability should be achieved.
Cite this paper
Kawther K.Al-Swat,Said G.Ahmed, (2015) A New Iterative Method for Multi-Moving Boundary Problems Based Boundary Integral Method. Journal of Applied Mathematics and Physics,03,1126-1137. doi: 10.4236/jamp.2015.39140
: Fundamental solution
: Thermal diffusivity
: Input heat flux
: Liquid-solid interface
: Vapor-liquid interface
: Truncated long enough boundary length
: Latent heat