NEW IMPLEMENTATION OF MLBIE METHOD FOR HEAT CONDUCTION ANALYSIS IN FUNCTIONALLY GRADED MATERIALS
|
|
- James Martin
- 5 years ago
- Views:
Transcription
1 NEW IMPLEMENTATION OF MLBIE METHOD FOR HEAT CONDUCTION ANALYSIS IN FUNCTIONALLY GRADED MATERIALS DAVOUD MIRZAEI,,, MEHDI DEHGHAN Abstract. The meshless local boundary integral equation (MLBIE) method with an efficient technique to deal with the time variable are presented in this article to analyze the transient heat conduction in continuously nonhomogeneous functionally graded materials (FGMs). In space, the method is based on the local boundary integral equations and the moving least squares (MLS) approximation of the temperature and heat flux. In time, again the MLS approximates the equivalent Volterra integral equation derived from the heat conduction problem. It means that, the MLS is used for approximation in both time and space domains, and we avoid using the finite difference discretization or Laplace transform methods to overcome the time variable. Finally the method leads to a single generalized Sylvester equation rather than some (many) linear systems of equations. The method is computationally attractive, which is shown in couple of numerical examples for a finite strip and a hollow cylinder with an exponential spatial variation of material parameters. Keywords. Meshless local boundary integral equation (MLBIE) method; Moving least squares (MLS) approximation; Heat conduction problem; Functionally graded materials (FGM).. Introduction Functionally graded materials (FGMs) are those in which the composition and the volume fraction of FGM constituents vary gradually, giving a non-uniform microstructure with continuously graded macro-properties (e.g. specific heat, conductivity, density). For instance, one face of a structural component (plate or cylinder) may be an engineering ceramic that can resist severe thermal loading, and the other face may be a metal to maintain structural strength and toughness [, 2]. In recent Date: October 6, 2. Corresponding author. address: d.mirzaei@sci.ui.ac.ir (D. Mirzaei). address: mdehghan@aut.ac.ir (M. Dehghan). Department of Mathematics, University of Isfahan, , Isfahan, Iran. Tel: Department of Applied Mathematics, Faculty of Mathematics and Computer Science, Amirkabir University of Technology, No. 424, Hafez Ave.,594, Tehran, Iran.
2 2 DAVOUD MIRZAEI,,, MEHDI DEHGHAN years FGMs have drawn considerable interest [3]. These FGMs are expected to be highly heat resistant materials that can be used under high temperature and high temperature gradient conditions. In this case, it is important to investigate the temperature distributions in the FGMs. Due to the high mathematical complexity of the initial boundary value problems, analytical approaches for the thermo-mechanics of FGMs are restricted to simple geometry and boundary conditions. Thus, the transient heat conduction analysis in FGM demands accurate and efficient numerical methods [4]. There are plenty of papers considered the heat conduction problem in FGMs. Here we present some of those papers which are more related to the meshless methods. Sutradhar et.al [] presented the Galerkin boundary element method with Laplace transform for heat conduction problem in homogeneous and nonhomogeneous FGMs. J. Sladek et.al [5] proposed a meshless method for transient heat conduction analysis in FGMs using local boundary integral equation. In [6] they used this technique for heat conduction with nonlinear source term in FGMs. V. Sladek et.al [7] developed a local boundary integral method for transient heat conduction in anisotropic and functionally graded media. Again, Sladek brothers and their collaborates employed meshless local Petrov-Galerkin (MLPG) method to analyze heat conduction in FGMs [8, 4]. Wang et.al [9] developed a meshless numerical method based on RBF interpolation and time discretization to analyze this problem. Hamza-Cherif et.al [] solved transient temperature fields in axisymmetric FGM cylinders under various boundary conditions using an h-p finite element method. There are also many works on meshless methods for solving some other timedependent problems. For example see [, 2, 3, 4, 5] and etc. Meshless methods are based on approximation in terms of scattered data [6]. These methods are becoming popular, due to their high adaptivity and a low cost to prepare input data for numerical analysis. The moving least squares (MLS) is one of the scattered data approximation methods, which has been used successfully to approximate the trial space in many meshless methods. For instance we can mention the MLBIE method [7, 8]. This method is based on local sub-domains, rather than a global domain, and requires neither domain elements nor background cells in either the approximation or the integration. Hence it is a truly weak-based meshless method. During these recent years, MLBIE methods have been widely used to solve many steady-state and time-dependent engineering problems. The authors of all above-mentioned papers employed either the Laplace transform or the time discretization method to overcome the time variable. In this paper, we propose an interesting meshless local boundary integral equation (ML- BIE) method which employs the MLS to approximate the quantities in both time
3 2 3 and space domains. Although the technique is presented for the MLBIE method, it can be easily extended to other MLS based meshless methods. This paper is organized as follows. In Section 2, MLS approximation is discussed. In Section 3, the governing equations of heat conduction problem in FGMs are converted to an algebraic system of equations using MLBIE method. In Section 4, numerical results for a finite stirp and a hallow cylinder are presented. 2. MLS approximation The MLS as an approximation method has been introduced by Lancaster and Salkauskas [9] inspired by the pioneer work of Shepard [2]. Since the numerical approximations of MLS are based on a cluster of scattered nodes instead of interpolation on elements, there have many meshless methods based on the MLS method for the numerical solution of differential equations in recent years. Given data values {u j } N j= at centers {x,..., x N }, the MLS method produces a function û that approximates u in a weighted square sense. Let P d m be the space of polynomial of degree m in R d, and let B = {p, p 2,..., p Q } be any polynomial basis of P d m, where Q = ( ) m+d d. The MLS approximation û(x) of u(x), x Ω, can be defined as û(x) = p T (x)a(x), x Ω, where p T (x) = [ p (x), p 2 (x),..., p Q (x) ] and a(x) is a vector with components a j (x), j =, 2,..., Q, which are functions of the space coordinates x. In this paper we suppose that the basis of P d m is a complete monomial basis of order Q. The coefficient vector a(x) is determined by minimizing a weight discrete square norm, which is defined as I(x) = N w(x, x j ) ( p T ) 2, (x j )a(x) u j (2.) j= where w(x, x j ) is the weight function associated with the node j. In MLS, we are interested in local continuous weight function w. To be more precise we choose a continuous function φ : [, ) [, ) with and define φ(r) >, r <, φ(r) =, r, ( ) x y 2 w δ (x, y) = φ, δ for δ > as a weight function. We further define the set of indices J(x) = J(x, δ, X) = {j {, 2,..., N} : x x j 2 δ}, of centers contained in the interior of closed ball B(x, δ) of radius δ around x.
4 4 DAVOUD MIRZAEI,,, MEHDI DEHGHAN If we define the matrices P and W as P := ( p l (x j ) ), l Q, j J(x), W (x) := ( δ jk w δ (x, x k ) ) j,k J(x), the stationarity point of I, in equation (2.), with respect to a(x) leads to and ϕ(x) = p(x)a B, (2.2) û(x) = N ϕ j (x)u(x j ), (2.3) j= where A(x) = P W (x)p T, B(x) = P W (x) and ϕ(x) = [ϕ (x),..., ϕ N (x)], providing ϕ j (x) = for j / J(x). The functions ϕ j (x) are called the shape functions of the MLS approximation, corresponding to nodal point x j. If w(x, x j ) C k (Ω), then ϕ j (x) C k (Ω) [2, 22]. The matrix A is of size Q Q. The MLS approximation is well-defined only when the moment matrix A is non-singular. It can be seen that this is the case if and only if the rank of P equals Q (P should be a full-rank matrix). A necessary condition is J(x) Q, i.e. there should be at least Q points in the support of weight function w δ (the domain of definition of MLS approximation), but this condition is not sufficient. A necessary and sufficient condition is the P d m-unisolvency condition which is defined on set {x,..., x N } and depends on the geometric shape of Ω. For more precise details see [2] and Chapters 3 and 4 of [22]. The partial derivatives of û(x) are obtained as: where ϕ j,i (x) = û,i (x) = Q k= N ϕ j,i (x)u(x j ), x Ω, (2.4) j= ( pk,i [A B] kj + p k [A B,i + A,i B] kj), where A,i = (A ),i represents the derivative of the inverse of A with respect to i-th component of x, which is given by A,i = A A,i A. In the MLS and other scattered data approximation methods two quantities fill distance and separation distance are important to measure the quality of centers and derive the rate of convergence. For a set of points X = {x, x 2,..., x N } in a bounded domain Ω R d the fill distance is defined to be h X,Ω = sup min x x j 2, j N x Ω
5 2 5 and the separation distance is defined by q X = 2 min i j x i x j 2. A set X of data sites is said to be quasi-uniform with respect to a constant c qu > if q X h X,Ω c qu q X. (2.5) The error analysis of MLS approximation can be found in [23, 2, 22]. If Ω satisfies an interior cone condition, then there exist constants C, such that for all u C m+ (Ω ) and all quasi-uniform X Ω, we have u û Ch m+ X,Ω u C m+ (Ω ), (2.6) where Ω is the closure of x Ω B(x, Ch ) with h X,Ω h and h depends on Ω. The semi-norm in the right hand side is defined by u Ck (Ω ) := max β =k Dβ u L (Ω). The order of convergence of MLS approximation mainly depends on the polynomial basis function B and the quality of point set X. It dose not depend on the weight function. The weight function is used to make a moving, local and consequently stable approximation. Here, the following Gaussian weight function is used: { exp[ (dj/c) 2 ] exp[ (δ/c) 2 ] w δ (x, x j ) = exp[ (δ/c) 2 ], d j δ,, d j > δ, where d j = x x j 2, c = c h X,Ω is a constant controlling the shape of the weight function and δ = δ h X,Ω is the size of the support domain, with constants c and δ. 3. Transient heat conduction in FGM The thermal heat conduction problem in the FGM is described by the following governing equation [5] u (x, t) = u(x, t) + κ(x) u(x, t) + f(x, t), (3.) α(x) t κ(x) κ(x) where x Ω and t t F denote the time and the space variables, respectively, and t F is the finial time. The initial and boundary conditions are u(x, ) = u (x), x Ω, (3.2) u(x, t) = u(x, t), x Γ u, t t F, (3.3) κ(x) u n (x, t) = q(x, t), x Γ q, t t F. (3.4) In (3.)-(3.4), represents the Laplacian operator, u(x, t) is the temperature field, κ(x) and α(x) stand for the thermal conductivity and diffusivity, respectively, and f(x, t) is the density of the body heat sources. Moreover n is the unit outward
6 6 DAVOUD MIRZAEI,,, MEHDI DEHGHAN normal to the boundary Γ, u and q are specified values on the Dirichlet boundary Γ u and Neumann boundary Γ q where Γ = Γ u Γ q. To employ the MLS approximation for the time domain, we do the following: from (3.) and (3.2) we have u(x, t) = α t [ u(x, τ) + κ κ u(x, τ) ] dτ + κ t f(x, τ)dτ + α u (x), (3.5) keeping in mind α := α(x) and κ := κ(x) are some functions of spatial variable. First the MLS approximation is written in respect to the time variable t. The idea comes up from the MLS based method for Fredholm and Volterra integral equations presented by Mirzaei and Dehghan [24]. Here the Volterra type is compatible. Consider F distinct points T = {t, t 2,..., t F } in the time domain [, t F ], with the fill distance h T,[,tF ]. It is better t F be included because in various cases we need the solution at the final time. According to (2.3), replacing x by t and N by F, the univariate MLS approximation to equation (3.5) in respect to the time variable t, after imposing at t = t k for k F, is α F F ( t k ϕ l (t k )u(x, t l ) l= l= )( ϕ l (τ)dτ u(x, t l ) + ) κ κ u(x, t l) = κ tk f(x, τ)dτ + α u (x). The integrations over [, t k ] can be done by converting this interval to interval [, ] using the following linear transformation τ(t, θ) = tθ. (3.6) Therefore, if we set E k,l = ϕ l (t k ), k, l F, G k,l = tk ϕ l (τ)dτ = t k ϕ l (τ(t k, θ))dθ M = t k ϕ l (t k θ j )ω j, k, l F, j= u(x) = [u(x, t ),..., u(x, t F )], f(x) = [f (x),..., f F (x)], f k (x) = t k M j= u (x) = [u (x),..., u (x)] F, f(x, t k θ j )ω j, k F, (3.7) where {θ j, ω j } M j= is a M-point Gauss quadrature (or others) formula in [, ], then we have ( α EuT (x) G u T (x) + [ ] ) T κ u(x) = κ κ f T (x) + α ut (x). (3.8)
7 2 7 In addition from Dirichlet boundary condition (3.3), we have F ϕ l (t k )u(x, t l ) = u(x, t k ), k F, x Γ u. l= If we set u(x) = [u(x, t ),..., u(x, t F )], then Eu T (x) = u T (x), x Γ u. (3.9) Moreover, from the Neumann boundary condition (3.4), we have F l= ϕ l (t k ) u n (x, t l) = κ q(x, t k), k F, x Γ q and if we set q(x) = [q(x, t ),..., q(x, t F )], then E ut n (x) = κ qt (x), x Γ q. (3.) Equations (3.8)-(3.) construct a time-free system of boundary value problems which can be solved by a proper method. We are trying to apply the MLBIE method, presented by Zhu et.al [7], to this system of equations. Thus the MLS approximation in respect to the space variable x should be applied. According to the MLBIE process, first, the local weak form is written on local sub-domains Ω y σ Ω where Ω y σ = L y σ Γ y σ, and Γ y σ is a part of sub-domain s boundary which has intersection with global boundary Γ over which the Neumann boundary conditions are applied. The local sub-domains could be of any geometric shape and size. For simplicity they are taken to be B(y, σ) Ω. It is clear that for every internal node y, Γ y σ = and L y σ = Ω y σ. In MLBIE method the solution of the governing equation (3.8) in two dimensional space can be found in a second weak form via a modified fundamental solution (companion solution [7]) corresponds to the Poisson equation u (x, y) = 2π ln σ r u (x, y) + δ(x, y) =, where δ(x, y) is the Dirac delta function, r is the distance of the field and source points, i.e., r = x y 2 and σ is the radius of sub-domains. By the use of the Gauss divergence theorem to the weak form of Equation (3.8) and then applying the Neumann boundary condition (3.), we obtain an integral
8 8 DAVOUD MIRZAEI,,, MEHDI DEHGHAN representation E Ω y σ α ut (x)u (x, y)dω [ + G c(y)u T (y) + G = Ω y σ L y σ u T (x) u n (x, y)dγ + Γ y σ [ ] T κ u(x) u (x, y)dω Ω y κ σ [ κ f T (x) + α u (x)]u (x, y)dω + GE Γ y σ ] u T (x) u (x, y)dγ n κ qt (x)u (x, y)dγ, (3.) which is held for every y int(ω) Γ q, providing E is an invertible matrix. Note that the value c(y) depends on the location of y in the domain Ω or its boundary. If y int(ω) then c(y) =. For y Γ it depends on whether the source point y lies on the smooth boundary /2 or at corners β(y)/(2π) where β(y) is the interior angle at point y. The symbol denotes the Cauchy Principle Value (CPV) integral. The right hand side of Equation (3.) is a known vector of size F, say r(y). Note that neither Lagrange multiplier nor penalty parameters are introduced in the local weak form, because the Dirichlet boundary conditions are imposed directly using the MLS approximation. Hence, for y Γ u, Equation (3.9) is repeated replacing x by y. Now consider a quasi-uniform set X = {x, x 2,..., x N } in Ω with the fill distance h X,Ω. The set X should also satisfy the unisolvency condition to have a well-defined MLS approximation [22]. For simplicity suppose that x j Γ u for j N and x j Γ q int(ω) for N + j N. Thanks to the MLS approximations (2.3) and (2.4), the vectors u and its derivatives can be written in terms of values of u at centers X. To make difference between the MLS approximation in time and space, we use the notation φ j (x) instead of ϕ j (x) for shape functions in (2.3) and (2.4). Thus equation (3.9) can be approximated by and equation (3.) by ( N ) E φ j (y)u T (x j ) = u T (y), y Γ u. (3.2) j= ( N [ E j= Ω y σ ( N + G j= ] ) α φ j(x)u (x, y)dω u T (x j ) [ c(y)φ j (y) + L y σ φ j (x) u n Ω y κ κ φ j(x)u (x, y)dω σ ] u T (x j ) (x, y)dγ + Γ y σ ) = r(y). φ j (x) u (x, y)dγ n (3.3)
9 2 9 Imposing (3.2) at y = x i, i N and (3.3) at y = x i, N + i N and setting H (i, j) = φ j (x i ), i N, j N, H 2 (i, j) = Ω x i α(x) φ j(x)u (x, x i )dω, N + i N, j N, σ K 2 (i, j) = c(x i )φ j (x i ) + φ j (x) u n (x, x i)dγ + φ j (x) u n (x, x i)dγ Ω x i σ L x i σ Γ x i σ κ κ φ j(x)u (x, x i )dω, N + i N, j N, R (i, :) = u(x i ), N + i N, R 2 (i, :) = r(x i ), N + i N, U(j, :) = u(x j ), j N, give { E(H U) T = R T. E(H 2 U) T + G(K 2 U) T = R2 T Let H = [H T H T 2 ] N N, K = [ K T 2 ] N N and R = [R T R T 2 ] F N, thus we have EU T H + GU T K = R, (3.4) which is a generalized Sylvester equation. This equation should be solved by some linear algebra techniques that we will discuss in the following. Note that U(j, l) = u(x j, t l ) for j N and l F. When U is derived from equation (3.4), the solution at every point (x, t) Ω [, t F ] can be approximated by u(x, t) û(x, t) = = N φ j (x)u(x j, t) j= N j= l= F φ j (x)ϕ l (t)u(x j, t l ). If we set φ(x) = [φ (x),..., φ N (x)] and ϕ(t) = [ϕ (t),..., ϕ F (t)] then u(x, t) φ(x)uϕ T (t). (3.5) Note that, E and G are time and H and K are space matrices. The method provides the time and space matrices separately. Therefore other meshless methods in space can be easily replaced by the MLBIE, presented in this paper. The advantages of this method over the finite difference discretization in time can be expressed as bellow. In time discretization method we should use a difference equation of order O( t) like [u(t + t) u(t)]/ t instead of u/ t. Thus a small time step t requires to get accurate results. In the method of this paper the
10 DAVOUD MIRZAEI,,, MEHDI DEHGHAN convergence rate in the time domain is of order h m+ T,[,t F ]. Therefore, using not so many small time fill distance h T,[,tF ], accurate results are obtained. Laplace transform method is another technique to treat the time dependent problems. Several linear system of equations should be solved to get the solution at a sample time t. Unfortunately, the accuracy of the space method (here ML- BIE) will finally lose due to the poor accuracy of Laplace inversion techniques. In addition Laplace transforms of the body force and the boundary conditions should be provided. However, the method presented in this paper needs a linear algebra technique to solve the generalized Sylvester equation (3.4). There are several strategies. This equation can be converted to T U T + U T S = C, (3.6) with the time matrix T = G E, the space matrix S = KH and the right-hand side C = G RH, providing G and H are invertible matrices. Equation (3.6) is a standard Sylvester equation with unknown U. In Matlab the function X = lyap(t,s,-c) gives the solution U = X T. Numerical results show, matrix H is invertible and if is not included in {t,..., t F }, matrix G is also invertible. Note that there is no force to include zero. Equation (3.4) can be also solved directly using the Bartels-Stewart (BS) or Hessenberg-Schur (HS) algorithms presented in [25] based on QZ factorization, without any assumption on invertibility. Gardiner et.al [25] derived: if N F (if is not just change N and F ), the computational costs of BS and HS algorithms are BS : 33N F 3 + 3N 2 F + 3F 2 N, HS : 8.8N F 3 + ( p)N 2 F + 3F 2 N, where p.5. In this method, F is the number of meshless points in onedimensional time domain, and in many cases can be ignored in contrast to the number of meshless points in multi-dimensional space domain, i.e. N. These costs are roughly the same as the costs of solving many linear system of equations which is required in time discretization technique. Note that, solving a single generalized sylvester equation (3.4) gives the results at all time and space points, at once. This means that solving a generalized Sylvester equation is not a noticeable problem in sense of computational costs. Anyway, since in time difference method, LU decomposition and backward and forward substitutions can be used in some cases, the method of this paper is slightly more expensive.
11 2 4. Numerical results In numerical results, the basis { } (x z) β B = h β β m is used where h can be q X, h X,Ω or an average of them. And z is a fixed evaluation point such as a test point or a Gaussian point for integration in weak-form based techniques. Here β = (β, β 2 ) N 2 is a multi-index and β = β + β 2. If x = (x, x 2 ) then x β = x β xβ2 2. This choice of basis function, instead of B = {x β } β m, leads to a well-conditioned matrix A = P W P T in the MLS approximation (Equation (2.2)) and A is computed accurately, especially when the mesh-size is reduced. The effect of this variation on the conditioning of the MLS algorithm has been analytically investigated in [26]. We take m = 2 to span the trial spaces by the MLS approximation in both time and space domains. The local sub-domain s size is chosen as σ =.2. Numerical results in [27] show that it is a good choice. For integration over the local sub-domains in MLBIE (space) we use the Gauss integration formulas [2]. A fifteen points Gauss formula is employed for regular integrals and a seven points special Gauss formula for integrations with singular kernel ln r in [, ]. For the time integration (to generate the matrix G), if the time domain is large it is better to use some more accurate quadratures. Here we use a fifty points Gauss formula. If the time domain is small we can use some fewer points quadrature formulas. Note that the integration in the time domain is based on calling the one-dimensional MLS shape functions which has extremely low cost in comparison to the integration in multi-dimensional space domain. Therefore, the cost of integration in time can be really ignored. 4.. Problem. Consider a finite strip with a unidirectional variation of the thermal conductivity and diffusivity. The same exponential spatial variation for both parameters is taken κ(x) = κ exp(γx ), α(x) = α exp(γx ), (4.) with α =.7 4 m 2 s and κ = 7 Wm C. Three different exponential parameters γ =.2,.5 and. cm (2, 5 and m ) are selected in numerical calculation. On both opposite sides parallel to the x 2 -axis two different temperatures are prescribed. One side is kept to zero temperature and the other has the Heaviside step time variation i.e., u = T H(t) with T = C. On the lateral sides of the strip the heat flux vanishes. See Fig.. Such problem has been considered in [5] using MLBIE method with Laplace transform in time, and in [9] using an RBF based meshless collocation method with time difference approximation.
12 2 DAVOUD MIRZAEI,,, MEHDI DEHGHAN In numerical calculations, a square with a side-length a =.4 m and a regular node distribution is used for the MLS approximation in the space domain (Fig. ). Also a regular node distribution in the time domain t [, 6] with fill distance (mesh-size) h T,[,tF ] = 2 sec. is used. The values of unknown function u at arbitrary point (x, t) Ω [, t F ] is approximated using equation (3.5), after solving U from (3.4). The special case with an exponential parameter γ = corresponds to a homogeneous material. In such a case an analytical solution is available u(x, t) = T x a + 2 π n= T cos nπ n sin nπx a ( exp α n 2 π 2 ) t a 2, which can be used to check the accuracy of the present numerical method. Numerical results are computed at three locations along the x -axis with x /a =.25,.5 and.75. Results are depicted in Fig. 2. An excellent agreement between numerical and analytical results is obtained. It is known that the numerical results are rather inaccurate at very early time instants and at points close to the application of thermal shocks. Therefore in Fig. 3 we have compared the numerical and analytical solutions at very early time instants. Numerical results are obtained for x /a =.5, t F = 2 and h T,[,tF ] =.2 and then given for time-lengths 4 3 (using Equation (3.5)) by full-rectangles in this figure. The absolute error is approximately bounded by 2 4. The results are more accurate for points closer to the starting time. Besides, in Fig. 4 the numerical and analytical solutions at points very close to the application of thermal shocks are given and compared for t = sec. The maximum absolute error is Consequently, the method gives accurate results in these cases. The discussion above concerns heat conduction in homogeneous materials only since analytical solutions can be used for verification. To illustrate the application of the proposed algorithm, consider now the cases γ = 2, and 5 m, respectively. The variation of temperature with time for three γ-values and at position x =. and.2 m are presented in Figs. 5 and 6, respectively. Fig. 6 is in good agreement with Figure presented in [9], but some differences are evident between Fig. 5 of this paper and Figure 4 of [5]. This disagreement has also been reported in Table of [9]. Note that in this case the results of this paper are in agreement with Table of [9]. For final steady state an analytical solution can be obtained as u(x, t ) = T exp( γx ) ( exp( γa), u T x ) a, as γ. Analytical and numerical results computed at time t = 6 sec. corresponding to stationary or static loading conditions are presented in Fig. 7. Numerical results are in good agreement with analytical ones for the steady state temperatures.
13 2 3 As expected, it is found from Figs. 5, 6 and 7 that the temperature increases along with an increase in γ-values (or equivalently in thermal conductivity), and the temperature approaches a steady state when t > 2 sec Problem 2. In this example, an infinitely long and functionally graded thickwalled hollow cylinder is considered, where the following radii R = 8 2 m and R 2 = 2 m are selected. This problem has been considered in [5] using MLBIE method, where Laplace transform employed to eliminate the time variable. Heaviside step boundary condition is prescribed on the external surface of the hollow cylinder for the time-dependent thermal loading as a thermal shock with T = C. The inner surface is kept at zero temperature. Due to the symmetry in geometry and boundary conditions it is sufficient to analyze only a quarter of the cross section of the hollow cylinder. The boundary conditions and the total nodes lying in the domain and on the global boundary are depicted in Fig. 8. The thermal conductivity and diffusivity functions are chosen such as (4.) with same constants α and κ. For comparison purpose, a homogeneous material (γ = ) is first considered. In this homogeneous case, an analytical solution is known as: where u(r, t) = T ln(r/r ) ln(r 2 /R ) π J 2 (R β n )U (rβ n ) T J 2(R β n ) J 2(R 2β n ) exp ( α βnt 2 ), n= U (rβ n ) = J (rβ n )Y (R 2 β n ) J (R 2 β n )Y (rβ n ), and β n are the roots of the following transcendental equation J (r)y (rr 2 /R ) Y (r)j (rr 2 /R ) =, with J (r) and Y (r) being Bessel functions of first and second kinds and zeroth order. In numerical results we use time step h T,[,tF ] =.5 sec. in the time domain [, 6]. Fig. 9 depicts the numerical and exact solutions for various times at r =.9 m. The results show an excellent agreement with the analytical solution. Fig. shows the variation of the temperature with radial coordinate at four different time instants. Finally, we consider the functionally graded hollow cylinder with the thermal diffusivity and conductivity being graded in the radial direction. Numerical results for the time variation of the temperature are shown in Fig.. Similar to the results of a finite strip in the first example, the temperature level at interior points in the steady state increases with increasing γ-value. Here, again there are some differences between the results of this paper and the results presented in Figure 9 of [5] using MLBIE method and Laplace transform. It can be seen that, before getting the steady state solution, the results of this paper are larger than those obtained by MLBIE and Laplace transform. This disagreement also observed in the previous example. For γ = 5 m and m the temperature grows rapidly and gets
14 4 DAVOUD MIRZAEI,,, MEHDI DEHGHAN the steady-state in lower times. To deal with this argument we use very fine time steps h T,[,tF ] =. and h T,[,tF ] =.5, respectively. The results are presented in Fig. 2 and 3. As we can see, in the fist case for t > sec., and in the second case for t >.5, we get the steady state solutions. Finally Fig. 4 illustrates the variation of the temperature with the normalized radial coordinate r/r for different choices of γ at t = 2 sec. 5. Conclusion Heat conduction problem in functionally graded materials (FGMs) is investigated in this paper using the truly meshless local boundary integral equation (MLBIE) method. This method is based on the moving least squares (MLS) approximation in trial space. We avoid using Laplace transform and time discretization methods to eliminate the time variable. We employ again MLS as an approximation method in the time domain. Up to the best of our knowledge, it is a new technique to treat the time variable (at least) in the field of meshless methods. Numerical results which are provided for a finite strip and a hallow cylinder, show the accuracy and efficiency of this technique. Although the process has been presented for the MLBIE method, it can be extend to other MLS based meshless methods, similarly. Acknowledgements The first author was partially supported by the Center of Excellence for Mathematics, University of Isfahan. References [] A. Sutradhar, G.H. Paulino, and L.G. Gray. Transinet heat conduction in homogeneous and nonhomogeneous materials by the Laplace transfom Galerkin boundary element method. Engineering Analysis with Boundary Elements, 26:9 32, 22. [2] A. Sutradhar, G.H. Paulino, and L.G. Gray. On hypersingular surface integrals in the symmetric Galerkin boundary element method: application to heat conduction in exponentially graded materials. International Journal for Numerical Methods in Engineering, 62:22 57, 25. [3] Y. Miyamoto, W.A. Kaysser, B.H. Rabin, A. Kawasaki, and R.G. Ford. Functionally Graded Materials; Design, Processing and Applications. Kluwer, Dordrecht, 999. [4] J. Sladek, V. Sladek, Ch. Hellmich, and J. Eberhardsteiner. Heat conduction analysis of 3-D axisymmetric and anisotropic FGM bodies by meshless local Petrov-Galerkin method. Computational Mechanics, 39: , 27. [5] J. Sladek, V. Sladek, and Ch. Zhang. Transient heat conduction analysis in functionally graded materials by the meshless local boundary integral equation method. Computational Materials Science, 28:494 54, 23. [6] J. Sladek, V. Sladek, and Ch. Zhang. A local BIEM for analysis of transient heat conduction with nonlinear source terms in FGMs. Engineering Analysis with Boundary Elements, 28:, 24.
15 2 5 [7] V. Sladek, J. Sladek, M. Tanaka, and Ch. Zhang. Transient heat conduction in anisotropic and functionally graded media by local integral equations. Engineering Analysis with Boundary Elements, 29:47 65, 25. [8] J. Sladek, V. Sladek, and S.N. Atluri. Meshless local Petrov-Galerkin method for heat conduction problem in an anisotropic medium. CMES Computer Modeling in Engineering and Sciences, 6:39 38, 24. [9] H. Wang, Q.-H. Qin, and Y.-L. Kang. A meshless model for transient heat conduction in functionally graded materials. Computational Mechanics, 38:5 6, 26. [] S. Hamaza-Cherif, A. Houmat, and A. Hadjou. Transinet heat conduction in functionally graded materials. International Journal of Computational Methods, 4:63 69, 27. [] M. Dehghan and D. Mirzaei. The meshless local Petrov-Galerkin (MLPG) method for the generalized two-dimensional non-linear Schrödinger equation. Engineering Analysis with Boundary Elements, 32: , 28. [2] D. Mirzaei and M. Dehghan. Implementation of meshless LBIE method to the 2D non-linear SG problem. International Journal for Numerical Methods in Engineering, 79: , 29. [3] D. Mirzaei and M. Dehghan. Meshless local Petrov-Galerkin (MLPG) approximation to the two dimensional sine-gordon equation. Journal of Computational and Applied Mathematics, 233: , 2. [4] S. Abbasbandy and A. Shirzadi. A meshless method for two-dimensional diffusion equation with an integral condition. Engineering Analysis with Boundary Elements, 34:3 37, 2. [5] S. Abbasbandy and A. Shirzadi. MLPG method for two-dimensional diffusion equation with neumann s and non-classical boundary conditions. Applied Numerical Mathematics, 6:79 8, 2. [6] T. Belytschko, Y. Krongauz, D.J. Organ, M. Fleming, and P. Krysl. Meshless methods: an overview and recent developments. Computer Methods in Applied Mechanics and Engineering, special issue, 39:3 47, 996. [7] T. Zhu, J.D. Zhang, and S.N. Atluri. A local boundary integral equation (LBIE) method in computational mechanics, and a meshless discretization approach. Computational Mechanics, 2: , 998. [8] T. Zhu, J.D. Zhang, and S.N. Atluri. A meshless local boundary integral equation (LBIE) method for solving nonlinear problems. Computational Mechanics, 22:74 86, 998. [9] P. Lancaster and K. Salkauskas. Surfaces generated by moving least squares methods. Mathematics of Computation, 37:4 58, 98. [2] D. Shepard. A two-dimensional interpolation function for irregularly-spaced data. In Proceedings of the 23th National Conference ACM, pages , 968. [2] H. Wendland. Local polynomial reproduction and moving least squares approximation. 2:285 3, 2. [22] H. Wendland. Scattered Data Approximation. Cambridge University Press, 25. [23] D. Levin. The approximation power of moving least-squares. Mathematics of Computation, 67:57 53, 998. [24] D. Mirzaei and M. Dehghan. A meshless based method for solution of integral equations. Applied Numerical Mathematics, 6: , 2. [25] G.D. Gardiner, A.J. Laub, J.J. Amato, and C.B. Moler. Solution of the Sylvester matrix equation AXB T + CXD T = E. AMS Transactions on Mathematical Software, 8(2):223 23, 992.
16 6 DAVOUD MIRZAEI,,, MEHDI DEHGHAN [26] D. Mirzaei, R. Schaback, and M. Dehghan. On generalized moving least squares and diffuse derivatives. IMA Journal of Numerical Analysis, 2. doi:.98/imanum/drr3. [27] M. Dehghan and D. Mirzaei. Numerical solution to the unsteady two-dimensional Schrödinger equation using meshless local boundary integral equation method. International Journal for Numerical Methods in Engineering, 76:5 52, 28.
17 2 7 x 2 (cm) 4 q = 3 u = 2 u = T.H(t) q = x (cm) Figure. Boundary conditions and node distribution for a finite strip. Temperature u (C ) x/a =.25, Numerical. " Analytical x/a =.5, Numerical " Analytical -. x/a =.75, Numerical " Analytical Time t (sec.) Figure 2. Time variation of the temperature in a finite strip at three positions with γ =.
18 8 DAVOUD MIRZAEI,,, MEHDI DEHGHAN Temperature u (C ) 5 x Analytical Numerical Time t (sec.) Figure 3. Accuracy of method for early time instants at position x /a =.5 of the functionally graded finite strip. Temprature u (C ) Analytical Numerical x /a (m) Figure 4. Accuracy of method for points close to the thermal shock at time t = sec. of the functionally graded finite strip.
19 2 9.4 Temprature u (C ).3.2. = m - = 2 m - = 5 m Time t (sec.) Figure 5. Time variation of the temperature at position x /a =.25 of the functionally graded finite strip Temperature u (C ) = m -. = 2 m - = 5 m Time t (sec.) Figure 6. Time variation of the temperature at position x /a =.5 of the functionally graded finite strip.
20 2 DAVOUD MIRZAEI,,, MEHDI DEHGHAN Temperature u (C ) = m -, Num. " Anal. = 2 m -, Num. " Anal. = 5 m -, Num. " Anal. = m -, Num. " Anal x /a (m) Figure 7. Distribution of temperature along x -axis for a functionally graded finite square strip under steady-state loading conditions. q = 8 u = T.H(t) x 2 (cm) 6 4 R 2 u = 2 R q = x (cm) Figure 8. Boundary conditions and node distribution for a hollow cylinder.
21 2 2.5 Temperature u (C ) Analytical solution Numerical solution Time t (sec.) Figure 9. Exact and numerical solution at r =.9 m in hallow cylinder with homogeneous material properties (γ = )..8 Temperature u (C ).6.4 t = 2 sec..2 t = 4 sec. t = 6 sec. t = sec r (m) Figure. Temperature variation with radial coordinate r in hallow cylinder at different time levels with homogeneous material properties (γ = ).
22 22 DAVOUD MIRZAEI,,, MEHDI DEHGHAN.6 Temperature u (C ) = m - = m - = 2 m Time t (sec.) Figure. Time variation of the temperature at r =.9 m in a functionally graded hollow cylinder..7.6 Temperature u (C ) = 5 m Time t (sec.) Figure 2. Time variation of the temperature at r =.9 m for γ = 5 m for a functionally graded hollow cylinder.
23 Temperature u (C ) = m Time t (sec.) Figure 3. Time variation of the temperature at r =.9 m for γ = m for a functionally graded hollow cylinder..2 Temperature u (C ) = m - = 2 m - = 5 m - = m r (m) Figure 4. Temperature variation with radial coordinate r at t = 2 sec. for a functionally graded hollow cylinder
Solving Heat Conduction Problems by the Direct Meshless Local Petrov-Galerkin (DMLPG) method
Solving Heat Conduction Problems by the Direct Meshless Local Petrov-Galerkin (DMLPG) method Davoud Mirzaei, Robert Schaback Department of Mathematics, University of Isfahan, 8745-63 Isfahan, Iran. Institut
More informationON GENERALIZED MOVING LEAST SQUARES AND DIFFUSE DERIVATIVES
ON GENERALIZED MOVING LEAST SQUARES AND DIFFUSE DERIVATIVES DAVOUD MIRZAEI, ROBERT SCHABACK,, MEHDI DEHGHAN Abstract. The Moving Least Squares (MLS) method provides an approximation û of a function u based
More informationPoint interpolation method based on local residual formulation using radial basis functions
Structural Engineering and Mechanics, Vol. 14, No. 6 (2002) 713-732 713 Point interpolation method based on local residual formulation using radial basis functions G.R. Liu, L. Yan, J.G. Wang and Y.T.
More informationNumerical solution of Maxwell equations using local weak form meshless techniques
Journal of mathematics and computer science 13 2014), 168-185 Numerical solution of Maxwell equations using local weak form meshless techniques S. Sarabadan 1, M. Shahrezaee 1, J.A. Rad 2, K. Parand 2,*
More informationAnalysis of Steady State Heat Conduction Problem Using EFGM
International Journal of Engineering and Management Research, Vol.-2, Issue-6, December 2012 ISSN No.: 2250-0758 Pages: 40-47 www.ijemr.net Analysis of Steady State Heat Conduction Problem Using EFGM Manpreet
More informationA truly meshless Galerkin method based on a moving least squares quadrature
A truly meshless Galerkin method based on a moving least squares quadrature Marc Duflot, Hung Nguyen-Dang Abstract A new body integration technique is presented and applied to the evaluation of the stiffness
More informationDavoud Mirzaei s CV. University of Shahid Chamran, Ahvaz, Iran, September 2001 July 2005, Bachelor of Science, Applied Mathematics.
Davoud Mirzaei s CV Davoud Mirzaei, Associate Professor of Applied Mathematics, (Numerical Analysis) Department of Mathematics, University of Isfahan, 81745-163 Isfahan, Iran. Tell: 0098-31-37934635 Email:
More informationSolving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels
Solving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels Y.C. Hon and R. Schaback April 9, Abstract This paper solves the Laplace equation u = on domains Ω R 3 by meshless collocation
More informationSolving Elastic Problems with Local Boundary Integral Equations (LBIE) and Radial Basis Functions (RBF) Cells
Copyright 2010 Tech Science Press CMES, vol.57, no.2, pp.109-135, 2010 Solving Elastic Problems with Local Boundary Integral Equations (LBIE) and Radial Basis Functions (RBF) Cells E. J. Sellountos 1,
More informationSolving the One Dimensional Advection Diffusion Equation Using Mixed Discrete Least Squares Meshless Method
Proceedings of the International Conference on Civil, Structural and Transportation Engineering Ottawa, Ontario, Canada, May 4 5, 2015 Paper No. 292 Solving the One Dimensional Advection Diffusion Equation
More informationIn this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots,
Chapter 8 Elliptic PDEs In this chapter we study elliptical PDEs. That is, PDEs of the form 2 u = lots, where lots means lower-order terms (u x, u y,..., u, f). Here are some ways to think about the physical
More informationAnalyzing Plane-plate Bending with EFGM
Journal of Mathematics Research March, 2009 Analyzing Plane-plate Bending with EFGM Yajing Zhang, Maohui Xia & Yun Zhai Department of Science, YanShan University Qinhuangdao 066004, China E-mail: zhang
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 informationDIRECT ERROR BOUNDS FOR SYMMETRIC RBF COLLOCATION
Meshless Methods in Science and Engineering - An International Conference Porto, 22 DIRECT ERROR BOUNDS FOR SYMMETRIC RBF COLLOCATION Robert Schaback Institut für Numerische und Angewandte Mathematik (NAM)
More informationNUMERICAL MLPG ANALYSIS OF PIEZOELECTRIC SENSOR IN STRUCTURES
DOI: 10.2478/sjce-2014-0009 NUMERICAL MLPG ANALYSIS OF PIEZOELECTRIC SENSOR IN STRUCTURES Peter STAŇÁK 1*, Ján SLÁDEK 1, Vladimír SLÁDEK 1, Slavomír KRAHULEC 1 Abstract The paper deals with a numerical
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 14: The Power Function and Native Space Error Estimates Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu
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 informationPSEUDO ELASTIC ANALYSIS OF MATERIAL NON-LINEAR PROBLEMS USING ELEMENT FREE GALERKIN METHOD
Journal of the Chinese Institute of Engineers, Vol. 27, No. 4, pp. 505-516 (2004) 505 PSEUDO ELASTIC ANALYSIS OF MATERIAL NON-LINEAR PROBLEMS USING ELEMENT FREE GALERKIN METHOD Raju Sethuraman* and Cherku
More informationThe purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and engineering.
Lecture 16 Applications of Conformal Mapping MATH-GA 451.001 Complex Variables The purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and
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 informationMeshless Local Petrov-Galerkin (MLPG) Methods in
Meshless Local Petrov-Galerkin MLPG) Methods in Quantum Mechanics *Williams L Nicomedes Renato C Mesquita *Ferno J S Moreira * Federal University of Minas Gerais Dept of Electronic Engineering Belo Horizonte
More informationApplication of meshless EFG method in fluid flow problems
Sādhan ā Vol. 29, Part 3, June 2004, pp. 285 296. Printed in India Application of meshless EFG method in fluid flow problems I V SINGH Mechanical Engineering Group, Birla Institute of Technology and Science,
More informationExistence Theory: Green s Functions
Chapter 5 Existence Theory: Green s Functions In this chapter we describe a method for constructing a Green s Function The method outlined is formal (not rigorous) When we find a solution to a PDE by constructing
More informationJournal of Theoretical and Applied Mechanics, Sofia, 2012, vol. 42, No. 1, pp
Journal of Theoretical and Applied Mechanics, Sofia, 2012, vol. 42, No. 1, pp. 23 40 SOLID MECHANICS CALCULATION OF STRESS INTENSITY FACTORS BASED ON FORCE-DISPLACEMENT CURVE USING ELEMENT FREE GALERKIN
More informationLeast-squares collocation meshless method
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2001; 51:1089 1100 Least-squares collocation meshless method Xiong Zhang ;, Xiao-HuLiu, Kang-ZuSong and Ming-Wan Lu
More informationStability of Kernel Based Interpolation
Stability of Kernel Based Interpolation Stefano De Marchi Department of Computer Science, University of Verona (Italy) Robert Schaback Institut für Numerische und Angewandte Mathematik, University of Göttingen
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationRadial point interpolation based finite difference method for mechanics problems
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 26; 68:728 754 Published online 11 April 26 in Wiley InterScience (www.interscience.wiley.com). DOI: 1.12/nme.1733
More informationNumerical solution of surface PDEs with Radial Basis Functions
Numerical solution of surface PDEs with Radial Basis Functions Andriy Sokolov Institut für Angewandte Mathematik (LS3) TU Dortmund andriy.sokolov@math.tu-dortmund.de TU Dortmund June 1, 2017 Outline 1
More informationA Posteriori Error Bounds for Meshless Methods
A Posteriori Error Bounds for Meshless Methods Abstract R. Schaback, Göttingen 1 We show how to provide safe a posteriori error bounds for numerical solutions of well-posed operator equations using kernel
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 informationNumerical solution of nonlinear sine-gordon equation with local RBF-based finite difference collocation method
Numerical solution of nonlinear sine-gordon equation with local RBF-based finite difference collocation method Y. Azari Keywords: Local RBF-based finite difference (LRBF-FD), Global RBF collocation, sine-gordon
More informationIndex. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2
Index advection equation, 29 in three dimensions, 446 advection-diffusion equation, 31 aluminum, 200 angle between two vectors, 58 area integral, 439 automatic step control, 119 back substitution, 604
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 informationMeshless Local Petrov-Galerkin Mixed Collocation Method for Solving Cauchy Inverse Problems of Steady-State Heat Transfer
Copyright 2014 Tech Science Press CMES, vol.97, no.6, pp.509-533, 2014 Meshless Local Petrov-Galerkin Mixed Collocation Method for Solving Cauchy Inverse Problems of Steady-State Heat Transfer Tao Zhang
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 34: Improving the Condition Number of the Interpolation Matrix Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu
More informationSpotlight on Laplace s Equation
16 Spotlight on Laplace s Equation Reference: Sections 1.1,1.2, and 1.5. Laplace s equation is the undriven, linear, second-order PDE 2 u = (1) We defined diffusivity on page 587. where 2 is the Laplacian
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 informationThe Nottingham eprints service makes this work by researchers of the University of Nottingham available open access under the following conditions.
Power, H. and Caruso, N. and Portapila, M. (2017) A note on the use of the Companion Solution (Dirichlet Green's function) on meshless boundary element methods. Engineering Analysis with Boundary Elements,
More informationCLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE
CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that
More informationNumerical and Analytical Approach of Thermo-Mechanical Stresses in FGM Beams
Numerical and Analytical Approach of Thermo-Mechanical Stresses in FGM Beams Fatemeh Farhatnia, Gholam-Ali Sharifi and Saeid Rasouli Abstract In this paper, thermo-mechanical stress distribution has been
More informationMeshless equilibrium on line method (MELM) for linear elasticity
Structural Engineering and Mechanics, Vol. 35, No. 4 (2010) 000-000 1 Meshless equilibrium on line method (MELM) for linear elasticity A. Sadeghirad*, S. Mohammadi and I. Mahmoudzadeh Kani School of Civil
More informationD. Shepard, Shepard functions, late 1960s (application, surface modelling)
Chapter 1 Introduction 1.1 History and Outline Originally, the motivation for the basic meshfree approximation methods (radial basis functions and moving least squares methods) came from applications in
More information13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)
13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs) A prototypical problem we will discuss in detail is the 1D diffusion equation u t = Du xx < x < l, t > finite-length rod u(x,
More informationLecture 13: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay. Poisson s and Laplace s Equations
Poisson s and Laplace s Equations Lecture 13: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay We will spend some time in looking at the mathematical foundations of electrostatics.
More informationError estimates for moving least square approximations
Applied Numerical Mathematics 37 (2001) 397 416 Error estimates for moving least square approximations María G. Armentano, Ricardo G. Durán 1 Departamento de Matemática, Facultad de Ciencias Exactas y
More informationA Meshless Radial Basis Function Method for Fluid Flow with Heat Transfer
Copyright c 2008 ICCES ICCES, vol.6, no.1, pp.13-18 A Meshless Radial Basis Function Method for Fluid Flow with Heat Transfer K agamani Devi 1,D.W.Pepper 2 Summary Over the past few years, efforts have
More informationAnalysis of free vibration problems with the Element-Free Galerkin method
Analysis of free vibration problems with the Element-Free Galerkin method Carlos M. Tiago, Vitor M. A. Leitão Universidade Técnica de Lisboa, Instituto Superior Técnico Av. Rovisco Pais, 49- Lisboa, Portugal
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationSSPH basis functions for meshless methods, and comparison of solutions with strong and weak formulations
Comput Mech 8 4:57 545 DOI.7/s466-7-9- ORIGINAL PAPER SSPH basis functions for meshless methods, and comparison of solutions with strong and weak formulations R. C. Batra G. M. Zhang Received: 7 February
More informationGeneralized Finite Element Methods for Three Dimensional Structural Mechanics Problems. C. A. Duarte. I. Babuška and J. T. Oden
Generalized Finite Element Methods for Three Dimensional Structural Mechanics Problems C. A. Duarte COMCO, Inc., 7800 Shoal Creek Blvd. Suite 290E Austin, Texas, 78757, USA I. Babuška and J. T. Oden TICAM,
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 informationIntroduction to Finite Element Method
Introduction to Finite Element Method Dr. Rakesh K Kapania Aerospace and Ocean Engineering Department Virginia Polytechnic Institute and State University, Blacksburg, VA AOE 524, Vehicle Structures Summer,
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationScattered Data Interpolation with Polynomial Precision and Conditionally Positive Definite Functions
Chapter 3 Scattered Data Interpolation with Polynomial Precision and Conditionally Positive Definite Functions 3.1 Scattered Data Interpolation with Polynomial Precision Sometimes the assumption on the
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 informationConsistency Estimates for gfd Methods and Selection of Sets of Influence
Consistency Estimates for gfd Methods and Selection of Sets of Influence Oleg Davydov University of Giessen, Germany Localized Kernel-Based Meshless Methods for PDEs ICERM / Brown University 7 11 August
More informationFundamental Solutions and Green s functions. Simulation Methods in Acoustics
Fundamental Solutions and Green s functions Simulation Methods in Acoustics Definitions Fundamental solution The solution F (x, x 0 ) of the linear PDE L {F (x, x 0 )} = δ(x x 0 ) x R d Is called the fundamental
More informationQuintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationNumerical Analysis of Electromagnetic Fields
Pei-bai Zhou Numerical Analysis of Electromagnetic Fields With 157 Figures Springer-Verlag Berlin Heidelberg New York London Paris Tokyo Hong Kong Barcelona Budapest Contents Part 1 Universal Concepts
More informationFINITE ELEMENT METHOD: APPROXIMATE SOLUTIONS
FINITE ELEMENT METHOD: APPROXIMATE SOLUTIONS I Introduction: Most engineering problems which are expressed in a differential form can only be solved in an approximate manner due to their complexity. The
More informationAPPLIED PARTIM DIFFERENTIAL EQUATIONS with Fourier Series and Boundary Value Problems
APPLIED PARTIM DIFFERENTIAL EQUATIONS with Fourier Series and Boundary Value Problems Fourth Edition Richard Haberman Department of Mathematics Southern Methodist University PEARSON Prentice Hall PEARSON
More informationA Continuation Approach to a Quadratic Matrix Equation
A Continuation Approach to a Quadratic Matrix Equation Nils Wagner nwagner@mecha.uni-stuttgart.de Institut A für Mechanik, Universität Stuttgart GAMM Workshop Applied and Numerical Linear Algebra September
More informationPositive Solution to a Generalized Lyapunov Equation via a Coupled Fixed Point Theorem in a Metric Space Endowed With a Partial Order
Filomat 29:8 2015, 1831 1837 DOI 10.2298/FIL1508831B Published by Faculty of Sciences Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Positive Solution to a Generalized
More informationA Meshless Local Petrov-Galerkin Method for Solving the Bending Problem of a Thin Plate
Copyright cfl 2002 Tech Science Press CMES, vol.3, no.1, pp.53-63, 2002 A Meshless Local Petrov-Galerkin Method for Solving the Bending Problem of a Thin Plate Shuyao Long 1, S. N. Atluri 2 Abstract: Meshless
More informationLecture notes: Introduction to Partial Differential Equations
Lecture notes: Introduction to Partial Differential Equations Sergei V. Shabanov Department of Mathematics, University of Florida, Gainesville, FL 32611 USA CHAPTER 1 Classification of Partial Differential
More informationAn Introduction to Partial Differential Equations in the Undergraduate Curriculum
An Introduction to Partial Differential Equations in the Undergraduate Curriculum John C. Polking LECTURE 12 Heat Transfer in the Ball 12.1. Outline of Lecture The problem The problem for radial temperatures
More informationMeshless Method for Crack Analysis in Functionally Graded Materials with Enriched Radial Base Functions
Copyright c 2008 Tech Science Press CMES, vol.30, no.3, pp.133-147, 2008 Meshless Method for Crack Analysis in Functionally Graded Materials with Enriched Radial Base Functions P.H. Wen 1, M.H. Aliabadi
More informationStability constants for kernel-based interpolation processes
Dipartimento di Informatica Università degli Studi di Verona Rapporto di ricerca Research report 59 Stability constants for kernel-based interpolation processes Stefano De Marchi Robert Schaback Dipartimento
More informationCIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen
CIV-E16 Engineering Computation and Simulation Examination, December 12, 217 / Niiranen This examination consists of 3 problems rated by the standard scale 1...6. Problem 1 Let us consider a long and tall
More informationConnection to Laplacian in spherical coordinates (Chapter 13)
Connection to Laplacian in spherical coordinates (Chapter 13) We might often encounter the Laplace equation and spherical coordinates might be the most convenient 2 u(r, θ, φ) = 0 We already saw in Chapter
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationPARTIAL DIFFERENTIAL EQUATIONS. Lecturer: D.M.A. Stuart MT 2007
PARTIAL DIFFERENTIAL EQUATIONS Lecturer: D.M.A. Stuart MT 2007 In addition to the sets of lecture notes written by previous lecturers ([1, 2]) the books [4, 7] are very good for the PDE topics in the course.
More informationThe first order quasi-linear PDEs
Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.
More informationMixed discrete least square meshless method for solution of quadratic partial di erential equations
Scientia Iranica A (014) 1(3), 49{504 Sharif University of Technology Scientia Iranica Transactions A: Civil Engineering www.scientiairanica.com Mixed discrete least square meshless method for solution
More informationRegularized Meshless Method for Solving Acoustic Eigenproblem with Multiply-Connected Domain
Copyright c 2006 Tech Science Press CMES, vol16, no1, pp27-39, 2006 Regularized Meshless Method for Solving Acoustic Eigenproblem with Multiply-Connected Domain KH Chen 1,JTChen 2 and JH Kao 3 Abstract:
More informationPartial Differential Equations
M3M3 Partial Differential Equations Solutions to problem sheet 3/4 1* (i) Show that the second order linear differential operators L and M, defined in some domain Ω R n, and given by Mφ = Lφ = j=1 j=1
More informationLinear Regression and Its Applications
Linear Regression and Its Applications Predrag Radivojac October 13, 2014 Given a data set D = {(x i, y i )} n the objective is to learn the relationship between features and the target. We usually start
More informationMAT389 Fall 2016, Problem Set 4
MAT389 Fall 2016, Problem Set 4 Harmonic conjugates 4.1 Check that each of the functions u(x, y) below is harmonic at every (x, y) R 2, and find the unique harmonic conjugate, v(x, y), satisfying v(0,
More informationA Meshless Simulations for 2D Nonlinear Reaction-diffusion Brusselator System
Copright 3 Tech Science Press CMES, vol.95, no.4, pp.59-8, 3 A Meshless Simulations for D Nonlinear Reaction-diffusion Brusselator Sstem Ahmad Shirzadi, Vladimir Sladek, Jan Sladek 3 Abstract: This paper
More informationOutline. 1 Boundary Value Problems. 2 Numerical Methods for BVPs. Boundary Value Problems Numerical Methods for BVPs
Boundary Value Problems Numerical Methods for BVPs Outline Boundary Value Problems 2 Numerical Methods for BVPs Michael T. Heath Scientific Computing 2 / 45 Boundary Value Problems Numerical Methods for
More informationMathematics Qualifying Exam Study Material
Mathematics Qualifying Exam Study Material The candidate is expected to have a thorough understanding of engineering mathematics topics. These topics are listed below for clarification. Not all instructors
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 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 informationHEAT CONDUCTION USING GREEN S FUNCTIONS
HEAT CONDUCTION USING GREEN S FUNCTIONS Preface to the first edition Preface to the second edition Author Biographies Nomenclature TABLE OF CONTENTS FOR SECOND EDITION December 2009 Page viii x xii xiii
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 informationApproximation by Conditionally Positive Definite Functions with Finitely Many Centers
Approximation by Conditionally Positive Definite Functions with Finitely Many Centers Jungho Yoon Abstract. The theory of interpolation by using conditionally positive definite function provides optimal
More informationMatrix Theory. A.Holst, V.Ufnarovski
Matrix Theory AHolst, VUfnarovski 55 HINTS AND ANSWERS 9 55 Hints and answers There are two different approaches In the first one write A as a block of rows and note that in B = E ij A all rows different
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 information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
More informationA mesh-free minimum length method for 2-D problems
Comput. Mech. (2006) 38: 533 550 DOI 10.1007/s00466-005-0003-z ORIGINAL PAPER G. R. Liu K. Y. Dai X. Han Y. Li A mesh-free minimum length method for 2-D problems Received: 7 March 2005 / Accepted: 18 August
More informationThe DMLPG Meshless Technique for. Poisson Problems
Applied Mathematical Sciences, Vol. 8, 2014, no. 164, 8233-8250 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.49739 The DMLPG Meshless Technique for Poisson Problems Flavio Sartoretto
More informationA gradient recovery method based on an oblique projection and boundary modification
ANZIAM J. 58 (CTAC2016) pp.c34 C45, 2017 C34 A gradient recovery method based on an oblique projection and boundary modification M. Ilyas 1 B. P. Lamichhane 2 M. H. Meylan 3 (Received 24 January 2017;
More informationPreliminary Examination in Numerical Analysis
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify
More informationMeshless analysis of three-dimensional steady-state heat conduction problems
Cin. Pys. B ol. 19, No. 9 (010) 09001 Mesless analysis of tree-dimensional steady-state eat conduction problems Ceng Rong-Jun( 程荣军 ) a) and Ge Hong-Xia( 葛红霞 ) b) a) Ningbo Institute of ecnology, Zejiang
More informationDiscrete Projection Methods for Integral Equations
SUB Gttttingen 7 208 427 244 98 A 5141 Discrete Projection Methods for Integral Equations M.A. Golberg & C.S. Chen TM Computational Mechanics Publications Southampton UK and Boston USA Contents Sources
More informationDepartment of Structural, Faculty of Civil Engineering, Architecture and Urban Design, State University of Campinas, Brazil
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm A SIMPLIFIED FORMULATION FOR STRESS AND TRACTION BOUNDARY IN- TEGRAL EQUATIONS USING THE
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationError Analyses of Stabilized Pseudo-Derivative MLS Methods
Department of Mathematical Sciences Revision: April 11, 2015 Error Analyses of Stabilized Pseudo-Derivative MLS Methods Don French, M. Osorio and J. Clack Moving least square (MLS) methods are a popular
More information