CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

Internal mesh

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 18, 2016, 14:13
Default Internal mesh
  #1
Member
 
anonymous
Join Date: Mar 2016
Location: Canada
Posts: 93
Rep Power: 10
cute is on a distinguished road
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]);
  }
Thanks.
cute is offline   Reply With Quote

Old   June 20, 2016, 18:55
Default
  #2
Member
 
anonymous
Join Date: Mar 2016
Location: Canada
Posts: 93
Rep Power: 10
cute is on a distinguished road
Anyone has any comment or idea.
cute is offline   Reply With Quote

Old   June 21, 2016, 02:20
Default
  #3
Member
 
Rodrigo
Join Date: Mar 2010
Posts: 98
Rep Power: 16
guin is on a distinguished road
What about
FC.internalField()=k1.value()*C.internalField()/(k2.value()+C.internalFied()) ? If this compiles it shall do what you want
guin is offline   Reply With Quote

Old   June 21, 2016, 03:34
Default
  #4
Member
 
anonymous
Join Date: Mar 2016
Location: Canada
Posts: 93
Rep Power: 10
cute is on a distinguished road
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.
cute is offline   Reply With Quote

Old   June 21, 2016, 03:43
Default
  #5
Senior Member
 
akidess's Avatar
 
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,377
Rep Power: 29
akidess will become famous soon enough
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.
akidess is offline   Reply With Quote

Old   June 21, 2016, 03:51
Default
  #6
Member
 
anonymous
Join Date: Mar 2016
Location: Canada
Posts: 93
Rep Power: 10
cute is on a distinguished road
Internal field is zero and all side walls has zeroGradient.
What I don't understand is that, why the values at walls are changing?
cute is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


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


All times are GMT -4. The time now is 11:21.