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

How about I write U by reconstruct(phi)?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 3, 2013, 21:14
Default How about I write U by reconstruct(phi)?
  #1
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 675
Rep Power: 8
sharonyue is on a distinguished road
Hi Foamers!

Wonderful day!

lol, AFAWK, in the algorithm, we can get the U after solving possion equation via this equation: U = HbyA - rAU*fvc::grad(p);

In OpenFOAM, We all know that phi is divergence free indtead of U, Well, how about writing U by myU = fvc::reconstruct(phi)? Shouldn't it be more precise?

I made a test on cavity, there is a little difference but not much.
Attached Images
File Type: jpg 1.jpg (56.1 KB, 30 views)
File Type: jpg 2.jpg (20.0 KB, 30 views)
sharonyue is offline   Reply With Quote

Old   September 5, 2013, 05:21
Default
  #2
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 675
Rep Power: 8
sharonyue is on a distinguished road
Hello guys??
sharonyue is offline   Reply With Quote

Old   September 11, 2013, 04:49
Default
  #3
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 675
Rep Power: 8
sharonyue is on a distinguished road
Any ideas?
sharonyue is offline   Reply With Quote

Old   September 17, 2013, 09:55
Default
  #4
Member
 
A. Bernath
Join Date: Jun 2011
Location: Karlsruhe, Germany
Posts: 39
Rep Power: 6
derkermit is on a distinguished road
Shouldn't reconstruct just (linear) interpolate the face centered values of (in this case) phi to the cell centers?
Do you get the same dimensions for myU and U?
derkermit is offline   Reply With Quote

Old   September 22, 2013, 21:30
Default
  #5
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 675
Rep Power: 8
sharonyue is on a distinguished road
Quote:
Originally Posted by derkermit View Post
Shouldn't reconstruct just (linear) interpolate the face centered values of (in this case) phi to the cell centers?
Do you get the same dimensions for myU and U?
Yep, its the same. And it surely interpolate the phi to the cell centers. Just as OpenFOAM make sure phi is divergence free. So I think myU would be more precise than U. but as of the algorithm, U can be get by another equation. So I'm wondering which method would be more precise.
sharonyue is offline   Reply With Quote

Old   September 23, 2013, 06:28
Default
  #6
Member
 
A. Bernath
Join Date: Jun 2011
Location: Karlsruhe, Germany
Posts: 39
Rep Power: 6
derkermit is on a distinguished road
How can they both have the same dimensions?
As an example: reconstruct(T) would give me the linear interpolated cell centered temperature values from face centered values. That will not affect the units or am I wrong with this?
If not, why would you get the same dimensions for myU and U?
derkermit is offline   Reply With Quote

Old   September 23, 2013, 20:41
Default
  #7
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 675
Rep Power: 8
sharonyue is on a distinguished road
Quote:
Originally Posted by derkermit View Post
How can they both have the same dimensions?
As an example: reconstruct(T) would give me the linear interpolated cell centered temperature values from face centered values. That will not affect the units or am I wrong with this?
If not, why would you get the same dimensions for myU and U?
I dont know if you can do reconstruct(t), but you can do reconstruct(phi) then you would get U at the cell center. more details see here and welcome comments: u_f , u, are these two the same?. so myU 's unit,which is get by reconstructing, would be the same with U.

Actually, I have to make a new file in time file. I just copy U and rename it by myU.
sharonyue is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Helyx-OS (GUI for SnappyHexMesh elvis OpenFOAM Native Meshers: snappyHexMesh and Others 176 July 10, 2015 16:50
mesh file for flow over a circular cylinder Ardalan Main CFD Forum 6 April 17, 2010 23:40
How to write UDF programme inside FLUENT? coolyihao Fluent UDF and Scheme Programming 4 April 1, 2009 02:07
Phase locked average in run time panara OpenFOAM 2 February 20, 2008 15:37
Problem when I try to write the geometry file massimo CD-adapco 0 January 14, 2003 05:58


All times are GMT -4. The time now is 19:42.