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

non uniform initial conditions

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 19, 2014, 15:09
Default non uniform initial conditions
  #1
New Member
 
Join Date: Feb 2014
Posts: 12
Rep Power: 3
ajnewman is on a distinguished road
Hello all,

Simple (?) question: If I want to specify a nonuniform initial condition using the internalField nonuniform List<vector> approach, like in the boxTurb tutorial, do I specify my initial velocity at all the points listed in the "points" file and in that order? That is, if the first entry in points is (0,0,0) should the first entry in U be the velocity at (0,0,0)? This simulation has multiple blocks.

I have seen the setFields as used in the dam break tutorial but am not currently savvy enough to use this to do what I want.

The background is I am attempting to run a simulation in OpenFOAM starting with an initial condition generated by a fortran code on a mesh that is not the same as my OpenFOAM mesh.

Any input is appreciated.

-Jensen
ajnewman is offline   Reply With Quote

Old   February 19, 2014, 15:49
Default
  #2
New Member
 
Join Date: Feb 2014
Posts: 12
Rep Power: 3
ajnewman is on a distinguished road
I think I answered a bit of my own question since the points file has (no_cells_x+1)*(no_cells_y+1)*(no_cells_z+1) entries whereas a nonuniform velocity entry should have only (no_cells_x)*(no_cells_y)*(no_cells_z) entries. Based on this I am guessing that velocity needs to be spec'd at cell centers.
ajnewman is offline   Reply With Quote

Old   February 19, 2014, 18:03
Default
  #3
Member
 
Jace
Join Date: Oct 2012
Posts: 77
Rep Power: 5
zhengzh5 is on a distinguished road
Quote:
Originally Posted by ajnewman View Post
I think I answered a bit of my own question since the points file has (no_cells_x+1)*(no_cells_y+1)*(no_cells_z+1) entries whereas a nonuniform velocity entry should have only (no_cells_x)*(no_cells_y)*(no_cells_z) entries. Based on this I am guessing that velocity needs to be spec'd at cell centers.
you're right on. it's finite-volume approach. So all volScalarFields/volVectorFields are specified at the cell-center. Surface and point values are usually interpolated based on cell-center values.
zhengzh5 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
a problem with convergence in buoyantSimpleFoam skuznet OpenFOAM Running, Solving & CFD 5 February 19, 2014 05:30
Extremely slow simulation with interDyMFoam jrrygg OpenFOAM Running, Solving & CFD 9 April 23, 2013 10:14
pisoFoam with k-epsilon turb blows up - Some questions Heroic OpenFOAM Running, Solving & CFD 26 December 17, 2012 04:34
High Courant Number @ icoFoam Artex85 OpenFOAM Running, Solving & CFD 9 January 3, 2012 09:06
Velocity blows up suddenly after 30,000+ iterations lordvon OpenFOAM 14 October 27, 2010 11:13


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