Dear Label_Name,

Greetings! I hope you are well.

I am writing on behalf of an international publishing house, Science and Education Publishing, USA (http://www.sciepub.com).

Science and Education Publishing (SciEP) is a multidisciplinary open access publisher of Journals covering the fields of Science, Technology and Medicine. All SciEP publications are subject to high-quality peer review, editorial and production processes and are freely available online.

Welcome to the April, 2013 edition of SciEP newsletter, providing monthly summaries of articles on science & education related topics of interest to scientist and professionals. The mission of our journals is to provide an open access forum for publishing new findings on issues pertinent to scientific and academic field. We also invite you to submit original research papers, review articles, short reports, or case studies are considered for publication in our Multidisciplinary Journals.

American Journal of Public Health Research
DOI: 10.12691/ajphr-1-1-3

The Implementation of a Tele-Homecare System with Service Oriented Architecture and HL7 Message Transmission Standard
Ching-Sung Wang

American Journal of Medical and Biological Research
DOI: 10.12691/ajmbr-1-1-2

The Influence of Ethanol on Pyruvate Kinases Activity in Vivo, in Vitro, in Silico
Lelevich Sergey Vladimirovich, Khrustalev Vladislav Victorovich, Barkovsky Eugene Victorovich, Shedogubova Tatyana Aleksandrovna

American Journal of Medical Sciences and Medicine
DOI: 10.12691/ajmsm-1-1-1

Clinical and Histological Features in Wegener Granulomatosis
Li Cavoli Gioacchino, Passantino Rita, Bono Luisa, Tortorici Calogera, Ferrantelli Angelo, Giammarresi Carlo, Rotolo Ugo

American Journal of Infectious Diseases and Microbiology
DOI: 10.12691/ajidm-1-1-2

The Curse of Lymphatic Filariasis: Would the Continual Use of Diethylcarbamazine Eliminate this Scourge in Papua New Guinea?
Anne Lanham, Lillian Mwanri

We have released issues of some journals on February 2013.
American Journal of Public Health Research
American Journal of Medical Sciences and Medicine
American Journal of Medical and Biological Research
American Journal of Pharmacological Sciences
American Journal of Clinical Medicine Research
American Journal of Nursing Research
American Journal of Cancer Prevention
American Journal of Medicine Studies
American Journal of Infectious Diseases and Microbiology
American Journal of Cardiovascular Disease Research
American Journal of Biomedical Research
American Journal of Mechanical Engineering
American Journal of Microbiological Research
Journal of Preventive Medicine
American Journal of Food and Nutrition
American Journal of Hypertension Research
American Journal of Sports Science and Medicine
Journal of Cancer Research and Treatment
American Journal of Zoological Research
Journal of Biomedical Engineering and Technology
American Journal of Microbiological Research
Biomedical Science and Engineering

We are pleased to announce that submissions are now being accepted for peer review for all the journals across all our subject areas in SciEP.

Please submit your manuscript via our online submission system:
http://www.sciepub.com/author/signin

To unsubscribe from our mailing list, please click on the following link: unsubscribe

I am looking forward to hearing from you.

Kind regards,

David Rain
Email:join@sciepub.com
Science and Education Publishing (http://www.sciepub.com)

Analytical and FDM Solutions for Anisotropic Heat Conduction in an Orthotropic Rectangular

Mohammad Sameti, Fatemeh Razi Astaraie, Fathollah Pourfayaz*, Alibakhsh Kasaeian

Department of Renewable Energies, Faculty of New Sciences & Technologies, University of Tehran, Tehran, Iran

Abstract: Anisotropic materials have physical properties such as thermal conductivity which vary with orientation throughout the structure. These properties make the analysis of heat transfer much more complex than standard isotropic materials. Therefore, the heat flux equation and Fourier’s Law should be expanded in order to properly govern the thermal conduction. The analytical exact solutions and numerical finite difference solutions of a fundamental heat conduction problem in an anisotropic medium are presented in this study. The problem is solved for steady state condition in a rectangular thin film of quartz. The key for analytical solution is to use two variables to convert the anisotropic equations into isotropic ones. The temperatures of 35 points on the rectangular are obtained and compared by the two methods. The results show an excellent accuracy for the numerical method. The temperature and the heat flux profile are also illustrated to show how the anisotropicity influences the thermal behavior of the media. MAPLE and EES were used to numerically integrate or solve the system of the eqautions.

Keywords: anisotropic materials, heat conduction, Finite Difference Method (FDM)

Received February 21, 2014; Revised March 12, 2014; Accepted March 14, 2014

Citation: Mohammad Sameti, Fatemeh Razi Astaraie, Fathollah Pourfayaz, Alibakhsh Kasaeian, Analytical and FDM Solutions for Anisotropic Heat Conduction in an Orthotropic Rectangular. American Journal of Numerical Analysis, 2 (2). 65-68. January 2014. doi:10.12691/ajna-2-2-6

1. Introduction

There are many natural and synthetic materials in which the thermal conductivity varies with direction and are called anisotropic materials. For example, crystals, wood, sedimentary rocks, metals that have undergone heavy cold pressing, laminated sheets, cables, heat shielding materials for space vehicles, fiber reinforced composite structures and many others are anisotropic materials. Ma and Chang [1] considered a two-dimensional steady state thermal conduction problem of anisotropic multi-layered media. They used the Fourier transform and the series expansion technique to find the explicit closed form solutions of the problem. Maruyama et al. investigated the anisotropic heat transfer in a single-wall carbon nanotube using molecular dynamic simulation [2]. Mera et al. solved the steady state anisotropic heat conduction problems by the implementation of the boundary element method (BEM) [3]. Sladek et al. used the meshless methods based on the Petrov-Galerkin approach to solve the steady and transient heat conduction problem in a continuously nonhomogeneous anisotropic medium [4]. Materials used for renewable heating and cooling of buildings can also have the anisotropic behavior [5].

In this study the steady state heat conduction is of interest and the finite difference method is implemented as the numerical solution. Two methods were combined to solve the problem, analytically. Then, the energy balance formula is used as in to numerically model the equation with finite difference method. In rectangular coordinates for two dimensional conduction, the magnitudes of heat flux qx and qy in the x and y directions can be determined from Fourier’s low:

[picture] | (1)

[picture] | (2)

Where T is the temperature and kij is the thermal conductivity. From thermodynamics, the determinant and the elements of the matrix M defined as:

[picture] | (3)

Should be positive. This known as Onsager’s reciprocity relation [6].

2. Formulation of the Problem

A rectangular thin film which is made of quartz is considered in this study. As Figure 1 shows, owing to its low height l, instead of solving the three dimensional problem, the problem is assumed to be two dimensional. Temperature variations of the medium is negligible in the z direction, however it varies in the directions of x and y with different rates. The thermal and geometric specifications of the rectangular can be found in Table 1. The problem can be expressed as follows:

Figure 1. Orthotropic rectangular considered for the anisotropic heat conduction

[picture] | (4)

Where g is the heat generation, ρ is the density and C is the specific heat. For steady state condition with no heat source, the equation is simplified to:

[picture] | (5)

Where 0<x<6 and 0<y<4. The boundary conditions for the temperature distribution on the sides are given by:

[picture] | (6)

[picture] | (7)

2.1. Analytical Solution

The heat conduction equation for an orthotropic medium can be transformed to a standard heat conduction equation for an isotropic medium [7]. New independent variables X and Y are defined as:

[picture] | (8)

[picture] | (9)

Where k is a reference conductivity. Equation (1) then becomes:

[picture] | (10)

The choice of the reference conductivity is not arbitrary. The reason for this is that a volume element in the original space dxdy transforms using equations (8) and (9) into equation (11):

[picture] | (11)

If the quantities ρCp and the generation term g defined on the basis of unit volume should have the same physical significance, we should have:

[picture] | (12)

Now, the equation (10) which is the 2D potential equation can be solved using the separation of variables. The method is well-known and applicable for such problems [8]. Assuming:

[picture] | (13)

Where T₁(x) and T₂(y) are independent functions. By introducing equation (13) into equation (10), it is obtained that:

[picture] | (14)

[picture] | (15)

Where is the second derivative of the function . The boundary conditions then become:

[picture] | (16)

[picture] | (17)

The length a and the width a of the rectangular is transformed to A and B, respectively.

[picture] | (18)

[picture] | (19)

Solving equation (14), it is obtained that the problem is the generator for the sinusoidal Fourier series where the eigenfunctions and eigenvalues are:

[picture] | (20)

[picture] | (21)

Where n is a natural number. Substituting equation (20) in equation (15) and solving it with the condition (17), yields:

[picture] | (22)

Therefore, for each particular number n, the solution for the temperature is:

[picture] | (23)

The final answer is the summation of al , so

[picture] | (24)

The boundary condition (6) should now be used, but it is necessary to transform it to a new function using equations (8) and (9):

[picture] | (25)

Applying the boundary condition we have:

[picture] | (26)

Which is solved to find the Fourier coefficients as:

[picture] | (27)

Table 1. Specifications of the orthotropic media

The temperature at each point can now be found using MAPLE to numerically integrate equation (27) and to find the series answer in equation (26). Furthermore, the temperature profile and the heat flux profiles can also be illustrated.

2.2. Numerical Solution

A mesh can be made with 35 nodes in which the horizontal step and the vertical step are both equal to 1 cm. The grid is illustrated in Figure 2. In a steady state condition, the energy balance can be written for each of the nodes [9]. For example, for the node number 9 the suitable equation is:

[picture] | (28)

The general equation for all of the nodes can be written as follows:

[picture] | (29)

Where m falls within the following interval.

[picture] | (30)

For the other nodes, the temperature is given in equations (2) and (3) as:

[picture] | (31)

Now, 15 equations are created which can be solved using the software EES (Engineering Equation Solver).

Figure 2. Surface temperature profile (blue line) and position of the nodes

3. Results and Discussions

The temperature for the 35 nodes is calculated by the both methods and can be found in Table 2. As it can be seen, the results are very similar and the numerical method is accurate enough. The isothermal curves and the heat flux vectors are also illustrated in Figure 3 and Figure 4. It can be understood from the results that:

1. The heat flux vectors are not perpendicular to isothermal curves (despite the isotropic heat conduction), because the heat flows in an anisotropic way in the medium.

2. The slope of all vectors is more than 1 or the angles are more than 45°. This happens because the conductivity in the horizontal direction is less than the vertical direction.

3. The variations in the surface temperature in equation (2) have led to the accumulation of isothermal curves in the left side and top of the rectangular (edge effects).

4. The temperature difference between the two methods reach its maximum level at node 9, because the surface temperature described by the equation (2) suddenly increases from 0 to 50 in that region.

Table 2. Results and comparison of nodes temperatures

Figure 3. Isothermal profiles in the orthotropic rectangular

Figure 4. Heat flow directions in the orthotropic rectangular

4. Conclusions

Mathematically, it is understood that the most important difference between the isotropic and anisotropic heat conduction is that the heat flux vectors are not perpendicular to isothermal curves. Steady state anisotropic heat conduction problem is solved numerically and analytically for an orthotropic rectangular thin film. The analytical method is based on the introducing two new variables to transform the equations to the well-known 2D potential equation which can be solved by the separation of variables. The numerical method is based on the energy balance and the finite difference method. The results are compared and the results show an excellent accuracy.

References

[1]   | Ma, C. C., & Chang, S. W. (2004). Analytical exact solutions of heat conduction problems for anisotropic multi-layered media. International journal of heat and mass transfer, 47 (8), 1643-1655.
[2]   | Maruyama, S., Igarashi, Y., Taniguchi, Y., & Shiomi, J. (2006). Anisotropic heat transfer of single-walled carbon nanotubes. Journal of Thermal Science and Technology, 1, 138-148.
[3]   | Mera, N. S., Elliott, L., Ingham, D. B., & Lesnic, D. (2001). A comparison of boundary element method formulations for steady state anisotropic heat conduction problems. Engineering analysis with boundary elements, 25 (2), 115-128.
[4]   | Sladek, J., Sladek, V., & Atluri, S. N. (2004). Meshless local Petrov-Galerkin method for heat conduction problem in an anisotropic medium. Computer Modeling in Engineering and Sciences, 6, 309-318.
[5]   | Kasaeian, A., Sameti, M., & Eshghi, A. T. (2014). Simplified Method for Night Sky Radiation Analysis in a Cool-Pool System. Science and Education, 2(1), 29-34.
[6]   | Ozisik, N. (1993). Heat conduction. John Wiley & Sons.
[7]   | Su, S., Chen, J., & Zhang, C. (2011). Study on Performance of Anisotropic Materials of Thermal Conductivity. Open Civil Engineering Journal, 5.
[8]   | Kreyszig, E. (2010). Advanced engineering mathematics. John Wiley & Sons.
[9]   | Incropera, F. P., Lavine, A. S., Bergman, T. L., & DeWitt, D. P. (2013). Principles of heat and mass transfer. Wiley.