Parameter-Uniform Numerical Methods for a Class of Singularly Perturbed Problems with a Neumann Boundary Condition
|
|
- Norma Gilmore
- 5 years ago
- Views:
Transcription
1 Parameter-Uniform Numerical Methods for a Class of Singularly Perturbed Problems with a Neumann Boundary Condition P. A. Farrell 1,A.F.Hegarty 2, J. J. H. Miller 3, E. O Riordan 4,andG.I. Shishkin 5 1 Department of Mathematics and Computer Science, Kent State University Kent, Ohio 44242, USA 2 Department of Mathematics and Statistics, University of Limerick, Ireland 3 Department of Mathematics, Trinity College, Dublin, Ireland 4 School of Mathematical Sciences, Dublin City University, Ireland 5 Institute for Mathematics and Mechanics, Russian Academy of Sciences Ekaterinburg, Russia Abstract. The error generated by the classical upwind finite difference method on a uniform mesh, when applied to a class of singularly perturbed modelordinary differentialequations with a singularly perturbed Neumann boundary condition, tends to infinity as the singular perturbation parameter tends to zero. Note that the exact solution is uniformly bounded with respect to the perturbation parameter. For the same classicalfinite difference operator on an appropriate piecewise uniform mesh, it is shown that the numerical solutions converge, uniformly with respect to the perturbation parameter, to the exact solution of any problem from this class. 1 Introduction Consider the following class of linear one dimensional convection diffusion problems L ε u ε εu ε + a(xu ε = f(x, x Ω =(, 1, (1a εu ε ( = A, u ε(1 = B, (1b a, f C 2 (Ω, a(x α>, x Ω. (1c Note that a Neumann boundary condition has been specified at x =. We recall the comparison principle for this problem (see [2], for example. Theorem 1. Assume that v C 2 (Ω. Then, if v (, v(1 and L ε v(x for all x Ω, it follows that v(x for all x Ω. This research was supported in part by the Russian Foundation for Basic Research under grant No , by the NationalScience Foundation grant DMS and by the Enterprise Ireland grant SC L. Vulkov, J. Waśniewski, and P. Yalamov (Eds.: NAA 2, LNCS 1988, pp , 21. c Springer-Verlag Berlin Heidelberg 21
2 Parameter-Uniform NumericalMethods 293 From this we can easily establish the following stability bound on the solution u ε (x u ε (1 + ε α u ε ( e αx/ε + 1 f (1 x. α Lemma 1. [2] The derivatives u (k ε of the solution of (1 satisfy the bounds u (k u (3 ε where C depends only on a and a. ε Cε k max{ f, u ε }, k =1, 2 Cε 3 max{ f, f, u ε } Consider the following decomposition of the solution u ε u ε = v ε + w ε,v ε = v + εv 1 + ε 2 v 2 ; w ε = w + εw 1 (2a where the components v,v 1 and v 2 are the solutions of the problems av = f, v (1 = u ε (1, (2b av 1 = v, v 1(1 =, (2c L ε v 2 = v 1, εv 2( =, v 2 (1 = (2d and the components w,w 1 are the solutions of L ε w =, εw ( = εu ε(, w (1 =. (2e L ε w 1 =, εw 1( = v ε(, w 1 (1 =. (2f Thus the components v ε and w ε are the solutions of the problems L ε v ε = f, v ε( = v ( + εv 1(, v ε (1 = u ε (1 L ε w ε =, w ε ( = u ε ( v ε (, w ε(1 =. Also, they satisfy the bounds given in the following lemma. Lemma 2. The components v ε,w ε and their derivatives satisfy v (k ε C(1 + ε 2 k, k =, 1, 2, 3 w (k ε (x C(ε u ε( + εε k e αx/ε, k =, 1, 2, 3. Proof. UseLemma1andthefactthat εψ ε (x = is the exact solution of 1 x e A(t/ε dt, where A(t = t L ε ψ ε =, εψ ε( = 1, ψ ε (1 =. a(sds
3 294 P. A. Farrell et al. 2 Upwinding on a Uniform Mesh In this section we examine the convergence behaviour of standard upwinding on a uniform mesh. The Neumann boundary condition is discretized by the scaled discrete derivative εd + U ε (. L N ε U ε εδ 2 U ε + a(x i D + U ε = f(x i, x i Ω N, (3a εd + U ε ( = εu ε (, U ε(1 = u ε (1, (3b where Ω N is an arbitrary mesh. We nowstate a discrete comparison principle. Theorem 2. [2] LetL N ε be the upwind finite difference operator defined in (3 and let Ω N be an arbitrary mesh of N +1 mesh points. If V is any mesh function defined on this mesh such that D + V (x, V(x N and L N ε V in Ω N, then V (x i for all x i Ω N. Hence, on an arbitrary mesh, the discrete solution U ε satisfies the bound U ε (x i u ε (1 + ε u ε( Φ i + 1 α f (1 x i. where Φ i is the solution of the constant coefficient problem εδ 2 Φ i + αd + Φ i =, εd + Φ = 1, Φ N =. Theorem 3. Let u ε be the continuous solution of (1 andletu ε be the numerical solution generated from the upwind finite difference scheme (3 on a uniform mesh Ωu N.Then, (a if u ε ( C, we have U ε u ε Ω CN 1 where U ε is the linear interpolant of U ε and C is a constant independent of N and ε. Also, (b if ε u ε( = C, then for any fixed N, U ε as ε Proof. (aconsiderfirstthecaseof u ε( C. The discrete solution U ε can be decomposed into the sum U ε = V ε + W ε where V ε and W ε are respectively the solutions of the problems L N ε V ε = f(x i, x i Ωu N, εd + V ε ( = εv ε(, V ε (1 = v ε (1 L N ε W ε =,x i Ωu N, εd+ W ε ( = εw ε (, W ε(1 =.
4 Parameter-Uniform NumericalMethods 295 We estimate the errors V ε v ε and W ε w ε separately. By standard local truncation error estimates, we obtain L N ε (V ε v ε (x i ε 3 (x i+1 x i 1 v ε (3 + a(x i 2 (x i+1 x i v ε (2 CN 1. Note also that D + (V ε v ε ( = v ε( D + v ε ( = 1 h h (s hv ε (s ds CN 1. With the two functions ψ ± (x i =CN 1 (1 x i ± (V ε v ε (x i, and the discrete minimum principle for L N ε we easily derive Note that if u ε ( C then (V ε v ε (x i CN 1. w (k ε (x Cε 1 k e αx/ε, k =, 1, 2, 3. The local truncation error for the layer component is given by and L N ε (W ε w ε (x i Cε 1 (x i+1 x i 1 e αxi 1/ε Cε 1 N 1 e αxi 1/ε εd + (W ε w ε ( = ε h h (s hw ε (s ds CN 1. Introduce the two mesh functions Ψ ± i = Cλ2 γ(α γ N 1 Y i ± (W ε w ε (x i where γ is any constant satisfying <γ<αand Y i = λn i 1 λ N 1, λ =1+γh ε, h =1/N. Note that Y i λ i. It is easy to see that λ 2 D + Y i γ ε e γxi 1/ε and so Y i decreases monotonically with Y i 1. We then have εd + Ψ ±, Ψ ± N =andusing( εδ 2 + γd + Y i =,weobtain L N ε Ψ ± i = Cλ2 γ(α γ N 1 (a(x i γd + Y i ± L N ε (W ε w ε (x i Cε 1 N 1( a(x i γ α γ e γxi 1/ε e αxi 1/ε <.
5 296 P. A. Farrell et al. By the discrete minimum principle we conclude that Ψ ± i andsoforallx i Ω N u ( W ε w ε (xi Cλ2 γ(α γ N 1 Y i Cλ 2 N 1. Thus, we have that ( W ε w ε (xi CN 1, when h ε. From an integral representation of the truncation error, we have xi+1 L N 1 ε (W ε w ε (x i C x i 1 ε e αt/ε dt Ce αxi 1/ε. As before we can establish Hence, for 1 i N and ε h, Note that which implies that ( W ε w ε (xi Cελ 2 i. ( W ε w ε (xi Cελ Ch. D + (W ε w ε ( = 1 h h (s hw (s ds C ( W ε w ε ( ( Wε w ε (x1 + Ch Ch. On the interval [x i,x i+1 ]wehave x (w ε w ε (x = w ε (tdt x x i w ε x i x i+1 x (tdt CN 1 w ε i x i Combining this with the argument in [2] completes part (a. (b If we discretize the Neumann boundary condition εu ε( = C by the standard discrete derivative εd + U ε ( = C on a uniform mesh then For a fixed distance h, ε xi+1 U ε (h =U ε ( + C h ε. lim U ε(h U ε (. ε On a uniform mesh, the discrete solution is not bounded independently of ε.
6 Parameter-Uniform NumericalMethods 297 Remarks. (i We define a a weak boundary layer by u ε (x C, thatis,the derivative is uniformly bounded with respect to ε. In this case, Theorem 3 states that the solution u ε can be approximated ε-uniformly on a uniform mesh. However, the first derivative u ε (x is not approximated ε-uniformly by the discrete derivative D + U ε (x i on a uniform mesh. This can be checked by solving a nontrivial constant coefficient continuous problem and its corresponding discrete problem directly, setting εn = 1 and then taking the limit as N. (ii In the case of the constant coefficient problem (1, we observe that lim ε εd+ U ε (x i εu ε (x i =. Thus, although the discrete solutions are unbounded as ε, the scaleddiscrete derivatives are at least ε uniformly bounded and, moreover, converge as N to εu (x i for each fixed ε. However, the scaled discrete derivatives εd + U ε ( are not ε uniformly convergent to εu ε (. In contrast, for the problem L ε u ε εu ε + a(xu ε = f(x, x Ω, (4a u ε ( = A, u ε (1 = B, (4b with Dirichelet boundary conditions, we have that lim U ε(x i u ε (x i =, ε and that lim ε εd+ U ε ( εu ε ( = O(1. This can be seen easily from the explicit solutions to the constant coefficient continuous problem. u ε (x =u ε ( + f α x (u ε( u ε (1 + f ( 1 e αx/ε α and the discrete problem 1 e α/ε U ε (x i =u ε ( + f α x i (u ε ( u ε (1 + f α ( 1 λ i 1 λ N,λ=1+αh/ε. 3 Upwinding on a Piecewise Uniform Mesh Consider the same upwind finite difference scheme (3 on the piecewise uniform mesh Ω N ε = {x i x i =2iσ/N, i N/2; x i = x i 1 +2(1 σ/n, N/2 <i} (5a where the transition parameter σ is fitted to the boundary layer by taking σ =min{ 1 2, 1 ε ln N}. α (5b The next result shows that upwinding on this mesh produces an ε uniform numerical method.
7 298 P. A. Farrell et al. Theorem 4. Let u ε be the continuous solution of (1 andletu ε be the numerical solution generated from an upwind finite difference scheme (3 on the piecewiseuniform mesh (5. Then, for all N 4, we have U ε u ε N Ω CN 1 ln N ε where C is a constant independent of N and ε. Proof. As for the uniform mesh we derive (V ε v ε (x i CN 1. When σ =1/2, the mesh is uniform and applying the argument of the previous theorem, we get ( W ε w ε (xi CN 1 ε 1 CN 1 ln N. When σ<1/2, the argument is divided between the coarse mesh and fine mesh regions. Consider first the coarse mesh region [σ, 1], where w ε (x Ce ασ/ε CN 1. Using the discrete comparison principle, we get W ε (x i ε w ε( Φ i where Φ i is the solution of the constant coefficient problem εδ 2 Φ i + αd + Φ i =, εd + Φ = 1, Φ N =. From an explicit representation of Φ i one can showthat Φ N/2 CN 1 Hence, for x i σ W ε (x i w ε (x i W ε (x i + w ε (x i CN 1 Consider nowthe fine mesh region, using the same argument as in the previous theorem we get ( W ε w ε (xi Cλ 2 N 1 ln N CN 1 ln N. This completes the proof. In [1], an essentially second order scheme is constructed on a piecewiseuniform mesh, using a more complicated finite difference operator. As in [2], the nodal error estimate for the simpler scheme presented here can easily be extended to a global error estimate by simple linear interpolation. That is, we have Ūε u ε Ω CN 1 ln N where Ūε is the linear interpolant of U ε. Also, using the techniques in [2], one can deduce that ε D + U ε u ε Ω\{1} CN 1 ln N.
8 4 Parabolic Boundary Layers Parameter-Uniform NumericalMethods 299 In this section, we introduce a new class of problems. Let Ω = (, 1,D = Ω (,T]andΓ = Γ l Γ b Γ r where Γ l and Γ r are the left and right sides of the box D and Γ b is its base. Consider the following linear parabolic partial differential equation in D with Dirichlet-Neumann boundary conditions on Γ L ε u ε (x, t ε 2 u ε x 2 + b(x, tu ε + d(x, t u ε t = f(x, t, (x, t D, (6a u ε = ϕ b on Γ b, ε u ε x = ϕ l on Γ l, u ε = ϕ r on Γ r (6b d(x, t >δ> and b(x, t β>, (x, t D (6c ϕ l ( = εϕ b(, ϕ b (1 = ϕ r (. (6d We have the comparison principle Lemma 3. Assume b, d C (D and ψ C 2 (D C 1 (D. Suppose that ψ on Γ b Γ r and ψ x on Γ r.thenl ε ψ in D implies that ψ in D. and the following stability bound Theorem 5. Let v be any function in the domain of the differential operator L ε. Then ε v (1 + αt max{ L ε v, v Γb Γ r } + β v x Γl e βx/ ε where α =max D {, (1 b/d} 1/δ. Assume that the data b, d, f, ϕ satisfy sufficient regularity and compatibility conditions so that the problem has a unique solution u ε and u ε Cλ 4 (D and, furthermore, such that the derivatives of the solution u ε satisfy, for all nonnegative integers i, j, i +2j 4 i+j u ε x i t j Cε i/2 D where the constant C is independent of ε. We write the solution as the sum u ε = v ε + w ε where v ε,w ε are smooth and singular components of u ε defined in the following way. The smooth component is further decomposed into the sum where v,v 1 are defined by v ε = v + εv 1 bv + d v t = f in D, v = u ε on Γ b (7a L ε v 1 = 2 v x 2 in D, v 1 =onγ \ Γ l, v 1 x =onγ l. (7b
9 3 P. A. Farrell et al. The singular component is decomposed into the sum where w l and w r are defined by w ε = w l + w r L ε w r =ind, w r = u ε v on Γ r,w r =onγ b Γ l (8a L ε w l =ind, (8b w l x = u ε x v x w r x on Γ l,w l =onγ r Γ b. (8c It is clear that w l,w r correspond respectively to the boundary layer functions on Γ l and Γ r. Assume that the data satisfy sufficient regularity and compatibility conditions so that v ε,w ε Cλ 4(D. Theorem 6. [3] For all non-negative integers i, j, such that i +2j 4 i+j v ε x i t j C(1 + ε 1 i/2 D and for all (x, t D, i+j w l (x, t x i t j Cε i/2 e x/ ε, i+j w r (x, t x i t j Cε i/2 e (1 x/ ε where C is a constant independent of ε. Problem (6 is discretized using a standard numerical method composed of a standard finite difference operator on a fitted piecewise uniform mesh. L N ε U ε = εδxu 2 ε + bu ε + ddt U ε = f,(x, t Dσ N U ε = u ε on Γb,σ N Γ r,σ N, D+ x U ε = u ε x on Γ l,σ N (9a (9b where Dσ N = ΩNx σ Ω Nt, and Γσ N = DN σ Γ. (9c A uniform mesh Ω Nt with N t mesh elements is used on (,T. A piecewise uniform mesh Ωσ Nx on Ω with N x mesh elements is obtained by putting a uniform mesh with N x /4meshelementsonboth(,σand(1 σ, 1 and one with N x /2 mesh elements on (σ, 1 σ, with the transition parameter { 1 σ =min 4, 2 } ε ln N x. (9d We have the following discrete comparison principle Lemma 4. Assume that the mesh function Ψ satisfies Ψ on Γ N b,σ Γ N r,σ and D + x Ψ on Γ N l,σ.thenln ε Ψ on DN σ implies that Ψ on DN σ.
10 Parameter-Uniform NumericalMethods 31 The ε uniform error estimate is contained in Theorem 7. Let u ε be the continuous solution of (6 andletu ε be the numerical solution generated from (9. Assume that v ε,w ε Cλ 4 (D. Then,forallN 4, we have sup U ε u ε N D CNx 1 ln N x + CNt 1 <ε 1 σ where C is a constant independent of N x,n t and ε. Proof. The argument follows [3]. The discrete solution U ε is the sum U ε = V ε + W ε where V ε and W ε are the obvious discrete counterparts to v ε and w ε.the classical truncation error estimate yields L N ε (V ε v ε C εn 1 x It follows that V ε v ε CN 1 x + CN 1 t and D + x (V ε v ε (,t CN 1 x. + CNt 1. Note also that L N ε (W l w l CNx 1 ln N x + CNt 1 and ε D + x (W l w l (,t CNx 1 ln N x. The proof is completed as in [3]. 5 Numerical Results In this section we present numerical results for the following specific elliptic problem ε u ε x =, (x, y Γ L, ε u ε + u ε x =16x(1 x(1 yy, (x, y (, 12 (1a u ε =1, (x, y Γ R Γ T (1b ε u ε y = 16x2 (1 x 2, (x, y Γ B (1c whose solution has a parabolic boundary layer near Γ B. The nature of the boundary layer function associated with this layer is related to the solutions of the parabolic problems examined in the previous section. In Figure 1 we present the numerical solution generated by applying standard upwinding on a uniform mesh. The numerical solutions are not bounded uniformly with respect to ε as ε. This should be compared with the accurate approximation given in Figure 2, which was generated by applying standard upwinding on the piecewise uniform mesh Ωσ N,N 2 = Ωu N Ωτ N,
11 32 P. A. Farrell et al x y Fig. 1. Numerical solution generated by upwinding on a uniform mesh with N=32, ε=1 8 for problem (1 where { 1 τ =min 2, } ε ln N. Note the significant difference in the vertical scale in these two figures. In Table 1 we present the computed orders of convergence (see [2] generated by applying standard upwinding on this piecewise uniform mesh. These indicate that the method is ε uniformly convergent for problem (1. References 1. Andreyev V. B. and Savin I. A. (1996. The computation of boundary flow with uniform accuracy with respect to a small parameter. Comput. Maths. Math. Phys., 36 ( Farrell, P. A., Hegarty A.F, Miller, J. J. H., O Riordan, E.,Shishkin G. I., Robust Computational techniques for boundary layers, Chapman and Hall/CRC Press, , 293, 294, 296, 298, Miller, J. J. H.,O Riordan, E.,Shishkin G. I.and Shishkina,L. P. (1998 Fitted mesh methods for problems with parabolic boundary layers, Math. Proc. RoyalIrish Academy, 98A ( , 31
12 Parameter-Uniform NumericalMethods x y Fig. 2. Numerical solution generated by upwinding on a piecewise uniform mesh with N=32, ε=1 8 for problem (1 Table 1. Computed orders of convergence generated by upwinding on a piecewise uniform mesh applied to problem (1 Number of intervals N ε p N
A Singularly Perturbed Convection Diffusion Turning Point Problem with an Interior Layer
A Singularly Perturbed Convection Diffusion Turning Point Problem with an Interior Layer E. O Riordan, J. Quinn School of Mathematical Sciences, Dublin City University, Glasnevin, Dublin 9, Ireland. Abstract
More informationA PARAMETER ROBUST NUMERICAL METHOD FOR A TWO DIMENSIONAL REACTION-DIFFUSION PROBLEM
A PARAMETER ROBUST NUMERICAL METHOD FOR A TWO DIMENSIONAL REACTION-DIFFUSION PROBLEM C. CLAVERO, J.L. GRACIA, AND E. O RIORDAN Abstract. In this paper a singularly perturbed reaction diffusion partial
More informationInterior Layers in Singularly Perturbed Problems
Interior Layers in Singularly Perturbed Problems Eugene O Riordan Abstract To construct layer adapted meshes for a class of singularly perturbed problems, whose solutions contain boundary layers, it is
More informationIterative Domain Decomposition Methods for Singularly Perturbed Nonlinear Convection-Diffusion Equations
Iterative Domain Decomposition Methods for Singularly Perturbed Nonlinear Convection-Diffusion Equations P.A. Farrell 1, P.W. Hemker 2, G.I. Shishkin 3 and L.P. Shishkina 3 1 Department of Computer Science,
More informationParameter robust methods for second-order complex-valued reaction-diffusion equations
Postgraduate Modelling Research Group, 10 November 2017 Parameter robust methods for second-order complex-valued reaction-diffusion equations Faiza Alssaedi Supervisor: Niall Madden School of Mathematics,
More informationACCURATE NUMERICAL METHOD FOR BLASIUS PROBLEM FOR FLOW PAST A FLAT PLATE WITH MASS TRANSFER
FDS-2000 Conference, pp. 1 10 M. Sapagovas et al. (Eds) c 2000 MII ACCURATE UMERICAL METHOD FOR BLASIUS PROBLEM FOR FLOW PAST A FLAT PLATE WITH MASS TRASFER B. GAHA 1, J. J. H. MILLER 2, and G. I. SHISHKI
More informationA uniformly convergent numerical method for a coupled system of two singularly perturbed linear reaction diffusion problems
IMA Journal of Numerical Analysis (2003 23, 627 644 A uniformly convergent numerical method for a coupled system of two singularly perturbed linear reaction diffusion problems NIALL MADDEN Department of
More informationA Linearised Singularly Perturbed Convection- Diffusion Problem with an Interior Layer
Dublin Institute of Technology ARROW@DIT Articles School of Mathematics 201 A Linearised Singularly Perturbed Convection- Diffusion Problem with an Interior Layer Eugene O'Riordan Dublin City University,
More informationA coupled system of singularly perturbed semilinear reaction-diffusion equations
A coupled system of singularly perturbed semilinear reaction-diffusion equations J.L. Gracia, F.J. Lisbona, M. Madaune-Tort E. O Riordan September 9, 008 Abstract In this paper singularly perturbed semilinear
More informationOn the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations
On the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations Lutz Tobiska Institut für Analysis und Numerik Otto-von-Guericke-Universität
More informationRemarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable?
Remarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable? Thomas Apel, Hans-G. Roos 22.7.2008 Abstract In the first part of the paper we discuss minimal
More informationSuperconvergence analysis of the SDFEM for elliptic problems with characteristic layers
Superconvergence analysis of the SDFEM for elliptic problems with characteristic layers S. Franz, T. Linß, H.-G. Roos Institut für Numerische Mathematik, Technische Universität Dresden, D-01062, Germany
More informationA COMPUTATIONAL TECHNIQUE FOR SOLVING BOUNDARY VALUE PROBLEM WITH TWO SMALL PARAMETERS
Electronic Journal of Differential Equations, Vol. 2013 (2013), No. 30, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu A COMPUTATIONAL
More informationAn exponentially graded mesh for singularly perturbed problems
An exponentially graded mesh for singularly perturbed problems Christos Xenophontos Department of Mathematics and Statistics University of Cyprus joint work with P. Constantinou (UCY), S. Franz (TU Dresden)
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 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 Element Superconvergence Approximation for One-Dimensional Singularly Perturbed Problems
Finite Element Superconvergence Approximation for One-Dimensional Singularly Perturbed Problems Zhimin Zhang Department of Mathematics Wayne State University Detroit, MI 48202 Received 19 July 2000; accepted
More informationOn solving linear systems arising from Shishkin mesh discretizations
On solving linear systems arising from Shishkin mesh discretizations Petr Tichý Faculty of Mathematics and Physics, Charles University joint work with Carlos Echeverría, Jörg Liesen, and Daniel Szyld October
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 information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
More informationA New Improvement of Shishkin Fitted Mesh Technique on a First Order Finite Difference Method
A New Improvement of Shishkin Fitted Mesh Technique on a First Order Finite Difference Method With applications on singularly perturbed boundary value problem ODE Rafiq S. Muhammad College of Computers
More informationOn a Two-Point Boundary-Value Problem with Spurious Solutions
On a Two-Point Boundary-Value Problem with Spurious Solutions By C. H. Ou and R. Wong The Carrier Pearson equation εü + u = with boundary conditions u ) = u) = 0 is studied from a rigorous point of view.
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationThe method of lines (MOL) for the diffusion equation
Chapter 1 The method of lines (MOL) for the diffusion equation The method of lines refers to an approximation of one or more partial differential equations with ordinary differential equations in just
More informationON SOME ELLIPTIC PROBLEMS IN UNBOUNDED DOMAINS
Chin. Ann. Math.??B(?), 200?, 1 20 DOI: 10.1007/s11401-007-0001-x ON SOME ELLIPTIC PROBLEMS IN UNBOUNDED DOMAINS Michel CHIPOT Abstract We present a method allowing to obtain existence of a solution for
More informationSPLINE COLLOCATION METHOD FOR SINGULAR PERTURBATION PROBLEM. Mirjana Stojanović University of Novi Sad, Yugoslavia
GLASNIK MATEMATIČKI Vol. 37(57(2002, 393 403 SPLINE COLLOCATION METHOD FOR SINGULAR PERTURBATION PROBLEM Mirjana Stojanović University of Novi Sad, Yugoslavia Abstract. We introduce piecewise interpolating
More informationOn Multigrid for Phase Field
On Multigrid for Phase Field Carsten Gräser (FU Berlin), Ralf Kornhuber (FU Berlin), Rolf Krause (Uni Bonn), and Vanessa Styles (University of Sussex) Interphase 04 Rome, September, 13-16, 2004 Synopsis
More informationANISOTROPIC EQUATIONS: UNIQUENESS AND EXISTENCE RESULTS
ANISOTROPIC EQUATIONS: UNIQUENESS AND EXISTENCE RESULTS STANISLAV ANTONTSEV, MICHEL CHIPOT Abstract. We study uniqueness of weak solutions for elliptic equations of the following type ) xi (a i (x, u)
More informationApproximation of interior-layer. solutions. in a singularly perturbed semilinear. reaction-diffusion problem
Numerische Mathematik manuscript No. (will be inserted by the editor) Natalia Kopteva Martin Stynes Approximation of interior-layer solutions in a singularly perturbed semilinear reaction-diffusion problem
More informationCentrum voor Wiskunde en Informatica
Centrum voor Wiskunde en Informatica Modelling, Analysis and Simulation Modelling, Analysis and Simulation High-order time-accurate schemes for singularly perturbed parabolic convection-diffusion problems
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 informationSUPERCONVERGENCE OF DG METHOD FOR ONE-DIMENSIONAL SINGULARLY PERTURBED PROBLEMS *1)
Journal of Computational Mathematics, Vol.5, No., 007, 185 00. SUPERCONVERGENCE OF DG METHOD FOR ONE-DIMENSIONAL SINGULARLY PERTURBED PROBLEMS *1) Ziqing Xie (College of Mathematics and Computer Science,
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
More informationarxiv:math/ v1 [math.ap] 28 Oct 2005
arxiv:math/050643v [math.ap] 28 Oct 2005 A remark on asymptotic completeness for the critical nonlinear Klein-Gordon equation Hans Lindblad and Avy Soffer University of California at San Diego and Rutgers
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationON SINGULAR PERTURBATION OF THE STOKES PROBLEM
NUMERICAL ANALYSIS AND MATHEMATICAL MODELLING BANACH CENTER PUBLICATIONS, VOLUME 9 INSTITUTE OF MATHEMATICS POLISH ACADEMY OF SCIENCES WARSZAWA 994 ON SINGULAR PERTURBATION OF THE STOKES PROBLEM G. M.
More informationBoundary layers in a two-point boundary value problem with fractional derivatives
Boundary layers in a two-point boundary value problem with fractional derivatives J.L. Gracia and M. Stynes Institute of Mathematics and Applications (IUMA) and Department of Applied Mathematics, University
More informationSplitting methods with boundary corrections
Splitting methods with boundary corrections Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Strang s paper, SIAM J. Numer. Anal., 1968 S (5)
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationEverywhere differentiability of infinity harmonic functions
Everywhere differentiability of infinity harmonic functions Lawrence C. Evans and Charles K. Smart Department of Mathematics University of California, Berkeley Abstract We show that an infinity harmonic
More informationRenormalized Solutions of a Nonlinear Parabolic Equation with Double Degeneracy
Electronic Journal of Qualitative Theory of Differential Equations 26, No. 5, -2; http://www.math.u-szeged.hu/ejqtde/ Renormalized Solutions of a Nonlinear Parabolic Equation with Double Degeneracy Zejia
More informationTHE MULTIPLICATIVE ERGODIC THEOREM OF OSELEDETS
THE MULTIPLICATIVE ERGODIC THEOREM OF OSELEDETS. STATEMENT Let (X, µ, A) be a probability space, and let T : X X be an ergodic measure-preserving transformation. Given a measurable map A : X GL(d, R),
More informationUNIFORMLY CONVERGENT 3-TGFEM VS LSFEM FOR SINGULARLY PERTURBED CONVECTION-DIFFUSION PROBLEMS ON A SHISHKIN BASED LOGARITHMIC MESH
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING, SERIES B Volume 4, Number 4, Pages 35 334 c 23 Institute for Scientific Computing and Information UNIFORMLY CONVERGENT 3-TGFEM VS LSFEM FOR SINGULARLY
More informationMeasure and Integration: Solutions of CW2
Measure and Integration: s of CW2 Fall 206 [G. Holzegel] December 9, 206 Problem of Sheet 5 a) Left (f n ) and (g n ) be sequences of integrable functions with f n (x) f (x) and g n (x) g (x) for almost
More informationNONLOCAL DIFFUSION EQUATIONS
NONLOCAL DIFFUSION EQUATIONS JULIO D. ROSSI (ALICANTE, SPAIN AND BUENOS AIRES, ARGENTINA) jrossi@dm.uba.ar http://mate.dm.uba.ar/ jrossi 2011 Non-local diffusion. The function J. Let J : R N R, nonnegative,
More informationLECTURE 3: DISCRETE GRADIENT FLOWS
LECTURE 3: DISCRETE GRADIENT FLOWS Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA Tutorial: Numerical Methods for FBPs Free Boundary Problems and
More informationSYMMETRY RESULTS FOR PERTURBED PROBLEMS AND RELATED QUESTIONS. Massimo Grosi Filomena Pacella S. L. Yadava. 1. Introduction
Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 21, 2003, 211 226 SYMMETRY RESULTS FOR PERTURBED PROBLEMS AND RELATED QUESTIONS Massimo Grosi Filomena Pacella S.
More informationA first-order system Petrov-Galerkin discretisation for a reaction-diffusion problem on a fitted mesh
A first-order system Petrov-Galerkin discretisation for a reaction-diffusion problem on a fitted mesh James Adler, Department of Mathematics Tufts University Medford, MA 02155 Scott MacLachlan Department
More informationON A DIFFERENCE EQUATION WITH MIN-MAX RESPONSE
IJMMS 2004:55, 295 2926 PII. S067204402270 http://ijmms.hindawi.com Hindawi Publishing Corp. ON A DIFFERENCE EQUATION WITH MIN-MAX RESPONSE GEORGE L. KARAKOSTAS and STEVO STEVIĆ Received 25 February 2004
More informationResearch Article Almost Periodic Viscosity Solutions of Nonlinear Parabolic Equations
Hindawi Publishing Corporation Boundary Value Problems Volume 29, Article ID 873526, 15 pages doi:1.1155/29/873526 Research Article Almost Periodic Viscosity Solutions of Nonlinear Parabolic Equations
More informationA Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 41, pp. 376-395, 2014. Copyright 2014,. ISSN 1068-9613. A ROBUST NUMERICAL SCHEME FOR SINGULARLY PERTURBED DELAY PARABOLIC INITIAL-BOUNDARY-VALUE PROBLEMS
More informationExplosive Solution of the Nonlinear Equation of a Parabolic Type
Int. J. Contemp. Math. Sciences, Vol. 4, 2009, no. 5, 233-239 Explosive Solution of the Nonlinear Equation of a Parabolic Type T. S. Hajiev Institute of Mathematics and Mechanics, Acad. of Sciences Baku,
More informationA Ginzburg-Landau approach to dislocations. Marcello Ponsiglione Sapienza Università di Roma
Marcello Ponsiglione Sapienza Università di Roma Description of a dislocation line A deformed crystal C can be described by The displacement function u : C R 3. The strain function β = u. A dislocation
More informationGoal. Robust A Posteriori Error Estimates for Stabilized Finite Element Discretizations of Non-Stationary Convection-Diffusion Problems.
Robust A Posteriori Error Estimates for Stabilized Finite Element s of Non-Stationary Convection-Diffusion Problems L. Tobiska and R. Verfürth Universität Magdeburg Ruhr-Universität Bochum www.ruhr-uni-bochum.de/num
More informationON A PRIORI ERROR ANALYSIS OF FULLY DISCRETE HETEROGENEOUS MULTISCALE FEM
MULTISCALE MODEL. SIMUL. Vol. 4, No. 2, pp. 447 459 c 2005 Society for Industrial and Applied Mathematics ON A PRIORI ERROR ANALYSIS OF FULLY DISCRETE HETEROGENEOUS MULTISCALE FEM ASSYR ABDULLE Abstract.
More informationRobust exponential convergence of hp-fem for singularly perturbed systems of reaction-diffusion equations
Robust exponential convergence of hp-fem for singularly perturbed systems of reaction-diffusion equations Christos Xenophontos Department of Mathematics and Statistics University of Cyprus joint work with
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationA uniformly convergent difference scheme on a modified Shishkin mesh for the singularly perturbed reaction-diffusion boundary value problem
Journal of Modern Methods in Numerical Mathematics 6:1 (2015, 28 43 A uniformly convergent difference scheme on a modified Shishkin mesh for the singularly perturbed reaction-diffusion boundary value problem
More informationComplex geometrical optics solutions for Lipschitz conductivities
Rev. Mat. Iberoamericana 19 (2003), 57 72 Complex geometrical optics solutions for Lipschitz conductivities Lassi Päivärinta, Alexander Panchenko and Gunther Uhlmann Abstract We prove the existence of
More informationNUMERICAL ANALYSIS OF THE CAHN-HILLIARD EQUATION AND APPROXIMATION FOR THE HELE-SHAW PROBLEM, PART II: ERROR ANALYSIS AND CONVERGENCE OF THE INTERFACE
NUMERICAL ANALYSIS OF THE CAHN-HILLIARD EQUATION AND APPROXIMATION FOR THE HELE-SHAW PROBLEM, PART II: ERROR ANALYSIS AND CONVERGENCE OF THE INTERFACE XIAOBING FENG AND ANDREAS PROHL Abstract. In this
More informationUNIFORMLY CONVERGENT APPROXIMATION ON SPECIAL MESHES
UNIFORMLY CONVERGENT APPROXIMATION ON SPECIAL MESHES A Thesis Submitted to the Graduate School of Engineering and Sciences of İzmir Institute of Technology in Partial Fulfillment of the Requirements for
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 informationDiscrete Maximum Principle for a 1D Problem with Piecewise-Constant Coefficients Solved by hp-fem
The University of Texas at El Paso Department of Mathematical Sciences Research Reports Series El Paso, Texas Research Report No. 2006-10 Discrete Maximum Principle for a 1D Problem with Piecewise-Constant
More informationCLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE
CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that
More informationThree remarks on anisotropic finite elements
Three remarks on anisotropic finite elements Thomas Apel Universität der Bundeswehr München Workshop Numerical Analysis for Singularly Perturbed Problems dedicated to the 60th birthday of Martin Stynes
More informationSpace-Time Nonconforming Optimized Schwarz Waveform Relaxation for Heterogeneous Problems and General Geometries
Space-Time Nonconforming Optimized Schwarz Waveform Relaxation for Heterogeneous Problems and General Geometries Laurence Halpern, Caroline Japhet, and Jérémie Szeftel 3 LAGA, Université Paris XIII, Villetaneuse
More informationNonlinear stabilization via a linear observability
via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result
More informationRelaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni
Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)
More informationPartial Differential Equations
Partial Differential Equations Xu Chen Assistant Professor United Technologies Engineering Build, Rm. 382 Department of Mechanical Engineering University of Connecticut xchen@engr.uconn.edu Contents 1
More informationSuggested Solution to Assignment 7
MATH 422 (25-6) partial diferential equations Suggested Solution to Assignment 7 Exercise 7.. Suppose there exists one non-constant harmonic function u in, which attains its maximum M at x. Then by the
More informationfor all subintervals I J. If the same is true for the dyadic subintervals I D J only, we will write ϕ BMO d (J). In fact, the following is true
3 ohn Nirenberg inequality, Part I A function ϕ L () belongs to the space BMO() if sup ϕ(s) ϕ I I I < for all subintervals I If the same is true for the dyadic subintervals I D only, we will write ϕ BMO
More informationNEWMAN S INEQUALITY FOR INCREASING EXPONENTIAL SUMS
NEWMAN S INEQUALITY FOR INCREASING EXPONENTIAL SUMS Tamás Erdélyi Dedicated to the memory of George G Lorentz Abstract Let Λ n := {λ 0 < λ < < λ n } be a set of real numbers The collection of all linear
More informationKonstantinos Chrysafinos 1 and L. Steven Hou Introduction
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher ANALYSIS AND APPROXIMATIONS OF THE EVOLUTIONARY STOKES EQUATIONS WITH INHOMOGENEOUS
More informationGreen s function estimates for a singularly perturbed convection-diffusion problem
Green s function estimates for a singularly perturbed convection-diffusion problem Sebastian Franz Natalia Kopteva Abstract We consider a singularly perturbed convection-diffusion problem posed in the
More informationParameter Dependent Quasi-Linear Parabolic Equations
CADERNOS DE MATEMÁTICA 4, 39 33 October (23) ARTIGO NÚMERO SMA#79 Parameter Dependent Quasi-Linear Parabolic Equations Cláudia Buttarello Gentile Departamento de Matemática, Universidade Federal de São
More informationEXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN EQUATION WITH DEGENERATE MOBILITY
Electronic Journal of Differential Equations, Vol. 216 216), No. 329, pp. 1 22. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN
More informationWe denote the space of distributions on Ω by D ( Ω) 2.
Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study
More informationWeak Convergence Methods for Energy Minimization
Weak Convergence Methods for Energy Minimization Bo Li Department of Mathematics University of California, San Diego E-mail: bli@math.ucsd.edu June 3, 2007 Introduction This compact set of notes present
More informationAnalysis of a Streamline-Diffusion Finite Element Method on Bakhvalov-Shishkin Mesh for Singularly Perturbed Problem
Numer. Math. Theor. Meth. Appl. Vol. 10, No. 1, pp. 44-64 doi: 10.408/nmtma.017.y1306 February 017 Analysis of a Streamline-Diffusion Finite Element Method on Bakhvalov-Shishkin Mesh for Singularly Perturbed
More informationThe Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition
The Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition Sukjung Hwang CMAC, Yonsei University Collaboration with M. Dindos and M. Mitrea The 1st Meeting of
More informationOn Some Mean Value Results for the Zeta-Function and a Divisor Problem
Filomat 3:8 (26), 235 2327 DOI.2298/FIL6835I Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat On Some Mean Value Results for the
More informationSingular Integrals. 1 Calderon-Zygmund decomposition
Singular Integrals Analysis III Calderon-Zygmund decomposition Let f be an integrable function f dx 0, f = g + b with g Cα almost everywhere, with b
More informationBoundary Layer Solutions to Singularly Perturbed Problems via the Implicit Function Theorem
Boundary Layer Solutions to Singularly Perturbed Problems via the Implicit Function Theorem Oleh Omel chenko, and Lutz Recke Department of Mathematics, Humboldt University of Berlin, Unter den Linden 6,
More informationMath 46, Applied Math (Spring 2008): Final
Math 46, Applied Math (Spring 2008): Final 3 hours, 80 points total, 9 questions, roughly in syllabus order (apart from short answers) 1. [16 points. Note part c, worth 7 points, is independent of the
More informationNONTRIVIAL SOLUTIONS TO INTEGRAL AND DIFFERENTIAL EQUATIONS
Fixed Point Theory, Volume 9, No. 1, 28, 3-16 http://www.math.ubbcluj.ro/ nodeacj/sfptcj.html NONTRIVIAL SOLUTIONS TO INTEGRAL AND DIFFERENTIAL EQUATIONS GIOVANNI ANELLO Department of Mathematics University
More informationAn Exponentially Fitted Non Symmetric Finite Difference Method for Singular Perturbation Problems
An Exponentially Fitted Non Symmetric Finite Difference Method for Singular Perturbation Problems GBSL SOUJANYA National Institute of Technology Department of Mathematics Warangal INDIA gbslsoujanya@gmail.com
More informationFitted operator finite difference scheme for a class of singularly perturbed differential difference turning point problems exhibiting interior layer
Fitted operator finite difference scheme for a class of singularly perturbed differential difference turning point problems exhibiting interior layer Dr. Pratima Rai Department of mathematics, University
More information1. Statement of the problem.
218 Нелинейные граничные задачи 18, 218-229 (2008 c 2008. M. A. Borodin THE STEFAN PROBLEM The Stefan problem in its classical statement is a mathematical model of the process of propagation of heat in
More information13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)
13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs) A prototypical problem we will discuss in detail is the 1D diffusion equation u t = Du xx < x < l, t > finite-length rod u(x,
More informationExistence, stability and instability for Einstein-scalar field Lichnerowicz equations by Emmanuel Hebey
Existence, stability and instability for Einstein-scalar field Lichnerowicz equations by Emmanuel Hebey Joint works with Olivier Druet and with Frank Pacard and Dan Pollack Two hours lectures IAS, October
More informationConvergence on Layer-adapted Meshes and Anisotropic Interpolation Error Estimates of Non-Standard Higher Order Finite Elements
Convergence on Layer-adapted Meshes and Anisotropic Interpolation Error Estimates of Non-Standard Higher Order Finite Elements Sebastian Franz a,1, Gunar Matthies b, a Department of Mathematics and Statistics,
More informationNUMERICAL SOLUTION OF GENERAL SINGULAR PERTURBATION BOUNDARY VALUE PROBLEMS BASED ON ADAPTIVE CUBIC SPLINE
TWMS Jour. Pure Appl. Math., V.3, N.1, 2012, pp.11-21 NUMERICAL SOLUTION OF GENERAL SINGULAR PERTURBATION BOUNDARY VALUE PROBLEMS BASED ON ADAPTIVE CUBIC SPLINE R. MOHAMMADI 1 Abstract. We use adaptive
More informationSome recent results on controllability of coupled parabolic systems: Towards a Kalman condition
Some recent results on controllability of coupled parabolic systems: Towards a Kalman condition F. Ammar Khodja Clermont-Ferrand, June 2011 GOAL: 1 Show the important differences between scalar and non
More informationResearch Article On Behavior of Solution of Degenerated Hyperbolic Equation
International Scholarly Research Network ISRN Applied Mathematics Volume 2012, Article ID 124936, 10 pages doi:10.5402/2012/124936 Research Article On Behavior of Solution of Degenerated Hyperbolic Equation
More informationInverse Transport Problems and Applications. II. Optical Tomography and Clear Layers. Guillaume Bal
Inverse Transport Problems and Applications II. Optical Tomography and Clear Layers Guillaume Bal Department of Applied Physics & Applied Mathematics Columbia University http://www.columbia.edu/ gb23 gb23@columbia.edu
More informationNitsche XFEM with Streamline Diffusion Stabilization for a Two Phase Mass Transport Problem
Nitsche XFEM with Streamline Diffusion Stabilization for a Two Phase Mass Transport Problem Christoph Lehrenfeld and Arnold Reusken Bericht Nr. 333 November 2011 Key words: transport problem, Nitsche method,
More informationFEM Convergence for PDEs with Point Sources in 2-D and 3-D
FEM Convergence for PDEs with Point Sources in -D and 3-D Kourosh M. Kalayeh 1, Jonathan S. Graf, and Matthias K. Gobbert 1 Department of Mechanical Engineering, University of Maryland, Baltimore County
More informationBoundary Value Problems and Iterative Methods for Linear Systems
Boundary Value Problems and Iterative Methods for Linear Systems 1. Equilibrium Problems 1.1. Abstract setting We want to find a displacement u V. Here V is a complete vector space with a norm v V. In
More information