THE METHOD OF FUNDAMENTAL SOLUTION FOR SOLVING MULTIDIMENSIONAL INVERSE HEAT CONDUCTION PROBLEMS
|
|
- Laurel Scott
- 6 years ago
- Views:
Transcription
1 THE METHOD OF FUNDAMENTAL SOLUTION FOR SOLVING MULTIDIMENSIONAL INVERSE HEAT CONDUCTION PROBLEMS Y.C. HON AND T.WEI Abstract. We propose in this paper an effective meshless and integration-free method for the numerical solution of multidimensional inverse heat conduction problems. Due to the use of fundamental solutions as basis functions, the method leads to a global approximation scheme in both the spatial and time domains. To tackle the ill-conditioning problem of the resultant linear system of equations, we apply the Tikhonov regularization method by using the L-curve criterion for the regularization parameter to obtain a stable approximation to the solution. The effectiveness of the algorithm is illustrated by several numerical two- and three-dimensional examples.. Introduction A standard inverse heat conduction problem (IHCP) is to compute the unknown temperature and heat flux at an unreachable boundary from scattered temperature measurements at reachable interior or boundary of the domain. In solving direct heat conduction problems, the errors induced from boundary or interior measurements are reduced due to the diffusive characteristic of heat conduction process. These errors, however, are extrapolated and amplified due to the extremely illposedness of the inverse heat conduction problems. In other words, a small error in the measurement can induce enormous error in computing the unknown solution at the unreachable boundary. Several techniques have been proposed for solving a one-dimensional IHCP [4, 6, 32, 33, 34, 4]. Among the methods proposed for higher dimensional IHCP, boundary element [7, 3], finite difference [3, 27] and finite element [24, 4] have Key words and phrases. Inverse heat conduction problem, fundamental solution method, Tikhonov regularization. The work described in this paper was partially supported by a grant from CityU (Project No. 746).
2 2 Y.C. HON AND T.WEI been widely adopted for problems in two-dimension. Besides, the sequential function specification method [4, 7] and mollification method [36] have also been used in solving the IHCP. There is, however, still a need on numerical scheme for multidimensional IHCP. The traditional mesh-dependent finite difference and finite element methods require dense meshes, and hence tedious computational time, to give a reasonable approximation to the solution of IHCP and suffers from numerical instability problem. The use of boundary element method (BEM) reduces the computational time and storage requirement but the problem of numerical stability still persists. Since there is no need on domain discretization in the BEM, the location of interior points, where the temperature data are collected, can be chosen in a quite arbitrary way [7]. In this paper a new meshless computational method is proposed to approximate the solution of a multidimensional IHCP under arbitrary geometry. The proposed method uses the fundamental solution of the corresponding heat equation to generate a basis for approximating the solution of the problem. Comparing with the mesh-dependent methods like FEM and BEM, the proposed method does not require any domain or boundary discretization. This meshless advantage makes the method feasible to solve high dimensional IHCPs under arbitrary geometry. Furthermore, unlike BEM which requires an interpolation of the unknown function on the boundary for complicated boundary integrations, the proposed method gives a global numerical solution on the entire time-spatial domain. The required surface temperature and heat flux on the unreachable boundary can easily be computed without any extra quadratures. To tackle the numerical instability in the resultant ill-conditioned linear system, we use Tikhonov regularization technique and L-curve method [4, 5, 6, 7, 3] to obtain a stable numerical solution for the multi-dimensional IHCP problem. In the recent rapid development of meshless computational schemes, the radial basis function (RBF) has been successfully developed as an efficient meshless scheme for solving various kinds of partial differential equations (PDEs). Most of the problems are still confined to direct problems [8, 9, 2, 23, 25]. In [2], Hon and Wu first applied the meshless RBF to solve a Cauchy problem for Laplace equation,
3 METHOD OF FUNDAMENTAL SOLUTION FOR IHCP 3 which is a severely ill-posed problem. In [22], Hon and Wei combined the meshless RBF with the method of fundamental solution (MFS) to successfully solve an onedimensional inverse heat conduction problem. This approach can be interpreted as an extension of the MFS for treating elliptic problems, for examples, the Laplace equation [5, 35], the biharmonic equations [26], elastostatics problems [38], and wave scattering problems [28, 29]. The MFS was also applied to solve nonhomogeneous linear and nonlinear Poisson equations [, 2, 3,, 37]. More details of the MFS can be found in the review papers of Fairweather and Karageorghis [] and Golberg and Chen [2]. It is noted here that most of these research works focus on wellposed problems in which the Dirichlet or Neumann data are given on the whole boundary. In the IHCP problem, the boundary conditions are usually complicated and incomplete. We expect that the combination of the meshless RBF and MFS will extend its application to solve higher dimensional inverse problems under irregular geometry. In [9], Frankel et al. gave an unified treatment for the IHCP by using the Chebyshev polynomials as basis functions. The use of fundamental solution as the basis functions in this paper provides a global approximation to the solution in both the spatial and time variables. Since IHCP problem is severely ill-posed and its solution is extremely sensitive to any perturbation of given data, enormous error will be accumulated from using any finite difference scheme for discretizing the time variable [8]. The proposed method has a definite advantage over most of the existing numerical methods in solving these kinds of time-dependent inverse problems. The use of Laplace transform for the time variable will reduce this time discretization error but result in solving a Cauchy problem for inhomogeneous modified Helmholtz equation with an extra parameter in the resulting equation [8]. For a slightly larger parameter, the Laplace transform method fails to give a stable and accurate approximation. In fact, it is even more difficult to solve a Cauchy problem for modified Helmholtz than solving the original IHCP. For large scale problems, an efficient numerical method for treating the ill-conditioning discrete problem is still needed for improvement. The recently developed iterative algorithms based on Lanczos bidiagonalization will be a consideration for future work.
4 4 Y.C. HON AND T.WEI 2. Methodology Let Ω be a simply connected domain in R d, d = 2, 3 and Γ, Γ 2, Γ 3 be three parts of the boundary Ω. Suppose that Γ Γ 2 Γ 3 = Ω, Γ or Γ 2 can be empty set. The IHCP to be investigated in this paper is to determine the temperature and heat flux on boundary Γ 3 from given Dirichlet data on Γ, Neumann data on Γ 2 and scattered measurement data at some interior points. Consider the following heat equation: (2.) u t (x, t) = a 2 u(x, t), x Ω, t (, t max ), under the initial condition (2.2) u(x, ) = ϕ(x), x Ω, and the boundary conditions (2.3) u(x, t) = f(x, t), x Γ, t (, t max ], and (2.4) ν (x, t) = g(x, t), x Γ 2, t (, t max ], where ν is the outer unit normal with respect to Γ 2. Let {x i } M i= of exact temperature u(x i, t (k) i Ω be a set of locations with noisy measured data h (k) i ) = h (k) i, i =, 2,, M, k =, 2,, I i where t (k) i (, t max ] are discrete times. The absolute error between the noisy measurement and exact data is assumed to be bounded, i.e. h (k) i h (k) i δ for all measurement points at all measured times. Here, the constant δ is called the noisy level of input data. The IHCP is now formulated as: Reconstruct u Γ3 and ν Γ 3 from (2.) (2.4) (k) and the scattered noisy measurements h i, i =, 2,, M, k =, 2,, I i. Remark 2.. Most of the existing numerical methods for solving the IHCP problem require a continuous time measurement u(x i, t) = h i (t), which is not realistic in real life. The fundamental solution of (2.) is given by (2.5) F (x, t) = (4πa 2 t) d 2 e x 2 4a 2 t H(t),
5 METHOD OF FUNDAMENTAL SOLUTION FOR IHCP 5 where H(t) is the Heaviside function. Assuming that T > t max is a constant, the following function (2.6) φ(x, t) = F (x, t + T ), is a general solution of (2.) in the solution domain Ω [, t max ]. We denote the measurement points to be {(x j, t j )} m j=, m = M i= I i so that a point at the same location but with different time is treated as two distinct points. The corresponding measured noisy data and exact data are denoted by h j h j. The collocation points are then chosen as {(x j, t j )} m+n j=m+ on the initial region Ω {}, {(x j, t j )} m+n+p j=m+n+ on surface Γ (, t max ] and {(x j, t j )} m+n+p+q j=m+n+p+ on surface Γ 2 (, t max ]. Here, n, p, q denote the total number of collocation points for the initial condition (2.2), Dirichlet boundary condition (2.3) and Neumann boundary condition (2.4) respectively. and The only requirement on the collocation points are pairwisely distinct in the (d + )-dimensional space (x, t). Following the idea of MFS, an approximation ũ to the solution of the IHCP under the conditions (2.2)-(2.4) with the noisy measurements h j can be expressed by the following linear combination: (2.7) ũ(x, t) = n+m+p+q j= λ j φ(x x j, t t j ), where φ(x, t) is given by (2.6) and λ j are unknown coefficients to be determined. For this choice of basis function φ, the approximated solution ũ automatically satisfies the original heat equation (2.). Using the initial condition (2.2) and collocating at the boundary conditions (2.3) and (2.4), we then obtain the following system of linear equations for the unknown coefficients λ j : (2.8) A λ = b, where (2.9) A = φ(x i x j, t i t j ) φ ν (x k x j, t k t j )
6 6 Y.C. HON AND T.WEI and (2.) b = hi ϕ(x i, t i ) f(x i, t i ) g(x k, t k ) where i =, 2,, (m + n + p), k = (m + n + p + ),, (m + n + p + q), j =, 2,, (n + m + p + q) respectively. The solvability of the system (2.8) depends on the non-singularity of the matrix A, which is still an open research problem. It is not surprise that the resultant matrix A is extremely ill-conditioned due to the ill-posed nature of the IHCP. The following sections shows that the use of regularization technique can produce a stable and accurate solution for the unknown parameters λ of the matrix equation (2.8), and hence the solution for the IHCP is obtained. 3. Regularization Techniques The ill-conditionedness of the coefficient matrix A indicates that the numerical result is sensitive to the noise of right hand side b and the number of collocation points. In fact, the condition number of the matrix A increases dramatically with respect to the total number of collocation points. The singular value decomposition (SVD) usually works well for direct problem [39] but still fails to provide a stable and accurate solution to the system (2.8). A number of regularization methods have been developed for solving this kind of ill-conditioning problem [4, 5, 7]. In our computation we apply the Tikhonov regularization technique [42] based on SVD with L-curve criterion for choosing a good regularization parameter to solve the matrix equation (2.8). Denote N = n + m + p + q. The SVD of the N N matrix A is a decomposition of the form: N (3.) A = W ΣV = w i σ i v i where W = (w, w 2,, w N ) and V = (v, v 2,, v N ) satisfying W W = V V = I N. Here, the superscript represents the transpose of a matrix. It is known that i=
7 METHOD OF FUNDAMENTAL SOLUTION FOR IHCP 7 Σ = diag(σ, σ 2,, σ N ) has non-negative diagonal elements satisfying: (3.2) σ σ N. The values σ i are called the singular values of A and the vectors w i and v i are called the left and right singular vectors of A respectively. For the matrix A arising from the discretization of MFS, the singular values decay rapidly to zero and the ratio between the largest and the smallest nonzero singular values is often huge. Thus the linear system (2.8) is a discrete ill-posed problem in the sense defined in Hansen s paper [4]. Based on the singular value decomposition, it is easy to know that the solution for the linear equations (2.8) is given by (3.3) λ = N i= w i b σ i v i. The corresponding solution with exact data is (3.4) λ = N i= w i b σ i v i where b is calculated by (2.) from the exact data h j. these two solutions is then given by (3.5) λ λ = N i= w i e σ i v i, The difference between where e = b b, e δ. The terms in the difference with small values of σ i will be large since the noisy level is generally greater than the smallest nonzero singular value. proposed. This is the reason why the following Tikhonov regularization method is The Tikhonov regularized solution λ α for equation (2.8) is defined to be the solution to the following least square problem: (3.6) min { A λ b λ e 2 + α 2 λ 2 }, where denotes the usual Euclidean norm and α is called the regularization parameter. The Tikhonov regularized solution based on SVD can the be expressed as: (3.7) λα = N i= f i w i b σ i v i,
8 8 Y.C. HON AND T.WEI where f i = σ 2 i /(σ2 i + α2 ) are called the filter factors, i =, 2,, N. The difference between λ α and λ is then given by (3.8) λα λ = The norm of this error vector is (3.9) λ α λ = N i= ( N i= (f i ) w i b σ i v i + N i= f i w i e σ i v i. ( (f i ) w i b ) 2 ( w i + f e ) ) /2 2 i. σ i σ i Hence, the approximated solution ũ α with noisy measurement for the IHCP is given by (3.) ũ α (x, t) = N ( λ α ) j φ(x x j, t t j ), j= where ( λ α ) j is the i-th entry of vector λ α. Let u N = N i= λ jφ(x x j, t t j ) be an approximated solution with exact data given in the IHCP. An error estimation is then given as follow: (3.) ũ α u N L 2 (D) λ α λ N φ(x x j, t t j ) L 2 (D), where D = Ω (, t max ) and the value λ α λ is given by formula (3.9). j= The determination of a suitable value for the regularization parameter α is crucial and is still under intensive research [42]. In this paper, we use the L-curve criterion to choose a good regularization parameter. This method was firstly developed by Lawson and Hanson [3] and more recently by Hansen and O Leary [5]. L-curve method is sketched in the following. Define a curve L by: The (3.2) L = { (log( λ α 2 ), log( A λ α b 2 )), α > }. The curve is known as L-curve and a suitable regularization parameter α is the one near the corner of the L-curve [4, 7]. In the numerical computation, a point with maximum curvature is seek to be the corner as given in [7] by (3.3) K(α) = 2 ηρ η α 2 ηρ + 2αηρ + α 4 ηη (α 2 η 2 + ρ 2 ) 3/2
9 METHOD OF FUNDAMENTAL SOLUTION FOR IHCP 9 where and η = λ α 2 = N w (f i b i ) 2, ρ = A λ α σ b N 2 = (( f i )w i b) 2, i i= η = 4 α N ( f i )fi 2 ( w i b ) 2. σ i i= In our computation, the Matlab code developed by Hansen [6] was used to obtain the optimal choice of regularization parameter α for solving the discrete illconditioned system (2.8). The corresponding approximated solution for the problem (2.) (2.4) with noisy measurement data is then given by (3.4) ũ α (x, t) = i= N ( λ α ) j φ(x x j, t t j ). j= The temperature and heat flux at surface Γ 3 can also be calculated respectively. 4. Numerical Examples For numerical verification, we assume that the heat conduction coefficient a = and t max = for all the following examples. In the cases when the input data contain noises, we use the function rand given in Matlab to generate the noisy data hi = h i + δrand(i) where h i is the exact data and rand(i) is a random number in [, ]. The magnitude δ indicates the noisy level of the measurement data. To test the accuracy of the approximated solution, we compute the Root Mean Square Error (RMSE) by (4.) E(u) = N t ((ũ α ) i u i ) N 2, t i= at sufficiently total N t number of testing points in the domain Γ 3 [, ] where(ũ α ) i and u i are respectively the approximated and exact temperature at a test point. The RMSE for the heat flux E( ν ) is also similarly defined. Consider the following three examples: Two-dimensional IHCPs: Example : The exact solution of (2.) is given by (4.2) u(x, x 2, t) = 2t + 2 (x2 + x 2 2).
10 Y.C. HON AND T.WEI Example 2: The exact solution of (2.) is given by (4.3) u(x, x 2, t) = e 4t (cos(2x ) + cos(2x 2 )). Both examples are computed under the following three different domains and boundary conditions: Case : Let Ω = { (x, x 2 ) < x <, < x 2 < }, Γ = { (x, x 2 ) x =, < x 2 < }, Γ 2 = { (x, x 2 ) < x <, x 2 = }, Γ 3 = Ω \ {Γ Γ 2 }. Case 2: Let Ω = { (x, x 2 ) x 2 + x 2 2 < }, Γ =, Γ 2 = { (x, x 2 ) x 2 + x 2 2 =, x >, }. Γ 3 = Ω \ {Γ Γ 2 }. Case 3: Let Ω be the same as Case, Γ = Γ 2 =, Γ 3 = Ω \ {Γ Γ 2 }. Locations of the internal measurements and collocation points over Ω for the three cases are shown in Figure. Case. Case 2. Case 3..9 Γ 2.9 Γ x 2.6 Γ.5 3 Γ Ω Γ x Ω Γ 3 Γ 2.6 Γ 3 Γ.5 3 Ω Γ Figure. Distribution of measurement points and collocation points. Here, star represents collocation point matching Dirichlet data, square represents collocation point matching Neumann data, dot represents collocation point matching initial data, and circle represents point with sensor for internal measurement. The boundary data f(x, t), g(x, t) and initial temperature ϕ(x, t) are obtained from the given exact solutions. The values of temperature at measurement points
11 METHOD OF FUNDAMENTAL SOLUTION FOR IHCP will be used respectively by the exact and noisy data. Numerical results are obtained by taking the constant T =.8 for Example and T =.6 for Example 2 in all different cases. The total numbers of various collocation points and measurement points are n = 36, m =, p = 55, q = 5 in Case for Examples -2. In this situation the total number of points and testing points are N = 24 and N t = 882. For Case 2, the total numbers are n = 96, m = 85, p =, q = 85, N = 266, N t = 693 and for Case 3, n = 36, m = 3, p =, q =, N = 66, N t = 764 respectively. Numerical results by using only SVD for the solution (3.4) and u N are presented in Table. The computed RMSEs from the experiments show that even for exact input data the direct method cannot produce an acceptable solution to this kind of ill-conditioned linear systems. In fact, the condition numbers in Case and Case 2 lie between which are too large to compute an accurate solution without the use of any regularization. For Case 3, the condition number in Example and in Example 2 are much smaller and hence the RMSEs in Case 3 looks better. The numerical results given in Table indicate that the IHCP is an severely ill-posed problem. The discretization by using MFS also leads to a highly ill-conditioned discrete problem. The use of the regularization technique as shown in the following numerical results will give a stable and much more accurate approximation to the solution of the IHCP. Table 2 gives the RMSEs on the temperature and heat flux in domain Γ 3 [, ] with exact data by using Tikhonov regularization method with L-curve choice α for the regularization parameter. It can be observed from Table 2 that the RMSEs have much been reduced compared to Table. It is remarked here that the L-curve method works well in searching the crucial regularization parameter but can further be improved if a scaling factor is used. Numerical results for the three examples with noisy data (noisy level δ =. in all cases) are shown in Table 3. It can be observed that the regularization method with L-curve technique provides an acceptable approximation to the solution of the IHCP whilst the direct method completely fails. The numerical result for the Example with Case 2 looks bad but can be much improved if a scaling factor. is multiplied to the L-curve regularized parameter α. In this case the
12 2 Y.C. HON AND T.WEI RMSEs for E(u) and E( ν ) are.248 and.6 respectively. The problem to choose an optimal regularization parameter is still an open question to researchers. It is noted here that the setting for Case 3 comes from a real-life problem produced by a steel company. In fact, the solution for the IHCP may not unique but our computational results demonstrate that the proposed method is flexible enough to give a reasonable approximation to the solution of the IHCP under insufficient information. Table. RMSEs in domain Γ 3 [, ] with exact data. No regularization technique Case Case 2 Case 3 Example Example 2 E(u)=5.922 E(u)=.499 E( ν )=2.255 E(u)= E( ν )=.54 E( ν )=.2325 E(u)=7.966 E( ν )=.7736 E(u)=.2 E(u)=. E( ν )=.543 E( ν )=.63 Table 2. RMSEs in domain Γ 3 [, ] with exact data by using Tikhonov regularization technique with L-curve parameter α Case Case 2 Case 3 Example Example 2 α = 2.744e-2 E(u)= 2.787e-4 α =3.4429e- E(u)=9.35e-5 E( ν )=6.692e-4 E( ν )=3.54e-4 α =4.55e-3 E(u)=. α = 8.663e-6 E(u)=. E( ν )=.29 E( ν )=.33 α =.37e-8 E(u)= 8.992e-4 α =2.3e-2 E(u)=4.82e-4 E( ν )=.3 E( ν )=.25 The relationship between the RMSEs and the value of constant T for Examples -2 under Case 2 with noisy data (δ =.) is displayed in Figure 2. Here, the regularization parameter α is obtained again by using the L-curve method. The
13 METHOD OF FUNDAMENTAL SOLUTION FOR IHCP 3 Table 3. RMSEs in domain Γ 3 [, ] with noisy data (δ =.) by using Tikhonov regularization with L-curve parameter α Case Case 2 Case 3 Example Example 2 α =2.77e-3 E(u)=.729 α =3.668e-4 E(u)=.46 E( ν )=.348 E( ν )=.45 α =.28 E(u)=.2778 α =.23e-5 E(u)=.4 E( ν )=.465 E( ν )=.38 α =2.7527e-5 E(u)=.22 α = 4.982e-6 E(u)=.2 E( ν )=.494 E( ν )=.578 numerical results indicate that the choice of the parameter T plays a minor but strange role for obtaining an accurate approximation. For Example, the RMSEs decrease with respect to the increasing value of T but decrease in Example 2. This interesting behavior is one of the focuses of our future research works..7 Example for case 2. σ=..4 Example 2 for case 2. σ= RMS(u) and RMS(du).4.3 RMS(u) and RMS(du) parameter T parameter T Figure 2. RMSEs of temperature and heat flux on Γ 3 [, ] with respect to parameter T To further extend the application of the proposed method, we investigate the following sample problem given in [7]: Example 3: Let Ω be defined as in Case. Γ =, Γ 2 = { (x, x 2 ) x =, < x 2 < } { (x, x 2 ) x 2 =, < x < }, Γ 3 = Ω \ {Γ Γ 2 }.
14 4 Y.C. HON AND T.WEI The locations of measurement points are shown in Figure 3. The temperature distribution is given by U(x, x 2, t), for t.5, u(x, x 2, t) = U(x, x 2, t) 2U(x, x 2, t.5), for.5 < t., U(x, x 2, t) 2U(x, x 2, t.5) + U(x, x 2, t.), for t >., where U(x, x 2, t) =.5t 2 + t(.5x 2 x + x 2 2 2x 2 + ) 4 j= (jπ) (.5cos(jπx 4 ) + cos(jπx 2 ))( e j2 π 2t ). The exact temperature data at the sensor locations are given by u(x, x 2, t) and the noisy data are randomly generated as before. In this computation, we take n = 2, m = 89, p =, q = 89, N = 499, N t = 882. The plots of temperature and heat flux with respect to time at points (, ), (, ), (, ) are displayed in Figure 4 in which T =.2 and δ =.. As shown in Figure 4, although the basis functions generated by MFS are sufficient smooth over the solution domain and the solution to Example 3 does not have a continuous second order time derivative, the computed temperature match the exact data excellently. The computed heat flux looks less accurate but is also comparable to the results given in paper [7]. It is also noted that the measurement points in this example are far away from the unspecified boundary whereas these points are close to the boundary in [7]. Finally, we consider the following three-dimensional IHCP: Example 4: Let Ω = { (x, x 2, x 3 ) < x i <, i =, 2, 3}, Γ = { (x, x 2, x 3 ) < x <, < x 2 <, x 3 = }, Γ 2 = { (x, x 2, x 3 ) < x <, < x 2 <, x 3 = }, Γ 3 = Ω \ {Γ Γ 2 }. The locations of measurement points and collocation points in the domain Ω are shown in Figure 5. In this computation, we take n = 245, m = 25, p = 8, q = 8, N = 855, N t = The exact solution for Example 4 is given by (4.4) u(x, x 2, x 3, t) = e ( 4t) (cos(2x ) + cos(2x 2 ) + cos(2x 3 )).
15 METHOD OF FUNDAMENTAL SOLUTION FOR IHCP 5.9 Example 3. Γ x Γ3 Ω Γ Γ x Figure 3. Distribution of measurement and collocation points. Square represents point with Neumann data, dot represents collocation point for initial temperature, and circle represents point with sensor. In the computation, the value of the parameter T is 2.3 and the noisy level is set to be δ =.. The errors between the exact solution and the approximated solution for the temperature and heat flux on boundary Γ 3 shown in Figure 6 and Figure 7 respectively. ν at time t = are Note that the L 2 norm of u and over Γ 3 [, ] are about.62 and.52 respectively. Their relative errors are approximately double the values given in Figures 6 and Figure 7. These small relative errors show that the proposed scheme is effective for solving the threedimensional inverse heat conduction problem of which very little numerical result has been given so far. 5. Conclusions The universal approach for solving time-dependent problems involves a timemarching procedure, i.e., advancing each time step and solving the remained problem in the spatial domain. The approach proposed in this paper gives a global approximated solution in both time and spatial domain. Numerical results indicate that the method of fundamental solution MFS combined with the regularization technique provides an efficient and accurate approximation to this highly ill-posed inverse heat conduction. The use of L-curve criterion for a suitable regularization
16 6 Y.C. HON AND T.WEI u(,,t).6.4 dudy(,,t) t t u(,,t) dudy(,,t) t t u(,,t).4.3 dudx(,,t) t t Figure 4. Plots of temperature and heat flux versus time at points (, ), (, ), (, ). Solid line represents the exact value and dotted-line represents the computed value parameter stabilizes the resultant ill-conditioned system but still needed to be further investigated for optimal convergence. The proposed approach is potentially valuable to solve the inverse problems in multidimensional space but for large-scale problems, other iterative algorithms based on Lanczos bidiagonalization may be
17 METHOD OF FUNDAMENTAL SOLUTION FOR IHCP 7 Example 4. Γ.8.6 x 3.4 Ω.2.5 x 2.5 Γ 2 x Figure 5. Locations of measurement and collocation points. Star represents measurement point with Dirichlet data, square represents measurement point with Neumann data, and circle represents point with sensor. more suitable. This will be our future work. Acknowledgment We wish to thank the referees for their constructive comments which helped us much in improving the presentation of this paper.
18 8 Y.C. HON AND T.WEI x 3 x 3 Error of u(x,,x 3,) x 3 x.5 Error of u(,x 2,x 3,) x 3 x 2.5 x 3 x 3 Error of u(x,,x 3,) x 3 x.5 Error of u(,x 2,x 3,) x 3 x 2.5 Figure 6. Surface plots of errors to temperature on boundary Γ 3 References [] Balakrishnan, K. and Ramachandran, P.A., A particular solution Trefftz method for nonlinear Poisson problems in heat and mass transfer, Journal of computational physics, 5(999), [2] Balakrishnan, K. and Ramachandran, P.A., The method of fundamental solutions for linear diffusion-reaction equations, Mathematical and Computer Modelling, 3(2), [3] Balakrishnan, K. and Ramachandran, P.A., Osculatory Interpolation in the Method of Fundamental Solution for Nonliear Possion Problems, Journal of Computational Physics, 72(2), 8. [4] Beck, J.V., Blackwell, B. and Clair, Ch.R.St., Inverse heat conduction, Ill-posed Problems, Wiley-Interscience Publication, New York, 985. [5] Bogomonlny, A., Fundamental solutions method for elliptic boundary value problems, SIAM Journal on Numerical Analysis, 22 (985), [6] Chantasiriwan, S., Comparison of three sequential function specification algorithms for the inverse heat conduction problem, International Communications in Heat and Mass Transfer, 26(999), no., 5-24.
19 METHOD OF FUNDAMENTAL SOLUTION FOR IHCP 9 x 3 Error of dudy(x,,x 3,)...5 x 3 x.5 Error of dudx(,x 2,x 3,) x 3 x 2.5 Error of dudy(x,,x 3,)...5 x 3 x.5 Error of dudx(,x 2,x 3,)..2.5 x 3 x 2.5 Figure 7. Surface plots of errors to heat flux on boundary Γ 3 [7] Chantasiriwan, S., An algorithm for solving multidimensional inverse heat conduction problem, International Journal of Heat and Mass Transfer, 44 (2), [8] Cho, H.A., Golberg, M.A., Muleshkov, A.S. and Li. X, Trefftz methods for time dependent partial differential equations (preprint). [9] Frankel, J.I. and Keyhani, M., A global time treatment for inverse heat conduction problems, Journal of Heat Transfer, 9 (997), [] Fairweather, G. and Karageorghis, A., The method of fundamental solutions for elliptic boundary value problems, Advances in Computational Mathematics, 9(998), [] Golberg, M.A., The method of fundamental solutions for Poisson s equation, Engineering Analysis with Boundary Elements, 6(995), [2] Golberg, M.A. and Chen, C.S., The method of fundamental solutions for potential, Helmholtz and diffusion problems. In: M. A. Golberg editor. Boundary integral methods-numerical and mathematical aspects, Sounthampton: Computational Mechanic Publications, 998, [3] Guo, L. and Murio, D.A, A mollified space-marching finite-difference algorithm for the twodimensional inverse heat conduction problem with slab symmetry, Inverse problem, 7 (99),
20 2 Y.C. HON AND T.WEI [4] Hansen, P.C., Analysis of discrete ill-posed problems by means of the L-curve, SIAM Review, 34 (992), [5] Hansen, P.C. and O Leary, D.P., The use of the L-curve in the regularization of discrete ill-posed problems. SIAM Journal of Scientific Computing, 4 (993), [6] Hansen, P.C., Regularization tools: a Matlab package for analysis and solution of discrete ill-posed problems, Numerical Algorithms, 6 (994), 35. [7] Hansen, P.C., The L-curve and its use in the numerical treatment of inverse problems, Computational Inverse Problems in Electrocardiology, 4, edited by Johnston, P., Advances in Computational Bioengineering Series, WIT Press, Southampton, 2. [8] Hon, Y.C. and Mao, X.Z., A radial basis function method for solving options pricing models, The Journal of Financial Engineering, 8(999), [9] Hon, Y.C., Cheung, K.F., Mao, X.Z. and Kansa, E.D., Multiquadric solution for shallow water equations, ASCE Journal of Hydraulic Engineering, 25(999), [2] Hon, Y.C., Lu, M.W., Xue, W.M. and Zhou, X, A new formulation and computation of the triphasic model for mechano-electrochemical mixtures, Computational Mechanics, 24 (999), [2] Hon, Y. C. and Wu, Z.M., A numerical computation for inverse boundary determination problem, Engineering Analysis with Boundary Elements, 24(2), [22] Hon, Y.C. and Wei, T., A Meshless Computational Method for Solving Inverse Heat Conduction Problem, International Series on Advances in Boundary Elements, 3(22), [23] Hon, Y.C. and Chen, W., Boundary knot method for 2D and 3D Helmholtz and convectiondiffusion problems with complicated geometry, International Journal for Numerical Methods in Engineering, 56(23), [24] Hsu, T.R., Sun, N.S., Chen, G.G. and Gong, Z.L., Finite element formulation for twodimensional inverse heat conduction analysis, Transactions of the ASME, Journal of Heat Transfer, 4 (992), [25] Kansa, E.J. and Hon, Y.C., Circumventing the ill-conditioning problem with multiquadric radial basis functions: applications to elliptic partial differential equations, Computers and Mathematics with Applications, 39(2), [26] Karageorghis A. and Fairweather G., The method of fundamental solutions for the numerical solution of the biharmonic equation, J. Comput. Phys., 69(987) [27] Khalidy, N.A., A general space marching algorithm for the solution of two-dimensional boundary inverse heat conduction problems, Numerical Heat Transfer, Part B, 34 (998), [28] Kondapalli, P.S., Shippy, D.J. and Fairweather, G., Analysis of acoustic scattering in fluids and solids by the method of fundamental solutions, J. Acoust. Soc. Amer., 9(992), [29] Kondapalli, P.S., Shippy, D.J. and Fairweather, G., The method of fundamental solutions for transmission and scattering of elastic waves, Comput. Methods Appl. Mech. Engrg., 96(992),
21 METHOD OF FUNDAMENTAL SOLUTION FOR IHCP 2 [3] Kurpisz, K. and Nowak, A.J., BEM approach to inverse heat conduction problems, Engineering Analysis with Boundary Elements, (992), [3] Lawson, C.L. and Hanson, R.J., Solving Least Squares Problems, SIAM, Philadelphia, PA, 995. First published by Prentice-Hall, 974. [32] Lesnic, D., Elliott, L. and Ingham, D.B., Application of the boundary element method to inverse heat conduction problems, International Communications in Heat and Mass Transfer 39 (996), No. 7, [33] Lesnic, D. and Elliott, L, The decomposition approach to inverse heat conduction, Journal of Mathematical Analysis and Application, 232(999), [34] Jonas, P and Louis, A.k., Approximate inverse for a one-dimensional inverse heat conduction problem, Inverse Problems, 6 (2), [35] Mathon, R. and Johnston, R.L., The approximatical solution of elliptic boundary-value problems by fundamental solutions, SIAM J. Numer. Anal., 4(977), [36] Murio, D.A., The mollification method and the numerical solution of ill-posed problems, A Wiley-Interscience Publication, 993. [37] Partridge, P.W. and Sensale, B., The method of fundamental solutions with dual reciprocity for diffusion and diffusion-convection using subdomains, Engineering Analysis with Boundary Elements, 24(22), [38] Poullikkas, A., Karageorghis, A. and Georgiou, G., The method of fundamental solutions for three-dimensional elastostatics problems, Computers and Structures, 8(22) [39] Ramachandran, P.A., Method of fundamental solutions: singular value decomposition analysis, Communications in numerical methods in engineering, 8(22), [4] Reinhardt, H.J., A numerical method for the solution of two-dimensional inverse heat conduction problems, International Journal for Numerical Methods in Enginerring, 32 (99), [4] Shen, S.Y., A numerical study of inverse heat conduction problems, Computers and Mathematics with Applications, 38(999), [42] Tikhonov, A. N. and Arsenin, V.Y., On the solution of ill-posed problems, John Wiley and Sons, New York, 977. Department of Mathematics, City University of Hong Kong, Tat Chee Avenue, Kowloon, HongKong address: maychon@cityu.edu.hk Department of Mathematics, City University of Hong Kong & Department of Mathematics, Lanzhou University, Lanzhou, Gansu Province 73, China address: t.wei@student.cityu.edu.hk
A Meshless Method for the Laplace and Biharmonic Equations Subjected to Noisy Boundary Data
Copyright c 4 Tech Science Press CMES, vol.6, no.3, pp.3-61, 4 A Meshless Method for the Laplace and Biharmonic Equations Subjected to Noisy Boundary Data B. Jin 1, Abstract: In this paper, we propose
More informationDetermination of a source term and boundary heat flux in an inverse heat equation
ISSN 1746-7659, England, UK Journal of Information and Computing Science Vol. 8, No., 013, pp. 3-114 Determination of a source term and boundary heat flux in an inverse heat equation A.M. Shahrezaee 1
More information1 Comparison of Kansa s method versus Method of Fundamental Solution (MFS) and Dual Reciprocity Method of Fundamental Solution (MFS- DRM)
1 Comparison of Kansa s method versus Method of Fundamental Solution (MFS) and Dual Reciprocity Method of Fundamental Solution (MFS- DRM) 1.1 Introduction In this work, performances of two most widely
More informationBOUNDARY PARTICLE METHOD FOR INVERSE CAUCHY PROBLEMS OF INHOMOGENEOUS HELMHOLTZ EQUATIONS
Journal of Marine Science and Technology, Vol. 7, No. 3, pp. 57-63 (9) 57 BOUNDARY PARTICLE METHOD FOR INVERSE CAUCHY PROBLEMS OF INHOMOGENEOUS HELMHOLTZ EQUATIONS Wen Chen* and Zhuojia Fu* Key words:
More informationNEW INVESTIGATIONS INTO THE MFS FOR 2D INVERSE LAPLACE PROBLEMS. Received June 2015; revised October 2015
International Journal of Innovative Computing, Information and Control ICIC International c 2016 ISSN 1349-4198 Volume 12, Number 1, February 2016 pp. 55 66 NEW INVESTIGATIONS INTO THE MFS FOR 2D INVERSE
More informationAn efficient approximation method for the nonhomogeneous backward heat conduction problems
International Journal o Engineering and Applied Sciences IJEAS) ISSN: 394-3661 Volume-5 Issue-3 March 018 An eicient appromation method or the nonhomogeneous bacward heat conduction problems Jin Wen Xiuen
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 informationSOLVING INHOMOGENEOUS PROBLEMS BY SINGULAR BOUNDARY METHOD
8 Journal of Marine Science and Technology Vol. o. pp. 8-4 03) DOI: 0.69/JMST-0-0704- SOLVIG IHOMOGEEOUS PROBLEMS BY SIGULAR BOUDARY METHOD Xing Wei Wen Chen and Zhuo-Jia Fu Key words: singular boundary
More informationA novel method for estimating the distribution of convective heat flux in ducts: Gaussian Filtered Singular Value Decomposition
A novel method for estimating the distribution of convective heat flux in ducts: Gaussian Filtered Singular Value Decomposition F Bozzoli 1,2, L Cattani 1,2, A Mocerino 1, S Rainieri 1,2, F S V Bazán 3
More informationA Truly Boundary-Only Meshfree Method Applied to Kirchhoff Plate Bending Problems
Advances in Applied Mathematics and Mechanics June 2009 Adv. Appl. Math. Mech., Vol. 1, No. 3, pp. 341-352 (2009) A Truly Boundary-Only Meshfree Method Applied to Kirchhoff Plate Bending Problems Zhuojia
More informationA method of fundamental solutions for two-dimensional heat conduction
A method of fundamental solutions for two-dimensional heat conduction B. Tomas Johansson, D. Lesnic, Thomas Henry Reeve To cite this version: B. Tomas Johansson, D. Lesnic, Thomas Henry Reeve. A method
More informationREGULARIZATION PARAMETER SELECTION IN DISCRETE ILL POSED PROBLEMS THE USE OF THE U CURVE
Int. J. Appl. Math. Comput. Sci., 007, Vol. 17, No., 157 164 DOI: 10.478/v10006-007-0014-3 REGULARIZATION PARAMETER SELECTION IN DISCRETE ILL POSED PROBLEMS THE USE OF THE U CURVE DOROTA KRAWCZYK-STAŃDO,
More informationInverse Source Identification for Poisson Equation
Inverse Source Identification for Poisson Equation Leevan Ling, Y. C. Hon, and M. Yamamoto March 1, 2005 Abstract A numerical method for identifying the unknown point sources for a two-dimensional Poisson
More informationJournal Articles (Refereed. Total 127) 1. Hon Y.C., Kong Y., and Xu C., Identification of gifted children from the nature of giftedness, draft paper.
Journal Articles (Refereed. Total 127) 1. Hon Y.C., Kong Y., and Xu C., Identification of gifted children from the nature of giftedness, draft paper. 2. Hon Y.C. and Schaback R., Kernel techniques for
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 informationThe Method of Fundamental Solutions applied to the numerical calculation of eigenfrequencies and eigenmodes for 3D simply connected domains
ECCOMAS Thematic Conference on Meshless Methods 2005 C34.1 The Method of Fundamental Solutions applied to the numerical calculation of eigenfrequencies and eigenmodes for 3D simply connected domains Carlos
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 informationCompactly supported radial basis functions for solving certain high order partial differential equations in 3D
Compactly supported radial basis functions for solving certain high order partial differential equations in 3D Wen Li a, Ming Li a,, C.S. Chen a,b, Xiaofeng Liu a a College of Mathematics, Taiyuan University
More informationRECONSTRUCTION OF NUMERICAL DERIVATIVES FROM SCATTERED NOISY DATA. 1. Introduction
RECONSTRUCTION OF NUMERICAL DERIVATIVES FROM SCATTERED NOISY DATA T. WEI, Y. C. HON, AND Y. B. WANG Abstract. Based on the thin plate spline approximation theory, we propose in this paper an efficient
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 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 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 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 informationTransient heat conduction by different versions of the Method of Fundamental Solutions a comparison study
Computer Assisted Mechanics and Engineering Sciences, 17: 75 88, 2010. Copyright c 2010 by Institute of Fundamental Technological Research, Polish Academy of Sciences Transient heat conduction by different
More informationRecent developments in the dual receiprocity method using compactly supported radial basis functions
Recent developments in the dual receiprocity method using compactly supported radial basis functions C.S. Chen, M.A. Golberg and R.A. Schaback 3 Department of Mathematical Sciences University of Nevada,
More informationThis is a repository copy of The method of fundamental solutions for an inverse boundary value problem in static thermo-elasticity.
This is a repository copy of The method of fundamental solutions for an inverse boundary value problem in static thermo-elasticity. White Rose Research Online URL for this paper: http://eprints.whiterose.ac.uk/816/
More informationDiscrete ill posed problems
Discrete ill posed problems Gérard MEURANT October, 2008 1 Introduction to ill posed problems 2 Tikhonov regularization 3 The L curve criterion 4 Generalized cross validation 5 Comparisons of methods Introduction
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 informationA MODIFIED TSVD METHOD FOR DISCRETE ILL-POSED PROBLEMS
A MODIFIED TSVD METHOD FOR DISCRETE ILL-POSED PROBLEMS SILVIA NOSCHESE AND LOTHAR REICHEL Abstract. Truncated singular value decomposition (TSVD) is a popular method for solving linear discrete ill-posed
More informationAn eigenvalue method using multiple frequency data for inverse scattering problems
An eigenvalue method using multiple frequency data for inverse scattering problems Jiguang Sun Abstract Dirichlet and transmission eigenvalues have important applications in qualitative methods in inverse
More informationarxiv: v1 [physics.comp-ph] 29 Sep 2017
Meshfree Local Radial Basis Function Collocation Method with Image Nodes Seung Ki Baek and Minjae Kim arxiv:179.124v1 [physics.comp-ph] 29 Sep 217 Department of Physics, Pukyong National University, Busan
More informationFast Boundary Knot Method for Solving Axisymmetric Helmholtz Problems with High Wave Number
Copyright 2013 Tech Science Press CMES, vol.94, no.6, pp.485-505, 2013 Fast Boundary Knot Method for Solving Axisymmetric Helmholtz Problems with High Wave Number J. Lin 1, W. Chen 1,2, C. S. Chen 3, X.
More informationA study on regularization parameter choice in Near-field Acoustical Holography
Acoustics 8 Paris A study on regularization parameter choice in Near-field Acoustical Holography J. Gomes a and P.C. Hansen b a Brüel & Kjær Sound and Vibration Measurement A/S, Skodsborgvej 37, DK-285
More informationDetermination of a space-dependent force function in the one-dimensional wave equation
S.O. Hussein and D. Lesnic / Electronic Journal of Boundary Elements, Vol. 12, No. 1, pp. 1-26 (214) Determination of a space-dependent force function in the one-dimensional wave equation S.O. Hussein
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 informationWe first repeat some well known facts about condition numbers for normwise and componentwise perturbations. Consider the matrix
BIT 39(1), pp. 143 151, 1999 ILL-CONDITIONEDNESS NEEDS NOT BE COMPONENTWISE NEAR TO ILL-POSEDNESS FOR LEAST SQUARES PROBLEMS SIEGFRIED M. RUMP Abstract. The condition number of a problem measures the sensitivity
More informationNational Taiwan University
National Taiwan University Meshless Methods for Scientific Computing (Advisor: C.S. Chen, D.L. oung) Final Project Department: Mechanical Engineering Student: Kai-Nung Cheng SID: D9956 Date: Jan. 8 The
More informationAn improved subspace selection algorithm for meshless collocation methods
An improved subspace selection algorithm for meshless collocation methods Leevan Ling 1 and Robert Schaback 2 1 Department of Mathematics, Hong Kong Baptist University, Kowloon Tong, Hong Kong. 2 Institut
More informationCubic spline Numerov type approach for solution of Helmholtz equation
Journal of Linear and Topological Algebra Vol. 03, No. 01, 2014, 47-54 Cubic spline Numerov type approach for solution of Helmholtz equation J. Rashidinia a, H. S. Shekarabi a and M. Aghamohamadi a a Department
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 informationResearch Article Numerical Solution of the Inverse Problem of Determining an Unknown Source Term in a Heat Equation
Applied Mathematics Volume 22, Article ID 39876, 9 pages doi:.55/22/39876 Research Article Numerical Solution of the Inverse Problem of Determining an Unknown Source Term in a Heat Equation Xiuming Li
More informationPerformance of Multiquadric Collocation Method in Solving Lid-driven Cavity Flow Problem with Low Reynolds Number
Copyright c 2006 Tech Science Press CMES, vol.15, no.3, pp.137-146, 2006 Performance of Multiquadric Collocation Method in Solving Lid-driven Cavity Flow Problem with Low Reynolds umber S. Chantasiriwan
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 informationSolving discrete ill posed problems with Tikhonov regularization and generalized cross validation
Solving discrete ill posed problems with Tikhonov regularization and generalized cross validation Gérard MEURANT November 2010 1 Introduction to ill posed problems 2 Examples of ill-posed problems 3 Tikhonov
More informationAN EVOLUTIONARY ALGORITHM TO ESTIMATE UNKNOWN HEAT FLUX IN A ONE- DIMENSIONAL INVERSE HEAT CONDUCTION PROBLEM
Proceedings of the 5 th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15 th July 005. AN EVOLUTIONARY ALGORITHM TO ESTIMATE UNKNOWN HEAT FLUX IN A
More informationBoundary knot method: A meshless, exponential convergence, integration-free, Abstract
Boundary knot method: A meshless, exponential convergence, integration-free, and boundary-only RBF technique W. Chen Department of Mechanical System Engineering, Shinshu University, Wakasato 4-17-1, Nagano
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 informationHIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS
HIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS JASON ALBRIGHT, YEKATERINA EPSHTEYN, AND QING XIA Abstract. Highly-accurate numerical methods that can efficiently
More informationCompact Local Stencils Employed With Integrated RBFs For Fourth-Order Differential Problems
Copyright 2011 Tech Science Press SL, vol.6, no.2, pp.93-107, 2011 Compact Local Stencils Employed With Integrated RBFs For Fourth-Order Differential Problems T.-T. Hoang-Trieu 1, N. Mai-Duy 1 and T. Tran-Cong
More informationThe method of fundamental solutions for Poisson's equation M.A. Golberg Czrc/e
The method of fundamental solutions for Poisson's equation M.A. Golberg Czrc/e ABSTRACT We show how to extend the method of fundamental solutions (MFS) to solve Poisson's equation in 2D without boundary
More informationA Simple Compact Fourth-Order Poisson Solver on Polar Geometry
Journal of Computational Physics 182, 337 345 (2002) doi:10.1006/jcph.2002.7172 A Simple Compact Fourth-Order Poisson Solver on Polar Geometry Ming-Chih Lai Department of Applied Mathematics, National
More informationA FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS
Proceedings of ALGORITMY 2005 pp. 222 229 A FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS ELENA BRAVERMAN, MOSHE ISRAELI, AND ALEXANDER SHERMAN Abstract. Based on a fast subtractional
More informationSolving Pure Torsion Problem and Modelling Radionuclide Migration Using Radial Basis Functions
International Workshop on MeshFree Methods 3 1 Solving Pure Torsion Problem and Modelling Radionuclide Migration Using Radial Basis Functions Leopold Vrankar (1), Goran Turk () and Franc Runovc (3) Abstract:
More informationPreliminary Examination, Numerical Analysis, August 2016
Preliminary Examination, Numerical Analysis, August 2016 Instructions: This exam is closed books and notes. The time allowed is three hours and you need to work on any three out of questions 1-4 and any
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 informationAn Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations
An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations Moshe Israeli Computer Science Department, Technion-Israel Institute of Technology, Technion city, Haifa 32000, ISRAEL Alexander
More informationNon-Singular Method of Fundamental Solutions and its Application to Two Dimensional Elasticity Problems
on-singular Method of Fundamental Solutions and its Application to Two Dimensional Elasticit Problems Qingguo Liu Mentor: Prof. Dr. Božidar Šarler Scope Motivation Introduction Previous Applications of
More informationInverse Source Identification based on Acoustic Particle Velocity Measurements. Abstract. 1. Introduction
The 2002 International Congress and Exposition on Noise Control Engineering Dearborn, MI, USA. August 19-21, 2002 Inverse Source Identification based on Acoustic Particle Velocity Measurements R. Visser
More informationHIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS
HIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS JASON ALBRIGHT, YEKATERINA EPSHTEYN, AND QING XIA Abstract. Highly-accurate numerical methods that can efficiently
More informationHybrid Chebyshev Polynomial Scheme for the Numerical Solution of Partial Differential Equations
The University of Southern Mississippi The Aquila Digital Community Dissertations Summer 8-2016 Hybrid Chebyshev Polynomial Scheme for the Numerical Solution of Partial Differential Equations Balaram Khatri
More informationFirst- and second-order sensitivity analysis schemes by collocation-type Trefftz method
Computer Assisted Mechanics and Engineering Sciences, 4: 477 489, 1997 Copyright c 1997 by Polska Akademia Nauk First- and second-order sensitivity analysis schemes by collocation-type Trefftz method Eisuke
More informationA Two-Side Equilibration Method to Reduce the Condition Number of an Ill-Posed Linear System
Copyright 2013 Tech Science Press CMES, vol.91, no.1, pp.17-42, 2013 A Two-Side Equilibration Method to Reduce the Condition Number of an Ill-Posed Linear System Chein-Shan Liu 1 Abstract: In the present
More informationINTRODUCTION TO PDEs
INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial
More informationOPTIMAL CONVERGENCE RATES FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH POTENTIAL FORCES
OPTIMAL CONVERGENCE RATES FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH POTENTIAL FORCES RENJUN DUAN Department of Mathematics, City University of Hong Kong 83 Tat Chee Avenue, Kowloon, Hong Kong,
More informationEstimation of transmission eigenvalues and the index of refraction from Cauchy data
Estimation of transmission eigenvalues and the index of refraction from Cauchy data Jiguang Sun Abstract Recently the transmission eigenvalue problem has come to play an important role and received a lot
More informationSIMULTANEOUS DETERMINATION OF STEADY TEMPERATURES AND HEAT FLUXES ON SURFACES OF THREE DIMENSIONAL OBJECTS USING FEM
Proceedings of ASME IMECE New ork, N, November -6, IMECE/HTD-43 SIMULTANEOUS DETERMINATION OF STEAD TEMPERATURES AND HEAT FLUES ON SURFACES OF THREE DIMENSIONAL OBJECTS USING FEM Brian H. Dennis and George
More informationMeshfree Approximation Methods with MATLAB
Interdisciplinary Mathematical Sc Meshfree Approximation Methods with MATLAB Gregory E. Fasshauer Illinois Institute of Technology, USA Y f? World Scientific NEW JERSEY LONDON SINGAPORE BEIJING SHANGHAI
More informationMethod of Fundamental Solutions for Scattering Problems of Electromagnetic Waves
Copyright c 2005 Tech Science Press CMES, vol.7, no.2, pp.223-232, 2005 Method of Fundamental Solutions for Scattering Problems of Electromagnetic Waves D.L. Young 1,2 and J.W. Ruan 2 Abstract: The applications
More informationThe method of fundamental solutions for solving Poisson problems
The method of fundamental solutions for solving Poisson problems C.J.S. Alves', C.S. Chen2and B. %der3 'Centro de Matematica Aplicada, Instituto Superior Thaico,Portugal 2Departmentof Mathematical Sciences,
More informationdifferential problems
A numerical study of 2D integrated RBFNs incorporating Cartesian grids for solving 2D elliptic differential problems N. Mai-Duy and T. Tran-Cong Computational Engineering and Science Research Centre Faculty
More informationA numerical solution of a Kawahara equation by using Multiquadric radial basis function
Mathematics Scientific Journal Vol. 9, No. 1, (013), 115-15 A numerical solution of a Kawahara equation by using Multiquadric radial basis function M. Zarebnia a, M. Takhti a a Department of Mathematics,
More informationOptimal Interface Conditions for an Arbitrary Decomposition into Subdomains
Optimal Interface Conditions for an Arbitrary Decomposition into Subdomains Martin J. Gander and Felix Kwok Section de mathématiques, Université de Genève, Geneva CH-1211, Switzerland, Martin.Gander@unige.ch;
More informationFinite Elements. Colin Cotter. January 15, Colin Cotter FEM
Finite Elements January 15, 2018 Why Can solve PDEs on complicated domains. Have flexibility to increase order of accuracy and match the numerics to the physics. has an elegant mathematical formulation
More informationSolving the Laplace Equation by Meshless Collocation Using Harmonic Kernels
Solving the Laplace Equation by Meshless Collocation Using Harmonic Kernels Robert Schaback February 13, 2008 Abstract We present a meshless technique which can be seen as an alternative to the Method
More informationc 1999 Society for Industrial and Applied Mathematics
SIAM J. MATRIX ANAL. APPL. Vol. 21, No. 1, pp. 185 194 c 1999 Society for Industrial and Applied Mathematics TIKHONOV REGULARIZATION AND TOTAL LEAST SQUARES GENE H. GOLUB, PER CHRISTIAN HANSEN, AND DIANNE
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 informationNumerical analysis of heat conduction problems on 3D general-shaped domains by means of a RBF Collocation Meshless Method
Journal of Physics: Conference Series PAPER OPEN ACCESS Numerical analysis of heat conduction problems on 3D general-shaped domains by means of a RBF Collocation Meshless Method To cite this article: R
More informationA collocation method for solving some integral equations in distributions
A collocation method for solving some integral equations in distributions Sapto W. Indratno Department of Mathematics Kansas State University, Manhattan, KS 66506-2602, USA sapto@math.ksu.edu A G Ramm
More informationNumerical Analysis and Methods for PDE I
Numerical Analysis and Methods for PDE I A. J. Meir Department of Mathematics and Statistics Auburn University US-Africa Advanced Study Institute on Analysis, Dynamical Systems, and Mathematical Modeling
More informationSolution of inverse heat conduction equation with the use of Chebyshev polynomials
archives of thermodynamics Vol. 372016, No. 4, 73 88 DOI: 10.1515/aoter-2016-0028 Solution of inverse heat conduction equation with the use of Chebyshev polynomials MAGDA JOACHIMIAK ANDRZEJ FRĄCKOWIAK
More informationThe double layer potential
The double layer potential In this project, our goal is to explain how the Dirichlet problem for a linear elliptic partial differential equation can be converted into an integral equation by representing
More informationSyllabus for Applied Mathematics Graduate Student Qualifying Exams, Dartmouth Mathematics Department
Syllabus for Applied Mathematics Graduate Student Qualifying Exams, Dartmouth Mathematics Department Alex Barnett, Scott Pauls, Dan Rockmore August 12, 2011 We aim to touch upon many topics that a professional
More informationDUAL REGULARIZED TOTAL LEAST SQUARES SOLUTION FROM TWO-PARAMETER TRUST-REGION ALGORITHM. Geunseop Lee
J. Korean Math. Soc. 0 (0), No. 0, pp. 1 0 https://doi.org/10.4134/jkms.j160152 pissn: 0304-9914 / eissn: 2234-3008 DUAL REGULARIZED TOTAL LEAST SQUARES SOLUTION FROM TWO-PARAMETER TRUST-REGION ALGORITHM
More informationNumerical comparison of two boundary meshless methods for water wave problems
Boundary Elements and Other Mesh Reduction Methods XXXVI 115 umerical comparison of two boundary meshless methods for water wave problems Zatianina Razafizana 1,2, Wen Chen 2 & Zhuo-Jia Fu 2 1 College
More informationQ. Zou and A.G. Ramm. Department of Mathematics, Kansas State University, Manhattan, KS 66506, USA.
Computers and Math. with Applic., 21 (1991), 75-80 NUMERICAL SOLUTION OF SOME INVERSE SCATTERING PROBLEMS OF GEOPHYSICS Q. Zou and A.G. Ramm Department of Mathematics, Kansas State University, Manhattan,
More informationThree-Phase Inverse Design Stefan Problem
Three-Phase Inverse Design Stefan Problem Damian S lota Institute of Mathematics Silesian University of Technology Kaszubska 23, 44-1 Gliwice, Poland d.slota@polsl.pl Abstract. The method of the convective
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 informationCalibration and Rescaling Principles for Nonlinear Inverse Heat Conduction and Parameter Estimation Problems
University of Tennessee, Knoxville Trace: Tennessee Research and Creative Exchange Doctoral Dissertations Graduate School 5-2015 Calibration and Rescaling Principles for Nonlinear Inverse Heat Conduction
More informationSolving the complete-electrode direct model of ERT using the boundary element method and the method of fundamental solutions
Solving the complete-electrode direct model of ERT using the boundary element method and the method of fundamental solutions T. E. Dyhoum 1,2, D. Lesnic 1, and R. G. Aykroyd 2 1 Department of Applied Mathematics,
More informationLecture 8: Boundary Integral Equations
CBMS Conference on Fast Direct Solvers Dartmouth College June 23 June 27, 2014 Lecture 8: Boundary Integral Equations Gunnar Martinsson The University of Colorado at Boulder Research support by: Consider
More information3rd Int. Conference on Inverse Problems in Engineering AN INVERSE PROBLEM FOR SLOW VISCOUS INCOMPRESSIBLE FLOWS. University of Leeds
Inverse Problems in Engineering: Theory and Practice 3rd Int. Conference on Inverse Problems in Engineering June 13-18, 1999, Port Ludlow, WA, USA ME06 AN INVERSE PROBLEM FOR SLOW VISCOUS INCOMPRESSIBLE
More informationNonlinear Diffusion. Journal Club Presentation. Xiaowei Zhou
1 / 41 Journal Club Presentation Xiaowei Zhou Department of Electronic and Computer Engineering The Hong Kong University of Science and Technology 2009-12-11 2 / 41 Outline 1 Motivation Diffusion process
More informationA technique based on the equivalent source method for measuring the surface impedance and reflection coefficient of a locally reacting material
A technique based on the equivalent source method for measuring the surface impedance and reflection coefficient of a locally reacting material Yong-Bin ZHANG 1 ; Wang-Lin LIN; Chuan-Xing BI 1 Hefei University
More information8 The SVD Applied to Signal and Image Deblurring
8 The SVD Applied to Signal and Image Deblurring We will discuss the restoration of one-dimensional signals and two-dimensional gray-scale images that have been contaminated by blur and noise. After an
More information8 The SVD Applied to Signal and Image Deblurring
8 The SVD Applied to Signal and Image Deblurring We will discuss the restoration of one-dimensional signals and two-dimensional gray-scale images that have been contaminated by blur and noise. After an
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 informationENERGY METHODS IN IMAGE PROCESSING WITH EDGE ENHANCEMENT
ENERGY METHODS IN IMAGE PROCESSING WITH EDGE ENHANCEMENT PRASHANT ATHAVALE Abstract. Digital images are can be realized as L 2 (R 2 objects. Noise is introduced in a digital image due to various reasons.
More informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationRBF Collocation Methods and Pseudospectral Methods
RBF Collocation Methods and Pseudospectral Methods G. E. Fasshauer Draft: November 19, 24 Abstract We show how the collocation framework that is prevalent in the radial basis function literature can be
More informationOn the Solution of the Elliptic Interface Problems by Difference Potentials Method
On the Solution of the Elliptic Interface Problems by Difference Potentials Method Yekaterina Epshteyn and Michael Medvinsky Abstract Designing numerical methods with high-order accuracy for problems in
More information