5
$\begingroup$

I have question connected with finite volume method. Consider equation

$$\frac{\partial u}{\partial t}=\operatorname{div}A\nabla u +f . \quad x\in \overline{B}_{1} (0)\subset \mathbb{R}^3 -\text{unit ball with its boundary}$$ I solved IBVP for this equation with a help of finite volume method, for grid I took uniform spherical coordinates grid,. (picture lower is related to topic) enter image description here

I have a problem with finding flows for the faces of cell: how to approximate the value $$\oint_{\partial V_{i,j,k}}(A\nabla u, n)dS $$ (where $A$ is a variable matrix)? When $A$ is an identity matrix everything is clear. Could you help me with this issue or to advise literature where this issue is represented?

| cite | improve this question | |
$\endgroup$
  • $\begingroup$ Product of A and grad u is a vector quantity - I would try to find value of each component of that vector at cell face. I see that as the main difficulty. The way you will do that will probably be similar to the procedure for discretization of diffusion in general transport eq. I'm not sure you are probably aware of that and look for more specific answer. $\endgroup$ – Johntra Volta Feb 4 '14 at 16:14
  • $\begingroup$ Actually, I thought almost in this way: to find the value of gradient, normal vector and matrix in the center of the cells' face and just write scalar product, and then just add these values for all faces. But also I thought that there are some special methods for doing it for my case. $\endgroup$ – cool Feb 4 '14 at 17:11
  • $\begingroup$ Or to find the product of the matrix and grad u, and then interpolate. Finding cell face values of any scalar (including individual components of your vector) using least squares gradients at cell centers may be a good idea... $\endgroup$ – Johntra Volta Feb 4 '14 at 17:25
  • $\begingroup$ Isn't the governing equation in Cartesian coordinate system? First you will need to write the equation in spherical coordinate system. $\endgroup$ – Discrete_Reynolds Jan 28 '17 at 9:27

Your Answer

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

Browse other questions tagged or ask your own question.