Numerical solutions of Dirichlet problem for elliptic operator in divergence form with a right-hand side measure

  • Authors:

  • Nedžad Limić;Mladen Rogina

  • Affiliations:

  • Department of Mathematics, University of Zagreb, Bijenicka 30, P.O. Box 335, 10002 Zagreb, Croatia;Department of Mathematics, University of Zagreb, Bijenicka 30, P.O. Box 335, 10002 Zagreb, Croatia

  • Venue:
  • Journal of Computational and Applied Mathematics - Special Issue: Proceedings of the 10th international congress on computational and applied mathematics (ICCAM-2002)
  • Year:
  • 2004

Quantified Score

Hi-index 0.00

Visualization

Abstract

We consider a second-order elliptic operator A=A(x)=-Σi,j=1d ∂iaij(x)∂j + Σj=1d b'j(x)∂j + Σj=1d∂j(b''j(x).) + c(x) on Rd from the point of view of its numerical approximations in terms of matrices An having compartmental structure, that is (An)ii 0, (An)ij ≤ 0, i ≠ j, Σi(An)ij ≥ 0. We solve numerically the corresponding Dirichlet problem on a bounded domain D ⊂ Rd (d = 2, 3), for which the right-hand side is a probability measure with support in D. Numerical solutions on grids are nonpositive, and can be naturally embedded into linear spaces of 'hat' functions approximating the original solution in W11 (D). Numerical solutions converge in L1(D).The construction of our approximations is valid for general dimensions, but we give the convergence proof only for d = 2, 3. We end by a nontrivial example that illustrates the obtained results.