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

Error due to Unstructured Mesh with custom solver

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

Like Tree5Likes
  • 5 Post By fredo490

Reply
 
LinkBack Thread Tools Display Modes
Old   February 22, 2013, 08:11
Default Error due to Unstructured Mesh with custom solver
  #1
Senior Member
 
HECKMANN Frédéric
Join Date: Jul 2010
Posts: 236
Rep Power: 7
fredo490 is on a distinguished road
Edit, problem solved

Dear all,

I'm facing some issues with unstructured Mesh (quality is ok and Skewness is less than 0.45) with my custom solver. My "density variable" has some strange error lines that follows the flow streamlines as shown bellow. The picture shows the normalized density (1 is the freestream value) when the case reaches steadyness. I tried different scheme, the linearUpwindV scheme improve/smooth a little the results but it's not perfect yet.

My case computes the behavior of highly dispersed particles in a fluid (only the fluid impact the particles motion through the drag) using an Eulerian approach. The goal is to compute the particles collection over a cylinder (some particles impact the cylinder and some are moved away by the air flow via drag)


My solver:
An unsteady Eulerian solver to compute the behavior of highly dispersed tiny particles in a fluid (only the fluid impact the particles motion through the drag).

Nomenclature :
- \alpha volume fraction of the particles
- \vec{u} velocity of the particles
- f(Re) function of the relative Reynolds number
- \tau time respond of the particle
- \vec{V} velocity of the fluid


Continuity equation :
\frac{\partial \alpha}{\partial t} + {\nabla \alpha \vec{u}} = 0

Momentum equation
\frac{\partial \alpha \vec{u}}{\partial t} + {\nabla \alpha \vec{u} \vec{u}} = \frac{f(Re)}{\tau } \alpha (\vec{V} -\vec{u})

My scheme
Code:
ddtSchemes
{
    default         backward;
}

divSchemes
{
    default         no;
    div(phiR,U)     Gauss linearUpwind grad(U);
    div(phi,rho)    Gauss upwind grad(rho);
}

gradSchemes
{
    default         Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear;
}

interpolationSchemes
{
    default linear;
}

snGradSchemes
{
    default         fourth;
}

fluxRequired
{
    default         no;
    rho               ;
}
EDIT: I've found a solution
To solve my problem, I used a skewCorrected

My new fvSchemes is:
Code:
ddtSchemes
{
    default         backward;
}

divSchemes
{
    default         no;
    div(phiR,U)     Gauss skewCorrected linearUpwind grad(U);
    div(phi,rho)    Gauss skewCorrected upwind grad(rho);
}

gradSchemes
{
    default         Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear;
}

interpolationSchemes
{
    default linear;
}

snGradSchemes
{
    default         fourth;
}

fluxRequired
{
    default         no;
    rho               ;
}
chegdan, wyldckat, Pat84 and 2 others like this.

Last edited by fredo490; February 24, 2013 at 08:22. Reason: problem solved
fredo490 is offline   Reply With Quote

Old   February 23, 2013, 04:04
Default
  #2
Senior Member
 
HECKMANN Frédéric
Join Date: Jul 2010
Posts: 236
Rep Power: 7
fredo490 is on a distinguished road
I've checked my problem with a structured grid having the same y+ and the small disturbances are gone... The problem is linked to the unstructured but I cannot find any scheme able to smooth/correct those errors.
fredo490 is offline   Reply With Quote

Old   March 6, 2015, 03:46
Default
  #3
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 592
Rep Power: 6
sharonyue is on a distinguished road
Hello HECKMANN!

Could u pls upload a comparison of contour plot of the right result? Using skewCorrected scheme.

Thanks.
sharonyue 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
how to set periodic boundary conditions Ganesh FLUENT 13 January 22, 2014 05:11
No layers in a small gap bobburnquist OpenFOAM Native Meshers: snappyHexMesh and Others 2 November 25, 2012 09:54
[ICEM] Why do I have to convert my mesh to an unstructured mesh? iznish ANSYS Meshing & Geometry 1 November 2, 2012 09:29
[ICEM] Just started unstructured mesh ! diamondx ANSYS Meshing & Geometry 6 September 5, 2012 07:15
unstructured vs. structured grids Frank Muldoon Main CFD Forum 1 January 5, 1999 11:09


All times are GMT -4. The time now is 16:48.