Find $u\in H_D^1(\Omega)$ such that

  • $-\nabla\cdot(a\nabla u)=0$ in $\Omega$,
  • $\dfrac{\partial u}{\partial n}=g$ on $\Gamma_N$,
  • $u=0$ on $\Gamma_D$.

The function $a(x,y)$ is piecewise constant on each of the $9$ subdomains

$$\begin{array}{|c|c|c|} \hline \Omega_7&\Omega_8&\Omega_9\\ \hline \Omega_4&\Omega_5&\Omega_6\\ \hline \Omega_1&\Omega_2&\Omega_3\\ \hline \end{array}$$

(or see the picture here).

I tried multiplying by $v$ then integrating (the standard approach). I broke down the integral to have $9$ integrals, one for each subdomain, to use Green’s formula, but I’m stuck since now I have border integrals over all the subdomains...

Edit : $\Gamma_D$ is the top edge of the square, $\Gamma_N$ is the the union of the three other edges.

Edit 2 : I tried another approach, can someone JUST check if this is correct ? $$ \int_{\Omega} a\nabla u \cdot \nabla v = \int_{\partial \Omega }\frac{\partial u}{\partial n}}$$

  • $\begingroup$ Welcome to math.SE! This might help you learning the latex-based formatting program used in this site. $\endgroup$ – user228113 Jun 10 '15 at 1:35
  • $\begingroup$ Thanks for the edit it is much readable now. $\endgroup$ – mdrlol Jun 10 '15 at 1:46

I think that your first approach is correct.

You have to keep the weighted Neumann jumps $[a \frac{\partial u}{\partial n}]$ on the internal edges. Usually, these jumps are explicitly specified in the strong formulation.

I don't understand the formulation that you provide in your second edit. If $a$ is a function, why is it outside the integral? And where is the test function $v$ in the boundary part?

  • $\begingroup$ Thanks for answering, the problem is that the jumps you talk about are not specified in the strong formulation ! As for the last edit : I thought I could apply the formula grad(a(grad(u) (where a is a function) to obtain a.Laplacian(u)+grad(a).grad(u) and say that this is equal to a.Laplacian(u) then multiply by v to apply Green formula. Is this correct ? $\endgroup$ – mdrlol Jun 11 '15 at 16:59
  • $\begingroup$ The jumps are only specified on the 4 edges of the big square $\endgroup$ – mdrlol Jun 11 '15 at 17:00
  • $\begingroup$ I have found something very similar to my problem. How did they get rid of the internal edges ? fenicsproject.org/documentation/dolfin/1.2.0/python/demo/pde/… $\endgroup$ – mdrlol Jun 11 '15 at 17:57
  • $\begingroup$ I don't think that you are allowed to expand $\nabla (a\nabla u)$ $\endgroup$ – Albe Jun 12 '15 at 9:45
  • $\begingroup$ If $a$ is only piecewise constant, you have to rely on the notion of weak derivative and perform integration by parts. As for the link, you may assume that you have homogeneous weighted Neumann jumps. However, in general, these jump conditions are provided by the physical properties of the problem under consideration. $\endgroup$ – Albe Jun 12 '15 at 9:56

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.