CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Discontinuity at the blunt trailing edge of an airfoil

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

LinkBack Thread Tools Search this Thread Display Modes
Old   October 28, 2020, 04:42
Default Discontinuity at the blunt trailing edge of an airfoil
New Member
Join Date: Sep 2020
Posts: 1
Rep Power: 0
klkl is on a distinguished road
Hey Foamers,

I am a beginner to openFoam and CFD. My research project is to study the airframe noise produced near the blunt trailing edge of a NACA 0012 airfoil.

I am trying to run a transient compressible simulation on openFoam, using rhoPimpleFoam as the solver and kOmegaSST as the turbulence model.

There is some discontinuity of pressure at the trailing edge, (likely induced by numerical instability). (See pics). In the early stage of computation, there is still some vortex shedding pattern near the trailing edge. In the later stage of computation, the vortex shedding disappears and there are only lines of discontinuity of pressure.

I am not sure what have I done wrong that causes such a discontinuity. I defined the initial alpha, k, nut, omega, p, T, U based on the rhoPimpleFoam/RAS/aerofoilNACA0012 tutorial. The value of k/nut/omega should be alright since I already have a successful transient incompressible simulation using pisoFoam.

The structured mesh should not be an issue, since it is built based on a reference article. I have turned on the adjustTimeStep and fixed the maxco to 0.5, so it should not be the problem of CFL. I guess it may be related to the fvSchemes or fvSolutions, but I am not sure.

Here are the numerical schemes I used.
default backward;

default Gauss linear;

default none;

div(phi,U) Gauss linearUpwind grad(U);

turbulence Gauss upwind;
energy Gauss upwind;

div(phi,k) $turbulence;
div(phi,omega) $turbulence;

div(phi,e) $energy;
div(phi,K) $energy;
div(phi,Ekp) $energy;

div(phiv,p) Gauss upwind;
div((phi|interpolate(rho)),p) Gauss upwind;

div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;

default Gauss linear corrected;

default linear;

default corrected;

method meshWave;

Here are the solvers that I used.
solver GAMG;
smoother GaussSeidel;
tolerance 1e-6;
relTol 0.01;

relTol 0;

solver smoothSolver;
smoother GaussSeidel;
nSweeps 2;
tolerance 1e-08;
relTol 0.1;

relTol 0;

nCorrectors 2;
nNonOrthogonalCorrectors 1;
nOuterCorrectors 1;

I am wondering if any of you encounter a similar problem before or have any idea what should I do to get rid of this discontinuity?

Many thanks!!!

Attached Images
File Type: jpg discontiuity-1.jpg (32.5 KB, 19 views)
File Type: jpg discontiuity-2.jpg (41.0 KB, 14 views)
File Type: jpg flow.jpg (27.6 KB, 10 views)
File Type: jpg mesh.jpg (100.9 KB, 10 views)
klkl is offline   Reply With Quote


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
SU2 AOA optimization SU2 9 March 7, 2022 16:17
how to set periodic boundary conditions Ganesh FLUENT 15 November 18, 2020 06:09
[DesignModeler] 3mm Trailing edge radius on closed airfoil henrikdahlberg ANSYS Meshing & Geometry 1 June 10, 2015 05:04
Leading edge and trailing edge of an airfoil rotor blade Jun_Milan CFX 0 July 10, 2014 23:07
Problem with restart solution in robyTKD SU2 Shape Design 21 May 29, 2013 09:26

All times are GMT -4. The time now is 05:29.