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

Some problems in naca 0012 V&V case of NASA TMR and DPW using OpenFOAM

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 10, 2017, 10:08
Question Some problems in naca 0012 V&V case of NASA TMR and DPW using OpenFOAM
  #1
Member
 
Di Cheng
Join Date: May 2010
Location: Beijing, China
Posts: 47
Rep Power: 16
chengdi is on a distinguished road
Hi, everyone

I am working on naca0012 validation case in NASA TMR using the grid (the most coarse one for now) and model provided by TMR. It is required to join DPW 6.

2D case with Mach 0.15 and Re=6e6. Seems easy at first glance.
However, there are a lot of obstacles and traps!!! Including:
1. max aspect ratio of grid (most coarse one, familty-II) is 2.4E7! It did not past the checkMesh. However it can be run.
2. SA model in OpenFOAM is not a standard SA.
3. The effective Prandtl number of original sutherlandTransport model of OpenFOAM is not 0.72 but 0.69. (sutherlandTransport is using Eugene's equation rather than Prandtl number approach)
4. The turbulent Prandtl number of EddyViscosity is 1.0 by default.
5. No Riemann BC in OpenFOAM.
6. The 500*cord far field boundary is so large that I need to run too many timesteps to reach steady state with an unsteady solvers. However, there is only one compressible steady state solver available in OpenFOAM 201612+ (the rhoSimpleFoam) which seems not very good for highly compressible flow.

Now I am running into turbulence. I just found that there is no example case of compressible SA model in OpenFOAM 1612+ and foam extend 4.0.

And I found that the initial condition and solver configuration existing in incompressible SA model case does not work in compressible SA case becasue of lacking `0/alphat` field.

And the alphat field is not thermal diffsivity. It seems like `density*thermal diffusivity` because its dimension must be [1 -1 -1 0 0 0 0]. (confirmed by the implementation of SA model code)

I am totally frustrated by this validation case. Is there anyone have any experience on doing DPW-like verification with naca 0012 using OpenFOAM?
chengdi is offline   Reply With Quote

Old   May 12, 2017, 09:13
Default
  #2
New Member
 
Ajeje Brazov
Join Date: Apr 2017
Posts: 14
Rep Power: 9
Bernoulli666 is on a distinguished road
I run the simulation with the parameters in the folder and Mesh size 897 and seems that results are very similar to the nlrc ones.
Attached Files
File Type: zip SANACA0012.zip (5.9 KB, 36 views)
Bernoulli666 is offline   Reply With Quote

Old   May 13, 2017, 21:12
Default
  #3
Member
 
Di Cheng
Join Date: May 2010
Location: Beijing, China
Posts: 47
Rep Power: 16
chengdi is on a distinguished road
NASA TMR requires that the case should be calculated using compressible code rather than simpleFoam (an incompressible code). I tried to used rhoCentralFoam in unsteady manner, it works. However the time is too long. I also tried to use rhoSimpleFoam, it is not working.
chengdi is offline   Reply With Quote

Old   May 14, 2017, 07:31
Default
  #4
New Member
 
Ajeje Brazov
Join Date: Apr 2017
Posts: 14
Rep Power: 9
Bernoulli666 is on a distinguished road
You can also use compressibile solvers for uncompressible flows. If you look in the tutorial folder "incompressibile" there are cases solved with simple algorithm.


Sent from my iPad using CFD Online Forum mobile app
Bernoulli666 is offline   Reply With Quote

Old   May 15, 2017, 00:52
Default
  #5
Member
 
Di Cheng
Join Date: May 2010
Location: Beijing, China
Posts: 47
Rep Power: 16
chengdi is on a distinguished road
Yes, I can solve it even with potentialFoam. However, according to NASA's comment, the solver is required to be an incompressible solver.
chengdi is offline   Reply With Quote

Old   May 15, 2017, 02:29
Default
  #6
Senior Member
 
Mikko
Join Date: Jul 2014
Location: The Hague, The Netherlands
Posts: 243
Rep Power: 13
Flowkersma is on a distinguished road
Hi,

My first question is, why are you using compressible solver? I have run the same case with simpleFoam and kOmegaSST turbulence model and I got very similar results to NASA's results.

1. High aspect ratios does not seem to be a problem for OpenFOAM even though checkMesh gives an error.
2. This depends on which OF version you are using? Older releases are using the not recommended one with fv3 term. I think that OF foundation is using now the "good" implementation without fv3 term. However, I am not sure about OF+. Unfortunately, foam-extend is probably still using the fv3 term implementation.
3-4. If the compressibility effects are almost negligible, do we care?
6. This flow is not highly compressible so I would use simpleFoam or if you want to stick to compressible solver then rhoSimpleFoam.

Regards,
Mikko
Flowkersma is offline   Reply With Quote

Old   May 16, 2017, 05:18
Default
  #7
Member
 
Di Cheng
Join Date: May 2010
Location: Beijing, China
Posts: 47
Rep Power: 16
chengdi is on a distinguished road
Hi, Mikko

For your questions:

1. Not just high aspect ratio problem. The highest AR element is located at the the wake baffle and it is very orthogonal rectangular element. The non-orthogonality might be the worst factor. The max angle is > 80 degree.

2. I am using compressible solver because the verification case requires that. And I just want to make sure the OpenFOAM 201612+ can be used to produce a result almost exactly as CFL3D/FUN3D and TAU of DLR. So I can assure that the program is correctly coded and can be used to deal with external aerodynamic problems such as DPW6.

3. For application, it is negligible and I will definitely use incompressible solver, however, this is a verification case. I must follow the guide as closely as I can.

4. rhoSimpleFoam is very difficult to run without divergence or thermophysical problem.
chengdi is offline   Reply With Quote

Old   October 5, 2019, 14:20
Default NACA0012 OF validation for Re = 6Million
  #8
New Member
 
Join Date: Dec 2016
Posts: 5
Rep Power: 9
Iain_symthe is on a distinguished road
Quote:
Originally Posted by Bernoulli666 View Post
I run the simulation with the parameters in the folder and Mesh size 897 and seems that results are very similar to the nlrc ones.

Thanks for attaching the zip file with your values. Can you also share your validation results (OF simulations with data) in text files. How many steps did you need to run for convergence? Looking at your 0/U and 0/nut files it looks like you are running a Re = 3M case (87/3.e-6) for Lref = 1



I am trying to run the Re = 10^6 case as outlined in

https://www.openfoam.com/documentati...irfoil-2d.html
using the n0012_897-257.p3dfmt file.


I had to do

plot3dToFoam -noBlank n0012_897-257.p3dfmt

followed by
autoPatch -overwrite 15


to generate a proper boundary file (with inlet, outlet, front, back and airfoil) in it. I have used values for U_alpha, nu and nuT as prescribed in the OF validation page and I have run the simulation with
simpleFoam for 10000 steps but it has still not converged.


It will be helpful if someone has run this case (Re = 1e6) to completion and compared their results to data as in the OF validation page.
I have attached my setup file (without polyMesh as it is too large).



Thanks
Iain
Attached Files
File Type: gz setup1.tar.gz (2.3 KB, 15 views)
Iain_symthe is offline   Reply With Quote

Reply

Tags
naca0012, tmr


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
How to run the naca 0012 V&V case of NASA TMR and DPW using OpenFOAM chengdi OpenFOAM Running, Solving & CFD 0 May 9, 2017 08:45


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