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

rhoPimpleFoam Boundary Condition Problem

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 10, 2012, 23:21
Default rhoPimpleFoam Boundary Condition Problem
  #1
Senior Member
 
Daniel
Join Date: Jul 2009
Location: Montreal, Canada
Posts: 151
Rep Power: 8
dancfd is on a distinguished road
Hello all,

I have run my 2D airfoil case successfully in pimpleFoam, but now that I try to migrate my case to rhoPimpleFoam, the boundary conditions misbehave and the simulation stops with a floating point exception. I am running a NACA 0012 airfoil at M 0.3 with Re 6e6. Using MaxCo of 0.5. I have attached the boundary conditions and thermophysical properties, as well as the checkMesh results. I have also attached three images of pressure distributions across the domain: 1) from pimpleFoam that looks the way it should, 2) from rhoPimpleFoam using the same BCs, and 3) from rhoPimpleFoam defining P at the inlet instead of the outlet. It is clear that wherever I define P leads to anomalous results. Note that in this mesh, the outlet is the rightmost limit of the domain, and the inlet is the top, bottom and curved surfaces, as in the nacaAirfoil tutorial. InletOutlet is used as a BC for U since I will be running the simulation at various angles of attack by changing the BCs, not the mesh. I want to use a compressible solver since I will be going as high as M 0.8, and I expect shocks on the airfoil.

Here is the error:

Code:
[2] #0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[2] #1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[2] #2   in "/lib/x86_64-linux-gnu/libc.so.6"
[2] #3  Foam::hPsiThermo<Foam::pureMixture<Foam::sutherlandTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::calculate() in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
[2] #4  Foam::hPsiThermo<Foam::pureMixture<Foam::sutherlandTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::correct() in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
Can someone please suggest what the issue might be, or if the boundary conditions can be improved?

Many thanks in advance,

Dan
Attached Images
File Type: jpg 1-pimplefoam_P_def_at_outlet.jpg (31.2 KB, 51 views)
File Type: jpg 2-P_def_at_outlet_like_pimplefoam_case.jpg (65.9 KB, 69 views)
File Type: jpg 3-P_def_at_inlet.jpg (25.3 KB, 57 views)
Attached Files
File Type: zip case.zip (5.1 KB, 26 views)
dancfd is offline   Reply With Quote

Old   May 11, 2012, 02:27
Default
  #2
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
I'd consider rhoCentralFoam.

Best,
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   May 11, 2012, 06:31
Default
  #3
Senior Member
 
Daniel
Join Date: Jul 2009
Location: Montreal, Canada
Posts: 151
Rep Power: 8
dancfd is on a distinguished road
Hello Alberto,

Thanks for your suggestion. I ran a case with P defined at the inlet, and although I did not receive a floating point error in the time that I let the case run, the pressure distribution had the same issue as picture #2 above.

I am currently running a case with P defined at the outlet, and I will post the results when they are ready.

Regards,

Dan
dancfd is offline   Reply With Quote

Old   May 11, 2012, 18:20
Default
  #4
Senior Member
 
Daniel
Join Date: Jul 2009
Location: Montreal, Canada
Posts: 151
Rep Power: 8
dancfd is on a distinguished road
It looks like the boundary condition issue has been resolved by defining P at the outlet instead of the inlet and by using rhoCentralFoam. The p distribution is now free of anomalies. However, the Cp plot across the airfoil now has a wavy characteristic and the data points are spread out much more than they should be. Angle of attack is <1 deg for both plots. Any ideas what could cause this disruption to the Cp plot?

Thanks,

Dan
dancfd is offline   Reply With Quote

Old   May 14, 2012, 21:56
Default
  #5
Senior Member
 
Daniel
Join Date: Jul 2009
Location: Montreal, Canada
Posts: 151
Rep Power: 8
dancfd is on a distinguished road
I gave up on rhoCentralFoam. SonicFoam looks promising, however when I change from Euler to backward differencing for the temporal terms, I get strange results (see pics). I think it is necessary to use a second-order scheme such as backward if I want to publish, therefore this is a concern for me. Has anyone else encountered this problem?

Note that CrankNicholson 1 does not converge (floating point error) and CrankNicholson 0.5 produced a similar effect to that shown below.

Thank you,

Dan
Attached Images
File Type: jpg p_dist.jpg (43.6 KB, 50 views)
File Type: jpg p_dist_euler.jpg (80.5 KB, 35 views)

Last edited by dancfd; May 15, 2012 at 20:04.
dancfd is offline   Reply With Quote

Old   May 15, 2012, 00:24
Default
  #6
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
How do you set your case up in rhoCentralFoam? Take a look at the tutorials. Also, do you have a small case (it should run on 1 CPU in a short time) that reproduces your problem?

Best,
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   May 16, 2012, 23:09
Default
  #7
Senior Member
 
Daniel
Join Date: Jul 2009
Location: Montreal, Canada
Posts: 151
Rep Power: 8
dancfd is on a distinguished road
Hello Alberto,

I have tried to apply the case setups from the tutorial on my case, which is how I got as far as I did. I have attached my case, though I had to make a couple of changes for it to be within the size limits of this board. I am afraid it takes quite some time to run, unfortunately - blame the low courant number limit for that one!

I also tried changing the discretization of the time terms from backward to Euler, since that worked for sonicFoam. Unfortunately, it made no visible difference to the results with rhoCentralFoam - they still are not representative of experimental data.

I would appreciate any assistance you could offer.

Thanks,
Dan
Attached Files
File Type: zip basic_tadmor2_OFpost_nmsh.zip (38.2 KB, 21 views)
dancfd is offline   Reply With Quote

Old   May 16, 2012, 23:32
Default
  #8
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
I am looking into it. At what time do you start seeing the problem? Also, what version of OF are you using?
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods

Last edited by alberto; May 17, 2012 at 00:17. Reason: Added question on version
alberto is offline   Reply With Quote

Old   May 17, 2012, 12:47
Default
  #9
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Sorry, in my previous post I think I mixed the issues you have in rhoCentralFoam and those in sonicFoam.
I confirm the scatter in values using rhoCentralFoam. Have you checked that the solution stops changing?
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   May 17, 2012, 17:13
Default
  #10
Senior Member
 
Mahdi Hosseinali
Join Date: Apr 2009
Posts: 124
Rep Power: 8
anishtain4 is on a distinguished road
I tried to reproduce the case but it takes too much time, however there are two points that came to my mind:
1. change the velocity outlet from inletoutlet to zerogradient, after all your grid is far enough to be free of any recirculation, isn't it?
2. change the interpolation method, that maybe is the case for a non-smooth pressure on airfoil
anishtain4 is offline   Reply With Quote

Old   May 17, 2012, 20:18
Default
  #11
Senior Member
 
Daniel
Join Date: Jul 2009
Location: Montreal, Canada
Posts: 151
Rep Power: 8
dancfd is on a distinguished road
I am not sure if it appears any sooner that the time in the file. I had set the transient sim to output data at that interval. I agree that it takes a long time - is there any way to accelerate this sim? The Co must stay low or it will diverge. I am using OF 2.0.1. The solution should stop, since the convergence criteria are very demanding and I normally run the simulation with a time-varying flow speed. Thank you Mahdi, Alberto for your assistance - I will try the zeroGradient condition and another interpolation method and write back.

Regards,
Dan
dancfd is offline   Reply With Quote

Old   May 18, 2012, 04:45
Default
  #12
Senior Member
 
Mahdi Hosseinali
Join Date: Apr 2009
Posts: 124
Rep Power: 8
anishtain4 is on a distinguished road
Well you have a too dense grid, I'm not sure if it is necessary or not? This confines your Co to very low values, and you are solving a viscose flow with turbulence model.

Maybe it is better to first run this case as inviscid to come along with a good farfield boundary condition, and then take care of wall anomalities
anishtain4 is offline   Reply With Quote

Old   May 18, 2012, 21:47
Default
  #13
Senior Member
 
Daniel
Join Date: Jul 2009
Location: Montreal, Canada
Posts: 151
Rep Power: 8
dancfd is on a distinguished road
The zeroGradient condition did not have a noticeable effect on the results - the Cp plot is still oscillatory. I changed all of the interpolation schemes to linear and the simulation did not converge. Regarding the mesh, I believe the y+ is reasonable for the SST (~15) and the aspect ratios are acceptable to checkMesh. In which direction do you think I could coarsen the mesh?

The mesh notwithstanding (since that would allow me to run a faster sim, but I do not believe it would remove the oscillations), any other ideas on how to make rhoCentralFoam work with this case?

Thanks,
Dan
dancfd 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
Impinging Jet Boundary Condition Problem mattyg101 FLUENT 6 January 29, 2013 11:35
CFX two-phase cyclic boundary condition problem ukbid CFX 1 May 2, 2012 04:09
Transient Simulation: Boundary Condition Problem Shafiul CFX 7 January 11, 2011 17:40
problem about periodic boundary condition in Fluent winnawinna FLUENT 0 December 29, 2010 00:32
RPM in Wind Turbine Pankaj CFX 9 November 23, 2009 05:05


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