1-D spherical finite volume discretisation

lost.identity · October 6, 2010, 09:54

I'm trying to discretise the conservative equations using FVM for a 1-D spherical geometry (i.e. spherical symmetry).

After applying the Gaus's theorem I get the following for spatial discretisation of the diffusive terms

$\boldsymbol{\Gamma_f}\cdot\boldsymbol{A}_f$

I know that for 1-D Cartesian the face area vector is simply $\boldsymbol{A}=\boldsymbol{i}$ where i is the unit vector in the x-direction.

However, what will be the face area vector for a spherical geometry? Surely it can't just be the unit vector.

I know that the area element in spherical coordinates is given by

$d\boldsymbol{A} = r^2\sin\theta{d}\theta{d}\phi\boldsymbol{e}_r$ symbol. So in 1-D can I assume it to be just r^2?

Thanks.

david_e · October 7, 2010, 15:36

You can take it to be

, because to get the conservative formulation one usually integrates over $\theta$ and $\phi$ .

October 6, 2010, 09:54	1-D spherical finite volume discretisation	#1
lost.identity Senior Member Join Date: Apr 2009 Posts: 118 Rep Power: 10	I'm trying to discretise the conservative equations using FVM for a 1-D spherical geometry (i.e. spherical symmetry). After applying the Gaus's theorem I get the following for spatial discretisation of the diffusive terms $\boldsymbol{\Gamma_f}\cdot\boldsymbol{A}_f$ I know that for 1-D Cartesian the face area vector is simply $\boldsymbol{A}=\boldsymbol{i}$ where i is the unit vector in the x-direction. However, what will be the face area vector for a spherical geometry? Surely it can't just be the unit vector. I know that the area element in spherical coordinates is given by $d\boldsymbol{A} = r^2\sin\theta{d}\theta{d}\phi\boldsymbol{e}_r$ symbol. So in 1-D can I assume it to be just r^2? Thanks.

Thread Tools
Show Printable Version Email this Page
Display Modes
Linear Mode Switch to Hybrid Mode Switch to Threaded Mode

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
BlockMesh FOAM warning	gaottino	OpenFOAM Native Meshers: blockMesh	7	July 19, 2010 14:11
Control Volume , Finite Volume, Finite Control Volume, Finite Element Method	technocrat.prakash	Main CFD Forum	1	April 24, 2010 19:24
On the damBreak4phaseFine cases	paean	OpenFOAM Running, Solving & CFD	0	November 14, 2008 22:14
fluent add additional zones for the mesh file	SSL	FLUENT	2	January 26, 2008 12:55
Axisymmetrical mesh	Rasmus Gjesing (Gjesing)	OpenFOAM Native Meshers: blockMesh	10	April 2, 2007 14:00

October 7, 2010, 15:36		#2
david_e New Member Join Date: May 2010 Posts: 8 Rep Power: 9	You can take it to be , because to get the conservative formulation one usually integrates over $\theta$ and $\phi$ .