NOMENCLATURE. H j a r s u o. u X. U x o
|
|
- Daniel Horton
- 5 years ago
- Views:
Transcription
1 Numerical Heat Transfer, Part A, 50: , 2006 Copyright # Taylor & Francis Group, LLC ISSN: print= online DOI: / INVERSE RADIATION DESIGN PROBLEM IN A TWO-DIMENSIONAL RADIATIVELY ACTIVE CYLINDRICAL MEDIUM USING AUTOMATIC DIFFERENTIATION AND BROYDEN COMBINED UPDATE Ki Wan Kim and Seung Wook Baek Division of Aerospace Engineering, Department of Mechanical Engineering, Korea Advanced Institute of Science and Technology, Yuseong-Gu, Taejon, Republic of Korea Inverse radiation design analysis in a two-dimensional, concentric, cylindrical, absorbing, emitting and scattering medium has been conducted, given desired boundary conditions on the design surface. The finite-volume method was adopted to solve the radiative transfer equation and the energy conservation equation in the direct problem, while the Levenberg- Marquardt method was used to solve a set of equations, which are expressed by errors between estimated and desired radiative heat fluxes on the design surface. In order to diminish the computational time required for calculating sensitivity matrix, automatic differentiation as well as the Broyden combined update were utilized. INTRODUCTION Inverse radiation analysis has been utilized to estimate radiative properties when the other measured radiative data were available or to deduce the optimal condition satisfying a design goal [1]. In the former case, medium properties such as extinction coefficient, absorption coefficient, single scattering albedo, phase function, optical depth, and temperature, as well as surface properties such as boundary emissivity and temperature, have been estimated from measured intensities exiting from the boundary or temperature inside the medium [2 6]. This concept has also been applied to optical tomography for examining cross-sectional images of the human body in medical science [7]. In the latter case, boundary or source design problems have been solved to find optimal distribution of boundary condition or source strength satisfying the desired condition on design surfaces in irregular as well as regular geometries containing either a nonparticipating or a gray participating medium [8 12]. However, cylindrical geometry problems have not attracted as much Received 1 October 2005; accepted 23 November The financial assistance by the Combustion Engineering Research Center at the Korea Advanced Institute of Science and Technology is gratefully acknowledged. Address correspondence to Seung Wook Baek, Division of Aerospace Engineering, Department of Mechanical Engineering, Korea Advanced Institute of Science and Technology, Guseong-Dong, Yuseong-Gu, Daejeon , Republic of Korea. swbaek@kaist.ac.kr 525
2 526 K. W. KIM AND S. W. BAEK NOMENCLATURE D m ci direction weights ~e x ;~e y ;~e z unit vector in x, y, z directions E b emissive power ð¼pi b Þ; W=m 2 I radiation intensity, W=m 2 sr I b blackbody radiation intensity, W=m 2 sr L cylinder height, m ~n i unit vector normal to control-volume surface i n e number of elements on design surface n h number of heaters n p number of unknown parameters P a vector of unknown parameters q R radiative heat flux, W=m 2 ~r position vector of intensity r in inner cylinder radius, m r out outer cylinder radius, m _ s direction vector of intensity S equations set T temperature of participating medium, K T i boundary temperature at the surface i X sensitivity matrix b o extinction coefficient ð¼j a þ r s Þ; m 1 DA i ; DV ith surface area and volume of the control volume DX control angle boundary emissivity at the surface i e i h H j a r s u o u X U x o polar angle measured from the z direction, rad nondimensional medium temperature absorption coefficient, m 1 scattering coefficient, m 1 space variable in the azimuthal direction measured from x axis, rad angular variable in the azimuthal direction measured from x 0 axis, rad scattering phase function single scattering albedo Subscripts d design surface e estimated value e; w; n; s; t; b index of control-volume surfaces E; W; N; S; T; B index of neighboring nodal points h heater surface P calculating nodal point R radiation ref reference value w wall Superscripts k iteration number m; m 0 radiation direction m þ ; m boundaries of the control volume * nondimensional quantity attention [12], while the scattering term in the radiative transfer equation was not taken into account either [10, 11]. Various inverse methods have been adopted to obtain a stable solution in spite of the ill-posed characteristic of inverse problems. Even though matrix solver-based methods do not need an iteration procedure, special treatment is required to deal with the ill-conditioned nature of the equations [8, 9]. On the other hand, optimization technique-based methods are emerging as a stable inverse method due to multidimensionality as well as complexity in inverse problems [2 7, 10, 11]. A successful estimation of unknown properties depends on the magnitude of the sensitivity coefficient. The bigger its magnitude, the easier it is to estimate the inverse solutions. Furthermore, a reduction in computational time for calculating those values is one of the main issues in gradient-based optimization technique. The sensitivity method deals with the calculation of sensitivity coefficients. So far, equation-based sensitivity methods, such as sensitivity problems [1, 4, 5], boundary-value problems [1, 10], or adjoint problems [1], have been prevalent. Recently, code-based methods have received much attention due to their advantages [7]. These methods do not need to derive an adjoint problem or a sensitivity problem, so a discretization procedure of relevant equations is not necessary. Computational code
3 INVERSE RADIATION DESIGN PROBLEM 527 for the sensitivity matrix can be extracted directly from the code of the direct problem. Automatic differentiation belongs to this category. The goal of this study is to reduce the computational time for calculating the sensitivity matrix. Our previous study shows that the effect of the precision of the sensitivity matrix on estimation accuracy is not crucial [13]. Because of this, the Broyden combined update, which is one of the secant methods, is used with automatic differentiation to further shorten the computational time required to calculate the sensitivity matrix. Therefore, the major contribution of this work is to propose automatic differentiation with Broyden combined update as an alternative sensitivity method and to extend the inverse radiation design problem to a two-dimensional cylindrical enclosure as a successful application. ANALYSIS AND MODELING Model Description Figure 1 shows a concentric cylindrical enclosure which is filled with an absorbing, emitting, and scattering gray medium, and its simplified geometry with boundary conditions. The walls are diffusely emitting and reflecting gray walls. The direct problem is to calculate the radiative heat flux distribution on the design surface, given boundary conditions, such as temperature, wall emissivity, and medium properties such as absorption and scattering coefficient. The temperature distribution inside the medium is determined by solving the energy conservation equation in radiative equilibrium, since radiation is the only dominant heat transfer in this work. Radiative heat flux over the wall including the heater surface is obtained by integrating the dot product of intensity and the unit normal vector over all solid angles, while the intensity distribution is acquired by solving the radiative transfer equation. In this work, the finite-volume method is adopted for solving the radiative transfer equation. In the inverse problem, the goal is to find the radiative heat flux distribution over the heater surface which satisfies the desired radiative heat flux as well as temperature on the design surface. The Levenberg-Marquardt method is used to solve a set of equations expressed by errors between the estimated and desired radiative heat flux on the design surface. In order to reduce the computing time for the sensitivity matrix, a combined method of automatic differentiation and Broyden combined update is adopted in this study. Detailed description of mathematical formulations and numerical methods used in this work is given below. Radiative Transfer Equation The nondimensionalized radiative transfer equation (RTE) governing radiation intensity for a gray medium at any position ~r along a path _ s through an absorbing, emitting, and scattering medium is given by di ð~r; s _ Þ ds þ I ð~r; _ s Þ¼ð1 x o ÞH 4 ð~rþþ x Z o Uð _0 s ; _ s ÞI ð~r; _ s ÞdX 0 4p 4p ð1þ
4 528 K. W. KIM AND S. W. BAEK Figure 1. Schematic of the physical system and its simplified geometry: (a) physical system; (b) simplified geometry (L ¼ 4, r out ¼ 1, r in ¼ 0.5). where H ¼ T T ref I ¼ I I b ðt ref Þ x o ¼ r s b o s ¼ b o s b o ¼ j a þ r s Here, Uð _0 s ; _ s Þ is the scattering phase function for radiation from incoming direction _0 _ s to scattered direction s. It is approximated by a finite series of Legendre polynomials such that Uð s _0 ; s _ Þ¼Uðcos WÞ ¼ XJ j¼0 C j P j ðcos WÞ ð2þ where the C j s are the expansion coefficients, and J is the order of the phase function.
5 INVERSE RADIATION DESIGN PROBLEM 529 The boundary condition for a diffusely emitting and reflecting wall can be written as I ð~r w ; _ s Þ¼e w ð~r w ÞH 4 ð~r w Þþ 1 e Z wð~r w Þ I ð~r w ; _0 s Þj _0 s ~n w j dx 0 ð3þ p _0 s ~nw <0 where e w is the wall emissivity and ~n w, which has a positive value when the ray travels from the wall to medium, is the unit normal vector to the wall. Equation (3) can also be expressed as a form of prescribed radiative heat flux: Z I ð~r w ; _0 s Þj _0 s ~n w j dx 0 ð4þ I ð~r w ; s _ Þ¼q R ð~r w Þþ 1 p _0 s ~nw <0 Here, nondimensional heat flux is q R ¼ q R=rTref 4. The radiative heat flux on the design surface is calculated as follows: q R ð~r w;d Þ¼ 1 Z I ð~r w ; _ s Þð~n w;d _ s Þ dx p X¼4p ð5þ Finite-Volume Method for Radiation In solving the RTE, the finite-volume method (FVM) is adopted for its convenience in selecting the solid angle while guaranteeing an exact global conservation of radiative energy. Chui s way to treat the axisymmetric geometry is employed [14, 15]. To obtain the discretization equation, Eq. (1) is integrated over a control volume, DV, and a control angle, DX m, in the axisymmetric orthogonal grid as shown in Figure 2. By assuming that the magnitude of the intensity in a control volume and control angle is constant, the following finite-volume formulation can be obtained: X I ;m i DA i D m ;m ci ¼ð IP þ S ;m r;p ÞDVDXm ð6aþ where i¼e;w;n;s;t;b Z u mþ D m ci ¼ X u m X Z h mþ h m ð s _ ~n i Þ sin h dh du X ð6bþ _ s ¼ sin h cos ux ~e x þ sin h sin u X ~e y þ cos h~e z ð6cþ ~n e ¼ sin u o;p ~e x cos u o;p ~e y ð6dþ ~n w ¼ sin u o;p þ~e x þ cos u o;pþ ~e y ð6eþ ~n n ¼ ~n s ¼ cos u o;p ~e x þ sin u o;p ~e y ð6f Þ
6 530 K. W. KIM AND S. W. BAEK Figure 2. Angular control angle and spatial control volume: (a) angular control angle; (b) spatial control volume. ~n t ¼ ~n b ¼ ~e z ð6gþ S ;m r ¼ð1 x o ÞH 4 þ x Z o I ;m0 U m 4p 0!m dx 0 X 0 ¼4p ð6hþ DX m ¼ Z u mþ X u m X Z h mþ h m sin h dh du X ð6iþ To relate the intensities on the control-volume surfaces to a nodal one, the step scheme, which is not only simple and convenient but also ensures positive intensity, is
7 INVERSE RADIATION DESIGN PROBLEM 531 adopted. Then, the final discretized equation for the FVM becomes a m P I X ;m P ¼ a m I I ;m I þ b m P ð7aþ I¼E;W;S;N;T;B a m I ¼ DA i D m ci;in ð7bþ a m P ¼ X DA i D m ci;out þ DVDXm i¼e;w;s;n;t;b ð7cþ b m P ¼ðS;m r;p ÞDVDXm ð7dþ where D m ci;out ¼ ZDX m ð~n i s _ ÞdX ~n i s _ > 0 ð7eþ D m ci;in ¼ ZDX m ð~n i s _ ÞdX ~n i s _ < 0 ð7f Þ The spatial and angular domains are discretized into control volumes and 8 20 control angles, and temperature inside the domain is calculated by solving the energy conservation equation in radiative equilibrium as follows: rq R ¼ 4b oð1 x o Þ H 4 1 Z I dx ¼ 0 ð8þ 4p 4p Levenberg-Marquardt Method (LMM) In applying the quasi-newton method to inverse analysis, a set of equations has to be formulated using the difference between the estimated and desired heat flux on the design surface as follows: S i ðpþ ¼q R;e;i ðpþ q R;d;i for i ¼ 1;...; n e ð9þ Here, P is an unknown parameter, which is radiative heat flux on the heater surface, q R;h, in this work. In the above equation, S can be expressed in Taylor series, and should be zero to find appropriate parameters such that SðPÞ ffis 0 ðp 0 ÞþrSðPÞðP P 0 Þ¼0 ð10þ where rsðpþ ¼ qsðpþ qp which is the matrix of sensitivity coefficients, X. ¼ qq R;e ðpþ qp
8 532 K. W. KIM AND S. W. BAEK After some manipulation of Eq. (10), we can derive the following relation: P kþ1 ¼ P k ðx k Þ 1 S k ð11þ Here, it must be noted that the above relation may be applied only for the case that the number of equations is the same as that of the unknown variables. However, based on the fact that the number of equations can be greater than that of the unknown variables in the inverse analysis, Eq. (11) is modified as follows: P kþ1 ¼ P k ðx k Þ T X k 1 ðx k Þ T S k ð12þ The Levenberg-Marquardt method has a damping term to cope with ill-posed characteristics such that P kþ1 ¼ P k ðx k Þ T X k þ m k X k 1 ðx k Þ T S k ð13þ where m is a positive scalar called the damping parameter, and X k ¼ diag½ðx k Þ T X k Š is a diagonal matrix [1]. Automatic Differentiation (AD) No matter how complicated the objective function, it is executed on a computer as a sequence of code events containing only elementary operations such as addition, subtraction, multiplication, or division. By applying the chain rule on every code event, sensitivity coefficients can be calculated. Based on this fact, a program code for the calculation of sensitivity coefficients can be extracted from that of the direct problem. For this reason, this method does not need to derive an adjoint problem or a sensitivity problem, so a discretization procedure for these problems is not necessary either. Rather, we need the following two steps for automatic differentiation. First, the objective function for differentiation should be broken into its most elementary forms. Each elementary form permits only one computational operation, such as addition, subtraction, multiplication, or division. Second, the code event is sequentially differentiated in the form of total derivatives. There are two approaches to calculate a sensitivity matrix in AD, i.e., source transformation and operator overloading. The former is to change the source code following the previously described manner such that the source code is decomposed into elementary code events. The latter does not decompose the source code, but instead uses a derived type of modern complier. It needs a modern complier to overload elemental operators such as intrinsic functions in FORTRAN. Since the variables for differentiation are a derived type comprised of the function value, the first, and the second-order derivatives, redefinition of operators is needed. In this approach, all we have to do is define all variables involved in the calculation of differentiation in the derived type. This requires, however, additional computational resources, since all variables involved in the calculation of differentiation should be defined in derived type. We adopted ADIFOR as the source transformation type in this work [16].
9 There are two modes, i.e., forward and backward, according to the direction for applying the chain rule to code events. In this study, forward mode is adopted, and the general formulation for differentiation by the chain rule in forward mode is written by qs i qe j ¼ XN k¼n p þ1 where N is the number of code events. INVERSE RADIATION DESIGN PROBLEM 533 qs i qe k i ¼ 1;...; n e and j ¼ 1;...; n p ð14þ qe k qe j Broyden Combined Update (BC) The various secant methods have been studied by many researchers for their convenience in calculating the gradient value [17, 18]. They utilize the function values of present and previous iterations to obtain the gradient information. Broyden s method for nonlinear systems belongs to the present type. The Levenberg- Marquardt method finds the solution with the following iteration procedure: P kþ1 ¼ P k þ D k ð15þ D k ¼ ðx k Þ T X k þ m k X k 1 ðx k Þ T S k ð16þ The Broyden method approximates X k by B k. It does not involve computing derivatives at all. Moreover, B k is obtained from B k 1 using simple procedures, so called, update. Broyden s good method (BGM) is defined by B kþ1 ¼ B k þ ðy k B k D k ÞD T k D T k D k ð17þ Here, y k ¼ SðP kþ1 Þ SðP k Þ is the difference in function values of the previous and present iterations. Additionally, the definition of Broyden s bad method (BBM) is B kþ1 ¼ B k þ ðy k B k D k Þy T k B k y T k B kd k ð18þ The name of the method does not indicate that BGM is superior to BBM. In [17], a combined method was devised to choose BGM or BBM according to the test rule as defined by jd T k D k 1j jd T k ðx kþ 1 y k j < jyt k y k 1j y T k y k ð19þ If the above condition is satisfied, BGM is applied, otherwise, BBM is selected.
10 534 K. W. KIM AND S. W. BAEK The Broyden combined update turned out to be superior to BGM or BBM alone. In order to use this method, the first jacobian value, B 0, should be provided using finite-difference (FD) approximation or AD. Inverse Analysis Procedure In order to estimate the radiative heat flux distribution over the heater surface satisfying the desired radiative heat flux and temperature on the design surface, all properties are assumed to be known except the radiative heat flux distribution over the heater surface. The computational algorithm for the inverse problem can be summarized as follows. Suppose that the desired radiative heat flux q R;d ¼ðq R;d;1 ; q R;d;2 ;...; q R;d;n e Þ and temperature are given at each element on the design surface and an initial guess P 0 is available. The solution is sought by the following iterative steps. 1. Knowing P k, compute I ð~r; _ s Þ in Eq. (1) by solving the direct problem with Eq. (4) for the heater surface and Eq. (3) for other boundaries. 2. Compute internal temperature distribution H from Eq. (8) and radiative heat flux on the design surface q R;d given by Eq. (5). 3. Check the stopping criterion. If max½ðp k P k 1 Þ=P k Š10 6, all steps are terminated, otherwise, continue the inverse procedure. 4. Compute sensitivity matrix X with methods of AD from Eq. (14) or BC from Eqs. (17), (18), and (19). 5. Knowing X k, compute P kþ1 using Eq. (13). Replace k by k þ 1 and return to step 1. RESULTS AND DISCUSSION Validation of FVM Code for Radiation FVM code was applied for validation to infinite concentric cylinders with isothermal black boundaries in radiative equilibrium, since results obtained by the Monte Carlo method [19] and the discrete ordinates method (DOM) [20] are available in the literature. At radiative equilibrium, the temperature field is computed from Eq. (8), given the ratio of inner to outer radius. In order to simulate infinite concentric cylinders, cyclic condition was imposed on the east and west boundaries, that is, outward intensity on the east boundary was set as inward intensity on the west boundary if ray direction is the same, and vice versa, respectively. Emissive power distribution for various ratios of inner to outer radius was computed in comparison with other results as follows: E b ¼ H4 ðrþ H 4 ðr out Þ H 4 ðr in Þ H 4 ðr out Þ ð20þ Figure 3 shows the results acquired by various methods: Monte Carlo, DOM, and the present FVM. The best agreement between FVM and Monte Carlo is
11 INVERSE RADIATION DESIGN PROBLEM 535 Figure 3. Emissive power distribution inside medium for various ratios of inner to outer radius (s ¼ 2, black walls). obtained at radius ratios of 0.01 and 0.1, while FVM is better than DOM. Through this comparison, the computational code made here for radiation could be validated for future analysis. Inverse Problem In order to simulate desired radiative heat flux, the direct problem was solved using Eqs. (1), (3), (5), and (8), given boundary conditions and medium properties as shown in Figure 1. After obtaining simulated desired radiative heat flux, inversely we tried to estimate boundary temperature over the heater surface to see if the inverse method works properly, and then radiative heat flux over the heater surface was also estimated with boundary condition of Eq. (4) over the heater surface. As shown in Figure 4, the boundary temperature over the heater surface is estimated very accurately, while the radiative heat flux, which satisfies the desired heat flux on the design surface, could also be found. The result has the same trend as in [10]. However, in this study the downward direction is defined as positive for radiative heat flux over the design and heater surface. Comparison of computational time for the sensitivity matrix was carried out to investigate the efficiency of AD and BC. Finite-difference approximation was adopted as reference, while the forward approximation was used as follows
12 536 K. W. KIM AND S. W. BAEK Figure 4. Estimated temperature and radiative heat flux over heater surface and desired radiative heat flux on design surface. (e ¼ 10 6 is used): qs i qp j ¼ qq e;i qp j q e;i ðp 1;...; P j þ ep j ;...; P np Þ q e;i ðp 1;...; P j ;...; P np Þ ep j for i ¼ 1;...; n e and j ¼ 1;...; n p ð21þ Maximum optimization level was used in the FORTRAN compiler, using a Pentium IV 3.0-GHz processor. All sensitivity methods satisfied the stopping criterion at the same iteration number, while BC needed the shortest computational time as shown in Table 1. The result shows that BC is efficient as well as accurate enough to estimate unknown parameters, so we adopted BC as the sensitivity method for future inverse analysis, while AD is used for B 0. If iteration number increases, the difference in computational time among sensitivity methods will also increase. Table 1. Comparison of iteration number and computational time for three cases of LMM þ FD, LMM þ AD, and LMM þ BC LMM þ FD LMM þ AD LMM þ BC (AD used for B 0 ) Initial value Iteration no. CPU time (s) Iteration no. CPU time (s) Iteration no. CPU time (s) , , ,061.39
13 INVERSE RADIATION DESIGN PROBLEM 537 Figure 5. Comparison of estimated radiative heat flux over heater surface satisfying desired condition for various damping parameters at k ¼ 2 and radiative heat flux on design surface in case of L=r out ¼ 2. Effect of Ratio of Cylinder Length to Outer Radius The effect of the ratio of cylinder length to outer radius on the radiative heat flux distribution was considered at various dimensionless positions defined as z ¼ z=l over the heater surface. When the ratio was smaller than the reference value, L=r out ¼ 4, some damping was necessary to obtain a smooth and stable solution. For this reason, in order to select an appropriate damping parameter, various values were applied, especially to the case of L=r out ¼ 2, as shown in Figure 5. The results for values of 0.01 and 0.1 satisfied the desired radiative heat flux over the design surface. However, since design engineers usually prefer a more uniform distribution, the value of 0.1 is selected as damping parameter. Figure 6 shows radiative heat flux distributions over the heater and design surface for the various ratios of cylinder length to outer radius. As the ratio increases, the radiative heat flux distribution over the heater surface becomes more uniform. This result indicates that the desired radiative heat flux over the design surface can be achieved with a heater with more uniform heat flux if the length of the heater surface is extended. Effect of Emissivity at Side Walls When emissivity at side walls 3 and 4 is changed, the variation of radiative heat flux over the heater surface which minimizes the difference between estimated and desired heat flux over the design surface, is shown in Figure 7. The result shows that the magnitude of radiative heat flux over the heater surface near the side walls becomes bigger than that for the reference value, e 3;4 ¼ 0:8, when the wall emissivity
14 538 K. W. KIM AND S. W. BAEK Figure 6. Estimated radiative heat flux distributions over heater surface satisfying desired condition for various ratios of cylinder length to outer radius and resulting radiative heat flux on design surface. Figure 7. Estimated radiative heat flux distributions over heater surface satisfying desired condition for various wall emissivities and resulting radiative heat flux on design surface.
15 INVERSE RADIATION DESIGN PROBLEM 539 increases. Even though the emissivity was significantly changed, corresponding inverse solutions satisfying the desired heat flux were found. Effect of Temperature at Side Walls When temperature at side walls 3 and 4 increases, the magnitude of radiative heat flux over the heater surface near the side walls becomes smaller than that for the reference value, H 3;4 ¼ 0:5, as shown in Figure 8. The effect is smaller than that of wall emissivity, because the sensitivity of wall emissivity is greater than that of wall temperature. In the case of H 3;4 ¼ 0:7, damping is needed to obtain the stable inverse solution at k ¼ 2. It is also found that inverse solutions which satisfy the desired condition are obtained even when side wall temperatures are varied. Effect of Phase Function If intensities are scattered toward the forward direction, more energy is emitted to the design surface. For this reason, less radiative heat flux from the heater will be necessary to satisfy the desired condition for forward scattering, and vice versa. The effect of phase function on radiative heat flux over the heater surface is the smallest compared with previous cases, as shown in Figure 9. This means that the sensitivity of phase function with respect to radiative heat flux over the design surface is the smallest. Figure 8. Estimated radiative heat flux distributions over heater surface satisfying desired condition for various side wall temperatures and resulting radiative heat flux on design surface.
16 540 K. W. KIM AND S. W. BAEK Figure 9. Estimated radiative heat flux distributions over heater surface satisfying desired condition for various phase functions and resulting radiative heat flux on design surface. Effect of Single Scattering Albedo Additionally, the effect of variation of single scattering albedo was negligible, because in the case of radiative equilibrium the RTE with isotropic scattering does not contain single scattering albedo. This implies that the absorption and scattering coefficients do not affect the radiative heat flux distribution over the heater surface. Example of a Real Design Problem The methodology used in this study was applied to a real design problem in which an inverse solution is not known a priori. On the design surface, two boundary conditions of uniform dimensionless temperature, H 1 ¼ 0:5, and radiative heat flux, q R;d ¼ 0:4 or 0.7, are imposed. In order to inspect convergence more carefully, a damping parameter of m ¼ 0:1 was applied. Figure 10 shows the inverse solution satisfying each imposed condition. At early iteration stage of k ¼ 2, the distribution is smooth owing to damping. However, as iteration number increases, nonsmooth variation occurs near the edge position. Resulting radiative heat flux distributions over the design surface calculated with estimated heat flux over the heater surface are shown in Figure 11. Desired uniform heat fluxes are achieved at the final stage. However, even with the early-staged smooth solutions for k ¼ 2, desired heat flux distributions could be achieved except at both end regions. This minor shortcoming can be easily eliminated though the extension of heater length, as discussed previously. Additionally, the increment in boundary temperature would reduce the dramatic increase of radiative heat flux at end regions over the heat surface.
17 INVERSE RADIATION DESIGN PROBLEM 541 Figure 10. Estimated radiative heat flux distributions over heater surface satisfying various desired conditions on design surface. Figure 11. Resulting radiative heat fluxes on design surface.
18 542 K. W. KIM AND S. W. BAEK CONCLUSIONS Inverse radiation design analysis in a two-dimensional, concentric, cylindrical, absorbing, emitting, and scattering medium has been conducted, given desired boundary conditions on the design surface. By adopting the automatic differentiation and the Broyden combined update as sensitivity method, the computational time for calculating the sensitivity matrix could be reduced up to one-third compared with that required for the finite-difference approximation. Various parametric analyses have shown that in most cases the radiative heat flux over the heater surface satisfying the desired condition over the design surface could be well estimated using the current inverse analysis. Stable and smooth solution could be obtained by controlling the damping parameter used in the Levenberg-Marqurdt method, even for ill-posed cases. Finally, the inverse method was applied to a real design problem to get the radiative heat flux over a heater surface satisfying uniform heat flux and temperature distribution on the design surface. Even with the early-staged smooth solutions, the desired heat flux distributions could be achieved except at both end regions. This minor trouble could be successfully eliminated though extension of the heater surface length. REFERENCES 1. M. N. Özisik and H. R. B. Orlande, Inverse Heat Transfer, chap. 1, Taylor & Francis, New York, H. Y. Li and C. Y. Yang, A Genetic Algorithm for Inverse Radiation Problems, Int. J. Heat Mass Transfer, vol. 40, pp , N. R. Ou and C. H. Wu, Simultaneous Estimation of Extinction Coefficient Distribution, Scattering Albedo and Phase Function of a Two-Dimensional Medium, Int. J. Heat Mass Transfer, vol. 45, pp , H. Y. Li, An Inverse Source Problem in Radiative Transfer for Spherical Media, Numer. Heat Transfer B, vol. 32, pp , H. M. Park and T. Y. Yoon, Solution of the Inverse Radiation Problem Using a Conjugate Gradient Method, Int. J. Heat Mass Transfer, vol. 43, pp , K. W. Kim, S. W. Baek, M. Y. Kim, and H. S. Ryu, Estimation of Emissivities in a Two- Dimensional Irregular Geometry by Inverse Radiation Analysis Using Hybrid Genetic Algorithm, J. Quant. Spectrosc. Radiat. Transfer, vol. 87, pp. 1 14, D. K. Alexander and H. H. Andreas, Optical Tomography Using the Time-Independent Equation of Radiative Transfer Part 2: Inverse Model, J. Quant. Spectrosc. Radiat. Transfer, vol. 72, pp , M. R. Jones, Inverse Analysis of Radiative Heat Transfer Systems, J. Heat Transfer, vol. 121, pp , J. R. Howell, O. A. Ezekoye, and J. C. Morales, Inverse Design Model for Radiative Heat Transfer, J. Heat Transfer, vol. 122, pp , S. S. M. Hosseini, S. H. Mansouri, and J. R. Howell, Inverse Boundary Design Radiation Problem in Absorbing-Emitting Media with Irregular Geometry, Numer. Heat Transfer A, vol. 43, pp , S. S. M. Hosseini and S. H. Mansouri, Inverse Design for Radiative Heat Source in Two-Dimensional Participating Media, Numer. Heat Transfer B, vol. 46, pp , 2004.
19 INVERSE RADIATION DESIGN PROBLEM K. W. Kim and S. W. Baek, Inverse Surface Radiation Analysis in an Axisymmetric Cylindrical Enclosure Using a Hybrid Genetic Algorithm, Numer. Heat Transfer A, vol. 46, pp , K. W. Kim and S. W. Baek, Comparison of Sensitivity Methods for the Function Estimation of Boundary Conditions by Inverse Radiation Analysis, 4th International Conference on Computational Heat and Mass Transfer, Paris-Cachan, France, vol. 2, pp , E. H. Chui, G. D. Raithby, and P. M. J. Hughes, Prediction of Radiative Transfer in Cylindrical Enclosure with the Finite Volume Method, J. Thermophys. Heat Transfer, vol. 6, pp , S. W. Baek and M. Y. Kim, Analysis of Radiative Heating of a Rocket Plume Base with the Finite-Volume Method, Int. J. Heat Mass Transfer, vol. 40, pp , C. Bischof, A. Carle, P. Khademi, A. Mauer, and P. Hovland, ADIFOR2.0 User s Guide (Revision D), Argonne National Lab. Tech. Memo. ANL=MCS-TM-192, and CRPC Tech. Rep. CRPC-TR95516-S, Argonne Natl. Lab., Argonne, IL, J. M. Martínez, Practical Quasi-Newton Methods for Solving Nonlinear Systems, J. Comput. Appl. Math, vol. 124, pp , J. E. Dennis and R. B. Schnabel, Numerical Methods for Unconstrained Optimization and Nonlinear Equations, chap. 8, Prentice-Hall, Englewood Cliffs, NJ, M. Perlmuttter and J. R. Howell, Radiant Transfer through a Gray Gas Between Concentric Cylinders Using Monte Carlo, J. Heat Transfer, vol. 86, pp , R. Vaillon, M. Lallemand, and D. Lemonnier, Radiative Heat Transfer in Orthogonal Curvilinear Coordinates Using the Discrete Ordinates Method, J. Quant. Spectrosc. Radiat. Transfer, vol. 55, pp. 7 17, 1996.
ARTICLE IN PRESS. Received 20 June 2007; received in revised form 14 November 2007; accepted 20 November 2007
Journal of Quantitative Spectroscopy & Radiative Transfer 09 (008) 607 6 www.elsevier.com/locate/jqsrt Radiative heat transfer between two concentric spheres separated by a two-phase mixture of non-gray
More informationSOLVING TRANSIENT CONDUCTION AND RADIATION USING FINITE VOLUME METHOD
SOLVING TRANSIENT CONDUCTION AND RADIATION USING FINITE VOLUME METHOD 1 PRERANA NASHINE, 2 ASHOK K SATAPATHY 1,2 National Institute of Technology Rourkela, Mechanical Engineering Department, India E-mail:
More informationHigh Altitude Rocket Plume and Thermal Radiation Analysis
High Altitude Rocket Plume and Thermal Radiation Analysis [ Woo Jin Jeon, Seung Wook Baek, Jae Hyun Park and Dong Sung Ha ] Abstract In this study, rocket plume behavior at various altitudes and radiative
More informationMasato Akamatsu a & Zhixiong Guo b a Graduate School of Science and Engineering, Yamagata University,
This article was downloaded by: [Rutgers University] On: 07 February 2013, At: 12:09 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office:
More informationInternational Journal of Heat and Mass Transfer
nternational Journal of Heat and Mass Transfer 52 (29) 15 111 Contents lists available at ScienceDirect nternational Journal of Heat and Mass Transfer journal homepage:.elsevier.com/locate/ijhmt Transient
More informationTechnology, Chonbuk National University, Jeonju, Chonbuk, Republic of Korea 3 Institute of Fluid Science, Tohoku University, Sendai, Japan
Numerical Heat Transfer, Part A, 60: 254 287, 2011 Copyright # Taylor & Francis Group, LLC ISSN: 1040-7782 print=1521-0634 online DOI: 10.1080/10407782.2011.588581 ANALYSIS OF CONDUCTION AND RADIATION
More informationAdaptability analysis of radiative transport diffusion approximation in planar-graded-index media
Research Article Adaptability analysis of radiative transport diffusion approximation in planar-graded-index media Advances in Mechanical Engineering 2018, ol. 10(11) 1 6 Ó The Author(s) 2018 DOI: 10.1177/1687814018809613
More informationULTRAFAST LASER PULSE TRAIN RADIATION TRANSFER IN A SCATTERING-ABSORBING 3D MEDIUM WITH AN INHOMOGENEITY
Heat Transfer Research 46(9), 861 879 (2015) ULTRAFAST LASER PULSE TRAIN RADIATION TRANSFER IN A SCATTERING-ABSORBING 3D MEDIUM WITH AN INHOMOGENEITY Masato Akamatsu 1,* & Zhixiong Guo 2 1 Graduate School
More informationA New and Simple Technique to Normalize the HG Phase Function for Conserving Scattered Energy and Asymmetry Factor
This article was downloaded by: [Rutgers University] On: 24 February 2014, At: 06:10 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office:
More informationPII S (99) SCALING ANISOTROPIC SCATTERING IN RADIATIVE TRANSFER IN THREE-DIMENSIONAL NONHOMOGENEOUS MEDIA
Pergamon Int. Comm. Heat Mass Transfer, Vol. 26, No. 7, pp. 997-1007, 1999 Copyright 1999 Elsevier Science Ltd Printed in the USA. All rights reserved 0735-1933/99IS-see front matter PII S0735-1933(99)00089-5
More informationOptimal heliostat system using hybrid type solar tracking device
Optimal heliostat system using hybrid type solar tracking device Dong Il Lee 1,a, Woo Jin Jeon 1,b, Seung Wook Baek 1*,c and Soon eung Chang,d 1 Propulsion and Combustion Laboratory, School of Mechanical,
More informationFull terms and conditions of use:
This article was downloaded by: [Rutgers University] On: 27 September 2012, At: 09:51 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office:
More informationAn Inverse Boundary Design Problem in a Radiant Enclosure
he 6 th. International Chemical Engineering Congress & Exhibition IChEC 009 6 0 November 009, ish Island, Iran An Inverse Boundary Design Problem in a Radiant Enclosure A. arimipour *, S.M.H. Sarvari,
More informationAn inverse problem in estimating simultaneously the effective thermal conductivity and volumetric heat capacity of biological tissue
Applied Mathematical Modelling 31 (7) 175 1797 www.elsevier.com/locate/apm An inverse problem in estimating simultaneously the effective thermal conductivity and volumetric heat capacity of biological
More informationIndo-German Winter Academy
Indo-German Winter Academy - 2007 Radiation in Non-Participating and Participating Media Tutor Prof. S. C. Mishra Technology Guwahati Chemical Engineering Technology Guwahati 1 Outline Importance of thermal
More informationImprovement of computational time in radiative heat transfer of three-dimensional participating media using the radiation element method
Journal of Quantitative Spectroscopy & Radiative Transfer 73 (2002) 239 248 www.elsevier.com/locate/jqsrt Improvement of computational time in radiative heat transfer of three-dimensional participating
More informationLevenberg-Marquardt Method for Solving the Inverse Heat Transfer Problems
Journal of mathematics and computer science 13 (2014), 300-310 Levenberg-Marquardt Method for Solving the Inverse Heat Transfer Problems Nasibeh Asa Golsorkhi, Hojat Ahsani Tehrani Department of mathematics,
More informationNumerical solution of hyperbolic heat conduction in thin surface layers
International Journal of Heat and Mass Transfer 50 (007) 9 www.elsevier.com/locate/ijhmt Numerical solution of hyperbolic heat conduction in thin surface layers Tzer-Ming Chen * Department of Vehicle Engineering,
More informationSpectral gradient projection method for solving nonlinear monotone equations
Journal of Computational and Applied Mathematics 196 (2006) 478 484 www.elsevier.com/locate/cam Spectral gradient projection method for solving nonlinear monotone equations Li Zhang, Weijun Zhou Department
More informationResearch Article A Two-Step Matrix-Free Secant Method for Solving Large-Scale Systems of Nonlinear Equations
Applied Mathematics Volume 2012, Article ID 348654, 9 pages doi:10.1155/2012/348654 Research Article A Two-Step Matrix-Free Secant Method for Solving Large-Scale Systems of Nonlinear Equations M. Y. Waziri,
More informationA COMPARISON OF THE DISCRETE ORDINATES METHOD AND FINITE VOLUME METHOD FOR RADIATIVE HEAT TRANSFER ANALYSIS
Proceedings of the ASME 2011 International Mechanical Engineering Congress & Exposition IMECE2011 November 11-17, 2011, Denver, Colorado, USA IMECE2011-63658 A COMPARISON OF THE DISCRETE ORDINATES METHOD
More informationAnalysis of Scattering of Radiation in a Plane-Parallel Atmosphere. Stephanie M. Carney ES 299r May 23, 2007
Analysis of Scattering of Radiation in a Plane-Parallel Atmosphere Stephanie M. Carney ES 299r May 23, 27 TABLE OF CONTENTS. INTRODUCTION... 2. DEFINITION OF PHYSICAL QUANTITIES... 3. DERIVATION OF EQUATION
More informationDiscrete-ordinates solution of short-pulsed laser transport in two-dimensional turbid media
Discrete-ordinates solution of short-pulsed laser transport in two-dimensional turbid media Zhixiong Guo and Sunil Kumar The discrete-ordinates method is formulated to solve transient radiative transfer
More informationResearch Article Comparison of Two Models for Radiative Heat Transfer in High Temperature Thermal Plasmas
Modelling and Simulation in Engineering Volume 2, Article ID 2858, 7 pages doi:.55/2/2858 Research Article Comparison of Two Models for Radiative Heat Transfer in High Temperature Thermal Plasmas Matthieu
More informationA new modified Halley method without second derivatives for nonlinear equation
Applied Mathematics and Computation 189 (2007) 1268 1273 www.elsevier.com/locate/amc A new modified Halley method without second derivatives for nonlinear equation Muhammad Aslam Noor *, Waseem Asghar
More informationIntroduction to modeling of thermal radiation in participating gases
Project Report 2008 MVK 160 Heat and Mass Transfer May 07, 2008, Lund, Sweden Introduction to modeling of thermal radiation in participating gases Eric Månsson Dept. of Energy Sciences, Faculty of Engineering,
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS NUMERICAL FLUID MECHANICS FALL 2011
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 02139 2.29 NUMERICAL FLUID MECHANICS FALL 2011 QUIZ 2 The goals of this quiz 2 are to: (i) ask some general
More informationImproved Treatment of Anisotropic Scattering in Radiation Transfer Analysis Using the Finite Volume Method
Heat Transfer Engineering ISSN: 0145-7632 (Print) 1521-0537 (Online) Journal homepage: http://www.tandfonline.com/loi/uhte20 Improved Treatment of Anisotropic Scattering in Radiation Transfer Analysis
More informationTo remain at 0 K heat absorbed by the medium must be removed in the amount of. dq dx = q(l) q(0) dx. Q = [1 2E 3 (τ L )] σ(t T 4 2 ).
294 RADIATIVE HEAT TRANSFER 13.7 Two infinite, isothermal plates at temperatures T 1 and T 2 are separated by a cold, gray medium of optical thickness τ L = L (no scattering). (a) Calculate the radiative
More informationInverse Heat Flux Evaluation using Conjugate Gradient Methods from Infrared Imaging
11 th International Conference on Quantitative InfraRed Thermography Inverse Heat Flux Evaluation using Conjugate Gradient Methods from Infrared Imaging by J. Sousa*, L. Villafane*, S. Lavagnoli*, and
More informationHT NORMALIZATION FOR ULTRAFAST RADIATIVE TRANSFER ANALYSIS WITH COLLIMATED IRRADIATION
Proceedings of the ASME 2012 Summer Heat Transfer Conference HT2012 July 8-12, 2012, Rio Grande, Puerto Rico HT2012-58307 NORMALIZATION FOR ULTRAFAST RADIATIVE TRANSFER ANALYSIS WITH COLLIMATED IRRADIATION
More informationANALYSIS OF NATURAL CONVECTION IN HORIZONTAL CONCENTRIC ANNULI OF VARYING INNER SHAPE
Numerical Heat Transfer, Part A, 68: 1155 1174, 2015 Copyright # Taylor & Francis Group, LLC ISSN: 1040-7782 print=1521-0634 online DOI: 10.1080/10407782.2015.1032016 ANALYSIS OF NATURAL CONVECTION IN
More informationEquivalent isotropic scattering formulation for transient
Equivalent isotropic scattering formulation for transient short-pulse radiative transfer in anisotropic scattering planar media Zhixiong Guo and Sunil Kumar An isotropic scaling formulation is evaluated
More informationSimulation of Nanoscale Multidimensional Transient Heat Conduction Problems Using Ballistic-Diffusive Equations and Phonon Boltzmann Equation
Ronggui Yang Mem. ASME Gang Chen Mem. ASME Marine Laroche 1 Mechanical Engineering Department, Massachusetts Institute of Technology, Cambridge, MA 02139-4307 Yuan Taur Department of Electrical & Computer
More informationEFFECTS OF GRAIN SIZE ON THE SPECTRAL ENERGY DISTRIBUTION OF DUSTY CIRCUMSTELLAR ENVELOPES
The Astrophysical Journal, 604:238 251, 2004 March 20 # 2004. The American Astronomical Society. All rights reserved. Printed in U.S.A. EFFECTS OF GRAIN SIZE ON THE SPECTRAL ENERGY DISTRIBUTION OF DUSTY
More informationRadiative Equilibrium Models. Solar radiation reflected by the earth back to space. Solar radiation absorbed by the earth
I. The arth as a Whole (Atmosphere and Surface Treated as One Layer) Longwave infrared (LWIR) radiation earth to space by the earth back to space Incoming solar radiation Top of the Solar radiation absorbed
More informationSOLUTION OF THE INTEGRATED RADIATIVE TRANSFER EQUATION FOR GRAY AND NONGRAY MEDIA
Numerical Heat Transfer, Part A, 50: 205 228, 2006 Copyright # Taylor & Francis Group, LLC ISSN: 1040-7782 print=1521-0634 online DOI: 10.1080/10407780600602465 SOLUTION OF THE INTEGRATED RADIATIVE TRANSFER
More informationEffective Borehole Thermal Resistance of A Single U-Tube Ground Heat Exchanger
Numerical Heat Transfer, Part A: Applications An International Journal of Computation and Methodology ISSN: 1040-7782 (Print) 1521-0634 (Online) Journal homepage: http://www.tandfonline.com/loi/unht20
More informationHakwoon Kim Gunhee Jang Sanghoon Lee. 1 Introduction
Microsyst Technol (2011) 17:749 759 DOI 10.1007/s00542-010-1188-4 TECHNICAL PAPER Complete determination of the dynamic coefficients of coupled journal and thrust bearings considering five degrees of freedom
More informationDevelopment of New Solid Angle Quadrature Sets to Satisfy Even- and Odd-Moment Conditions
Journal of Nuclear Science and Technology ISSN: 22-11 (Print) 1881-1248 (Online) Journal homepage: http://www.tandfonline.com/loi/tnst2 Development of New Solid Angle Quadrature Sets to Satisfy Even- and
More informationDevelopment of the numerical method for calculating sound radiation from a rotating dipole source in an opened thin duct
Journal of Sound and Vibration 295 (2006) 739 752 JOURNAL OF SOUND AND VIBRATION www.elsevier.com/locate/jsvi Development of the numerical method for calculating sound radiation from a rotating dipole
More informationIn this method, one defines
1 Feautrier s Method for Radiative Transfer This project is to use Feautrier s method to solve the Gray Atmosphere problem for a plane-parallel atmosphere Although we have developed alternate techniques
More informationRigid body modeling issue in acoustical topology optimization
Purdue University Purdue e-pubs Birck and NCN Publications Birck Nanotechnology Center 2009 Rigid body modeling issue in acoustical topology optimization Jin Woo Lee Purdue University - Main Campus, lee04@purdue.edu
More informationPROBLEM (a) Long duct (L): By inspection, F12. By reciprocity, (b) Small sphere, A 1, under concentric hemisphere, A 2, where A 2 = 2A
PROBLEM 3. KNON: Various geometric shapes involving two areas and. FIND: Shape factors, F and F, for each configuration. SSUMPTIONS: Surfaces are diffuse. NLYSIS: The analysis is not to make use of tables
More informationExperimental designs for precise parameter estimation for non-linear models
Minerals Engineering 17 (2004) 431 436 This article is also available online at: www.elsevier.com/locate/mineng Experimental designs for precise parameter estimation for non-linear models Z. Xiao a, *,
More informationLight transport in biological tissue based on the simplified spherical harmonics equations
Journal of Computational Physics 220 (2006) 44 470 www.elsevier.com/locate/jcp Light transport in biological tissue based on the simplified spherical harmonics equations Alexander D. Klose a, *, Edward
More informationPHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.
PHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.4) We will consider two cases 1. f(x) = 0 1-dimensional 2. f(x) = 0 d-dimensional
More informationMaximum time-resolved hemispherical reflectance of absorbing and isotropically scattering media
Journal of Quantitative Spectroscopy & Radiative Transfer 14 (27) 384 399 www.elsevier.com/locate/jqsrt Maximum time-resolved hemispherical reflectance of absorbing and isotropically scattering media Kyle
More information12 Thermal Radiation including Participating Media
Thermal Radiation including Participating Media Nomenclature A area [m ] c speed of light c specific heat [J/(kgK)] p E emissive power [W/m ] F view factor [-] ij G incident radiation [W/m ] h Planck's
More informationImproved Damped Quasi-Newton Methods for Unconstrained Optimization
Improved Damped Quasi-Newton Methods for Unconstrained Optimization Mehiddin Al-Baali and Lucio Grandinetti August 2015 Abstract Recently, Al-Baali (2014) has extended the damped-technique in the modified
More informationHandout on Newton s Method for Systems
Handout on Newton s Method for Systems The following summarizes the main points of our class discussion of Newton s method for approximately solving a system of nonlinear equations F (x) = 0, F : IR n
More informationPROBLEM (a) Long duct (L): By inspection, F12. By reciprocity, (b) Small sphere, A 1, under concentric hemisphere, A 2, where A 2 = 2A
PROBLEM 3. KNON: Various geometric shapes involving two areas and. FIND: Shape factors, F and F, for each configuration. SSUMPTIONS: Surfaces are diffuse. NLYSIS: The analysis is not to make use of tables
More informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationMCRT: L4 A Monte Carlo Scattering Code
MCRT: L4 A Monte Carlo Scattering Code Plane parallel scattering slab Optical depths & physical distances Emergent flux & intensity Internal intensity moments Constant density slab, vertical optical depth
More informationHot Air Gun Identification by Inverse Heat Transfer
Vol. 8/ No. 4/ Winter 2016 pp. 29-34 Hot Air Gun Identification by Inverse Heat Transfer A. M. Tahsini *1 and S. Tadayon Mousavi 2 1, 2. Aerospace Research Institute; Ministry of Science, Research and
More informationINVERSE DETERMINATION OF SPATIAL VARIATION OF DIFFUSION COEFFICIENTS IN ARBITRARY OBJECTS CREATING DESIRED NON- ISOTROPY OF FIELD VARIABLES
Contributed Papers from Materials Science and Technology (MS&T) 05 October 4 8, 05, Greater Columbus Convention Center, Columbus, Ohio, USA Copyright 05 MS&T5 INVERSE DETERMINATION OF SPATIAL VARIATION
More informationInverse Problems and Optimal Design in Electricity and Magnetism
Inverse Problems and Optimal Design in Electricity and Magnetism P. Neittaanmäki Department of Mathematics, University of Jyväskylä M. Rudnicki Institute of Electrical Engineering, Warsaw and A. Savini
More informationMath 302 Outcome Statements Winter 2013
Math 302 Outcome Statements Winter 2013 1 Rectangular Space Coordinates; Vectors in the Three-Dimensional Space (a) Cartesian coordinates of a point (b) sphere (c) symmetry about a point, a line, and a
More informationNumerical Solution of an Inverse Diffusion Problem
Applied Mathematical Sciences, Vol. 1, 2007, no. 18, 863-868 Numerical Solution of an Inverse Diffusion Problem Y. Ahmadizadeh Physics Dept., Faculty of Sciences Air University, Postal Code 1384663113,
More informationCost-accuracy analysis of a variational nodal 2D/1D approach to pin resolved neutron transport
Cost-accuracy analysis of a variational nodal 2D/1D approach to pin resolved neutron transport ZHANG Tengfei 1, WU Hongchun 1, CAO Liangzhi 1, LEWIS Elmer-E. 2, SMITH Micheal-A. 3, and YANG Won-sik 4 1.
More information5. FVM discretization and Solution Procedure
5. FVM discretization and Solution Procedure 1. The fluid domain is divided into a finite number of control volumes (cells of a computational grid). 2. Integral form of the conservation equations are discretized
More informationA New Approach for Solving Dual Fuzzy Nonlinear Equations Using Broyden's and Newton's Methods
From the SelectedWorks of Dr. Mohamed Waziri Yusuf August 24, 22 A New Approach for Solving Dual Fuzzy Nonlinear Equations Using Broyden's and Newton's Methods Mohammed Waziri Yusuf, Dr. Available at:
More informationPOTENTIAL THEORY AND HEAT CONDUCTION DIRICHLET S PROBLEM
Chapter 6 POTENTIAL THEORY AND HEAT CONDUCTION DIRICHLET S PROBLEM M. Ragheb 9/19/13 6.1 INTRODUCTION The solution of the Dirichlet problem is one of the easiest approaches to grasp using Monte Carlo methodologies.
More informationRESEARCH ARTICLE. A strategy of finding an initial active set for inequality constrained quadratic programming problems
Optimization Methods and Software Vol. 00, No. 00, July 200, 8 RESEARCH ARTICLE A strategy of finding an initial active set for inequality constrained quadratic programming problems Jungho Lee Computer
More informationNUMERICAL SIMULATION OF THREE-DIMENSIONAL COMBINED CONVECTIVE RADIATIVE HEAT TRANSFER IN RECTANGULAR CHANNELS
NUMERICAL SIMULATION OF THREE-DIMENSIONAL COMBINED CONVECTIVE RADIATIVE HEAT TRANSFER IN RECTANGULAR CHANNELS A Dissertation by MIN SEOK KO Submitted to the Office of Graduate Studies of Texas A&M University
More informationModeling of Advanced Melting Zone for Manufacturing of Optical Fibers*
Zhiyong Wei e-mail: gte384w@prism.gatech.edu Kok-Meng Lee kokmeng.lee@me.gatech.edu The George W. Woodruff School of Mechanical Engineering, Georgia Institute of Technology, Atlanta, GA 30332-0405 Zhi
More informationA novel technique of friction and material property measurement by tip test in cold forging
A novel technique of friction and material property measurement by tip test in cold forging Y T Im*, S H Kang, and J S Cheon Department of Mechanical Engineering, Korea Advanced Institute of Science and
More informationA Boundary Condition for Porous Electrodes
Electrochemical Solid-State Letters, 7 9 A59-A63 004 0013-4651/004/79/A59/5/$7.00 The Electrochemical Society, Inc. A Boundary Condition for Porous Electrodes Venkat R. Subramanian, a, *,z Deepak Tapriyal,
More informationA derivative-free nonmonotone line search and its application to the spectral residual method
IMA Journal of Numerical Analysis (2009) 29, 814 825 doi:10.1093/imanum/drn019 Advance Access publication on November 14, 2008 A derivative-free nonmonotone line search and its application to the spectral
More informationA Second Order Radiative Transfer Equation and Its Solution by Meshless. Method with Application to Strongly Inhomogeneous Media
A Second Order Radiative Transfer Equation and Its Solution by Meshless Method with Application to Strongly Inhomogeneous Media J.M. Zhao a, J.Y. Tan b, L.H. Liu a,b * a School of Energy Science and Engineering,
More informationTHERMOELASTIC ANALYSIS OF THICK-WALLED FINITE-LENGTH CYLINDERS OF FUNCTIONALLY GRADED MATERIALS
Journal of Thermal Stresses, 28: 391 408, 2005 Copyright # Taylor & Francis Inc. ISSN: 0149-5739 print/1521-074x online DOI: 10.1080/01495730590916623 THERMOELASTIC ANALYSIS OF THICK-WALLED FINITE-LENGTH
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat 16 (2011) 2730 2736 Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns Homotopy analysis method
More informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
More informationSpontaneous magnetization of the square 2D Ising lattice with nearest- and weak next-nearest-neighbour interactions
Phase Transitions Vol. 82, No. 2, February 2009, 191 196 Spontaneous magnetization of the square 2D Ising lattice with nearest- and weak next-nearest-neighbour interactions H.J.W. Zandvliet a * and C.
More informationp(θ,φ,θ,φ) = we have: Thus:
1. Scattering RT Calculations We come spinning out of nothingness, scattering stars like dust. - Jalal ad-din Rumi (Persian Poet, 1207-1273) We ve considered solutions to the radiative transfer equation
More informationEffects of Radiative Transfer Modeling on Transient Temperature Distribution in Semitransparent Glass Rod
Zhiyong Wei e-mail: gte384w@prism.gatech.edu Kok-Meng Lee e-mail: kokmeng.lee@me.gatech.edu Serge W. Tchikanda The George W. Woodruff School of Mechanical Engineering, Georgia Institute of Technology,
More information2 The Radiative Transfer Equation
9 The Radiative Transfer Equation. Radiative transfer without absorption and scattering Free space or homogeneous space I (r,,) I (r,,) r -r d da da Figure.: Following a pencil of radiation in free space
More informationA Trust-region-based Sequential Quadratic Programming Algorithm
Downloaded from orbit.dtu.dk on: Oct 19, 2018 A Trust-region-based Sequential Quadratic Programming Algorithm Henriksen, Lars Christian; Poulsen, Niels Kjølstad Publication date: 2010 Document Version
More informationPROBLEM L. (3) Noting that since the aperture emits diffusely, I e = E/π (see Eq ), and hence
PROBLEM 1.004 KNOWN: Furnace with prescribed aperture and emissive power. FIND: (a) Position of gauge such that irradiation is G = 1000 W/m, (b) Irradiation when gauge is tilted θ d = 0 o, and (c) Compute
More informationStars AS4023: Stellar Atmospheres (13) Stellar Structure & Interiors (11)
Stars AS4023: Stellar Atmospheres (13) Stellar Structure & Interiors (11) Kenneth Wood, Room 316 kw25@st-andrews.ac.uk http://www-star.st-and.ac.uk/~kw25 What is a Stellar Atmosphere? Transition from dense
More informationThe Radiative Transfer Equation
The Radiative Transfer Equation R. Wordsworth April 11, 215 1 Objectives Derive the general RTE equation Derive the atmospheric 1D horizontally homogenous RTE equation Look at heating/cooling rates in
More informationA far-field based T-matrix method for three dimensional acoustic scattering
ANZIAM J. 50 (CTAC2008) pp.c121 C136, 2008 C121 A far-field based T-matrix method for three dimensional acoustic scattering M. Ganesh 1 S. C. Hawkins 2 (Received 14 August 2008; revised 4 October 2008)
More informationFinite Element Method (FEM)
Finite Element Method (FEM) The finite element method (FEM) is the oldest numerical technique applied to engineering problems. FEM itself is not rigorous, but when combined with integral equation techniques
More informationRadiation in the atmosphere
Radiation in the atmosphere Flux and intensity Blackbody radiation in a nutshell Solar constant Interaction of radiation with matter Absorption of solar radiation Scattering Radiative transfer Irradiance
More informationRadiation from planets
Chapter 4 Radiation from planets We consider first basic, mostly photometric radiation parameters for solar system planets which can be easily compared with existing or future observations of extra-solar
More informationInvestigation of radiative heat transfer in fixed bed biomass furnaces
Available online at www.sciencedirect.com Fuel 87 (2008) 2141 2153 www.fuelfirst.com Investigation of radiative heat transfer in fixed bed biomass furnaces T. Klason a, X.S. Bai a, *, M. Bahador b, T.K.
More informationProgramming, numerics and optimization
Programming, numerics and optimization Lecture C-3: Unconstrained optimization II Łukasz Jankowski ljank@ippt.pan.pl Institute of Fundamental Technological Research Room 4.32, Phone +22.8261281 ext. 428
More informationA Novel and Precise Sixth-Order Method for Solving Nonlinear Equations
A Novel and Precise Sixth-Order Method for Solving Nonlinear Equations F. Soleymani Department of Mathematics, Islamic Azad University, Zahedan Branch, Zahedan, Iran E-mail: fazl_soley_bsb@yahoo.com; Tel:
More informationANALYSIS OF SOLIDIFICATION OF A SEMITRANSPARENT PLANAR LAYER USING THE LATTICE BOLTZMANN METHOD AND THE DISCRETE TRANSFER METHOD
Numerical Heat Transfer, Part A, 49: 279 299, 2006 Copyright# Taylor & Francis LLC ISSN: 1040-7782 print=1521-0634 online DOI: 10.1080/10407780500359828 ANALYSIS OF SOLIDIFICATION OF A SEMITRANSPARENT
More informationExperimental confirmation of the negentropic character of the diffraction polarization of diffuse radiation
Experimental confirmation of the negentropic character of the diffraction polarization of diffuse radiation V. V. Savukov In the course of analyzing the axiomatic principles on which statistical physics
More informationDocumentation of the Solutions to the SFPE Heat Transfer Verification Cases
Documentation of the Solutions to the SFPE Heat Transfer Verification Cases Prepared by a Task Group of the SFPE Standards Making Committee on Predicting the Thermal Performance of Fire Resistive Assemblies
More informationHT14 A COMPARISON OF DIFFERENT PARAMETER ESTIMATION TECHNIQUES FOR THE IDENTIFICATION OF THERMAL CONDUCTIVITY COMPONENTS OF ORTHOTROPIC SOLIDS
Inverse Problems in Engineering: heory and Practice rd Int. Conference on Inverse Problems in Engineering June -8, 999, Port Ludlow, WA, USA H4 A COMPARISON OF DIFFEREN PARAMEER ESIMAION ECHNIQUES FOR
More informationThe detection of subsurface inclusions using internal measurements and genetic algorithms
The detection of subsurface inclusions using internal measurements and genetic algorithms N. S. Meral, L. Elliott2 & D, B, Ingham2 Centre for Computational Fluid Dynamics, Energy and Resources Research
More informationFull terms and conditions of use:
This article was downloaded by:[rollins, Derrick] [Rollins, Derrick] On: 26 March 2007 Access Details: [subscription number 770393152] Publisher: Taylor & Francis Informa Ltd Registered in England and
More informationIII. TRANSFORMATION RELATIONS
III. TRANSFORMATION RELATIONS The transformation relations from cartesian coordinates to a general curvilinear system are developed here using certain concepts from differential geometry and tensor analysis,
More informationPartial differential equations
Partial differential equations Many problems in science involve the evolution of quantities not only in time but also in space (this is the most common situation)! We will call partial differential equation
More informationTests of relative earthquake location techniques using synthetic data
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 110,, doi10.1029/2004jb003380, 2005 Tests of relative earthquake location techniques using synthetic data Guoqing Lin and Peter Shearer Institute of Geophysics and
More informationNumerical Optimization Professor Horst Cerjak, Horst Bischof, Thomas Pock Mat Vis-Gra SS09
Numerical Optimization 1 Working Horse in Computer Vision Variational Methods Shape Analysis Machine Learning Markov Random Fields Geometry Common denominator: optimization problems 2 Overview of Methods
More information1 The formation and analysis of optical waveguides
1 The formation and analysis of optical waveguides 1.1 Introduction to optical waveguides Optical waveguides are made from material structures that have a core region which has a higher index of refraction
More informationNumerical Methods I Solving Nonlinear Equations
Numerical Methods I Solving Nonlinear Equations Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 October 16th, 2014 A. Donev (Courant Institute)
More information