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

Non-physcial solution in the case of triangular mesh from interFoam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 11, 2011, 07:43
Default
  #21
ata
Senior Member
 
ata's Avatar
 
ata kamyabi
Join Date: Aug 2009
Location: Kerman
Posts: 323
Rep Power: 17
ata is on a distinguished road
Hi Niels Gjoel Jacobsen
I know it was in OF1.6 and in the OF1.7 p_rgh has been used but as I know they used p_rgh to make BC more comfortable grad(p_rgh)=0 but as I said I think we can solve this problem by using p in the momentum equation and setting grad(p)=grad(rgh) on the walls.
I do not know do you understand me?

Good luck
Best regards

Ata
ata is offline   Reply With Quote

Old   October 11, 2011, 23:21
Default
  #22
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,273
Rep Power: 34
arjun will become famous soon enougharjun will become famous soon enough
@Ata, Does openFOAM able to run standard test problems with triangular or tet meshes. For example there must be simple tutorial that openFOAM can run. What happens if you keep everything the same and just change the mesh?? Are you able to run the solver.

A simple test could be a water column falling etc etc.
arjun is offline   Reply With Quote

Old   October 12, 2011, 03:54
Default
  #23
ata
Senior Member
 
ata's Avatar
 
ata kamyabi
Join Date: Aug 2009
Location: Kerman
Posts: 323
Rep Power: 17
ata is on a distinguished road
Hi
Yes I examine some different grids and in all triangular grids I had the same problem.
Good luck

Ata
ata is offline   Reply With Quote

Old   October 12, 2011, 06:25
Default
  #24
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,273
Rep Power: 34
arjun will become famous soon enougharjun will become famous soon enough
Quote:
Originally Posted by ata View Post
Hi
Yes I examine some different grids and in all triangular grids I had the same problem.
Thats too bad.


Quote:
Originally Posted by ata View Post
Good luck
Thank you. I need lots of luck.

In my case though, my solver works well irrespective of the type of grid (if no negative volume cells exist then usually no problems). Grid type is not much a problem for me.
arjun is offline   Reply With Quote

Old   October 12, 2011, 06:43
Default
  #25
ata
Senior Member
 
ata's Avatar
 
ata kamyabi
Join Date: Aug 2009
Location: Kerman
Posts: 323
Rep Power: 17
ata is on a distinguished road
Hi
In your cases how much is the density ratio?
Good luck

Ata
ata is offline   Reply With Quote

Old   October 12, 2011, 07:16
Default
  #26
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,273
Rep Power: 34
arjun will become famous soon enougharjun will become famous soon enough
Quote:
Originally Posted by ata View Post
Hi
In your cases how much is the density ratio?
Good luck

Ata

1000 : 1

Water air.


Edited to add: So far as per my experience with multiphase simulations with various softwares like Fluent, CFX, StarCCM+ etc, I think the viscosity ratio is much much more problematic than density ratio. 1000 : 1 density ratio is usually not much a problem.

Last edited by arjun; October 12, 2011 at 07:45.
arjun is offline   Reply With Quote

Old   March 6, 2013, 17:25
Default
  #27
New Member
 
Join Date: Jan 2010
Posts: 23
Rep Power: 16
jdiorio is on a distinguished road
Similar issue here. Test problem of a sphere entering a free-surface. Using OF-2.1.0 with sliding AMI interfaces. InterDyMFoam solver although issue occurs with others as well. Results for cavity shape/pinch-off etc. are fine with cartesian grids. When there are just a few non-cartesian cells however, the cavity closes unphysically (see image).

Have a copy of Wemmenhove's paper. Their analysis makes sense to me. I am going to start looking into the evaluation of the snGrad terms in the momentum equation (pd and rho), unless others have some ideas. Has there been anymore discussion on this topic? This is the only thread I could find dealing with the issue.

Also, changing the snGrad scheme from corrected to limited did not fix the issue (results shown here are with the limited scheme).
Attached Images
File Type: jpg snapped_v_castellated_t0p030.jpg (95.6 KB, 60 views)
jdiorio is offline   Reply With Quote

Old   July 8, 2013, 06:34
Default opposite experience with mesh geometry
  #28
Member
 
Join Date: Mar 2013
Posts: 98
Rep Power: 13
giack is on a distinguished road
Hi to all,
this is a very interesting thread. I simulated a multiphase case in which a horizontal closed pipe filled with water is opened at one side at time zero, allowing the gas to enter and liquid to exit the domain.
I used a hexahedral mesh with all possible scheme for gradient and divergence scheme (also cellMDLimited version) and the result are very bad.
So i tried to change mesh geometry and I used a tetrahedral mesh. In this case the results are very good.

So this experience is opposite to yours.
The reason is probabily that using a tetahedral mesh there isn't a preferential direction for gradient calculations. What is your opinion?

Thank to all
giack is offline   Reply With Quote

Old   September 4, 2016, 12:59
Default
  #29
New Member
 
Saumitra Joshi
Join Date: Dec 2012
Posts: 14
Rep Power: 13
Zaphod'sSecondHead is on a distinguished road
Quote:
Originally Posted by arjun View Post
In my case though, my solver works well irrespective of the type of grid (if no negative volume cells exist then usually no problems). Grid type is not much a problem for me.
Hey Arjun,

I'm working on building an Adaptive Mesh Refinement algorithm based on the current version, but that also deals with arbitrary meshes. I face a huge problem of spurious currents in my solution - the grid I'm using is mainly tetrahedral.

Could you share what corrections you used to make your code independent of cell shape and orientation?
Zaphod'sSecondHead 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
CFL Condition Matt Umbel Main CFD Forum 19 June 30, 2020 09:20
[snappyHexMesh] snappyHexMesh won't work - zeros everywhere! sc298 OpenFOAM Meshing & Mesh Conversion 2 March 27, 2011 22:11
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 03:58
Discussion about Mesh independant solution Seb Main CFD Forum 13 May 22, 2001 14:37
unstructured vs. structured grids Frank Muldoon Main CFD Forum 1 January 5, 1999 11:09


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