The Coal’s Layer Temperature Regime under Conditions of Heat Removal of Absorbed Microwave Energy by Radiation and Convection Simultaneously ()

V. A. Karelin^{1,2}, Vl. V. Salomatov^{1,2}, Vas. V. Salomatov^{1}

^{1}Institute of Thermal Physics, Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia.

^{2}Physics Department, Novosibirsk National Research State University, Novosibirsk, Russia.

**DOI: **10.4236/jcc.2018.61023
PDF HTML XML
556
Downloads
786
Views
Citations

Coal is one of the most popular sources of energy. However, it has a relatively low energy efficiency due to high humidity and a greater release of harmful substances during combustion. On the other hand, the coal reserve on earth is estimated at 500 years and the cost is relatively low. This causes the search for new ways of processing coal. One way to efficiently process coal, while reducing the humidity and content of harmful components is microwave treatment. The basic information for microwave exposure to coal is the temperature field. In this paper, an approximate-analytically nonlinear mathematical model for heating a flat coal mass is studied, provided that the absorbed microwave energy is removing by heat radiation and convection simultaneously.

Keywords

Coal, Microwave, Temperature Distribution, Stefan–Boltzmann Law, Convection

Share and Cite:

Karelin, V. , Salomatov, V. and Salomatov, V. (2018) The Coal’s Layer Temperature Regime under Conditions of Heat Removal of Absorbed Microwave Energy by Radiation and Convection Simultaneously. *Journal of Computer and Communications*, **6**, 229-238. doi: 10.4236/jcc.2018.61023.

1. Introduction

At present, there are many experimental confirmations [1] [2] of the prospects of using microwave treatment of coal material. Advantages over conventional processing methods (convective, infrared, conductive, etc.) include the uniformity and volume character of heating; the possibility of concentrating large capacities in certain places; obtaining the required temperature distribution; low inertia; absence of a heat-bearing medium, etc., which significantly increases the effectiveness of such an impact on energy coal. These specific features of microwave radiation make it possible to search for optimal heating regimes taking into account changes in the characteristics of coal during the process (humidity, content of other components, temperature, thermophysical and electromagnetic inhomogeneity, etc.) [3]. A comparison was made between the economics of drying by various methods and the microwave advantage was demonstrated [4]. A number of experiments were performed showing the possibility and efficiency of combining different types of drying with microwave drying [5].

Energy costs for evaporation under the influence of microwave radiation are reduced to an average of two times. It is shown that with evaporation of moisture, complete removal of nitrogen from wet coal samples occurs. In addition, there is a reduction of about 50% sulfur, 50% mercury, 30% ash, 50% chlorine [6].

Even a brief analysis of the application of microwave treatment of coal massifs demonstrates the prospects of its use and development in modern energy technologies, which will significantly improve the environmental safety and efficiency of burning coal fuel.

2. Problem of Microwave Heating of Plane Coal Massive

Coal is a dielectric material capable of absorbing electromagnetic radiation from the microwave range and thus forming an internal heat source. In this formulation, it was modeled according to Bouguer’s exponential law (Equation (1)). This internal heat source ensures warming up, ignition and burning out of the coal massif. The layer itself has two boundaries. At the top-the mode of heat removal using two mechanisms: radiation and convection (see Figure 1), on the bottom - thermal insulation. Radiation is most pronounced at high temperatures (~T^{4}), convection at lower temperatures (~T). However, in the general case, it is necessary to take into account two parallel heat fluxes.

The mathematical model contains the following assumptions:

1) Thermophysical characteristics of coal are constant and isotropic; 2) The initial temperature of the layer is constant along the coordinate; 3) Heat from the surface of the array is removed by heat radiation and convection simultaneously; 4) The non-isothermicity of a coal massif is taken into account.

Taking into account the above assumptions, the system of differential equations of the given problem will be written in the following form:

$\frac{\partial T(x,t)}{\partial t}=a\frac{{\partial}^{2}T(x,t)}{\partial {x}^{2}}+\frac{{q}_{{v}_{0}}}{c\rho}{e}^{-kx}$ (1)

$0\le t\le {t}_{s}$ , $0\le x\le l$

$T(x,0)={T}_{0}$ (2)

$-\lambda \frac{\partial T(l,t)}{\partial x}=\sigma \left[{T}^{4}(l,t)-{T}_{c}{}^{4}\right]+\alpha \left[T(l,t)-{T}_{c}\right]$ (3)

$\frac{\partial T(0,t)}{\partial x}=0$ (4)

Figure 1. The scheme of the problem of microwave heating with heat dissipation by convection and heat radiation simultaneously.

Here t is the current time, s; x―current Cartesian coordinate, m; t_{s}―the moment of ignition, s; l―thickness of the coal massif, m; a―coefficient of thermal diffusivity, m^{2}/s; T (x, t) is the current temperature, K; λ―coefficient of thermal conductivity, W/m∙K; c―heat capacity, kJ/kg∙K; α―Heat transfer coefficient at the upper surface, W/m^{2}∙K; ρ―density, kg/m^{3}; σ―Visible coefficient of thermal radiation, W/m^{2}・K^{4}; k is the absorption coefficient of microwave radiation, m^{−1}; T_{c}―ambient temperature, K.

In order to find a generalized solution, we turn to dimensionless variables, choosing as scales the quantities entering into the uniqueness conditions of the problem. To this end, we introduce new variables into problem (1)-(4):

$X=\frac{x}{l}$ ―dimensionless coordinate,

$Fo=\frac{t}{{l}^{2}/a}$ ―Fourier number,

$Po=\frac{{q}_{{v}_{0}}{l}^{2}}{\lambda {T}_{c}}$ ―number of Pomerantsev,

$Bu=kl$ ―Bouguer’s number,

$\theta (X,Fo)=\frac{T(x,t)}{{T}_{c}}$ ―dimensionless temperature,

$Sk=\frac{\sigma {T}_{c}{}^{3}l}{\lambda}$ ―Stark number,

$Bi=\frac{\alpha l}{\lambda}$ ―Bio number.

As a result, we have:

$\frac{\partial \theta (X,Fo)}{\partial Fo}=\frac{{\partial}^{2}\theta (X,Fo)}{\partial {X}^{2}}+Po\cdot {e}^{-BuX}$ (5)

$\theta (X,0)={\theta}_{0}$ (6)

$-\frac{\partial \theta (1,Fo)}{\partial X}=Sk\left[{\theta}^{4}(1,Fo)-1\right]+Bi\left[\theta (1,Fo)-1\right]\equiv Ki(Fo)$ (7)

$\frac{\partial \theta (0,Fo)}{\partial X}=0$ (8)

Problem (5)-(8) does not have a strict analytic solution, since the Stefan-Boltzmann law is nonlinear in (7). However, as shown in our paper [7], it is possible to construct sufficiently effective asymptotics for such problems. Next, we turn to the construction of such asymptotics for small and large values of Fo. To this end, we translate the system (5)-(8) from the region of the originals into the Laplace region of images.

The transformed system looks like this:

$s{\theta}_{L}(X,s)-{\theta}_{0}={{\theta}^{\u2033}}_{L}(X,s)+\frac{Po}{s}{e}^{-BuX}$ (9)

$-{{\theta}^{\prime}}_{L}(1,s)=K{i}_{L}(s)$ (10)

${{\theta}^{\prime}}_{L}(0,s)=0$ (11)

We find the solution of this system as the sum of the solution of a homogeneous problem and a solution with allowance for the inhomogeneous term:

${\theta}_{L}(X,s)={Y}_{1}(X,s)+{Y}_{2}(X,s)$

The solution of the homogeneous problem in the general form gives:

${Y}_{1}(X,s)=A\cdot {e}^{\sqrt{s}\cdot X}+B\cdot {e}^{-\sqrt{s}\cdot X}$

The nonhomogeneous part must be of the form:

${Y}_{2}(X,s)=C\cdot {e}^{-Bu\cdot X}+\frac{{\theta}_{0}}{s}$

Next we substitute the general form of the solution

${\theta}_{L}(X,s)=A\cdot {e}^{\sqrt{s}\cdot X}+B\cdot {e}^{-\sqrt{s}\cdot X}+C\cdot {e}^{-Bu\cdot X}+\frac{{\theta}_{0}}{s}$

for an inhomogeneous problem in Equations (9)-(11), we find the constants.

Then the general and particular solution for the Laplace-transformed temperature distribution can be written:

$\begin{array}{l}{\theta}_{L}(X,s)=\frac{Po\cdot Bu}{{s}^{\frac{3}{2}}\cdot (s-B{u}^{2})}\left[\left(1+\frac{{e}^{\sqrt{s}}-{e}^{-Bu}}{{e}^{-\sqrt{s}}-{e}^{\sqrt{s}}}\right)\cdot {e}^{\sqrt{s}\cdot X}+\frac{{e}^{\sqrt{s}}-{e}^{-Bu}}{{e}^{-\sqrt{s}}-{e}^{\sqrt{s}}}{e}^{-\sqrt{s}\cdot X}\right]\\ -K{i}_{L}(s)\cdot \Phi (s,X)+\frac{Po}{s(s-B{u}^{2})}e{}^{-BuX}+\frac{{\theta}_{0}}{s}\end{array}$ (12)

Here is Transmission function:

$\Phi (s,X)=\frac{ch(\sqrt{s}X)}{\sqrt{s}\cdot sh\sqrt{s}}$

Next, we turn to the search for asymptotic solutions for small times.

3. Small Fo (Large s)

We introduce the notation ${Y}_{1}(X,s)=-Ki(s)\cdot \Phi (s,X)$ ―the general solution.

$\begin{array}{l}{Y}_{2}(X,s)=\frac{Po\cdot Bu}{\sqrt{s}\cdot s\cdot (s-B{u}^{2})}\left[\left(1+\frac{{e}^{\sqrt{s}}-{e}^{-Bu}}{{e}^{-\sqrt{s}}-{e}^{\sqrt{s}}}\right)\cdot {e}^{\sqrt{s}\cdot X}+\frac{{e}^{\sqrt{s}}-{e}^{-Bu}}{{e}^{-\sqrt{s}}-{e}^{\sqrt{s}}}{e}^{-\sqrt{s}\cdot X}\right]+\\ +\frac{Po}{s(s-B{u}^{2})}{e}^{-BuX}+\frac{{\theta}_{0}}{s}\end{array}$

- A particular solution.

Then the system of basic Equations (5)-(8) after the Laplace transform looks as follows:

${{Y}^{\u2033}}_{1}(X,s)-s{Y}_{1}(X,s)=0$

$-{{Y}^{\prime}}_{1}(1,s)=Ki(\; s\; )$

${{Y}^{\prime}}_{1}(0,s)=0$

We represent the transfer function in the form of an expansion in the large parameter s:

$\u0424(s,X)=\frac{ch(\sqrt{s}X)}{\sqrt{s}\cdot sh\sqrt{s}}\approx {\Psi}_{1}(X,s)\frac{{e}^{-\sqrt{s}}}{\sqrt{s}}+{\Psi}_{2}(X,s)\frac{{e}^{-2\sqrt{s}}}{s}+\mathrm{...}$

As a result, the general solution in the images becomes:

${Y}_{1}(X,s)=-K{i}_{L}(s)\left[{\Psi}_{1}(X,s)\frac{{e}^{-\sqrt{s}}}{\sqrt{s}}+{\Psi}_{2}(X,s)\frac{{e}^{-2\sqrt{s}}}{s}+\mathrm{...}\right]$ (13)

Further we substitute the solution ${Y}_{1}(X,s)$ into the system (9)-(11) and equate the terms with the same degrees of expansion.

As a result, for ${\Psi}_{1}(X,s)$ we get

${{\Psi}^{\u2033}}_{1}(X,s)-s{\Psi}_{1}(X,s)=0$ (14)

From Equation (10) it follows that

$-K{i}_{L}(s)\left[{{\Psi}^{\prime}}_{1}(1,s)\frac{{e}^{-\sqrt{s}}}{\sqrt{s}}\right]=-K{i}_{L}(\; s\; )$

${{\Psi}^{\prime}}_{1}(1,s)=\frac{\sqrt{s}}{{e}^{-\sqrt{s}}}$ (15)

${{\Psi}^{\prime}}_{1}(0,s)=0$ (16)

The solution of (15) is represented in the form

${\Psi}_{1}(1,s)=A\cdot sh(\sqrt{s}X)+B\cdot ch(\sqrt{s}X)$ (17)

From the boundary condition (16) it follows that

$A=0$

As a result, we have:

$-K{i}_{L}(s)\left[{\Psi}_{1}(X,s)\frac{{e}^{-\sqrt{s}}}{\sqrt{s}}+\mathrm{...}\right]=B\cdot ch(\; s\; X\; )$

For X = 1, using (15), we find

$B=-\frac{Ki(s)}{\sqrt{s}\cdot sh(\sqrt{s})}$ (18)

Then

${\Psi}_{1}(X,s)=-\frac{Ki(s)\cdot ch(\sqrt{s}X)}{\sqrt{s}\cdot sh(\sqrt{s})}$ (19)

Further we substitute hyperbolic functions into (19):

${\Psi}_{1}(X,s)=-Ki(s)\cdot ({e}^{-\sqrt{s}(X+1)}+{e}^{\sqrt{s}(X-1)})$ (20)

As a result, the solution in images becomes

$\begin{array}{l}{\theta}_{L}(X,s)=\frac{Po\cdot Bu}{{s}^{\frac{3}{2}}\cdot (s-B{u}^{2})}\left[\left(1+\frac{{e}^{\sqrt{s}}-{e}^{-Bu}}{{e}^{-\sqrt{s}}-{e}^{\sqrt{s}}}\right)\cdot {e}^{\sqrt{s}\cdot X}+\frac{{e}^{\sqrt{s}}-{e}^{-Bu}}{{e}^{-\sqrt{s}}-{e}^{\sqrt{s}}}{e}^{-\sqrt{s}\cdot X}\right]\\ -\frac{K{i}_{L}(s)}{\sqrt{s}}\cdot \left({e}^{-\sqrt{s}(X+2)}+{e}^{\sqrt{s}(X-2)}\right)+\frac{Po}{s(s-B{u}^{2})}{e}^{-BuX}+\frac{{\theta}_{0}}{s}\end{array}$ (21)

Then we carry out the inverse Laplace transform. For large s it is permissible to assume that $s\gg Bu$ and

${L}^{-1}\left[\frac{Po}{{s}^{2}}\cdot {e}^{-BuX}\right]=Po{e}^{-BuX}Fo$ (22)

For large s, the expression

$\frac{Po\cdot Bu}{\sqrt{s}\cdot s\cdot (s-B{u}^{2})}\left[\left(1+\frac{{e}^{\sqrt{s}}-{e}^{-Bu}}{{e}^{-\sqrt{s}}-{e}^{\sqrt{s}}}\right)\cdot {e}^{\sqrt{s}\cdot X}+\frac{{e}^{\sqrt{s}}-{e}^{-Bu}}{{e}^{-\sqrt{s}}-{e}^{\sqrt{s}}}{e}^{-\sqrt{s}\cdot X}\right]$

goes into

$-\frac{Po\cdot Bu}{\sqrt{s}\cdot {s}^{2}}{e}^{-\sqrt{s}\cdot X}$

Let

$N(X,Fo)={L}^{-1}\left[-\frac{Po\cdot Bu}{\sqrt{s}\cdot {s}^{2}}{e}^{-\sqrt{s}\cdot X}\right]$

Then the final solution is written in convolution

$\begin{array}{l}\theta (X,Fo)=-{\displaystyle \underset{0}{\overset{Fo}{\int}}\frac{Ki(\eta )}{\sqrt{\pi (Fo-\eta )}}\cdot ({e}^{-{\frac{(X+2)}{4(Fo-\eta )}}^{2}}+{e}^{{\frac{(X-2)}{4(Fo-\eta )}}^{2}})d}\eta \\ +Po\cdot {e}^{-BuX}Fo+N(X,Fo)\end{array}$ (23)

To obtain an explicit calculation expression, we expand $Ki(\eta )$ in the Taylor series near $\eta \approx Fo$ : $Ki(\eta )=Ki(Fo)+(\eta -Fo)\cdot Ki(Fo)+\mathrm{...}$

Limited to the first term of the expansion, we have

$\begin{array}{l}\theta (X,Fo)=-Ki(Fo){\displaystyle \underset{0}{\overset{Fo}{\int}}\frac{1}{\sqrt{\pi (Fo-\eta )}}\cdot ({e}^{-{\frac{(X+2)}{4(Fo-\eta )}}^{2}}+{e}^{{\frac{(X-2)}{4(Fo-\eta )}}^{2}})d}\eta \\ +Po\cdot {e}^{-BuX}Fo+N(X,Fo)\end{array}$ (24)

After taking the integral, we get:

$\begin{array}{l}\theta (X,Fo)=\left(Sk\left({\theta}^{4}(1,Fo)-1\right)+Bi\left(\theta (1,Fo)-1\right)\right)\\ \cdot (\frac{X}{\sqrt{\pi}}-\frac{X+2}{\sqrt{\pi}\cdot \frac{X+2}{\sqrt{Fo}}}\cdot \left[\frac{\sqrt{\pi}(X+2)}{2\sqrt{Fo}}erf\left(\frac{X+2}{2\sqrt{Fo}}\right)+{e}^{-{\left(\frac{(X+2)}{2\sqrt{Fo}}\right)}^{2}}\right]\\ -\frac{X-2}{\sqrt{\pi}\cdot \frac{X-2}{\sqrt{Fo}}}\cdot \left[\frac{\sqrt{\pi}(X-2)}{2\sqrt{Fo}}erf\left(\frac{X-2}{2\sqrt{Fo}}\right)+{e}^{-{\left(\frac{(X-2)}{2\sqrt{Fo}}\right)}^{2}}\right])\\ +Po\cdot {e}^{-BuX}Fo+{\theta}_{0}+N(X,Fo)\end{array}$ (25)

The boundary temperature is found by substituting X = 1 in (30):

$\begin{array}{l}\theta (1,Fo)=\left(Sk\left({\theta}^{4}(1,Fo)-1\right)+Bi\left(\theta (1,Fo)-1\right)\right)\\ \cdot \left(\frac{1}{\sqrt{\pi}}-\frac{3}{2}erf\left(\frac{3}{2\sqrt{Fo}}\right)+\frac{1}{2}erf\left(-\frac{1}{2\sqrt{Fo}}\right)-\frac{\sqrt{Fo}}{\sqrt{\pi}}\left[{e}^{-{\left(\frac{3}{2\sqrt{Fo}}\right)}^{2}}+{e}^{-{\left(\frac{1}{2\sqrt{Fo}}\right)}^{2}}\right]\right)\\ +Po\cdot {e}^{-Bu}Fo+{\theta}_{0}+N(1,Fo)\end{array}$

As a result, the surface temperature is calculated by the following expression:

$\theta \left(1,Fo\right)=\frac{\sqrt{b+c}\pm \sqrt{-b-c+2\left(1-Bi\right){\left(Sk\sqrt{b+c}\right)}^{-1}}}{2}+a$ (26)

Here are new variables:

$a=+Po\cdot {e}^{-Bu}Fo+{\theta}_{0}+N(1,Fo)$

$b=\frac{4\cdot {2}^{\frac{1}{3}}a}{{\left(27\cdot {(Bi\cdot d-1)}^{2}Sk\cdot d+\sqrt{729{(Bi\cdot d-1)}^{4}{(Sk\cdot d)}^{2}-6912\cdot {a}^{3}{(Sk\cdot d)}^{3}}\right)}^{\frac{1}{3}}}$

$c=\frac{{\left(27\cdot {(Bi\cdot d-1)}^{2}Sk\cdot d+\sqrt{729{(Bi\cdot d-1)}^{4}{(Sk\cdot d)}^{2}-6912\cdot {a}^{3}{(Sk\cdot d)}^{3}}\right)}^{\frac{1}{3}}}{3\cdot {2}^{\frac{1}{3}}Sk\cdot d}$

$d=\left(\frac{1}{\sqrt{\pi}}-\frac{3}{2}erf\left(\frac{3}{2\sqrt{Fo}}\right)+\frac{1}{2}erf\left(-\frac{1}{2\sqrt{Fo}}\right)-\frac{\sqrt{Fo}}{\sqrt{\pi}}\left[{e}^{-{\left(\frac{3}{2\sqrt{Fo}}\right)}^{2}}+{e}^{-{\left(\frac{1}{2\sqrt{Fo}}\right)}^{2}}\right]\right)$

4. Large Fo (Small s)

We take the inverse transformation from the term with the source as we do on previous stage, but now we calculating expressions with assumption of small s. Let’s set new variable:

$M(X,Fo)={L}^{-1}\left[\frac{Po\cdot Bu}{{s}^{2}B{u}^{2}}\left(1-{e}^{-Bu}\right)\right]=\frac{\left(1-{e}^{-Bu}\right)Fo\cdot Po}{Bu}$

We shall construct the asymptotics for large times-for this purpose we expand the transfer function in a series with respect to the small parameter s:

$\u0424(s,X)={\phi}_{0}(s,X)+s{\phi}_{1}(s,X)+{s}^{2}{\phi}_{2}(s,X)+\cdots $ (27)

Substituting (27) into the basic Equation (9) and equating the terms with the same powers of s, we obtain a system for determining ${\phi}_{0}$ , ${\phi}_{1}$ , ${\phi}_{2}$

The system that is used for this:

${{Y}^{\u2033}}_{1}(X,s)-s\cdot {Y}_{1}(X,s)=0$

$-{{Y}^{\prime}}_{1}(1,s)=K{i}_{L}(\; s\; )$

${{Y}^{\prime}}_{1}(0,s)=0$

As a result, we get

${{\phi}^{\u2033}}_{0}(X,s)=0,{{\phi}^{\u2033}}_{1}(X,s)={\phi}_{0},{{\phi}^{\u2033}}_{2}(X,s)={\phi}_{1}$ (28)

Each equation in (28) requires 2 boundary conditions for finding two constants.

The boundary condition:

${{\phi}^{\prime}}_{0}|{}_{X=0}=0,{{\phi}^{\prime}}_{1}|{}_{X=0}=0,{{\phi}^{\prime}}_{2}|{}_{X=0}=0$ (29)

The second constant is found from the integral relations:

$\underset{0}{\overset{1}{\int}}{{\phi}^{\u2033}}_{0}}dX=0,{\displaystyle \underset{0}{\overset{1}{\int}}{{\phi}^{\u2033}}_{1}}dX={\displaystyle \underset{0}{\overset{1}{\int}}{\phi}_{0}}dX,{\displaystyle \underset{0}{\overset{1}{\int}}{{\phi}^{\u2033}}_{2}}dX={\displaystyle \underset{0}{\overset{1}{\int}}{\phi}_{1}}dX$ (30)

We solve the system, finding for the asymptotics in the originals with allowance for the two terms of the series (27). As a result, we have

$\begin{array}{l}{\theta}_{L}(X,Fo)\approx -{\displaystyle \underset{0}{\overset{Fo}{\int}}Ki(\eta )d\eta}+Ki(Fo)\frac{1-3{X}^{2}}{6}\\ +Po\cdot Fo\cdot {e}^{-BuX}+{\theta}_{0}+M(X,Fo)\end{array}$ (31)

The unknown surface temperature $\theta (1,Fo)$ is found from the Voltaire integral equation of the second kind, while denoting ${\theta}_{w}(Fo)=\theta (1,Fo)$ :

$\begin{array}{l}{\theta}_{w}(Fo)\approx -{\displaystyle \underset{0}{\overset{Fo}{\int}}\left[Bi({\theta}_{w}(\eta )-1)+Sk({\theta}_{w}{}^{4}(\eta )-1)\right]d\eta}\\ -\frac{1}{3}\left[Bi({\theta}_{w}(Fo)-1)+Sk({\theta}_{w}{}^{4}(Fo)-1)\right]+Po\cdot Fo\cdot {e}^{-Bu}+{\theta}_{0}+M(1,Fo)\end{array}$ (32)

The solution of (32) is

$3Sk\cdot Fo={F}_{1}\left[{\theta}_{w}{}^{*},\frac{Bi}{Sk}\right]-{F}_{1}\left[{\theta}_{w},\frac{Bi}{Sk}\right]$ (33)

where the function F_{1} has the form

$\begin{array}{l}{F}_{1}[{\theta}_{w}]=-\frac{A}{4\sqrt{2}\cdot {B}^{\frac{3}{4}}}(2\cdot \mathrm{arctan}\left[1-\frac{\sqrt{2}\cdot {\theta}_{w}}{{B}^{\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$4$}\right.}}\right]-2\cdot \mathrm{arctan}\left[1+\frac{\sqrt{2}\cdot {\theta}_{w}}{{B}^{\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$4$}\right.}}\right]\\ +\mathrm{ln}\left(\frac{{\theta}_{w}{}^{2}-\sqrt{2}\cdot {\theta}_{w}\cdot {B}^{\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$4$}\right.}+\sqrt{B}}{{\theta}_{w}{}^{2}+\sqrt{2}\cdot {\theta}_{w}\cdot {B}^{\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$4$}\right.}+\sqrt{B}}\right))+Sk\cdot \mathrm{ln}\left({\theta}_{w}{}^{4}+B\right)\end{array}$

Here ${\theta}_{w}{}^{*}$ is the initial value in the integral Equation (32). We find it substituting Fo = 0:

${\theta}_{w}{}^{*}=-\frac{1}{3}\left[Bi({\theta}_{w}{}^{*}-1)+Sk({\theta}_{w}{{}^{*}}^{4}-1)\right]+{\theta}_{0}$ (34)

The solution of (34) is:

${\theta}_{w}{}^{*}=\frac{1}{2}\sqrt{-A+B}-\frac{1}{2}\sqrt{A-B+C}$ (35)

where

$\u0410=\frac{4{\left(\frac{2}{3}\right)}^{\frac{1}{3}}\left(Bi+Sk+3{\theta}_{0}\right)}{{\left(9{\left(3+Bi\right)}^{2}\cdot Sk+\sqrt{3\left(27{\left(3+Bi\right)}^{4}\cdot S{k}^{2}+256\cdot S{k}^{3}{\left(Bi+Sk+3{\theta}_{0}\right)}^{3}\right)}\right)}^{\frac{1}{3}}}$

$B=\frac{{\left(9{\left(3+Bi\right)}^{2}\cdot Sk+\sqrt{3\left(27{\left(3+Bi\right)}^{4}\cdot S{k}^{2}+256\cdot S{k}^{3}{\left(Bi+Sk+3{\theta}_{0}\right)}^{3}\right)}\right)}^{\frac{1}{3}}}{{2}^{\frac{1}{3}}\cdot {3}^{\frac{2}{3}}\cdot Sk}$

$C=\frac{2\cdot \left(3+Bi\right)}{Sk\cdot \sqrt{-A+B}}$

Knowing the surface temperature, you can calculate the required temperature:

$\begin{array}{l}{\vartheta}_{L}(X,Fo)\approx {\theta}_{w}(Fo)-(Bi\cdot ({\theta}_{w}(Fo)-1)+Sk({\theta}_{w}{}^{4}(Fo)-1))\cdot \frac{1-3X}{6}\mathrm{...}\\ +Po\cdot Fo\cdot {e}^{-BuX}\end{array}$ (36)

A typical temperature distribution obtained from this solution is demonstrated on Figure 2.

The stationary solution for large times is found from the basic Equation (5) and looks as follows:

$\theta \left(X,\infty \right)=\frac{Po}{B{u}^{2}}{e}^{-BuX}$ (37)

5. Conclusions

This article is a construction of a nonlinear model of heating a coal layer by microwave energy. At the same time internal heat sources are generated inside, which are modeled according to Bouguer’s law. On the upper surface, a heat release is prescribed according to the law of heat radiation and convection. Since they are decisive in the process of heating the coal to the ignition temperature. The lower boundary is considered to be insulated in the problem. The task itself to simplify the recording of the solution has been translated into dimensionless variables. The problem was analyzed using the Laplace transform. Because of the Stefan-Boltzmann law, the problem is nonlinear and, generally speaking, does

Figure 2. Temperature distribution over the depth of the layer.

not have a rigorous analytical solution. To do this, we searched for asymptotic solutions for large and small times.

The resulting temperature fields will be used in future to calculate thermal stresses, find the moment of ignition, optimum conditions, conduct microwave processing, taking into account ecology, energy saving, safety, etc.

Funding

The reported study was funded by Russian Foundation for Basic Research (RFBR) according to the research project No. 17-08-00752.

Conflicts of Interest

The authors declare no conflicts of interest.

[1] | Salomatov, Vl.V., Sladkov, S.O., Pashchenko, S.E. and Salomatov, Vas.V. (2013) Phenomenological Models of Microwave Heating of a Flat Coal Mass with Release of Absorbed Heat by the Convection Law. Journal of Engineering Thermophysics, 22, 134-142. https://doi.org/10.1134/S1810232813020057 |

[2] | Salomatov, V.V., Sladkov, S.O. and Pashchenko, S.E. (2012) Microwave Technologies in Coal Power Engineering. IFH, 85. 535-549. |

[3] | Pokras, S.M., Kolesnichenko, Yu.F., Kuryshev, V.A. and Vasiliev, K.B. (1991) On the Possibility of Using Microwave Energy to Solve Problems Arising from the Combustion of Solid Fuels. Proceedings of MRTI, Radiophysics Series, 266-273. |

[4] | Salomatov, Vl.V., Sladkov, S.O., Paschenko, S.E. and Salomatov, You.V. (2012) Application of Microwave in Coal Heat Power Engineering. Collection of Scientific Works of III Vseros. Scientific and Practical Conference with International Participation “Thermophysical Fundamentals of Energy Technologies”, Publishing House of Tomsk University, Tomsk, 8-23. |

[5] | Khaydurova, A.A. (2010) Microwave Drying of Brown Coals and Increase of Their Technological Characteristics. Izv. Universities. Problems of Energy. Publishing House KGEU, Kazan, 31-35. |

[6] | Patent USA (2007) 0295.590 Methods and Systems for Improving the Properties of Solid Fuels. |

[7] | Salomatov, Vl.V., Karelin, V.A. and Salomatov, Vas.V. (2016) Mathematical Models of Microwave Heating of a Coal Mass with Release of Absorbed Energy. Journal of Engineering Thermophysics, 25. |

Journals Menu

Copyright © 2021 by authors and Scientific Research Publishing Inc.

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.