Dear deal.ii community, I am new to deal.ii and trying to solve the
following set of equations with homogeneous Neumann BC:
$$u_t = \Delta u - \nabla \cdot (\lambda(u) \nabla v)+ f(u,v)$$
$$v_t = \Delta v + g(u,v)$$
Due to radial symmetry it reduces to:
$$u_t = u_{rr}+ \frac{1}{r} u_r - \frac{1}{r} \frac{\partial}{\partial r}
(r \lambda(u) \frac{\partial v}{\partial r}) +f(u,v)$$
$$v_t = v_{rr} + \frac{1}{r} v_r + g(u,v) $$
But at $r=0$, it becomes :
$$u_t = 2u_{rr}-2uv_{rr}-u_r h_r +f(u,v)$$
$$v_t = 2v_{rr}+g(u,v)$$
So, my doubts are:
1. Since, its a time dependent, vector valued \& non-linear set of PDEs,
which tutorials should I follow to solve these? I found 21, 26, 31 \& 31
related but some of them are too complicated, is there any relatively
simpler way !! \\
2. How to handle the condition at $r=0$ with the homogeneous Neumann BC?
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
For more options, visit https://groups.google.com/d/optout.
