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

IDDES-SA boundary conditions

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 16, 2011, 22:13
Default IDDES-SA boundary conditions
  #1
New Member
 
Johannes N Theron
Join Date: Feb 2010
Location: Hamburg
Posts: 21
Rep Power: 7
jantheron is on a distinguished road
I am trying to run a simple case of flow over a square cylinder at Re=20,000 using Spalart's DES model implemented in OF2.0.1. Can anyone shed some light on the appropriate boundary conditions to use for k, nuTilda, U, omega and p? I have slip walls with one simple inlet and one outlet (kind of like a wind tunnel).

My grid was generated with gmsh and extruded in the Z direction. It consists of purely hexahedra with maximum orthogonality of 45. It passed all grid checks.

When I run with pisoFoam (keeping Co below 0.4), all the residuals decrease nicely but I can't seem to get rid of some unevenness in the velocity contours - these seem to come from the inlet.

At first I used Gauss linear interpolation but switched to TVD (limitedLinear) which got rid of extreme oscillations in the velocity contours. Yet, these oscilaltions emanating from the inlet corners persist in spite of a very low time step (5e-6 sec). My y+ values are max 14 and average 6.

The domain extends about 25 body diameters around the body, 18 diameters upstream and 40 diameters downstream. I am interested in the effect of shed vorticity on the freestream velocity in the +Y direction (about 20 diameters above the body)

Jan Theron

This is a snippet from my log file:

**********
Time = 0.0035

Courant Number mean: 0.00305846 max: 0.333807
DILUPBiCG: Solving for Ux, Initial residual = 0.0105967, Final residual = 3.03098e-12, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.012991, Final residual = 6.97586e-12, No Iterations 2
GAMGPCG: Solving for p, Initial residual = 0.280629, Final residual = 0.00207036, No Iterations 1
time step continuity errors : sum local = 5.72685e-09, global = -2.32239e-10, cumulative = 4.26148e-10
GAMGPCG: Solving for p, Initial residual = 0.0890806, Final residual = 0.00107882, No Iterations 1
time step continuity errors : sum local = 2.82245e-09, global = -3.78255e-11, cumulative = 3.88323e-10
GAMGPCG: Solving for p, Initial residual = 0.0790894, Final residual = 3.05308e-09, No Iterations 8
time step continuity errors : sum local = 7.97931e-15, global = 2.15325e-16, cumulative = 3.88323e-10
DILUPBiCG: Solving for nuTilda, Initial residual = 1.90357e-06, Final residual = 5.79013e-12, No Iterations 1
ExecutionTime = 4890.48 s ClockTime = 5531 s

forces output:
forces(pressure, viscous)((0.00524665 0.00263487 1.15591e-19) (0.000392818 2.52438e-06 -2.2827e-23))
moment(pressure, viscous)((-1.31744e-05 2.62332e-05 -5.77975e-06) (-1.26219e-08 1.96409e-06 1.62711e-08))
local:
forces(pressure, viscous)((0.00524665 0.00263487 1.15591e-19) (0.000392818 2.52438e-06 -2.2827e-23))
moment(pressure, viscous)((-1.31744e-05 2.62332e-05 -5.77975e-06) (-1.26219e-08 1.96409e-06 1.62711e-08))

forceCoeffs output:
Cd = 1.43864
Cl = 0.672805
Cm = -0.0367569


**********

My fvSolution and fvScheme files have been attached, as well as a screenshot of the velocity contours (flow is coming in from the left at 4mps) and one showing a line plot of the velocity across the body and the yplus values.
Attached Images
File Type: jpg IDESSA-0.004sec-lineplot.jpg (25.6 KB, 16 views)
File Type: jpg IDESSA-0.004sec.jpg (64.6 KB, 23 views)
Attached Files
File Type: zip fvSchemes.zip (1.1 KB, 8 views)
File Type: zip fvSolution.zip (1.1 KB, 8 views)
jantheron 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
symmetry boundary conditions in cfx lost.identity CFX 41 May 22, 2013 07:21
Impinging Jet Boundary Conditions Anindya Main CFD Forum 24 January 11, 2012 14:40
OpenFOAM Variable Velocity Boundary Conditions NickolasPl OpenFOAM Programming & Development 2 May 19, 2011 05:37
boundary conditions and mesh exporting vaina74 Open Source Meshers: Gmsh, Netgen, CGNS, ... 2 May 27, 2010 09:38
Trimmed cell and embedded refinement mesh conversion issues michele OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 2 July 15, 2005 04:15


All times are GMT -4. The time now is 08:53.