An Inverse Boundary Design Problem in a Radiant Enclosure

Transcription

1 he 6 th. International Chemical Engineering Congress & Exhibition IChEC November 009, ish Island, Iran An Inverse Boundary Design Problem in a Radiant Enclosure A. arimipour *, S.M.H. Sarvari, Ali arimipour 3, E. Abedini 4 Ph.D. candidate, Islamic Azad Uni., Najaf Abad Branch, Arasharimipour@gmail.com Associate professor, University of Sistan & Baluchestan, sarvari@hamoon.usb.ac.ir 3 B.S., Sima Esfahan Company, Aliarimipour@ymail.com 4 Ph.D student, University of Sistan & Baluchestan, Ehsan.abedini@gmail.com Abstract his wor presents an inverse radiative design problem in which the aim is to determine the distribution of heater strengths over the heater surface to produces a desired temperature and heat flux distributions over the design surface of a radiant oven. he walls are assumed to be diffuse-gray and the medium is considered to be transparent. he conjugate gradient method is used to minimize an objective function which is defined by the sum of square errors between estimated and desired heat fluxes over the design surface. he inverse estimation consists of three problems; the direct problem, the inverse problem and the sensitivity problem. he direct problem and the sensitivity problem are solved by the net radiation method, and the configuration factors are calculated by the Hottel s crossed-string method. Our investigation shows that the presented algorithm is able to estimate heater strengths accurately. eywords: Radiation, Inverse boundary design, Conjugate Gradient Method (CGM). Introduction Radiation is the important mode of heat transfer in high temperature devices such as combustion chambers and furnaces. Inverse radiation design problems are concerned with the determination of unnown heater settings over the heater surface or unnown heat source distribution in the medium to produce both desired heat flux and temperature distributions over some part of the radiant enclosure, namely the design surface. he desired conditions over the design surface are dependent on the process for which the thermal system is built. he design conditions can be reached by controlling the thermal conditions in other parts of the system, namely the heater surface. herefore, the ultimate goal in thermal design problems is to find the required powers of heaters over the heater surface to satisfy both specified boundary conditions over the design surface. he forward design procedure is based on the nowledge of one and only one thermal condition on each element of the system. herefore, the forward design process involves a trial-and-error procedure, which requires a great deal of effort and time. Using the inverse design methodology is in fact an approach to avoid the trial and error procedure to meet the design goal. A comprehensive review of radiative heat transfer in combustion systems has been given by Visanta and Menguc []. Inverse analysis of radiative transfer is concerned with the determination of the radiative properties, boundary conditions and the temperature field or source term distribution from different inds of radiation measurements. A thorough review of inverse radiation problems has been given by McCormic []. A lot of wor has been reported on the estimation of radiative properties. Many researchers have also dealt with the inverse problems for determining the temperature profile or source term in media. Li and Ozisi [3], and Liu [4] have reconstructed the temperature profiles or source terms in plane-parallel, spherical, and cylindrical media by the inverse analysis from the data of the radiation intensities exiting the boundaries. Sarvari et al. [5,6,7] reported a methodology for designing radiant enclosures containing absorbing-emitting media to find the appropriate heater settings. Sarvari and Mansouri [8] used the optimization procedure to determine the heat source distribution in participating media. his procedure was extended by Sarvari [9] to determine the heat source distribution in conductiveradiative media. In this paper an inverse radiation boundary design problem for an enclosure filled with transparent medium is investigated numerically. he configuration factors are calculated by the Hottel s crossed-string method and the CGM is used to minimize an objective function which is defined by the sum of square errors between estimated and desired heat fluxes over the design surface. he direct problem and the sensitivity problem are solved by the net radiation method. Problem Description he geometry and the coordinate system of the radiant oven with diffuse gray walls and two insulated side walls, is illustrated in Fig.. he product surface (design surface) is placed on the bottom wall, whereas the heaters are placed * corresponding author

2 on the heater surface at the top of the oven. he medium is transparent and all the boundary conditions over the oven s walls are specified except for the heater surface. In order to meet the design goal, both the temperature and the heat flux require having uniform distributions over the product. he goal of the design problem is to identify the heater flux, q h, that produce the desired uniform heat flux, q d, over the temperature specified design surface with uniform temperature 4 of d (or with uniform emissive power, Ed = εσd ). Inverse radiative design method is applied to determine the distribution of temperature and heat flux sources on the heater surface, according to the distribution of desired temperature and heat flux on the design surface. Fig.. he Schematic shape of the radiant oven with gray-diffuse walls and transparent medium In first part of this paper the values of Ed =.0 ( w/ m ) and q d has the profile is shown in Fig., so they are nown parameters and the unnowns are E h =?, q h =?, and in the second part, the values of unnowns E h =?, q h =? are calculated while Ed =.0 ( w/ m ) and qd Numerical Procedure he inverse estimation consists of three problems: =.7 ( w/ m ) are nown and constant. -Direct problem: Suppose an enclosure with discrete internal surfaces involving two types of boundary conditions, the surfaces with specified temperature and the surfaces with specified heat flux, where the objectives are to analyze the radiation exchange between the surfaces. he net radiation method is used to solve the radiation exchange in a radiant enclosure. In this method, the boundary is subdivided into surface elements. he equation of radiation exchange for surface elements with specified temperature (emissive power) and for other surface elements with specified heat flux can be described by the following equations: δj ( ε) F j J j = ε E, and ( δj j) j F J = q, + (a,b) where δ j is the ronicer delta. he set of Eq. () is solved to calculate the outgoing heat fluxes, J j,j =,,, then the unnown boundary condition (emissive power or heat flux) is determined by the following equation: ε ε q + J = E, () he configuration factors are calculated by the Hottel s crossed-string method. -Inverse problem: For the inverse problem considered here, the heat flux distribution over the heater surface, qh ( z ), is regarded as unnown, and the desired heat flux distribution over the design surface, qd ( x ), is available for the analysis. he heat flux distributions over the design and the heater surfaces may be expressed as vectors of discrete elemental values, such as qh = { q h,n n =,,N }, qd = { qd,m m=,,m}

3 where N and M are the number of heaters and the number of surface elements on the design surface, respectively. he estimated heat fluxes by the direct method over the design surface is expressed as qe (x d ) { qe,m m,,m} solution of the inverse problem is based on the minimization of an objective function given by: = =. he ( q ) = q q ( q ) q q ( q ) G h d e h d e h (3) he minimization procedure is performed using the conjugate gradient method. Iterations are built in the following manner: h h β + = + and ( h ) q q C G C = q + α C and β where C is the direction of descent vector and Here, S, G ( q h ) and which are defined as follows: ( ) S = S = q q h q α mn e,m h,n G ) G ) ) G ) β ( ) d e h SC q q q = SC SC is the step size. (4a,b,c) α are the sensitivity matrix, gradient direction vector and conjugate coefficient, respectively, m n 0 = with α = 0 G and ) = S qd qe ) G 3-Sensitivity problem: o minimize the objective function given by Eq. (3), we need to calculate the components of the sensitivity matrix, S, defined by Eq. (5a). he sensitivity problem is obtained by differentiating the direct mn problem given by Eqs. () with respect to the elemental heat fluxes over the heater surface, q h,n. Hence, for temperature specified surface elements we have Eq. (6a) and for surface elements with specified heat flux we have Eq. (6b) (5a,b) (5c) δj ( ε) F j J j qh,n = 0, and ( δj F j ) J j q h,n = δn,+ (6a,b) he magnitudes of J j qh,n are calculated by solving Eqs. (6), then the elements of sensitivity matrix can be obtained by differentiation of Eq. () with respect to the elemental heat fluxes over the heater surface, q h,n. ε S = q q = J q, m M m mn e,m h,n m h,n ε m (7) he n-th column of the sensitivity matrix is obtained by solving Eq. (7). he set of equations (6) and (7) must be solved N times in order to complete the sensitivity matrix. Because the right hand sides of equations (6a) and (6b) are constants, then the magnitudes of Jm qh,n is not a function of time; hence, the sensitivity matrix, S, can be obtained before starting the iterative procedure of the inverse problem. Results and Discussion In this paper an inverse radiation boundary design problem is investigated numerically for an enclosure with unnowns E h =?, q h =?, on the heater surface while Ed =.0 ( w/ m ) and q d (desired) are nown (Fig. ). Our investigation is done by writing a computer code with a useful computational algorithm to apply the conjugate gradient method. In Fig. the comparison between the estimated heat flux, q e, and desired heat flux, shown. q ( w/ m ), over the design surface is d

4 his figure shows a good agreement between the values of qe and q d, so that we can find the distribution of the emissive power and the heat flux on the heater surface which lead to this profile of q e. In Fig. 3 and Fig. 4, the distribution of the heat flux, q h, and the emissive power on the heater surface are shown. In the next part of this paper the profile value of unnowns E h =?, q h =? are calculated while Ed =.0 ( w/ m ) and qd =.7 ( w/ m ) are nown. In Fig. 5 the comparison between the estimated heat flux, q e, and desired heat flux, q d, over the design surface is shown. his figure shows when the E d and q d have the constant value which any relationship to the X, the inverse radiation method only can estimate the q e profile which is close to value of q d and so its accuracy is acceptable. Conclusion In this paper an inverse radiation boundary design problem in enclosure with diffuse gray walls was investigated numerically. he conjugate gradient method was used to minimize an objective function which is defined by the sum of square errors between estimated and desired heat fluxes over the design surface. It was shown that the method used in this paper for inverse problem would have been a good accuracy in order to apply in experimental and industrial problems q d q e q -.6 q h Fig..Comparison between the estimated heat flux X qe ( wm / ) and desired heat flux q ( wm / ) on the design surface d Fig. 3. Distribution of Z qh ( wm / ) on the heater surface q e q d -.69 Eh q Z X Fig. 5.Comparison between the estimated heat flux qe ( wm / ) and

5 Fig. 4. Distribution of Eh ( wm / ) on the heater surface desired heat flux qd ( wm / ) on the design surface References []- Visanta, R., Menguc, M.P., Radiation heat transfer in combustion Systems, Prog. Energy Combust. Sci. 3, pp , 987. []- McCormic, N. J., Inverse radiative transfer problems: a review, Nucl.Sci. Eng., pp , 99. [3]- Li, H. Y., Ozisi, M. N., Identification of the temperature profile in an absorbing, emitting, and isotropically scattering medium by inverse analysis, J. Heat ransfer 4, pp , 99. [4]- Liu, L. H., an, H. P., Yu, Q.Z., Simultaneous identification of temperature profile and wall emissivities in onedimensional semitransparent medium by inverse radiation analysis, Numer. Heat ransfer, Part A 36, pp. 5 55, 999. [5]-Sarvari, S. M. H., Mansouri, S. H., Howell, J. R., Inverse Boundary Design Radiation Problem in Absorbing- Emitting Media with Irregular Geometry, Num. Heat rans., Part A, Vol. 43, pp , 003. [6]-Sarvari, S. M. H., Howell, J. R., Mansouri, S. H., Inverse Boundary Design Conduction-Radiation Problem in Irregular wo-dimensional Domains, Num. Heat rans., Part B, Vol. 44, pp. 09-4, 003. [7]- Mehraban, S., Sarvari, S. M. H., Farahat, S., A QUASI-SEADY MEHOD FOR INVERSE DESIGN AND CONROL OF A WO-DIMENSIONAL RADIAN OVEN IN RANSIEN SAE, Proceedings of CH-08 ICHM International Symposium on Advances in Computational Heat ransfer, Marraech, Morocco, May -6, 008. [8]-Sarvari, S. M. H., Mansouri, S. H., Inverse Design for Radiative Heat Source in an Irregular -D Participating Media, Num.l Heat rans., Part B, Vol. 46, pp , 004. [9]-Sarvari, S. M. H., Inverse Determination of Heat Source Distribution in Conductive-Radiative Media with Irregular Geometry, J. Quant. Spect. and Rad. rans., Vol. 93, pp , 005.