|July 8, 2010, 11:36||
how to create a volScalarField of mesh.V() ???
Join Date: Feb 2010
Posts: 9Rep Power: 7
I use a volScalarField with the volume of cells.
It works with the cell center but not with cell volume.
- Why the following line works
volVectorField centres = Sj.mesh().C();
- Why the following line dosn't work
volScalarField volume= Sj.mesh().V();
Sj is a volVectorField defined as follow:
|July 2, 2011, 07:48||
Join Date: Jun 2010
Posts: 107Rep Power: 7
Looking at the doxygen documentation on the openfoam.com website, you can see that:
Did you try using:
|July 6, 2011, 10:17||
Join Date: Mar 2009
Location: Besigheim, Germany
Posts: 134Rep Power: 8
This is due that a volScalarField does store values on the boundary, what does not make a lot of sense for cell volumes.
So just the internalField of a volScalarField does have cell volumes
So I personally would not try to cast this into a volScalarField! What will you do on the boundaries? If you initialize the volScalarField with zero than there will be zero at the boundaries too. What happens if you divide at a point by these values?
I would just work on the internalField
volScalarField myWhatEverField =mag(U);
scalarField volumes = mesh.V();
volScalarField result(IOobject(...),mesh, 0);
result.internalField() = myWahateverField.internalField/volumes;
Do whatever you need to do on the boundaries
|Thread||Thread Starter||Forum||Replies||Last Post|
|How to create initiate a volScalarField p without reading from disk NO_READ does not seem to work||dbxmcf||OpenFOAM Running, Solving & CFD||12||August 22, 2013 07:32|
|Actuator disk model||audrich||FLUENT||0||September 21, 2009 07:06|
|Where's the singularity/mesh flaw?||audrich||FLUENT||3||August 4, 2009 01:07|
|fluent add additional zones for the mesh file||SSL||FLUENT||2||January 26, 2008 12:55|
|Meshing a Sphere||Ajay||FLUENT||9||March 29, 2004 09:14|