( ) ( ) ( ) ( ) ( ) ( ) ( )
|
|
- Horace Richardson
- 6 years ago
- Views:
Transcription
1 Please cite this article as: Ewa Machrzak, Łukasz Turchan, Boundary element method for 3D Fourier-Kirchhoff heat transfer equation, Scientific Research of the Institute of Mathematics and Computer Science, 00, Volume 9, Issue, pages -30. The website: Scientific Research of the Institute of Mathematics and Computer Science BOUNDARY ELEMENT METHOD FOR 3D FOURIER-KIRCHHOFF HEAT TRANSFER EQUATION Ewa Machrzak,, Łukasz Turchan Department of Strength of Materials and Computational Mechanics Silesian University of Technology, Poland Institute of Mathematics, Czestochowa University of Technology, Poland Abstract. The 3D heat transfer problem (steady state) is considered. The equation describing the thermal processes contains the convective term (substantial derivative). The problem is solved by means of the boundary element method. The numerical model for constant boundary elements and constant internal cells is presented. In the final part of the paper the examples of computations are shown. The numerical results obtained by means of the BEM are compared with analytical solution and the very good compatibility can be observed.. Formulation of the problem The following Fourier-Kirchhoff equation (steady state, 3D problem) is considered x : λ T x cu T x + Q x = 0 () where λ is the thermal conductivity and c is the volumetric specific heat, respectively u = [u, u, u 3 ] is the velocity, Q(x) is the source function, T denotes the temperature and x = [x, x, x 3 ] are the spatial co-ordinates. For u = [u, 0, 0] the equation () takes a form T x x T x cu + Q x = The boundary conditions in the form : λ 0 x x : T x = T T x x : q( x) = λ = q n are also known, at the same time T/ n denotes the normal derivative, n = = [cos α, cos α, cos α 3 ] is the normal outward vector, T b, q b are the known boundary temperature and boundary heat flux, respectively. To solve the problem (), (3) the boundary element method is proposed. b b () (3)
2 E. Machrzak, Ł. Turchan. Boundary element method At first, the weighted residual criterion [, ] for equation () is formulated T x λ T ( x) cu + Q( x) T ( ξ, x) d = 0 x (4) where ξ is the observation point and T*( ξ, x) is the fundamental solution. The integral (4) is substituted by a sum of two integrals, while the first of them is transformed using the nd Green formula λ ξ, d λ ξ, d T x T x = T x T x + ( ξ, ) T x T x λt ( ξ, x) λ T ( x) d n n The next integral in formula (4) is integrated by parts T x cu T ξ, x d cu T ξ, x T x cosα d = x T cu x ( ξ, x) T x d (5) (6) Introducing (5), (6)into (4) one has ( ξ, ) ( x) ( ξ, x) + T + λ T ( ξ, x) cu T ( x) d x T x T x λt ( ξ, x) λ T ( x) cut ( ξ, x) T ( x) cos α d + n n Q x T ξ, d = 0 Fundamental solution T*( ξ, x) should fulfil the following equation ( ξ, x) T x T x + cu = x : λ ( ξ, ) x δ ξ, (7) (8) where δ(ξ, x) is the Dirac function. Taking into account the property (8) the equation (7) takes a form
3 Boundary element method for 3D Fourier-Kirchhoff heat transfer equation 3 where T ξ + T ξ, x q x d = q ξ, x cu T ξ, x cos α T x d + Q x T ξ, x d ( ξ, ) T x T x q( x) = λ, q ( ξ, x) = λ n n (9) (0) For ξ the boundary integral equation is obtained B ξ T ξ + T ξ, x q x d = q ξ, x cu T ξ, x cos α T x d + Q x T ξ, x d () where B(ξ) (0, ) is the coefficient connected with the location of point ξ on the boundary. For the problem considered the fundamental solution is the following [3] cu T ( ξ, x) = exp r + ( x ξ ) 4 π λ r λ () where r is the distance between the points ξ = (ξ, ξ, ξ 3 ) and x = (x, x, x 3 ). Using formula () the heat flux q*(ξ, x) resulting from fundamental solution can be calculated cu q ( ξ, x) = exp 3 r + ( x ξ ) 4 π r λ cu cu d + r + r λ λ cosα (3) where d = x ξ cosα + x ξ cosα + x ξ cosα (4) Numerical realization To solve equation () the boundary is divided into N boundary elements and the interior is divided into L internals cells. Next, the integrals appearing in () are substituted by the sums of integrals.
4 4 E. Machrzak, Ł. Turchan So, for optional boundary point ξ i one has N N i i i B ξ T ξ + T ξ, x q x d = = L i i l i q ξ, x cu T ξ, x cosα T x d + Q x T ξ, x dl = l= l (5) When the constant boundary elements and constant internal cells are used, this means x : ( ) ( ) T x = T x = T q x = q x = q (6) and l x : Q x = Q x = Q (7) l l then equation (5) takes form i T + q T ( ξ, x) = N i d = L i i i N T q ξ, x cu T ξ, x cosα d + Ql T ξ, x dl = l= l (8) or where and while N N L Gi, q = ( Hi, Ui, ) T + Pi, lql, i =,,, N (9) = = l= H i, i Gi, = T ξ, x d (0) i q ξ, x d, i = i q ( ξ, x) d, i = i ( ξ, ) i, () U = c u T x cosα d ()
5 Boundary element method for 3D Fourier-Kirchhoff heat transfer equation 5 and i Pi, l = T ξ, x dl (3) l The system of equations (9) allows one to determine the missing boundary values T and q. Next, the temperatures at the internal points ξ i, i = N +, N +,...,N + L can be calculated using the formula N N L (,, ),, (4) T = H U T G q + P Q i i i i i l l = = l= 4. Discretization of boundary and interior of cuboid The cuboid of dimension l x l x l 3 is considered (Fig. ). The boundary is divided into N constant boundary elements. It is assumed that hs = l / n, n = l / hs, n 3 = l 3 / hs, and then N = n n + n n 3 + n n 3. Below the fragment of Delphi code in which the boundary nodes coordinates are defined is shown. {przod} for i:= to n*n3 do y[i]:=0; for k:=0 to n3- do for i:= to n do x[i+n*k]:=i*hs-hs/.0; z[i+n*k]:=k*hs+hs/.0; {prawa} for i:=n*n3+ to n*n3+n*n3 do x[i]:=n*hs; for k:=0 to n3- do for := to n do y[+n*k+n*n3]:=*hs-hs/.0; z[+n*k+n*n3]:=k*hs+hs/.0; {tyl} for i:=n*n3+n*n3+ to n*n3+*n*n3 do y[i]:=n*hs; for k:=0 to n3- do for i:= to n do x[i+n*k+n*n3+n*n3]:=n*hs-(i*hs-hs/.0); z[i+n*k+n*n3+n*n3]:=k*hs+hs/.0; {lewa} for i:=*n*n3+n*n3+ to *n*n3+*n*n3 do x[i]:=0.0; for k:=0 to n3- do for := to n do y[+n*k+*n*n3+n*n3]:=n*hs-(*hs-hs/.0);
6 6 E. Machrzak, Ł. Turchan z[+n*k+*n*n3+n*n3]:=k*hs+hs/.0; {dol} for i:=*n*n3+*n*n3+ to *n*n3+*n*n3+n*n do z[i]:=0.0; for :=0 to n- do for i:= to n do x[i+n*+*n*n3+*n*n3]:=i*hs-hs/.0; y[i+n*+*n*n3+*n*n3]:=*hs+hs/.0; {gora} for i:=*n*n3+*n*n3+n*n+ to N do z[i]:=n3*hs; for :=0 to n- do for i:= to n do x[i+n*+*n*n3+*n*n3+n*n]:=i*hs-hs/.0; y[i+n*+*n*n3+*n*n3+n*n]:=*hs+hs/.0; The coordinates of internal nodes can be determined as follows: {wewnetrzne} for k:=0 to n3- do for :=0 to n- do for i:= to n do x[n+i+n*+n*n*k]:=i*hs-hs/.0; y[n+i+n*+n*n*k]:=*hs+hs/.0; z[n+i+n*+n*n*k]:=k*hs+hs/.0; To calculate the integrals (0), () and () the Gaussian cubatures method is used [, ]. In this method the coordinates of vertexes of each boundary element (quadrilateral) should be known. If the vertexes of quadrilateral one denotes by (x, y, z ), (x, y, z ), (x 3, y 3, z 3 ), (x 4, y 4, z 4 ) then on the basis of boundary node coordinates x[], y[], z[] they can be defined as follows: if <= n*n3 then {przod} x:=x[]-hs/; x:=x[]+hs/; x3:=x[]+hs/; x4:=x[]-hs/; y:=0.0; y:=0.0; y3:=0.0; y4:=0.0; z:=z[]-hs/; z:=z[]-hs/; z3:=z[]+hs/; z4:=z[]+hs/; end else if <= n*n3+n*n3 then {prawa} x:=l; x:=l; x3:=l; x4:=l; y:=y[]-hs/; y:=y[]+hs/; y3:=y[]+hs/; y4:=y[]-hs/;
7 Boundary element method for 3D Fourier-Kirchhoff heat transfer equation 7 z:=z[]-hs/; z:=z[]-hs/; z3:=z[]+hs/; z4:=z[]+hs/; end else if <= *n*n3+n*n3 then {tyl} x:=x[]+hs/; x:=x[]-hs/; x3:=x[]-hs/; x4:=x[]+hs/; y:=l; y:=l; y3:=l; y4:=l; z:=z[]-hs/; z:=z[]-hs/; z3:=z[]+hs/; z4:=z[]+hs/; end else if <= *n*n3+*n*n3 then {lewa} x:=0.0; x:=0.0; x3:=0.0; x4:=0.0; y:=y[]+hs/; y:=y[]+hs/; y3:=y[]-hs/; y4:=y[]-hs/; z:=z[]+hs/; z:=z[]-hs/; z3:=z[]-hs/; z4:=z[]+hs/; end else if <= *n*n3+*n*n3+n*n then {dol} x:=x[]-hs/; x:=x[]-hs/; x3:=x[]+hs/; x4:=x[]+hs/; y:=y[]-hs/; y:=y[]+hs/; y3:=y[]+hs/; y4:=y[]-hs/; z:=0.0; z:=0.0; z3:=0.0; z4:=0.0; end else if <= N then {gora}\ x:=x[]-hs/; x:=x[]+hs/; x3:=x[]+hs/; x4:=x[]-hs/; y:=y[]-hs/; y:=y[]-hs/; y3:=y[]+hs/; y4:=y[]+hs/; z:=l3; z:=l3; z3:=l3; z4:=l3; After determining the integrals G i,, H i,, U i,, P i,l and taking into account the boundary conditions (3), the system of equations (9) can be solved by means of the Gaussian elimination method []. 5. Example of computations The following input data are introduced: thermal conductivity λ = 0 W/(mK), volumetric specific heat c = 0 6 J/(m 3 K), velocity u = m/s, source function Q = 0 W/m 3. The cuboid of dimension 0.05 x 0.05 x 0.05 m 3 is considered as shown in Figure. It s assumed that n = n = 0, n 3 = 5, so N = 400 boundary elements have been distinguished. It should be pointed out that in the case Q = 0 only the boundary should be discretized (c.f. equations (9) and (8)).
8 8 E. Machrzak, Ł. Turchan x 3 x l 3 l l Fig. Domain considered x The following boundary conditions are accepted: x = 0, 0 < x < l, 0 < x < l : T x = 50 x = l, 0 < x < l, 0 < x < l : T x = 00 0 < x < l, x = 0, 0 < x < l : q x = 0 0 < x < l, x = l, 0 < x < l : q x = 0 0 < x < l, 0 < x < l, x = 0 : q x = 0 0 < x < l, 0 < x < l, x = l : q x = 0 (5) In Figure the temperature distribution in the plane x = l / is shown. Taking into account the assumed boundary conditions, the results are the same for each plane x = s, where s [0, l ]. Fig.. Temperature distribution (plane x = l / ) The solution obtained can be compared with analytical solution for D problem uc Ta ( x ) = C + C exp x λ (6)
9 Boundary element method for 3D Fourier-Kirchhoff heat transfer equation 9 where C, C are the integral constants determined from boundary conditions T a (0) = 50 C and T a (l ) = 00 C. T [ C] Fig. 3. Comparison of numerical and analytical solutions x [m] The results are shown in Figure 3, where the solid line presents the analytical solution, while the symbols illustrate the numerical one. It is visible, that both solutions are practically the same. Conclusions Application of the boundary element method for numerical solution of 3D Fourier-Kirchhoff equation is presented. Both the theoretical and practical aspects of the problem solution are discussed. Among others, the way of discretization for cuboid using constant boundary elements is shown. Numerical results compared with analytical solution confirm the exactness and effectiveness of the algorithm proposed. The method presented can be applied for numerical modelling of heat transfer proceeding in domain of porous media, in particular the bioheat transfer equation basing on the theory of porous media [4] can be considered. Acknowledgement This paper is a part of Grant No N N
10 30 E. Machrzak, Ł. Turchan References [] Brebbia C.A., Dominguez J., Boundary elements, an introductory course, Computational Mechanics Publications, McGraw-Hill Book Company, London 99. [] Machrzak E., Metoda elementów brzegowych w przepływie ciepła, Wyd. Pol. Częstochowskie, Częstochowa 00. [3] Klekot J., Boundary element method for elliptic equation - determination of fundamental solution. Scientific Research of the Institute of Mathematics and Computer Science, Czestochowa University of Technology 00, (9), 6-70 [4] Nakayama A., Kuwanara F., A general bioheat transfer model based on the theory of porous media, Int. J. of Heat and Mass Transfer 008, 5,
DETERMINATION OF THE TEMPERATURE FIELD IN BURNED AND HEALTHY SKIN TISSUE USING THE BOUNDARY ELEMENT METHOD - PART I
Journal of Applied Mathematics and Computational Mechanics 03, (3), 39-46 DETERMINATION OF THE TEMPERATURE FIELD IN BURNED AND HEALTHY SKIN TISSUE USING THE BOUNDARY ELEMENT METHOD - PART I Katarzyna Freus,
More informationSensitivity of the skin tissue on the activity of external heat sources
Copyright c 23 Tech Science Press CMES, vol.4, no.3&4, pp.431-438, 23 Sensitivity of the skin tissue on the activity of external heat sources B. Mochnacki 1 E. Majchrzak 2 Abstract: In the paper the analysis
More informationTwo Posts to Fill On School Board
Y Y 9 86 4 4 qz 86 x : ( ) z 7 854 Y x 4 z z x x 4 87 88 Y 5 x q x 8 Y 8 x x : 6 ; : 5 x ; 4 ( z ; ( ) ) x ; z 94 ; x 3 3 3 5 94 ; ; ; ; 3 x : 5 89 q ; ; x ; x ; ; x : ; ; ; ; ; ; 87 47% : () : / : 83
More information' Liberty and Umou Ono and Inseparablo "
3 5? #< q 8 2 / / ) 9 ) 2 ) > < _ / ] > ) 2 ) ) 5 > x > [ < > < ) > _ ] ]? <
More informationLOWELL WEEKLY JOURNAL
Y -» $ 5 Y 7 Y Y -Y- Q x Q» 75»»/ q } # ]»\ - - $ { Q» / X x»»- 3 q $ 9 ) Y q - 5 5 3 3 3 7 Q q - - Q _»»/Q Y - 9 - - - )- [ X 7» -» - )»? / /? Q Y»» # X Q» - -?» Q ) Q \ Q - - - 3? 7» -? #»»» 7 - / Q
More informationScientific Research of the Institute of Mathematics and Computer Science DIFFERENT VARIANTS OF THE BOUNDARY ELEMENT METHOD FOR PARABOLIC EQUATIONS
Scieniic Research o he Insiue o Mahemaics and Compuer Science DIERENT VARIANTS O THE BOUNDARY ELEMENT METHOD OR PARABOLIC EQUATIONS Ewa Majchrzak,, Ewa Ładyga Jerzy Mendakiewicz, Alicja Piasecka Belkhaya
More informationEducjatipnal. L a d ie s * COBNWALILI.S H IG H SCHOOL. I F O R G IR L S A n B k i n d e r g a r t e n.
- - - 0 x ] - ) ) -? - Q - - z 0 x 8 - #? ) 80 0 0 Q ) - 8-8 - ) x ) - ) -] ) Q x?- x - - / - - x - - - x / /- Q ] 8 Q x / / - 0-0 0 x 8 ] ) / - - /- - / /? x ) x x Q ) 8 x q q q )- 8-0 0? - Q - - x?-
More informationM E M P H I S, T E N N., S A T U E D A Y, OCTOBER 8, 1870.
5 L V 8 5 x - L : L Q ) L - \ \ Q Q - V 84 z < L L 4 Y z ( (
More informationNumerical Solution of Dual Phase Lag Model of Bioheat Transfer Using the General Boundary Element Method
Copyright 2010 Tech Science Press CMES, vol.69, no.1, pp.43-60, 2010 Numerical Solution of Dual Phase Lag Model of Bioheat Transfer Using the General Boundary Element Method Ewa Majchrzak 1 Abstract: Heat
More informationP A L A C E P IE R, S T. L E O N A R D S. R a n n o w, q u a r r y. W WALTER CR O TC H, Esq., Local Chairman. E. CO O PER EVANS, Esq.,.
? ( # [ ( 8? [ > 3 Q [ ««> » 9 Q { «33 Q> 8 \ \ 3 3 3> Q»«9 Q ««« 3 8 3 8 X \ [ 3 ( ( Z ( Z 3( 9 9 > < < > >? 8 98 ««3 ( 98 < # # Q 3 98? 98 > > 3 8 9 9 ««««> 3 «>
More informationA DARK GREY P O N T, with a Switch Tail, and a small Star on the Forehead. Any
Y Y Y X X «/ YY Y Y ««Y x ) & \ & & } # Y \#$& / Y Y X» \\ / X X X x & Y Y X «q «z \x» = q Y # % \ & [ & Z \ & { + % ) / / «q zy» / & / / / & x x X / % % ) Y x X Y $ Z % Y Y x x } / % «] «] # z» & Y X»
More informationOWELL WEEKLY JOURNAL
Y \»< - } Y Y Y & #»»» q ] q»»»>) & - - - } ) x ( - { Y» & ( x - (» & )< - Y X - & Q Q» 3 - x Q Y 6 \Y > Y Y X 3 3-9 33 x - - / - -»- --
More informationINFLUENCE OF EMISSIVITY CHANGES ON THE BLOOD FLOW RATE DETERMINED ON THE BASIS OF HEAT BALANCE EQUATION
Scientific Research of the Institute of Mathematics and Computer Science INFLUENCE OF EMISSIVITY CHANGES ON THE BLOOD FLOW RATE DETERMINED ON THE BASIS OF HEAT BALANCE EQUATION Marek Jasiński Department
More information9 th International Conference on Quantitative InfraRed Thermography July 2-5, 2008, Krakow - Poland Application of infrared thermography for validation of numerical analyses results of a finned cross-flow
More informationChapter 2: Heat Conduction Equation
-1 General Relation for Fourier s Law of Heat Conduction - Heat Conduction Equation -3 Boundary Conditions and Initial Conditions -1 General Relation for Fourier s Law of Heat Conduction (1) The rate of
More informationLinear System. Lotfi A. Zadeh & Charles A. Desoer. The State Space Approach
Linear System The State Space Approach Lotfi A. Zadeh & Charles A. Desoer Department of Electrical Engineering University of California Berkeley, California McGraw-Hill Book Company New York / San Francisco
More informationr/lt.i Ml s." ifcr ' W ATI II. The fnncrnl.icniccs of Mr*. John We mil uppn our tcpiiblicnn rcprc Died.
$ / / - (\ \ - ) # -/ ( - ( [ & - - - - \ - - ( - - - - & - ( ( / - ( \) Q & - - { Q ( - & - ( & q \ ( - ) Q - - # & - - - & - - - $ - 6 - & # - - - & -- - - - & 9 & q - / \ / - - - -)- - ( - - 9 - - -
More informationW i n t e r r e m e m b e r t h e W O O L L E N S. W rite to the M anageress RIDGE LAUNDRY, ST. H E LE N S. A uction Sale.
> 7? 8 «> ««0? [ -! ««! > - ««>« ------------ - 7 7 7 = - Q9 8 7 ) [ } Q ««
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter Boundary-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or
More informationFast multipole boundary element method for the analysis of plates with many holes
Arch. Mech., 59, 4 5, pp. 385 401, Warszawa 2007 Fast multipole boundary element method for the analysis of plates with many holes J. PTASZNY, P. FEDELIŃSKI Department of Strength of Materials and Computational
More informationA. H. Hall, 33, 35 &37, Lendoi
7 X x > - z Z - ----»»x - % x x» [> Q - ) < % - - 7»- -Q 9 Q # 5 - z -> Q x > z»- ~» - x " < z Q q»» > X»? Q ~ - - % % < - < - - 7 - x -X - -- 6 97 9
More informationLOWHLL #WEEKLY JOURNAL.
# F 7 F --) 2 9 Q - Q - - F - x $ 2 F? F \ F q - x q - - - - )< - -? - F - - Q z 2 Q - x -- - - - 3 - % 3 3 - - ) F x - \ - - - - - q - q - - - - -z- < F 7-7- - Q F 2 F - F \x -? - - - - - z - x z F -
More informationCeiling mounted radiant panels calculations of heat output in heating and cooling application
Ceiling mounted radiant panels calculations of heat output in heating and cooling application Lawrence Drojetzki 1,*, and Janusz Wojtkowiak 1 1 Poznan University of Technology, Institute of Environmental
More informationTime-Dependent Conduction :
Time-Dependent Conduction : The Lumped Capacitance Method Chapter Five Sections 5.1 thru 5.3 Transient Conduction A heat transfer process for which the temperature varies with time, as well as location
More informationMANY BILLS OF CONCERN TO PUBLIC
- 6 8 9-6 8 9 6 9 XXX 4 > -? - 8 9 x 4 z ) - -! x - x - - X - - - - - x 00 - - - - - x z - - - x x - x - - - - - ) x - - - - - - 0 > - 000-90 - - 4 0 x 00 - -? z 8 & x - - 8? > 9 - - - - 64 49 9 x - -
More informationPanHomc'r I'rui;* :".>r '.a'' W"»' I'fltolt. 'j'l :. r... Jnfii<on. Kslaiaaac. <.T i.. %.. 1 >
5 28 (x / &» )»(»»» Q ( 3 Q» (» ( (3 5» ( q 2 5 q 2 5 5 8) 5 2 2 ) ~ ( / x {» /»»»»» (»»» ( 3 ) / & Q ) X ] Q & X X X x» 8 ( &» 2 & % X ) 8 x & X ( #»»q 3 ( ) & X 3 / Q X»»» %» ( z 22 (»» 2» }» / & 2 X
More informationCHAPTER 8: Thermal Analysis
CHAPER 8: hermal Analysis hermal Analysis: calculation of temperatures in a solid body. Magnitude and direction of heat flow can also be calculated from temperature gradients in the body. Modes of heat
More informationMEEN 673 ********************************* NONLINEAR FINITE ELEMENT ANALYSIS ********************************* J. N. Reddy
MEEN 673 ********************************* NONLINEAR FINITE ELEMENT ANALYSIS ********************************* J. N. Reddy e-mail: jnreddy@tamu.edu Texas A&M University College Station, TX 77843-3123 USA
More informationME5286 Robotics Spring 2017 Quiz 2
Page 1 of 5 ME5286 Robotics Spring 2017 Quiz 2 Total Points: 30 You are responsible for following these instructions. Please take a minute and read them completely. 1. Put your name on this page, any other
More informationModule 7: The Laplace Equation
Module 7: The Laplace Equation In this module, we shall study one of the most important partial differential equations in physics known as the Laplace equation 2 u = 0 in Ω R n, (1) where 2 u := n i=1
More informationConjugate gradient-boundary element method for a Cauchy problem in the Lam6 system
Conjugate gradient-boundary element method for a Cauchy problem in the Lam6 system L. Marin, D.N. Hho & D. Lesnic I Department of Applied Mathematics, University of Leeds, UK 'Hanoi Institute of Mathematics,
More informationIntroduction to Heat and Mass Transfer. Week 7
Introduction to Heat and Mass Transfer Week 7 Example Solution Technique Using either finite difference method or finite volume method, we end up with a set of simultaneous algebraic equations in terms
More informationLecture Note 3: Polynomial Interpolation. Xiaoqun Zhang Shanghai Jiao Tong University
Lecture Note 3: Polynomial Interpolation Xiaoqun Zhang Shanghai Jiao Tong University Last updated: October 24, 2013 1.1 Introduction We first look at some examples. Lookup table for f(x) = 2 π x 0 e x2
More informationLOWELL WEEKLY JOURNAL
G $ G 2 G ««2 ««q ) q «\ { q «««/ 6 «««««q «] «q 6 ««Z q «««Q \ Q «q «X ««G X G ««? G Q / Q Q X ««/«X X «««Q X\ «q «X \ / X G XX «««X «x «X «x X G X 29 2 ««Q G G «) 22 G XXX GG G G G G G X «x G Q «) «G
More informationSENSITIVITY OR TEMPERATURE FIELD IN THE SYSTEM PROTECTIVE CLOTHING FOREARM WITH RESPECT TO PERTURBATIONS OF EXTERNAL HEATING CONDITIONS
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
More informationA Green s function for steady-state two-dimensional isotropic heat conduction across a homogeneously imperfect interface
A Green s function for steady-state two-dimensional isotropic heat conduction across a homogeneously imperfect interface W. T. Ang*, K. K. Choo and H. Fan Division of Engineering Mechanics, School of Mechanical
More informationAnalytical Modeling of Laser Moving Sources
Analytical Modeling of Laser Moving Sources Contains: Heat flow equation Analytic model in one dimensional heat flow Heat source modeling Point heat source Line heat source Plane heat source Surface heat
More informationCOMPUTATIONAL COMPLEXITY OF PARAMETRIC LINEAR PROGRAMMING +
Mathematical Programming 19 (1980) 213-219. North-Holland Publishing Company COMPUTATIONAL COMPLEXITY OF PARAMETRIC LINEAR PROGRAMMING + Katta G. MURTY The University of Michigan, Ann Arbor, MI, U.S.A.
More informationPARAMETRIC INTEGRAL EQUATIONS SYSTEMS METHOD IN SOLVING UNSTEADY HEAT TRANSFER PROBLEMS FOR LASER HEATED MATERIALS
DOI 1.1515/ama-15-8 acta mechanica et automatica, vol.9 no.3 (15) PARAMETRIC INTEGRAL EQUATIONS SYSTEMS METHOD IN SOLVING UNSTEADY HEAT TRANSFER PROBLEMS FOR LASER HEATED MATERIALS Dominik SAWICKI *, Eugeniusz
More informationTransactions on Modelling and Simulation vol 12, 1996 WIT Press, ISSN X
Simplifying integration for logarithmic singularities R.N.L. Smith Department ofapplied Mathematics & OR, Cranfield University, RMCS, Shrivenham, Swindon, Wiltshire SN6 SLA, UK Introduction Any implementation
More informationModeling using conservation laws. Let u(x, t) = density (heat, momentum, probability,...) so that. u dx = amount in region R Ω. R
Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...) so that u dx = amount in region R Ω. R Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...)
More informationDepartment of mathematics MA201 Mathematics III
Department of mathematics MA201 Mathematics III Academic Year 2015-2016 Model Solutions: Quiz-II (Set - B) 1. Obtain the bilinear transformation which maps the points z 0, 1, onto the points w i, 1, i
More informationDeterministic solution of stochastic groundwater flow equations by nonlocalfiniteelements A. Guadagnini^ & S.P. Neuman^
Deterministic solution of stochastic groundwater flow equations by nonlocalfiniteelements A. Guadagnini^ & S.P. Neuman^ Milano, Italy; ^Department of. Hydrology & Water Resources, The University ofarizona,
More informationYour use of the JSTOR archive indicates your acceptance of the Terms & Conditions of Use, available at
Regression Analysis when there is Prior Information about Supplementary Variables Author(s): D. R. Cox Source: Journal of the Royal Statistical Society. Series B (Methodological), Vol. 22, No. 1 (1960),
More informationTHE FINITE ELEMENT METHOD 2017 Dept. of Solid Mechanics
THE FINITE ELEMENT METHOD 07 Dept. of Solid Mechanics EXAMINATION: 07-05-9 A maximum of 60 points can be achieved in this examination. To pass at least 0 points are required. Permitted aid: Pocket calculator.
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationMath 216 First Midterm 19 October, 2017
Math 6 First Midterm 9 October, 7 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material that
More informationFUZZY BOUNDARY ELEMENT METHODS: A NEW MULTI-SCALE PERTURBATION APPROACH FOR SYSTEMS WITH FUZZY PARAMETERS
MODELOWANIE INŻYNIERSKIE ISNN 896-77X 3, s. 433-438, Gliwice 6 FUZZY BOUNDARY ELEMENT METHODS: A NEW MULTI-SCALE PERTURBATION APPROACH FOR SYSTEMS WITH FUZZY PARAMETERS JERZY SKRZYPCZYK HALINA WITEK Zakład
More informationME 1401 FINITE ELEMENT ANALYSIS UNIT I PART -A. 2. Why polynomial type of interpolation functions is mostly used in FEM?
SHRI ANGALAMMAN COLLEGE OF ENGINEERING AND TECHNOLOGY (An ISO 9001:2008 Certified Institution) SIRUGANOOR, TIRUCHIRAPPALLI 621 105 Department of Mechanical Engineering ME 1401 FINITE ELEMENT ANALYSIS 1.
More informationElectro-Thermal-Mechanical MEMS Devices Multi-Physics Problems
Electro-Thermal-Mechanical MEMS Devices Multi-Physics Problems Presentation for ASEN 5519 Joseph Pajot Department of Aerospace Engineering CU Boulder November 15, 2004 Typeset by FoilTEX 1 The ETM Problems
More informationThis examination consists of 10 pages. Please check that you have a complete copy. Time: 2.5 hrs INSTRUCTIONS
THE UNIVERSITY OF BRITISH COLUMBIA Department of Electrical and Computer Engineering EECE 564 Detection and Estimation of Signals in Noise Final Examination 08 December 2009 This examination consists of
More informationENTROPY GENERATION OF CONVECTION HEAT TRANSFER IN AN ASYMMETRICALLY HEATED PACKED DUCT
University of Nebraska - Lincoln From the SelectedWorks of YASAR DEMIREL 1997 ENTROPY GENERATION OF CONVECTION HEAT TRANSFER IN AN ASYMMETRICALLY HEATED PACKED DUCT YASAR DEMIREL H.H. Ali B.A. Abu-Al-Saud
More informationFourier Variant Homogenization of the Heat Transfer Processes in Periodic Composites
Mechanics and Mechanical Engineering Vol. 22, No. 3 (2018) 719 726 c Lodz University of Technology Fourier Variant Homogenization of the Heat Transfer Processes in Periodic Composites Ewaryst Wierzbicki
More informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More informationME FINITE ELEMENT ANALYSIS FORMULAS
ME 2353 - FINITE ELEMENT ANALYSIS FORMULAS UNIT I FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PROBLEMS 01. Global Equation for Force Vector, {F} = [K] {u} {F} = Global Force Vector [K] = Global Stiffness
More informationTransactions on Modelling and Simulation vol 8, 1994 WIT Press, ISSN X
Model analysis of plates using the dual reciprocity boundary element method T.W. Davies & F.A. Moslehy Department of Mechanical and Aerospace Engineering, University of Central Florida, Orlando, Florida,
More informationCHAPTER II MATHEMATICAL BACKGROUND OF THE BOUNDARY ELEMENT METHOD
CHAPTER II MATHEMATICAL BACKGROUND OF THE BOUNDARY ELEMENT METHOD For the second chapter in the thesis, we start with surveying the mathematical background which is used directly in the Boundary Element
More informationAnalytical solutions of heat transfer for laminar flow in rectangular channels
archives of thermodynamics Vol. 35(2014), No. 4, 29 42 DOI: 10.2478/aoter-2014-0031 Analytical solutions of heat transfer for laminar flow in rectangular channels WITOLD RYBIŃSKI 1 JAROSŁAW MIKIELEWICZ
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 informationELECTROMAGNETIC FIELDS AND RELATIVISTIC PARTICLES
ELECTROMAGNETIC FIELDS AND RELATIVISTIC PARTICLES Emil J. Konopinski Professor of Physics Indiana University McGraw-Hill Book Company New York St. Louis San Francisco Auckland Bogota Hamburg Johannesburg
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationThe Finite Difference Method
Chapter 5. The Finite Difference Method This chapter derives the finite difference equations that are used in the conduction analyses in the next chapter and the techniques that are used to overcome computational
More informationOne dimensional steady state diffusion, with and without source. Effective transfer coefficients
One dimensional steady state diffusion, with and without source. Effective transfer coefficients 2 mars 207 For steady state situations t = 0) and if convection is not present or negligible the transport
More informationChapter 1 INTRODUCTION AND BASIC CONCEPTS
Heat and Mass Transfer: Fundamentals & Applications 5th Edition in SI Units Yunus A. Çengel, Afshin J. Ghajar McGraw-Hill, 2015 Chapter 1 INTRODUCTION AND BASIC CONCEPTS Mehmet Kanoglu University of Gaziantep
More informationCode No: RT41033 R13 Set No. 1 IV B.Tech I Semester Regular Examinations, November - 2016 FINITE ELEMENT METHODS (Common to Mechanical Engineering, Aeronautical Engineering and Automobile Engineering)
More informationFeature detection.
Feature detection Kim Steenstrup Pedersen kimstp@itu.dk The IT University of Copenhagen Feature detection, The IT University of Copenhagen p.1/20 What is a feature? Features can be thought of as symbolic
More informationIntroduction to Heat and Mass Transfer. Week 8
Introduction to Heat and Mass Transfer Week 8 Next Topic Transient Conduction» Analytical Method Plane Wall Radial Systems Semi-infinite Solid Multidimensional Effects Analytical Method Lumped system analysis
More informationCharacteristic finite-difference solution Stability of C C (CDS in time/space, explicit): Example: Effective numerical wave numbers and dispersion
Spring 015 Lecture 14 REVIEW Lecture 13: Stability: Von Neumann Ex.: 1st order linear convection/wave eqn., F-B scheme Hyperbolic PDEs and Stability nd order wave equation and waves on a string Characteristic
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 informationSalmon: Lectures on partial differential equations
6. The wave equation Of the 3 basic equations derived in the previous section, we have already discussed the heat equation, (1) θ t = κθ xx + Q( x,t). In this section we discuss the wave equation, () θ
More informationThermal Analysis. with SolidWorks Simulation 2013 SDC. Paul M. Kurowski. Better Textbooks. Lower Prices.
Thermal Analysis with SolidWorks Simulation 2013 Paul M. Kurowski SDC PUBLICATIONS Schroff Development Corporation Better Textbooks. Lower Prices. www.sdcpublications.com Visit the following websites to
More informationChapter 11. Taylor Series. Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27
Chapter 11 Taylor Series Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27 First-Order Approximation We want to approximate function f by some simple function. Best possible approximation
More informationIntroduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA
Introduction to Finite and Spectral Element Methods Using MATLAB Second Edition C. Pozrikidis University of Massachusetts Amherst, USA (g) CRC Press Taylor & Francis Group Boca Raton London New York CRC
More informationQ SON,' (ESTABLISHED 1879L
( < 5(? Q 5 9 7 00 9 0 < 6 z 97 ( # ) $ x 6 < ( ) ( ( 6( ( ) ( $ z 0 z z 0 ) { ( % 69% ( ) x 7 97 z ) 7 ) ( ) 6 0 0 97 )( 0 x 7 97 5 6 ( ) 0 6 ) 5 ) 0 ) 9%5 z» 0 97 «6 6» 96? 0 96 5 0 ( ) ( ) 0 x 6 0
More informationUniversity of Hertfordshire Department of Mathematics. Study on the Dual Reciprocity Boundary Element Method
University of Hertfordshire Department of Mathematics Study on the Dual Reciprocity Boundary Element Method Wattana Toutip Technical Report 3 July 999 Preface The boundary Element method (BEM) is now recognised
More informationThe minus sign indicates that the centroid is located below point E. We will relocate the axis as shown in Figure (1) and take discard the sign:
AOE 304: Thin Walled Structures Solutions to Consider a cantilever beam as shown in the attached figure. At the tip of the beam, a bending moment M = 1000 N-m is applied at an angle θ with respect to the
More informationLecture Note 3: Interpolation and Polynomial Approximation. Xiaoqun Zhang Shanghai Jiao Tong University
Lecture Note 3: Interpolation and Polynomial Approximation Xiaoqun Zhang Shanghai Jiao Tong University Last updated: October 10, 2015 2 Contents 1.1 Introduction................................ 3 1.1.1
More informationMULTI-SCALE PERTURBATION METHODS IN MECHANICS
006/3 4 PAGES 13 17 RECEIVED 18 9 006 ACCEPTED 7 11 006 J SKRZYPCZYK MULTI-SCALE PERTURBATION METHODS IN MECHANICS Dr Hab Eng J Skrzypczyk, Prof ul Krzywoustego 7, 44-100 Gliwice, Pol, Phone:+48 604 540
More informationMETHOD OF LINES FOR THE INCOMPRESSIBLE NAVIER-STOKES EQUATIONS IN THE STREAM-FUNCTION FORMULATION
TASKQUARTERLY9No1,17 35 METHOD OF LINES FOR THE INCOMPRESSIBLE NAVIER-STOKES EQUATIONS IN THE STREAM-FUNCTION FORMULATION ZBIGNIEW KOSMA Institute of Applied Mechanics, Radom University of Technology,
More informationUniversity of Illinois at Urbana-Champaign College of Engineering
University of Illinois at Urbana-Champaign College of Engineering CEE 570 Finite Element Methods (in Solid and Structural Mechanics) Spring Semester 2014 Quiz #2 April 14, 2014 Name: SOLUTION ID#: PS1.:
More informationIntroduction to Mathematical Modelling: Ordinary Differential Equations and Heat Equations
Introduction to Mathematical Modelling: Ordinary Differential Equations and Heat Equations Knut Andreas Lie SINTEF ICT, Dept Applied Mathematics INF2340 / Spring 2005 p Overview 1 Ordinary differential
More informationSimulation of the Temperature Profile During Welding with COMSOL Multiphysics Software Using Rosenthal s Approach
Simulation of the Temperature Profile During Welding with COMSOL Multiphysics Software Using Rosenthal s Approach A. Lecoanet*, D. G. Ivey, H. Henein Department of Chemical & Materials Engineering, University
More informationChapter 2 HEAT CONDUCTION EQUATION
Heat and Mass Transfer: Fundamentals & Applications 5th Edition in SI Units Yunus A. Çengel, Afshin J. Ghajar McGraw-Hill, 2015 Chapter 2 HEAT CONDUCTION EQUATION Mehmet Kanoglu University of Gaziantep
More informationMethod of Weighted Residual Approach for Thermal Performance Evaluation of a Circular Pin fin: Galerkin Finite element Approach
ISSN (Print): 2347-671 (An ISO 3297: 27 Certified Organization) Vol. 6, Issue 4, April 217 Method of Weighted Residual Approach for Thermal Performance Evaluation of a Circular Pin fin: Galerkin Finite
More informationChapter 2 HEAT CONDUCTION EQUATION
Heat and Mass Transfer: Fundamentals & Applications Fourth Edition Yunus A. Cengel, Afshin J. Ghajar McGraw-Hill, 2011 Chapter 2 HEAT CONDUCTION EQUATION Mehmet Kanoglu University of Gaziantep Copyright
More informationX/92/ $ by John Wiley & Sons, Ltd. INTERNATIONAL JOURNAL OF ENERGY RESEARCH, VOL. 16, (1992)
INTERNATIONAL JOURNAL OF ENERGY RESEARCH, VOL. 16,401-41 1 (1992) EXACT SOLUTION OF COUPLED HEAT AND MASS TRANSFER WITH DOUBLE MOVING INTERFACES IN A POROUS HALF-SPACE PART 2: MASS TRANSFER CONTROLLED
More informationThermal Analysis with SOLIDWORKS Simulation 2015 and Flow Simulation 2015
Thermal Analysis with SOLIDWORKS Simulation 2015 and Flow Simulation 2015 Paul M. Kurowski SDC PUBLICATIONS Better Textbooks. Lower Prices. www.sdcpublications.com Powered by TCPDF (www.tcpdf.org) Visit
More informationMA261-A Calculus III 2006 Fall Midterm 2 Solutions 11/8/2006 8:00AM ~9:15AM
MA6-A Calculus III 6 Fall Midterm Solutions /8/6 8:AM ~9:5AM. Find the it xy cos y (x;y)(;) 3x + y, if it exists, or show that the it does not exist. Assume that x. The it becomes (;y)(;) y cos y 3 + y
More informationUniversity of Illinois at Urbana-Champaign College of Engineering
University of Illinois at Urbana-Champaign College of Engineering CEE 570 Finite Element Methods (in Solid and Structural Mechanics) Spring Semester 03 Quiz # April 8, 03 Name: SOUTION ID#: PS.: A the
More informationOn the Three Dimensional Laplace Problem with Dirichlet Condition
Applied Mathematical Sciences, Vol. 8, 2014, no. 83, 4097-4101 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.45356 On the Three Dimensional Laplace Problem with Dirichlet Condition P.
More informationME 315 Final Examination Solution 8:00-10:00 AM Friday, May 8, 2009 CIRCLE YOUR DIVISION
ME 315 Final Examination Solution 8:00-10:00 AM Friday, May 8, 009 This is a closed-book, closed-notes examination. There is a formula sheet at the back. You must turn off all communications devices before
More informationqxbxg. That is, the heat rate within the object is everywhere constant. From Fourier s
PROBLEM.1 KNOWN: Steady-state, one-dimensional heat conduction through an axisymmetric shape. FIND: Sketch temperature distribution and explain shape of curve. ASSUMPTIONS: (1) Steady-state, one-dimensional
More informationLOWELL WEEKLY JOURNAL.
Y $ Y Y 7 27 Y 2» x 7»» 2» q» ~ [ } q q $ $ 6 2 2 2 2 2 2 7 q > Y» Y >» / Y» ) Y» < Y»» _»» < Y > Y Y < )»» >» > ) >» >> >Y x x )»» > Y Y >>»» }> ) Y < >» /» Y x» > / x /»»»»» >» >» >»» > > >» < Y /~ >
More informationACCEPTS HUGE FLORAL KEY TO LOWELL. Mrs, Walter Laid to Rest Yesterday
$ j < < < > XXX Y 928 23 Y Y 4% Y 6 -- Q 5 9 2 5 Z 48 25 )»-- [ Y Y Y & 4 j q - Y & Y 7 - -- - j \ -2 -- j j -2 - - - - [ - - / - ) ) - - / j Y 72 - ) 85 88 - / X - j ) \ 7 9 Y Y 2 3» - ««> Y 2 5 35 Y
More informationInternational Conference on Mechanical, Industrial and Energy Engineering December, 2014, Khulna, BANGLADESH
International Conference on Mechanical, Industrial and Energy Engineering 214 25-2 December, 214, Khulna, BANGLADESH ICMIEE-PI-13581 Analysis of Bio-heat Transfer Problem Using Finite Element Approach
More informationSECM Study of Permeability of a Thiolated Aryl Multilayer and Imaging of Single Nanocubes Anchored to It
Supporting Information SECM Study of Permeability of a Thiolated Aryl Multilayer and Imaging of Single Nanocubes Anchored to It Pierre-Yves Blanchard, Tong Sun, Yun Yu, Zengyan Wei, Hiroshi Matsui, # *
More informationPROBLEM 1.2 ( ) 25 C 15 C dx L 0.30 m Ambient air temperature, T2 (C)
PROBLEM 1.2 KNOWN: Inner surface temperature and thermal conductivity of a concrete wall. FIND: Heat loss by conduction through the wall as a function of ambient air temperatures ranging from -15 to 38
More information