Open access peer-reviewed chapter

# Modeling of the Two-Dimensional Thawing of Logs in an Air Environment

By Nencho Deliiski, Ladislav Dzurenda and Natalia Tumbarkova

Submitted: March 12th 2020Reviewed: June 12th 2020Published: July 9th 2020

DOI: 10.5772/intechopen.93177

## Abstract

A two-dimensional mathematical model has been created, solved, and verified for the transient nonlinear heat conduction in logs during their thawing in an air environment. For the numerical solution of the model, an explicit form of the finite-difference method in the computing medium of Visual FORTRAN Professional has been used. The chapter presents solutions of the model and its validation towards own experimental studies. During the validation of the model, the inverse task of the heat transfer has been solved for the determination of the logs’ heat transfer coefficients in radial and longitudinal directions. This task has been solved also in regard to the logs’ surface temperature, which depends on the mentioned coefficients. The results from the experimental and simulative investigation of 2D nonstationary temperature distribution in the longitudinal section of poplar logs with a diameter of 0.24 m, length of 0.48 m, and an initial temperature of approximately –30°C during their many hours thawing in an air environment at room temperature are presented, visualized, and analyzed.

### Keywords

• heat conduction
• modeling
• logs
• thawing
• heat transfer coefficients
• surface temperature

## 1. Introduction

The duration and the energy consumption of the thermal treatment of frozen logs aimed at their thawing and plasticizing for the production of veneer in winter are very high [1, 2, 3, 4, 5, 6, 7, 8, 9]. For example, thawing and plasticizing of poplar and pine logs with an initial temperature of –10°C and moisture content of 0.6 kg·kg−1 about 53 kWh·m−3 and 64 kWh·m−3 thermal energy, respectively, are needed .

In the specialized literature, there are few reports about the temperature fields subjected to thawing in agitated water or steam frozen logs [7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21], and there is very scarce information about research of the temperature distribution in frozen logs during their thawing in an air environment given by the authors only .

The computation of the temperature field in logs during their thawing in water or steam is carried out using mathematical models, which solve the so-called direct task of the heat transfer. This is the task when all variables in the model are known, and this allows computing the temperature field in the body [23, 24].

The computation of the temperature field in logs during their thawing in an air environment requires solving of the so-called inverse task of the heat transfer. This is the task when the model of the studied object and the experimentally obtained temperature field in it are known, but one or more variables in the model need to be determined during the solving and validation of the model .

The results from investigations of the temperature change subjected to thawing frozen logs only at conductive boundary conditions (i.e., at prescribed surface temperature) have been reported [2, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21].

The modeling and the multiparameter study of the thawing process of logs in air environment is of considerable scientific and practical interest. For example, as a result of such a study, it is possible to determine the real initial temperature of logs depending on their dimensions, wood species, moisture content, and the temperature of the air near the logs during their many days staying in an open warehouse before the thermal treatment in the production of veneer. The information about the real value of that immeasurable parameter can be used for scientifically based computing of the optimal, energy saving regimes for thermal treatment of each specific batch of logs.

This chapter presents the creation, numerical solving and validation of a two-dimensional nonlinear mathematical model of the transient heat conduction in frozen logs during their thawing at convective boundary conditions in an air environment. A validation of the models towards own experimentally determined 2D temperature distribution in poplar logs with a diameter of 0.24 m, length of 0.48 m, initial temperature of approximately –30°C, and moisture content above the hygroscopic range during their 70 h thawing at room temperature has been carried out.

During the validation of the model, the inverse task has been solved for the determination of the unknown logs’ heat transfer coefficients in radial and longitudinal directions. This task has been solved also in regard to the logs’ surface temperature, which depends on the mentioned coefficients.

## 2. Mechanism of 2D heat distribution in logs during thawing

### 2.1 Mathematical model of the 2D temperature distribution in frozen logs during their thawing in an air environment

In  the following common form of a model, which describes the 2D nonstationary temperature distribution subjected to thawing frozen logs in an air environment, has been suggested:

cwe1,2,3ρwTrzττ=λwr2Trzτr2+1r.Trzτr+λwrTTrzτr2++λwp2Trzτz2+λwpTTrzτz2E1

with an initial condition

Trz0=Tw0avgE2

and boundary conditions for convective heat transfer:

• Along the radial coordinate r on the logs’ frontal surface during thawing

Tr0τr=αwpr0τλwpr0τTr0τTmτE3

• Along the longitudinal coordinate z on the logs’ cylindrical surface during thawing

T0zτz=αwr0zτλwr0zτT0zτTmτE4

In  solutions of Eq. (1) only at conductive boundary conditions for the case of autoclave steaming of logs aimed at their plasticizing in the production of veneer have been realized and graphically presented.

An approach for solving Eq. (1) at much more complicated convective boundary conditions and verification of model (1) to (4) is considered below.

### 2.2 Mathematical description of the thermophysical characteristics of logs

In Figure 1 the three temperature ranges are presented, at which the process of the logs’ thawing above the hygroscopic range is carried out, i.e., when u >ufsp. Figure 1.Temperature ranges of the logs’ thawing process at u > u fsp and thermophysical characteristics of the wood and of the frozen and nonfrozen bound and free water in it.

There thermophysical characteristics of the logs and of both the frozen and nonfrozen free and bound water in them during the separate temperature ranges are also shown. The information on these characteristics is very important for the solving of the model given above.

The mathematical descriptions of the thermal conductivities of nonfrozen wood, λwnfr, and frozen wood, λwfr, and also of the specific heat capacities of nonfrozen wood, cwnfr, frozen wood, cwfr, and frozen free and bound water in the wood, cfwand cbwm, have been suggested in [8, 9, 18] using the experimentally determined data in the dissertations by Kanter  and Chudinov  for their change as a function of t and u. According to the suggested in [8, 9, 18] approach, the wood thermal conductivity during thawing of logs with moisture content u above the hygroscopic range can be calculated with the help of the following equations for λwTuρbufsp:

λw=λw0γ1+βT273.15,E5
v=0.12840.013u.E7

The coefficients γ and β in Eq. (5) are calculated using the next equations:

• For nonfrozen wood at u>ufsp272.15and 272.15K<T423.15K

γ=1.0,E8
β=3.65579ρb0.124103.E9

• For frozen wood at u>ufsp272.15and 213.15KT272.15K

γ=1+0.341.15uufsp,E10
β=0.002uufsp0.0038579ρb0.124.E11

The fiber saturation points of the wood species, u fsp and ufsp272.15, are calculated according to following Eqs. :

ufsp=ufsp293.150.001T293.15,E12

and consequently

ufsp272.15=ufsp293.15+0.021,E13

where ufsp293.15is the standardized fiber saturation point of the wood at T = 293.15 K (i.e., at t = 20°C), kg·kg−1, and ufsp272.15is the fiber saturation point at T = 272.15 K (i.e., at t = –1°C), kg·kg−1. At t = –1°C, the melting of the frozen bound water in the wood is fully completed, and the melting of the free water in the wood starts [22, 26].

The effective specific heat capacities of the logs during the pointed three ranges of the thawing process, cwe1,2,3, which participate in Eq. (1), are equal to the following:

First range:c1=cwfr+cbwm,E14
Second range:cwe2=cwnfr+cfw,E15
Third range:cwe3=cwnfr.E16

According to the suggested in [8, 9, 22] mathematical description, the effective specific heat capacities of the logs during their thawing can be calculated with the help of the following equations for cwe1,2,3Tuufsp:

cwe1,2,3=1.06+0.04u+0.00075T272.15ufsp272.15526+2.95T+0.0022T2+2261u+1976ufsp272.151+u++1.8938104ufsp272.150.12exp0.0567T272.151+u@u>ufsp272.15&213.15KT272.15K1strange2862u+5551+u+5.42u+2.951+uT+0.00361+uT2++3.34105uufsp272.151+u@u>ufsp272.15&272.15K<T<273.15K2ndrange2862u+5551+u+5.42u+2.951+uT+0.00361+uT2@u>ufsp272.15&273.15KT413.15K3rdrangeE17

The wood density, ρ w, which participates in Eq. (1), is determined above the hygroscopic range according to the following Equation [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 27]:

ρw=ρb1+u.E18

### 2.3 Mathematical description of the heat transfer coefficients of logs

For solving of the 2D mathematical model given above, it is needed to have values for the heat transfer coefficients of the logs in radial and longitudinal directions, α wr and α wp, respectively, which participate in Eqs. (3) and (4).

As it was mentioned in the Introduction, the values of α wr and α wp can be computed by solving the inverse task of the heat transfer between the logs and surrounding air environment.

The calculation of α wr and α wp can be carried out with the help of the following equations of the similarity theory, which are valid for the cases of heating of horizontally situated cylindrical bodies in conditions of free air convection :

αwr=NurλaL,E19
Nur=fGrrPra0.25PraPrs0.25,E21
Nup=fGrpPra0.25PraPrs0.25,E22
Grr=gβaL3wa2TmτTsτ,E23
Pra=waTaaaTa,E25
Prs=waTsaaTs,E26

For the usage of Eqs. (19)(26), it is needed to have a mathematical description of the thermophysical characteristics of the air, λ, β, w, and а, depending on T and φ. The temperature of the air near the logs subjected to thawing during our experiments described below changes in the range from 243.15 to 303.15 K (i.е., from –30°С to 30°С), and φ changes from 40–100% (see Figures 2 and 3). Figure 2.Experimentally determined change in tm, φm, and t in four points of the studied poplar log P1 during its 70 h thawing. Figure 3.Experimentally determined change in tm, φm, and t in four points of the studied poplar log P2 during its 70 h thawing.

For the calculation of λa, β a, w a, and а a, the temperature of the air, Т a, must be used, but for the calculation of w s and а s in Eq. (26), the temperature of the surface of the logs, T s, has to be used.

In the accessible specialized sources, we did not find suitable mathematical descriptions of λ, β, w, and а of the air, depending on T and φ, which could be applied for the precise determination of α wr and α wp according to Eqs. (19)(26).

Our further study has shown that for solving the inverse task of the heat transfer between the logs and surrounding air, i.e., for the calculation of the heat transfer coefficients of the logs, which participate in the boundary conditions (3) and (4) of the model, the following equations are suitable :

• In the radial direction on the cylindrical surface of the logs

αwr=1.123T0zτTmτxE27

• In the longitudinal direction on the frontal surface of the logs

αwp=2.56Tr0τTmτxE28

where x is an exponent, whose values are determined during the solving and validation of the model through the minimization of the root-square-mean error (RSME) between the calculated model and experimentally obtained results about the change of the temperature fields subjected to thawing logs.

## 3. Experimental research of the thawing process of logs

### 3.1 Experimental research of the 2D temperature distribution in poplar logs during their thawing

For solving the inverse task of the heat transfer aimed at validation of the suggested above mathematical model, it is necessary to have experimentally obtained data about the 2D temperature distribution in logs during their thawing. That is why we realized such experiments using poplar (Populus nigra L.) logs with D = 0.24 m, L = 0.48 m, and u > u fsp .

In Figure 4 the coordinates of four representative points of the logs, in which the 2D change in the temperature was measured and registered during the logs’ thawing, are shown. Figure 4.Radial (left) and longitudinal (right) coordinates of four characteristic points for the measurement of the temperature in logs subjected to thawing.

For the freezing of the logs before their thawing, a horizontal freezer was used with adjustable temperature range from –1 to –30°C. Sensors Pt100 with long metal casings were positioned in the drilled four holes of the logs. After 50 h separately freezing each logs, the freezer was switched off. Then its lid was opened, and 70 h thawing of the log at room temperature was carried out.

The automatic measurement and record of t m, φ m, and t in the representative points of the logs during the experiments was accomplished by Data Logger type HygroLog NT3 produced by the Swiss firm ROTRONIC AG.

In Figures 2 and 3, the change in the temperature of the processing air medium, t m, and in its humidity, φ m, and also in the temperature in four representative points of two poplar logs, named below as P1 and P2, respectively, during their separate 70 h thawing is presented.

All curves of the experimentally obtained data on these figures are drawn using the licensed software HW4 of the Data Logger. The left coordinate axis on the figures is graduated at % of φ m, and the right one is graduated at °C of t.

### 3.2 Mathematical description of the air medium temperature during logs’ thawing

The change shown in Figures 2 and 3 air medium temperature T m during the logs’ thawing with correlation 0.98 and root-square-mean error, σ < 1.5°C, has been approximated with the help of the software package Table Curve 2D by the following equation:

Тm=a+cτ0.51+bτ0.5,E29

whose coefficients are equal to:

• For log P1: а = 293.3637194, b = −0.00236425, c = −0.69281743.

• For log P2: a = 299.2738855, b = −0.00245303, c = −0.73047119.

and τ is the sum of the time of logs’ freezing, equal to 50 h = 180,000 s, and the current time of the subsequent thawing of the logs, s.

Eq. (29) was used for solving Eqs. (3) and (4) of the model.

## 4. Numerical solution of the mathematical model of the logs’ thawing process

The mathematical descriptions of the thermophysical characteristics of the logs and also of T m considered above were introduced in the mathematical model (1) to (4). An explicit form of the finite-difference method was used for solving of the model without any simplifications [7, 8].

### 4.1 Presentation of the model in a form suitable for programming

#### 4.1.1 Presentation of Eq. (1) of the model

The presentation of Eq. (1) of the model suitable for programming discrete analogue has been carried out using the given in Figures 5 and 6 coordinate system. These figures show the positioning of the knots of the calculation mesh and four representative points, in which the nonstationary 2D distribution of the temperature in the longitudinal section subjected to thawing log has been calculated. The mesh has been built on ¼ of the longitudinal section of the log due to the fact that this ¼ is mirrored symmetrical towards the remaining ¾ of the same section. Figure 5.Positioning of the knots of 2D calculation mesh on ¼ of longitudinal section of a log subjected to thawing (left) and calculation mesh for solving of the model (right). Figure 6.Calculation mesh and representative points T1, T2, T3, and T4 on ¼ of the longitudinal section subjected to thawing log.

Taking into consideration Eqs. (5) and (6), it can be written that

λw=λw0r,pγ1+βT273.15E30
λw0r=Kwrv0.165+1.39+3.8u3.3107ρb2+1.015103ρbE31
λw0p=Kwpv0.165+1.39+3.8u3.3107ρb2+1.015103ρbE32

and using the coefficient

Kwp/wr=KwpKwrE33

The discrete finite-difference analogue of the left-hand part of Eq. (1), which is suitable for programming in FORTRAN, has the following form [7, 30]:

cwe1,2,3ρwTrzττcwe1,2,3nρwTi,kn+1Ti,knΔτ.E34

Taking into account Eqs. (30), (31), and (32), the discrete analogue of the right-hand part of Eq. (1) has the following form:

λwr2Trzτr2+1r.Trzτr+λwrTTrzτr2++λwp2Trzτz2+λwpTTrzτz2=λw0rγ1+βTi,kn273.15Ti1,kn2Ti,kn+Ti+1,knΔr2+1i1ΔrTi1,knTi,knΔr+βTi1,knTi,knΔr22++λw0pγ1+βTi,kn273.15Ti,k1n2Ti,kn+Ti,k+1nΔz2+βTi,k1nTi,knΔz22.E35

After alignment of Eq. (34) with Eq. (35) and taking into account Eq. (33), at Δz = Δr, it is obtained that Eq. (1) is transformed into the following system of algebraic equations:

Ti,kn+1=Ti,kn+λw0rγΔτcwe1,2,3nρwΔr2..1+βTi,kn273.15.Ti1,kn+Ti+1,kn+Kwp/wrTi,k1n+Ti,k+1n2+2Kwp/wrTi,kn+1i1Ti1,knTi,kn++βTi1,knTi,kn2+Kwp/wrTi,k1nTi,kn2.E36

The term 1rin the right-hand part of Eq. (1) is represented as 1i1Δrin Eq. (35). According to the requirements of FORTRAN [7, 30], the knots of the calculation mesh, which are situated on the log’s surfaces, are denoted by numbers i = 1 and k = 1 along the coordinate axes r and z, respectively.

The temperature in these surface knots is calculated with the help of Eqs. (38) and (42) given below. Since Eq. (36) calculates the temperature in the knots, which are located inside the logs, i.e., in the knots with i ≥ 2 and k ≥ 2, the denominator of the term 1i1in this equation is always greater than zero.

It can be noted that the effective specific heat capacities of the log during the pointed above three ranges of its thawing process (see Figure 1), cwe1, cwe2, and cwe3, which are unitedly represented as cwe1,2,3in Eq. (36), are computed according to Eq. (17) separately for each knot of the calculation mesh.

#### 4.1.2 Presentation of the equation of the initial condition in the model

The initial condition (2) of the model of logs’ thawing process obtains the following discrete finite-difference form:

Тi,k0=Tw0avgE37

where T w0-avg is the experimentally determined average mass temperature of the log at the beginning of the thawing process, K.

#### 4.1.3 Presentation of the equation of boundary condition in the model along the radial coordinate

The boundary condition (3) of the logs’ thawing process obtains the following final form, suitable for programming in FORTRAN:

Тi,1n+1=Ti,2n+Gi,1nTmn+11+Gi,1n.E38

The variable Gi,1nin Eq. (38) is equal to

Gi,1n=Δrαwpnλw0pγ1+βTi,1n273.15,E39

where according to Eqs. (28) and (29)

αwpn=2.56Ti,1nTmnх,E40

where

Тmn=a+cnΔτ0.51+bnΔτ0.5.E41

#### 4.1.4. Presentation of the equation of boundary condition in the model along the longitudinal coordinate

Analogously, the boundary condition (4) of the logs’ thawing process obtains the following final form, suitable for programming in FORTRAN:

Т1,kn+1=T2,kn+G1,knTmn+11+G1,kn.E42

The variable G1,knin Eq. (42) is equal to

G1,kn=Δrαwrnλw0rγ1+βT1,kn273.15,E43

where according to Eq. (27)

αwrn=1.123T1,knTmnхE44

and Тmnis calculated according to Eq. (41).

### 4.2 Input data for solving of the model

The numerical solving and verification of the model (1) to (4) has been realized in the calculation environment of Visual FORTRAN Professional.

Using own software package in that environment, computations were carried out for the determination of the 2D nonstationary change of t in the representative points of the logs P1 and P2, whose experimentally registered temperature fields are presented in Figures 2 and 3, respectively.

The initial temperature, t w0-avg; basic density, ρ b; and moisture content, u, of the logs during the experiments were as follows:

• For log P1: t w0-avg = −29.7°C, ρ b = 359 kg·kg−1, and u = 1.44 kg·kg−1.

• For log P2: t w0-avg = −28.0°C, ρ b = 364 kg·kg−1, and u = 1.78 kg·kg−1.

As it was mentioned above, the duration of the freezing and duration of the subsequent thawing of the logs were equal to 50 and 70 h, respectively.

The model was solved with step Δr = Δz = 0.006 m along the coordinates r and z, with step Δτ = 6 s [8, 23], and with the same initial and boundary conditions, as they were during the experimental research.

During the solving of the model, mathematical descriptions of the thermophysical characteristics of poplar sapwood with ufsp293.15=0.35kg·kg−1, K wr = 1.48, and K wp = 2.88 [7, 9] have been used.

### 4.3 Inverse determination of the heat transfer coefficients during solving of the model

The model (1) to (4) was solved with various values of the exponent x in Eqs. (27) and (28). The computed by the model change of t in the four representative points of the logs with each of the tested values of the exponent x during the thawing was compared mathematically with the corresponding one experimentally registered change of t in these points with an interval of 15 min.

The aim of this comparison was to determine the values of x, which ensure the best compliance between the computed and experimentally registered temperature fields in subjected to thawing logs.

As a criterion of the best compliance, the minimum average value of RSME, σ avg, was used, which is equal to

σavg=n=1Np=1Ptp,ncomptp,nexp2PN1,E45

where tp,ncompand tp,nexpare the computed and experimentally registered temperatures in the representative points; p is the number of the representative points of the logs, р = 1, 2, 3, 4, i.e., Р = 4 was inputted into Eq. (45); n is the number of the moments of the thawing process, (n = 1, 2, 3, …, N = τthaw)/(150Δτ) = 252,000 s/900 s = 280, because of the circumstance that the comparison of the computed values of t with experimentally registered values in the same points was made with an interval of 15 min = 900 s = 150Δτ.

For the calculation of σ avg, a software program in the calculation environment of MS Excel was prepared. At τthaw = 70 h = 252,000 s, RSME has been calculated with the help of the program simultaneously for a total of N·P = 1120 temperature–time points during the thawing of each log.

It was determined that the minimum values of RSME overall for the studied four representative points are equal to σ avg = 1.37°C for log P1 and to σ avg = 1.34°C for log P2. These minimum values of σ avg correspond to the following values of the exponent x in Eqs. (27) and (28), which were obtained during the solving of the inverse task, x = 0.22 for log P1 and x = 0.20 for log P2.

Figures 7 and 8 present the calculated change in α wr and α wp during the studied thawing process of the logs P1 and P2, respectively. Figure 7.Calculated change in αwr and αwp of the log P1 during its 70 h thawing. Figure 8.Calculated change in αwr and αwp of the log P2 during its 70 h thawing.

Figures 9 and 10 present the calculated change in t m and also in the logs’ surface temperature t s and t of 4 representative points of the studied logs. Figure 9.Experimentally determined and calculated change in tm, ts, and t in four points of the log P1 during its 70 h thawing. Figure 10.Experimentally determined and calculated change in tm, ts, and t in four points of the log P2 during its 70 h thawing.

It can be seen that with the decrease of the difference between t m and t s during the logs’ thawing, the heat transfer coefficients on Figures 7 and 8 gradually decrease, as follows:

• At α wr: from 2.3 to 1.0 W·m−2·K−1 for P1 and from 1.9 to 1.2 W·m−2·K−1 for P2.

• At α wp: from 5.1 to 2.2 W·m−2·K−1 for P1 and from 4.5 to 2.8 W·m−2·K−1 for P2.

Using the obtained change in the heat transfer coefficients, the change in the logs’ surface temperature during the thawing, t s, has been calculated by the model (refer to Figures 9 and 10).

The comparison to each other of the analogical curves in Figures 2 and 9, and also in Figures 3 and 10, shows good conformity between the calculated and experimentally determined changes in the very complicated temperature fields of the studied logs during their thawing.

During our extensive simulations with the model (1) to (4), we established good qualitative and quantitative compliance between computed and experimentally determined temperature fields of logs from numerous wood species with different moisture content above the hygroscopic range .

The overall RSME for the studied four representative points in the logs does not exceed 5% of the temperature ranges between the minimal and maximal temperatures of each log during its thawing.

## 5. Conclusions

This chapter describes the creation, solving, and validation of a 2D nonlinear mathematical model for the transient heat conduction subjected to thawing frozen logs in an air environment.

The mechanism of the heat distribution in logs during their thawing has been described by a 2D equation of heat conduction at convective boundary conditions. For the numerical solving of the model with the help of explicit form of the finite-difference method, a software package has been prepared in the calculation medium of Visual FORTRAN Professional developed by Microsoft.

A validation of the model towards our own experimentally determined 2D temperature distribution in poplar logs with a diameter of 0.24 m, length of 0.48 m, and initial temperature about –30°C during their 70 h separate thawing at room temperature has been carried out.

During the validation of the model, the inverse problem has been solved for the determination of the logs’ heat transfer coefficients in radial and longitudinal directions. This problem has been solved also in regard to the logs’ surface temperature, which depends on the mentioned coefficients.

The following minimum values of the average RSME total for the temperature change in four representative points in each of the studied logs have been obtained:

• σ avg = 1.37°C for log P1 with ρ b = 359 kg·m−3 and u = 1.44 kg·kg−1.

• σ avg = 1.34°C for log P2 with ρ b = 364 kg·m−3 and u = 1.78 kg·kg−1.

During the solving of the inverse task, it was determined that the heat transfer coefficients subjected to thawing logs decrease gradually, as follows:

• At α wr: from 2.3 to 1.0 W·m−2·K−1 for P1 and from 1.9 to 1.2 W·m−2·K−1 for P2.

• At α wp: from 5.1 to 2.2 W·m−2·K−1 for P1 and from 4.5 to 2.8 W·m−2·K−1 for P2.

Good adequacy and precision of the model towards the results from extensive own experimental studies allow for the carrying out of various calculations with it, which are connected to the nonstationary temperature distribution in logs during their thawing in an air environment. For example, as a result of such calculations, it is possible to determine the real initial temperature of logs depending on their dimensions, wood species, moisture content, and the temperature of the air near the logs during their many days staying in an open warehouse before the thermal treatment in the production of veneer.

The information about the real value of that immeasurable parameter is needed for scientifically based computing of the optimal, energy saving regimes for thermal treatment of each specific batch of logs.

The model of the logs’ thawing process can be applied also in the software for controllers used for advanced model predictive automatic control [20, 21, 32] of this treatment. The approach for solving of the inverse task of the heat transfer in this chapter could be further applied in the development and solving of analogous models, for example, for the calculation of the temperature fields during freezing or thawing processes of different wooden and other capillary porous materials.

## Acknowledgments

This document was supported by the APVV Grant Agency as part of the project, APVV-17-0456, as a result of work of authors and the considerable assistance of the APVV agency.

## Nomenclature

a

temperature conductivity, m2·s−1

c

specific heat capacity, J·kg−1·K−1

D

diameter, m

g

acceleration of gravity, g = 9.81 m·s−2

Gr

Grashoff’s number of similarity

L

length, m

Nu

Nusselt’s number of similarity

Pr

Prandtl’s number of similarity

R

r

radial coordinate: 0 ≤ r ≤ R, m

T

temperature, K

t

temperature, oC

u

moisture content, kg·kg−1 = %/100

w

kinematic viscosity coefficient, m2·s−1

x

exponent, −

z

longitudinal coordinate: 0 ≤ z ≤ L/2, m

α

heat transfer coefficients between log’s surfaces and the surrounding air medium, W·m−2·K−1

β

coefficient of the volume expansion of the air, K−1

λ

thermal conductivity (for wood or air), W·m−1·K−1

ρ

density, kg·m−3

σ

root-square-mean error (RSME), °C

τ

time, s

φ

relative humidity, %

Δr

step along the coordinates r and z for solving of the model, m

Δτ

step along the time coordinate for solving of the model, s

### Subscripts

a

air

avg

average (for mass temperature of logs or for root-square-mean error)

b

basic (for wood density, based on dry mass divided to green volume)

bw

bound water

bwm

maximum possible amount of the bound water in the wood

comp

computed

exp.

experimental

fr

freezing

fre

end of freezing

fsp

fiber saturation point

fw

free water

i

current number of the knot of the calculation mesh in the direction along the log’s radius: i = 1, 2, 3,…, 21 = (Rr + 1)

k

current number of the knot of the calculation mesh in longitudinal direction of the logs: k = 1, 2, 3, …, 41 = (L/2/Δr + 1)

m

medium (for temperature of the air environment near the logs during their thawing process)

p

parallel to the wood fibers

r

s

surface

thaw

thawing

w

wood

we

wood effective (for specific heat capacity)

w-fr

wood with frozen water in it

w-nfr

wood with fully liquid water in it

w0p

parallel to the wood fibers at °C

w0r

radial direction of wood at °C

0

initial or at 0°C

1,2,3

1st, 2nd, 3rd (for temperature ranges of the logs’ thawing process)

@

at

&

and simultaneously with this

### Superscripts

n

current number of the step Δτ along the time coordinate during solving of the model: n = 1, 2, 3, …, N = τthaw/Δτ

272.15

at 272.15 K, i.e., at –1°C

293.15

at 293.15 K, i.e., at 20°C

chapter PDF
Citations in RIS format
Citations in bibtex format

## More

© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## How to cite and reference

### Cite this chapter Copy to clipboard

Nencho Deliiski, Ladislav Dzurenda and Natalia Tumbarkova (July 9th 2020). Modeling of the Two-Dimensional Thawing of Logs in an Air Environment, Modeling and Simulation in Engineering - Selected Problems, Jan Valdman and Leszek Marcinkowski, IntechOpen, DOI: 10.5772/intechopen.93177. Available from:

### Related Content

#### Modeling and Simulation in Engineering

Edited by Jan Valdman

Next chapter

#### Innovations in Heat Pump Design Using Computational Fluid Dynamics with Control Volume Method

By Cemil Koyunoğlu

#### Applications from Engineering with MATLAB Concepts

Edited by Jan Valdman

First chapter

#### Digital Image Processing with MATLAB

By Mahmut Sinecen

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.