|
[Sponsors] |
June 18, 2016, 14:13 |
Internal mesh
|
#1 |
Member
anonymous
Join Date: Mar 2016
Location: Canada
Posts: 93
Rep Power: 10 |
Hi,
I would like to perform some operations on the internal mesh only. I used the following code, but it also changes my boundary condition. How can I get access to internal mesh only. Need help please. Code:
forAll(FC, cellI) { FC[cellI] = k1.value() * C[cellI] / (k2.value() + C[cellI]); } |
|
June 20, 2016, 18:55 |
|
#2 |
Member
anonymous
Join Date: Mar 2016
Location: Canada
Posts: 93
Rep Power: 10 |
Anyone has any comment or idea.
|
|
June 21, 2016, 02:20 |
|
#3 |
Member
Rodrigo
Join Date: Mar 2010
Posts: 98
Rep Power: 16 |
What about
FC.internalField()=k1.value()*C.internalField()/(k2.value()+C.internalFied()) ? If this compiles it shall do what you want |
|
June 21, 2016, 03:34 |
|
#4 |
Member
anonymous
Join Date: Mar 2016
Location: Canada
Posts: 93
Rep Power: 10 |
Thanks, it compiled. But it did not solve my problem.
I am solving for a concentration in a box with dimension 1 x 1 x 0.5 mm^3. I set all sides and bottom with zeroGradient B.C. and top with uniform value of 1. Initially I set the initial internal field as ZERO. When I run my case it also changes the values in all sides, it supposed to be zero. I thought that may be it is changing due to my operation on every cell as mentioned above, therefore I wanted to try updating internal mesh only. But it did not work out. Need help. Any suggestion. Thanks. |
|
June 21, 2016, 03:43 |
|
#5 |
Senior Member
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,377
Rep Power: 29 |
Reread your post carefully. Especially the parts "it's supposed to be zero" and "I set it to zeroGRADIENT".
__________________
*On twitter @akidTwit *Spend as much time formulating your questions as you expect people to spend on their answer. |
|
June 21, 2016, 03:51 |
|
#6 |
Member
anonymous
Join Date: Mar 2016
Location: Canada
Posts: 93
Rep Power: 10 |
Internal field is zero and all side walls has zeroGradient.
What I don't understand is that, why the values at walls are changing? |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[mesh manipulation] how should we remove the Concave cells from the mesh? | s.m | OpenFOAM Meshing & Mesh Conversion | 20 | March 5, 2020 17:35 |
[snappyHexMesh] SnappyHexMesh no layers and no decent mesh for complex geometry | pizzaspinate | OpenFOAM Meshing & Mesh Conversion | 1 | February 25, 2015 08:05 |
fluent add additional zones for the mesh file | SSL | FLUENT | 2 | January 26, 2008 12:55 |
Convergence moving mesh | lr103476 | OpenFOAM Running, Solving & CFD | 30 | November 19, 2007 15:09 |
Icemcfd 11: Loss of mesh from surface mesh option? | Joe | CFX | 2 | March 26, 2007 19:10 |