|
[Sponsors] |
September 3, 2013, 21:14 |
How about I write U by reconstruct(phi)?
|
#1 |
Senior Member
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 838
Rep Power: 17 |
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. |
|
September 5, 2013, 05:21 |
|
#2 |
Senior Member
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 838
Rep Power: 17 |
Hello guys??
|
|
September 11, 2013, 04:49 |
|
#3 |
Senior Member
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 838
Rep Power: 17 |
Any ideas?
|
|
September 17, 2013, 09:55 |
|
#4 |
Member
A. Bernath
Join Date: Jun 2011
Location: Karlsruhe, Germany
Posts: 39
Rep Power: 14 |
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? |
|
September 22, 2013, 21:30 |
|
#5 |
Senior Member
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 838
Rep Power: 17 |
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.
|
|
September 23, 2013, 06:28 |
|
#6 |
Member
A. Bernath
Join Date: Jun 2011
Location: Karlsruhe, Germany
Posts: 39
Rep Power: 14 |
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? |
|
September 23, 2013, 20:41 |
|
#7 | |
Senior Member
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 838
Rep Power: 17 |
Quote:
Actually, I have to make a new file in time file. I just copy U and rename it by myU. |
||
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
mesh file for flow over a circular cylinder | Ardalan | Main CFD Forum | 7 | December 15, 2020 13:06 |
[Helyx OS] Helyx-OS (GUI for SnappyHexMesh | elvis | OpenFOAM Community Contributions | 210 | January 30, 2017 18:57 |
How to write UDF programme inside FLUENT? | coolyihao | Fluent UDF and Scheme Programming | 9 | June 5, 2016 21:08 |
Phase locked average in run time | panara | OpenFOAM | 2 | February 20, 2008 14:37 |
Problem when I try to write the geometry file | massimo | Siemens | 0 | January 14, 2003 04:58 |