Difference between revisions of "Creating geom.spj"

From PHASTA Wiki
Jump to: navigation, search
Line 44: Line 44:
  
 
Initial Scalar_1: Set to 1.8e-7
 
Initial Scalar_1: Set to 1.8e-7
 +
 +
=== NoSlip Walls ===
 +
comp3: Set the no-slip walls with comp3. For Boeing, the the only walls that do not get no-slip are the tunnel walls. Therefore, what you should do is go to Selection -> "Select All Model Faces" and then create a comp3. The vector magnitude is 0, and the direction '''must be non-zero''' (i.e. [1 0 0]).

Revision as of 11:22, 15 July 2013


  • create a directory for the smf and spj files (mkdir SMF_and_SPJ_files). cd to that directory and start SimulationMaker:
 vglrun /usr/local/simapps/sdk-090321/simapps  SimulationMaker
  • Default Solver: Phasta
  • Browse to the appropriate geometry file and select it (make sure it is not the SpaceClaim parasolid file; it should be the file you saved from ThreeDViewer - geom.xmt_txt)
  • Turn OFF sliver feature supression
  • Click Define Simulation Model Instance
  • The attribute list and model will appear in split windows. At the top of the list, double click on the name "problem definition" and set the name to "geom".
  • To set attributes, make sure the field with "geom", which you just changed, is selected. Using the model to the right, click on whichever faces you want to give attributes, and go to "Create" at the upper left hand corner to select the appropriate attribute.
  • The names of the attributes are arbitrary, but should be something meaningful to the user and whoever else is involved in the project.

For Incompressible Flow

Inflow

comp3: The name is arbitrary but should be set to something meaningful (i.e. Inflow). Click on the inflow face (which is the face which makes up the front of the wind tunnel) and go to Create -> comp3. Set vector magnitude to the free-stream velocity you want. Set the vector direction to the direction (i.e. if the free-stream velocity moves strictly along the x-axis, the direction would be [1 0 0]). Even if the vector magnitude is 0, the direction CANNOT be [0 0 0]. This will cause Phasta to crash. It must have at least one non-zero component. In this case, just set it to [1 0 0].

scalar_1: Click on the inflow face and go to Create -> scalar_1. Set the field "Type: scalar_1" to 1.8e-7. This is the eddy viscosity at the inlet.

Outflow

natural pressure: Set "Type: natural pressure" to 0.

scalar_1 flux: Set this to 0.

traction vector: Set this to [0 0 0]

TopBottomWalls

  • These are the walls of the wind tunnel that bound the model from above and below. The terms "above" and "below" are arbitrary, but we typically define them so that their normal points along the z-axis. However, this is not necessary.

scalar_1 flux: Set this to 0.

traction vector: Set this to [0 0 0].

comp1: The vector magnitude should be 0. The direction depends on which direction the normal of the wall points. If it is along the z-axis, then this will be [0 0 1]. If it is along the y-axis, this will be [0 1 0]. This attribute signifies a no-penetration condition through the tunnel walls.

SideWalls

  • These are also arbitrary. But if you already defined the TopBottomWalls, then you already know which walls they are.

scalar_1 flux: Set this to 0

traction vector: Set this to [0 0 0]

comp1: Set vector magnitude to 0. The direction should be set to whichever direction the normal of the walls point.

Initial Conditions

Initial Velocity: In the model window, go to Selection -> Select Model. In the attribute window, go to Create -> initial velocity. Set this to [1e-8 0 0]

Initial Pressure: Set to 0

Initial Scalar_1: Set to 1.8e-7

NoSlip Walls

comp3: Set the no-slip walls with comp3. For Boeing, the the only walls that do not get no-slip are the tunnel walls. Therefore, what you should do is go to Selection -> "Select All Model Faces" and then create a comp3. The vector magnitude is 0, and the direction must be non-zero (i.e. [1 0 0]).