Manoj Kumar 1, Pratibha Joshi 2. (Received 31 August 2011, accepted 26 December 2011)
|
|
- Arleen Lambert
- 5 years ago
- Views:
Transcription
1 ISSN (print), (online) International Journal of Nonlinear Science Vol.14(2012) No.1,pp A Mathematical Model and Numerical Solution of a One Dimensional Steady State Heat Conduction Problem by Using High Order Immersed Interface Method on Non- Uniform Mesh Manoj Kumar 1, Pratibha Joshi 2 1 Department of Mathematics, Motilal Nehru National Institute of Technology, Allahabad, India 2 Department of Mathematics, Graphic Era University, Dehradun, India (Received 31 August 2011, accepted 26 December 2011) Abstract: Elliptic interface problems arise in various areas of science and engineering. The steady state heat conduction in layered bodies is one of the largest areas of application of elliptic interface problems. In this paper, a mathematical model and solution of a one dimensional elliptic interface problem which represents a steady state heat conduction problem in composite medium have been discussed by using high order immersed interface method on non uniform mesh. Numerical results show good agreement with exact solution. Keywords: Elliptic interface problems; Steady state heat conduction; Immersed interface methods; Adaptive mesh refinement. 1 Introduction Elliptic interface problems are basically elliptic boundary value problems with discontinuous coefficients or solution across some interface. Standard finite difference methods may not give accurate results for these problems especially near the interface because of the discontinuity. From past few years many mathematicians and researchers are working on developing methods for interface problems. Initially the popular methods were Peskin s immersed boundary method, the smoothing method, the method of harmonic averaging. In 1994, Z. Li and LeVeque [15] developed original immersed interface method in which they incorporated the jump conditions into the finite difference method. The method was second order accurate. The method adds some additional nodes to the numerical stencil which may lead to a nonsymmetric matrix which can not be solved by standard numerical solvers. Since then many researchers have proposed variety of methods [1 14, 16]. Some of them are developed with the modifications in the standard methods so that they can deal with the discontinuities and the singularities. In [9], a decomposed immersed interface method is introduced where a correction term is introduced to standard difference stencil on the right hand side only so that the resulting linear system remains symmetric and diagonally dominant and can be solved by standard solvers. The correction term is decomposed in both Cartesian directions. Another method EJIIM (Explicit Jump Immersed Interface Method)is presented in [10], where high order jumps at the intersection of the interface and coordinate system are introduced as auxiliary unknowns.the interpolation equation is derived for the high order jumps by taking one sided local approximation and the information about jumps. In [10], the grid generation by correcting the finite difference scheme near the interface is avoided. EJIIM was actually motivated by fast iterative IIM introduced by Z. Li in which elliptic equation is preconditioned before using the immersed interface method. There is another efficient approach Matched interface and boundary method [11], in which the solution on each side of the interface is smoothly extended along the interface by some fictitious points. The fictitious values on the fictitious are evaluated from enforcing the jump conditions at the exact position of the interface. Many other approaches have also been developed like finite element method, integral equation method, discontinuous Galerkin approach etc. Elliptic Interface problems have various physical and engineering applications, such as heat conduction in layered medium, alloy solidification, crystal growth, cell and bubble formation, biochemical processing. Corresponding author. address: pratibha.joshi@gmail.com Copyright c World Academic Press, World Academic Union IJNS /633
2 12 International Journal of Nonlinear Science, Vol.14(2012), No.1, pp Heat conduction in composite medium is a wide area of application of elliptic interface problems. Mathematical modelling of heat conduction in composite walls, fins, layered conductor leads to an elliptic interface problem with discontinuous coefficient. In many engineering problems knowledge of temperature distribution and heat flux in materials composed of different layers with different thermal conductivities e.g. construction of furnaces, conductors, thermal protection applications, laminated and fibre reinforced materials, cooling problems is needed. In addition, steady state heat conduction is very helpful for the study of the structural integrity of layered devices in electronic applications. There have been many studies previously on developing techniques for analytical and numerical solution of problems related to steady state heat conduction. [1 3, 12]. In this paper we have first discussed the Zhong s high order IIM [5] and then generalised the method on non uniform mesh and derived the difference formula. Then we applied a one dimensional adaptive grid generaion algorithm to increase the accuracy. Finally we have solved some steady state heat conduction problems in one dimension. The paper has been divided in five sections. In first section we have discussed the mathematical model of a steady state heat conduction problem.in second section the high order IIM has been described. In third section we have obtained the discretization on non uniform mesh. In fourth section we have discussed the adaptive grid generation algorithm which we are going to apply on the scheme. In the last section we have numerically solved some one dimensional steady state heat conduction problems and illustrated the computational results. 2 Mathematical model of steady state heat conduction in one dimensional layered medium In this section, we will discuss the mathematical interpretation of steady state heat conduction equation and derive the heat conduction equation. Due to the existence of temperature gradient within a body, heat energy will flow from the region of high temperature to the region of low temperature. This phenomenon is known as conduction heat transfer, and is described by Fourier s Law q = k T. This equation determines the heat flux vector q for a given temperature profile T and thermal conductivity k. The minus sign ensures that heat flows down the temperature gradient. The heat equation follows from the conservation of energy for a small element within the body, Heat conducted in + Heat generated within = Heat conducted out + Change in energy stored within We can combine the heats conducted in and out into one net heat conducted out term to give, Net heat conducted out = Heat generated within - Change in energy stored within Mathematically, this equation is expressed as, q = q gen de dt where q gen is the power generated per unit volume and e is the internal energy. The change in internal energy e is related to the body s ability to store heat by raising its temperature, given by, (1) de dt dt = ρc dt Here c is the specific heat and ρ is the density. Also, the thermal diffusivity α is related to the thermal conductivity k,the specific heatc and the density ρ by, α = k ρc Now substituting for q using Fourier s Law of heat conduction from above to arrive at the Heat Equation, IJNS for contribution: editor@nonlinearscience.org.uk
3 M. Kumar, P. Joshi: A Mathematical Model and Numerical Solution of a One Dimensional Steady State Heat Conduction 13 In steady state, the heat equation takes the form ( k T ) = q gen ρc dt dt k 2 T + ρc dt = q gen 2 T 1 α dt dt dt = 1 k q gen ( k T ) = q gen Hence a one dimensional steady state heat conduction problem leads to an elliptic boundary value problem of the following form: (β(x)u (x)) = f(x) xϵω (2) where u(x) is the temperature distribution, β(x) is the heat conductivity and f(x) is the source. The domain is a layered medium i.e. it consists of two layers, say Ω 1 and Ω 2 with different heat conductivities, β 1 (x) and β 2 (x). Suppose their common boundary is at x = a. The above problem can be modeled into an elliptic interface problem with discontinuous coefficient β(x) across the interface x = a. Figure 1: One dimensional layered medium. { β1 (x) x Ω β(x) = 1 β 2 (x) x Ω 2 The domain is Ω = Ω 1 Ω 1. Suppose the temperature and flux have jump discontinuity across the interface which are defined as where [u] x=a = u + u [βu x ] x=a = β + u x + β u x u = lim x a u(x) and u+ = lim x a +u(x) β = lim β(x) and β + = lim β(x) x a x a + If the jumps defined above are equal to zero, the contact between the layers is called ideal. 3 High order immersed interface method to solve elliptic interface problem Xiaolin Zhong presented a high order immersed interface method for the simulation of multi phase flow [5]. The benefits of the method are the following: I. We can derive finite difference scheme of arbitrary order with the same methodology. II. It needs only two physical jump conditions of the variables and their first derivatives to achieve second or higher order accuracy. III. The method is applicable to variety of interface problems since the finite difference formulas can be expressed in explicit form. IJNS homepage:
4 14 International Journal of Nonlinear Science, Vol.14(2012), No.1, pp In the above scheme, standard finite difference formula has been used at regular points whereas at irregular points, special finite difference formulae have been derived by using Taylor expansion and matched polynomial interpolation. The discretization has been done on a uniform mesh. Suppose we have a one dimensional second order elliptic boundary value problem with following jump conditions [αu] = α + u + α u = A (3) [βu x ] = β + u x + β u x = B (4) The interface is at x Γ where x i and x i+1 are the irregular points i.e. x i < x Γ < x i+1. In [5], to derive standard finite Figure 2: Stencil near the interface. difference formulae by using matched polynomial interpolation, two matched polynomials are used to interpolate both sides of the interface. Due to the discontinuity of the solution across the interface single interpolation polynomial can not be used. For discretization, n points on left side of the interface and m points on the right side have been taken which will lead to a n + m points grid stencil. The polynomial interpolating left hand side of the interface can be written as where a n is an undetermined coefficient and P (x) = l k (x)u i+k + a n R(x) (5) n+1 k=0 R(x) = (x x i+k ) (6) n+1 l k (x) is the Lagrange polynomial interpolating n points on left side of the interface n+1 l=0,l k l k (x) = (x x i+l) n+1 l=0,l k (x i+k x i+l ) k=0 The polynomial for the right hand side of the interface can be written in similar manner as m P + (x) = h k (x)u i+k + b m Q(x) (8) where b m is an undetermined coefficient and k=1 Q(x) = m (x x i+k ) (9) k=1 (7) IJNS for contribution: editor@nonlinearscience.org.uk
5 M. Kumar, P. Joshi: A Mathematical Model and Numerical Solution of a One Dimensional Steady State Heat Conduction 15 h k (x) is the Lagrange polynomial interpolating m points on the right hand side m l=1,l k h k (x) = (x x i+l) m l=1,l k (x i+k x i+l ) (10) The undetermined coefficients a n and b m are decided by the jump conditions. Using the values of (5) and (8) in jump conditions (3) and (4) we have and ( m ) α + h k (x Γ )u i+k + b m Q(x Γ ) k=1 ( m ) β + h k(x)u i+k + b m Q (x) k=1 Solving equations (11) and (12) we get a n = m k= n+1 ( n+1 α k=0 ( n+1 β k=0 γ k u i+k + ξ A A + ξ B B l k (x Γ )u i+k + a n R(x Γ ) l k(x)u i+k + a n R (x) ) ) = A (11) = B (12) where and { 1 J γ k = (α β + Q (x Γ )l k (x Γ )) α + β Q(x Γ )l k (x Γ)) (k = n + 1,..., 0) α + β + J ( Q (x Γ )h k (x Γ ) + Q(x Γ )h k (x Γ)) (k = 1,..., m) ξ A = 1 J β+ Q (x Γ ) ξ B = 1 J α+ Q(x Γ )) J = α + β R (x Γ )Q(x Γ ) α β + R(x Γ )Q (x Γ ) Similarly coefficientb m can be evaluated. Value of b m is not determined here because it is not needed for evaluating difference scheme at irregular point i. Hence the r th derivative at irregular point i can be approximated as ( d r ) ( y d r P ) (x) dx r i dx r (13) x=x i Once we have determined the undetermined coefficients we get the full expression for the polynomials. We can derive expression for any order of derivative from the equation (13). Hence differentiating the polynomial r times we get ( d r ) y dx r i m k= n+1 d k u i+k + d A A + d B B (14) where d k = { l (r) k x(i) + γ kr (r) x(i) (k = n + 1,..., 0) γ k R (r) x(i) (k = 1,..., m) d A = ξ A R(r) x(i) d B = ξ B R(r) x(i) This is the finite difference formula for the irregular point i on the left hand side of the interface. We can get the difference formula at irregular point i + 1 by a coordinate transformation as described in [5]. At regular points we can use the standard central difference formula. u xx = u i 1 2u i + u i+1 h 2 IJNS homepage:
6 16 International Journal of Nonlinear Science, Vol.14(2012), No.1, pp High order immersed interface method on non uniform mesh In the scheme [5] the mesh is taken uniform. For n = m = 2 the difference scheme (14) gives the approximation at irregular point i as Here d k s are defined as where du = d 1u i 1 + d 0 u 0 + d 1 u i+1 + d 2 u i+2 + d A A + hd B B + u i u i 1 dx 2h h ( d 2 ) u dx 2 = d 1u i 1 + d 0 u 0 + d 1 u i+1 + d 2 u i+2 + d A A + hd B B i h 2 (15) d 1 = 1 ( ( cα 3σ 2σ 2 ) ( c β 2 + 3σ σ 2 )) D d 0 = 1 ( ( cα 3 σ + 2σ 2 ) ( c β 2 3σ + σ 2 )) D d 1 = 1 ( 4 4σ + σ 2 ) D d 2 = 1 ( 1 + 2σ σ 2 ) D d A = 1 α + ( 3 + 2σ) D d B = 1 ( 2 3σ + σ 2 ) β + D D = 1 2 c α = α α + c β = β β + ( cβ ( 2 + σ 5σ 2 + 2σ 3) c α ( 3σ σ 2 + 2σ 3)) Here h is the uniform mesh length. Now we propose the discretization on non uniform mesh. Suppose we have a grid {x i }, i = 1, 2,... with the grid lengths {h i }, i = 1, 2,... where h i = x i+1 x i and interface is located at x Γ = x i + σh. Then we can approximate the second order derivative of unknown as Figure 3: Four point stencil taking two points on both side of the interface. ( d 2 ) u dx 2 = t 1 u i 1 + t 0 u 0 + t 1 u i+1 + t 2 u i+2 + t A A + t B B i IJNS for contribution: editor@nonlinearscience.org.uk
7 M. Kumar, P. Joshi: A Mathematical Model and Numerical Solution of a One Dimensional Steady State Heat Conduction 17 The t k s are defined as T = 1 2 (c β (2σh i + h i 1 ) ((1 σ)h i ((1 σ)h i + h i+1 )) + c α (2(1 σ)h i + h i+1 ) (σh i (σh i + h i 1 ))) σh c α (2(1 σ)h i + h i+1 ) i 1 h i 1 + c β (((1 σ)h i )((1 σ)h i + h i+1 )) h i 1 t 1 = T c α ( 2(1 σ)h i h i+1 ) hi 1+σh i 1 h i 1 c β (((1 σ)h i )((1 σ)h i + h i+1 )) h i 1 t 0 = T (2(1 σ)h i + h i+1 ) ((1 σ)hi )+h i+1 1 h i+1 (((1 σ)h i )((1 σ)h i + h i+1 )) h i 1 t 1 = (16) T ( 2(1 σ)h i + h i+1 ) (1 σ)hi 1 h i+1 + (((1 σ)h i )((1 σ)h i + h i+1 )) h i 1 t 2 = T t A = β+ ( 2(1 σ)h i h i+1 ) T t B = α+ ((1 σ)h i ((1 σ)h i + h i+1 )) T We can find the approximation for the irregular point i + 1 by a coordinate transformation x = x As a result we get the finite difference formula as ( d 2 ) u dx 2 = t 2 u i 1 + t 1 u 0 + t 0 u i+1 + t 1 u i+2 t A A + t B B (17) i where t k s can be determined by puttingσ = 1 σ and interchanging places of h i+1 and h i 1 in equations (16). 5 Adaptive grid generation in the immersed interface method Adaptive grid generation is a popular and efficient method to increase accuracy of difference schemes. Many researchers have done significant work in generating and applying adaptive grids in the numerical solution of differential equations [7, 8]. It is basically special type of grid generation where we add some more grid points in the interval in which the solution has large variations such as peaks or boundary layers etc. In interface problems generally the interval having largest variation in the solution is the interval containing the interface. We have applied this concept in the difference scheme [5]. Since the intervals are not of uniform length after the adaptive grid generation; we have determined the scheme on non uniform mesh. There are various methods to generate adaptive grids. The common property of most of these methods is that they divide the solution domain into subintervals such that some positive weight function has approximately same value over each subinterval. The choice of the weight function depends on the first derivative or the second derivative or the truncation error. The algorithm we have applied in the difference scheme above is presented in [6]. The algorithm has the following steps: I. Generate a uniform mesh in the given solution domain. II. Compute the error in each subinterval. III. Determine the subinterval of maximum error. IV. Include the middle point into that subinterval. V. Repeat the procedure until the stopping criteria is satisfied. We have applied the algorithm on the interval containing the interface. Since the immersed interface method developed in [5] is built on uniform grids, it is not possible to apply the adaptive grid generation in the same scheme. Therefore, we have determined the difference scheme on non uniform mesh. At regular point we have used the standard finite difference formula on non uniform mesh. u xx = 2[h i 1u(x i+1 ) (h i + h i 1 )u(x i ) + h i u(x i 1 )] h i h i 1 (h i + h i 1 ) IJNS homepage:
8 18 International Journal of Nonlinear Science, Vol.14(2012), No.1, pp We will prove the efficiency of the proposed scheme in next section by solving one dimensional heat conduction problems. The results are very satisfactory. 6 Numerical results To illustrate our proposed scheme we consider two one dimensional steady state heat conduction problems. Suppose we have a very thin rod of unit length constructed with two materials having different conductivities. The contact of the two layers is at x = α.i.e. { k1 x < α k = x > α where k 1 and k 2 are constants. k 2 Problem 1 In this problem the temperature distribution is governed by the following elliptic problem (ku x ) x = 12x 2, 0 < x < 1 where the temperature at end points of the rod is given by u(0) = 0, u(1) = 1 ( ) α. k 2 k 1 k 2 The jump conditions for the temperature and the flux along the interface are The exact solution of the problem is u(x) = [u] x=α = 0, [ku x ] x=α = 0 x 4 k 1 x 4 k 2 + ( 1 k 1 1 k 2 )α 4 x < α x > α Now we apply the new immersed interface method described in section 4 and solve the problem by taking different values of.we have developed a computer program of the proposed method applied on the above problem in MATLAB 7. The algorithm for adaptive grid generation discussed in section 5 terminates when some stopping criterion is satisfied which can be either required level of accuracy or number of subintervals etc. In our MATLAB program, number of subintervals required has been chosen as the criterion to terminate. The order of the scheme can be evaluated by en log p = e 2n log(2) (18) where e n is the maximum norm defined by e n = max u(x i ) U(x i ) (19) i where u(x) is the approximate solution and U(x) is the exact solution of the problem. The error ratio can be evaluated by e n e 2n. (20) The grid refinement analysis is displayed in Table I. We have taken three different sets of constants and evaluated the maximum absolute error in each case. We have taken two meshes for N = 10 and N = 20. In the proposed method the mesh is unequally distributed whereas in IIM [5] the mesh has equal subintervals. The order of convergence and the error ratio can easily be evaluated by the formulas. In general the order of proposed method is O(h 2 ). The resulting IJNS for contribution: editor@nonlinearscience.org.uk
9 M. Kumar, P. Joshi: A Mathematical Model and Numerical Solution of a One Dimensional Steady State Heat Conduction 19 Figure 4: Comparison of numerical solution by proposed IIM and exact solution of Problem-1 for k 1 = 1, k 2 = 100, α = 2/3 and N = 20 unequal subintervals. Figure 5: Comparison of numerical solution by proposed IIM and exact solution of Problem-1 for k 1 = 1, k 2 = 1000, α = 2/3 and N = 20 unequal subintervals. Table 1: Maximum Absolute Error. Constants e n (Proposed IIM) N = 10 N = 20 k 1 = 1, k 2 = 100, α = 2/ k 1 = 1, k 2 = 500, α = 2/ k 1 = 1, k 2 = 1000, α = 2/ IJNS homepage:
10 20 International Journal of Nonlinear Science, Vol.14(2012), No.1, pp system of equations comes nonsymmetric in both cases because of the change of difference scheme at irregular points. However, the facts that the method is very flexible to apply in different interface problems and one can easily construct high order schemes from the given schemes (in Section 2 and 3) make the methods overweight the lose of symmetry in the resultant matrix. We have used the operator in MATLAB for solving the nonsymmetric coefficient matrix because it is very popular among the people working in MATLAB because of its efficiency. Recently many multigrid solvers and iterative methods are being used for nonsymmetric sparse system of equations which gives fast and accurate solution e.g. BiConjugate Gradient Stabilized (BiCGSTAB), Generalized Minimum Residual Method - an orthogonalization and least squares method (GMRES) etc which can also be used here instead of the direct solver (\ in MATLAB). The computational errors of the above methods depend on the location of the interface i.e. the value of σ.the value of σ changes with changes in number of subintervals. Problem 2 In this problem the temperature distribution is governed by the following equation with boundary conditions and thermal conductivities of the two materials are x (k x u(x)) = 0, x (0, 1) u(0) = 0, u(1) = 1 k = { 1 x < α 2 x > α The exact solution of the problem is 2x x < α u(x) = α + x 1 q x > α This problem is solved on two non uniform meshes with 21 and 41 intervals. Computational results are given in Table II. We have considered three cases of the problem with three different positions of the interface. In each situation the proposed method works well. Figure 6: Comparison of numerical solution by proposed IIM and exact solution of Problem-2 for N = 41 unequal subintervals at α = 6/7. Table 2: Maximum Absolute Error. Position of the interface e n (Proposed IIM) N = 21 N = 41 α = 3/ α = 6/ α = 7/ The computational results show good agreement with the exact solution for both the problems. Mostly the proposed method gives better result than the high order IIM [5]. The adaptive grid is constructed as to give the required accuracy. But one has to take care as sometimes the convergence could be lost while constructing it. IJNS for contribution: editor@nonlinearscience.org.uk
11 M. Kumar, P. Joshi: A Mathematical Model and Numerical Solution of a One Dimensional Steady State Heat Conduction 21 Figure 7: Comparison of numerical solution by proposed IIM and exact solution of Problem-2 for N = 41 unequal subintervals at α = 3/4. 7 Conclusion We have discussed the mathematical model of one dimensional steady state heat conduction problems and solved them by using high order IIM on non uniform grid applying the one dimensional adaptive grid generation algorithm described in [6]. The proposed method is well suited for the above problems.the proposed non uniform discretization can be used with many non uniform grids and adaptive grids which produces more accurate results mostly. In [5], Xiaolin Zhong was motivated to derive the high order IIM because of the need of an efficient method to solve problem in two phase flow. However the method is applicable to some steady state heat conduction problems too. The advantage of using these methods in steady state heat conduction problem is that they need only two jump condition and they are easy to apply on different types of problems. We can easily generate difference scheme of arbitrary order by the procedure described in 3. References [1] L.G. Stanley. A Senstivity Equation method for Molding Processes. Proceedings of the 2000 IEEE International Conference on Control Applications Anchorage, Alaska, USA,September [2] Z. M. Seyidmamedov and E. Ozbilge. A mathematical model and numerical solution of interface problems for steady state heat conduction. Mathematical Problems in Engineering., 2006(2006):1 18. [3] Okey Oseloka Onyejekwe. Heat conduction in composite media: a boundary integral approach. Computers and Chemical Engineering., 26(2002): [4] Z. Li. The Immersed Interface Method - A Numerical Approach for Partial Differential Equations with Interfaces. Ph.D. thesis, University of Washington [5] Xiaolin Zhong. A new high-order immersed interface method for solving elliptic equations with imbedded interface of discontinuity. Journal of Computational Physics., 225(2007): [6] M. M. Ahmed Khodier and Y. Adel Hassan. One dimensional adaptive grid generation. Internat. J. Math. & Math. Sci., 20(1997)(3): [7] J. F. Thompson. A Survey of Dynamically Adaptive Grid in the Numerical Solution of Partial Differential Equations. Appl. Numer. Math., (1985):3 27. [8] V. E. Denny and It. B. Landis. A New Method for Solving Two Point Boundary Value Problems Using Optimal Node Distribution. J. of comp. Phys., 9(1972): [9] Berthelsen,Andreas Petter. A decomposed immersed interface method for variable coefficient elliptic equations with non-smooth and discontinuous solutions. Journal of Computational Physics, 197(2004): IJNS homepage:
12 22 International Journal of Nonlinear Science, Vol.14(2012), No.1, pp [10] A. Wiegmann and K. Bube. The explicit-jump immersed interface method: finite difference methods for PDEs with piecewise smooth solutions. SIAM J. Numer. Anal., 37(2000): [11] Y. C. Zhou and G. W. Wei. On the fictitious -domain and interpolation formulations of the matched interface and boundary method. Journal of Computational Physics., 219(2006): [12] S. Naser Al-Huniti and M. A. Al-Nimr. Steady-State Thermoelastic Behavior of a Two-Anisotropic Layer Thick Plate Strip. International Journal for Computational Methods in Engineering Science and Mechanics., 7(2006): [13] Sumir Chandra and Manish Parashar. Addressing spatiotemporal and computational heterogeneity in structured adaptive mesh refinement. Computing and Visualization in Science., 9(2006): [14] R. Kornhuber, R. Krause, O. Sander, P. Deuflhard and S. Ertel. A monotone multigrid solver for two body contact problems in biomechanics. Computing and Visualization in Science., 11(2008):3 15. [15] R. J. LeVeque and Z. Li. The immersed interface method for elliptic equations with discontinuous coefficients and singular sources. SIAM J. Numer. Anal., 31(1994): [16] Manoj Kumar and Pratibha Joshi. Some Numerical Techniques for Solving Elliptic Interface Problems. Numerical Methods for Partial Differential Equations., 28(2012)(1): [17] Pratibha Joshi, Manoj Kumar. Mathematical Model and Computer Simulation of Three Dimensional Thin Film Elliptic Interface Problems. Computers and Mathematics with Applications,63(1)(2012): IJNS for contribution: editor@nonlinearscience.org.uk
ANALYSIS AND NUMERICAL METHODS FOR SOME CRACK PROBLEMS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING, SERIES B Volume 2, Number 2-3, Pages 155 166 c 2011 Institute for Scientific Computing and Information ANALYSIS AND NUMERICAL METHODS FOR SOME
More informationExplicit Jump Immersed Interface Method: Documentation for 2D Poisson Code
Eplicit Jump Immersed Interface Method: Documentation for 2D Poisson Code V. Rutka A. Wiegmann November 25, 2005 Abstract The Eplicit Jump Immersed Interface method is a powerful tool to solve elliptic
More informationLocal Mesh Refinement with the PCD Method
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 1, pp. 125 136 (2013) http://campus.mst.edu/adsa Local Mesh Refinement with the PCD Method Ahmed Tahiri Université Med Premier
More informationABSTRACT. coefficients. The main motivation is to get not only a second order accurate solution but
ABSTRACT CHEN, GUANYU. Accurate Gradient Computation for Elliptic Interface Problems with Discontinuous and Variable Coefficients. (Under the direction of Dr. Zhilin Li.) A new numerical method is proposed
More information(Received 05 August 2013, accepted 15 July 2014)
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.18(2014) No.1,pp.71-77 Spectral Collocation Method for the Numerical Solution of the Gardner and Huxley Equations
More informationAccurate Gradient Computation for Interface Problems
Accurate Gradient Computation for Interface Problems Nnenna Anoruo, Kevin Luna, William Reese, Justin Rivera Faculty Advsior: Zhilin Li North Carolina State University February 11, 016 Outline Problem
More informationA finite difference Poisson solver for irregular geometries
ANZIAM J. 45 (E) ppc713 C728, 2004 C713 A finite difference Poisson solver for irregular geometries Z. Jomaa C. Macaskill (Received 8 August 2003, revised 21 January 2004) Abstract The motivation for this
More informationHomotopy Perturbation Method for the Fisher s Equation and Its Generalized
ISSN 749-889 (print), 749-897 (online) International Journal of Nonlinear Science Vol.8(2009) No.4,pp.448-455 Homotopy Perturbation Method for the Fisher s Equation and Its Generalized M. Matinfar,M. Ghanbari
More informationOn Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities
On Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities Heiko Berninger, Ralf Kornhuber, and Oliver Sander FU Berlin, FB Mathematik und Informatik (http://www.math.fu-berlin.de/rd/we-02/numerik/)
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 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 informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
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 informationAN OVERVIEW OF THE IMMERSED INTERFACE METHOD AND ITS APPLICATIONS
. AN OVERVIEW OF THE IMMERSED INTERFACE METHOD AND ITS APPLICATIONS ZHILIN LI Abstract. Interface problems have many applications. Mathematically, interface problems usually lead to differential equations
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 9
Spring 2015 Lecture 9 REVIEW Lecture 8: Direct Methods for solving (linear) algebraic equations Gauss Elimination LU decomposition/factorization Error Analysis for Linear Systems and Condition Numbers
More informationA Numerical Algorithm for Finding Positive Solutions of Superlinear Dirichlet Problem
ISSN 1479-3889 (print), 1479-3897 (online) International Journal of Nonlinear Science Vol.3(2007) No.1,pp.27-31 A Numerical Algorithm for Finding Positive Solutions of Superlinear Dirichlet Problem G.A.Afrouzi
More informationA MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY
A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,
More informationAnalytical Solution of BVPs for Fourth-order Integro-differential Equations by Using Homotopy Analysis Method
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.9(21) No.4,pp.414-421 Analytical Solution of BVPs for Fourth-order Integro-differential Equations by Using Homotopy
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 informationNumerical Solution of Second Order Singularly Perturbed Differential Difference Equation with Negative Shift. 1 Introduction
ISSN 749-3889 print), 749-3897 online) International Journal of Nonlinear Science Vol.8204) No.3,pp.200-209 Numerical Solution of Second Order Singularly Perturbed Differential Difference Equation with
More informationAspects of Multigrid
Aspects of Multigrid Kees Oosterlee 1,2 1 Delft University of Technology, Delft. 2 CWI, Center for Mathematics and Computer Science, Amsterdam, SIAM Chapter Workshop Day, May 30th 2018 C.W.Oosterlee (CWI)
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 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 informationINTRODUCTION TO FINITE ELEMENT METHODS ON ELLIPTIC EQUATIONS LONG CHEN
INTROUCTION TO FINITE ELEMENT METHOS ON ELLIPTIC EQUATIONS LONG CHEN CONTENTS 1. Poisson Equation 1 2. Outline of Topics 3 2.1. Finite ifference Method 3 2.2. Finite Element Method 3 2.3. Finite Volume
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 informationA THEORETICAL INTRODUCTION TO NUMERICAL ANALYSIS
A THEORETICAL INTRODUCTION TO NUMERICAL ANALYSIS Victor S. Ryaben'kii Semyon V. Tsynkov Chapman &. Hall/CRC Taylor & Francis Group Boca Raton London New York Chapman & Hall/CRC is an imprint of the Taylor
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 informationCubic B-spline Collocation Method for Fourth Order Boundary Value Problems. 1 Introduction
ISSN 1749-3889 print, 1749-3897 online International Journal of Nonlinear Science Vol.142012 No.3,pp.336-344 Cubic B-spline Collocation Method for Fourth Order Boundary Value Problems K.N.S. Kasi Viswanadham,
More informationINTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR GRIDS. 1. Introduction
Trends in Mathematics Information Center for Mathematical Sciences Volume 9 Number 2 December 2006 Pages 0 INTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR
More informationICES REPORT A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws
ICES REPORT 7- August 7 A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws by Todd Arbogast, Chieh-Sen Huang, and Xikai Zhao The Institute for Computational Engineering and Sciences The
More informationChapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS
Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS 5.1 Introduction When a physical system depends on more than one variable a general
More informationFixed point iteration and root finding
Fixed point iteration and root finding The sign function is defined as x > 0 sign(x) = 0 x = 0 x < 0. It can be evaluated via an iteration which is useful for some problems. One such iteration is given
More informationFast Structured Spectral Methods
Spectral methods HSS structures Fast algorithms Conclusion Fast Structured Spectral Methods Yingwei Wang Department of Mathematics, Purdue University Joint work with Prof Jie Shen and Prof Jianlin Xia
More informationChapter 6. Finite Element Method. Literature: (tiny selection from an enormous number of publications)
Chapter 6 Finite Element Method Literature: (tiny selection from an enormous number of publications) K.J. Bathe, Finite Element procedures, 2nd edition, Pearson 2014 (1043 pages, comprehensive). Available
More informationA Fourth Order Accurate Discretization for the Laplace and Heat Equations on Arbitrary Domains, with Applications to the Stefan Problem
A Fourth Order Accurate Discretization for the Laplace and Heat Equations on Arbitrary Domains, with Applications to the Stefan Problem Frédéric Gibou Ronald Fedkiw April 27, 2004 Abstract In this paper,
More informationBeyond Wiener Askey Expansions: Handling Arbitrary PDFs
Journal of Scientific Computing, Vol. 27, Nos. 1 3, June 2006 ( 2005) DOI: 10.1007/s10915-005-9038-8 Beyond Wiener Askey Expansions: Handling Arbitrary PDFs Xiaoliang Wan 1 and George Em Karniadakis 1
More informationHigh-order di erence potentials methods for 1D elliptic type models. Yekaterina Epshteyn and Spencer Phippen
High-order di erence potentials methods for 1D elliptic type models Yekaterina Epshteyn and Spencer Phippen Department of Mathematics, The University of Utah, Salt Lake City, UT, 84112 Abstract Numerical
More informationIntegration, differentiation, and root finding. Phys 420/580 Lecture 7
Integration, differentiation, and root finding Phys 420/580 Lecture 7 Numerical integration Compute an approximation to the definite integral I = b Find area under the curve in the interval Trapezoid Rule:
More informationarxiv: v1 [math.na] 11 Jul 2011
Multigrid Preconditioner for Nonconforming Discretization of Elliptic Problems with Jump Coefficients arxiv:07.260v [math.na] Jul 20 Blanca Ayuso De Dios, Michael Holst 2, Yunrong Zhu 2, and Ludmil Zikatanov
More informationThe Conjugate Gradient Method
The Conjugate Gradient Method Classical Iterations We have a problem, We assume that the matrix comes from a discretization of a PDE. The best and most popular model problem is, The matrix will be as large
More information2.2. Methods for Obtaining FD Expressions. There are several methods, and we will look at a few:
.. Methods for Obtaining FD Expressions There are several methods, and we will look at a few: ) Taylor series expansion the most common, but purely mathematical. ) Polynomial fitting or interpolation the
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 informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationSolutions of Burgers Equation
ISSN 749-3889 (print, 749-3897 (online International Journal of Nonlinear Science Vol.9( No.3,pp.9-95 Solutions of Burgers Equation Ch. Srinivasa Rao, Engu Satyanarayana Department of Mathematics, Indian
More information3.1 Interpolation and the Lagrange Polynomial
MATH 4073 Chapter 3 Interpolation and Polynomial Approximation Fall 2003 1 Consider a sample x x 0 x 1 x n y y 0 y 1 y n. Can we get a function out of discrete data above that gives a reasonable estimate
More informationAn Immersed Boundary Method for Restricted Diffusion with Permeable Interfaces
An Immersed Boundary Method for Restricted Diffusion with Permeable Interfaces Huaxiong Huang Kazuyasu Sugiyama Shu Takagi April 6, 009 Keywords: Restricted diffusion; Permeable interface; Immersed boundary
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 23: GMRES and Other Krylov Subspace Methods; Preconditioning
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 23: GMRES and Other Krylov Subspace Methods; Preconditioning Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 18 Outline
More informationA NUMERICAL METHOD FOR SOLVING THE VARIABLE COEFFICIENT WAVE EQUATION WITH INTERFACE JUMP CONDITIONS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 6, Number, Pages 7 c 9 Institute for Scientific Computing and Information A NUMERICAL METHOD FOR SOLVING THE VARIABLE COEFFICIENT WAVE EQUATION
More informationA Differential Quadrature Algorithm for the Numerical Solution of the Second-Order One Dimensional Hyperbolic Telegraph Equation
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.13(01) No.3,pp.59-66 A Differential Quadrature Algorithm for the Numerical Solution of the Second-Order One Dimensional
More informationFinite difference method for elliptic problems: I
Finite difference method for elliptic problems: I Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationDiscretization of PDEs and Tools for the Parallel Solution of the Resulting Systems
Discretization of PDEs and Tools for the Parallel Solution of the Resulting Systems Stan Tomov Innovative Computing Laboratory Computer Science Department The University of Tennessee Wednesday April 4,
More informationZhenjiang, Jiangsu, , P.R. China (Received 7 June 2010, accepted xx, will be set by the editor)
ISSN 1749-3889 print), 1749-3897 online) International Journal of Nonlinear Science Vol.132012) No.3,pp.380-384 Fractal Interpolation Functions on the Stability of Vertical Scale Factor Jiao Xu 1, Zhigang
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 informationSolution Methods. Steady State Diffusion Equation. Lecture 04
Solution Methods Steady State Diffusion Equation Lecture 04 1 Solution methods Focus on finite volume method. Background of finite volume method. Discretization example. General solution method. Convergence.
More informationOn some numerical convergence studies of mixed finite element methods for flow in porous media
On some numerical convergence studies of mixed finite element methods for flow in porous media Gergina Pencheva Abstract We consider an expanded mixed finite element method for solving second-order elliptic
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 3: Finite Elements in 2-D Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods 1 / 18 Outline 1 Boundary
More informationHigh Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation
High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation Faheem Ahmed, Fareed Ahmed, Yongheng Guo, Yong Yang Abstract This paper deals with
More informationFinite Difference Methods (FDMs) 1
Finite Difference Methods (FDMs) 1 1 st - order Approxima9on Recall Taylor series expansion: Forward difference: Backward difference: Central difference: 2 nd - order Approxima9on Forward difference: Backward
More informationUpscaling Wave Computations
Upscaling Wave Computations Xin Wang 2010 TRIP Annual Meeting 2 Outline 1 Motivation of Numerical Upscaling 2 Overview of Upscaling Methods 3 Future Plan 3 Wave Equations scalar variable density acoustic
More informationCompact High Order Finite Difference Stencils for Elliptic Variable Coefficient and Interface Problems
Compact High Order Finite Difference Stencils for Elliptic Variable Coefficient and Interface Problems Daniel Ritter 1, Ulrich Rüde 1, Björn Gmeiner 1, Rochus Schmid 2 Copper Mountain, March 18th, 2013
More informationPartially Penalized Immersed Finite Element Methods for Parabolic Interface Problems
Partially Penalized Immersed Finite Element Methods for Parabolic Interface Problems Tao Lin, Qing Yang and Xu Zhang Abstract We present partially penalized immersed finite element methods for solving
More informationParametric Spline Method for Solving Bratu s Problem
ISSN 749-3889 print, 749-3897 online International Journal of Nonlinear Science Vol4202 No,pp3-0 Parametric Spline Metod for Solving Bratu s Problem M Zarebnia, Z Sarvari 2,2 Department of Matematics,
More informationHigh order, finite volume method, flux conservation, finite element method
FLUX-CONSERVING FINITE ELEMENT METHODS SHANGYOU ZHANG, ZHIMIN ZHANG, AND QINGSONG ZOU Abstract. We analyze the flux conservation property of the finite element method. It is shown that the finite element
More informationSolving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners
Solving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners Eugene Vecharynski 1 Andrew Knyazev 2 1 Department of Computer Science and Engineering University of Minnesota 2 Department
More informationMultigrid absolute value preconditioning
Multigrid absolute value preconditioning Eugene Vecharynski 1 Andrew Knyazev 2 (speaker) 1 Department of Computer Science and Engineering University of Minnesota 2 Department of Mathematical and Statistical
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Discretization of Boundary Conditions Discretization of Boundary Conditions On
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT 16-02 The Induced Dimension Reduction method applied to convection-diffusion-reaction problems R. Astudillo and M. B. van Gijzen ISSN 1389-6520 Reports of the Delft
More informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
More informationAdaptive Time Space Discretization for Combustion Problems
Konrad-Zuse-Zentrum fu r Informationstechnik Berlin Takustraße 7 D-14195 Berlin-Dahlem Germany JENS LANG 1,BODO ERDMANN,RAINER ROITZSCH Adaptive Time Space Discretization for Combustion Problems 1 Talk
More informationNumerical Solution I
Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions
More informationAdditive Manufacturing Module 8
Additive Manufacturing Module 8 Spring 2015 Wenchao Zhou zhouw@uark.edu (479) 575-7250 The Department of Mechanical Engineering University of Arkansas, Fayetteville 1 Evaluating design https://www.youtube.com/watch?v=p
More informationPreface to the Second Edition. Preface to the First Edition
n page v Preface to the Second Edition Preface to the First Edition xiii xvii 1 Background in Linear Algebra 1 1.1 Matrices................................. 1 1.2 Square Matrices and Eigenvalues....................
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More informationHIGH ACCURACY NUMERICAL METHODS FOR THE SOLUTION OF NON-LINEAR BOUNDARY VALUE PROBLEMS
ABSTRACT Of The Thesis Entitled HIGH ACCURACY NUMERICAL METHODS FOR THE SOLUTION OF NON-LINEAR BOUNDARY VALUE PROBLEMS Submitted To The University of Delhi In Partial Fulfillment For The Award of The Degree
More informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationSolution of the Coupled Klein-Gordon Schrödinger Equation Using the Modified Decomposition Method
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.4(2007) No.3,pp.227-234 Solution of the Coupled Klein-Gordon Schrödinger Equation Using the Modified Decomposition
More informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
More informationPart I. Discrete Models. Part I: Discrete Models. Scientific Computing I. Motivation: Heat Transfer. A Wiremesh Model (2) A Wiremesh Model
Part I: iscrete Models Scientific Computing I Module 5: Heat Transfer iscrete and Continuous Models Tobias Neckel Winter 04/05 Motivation: Heat Transfer Wiremesh Model A Finite Volume Model Time ependent
More informationNew Solutions for Some Important Partial Differential Equations
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.4(2007) No.2,pp.109-117 New Solutions for Some Important Partial Differential Equations Ahmed Hassan Ahmed Ali
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 37, pp. 166-172, 2010. Copyright 2010,. ISSN 1068-9613. ETNA A GRADIENT RECOVERY OPERATOR BASED ON AN OBLIQUE PROJECTION BISHNU P. LAMICHHANE Abstract.
More informationJournal of Computational and Applied Mathematics. Multigrid method for solving convection-diffusion problems with dominant convection
Journal of Computational and Applied Mathematics 226 (2009) 77 83 Contents lists available at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elsevier.com/locate/cam
More informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More informationRecovery-Based A Posteriori Error Estimation
Recovery-Based A Posteriori Error Estimation Zhiqiang Cai Purdue University Department of Mathematics, Purdue University Slide 1, March 2, 2011 Outline Introduction Diffusion Problems Higher Order Elements
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 informationMath 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework
Math 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework Jan Mandel University of Colorado Denver May 12, 2010 1/20/09: Sec. 1.1, 1.2. Hw 1 due 1/27: problems
More informationA cartesian grid finite volume method for the solution of the Poisson equation with variable coefficients and embedded interfaces
Konrad-Zuse-Zentrum für Informationstechnik Berlin Takustraße 7 D-495 Berlin-Dahlem Germany MICHAEL OEVERMANN, RUPERT KLEIN A cartesian grid finite volume method for the solution of the Poisson equation
More informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationFinite volume method for nonlinear transmission problems
Finite volume method for nonlinear transmission problems Franck Boyer, Florence Hubert To cite this version: Franck Boyer, Florence Hubert. Finite volume method for nonlinear transmission problems. Proceedings
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 24: Preconditioning and Multigrid Solver Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 5 Preconditioning Motivation:
More informationHigh-resolution finite volume methods for hyperbolic PDEs on manifolds
High-resolution finite volume methods for hyperbolic PDEs on manifolds Randall J. LeVeque Department of Applied Mathematics University of Washington Supported in part by NSF, DOE Overview High-resolution
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 informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationFinite Difference Solution of the Heat Equation
Finite Difference Solution of the Heat Equation Adam Powell 22.091 March 13 15, 2002 In example 4.3 (p. 10) of his lecture notes for March 11, Rodolfo Rosales gives the constant-density heat equation as:
More informationEfficient Solvers for Stochastic Finite Element Saddle Point Problems
Efficient Solvers for Stochastic Finite Element Saddle Point Problems Catherine E. Powell c.powell@manchester.ac.uk School of Mathematics University of Manchester, UK Efficient Solvers for Stochastic Finite
More informationSolving the Stochastic Steady-State Diffusion Problem Using Multigrid
Solving the Stochastic Steady-State Diffusion Problem Using Multigrid Tengfei Su Applied Mathematics and Scientific Computing Advisor: Howard Elman Department of Computer Science Sept. 29, 2015 Tengfei
More informationPreconditioning for Nonsymmetry and Time-dependence
Preconditioning for Nonsymmetry and Time-dependence Andy Wathen Oxford University, UK joint work with Jen Pestana and Elle McDonald Jeju, Korea, 2015 p.1/24 Iterative methods For self-adjoint problems/symmetric
More informationTwo-parameter regularization method for determining the heat source
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 13, Number 8 (017), pp. 3937-3950 Research India Publications http://www.ripublication.com Two-parameter regularization method for
More informationMotion Estimation (I) Ce Liu Microsoft Research New England
Motion Estimation (I) Ce Liu celiu@microsoft.com Microsoft Research New England We live in a moving world Perceiving, understanding and predicting motion is an important part of our daily lives Motion
More informationIntroduction to Multigrid Methods
Introduction to Multigrid Methods Chapter 9: Multigrid Methodology and Applications Gustaf Söderlind Numerical Analysis, Lund University Textbooks: A Multigrid Tutorial, by William L Briggs. SIAM 1988
More information