Nonlinear transient heat conduction problems for a class of inhomogeneous anisotropic materials by BEM
|
|
- Helen Boyd
- 5 years ago
- Views:
Transcription
1 Engineering Analysis with Boundary Elements 32 (28) Nonlinear transient heat conduction problems for a class of inhomogeneous anisotropic materials by BEM M.I. Azis a, D.L. Clements b, a Jurusan Matematika, Fakultas MIPA, Universitas Hasanuddin, Jl. P. Kemerdekaan KM.1 Tamalanrea, Makassar 9245, Indonesia b School of Mathematics, University of Adelaide, SA 55, Australia Received 2 March 27; accepted 26 April 27 Available online 2 February 28 Abstract A class of nonlinear transient heat conduction problems for anisotropic inhomogeneous materials is considered. The problems under consideration are reduced to a boundary integral euation which may be solved numerically using standard techniues. Some numerical examples are considered in order to test the accuracy of the numerical procedure. r 28 Elsevier Ltd. All rights reserved. Keywords: Boundary integral methods; Heat transfer; Nonlinear differential euations; Numerical methods 1. Introduction Boundary integral euation methods have been used by a number of authors to consider steady-state and transient heat conduction problems for inhomogeneous materials. For example Divo and Kassab in [1] considered the linear elliptic euation governing the steady-state temperature field in inhomogeneous isotropic media while Clements et al. [2] obtained a boundary integral euation procedure which may be used to determine the steady-state temperature field for a wide class of inhomogeneous anisotropic materials. For time dependent problems Abdullah et al. [3] and Sutradhar et al. [4] employed the Laplace transform and a boundary integral formulation to solve problems for particular classes of isotropic inhomogeneous materials. Nonlinear heat conduction problems have been investigated by, for example, Clements and Budhi [5] who considered steady-state problems for a class of materials for which the thermal conductivity varied with the temperature and the three spatial coordinates. In the case of transient nonlinear problems boundary integral euation methods have been employed by a number of authors Corresponding author. addresses: mazis@indosat.net.id (M.I. Azis), david.clements@adelaide.edu.au (D.L. Clements). including Bialecki and Kuhn [6], oto and Suzuki [7] and Kikuta et al. [8]. This paper is concerned with nonlinear transient heat conduction problems for a class of anisotropic inhomogeneous materials. An outline of the paper is as follows. The Kirchhoff transformation is employed to convert the governing nonlinear euation to a linear second order partial differential euation with variable coefficients. A further transformation is then applied in order to transform the variable coefficients euation to a linear euation with constant coefficients. A boundary integral euation formulation is obtained for the solution of the constant coefficients euation. Then use of the inverse transformation provides a boundary integral euation for the solution of the linear variable coefficients euation. nce numerical values have been obtained from this boundary integral euation use of the inverse of the Kirchhoff transformation provides the temperature field for particular boundary value problems. Some numerical examples are considered in order to assess the effectiveness of the numerical procedure. Essentially the present work may be considered to be a generalisation of the paper of Clements and Budhi [5] to include transient effects and, more particularly, a generalisation of the paper by oto and Suzuki [7] to include a restricted class of thermal coefficients which vary with position /$ - see front matter r 28 Elsevier Ltd. All rights reserved. doi:1.116/j.enganabound
2 M.I. Azis, D.L. Clements / Engineering Analysis with Boundary Elements 32 (28) Basic euations Referred to the Cartesian frame x 1 x 2 x 3 this paper is concerned with obtaining the solution of the heat conduction problems governed by the euation Tðx; tþ Tðx; tþ k ij ðx; TÞ ¼ rðx; TÞcðx; TÞ x i x j for i; j ¼ 1; 2; 3, (1) T is the temperature, k ij are the heat conductivities, r is the density, c denotes the heat capacity, x ¼ðx 1 ; x 2 ; x 3 Þ, t is the time variable and the repeated suffix summation convention is employed. The analysis employed by oto and Suzuki in [7] is used and extended to derive a boundary element method (BEM) for the solution of (1). The matrix of the heat conductivities k ij is assumed to be a real symmetric positive definite matrix. Also the heat conductivities k ij are assumed to take the form k ij ðx; TÞ ¼m ij ðxþhðtþ, (2) m ij are differentiable functions and h is an integrable function. The analysis is specially relevant to an anisotropic medium but it eually applies to the isotropic case. For isotropic media the coefficients k ij in (1) take the form k 11 ¼ k 22 ¼ k 33 and k 12 ¼ k 13 ¼ k 23 ¼ and use of these euations in the following analysis immediately yields the corresponding results for isotropic materials. 3. The initial-boundary value problem A solution to (1) is sought in a domain for time tx. The boundary of consists of a finite number of piecewise smooth closed surfaces. In the initial temperature Tðx; Þ is given (see Fig. 1). n the temperature Tðx; tþ is specified and on the flux P T ðx; tþ is specified Flux P T (x,t) specified Temperature T(x,) specified Fig. 1. The initial-boundary value problem. Temperature T(x,t) specified ¼ [, P T ¼ k ij ðt=x j Þn i with n i denoting the ith component of the outward pointing vector n normal to. Substitution of (2) into (1) yields Tðx; tþ Tðx; tþ m x ij ðxþhðtþ ¼ rðx; TÞcðx; TÞ. (3) i x j Use of the Kirchhoff s transformation Y ¼ hðtþ dt (4) in (3) gives Yðx; tþ Yðx; tþ m x ij ðxþ ¼½rðx; YÞcðx; YÞ=hðYÞŠ. i x j (5) The corresponding initial and boundary conditions for (5) are Yðx; Þ ¼YðTðx; ÞÞ for x 2, Yðx; tþ ¼YðTðx; tþþ for x 2, P Y ðx; tþ ¼P T ðx; tþ for x 2, (6) P Y ¼ m ij ðy=x j Þn i. (7) 4. Integral euation formulation The coefficients m ij ðxþ are reuired to take the form m ij ðxþ ¼k ij gðxþ, (8) the k ij are constants and gðxþ is reuired to satisfy 2 g 1=2 k ij ¼. (9) Substitution of (8) and use of the transformation cðx; tþ ¼g 1=2 ðxþyðx; tþ (1) in (5) yields 2 cðx; tþ cðx; tþ k ij Cðx; cþ ¼, (11) rðx; cþcðx; cþ Cðx; cþ ¼. (12) hðcþgðxþ Let ðx; tjn; tþ be the fundamental solution of the euation cðx; tþ 2 cðx; tþ C r k ij ¼, (13) C r is a constant reference value of Cðx; cþ. Hence ðx; tjn; tþ 2 ðx; tjn; tþ C r k ij ¼ dðx nþdðt tþ, (14) ðx; tjn; tþ 2 ðx; tjn; tþ C r þ k ij ¼ dðx nþdðt tþ. (15)
3 156 ARTICLE IN PRESS M.I. Azis, D.L. Clements / Engineering Analysis with Boundary Elements 32 (28) Multiply (11) by ðx; tjn; tþ and integrate 2 cðn; tþ ðx; tjn; tþ Cðn; cþ k ij dðnþ dt ¼ (16) and multiply (15) by cðn; tþ and integrate ðx; tjn; tþ cðx; tþ ¼ cðn; tþ C r 2 ðx; tjn; tþ þ k ij dðnþ dt. (17) Adding (16) to (17) yields cðx; tþ ¼ Cðn; cþðx; tjn; tþ ðx; tjn; tþ þ C r cðn; tþ dðnþ dt 2 ðx; tjn; tþ þ cðn; tþk ij 2 cðn; tþ ðx; tjn; tþk ij dðnþ dt. (18) Now ðx; tjn; tþ ½ðx; tjn; tþcðn; tþš ¼ cðn; tþ þ ðx; tjn; tþ so that ðx; tjn; tþ cðn; tþ ¼ ðx; tjn; tþ þ ½ðx; tjn; tþcðn; tþš. (19) Use of (19) in (18) provides cðx; tþ ¼ ½Cðn; cþ C r Šðx; tjn; tþ dðnþ dt 2 ðx; tjn; tþ þ cðn; tþk ij 2 cðn; tþ ðx; tjn; tþk ij dðnþ dt þ C r ½ðx; tjn; tþcðn; tþš dðnþ dt. (2) Now t C r ½cŠ dðnþ dt ¼ C r dðnþ ½cŠ dt ¼ C r ½cŠ t¼ dðnþ (21) since ðx; tjn; tþ¼ when t ¼ t. Also reen s theorem gives 2 2 c ck ij k ij dðnþ ¼ ½cP P c Š dðnþ, (22) ðx; tjn; tþ P ðx; tjn; tþ ¼k ij n j, x i (23) P c ðn; tþ ¼k ij n j. x i (24) Use of (21) and (22) in (2) gives ðxþcðx; tþ ¼ ½ðx; tjn; tþp c ðn; tþ cðn; tþp ðx; tjn; tþš dðnþ dt ½Cðn; cþ C r Šðx; tjn; tþ dðnþ dt þ C r ½ðx; tjn; tþcðn; tþš t¼ dðnþ, (25) ¼ ifxe [, ¼ 1ifx2, ¼ 1 2 if x 2 and has a continuously turning tangent at x. Consider the function Hðn; cþ ¼ ½Cðn; cþ C r Š dc, (26) ¼ c c ¼ðC C rþ c. (27) Substitution of (27) into the second integral in (25) yields ðc C r Þ c dðnþ dt ¼ dðnþ dt. (28) Use of integration by parts for the integral on the righthand side of (28) gives dt ¼½H Št¼t t¼ H dt. (29) The unknowns are now assumed to be constant, on each time element. Hence, following the work of oto and Suzuki in [7], if the unknown c (and hence the temperature T) is assumed to take on its value at the upper bound ðt ¼ tþ of the semi-open time interval ð; tš then H dt ½HŠ t¼t½š t¼t t¼. (3) Substitution of (3) into (29) provides dt ½H Št¼t t¼ ½HŠ t¼t½š t¼t t¼ ¼½HŠ t¼t t¼ ½Š t¼. (31)
4 M.I. Azis, D.L. Clements / Engineering Analysis with Boundary Elements 32 (28) So the integral euation (25) may be written in the form ðxþcðx; tþ ¼ ½C r cðn; Þ Hðn; cðn; tþþ þ Hðn; cðn; ÞÞŠðx; tjn; Þ dðnþ þ P c ðn; tþ ðx; tjn; tþ dt cðn; tþ P ðx; tjn; tþ dt dðnþ. (32) However, if the temperature is assumed to take on its values at the lower bound t ¼ of the semi-open time interval ½; tþ then dt ¼. (33) So the integral euation (25) may be written in the form ðxþcðx; tþ ¼ P c ðn; Þ ðx; tjn; tþ dt cðn; Þ P ðx; tjn; tþ dt dðnþ þ C r cðn; Þðx; tjn; Þ dðnþ. (34) In (32) and (34) the fundamental solution is given in Chang et al. [9]. For the case of two-dimensional problems it may be written in the form s ðx; tjn; tþ ¼ 4pzðt tþ exp C rr 2, (35) 4zðt tþ R and z are given by R 2 ¼ð _ x 1 _x 1 Þ 2 þð _ x 2 _x 2 Þ 2 and z ¼½k 11 þ 2_sk 12 þð_s 2 þ s 2 Þk 22 Š=2 with _x 1 ¼ x 1 þ _sx 2 ; _x 1 ¼ x 1 þ _sx 2, _x ¼ sx 2 ; _x 2 ¼ sx 2, _s and s are, respectively, the real and the positive imaginary parts of the complex root s of the uadratic k 11 þ 2k 12 s þ k 22 s 2 ¼. Now from (7), (8) and (1) it follows that the corresponding value of P Y is P Y ðn; tþ ¼ P g ðnþcðn; tþþp c ðn; tþg 1=2 ðnþ, (36) P g ðnþ ¼k ij g 1=2 x j n i. (37) Substitution of (1) and (36) in (32) and (34), respectively, gives ðxþg 1=2 ðxþyðx; tþ ¼ ½C r g 1=2 ðnþyðn; Þ Hðn; g 1=2 ðnþyðn; tþþ and þ Hðn; g 1=2 ðnþyðn; ÞÞŠðx; tjn; Þ dðnþ þ g 1=2 ðnþp Y ðn; tþ ðx; tjn; tþ dt Yðn; tþ g 1=2 ðnþ P ðx; tjn; tþ dt P g ðnþ ðx; tjn; tþ dt dðnþ (38) ðxþg 1=2 ðxþyðx; tþ ¼ C r g 1=2 ðnþyðn; Þðx; tjn; Þ dðnþ þ g 1=2 ðnþp Y ðn; Þ ðx; tjn; tþ dt Yðn; Þg 1=2 ðnþ P ðx; tjn; tþ dt P g ðnþ ðx; tjn; tþ dt dðnþ. (39) The integral euations (38) and (39) may be used to obtain the solution Y of E. (5) subject to the initial and boundary conditions (6). Solutions for the temperature T may then be obtained using (4). Es. (38) and (39) are both based on the assumption that for each time step the dependent variable c (and hence the temperature T) takes on a constant value on the time interval ð; tþ. In the case of E. (38) this constant value is taken to be the value at the upper bound of the semi-open time interval ð; tš and in the case of E. (39) the constant value is taken to be the value at the lower bound of the semi-open time interval ½; tþ. 5. Discretisation and numerical examples In this section only the integral euation (39) will be considered. It will be used to find numerical solutions to some two-dimensional problems. Let the domain be discretised into K cells k for k ¼ 1; 2;...; K, the boundary into J elements j for j ¼ 1; 2;...; J and the time t into M time intervals ½t m 1 ; t m Š for m ¼ 1; 2;...; M. If it is further assumed that the unknowns are constant along any boundary segment j ¼½ j 1 ; j Š and take on their values at the mid-point j, and also
5 158 ARTICLE IN PRESS M.I. Azis, D.L. Clements / Engineering Analysis with Boundary Elements 32 (28) constant in any domain cell k taking on their values at the center point r k of the cell. Then the discretised form of (39) for the time element ½t m 1 ; t m Š may be written as ðxþg 1=2 ðxþyðx; t m Þ ¼ XK Yðr k ; t m 1 Þ C r g 1=2 ðnþðx; t m jn; t m 1 Þ dðnþ k¼1 k ( þ XJ tm P Y ð j ; t m 1 Þ g 1=2 ðnþ ðx; t m jn; tþ dt dðnþ j¼1 j t m 1 " m Yð j ; t m 1 Þ g 1=2 ðnþ P ðx; t m jn; tþ dt dðnþ j t m 1 #) tm P g ðnþ ðx; t m jn; tþ dt dðnþ. (4) j t m 1 For each time interval m, (4) is treated as a new problem. It should be noted that the variable Y on the lefthand side of (4) is associated with the time t m as on the right-hand side it is associated with the time t m 1. Thus in order to solve (4) for the boundary unknowns (Y on and P Y on ) using an explicit scheme of linear euations the source point x has to be located outside the domain so that ¼ and the left-hand side of (4) is zero. Having solved for the boundary unknowns, solutions Y at the interior points r k for k ¼ 1; 2;...; K need to be found. These are then used as initial conditions for the time interval m þ 1. The time integrals of the fundamental solutions and P in (4) are evaluated analytically. Specifically m t m 1 ðx; t m jn; tþ dt ¼ s 4pz E 1 m C r R 2 4zðt m t m 1 Þ, (41) P ðx; t m jn; tþ dt t m 1 ¼ s 4pzR 2 exp C r R 2 R k 2 ij n j, (42) 4zðt m t m 1 Þ x i C r is assumed to be a positive number and E 1 is the exponential integral function defined by E 1 ðzþ 1 e s z s ds ¼ g ln z þ X1 n¼1 ð 1Þ n 1 z n nn!, (43) g ¼ : is Euler s constant. For all the problems considered the boundary of the domain is divided into a number of segments of eual length. Simpson s 3 8 rule is applied for the evaluation of the line integrals. For the area integrals, the domain is divided into a number of rectangular cells of eual area. The number of cells is increased to ensure the reuired level of accuracy for the domain integral. A uniform time step Dt ¼ t m t m 1 is used for all m. Also, when solving for the boundary unknowns, each source point is located outside the domain at a distance of one and a half times the length of the boundary segment under consideration. The distance is measured from the mid point of each segment in the direction of the outward pointing normal vector n Problem 1: Quadratic spatial variation of the conductivity matrix Consider the heat conduction problem governed by (1) for an inhomogeneous anisotropic material for which the coefficients of heat conduction are independent of temperature and are given by E. (8). The spatial variation of the coefficients is determined by the function gðxþ which is reuired to satisfy E. (9). Here gðxþ to taken to assume the uadratic form (Fig. 2) gðx 1 ; x 2 Þ¼ð1þ:1x 1 Þ2, (44) x i ¼ x i=d (for i ¼ 1; 2) with d a reference length. Also the coefficients occurring in E. (3) take the forms: ½m ij ðx 1 ; x 2 ÞŠ ¼ ½k ij Šgðx 1 ; x 2 Þ¼ 3 1 ð1 þ :1x Þ2, hðt Þ¼1, c ðx 1 ; x 2 ; T Þ¼ 4:5 T g1=2 ðx 1 ; x 2 Þ, r ðx 1 ; x 2 ; T Þ¼1, m ij ¼ m ij=m and k ij ¼ k ij=m with m a reference coefficient of heat conduction, T ¼ T=T with T a reference temperature, c ¼ c=c with c is a reference heat capacity and r ¼ r=r with r ¼ m h t =c d 2 a reference density. The geometry, initial and boundary conditions for the test problem under consideration are given in Fig. 3 T ðx 1 ; x 2 ; Þ ¼1 :25ðx 1 þ x 2 Þ2 ð1 þ :1x 1 Þ for px 1 p1 and px 2 p1, T ðx 1 ; 1; t Þ¼ 1 :25ð1 þ x 1 Þ2 ð1 þ t Þð1 þ :1x 1 Þ g(x 1,x 2 ) x Fig. 2. Spatial variation in the heat conductivities. x 1 1
6 M.I. Azis, D.L. Clements / Engineering Analysis with Boundary Elements 32 (28) for px 1 p1 and t X, T ðx 1 ; ; t 1 :25x 2 1 Þ¼ ð1 þ t Þð1 þ :1x 1 Þ for px 1 p1 and t X, P T ð1; x 2 ; t Þ¼ :3 2:2ð1 þ x 2 Þþ:75ð1 þ x 2 Þ2 1 þ t for px 2 p1 and t X, P T ð; x 2 ; t Þ¼ :3 þ 2x 2 :75x þ t for px 2 p1 and t X, t ¼ t=t, t is a reference time and P T ¼ P T=P with P ¼ k T =d a reference flux. It may be readily verified using the method of separation of variables that, for this particular problem, E. (3) admits the analytical solution T ðx 1 ; x 2 ; tþ ¼ 1 :25ðx 1 þ x 2 Þ2 ð1 þ t Þð1 þ :1x 1 Þ. (45) The corresponding heat flux is given by P T ðx 1; x 2 ; tþ ¼ :1ð3n 1 þ n 2 Þ 1 þ t ½1 :25ðx 1 þ x 2 Þ2 Š :5ð4n 1 þ 5n 2 Þ 1 þ t ðx 1 þ x 2 Þð1 þ :1x 1Þ. (46) Fig. 3. The boundary and the initial conditions for Problem 1. Since hðt Þ¼1, the variable Y is given by Y ¼ T Y ¼ Y=T. The function Cðx; cþ in (12) is given by C ðx; cþ ¼4:5=ðg 1=2 T Þ C ¼ C=C, C ¼ r c. As pg 1=2 T p1 it follows that C ðx; cþx4:5. It is convenient to take the reference value C r as C r ¼ 4:5 C r ¼ C r=c. The number of boundary segments and domain cells used is 16 (i.e. 4 4) and 16 (i.e. 4 4), respectively. Table 1 shows a comparison between the BEM and the analytical solutions at the interior point x ¼ð:5; :5Þ for some different values of the time step Dt. The results indicate the rate of convergence of the BEM solution to the analytical solution as the time step Dt decrerases Problem 2: Linear temperature variation of the conductivity matrix Consider the heat conduction problem for a homogeneous isotropic material with coefficients m ij ðx 1 ; x 2 Þ¼ 1 gðx 1 1 ; x 2 Þ, hðt Þ¼1 þ T, c ðx 1 ; x 2 ; T Þ¼1 þ :5T, r ðx 1 ; x 2 ; T Þ¼1, gðx 1 ; x 2 Þ¼1. The initial and boundary conditions are given in Fig. 4. Since hðt Þ¼1þT (4) gives Y ¼ T þ 1 2 T 2. Solving for T and taking the positive suare root yields p T ¼ 1þ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 2Y þ 1. The function Cðx; cþ is given as C ðx; cþ ¼ 1 þ :5T 1 þ T. In view of the boundary conditions the temperature T is expected to satisfy pt p1, so that 3 4 pc ðx; cþp1 and hence a reference value of C r ¼ :75 is appropriate. This particular problem is euivalent to the problem in Example 2 of the paper by oto and Suzuki [7] (for the case when k ¼ 1, b ¼ :5). oto and Suzuki showed that results they obtained using their boundary element techniue for t ¼ 1, x 2 ¼ :1 and values of x 1 given in Table 1 Solution T ð:5; :5; t Þ for Problem 1 Time BEM solution Analytical solution t Dt ¼ :1 Dt ¼ :2 Dt ¼ :4 Dt ¼ :1 Dt ¼ :
7 16 ARTICLE IN PRESS M.I. Azis, D.L. Clements / Engineering Analysis with Boundary Elements 32 (28) Conclusion Fig. 4. The initial and boundary conditions for Problem 2. Table 2 Solution T for Problem 2 A boundary element method for the solution of a certain class of non-linear parabolic initial-boundary value problems for a certain class of anisotropic inhomogeneous media has been derived. The method is generally easy to implement to obtain numerical values for particular problems. It can be applied to the time dependent heat conduction problems for inhomogeneous anisotropic materials and the analysis includes isotropy as a special case. The numerical results obtained using the method indicate that it can provide accurate numerical solutions for particular transient 2-D heat conduction problems. Position, x 1 oto and Suzuki, t ¼ 1 Present procedure, t ¼ 1 References Table 2 agreed with the numerical solution obtained for the same problem using several different numerical techniues. It is therefore reasonable to assume that the figures calculated by them provide a reference solution which is close enough to the exact solution to the problem for the values of x 1 ; x 2 and t indicated. By way of comparison the current techniue was employed to solve the same problem with the time step decreased and with the boundary elements and corresponding cells increased until convergence was achieved to the 1 values shown in Table 2 using a time step of 3, 96 boundary elements and 32 ð4 8Þ cells. The calculated values which are shown in Table 2 agree with the reference solution correct to two decimal places. For convergence to the given values the current explicit method based on the formula (39) reuires more boundary elements and cells than the implicit method of oto and Suzuki which is based on a special case of the formula (38). [1] Divo E, Kassab AJ. eneralized boundary integral euation for heat conduction in non-homogeneous media: recent developments on the sifting property. Eng Anal Bound Elem 1998;22: [2] Clements DL, Budhi WS, Azis MI. Boundary element methods for the solution of a class of elliptic boundary value problems for anisotropic inhomogeneous media. ANIAM J 23;44E: [3] Abdullah H, Clements DL, Ang WT. A boundary element method for the solution of a class of heat conduction problems for inhomogeneous media. Eng Anal Bound Elem 1993;11: [4] Sutradhar A, Paulino H, ray LJ. Transient heat conduction in homogeneous and non-homogeneous materials by the Laplace transform alerkin boundary element method. Eng Anal Bound Elem 22;26: [5] Clements DL, Budhi WS. A boundary element method for the solution of a class of steady-state problems for anisotropic media. ASME J Heat Transfer 1999;121: [6] Bialecki R, Kuhn. Boundary element solution of heat conduction prolems in multizone bodies of nonlinear material. Int J Numer Methods Eng 1993;36: [7] oto T, Suzuki M. A boundary integral euation method for nonlinear heat conduction problems with temperature-dependent material properties. Int J Heat Mass Transfer 1996;39: [8] Kikuta M, Togoh H, Tanaka M. Boundary element analysis of nonlinear transient heat conduction problems. Comput Methods Appl Math 1987;62: [9] Chang YP, Kang S, Chen DJ. The use of fundamental reen s functions for the solution of problems of heat conduction in anisotropic media. Int J Heat Mass Transfer 1973;16:
Numerical 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 informationA new meshless method for steady-state heat conduction problems in anisotropic and inhomogeneous media
Archive of Applied Mechanics 74 25 563--579 Springer-Verlag 25 DOI 1.17/s419-5-375-8 A new meshless method for steady-state heat conduction problems in anisotropic and inhomogeneous media H. Wang, Q.-H.
More informationExistence of positive periodic solutions for a periodic logistic equation
Applied Mathematics and Computation 139 (23) 311 321 www.elsevier.com/locate/amc Existence of positive periodic solutions for a periodic logistic equation Guihong Fan, Yongkun Li * Department of Mathematics,
More informationThis article was published in an Elsevier journal. The attached copy is furnished to the author for non-commercial research and education use, including for instruction at the author s institution, sharing
More informationInitial inverse problem in heat equation with Bessel operator
International Journal of Heat and Mass Transfer 45 (22) 2959 2965 www.elsevier.com/locate/ijhmt Initial inverse problem in heat equation with Bessel operator Khalid Masood *, Salim Messaoudi, F.D. Zaman
More informationDipartimento di Matematica, Politecnico di Torino, Torino 10129, Italy PLEASE SCROLL DOWN FOR ARTICLE
This article was downloaded by: [Camporesi, Roberto] On: 19 May 211 Access details: Access Details: [subscription number 936142537] Publisher Taylor & Francis Informa Ltd Registered in England and Wales
More informationTwo numerical methods for solving a backward heat conduction problem
pplied Mathematics and Computation 79 (6) 37 377 wwwelseviercom/locate/amc Two numerical methods for solving a backward heat conduction problem Xiang-Tuan Xiong, Chu-Li Fu *, hi Qian epartment of Mathematics,
More information1 Introduction Plane contact problems for homogeneous anisotropic materials have been studied quite widely in the past and analytical solutions have b
On the indentation of an inhomogeneous anisotropic elastic material by multiple straight rigid punches D. L. Clements a;λ, W.T. Ang b a School of Mathematics, University of Adelaide, SA 5005, Australia.
More informationInfiltration from irrigation channels into soil with impermeable inclusions
ANZIAM J. 46 (E) pp.c1055 C1068, 2005 C1055 Infiltration from irrigation channels into soil with impermeable inclusions Maria Lobo David L. Clements Nyoman Widana (Received 15 October 2004; revised 16
More informationTransient 3D Heat Conduction in Functionally. Graded Materials by the Method of Fundamental. Solutions
Transient 3D Heat Conduction in Functionally Graded Materials by the Method of Fundamental Solutions Ming Li*, C.S. Chen*, C.C. Chu, D.L.Young *Department of Mathematics, Taiyuan University of Technology,
More informationTESTS OF STATISTICAL SIGNIFICANCE AND BACKGROUND ESTIMATION IN GAMMA-RAY AIR SHOWER EXPERIMENTS
The Astrophysical Journal, 603:355 362, 2004 March 1 # 2004. The American Astronomical Society. All rights reserved. Printed in U.S.A. TESTS OF STATISTICAL SIGNIFICANCE AND BACKGROUND ESTIMATION IN GAMMA-RAY
More informationTransactions on Modelling and Simulation vol 8, 1994 WIT Press, ISSN X
Boundary element method for an improperly posed problem in unsteady heat conduction D. Lesnic, L. Elliott & D.B. Ingham Department of Applied Mathematical Studies, University of Leeds, Leeds LS2 9JT, UK
More informationExistence and uniqueness of solutions for a diffusion model of host parasite dynamics
J. Math. Anal. Appl. 279 (23) 463 474 www.elsevier.com/locate/jmaa Existence and uniqueness of solutions for a diffusion model of host parasite dynamics Michel Langlais a and Fabio Augusto Milner b,,1
More informationCharacteristics of forced convection heat transfer in vertical internally finned tube B
International Communications in Heat and Mass Transfer 32 (2005) 557 564 www.elsevier.com/locate/ichmt Characteristics of forced convection heat transfer in vertical internally finned tube B A. Al-Sarkhi*,
More informationAnalytical approximation of open-channel flow for controller design
Applied Mathematical Modelling 28 (2004) 677 695 www.elsevier.com/locate/apm Analytical approximation of open-channel flow for controller design Xavier Litrico a, *, Vincent Fromion b,1 a Cemagref, UR
More informationAn inverse vibration analysis of a tower subjected to wind drags on a shaking ground
Applied Mathematical Modelling 26 (2002) 517 528 www.elsevier.com/locate/apm An inverse vibration analysis of a tower subjected to wind drags on a shaking ground C.C. Kang a, C.Y. Lo b, * a Department
More informationHigh-order representation of Poincare maps
Nuclear Instruments and Methods in Physics Research A 558 (2006) 106 111 wwwelseviercom/locate/nima High-order representation of Poincare maps Johannes Grote, Martin Berz, Kyoko Makino Department of Physics
More informationA peak factor for non-gaussian response analysis of wind turbine tower
Journal of Wind Engineering and Industrial Aerodynamics 96 (28) 227 2227 www.elsevier.com/locate/jweia A peak factor for non-gaussian response analysis of wind turbine tower Luong Van Binh a,, Takeshi
More informationNumerical solution of a Fredholm integro-differential equation modelling _ h-neural networks
Available online at www.sciencedirect.com Applied Mathematics and Computation 195 (2008) 523 536 www.elsevier.com/locate/amc Numerical solution of a Fredholm integro-differential equation modelling _ h-neural
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 informationNON-HYPERELLIPTIC RIEMANN SURFACES OF GENUS FIVE ALL OF WHOSE WEIERSTRASS POINTS HAVE MAXIMAL WEIGHT1. Ryutaro Horiuchi. Abstract
R HORIUCHI KODAI MATH J 0 (2007), 79 NON-HYPERELLIPTIC RIEMANN SURFACES OF GENUS FIVE ALL OF WHOSE WEIERSTRASS POINTS HAVE MAXIMAL WEIGHT1 Ryutaro Horiuchi Abstract It is proved that any non-hyperelliptic
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat 15 (21) 3965 3973 Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns A similarity solution in
More informationHow do our representations change if we select another basis?
CHAPTER 6 Linear Mappings and Matrices 99 THEOREM 6.: For any linear operators F; G AðV Þ, 6. Change of Basis mðg FÞ ¼ mðgþmðfþ or ½G FŠ ¼ ½GŠ½FŠ (Here G F denotes the composition of the maps G and F.)
More informationMATH 260 Homework assignment 8 March 21, px 2yq dy. (c) Where C is the part of the line y x, parametrized any way you like.
MATH 26 Homework assignment 8 March 2, 23. Evaluate the line integral x 2 y dx px 2yq dy (a) Where is the part of the parabola y x 2 from p, q to p, q, parametrized as x t, y t 2 for t. (b) Where is the
More informationMFS with RBF for Thin Plate Bending Problems on Elastic Foundation
MFS with RBF for Thin Plate Bending Problems on Elastic Foundation Qing-Hua Qin, Hui Wang and V. Kompis Abstract In this chapter a meshless method, based on the method of fundamental solutions (MFS and
More information1 Basic tools. 1.1 The Fourier transform. Synopsis. Evenness and oddness. The Fourier transform of f(x) is defined as Z 1
1 Basic tools 1.1 The Fourier transform 1 Synopsis The Fourier transform of f(x) is defined as FðsÞ ¼ f ðxþe i2pxs dx The inverse Fourier transform is given by f ðxþ ¼ FðsÞe þi2pxs ds Evenness and oddness
More informationA two armed bandit type problem*
Int J Game Theory (2003) 32:3 16 DOI 10.1007/s001820300135 A two armed bandit type problem* Michel Benaïm1, Gerard Ben Arous2 1 Department of Mathematics, Université Cergy-Pontoise, 2 Avenue A. Chauvin,
More informationThe Eulerian Lagrangian method of fundamental solutions for two-dimensional unsteady Burgers equations
Engineering Analysis with Boundary Elements 32 (2008) 395 412 www.elsevier.com/locate/enganabound The Eulerian Lagrangian method of fundamental solutions for two-dimensional unsteady Burgers equations
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationSequences, Series, and the Binomial Formula
CHAPTER Sequences, Series, nd the Binomil Formul. SEQUENCES. ; ; ; ; 6 ; 6 6. ðþ ; ðþ ; ð Þ 6; ðþ ; ðþ 6; 6 ð6þ. ðþ ; ðþ : ðþ ; ðþ ; ðþ ; 6 ðþ 6 6 6. ; ; ; ; ; 6 6 6. 0 ; ; ; 8 ; 6 8 ; 6. 0; ; 6 ; ; 6
More informationTHE WAVE EQUATION. d = 1: D Alembert s formula We begin with the initial value problem in 1 space dimension { u = utt u xx = 0, in (0, ) R, (2)
THE WAVE EQUATION () The free wave equation takes the form u := ( t x )u = 0, u : R t R d x R In the literature, the operator := t x is called the D Alembertian on R +d. Later we shall also consider the
More informationApplied Mathematics and Computation
Applied Mathematics and Computation 245 (2014) 86 107 Contents lists available at ScienceDirect Applied Mathematics and Computation journal homepage: www.elsevier.com/locate/amc An analysis of a new family
More information1 Constant Pressure Inner Boundary Condition
Jim Lambers ENERGY 281 Spring Quarter 2007-08 Lecture 8 Notes These notes are based on Rosalind Archer s PE281 lecture notes, with some revisions by Jim Lambers. 1 Constant Pressure Inner Boundary Condition
More informationSimplified marginal effects in discrete choice models
Economics Letters 81 (2003) 321 326 www.elsevier.com/locate/econbase Simplified marginal effects in discrete choice models Soren Anderson a, Richard G. Newell b, * a University of Michigan, Ann Arbor,
More informationSolvability condition in multi-scale analysis of gyroscopic continua
Journal of Sound and Vibration 309 (2008) 338 342 Short Communication Solvability condition in multi-scale analysis of gyroscopic continua Li-Qun Chen a,b,, Jean W. Zu c JOURNAL OF SOUN AN VIBRATION a
More informationMathematics (JUN11MPC301) General Certificate of Education Advanced Level Examination June Unit Pure Core TOTAL
Centre Number Candidate Number For Examiner s Use Surname Other Names Candidate Signature Examiner s Initials Mathematics Unit Pure Core 3 Monday 13 June 2011 General Certificate of Education Advanced
More informationAn Explicit Formula for the Discrete Power Function Associated with Circle Patterns of Schramm Type
Funkcialaj Ekvacioj, 57 (2014) 1 41 An Explicit Formula for the Discrete Power Function Associated with Circle Patterns of Schramm Type By Hisashi Ando1, Mike Hay2, Kenji Kajiwara1 and Tetsu Masuda3 (Kyushu
More informationTransient conduction and convection heat transfer across a multi-layer floor subjected to multiple heat sources
ARTICLE IN PRESS Building and Environment 41 (26) 1299 131 www.elsevier.com/locate/buildenv Transient conduction and convection heat transfer across a multi-layer floor subjected to multiple heat sources
More informationThree Boundary Meshless Methods for Heat Conduction Analysis in Nonlinear FGMs with Kirchhoff and Laplace Transformation
Advances in Applied Mathematics and Mechanics Adv. Appl. Math. Mech., Vol. 4, No. 5, pp. 519-54 DOI: 10.408/aamm.10-m1170 October 01 Three Boundary Meshless Methods for Heat Conduction Analysis in Nonlinear
More informationENGIN 211, Engineering Math. Laplace Transforms
ENGIN 211, Engineering Math Laplace Transforms 1 Why Laplace Transform? Laplace transform converts a function in the time domain to its frequency domain. It is a powerful, systematic method in solving
More informationAsymptotics of minimax stochastic programs $
Statistics & Probability Letters 78 (2008) 150 157 www.elsevier.com/locate/stapro Asymptotics of minimax stochastic programs $ Alexander Shapiro School of Industrial and Systems Engineering, Georgia Institute
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat xxx ) xxx xxx Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: wwwelseviercom/locate/cnsns A note on the use of Adomian
More informationEnergy Conversion and Management
Energy Conversion and Management 58 (202) 63 70 Contents lists available at SciVerse ScienceDirect Energy Conversion and Management journal homepage: www.elsevier.com/locate/enconman Computation of effectiveness
More informationConvergence of a linear recursive sequence
int. j. math. educ. sci. technol., 2004 vol. 35, no. 1, 51 63 Convergence of a linear recursive sequence E. G. TAY*, T. L. TOH, F. M. DONG and T. Y. LEE Mathematics and Mathematics Education, National
More informationChapter 2 Green s Functions for Laplace and Wave Equations
Chapter 2 Green s Functions for Laplace and Wave Equations This chapter shows the solution method for Green s functions of, 2 and 3D Laplace and wave equations. Lengthy and detailed explanations are given
More informationThermal field in water pipe cooling concrete hydrostructures simulated with singular boundary method
Water Science and Engineering 017, 10(): 107e114 HOSTED BY Available online at www.sciencedirect.com Water Science and Engineering journal homepage: http://www.waterjournal.cn Thermal field in water pipe
More informationMathematical Model for Thermal Shock Problem of a Generalized Thermoelastic Layered Composite Material with Variable Thermal Conductivity
COMPUTATIONAL METHODS IN SCIENCE AND TECHNOLOGY ( 65-7 (6 Mathematical Model for Thermal Shock Problem of a Generalized Thermoelastic Layered Composite Material with Variable Thermal Conductivity H. M.
More informationA numerical analysis of a model of growth tumor
Applied Mathematics and Computation 67 (2005) 345 354 www.elsevier.com/locate/amc A numerical analysis of a model of growth tumor Andrés Barrea a, *, Cristina Turner b a CIEM, Universidad Nacional de Córdoba,
More informationSinusoidal Steady- State Circuit Analysis
Sinusoidal Steady- State Circuit Analysis 9. INTRODUCTION This chapter will concentrate on the steady-state response of circuits driven by sinusoidal sources. The response will also be sinusoidal. For
More informationTHREE-DIMENSIONAL FINITE DIFFERENCE MODEL FOR TRANSPORT OF CONSERVATIVE POLLUTANTS
Pergamon Ocean Engng, Vol. 25, No. 6, pp. 425 442, 1998 1998 Elsevier Science Ltd. All rights reserved Printed in Great Britain 0029 8018/98 $19.00 + 0.00 PII: S0029 8018(97)00008 5 THREE-DIMENSIONAL FINITE
More informationElectronic energy transfer in polymers labeled at both ends with fluorescent groups
Journal of Luminescence 96 (22) 269 278 Electronic energy transfer in polymers labeled at both ends with fluorescent groups E.. Bodunov, M.. Berberan-Santos, J.M.G. Martinho* Centro de Qu!ımica-F!ısica
More informationOn Behaviors of the Energy of Solutions for Some Damped Nonlinear Hyperbolic Equations with p-laplacian Soufiane Mokeddem
International Journal of Advanced Research in Mathematics ubmitted: 16-8-4 IN: 97-613, Vol. 6, pp 13- Revised: 16-9-7 doi:1.185/www.scipress.com/ijarm.6.13 Accepted: 16-9-8 16 cipress Ltd., witzerland
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationPUBLICATIONS. Journal of Geophysical Research: Oceans. On Adler space-time extremes during ocean storms RESEARCH ARTICLE 10.
PUBLICATIONS RESEARCH ARTICLE Key Points: The statistical properties of extreme sea waves in space-time domain are examined The exceedance probability of maximum wave elevation over an area is obtained
More informationSec. 7.4: Basic Theory of Systems of First Order Linear Equations
Sec. 7.4: Basic Theory of Systems of First Order Linear Equations MATH 351 California State University, Northridge April 2, 214 MATH 351 (Differential Equations) Sec. 7.4 April 2, 214 1 / 12 System of
More informationwe get y 2 5y = x + e x + C: From the initial condition y(0) = 1, we get 1 5 = 0+1+C; so that C = 5. Completing the square to solve y 2 5y = x + e x 5
Math 24 Final Exam Solution 17 December 1999 1. Find the general solution to the differential equation ty 0 +2y = sin t. Solution: Rewriting the equation in the form (for t 6= 0),we find that y 0 + 2 t
More informationSynchronization of self-sustained oscillators by common white noise
Physica A 351 (5) 16 13 www.elsevier.com/locate/physa Synchronization of self-sustained oscillators by common white noise D.S. Goldobin a,b, A.S. Pikovsky a, a Department of Physics, University of Potsdam,
More informationOn the torsion of functionally graded anisotropic linearly elastic bars
IMA Journal of Applied Mathematics (2007) 72, 556 562 doi:10.1093/imamat/hxm027 Advance Access publication on September 25, 2007 edicated with admiration to Robin Knops On the torsion of functionally graded
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 informationBiological population model and its solution by reduced differential transform method
Asia Pacific Journal of Engineering Science and Technology () (05) -0 Asia Pacific Journal of Engineering Science and Technology journal homepage: www.apjest.com Full length article Biological population
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
More informationFrequency and time domain dynamic structural analysis: convergence and causality
Computers and Structures 80 (2002) 1503 1509 www.elsevier.com/locate/compstruc Frequency and time domain dynamic structural analysis: convergence and causality F. Venancio Filho a, *, A.M. Claret b, F.S.
More informationA simple Galerkin boundary element method for three-dimensional crack problems in functionally graded materials
Materials Science Forum Vols. 492-493 (2005) pp 367-372 Online available since 2005/Aug/15 at www.scientific.net (2005) Trans Tech Publications, Switzerland doi:10.4028/www.scientific.net/msf.492-493.367
More informationARTICLE IN PRESS. Engineering Analysis with Boundary Elements
Engineering Analysis with Boundary Elements 33 (29) 561 571 Contents lists available at ScienceDirect Engineering Analysis with Boundary Elements journal homepage: www.elsevier.com/locate/enganabound BEM
More informationINTERIOR POINT ALGORITHMS FOR NONLINEAR CONSTRAINED LEAST SQUARES PROBLEMS
Inverse Problems in Engng., 00?, Vol. 00, No. 0, pp. 1 13 INTERIOR POINT ALGORITHMS FOR NONLINEAR CONSTRAINED LEAST SQUARES PROBLEMS JOSE HERSKOVITS a, *, VERANISE DUBEUX a, CRISTOVA O M. MOTA SOARES b
More informationThe Approximate Solution of Non Linear Fredholm Weakly Singular Integro-Differential equations by Using Chebyshev polynomials of the First kind
AUSTRALIA JOURAL OF BASIC AD APPLIED SCIECES ISS:1991-8178 EISS: 2309-8414 Journal home page: wwwajbaswebcom The Approximate Solution of on Linear Fredholm Weakly Singular Integro-Differential equations
More informationElectrostatics: Electrostatic Devices
Electrostatics: Electrostatic Devices EE331 Electromagnetic Field Theory Outline Laplace s Equation Derivation Meaning Solving Laplace s equation Resistors Capacitors Electrostatics -- Devices Slide 1
More informationOn the design of reactionless 3-DOF planar parallel mechanisms
Mechanism and Machine Theory 41 (2006) 70 82 Mechanism and Machine Theory www.elsevier.com/locate/mechmt On the design of reactionless 3-DOF planar parallel mechanisms Abbas Fattah *, Sunil K. Agrawal
More informationFINITE VOLUME METHOD: BASIC PRINCIPLES AND EXAMPLES
FINITE VOLUME METHOD: BASIC PRINCIPLES AND EXAMPLES SHRUTI JAIN B.Tech III Year, Electronics and Communication IIT Roorkee Tutors: Professor G. Biswas Professor S. Chakraborty ACKNOWLEDGMENTS I would like
More information35-959, Rzeszów, Poland b Institute of Computer Science, Jagiellonian University,
This article was downloaded by: [Polska Akademia Nauk Instytut Matematyczny] On: 07 March 03, At: 03:7 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 0795 Registered
More information(ii) Given that gðxþ ¼ð4x 1Þðax 2 þ bx þ cþ, find the values of the integers a, b and c. (3 marks)
Centre Number Candidate Number For Examiner s Use Surname Other Names Candidate Signature Examiner s Initials Mathematics Unit Pure Core 4 Tuesday 15 June 2010 General Certificate of Education Advanced
More informationMath 132, Spring 2003 Exam 2: Solutions
Math 132, Spring 2003 Exam 2: Solutions No calculators with a CAS are allowed. e sure your calculator is set for radians, not degrees if you do any calculus computations with trig functions. Part I, Multiple
More informationHigher Order Averaging : periodic solutions, linear systems and an application
Higher Order Averaging : periodic solutions, linear systems and an application Hartono and A.H.P. van der Burgh Faculty of Information Technology and Systems, Department of Applied Mathematical Analysis,
More informationA note on Two-warehouse inventory model with deterioration under FIFO dispatch policy
Available online at www.sciencedirect.com European Journal of Operational Research 190 (2008) 571 577 Short Communication A note on Two-warehouse inventory model with deterioration under FIFO dispatch
More informationEuropean Journal of Operational Research
European Journal of Operational Research 219 (2012) 719 726 Contents lists available at SciVerse ScienceDirect European Journal of Operational Research journal homepage: www.elsevier.com/locate/ejor imax
More informationNumerical Analysis In Short-Pulse Laser Heating Of Metals Based On The Dual-Phase-Lag Model Of Heat Conduction
Proceedings of the World Congress on New echnologies (Newech 5) Barcelona, Spain July 5-7, 5 Paper No. Numerical Analysis In Short-Pulse Laser Heating Of Metals Based On he Dual-Phase-Lag Model Of Heat
More informationNonlinear structural dynamical system identification using adaptive particle filters
Journal of Sound and Vibration 306 (2007) 524 563 JOURNAL OF SOUND AND VIBRATION www.elsevier.com/locate/jsvi Nonlinear structural dynamical system identification using adaptive particle filters Vikas
More informationCVBEM for a System of Second Order Elliptic Partial Differential Equations
CVBEM for a System of Second Order Elliptic Partial Differential Equations W. T. Ang and Y. S. Par Faculty of Engineering, Universiti Malaysia Sarawa, 94300 Kota Samarahan, Malaysia Abstract A boundary
More informationPhase-Field Simulation of the Effect of Elastic Inhomogeneity on Microstructure Evolution in Ni-Based Superalloys
Materials Transactions, Vol. 50, No. 4 (2009) pp. 744 to 748 #2009 The Japan Institute of Metals Phase-Field Simulation of the Effect of Elastic Inhomogeneity on Microstructure Evolution in Ni-Based Superalloys
More informationViscosity approximation methods for nonexpansive nonself-mappings
J. Math. Anal. Appl. 321 (2006) 316 326 www.elsevier.com/locate/jmaa Viscosity approximation methods for nonexpansive nonself-mappings Yisheng Song, Rudong Chen Department of Mathematics, Tianjin Polytechnic
More informationAiry Functions. Chapter Introduction
Chater 4 Airy Functions 4. Introduction Airy functions are named after the English astronomer George Biddell Airy (8 892). Airy s first mathematical work was on the diffraction henomenon, namely, the Airy
More informationElastic solutions for stresses in a transversely isotropic half-space subjected to three-dimensional buried parabolic rectangular loads
INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS Int. J. Numer. Anal. Meth. Geomech., 2002; 26:1449 1476 (DOI: 10.1002/nag.253) Elastic solutions for stresses in a transversely
More informationThe Method of Fundamental Solutions with Eigenfunction Expansion Method for Nonhomogeneous Diffusion Equation
The Method of Fundamental Solutions with Eigenfunction Expansion Method for Nonhomogeneous Diffusion Equation D. L. Young, 1 C. W. Chen, 1 C. M. Fan, 1 C. C. Tsai 2 1 Department of Civil Engineering and
More informationFinite Difference Solution of the Heat Equation
Finite Difference Solution of the Heat Equation Adam Powell 22.091 March 13 15, 2002 In example 4.3 (p. 10) of his lecture notes for March 11, Rodolfo Rosales gives the constant-density heat equation as:
More informationMathematics (JUN10MFP301) General Certificate of Education Advanced Level Examination June Unit Further Pure TOTAL
Centre Number Candidate Number For Examiner s Use Surname Other Names Candidate Signature Examiner s Initials Mathematics Unit Further Pure 3 Friday 11 June 2010 General Certificate of Education Advanced
More informationPractice Midterm Solutions
Practice Midterm Solutions Math 4B: Ordinary Differential Equations Winter 20 University of California, Santa Barbara TA: Victoria Kala DO NOT LOOK AT THESE SOLUTIONS UNTIL YOU HAVE ATTEMPTED EVERY PROBLEM
More information1 Solution to Homework 4
Solution to Homework Section. 5. The characteristic equation is r r + = (r )(r ) = 0 r = or r =. y(t) = c e t + c e t y = c e t + c e t. y(0) =, y (0) = c + c =, c + c = c =, c =. To find the maximum value
More informationMATH 425, FINAL EXAM SOLUTIONS
MATH 425, FINAL EXAM SOLUTIONS Each exercise is worth 50 points. Exercise. a The operator L is defined on smooth functions of (x, y by: Is the operator L linear? Prove your answer. L (u := arctan(xy u
More informationApproximate step response of a nonlinear hydraulic mount using a simplified linear model
Journal of Sound and Vibration 99 (007) 656 663 Short Communication JOURNAL OF SOUND AND VIBRATION Approximate step response of a nonlinear hydraulic mount using a simplified linear model Song He, Rajendra
More informationNAG Library Routine Document D02JAF.1
D02 Ordinary Differential Equations NAG Library Routine Document Note: before using this routine, please read the Users Note for your implementation to check the interpretation of bold italicised terms
More informationLearning for non-stationary Dirichlet processes. Czech Republic, Czech Republic SUMMARY
ACS : 4 PROD.TYPE: COM ED: KAVYASHREE pp.^ (col. g.: NIL) B2 PAGN: G.PRAKASH SCAN: INTERNATIONAL JOURNAL OF ADAPTIVE CONTROL AND SIGNAL PROCESSING Int. J. Adapt. Control Signal Process. 200; 20:000 000
More informationarxiv: v1 [math.ap] 11 Jan 2014
THE UNIFIED TRANSFORM FOR THE MODIFIED HELMHOLTZ EQUATION IN THE EXTERIOR OF A SQUARE A. S. FOKAS AND J. LENELLS arxiv:4.252v [math.ap] Jan 24 Abstract. The Unified Transform provides a novel method for
More informationStabilization of a solid propellant rocket instability by state feedback
INTERNATIONAL JOURNAL OF ROBUST AND NONLINEAR CONTROL Int. J. Robust Nonlinear Control 23; 13:483 495 (DOI: 1.12/rnc.732) Stabilization of a solid propellant rocket instability by state feedback Dejan
More informationInput output linearization with delay cancellation for nonlinear delay systems: the problem of the internal stability
INTERNATIONAL JOURNAL OF ROBUST AND NONLINEAR CONTROL Int J Robust Nonlinear Control 23; 13:99 937 (DOI: 112/rnc83) Input output linearization with delay cancellation for nonlinear delay systems: the problem
More informationJournal of Computational Physics
Journal of Computational Physics 8 (009) 7565 7595 Contents lists available at ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp An accurate and efficient method
More informationStability analysis for a class of Takagi Sugeno fuzzy control systems with PID controllers
International Journal of Approximate Reasoning 46 (2007) 109 119 www.elsevier.com/locate/ijar Stability analysis for a class of Takagi Sugeno fuzzy control systems with PID controllers Le Hung Lan * Institute
More informationMicroeconomics CHAPTER 2. THE FIRM
Chapter The Firm Exercise. Suppose that a unit of output can be produced by any of the following combinations of inputs z = 0: 0:5 ; z 0:3 = 0:. Construct the isouant for =. ; z 3 0:5 = 0:. Assuming constant
More informationSuction Potential and Water Absorption from Periodic Channels in Different Types of Homogeneous Soils
I. Solekhudin and K.C. Ang/ Electronic Journal of Boundary Elements, Vol. 1, No. 2, pp. 42-55 (212) Suction Potential and Water Absorption from Periodic Channels in Different Types of Homogeneous Soils
More informationOn the Cauchy problem for non linear PDEs in the Gevrey class with shrinkings
J. Math. Soc. Japan Vol. 54, No. 3, 22 On the Cauchy problem for non linear PDEs in the Gevrey class with shrinkings By Masaki Kawagishi Takesi Yamanaka (Received May 29, 2) (Revised Jan. 5, 2) Abstract.
More informationNumerical methods to compute optical errors due to stress birefringence
Numerical methods to compute optical errors due to stress birefringence Keith B. Doyle Optical Research Associates, 8 West Park Drive, Westborough, MA Victor L. Genberg & Gregory J. Michels Sigmadyne,
More information