PHASTA Wiki - User contributions [en]

PHASTA/Extract Wall Shear Stress

2012-08-04T04:34:12Z

Ychen:

[Compressible Flow]
How to extract the wall shear stress (wss) vectors from the ybar file containing time-averaged statistics or the restart file containing instantaneous solutions?

1. Go to the directory containing the reduced one-part geombc.dat.1 file and the ybar file, for example, ybar.352000-372000.0

2. Run the following command one by one

'''mkdir wss''' ''create a new directory to put the wss vtk file''

'''cd wss'''

'''ln -s ../geombc.dat.1''' ''create a softlink of the geombc.dat.1 file''

'''cp /users/ychen/develop/phasta/wallshearC/wallshearC/input.config .'''

'''cp /users/ychen/develop/phasta/wallshearC/wallshearC/solver.inp .''' ''input.config and solver.inp here have no real meanings, the only reason is that the code need to read them''

'''ln -s ../ybar.352000-372000.0 ybar.372000.0 ''' ''create a softlink of the ybar file, notice that the link only takes the final time step in the file name''

or '''ln -s ../restart.372000.0''' ''if we are extracting wss from the solution field of the restart file''

'''echo 372000 > numstart.dat'''

'''echo 1 > itagybar.dat''' ''tell the code to extract wss from the ybar file ''

'''mpirun /users/ychen/develop/phasta/wallshearC/wallshearC/bin/x86_64_linux-icc/wallshearC.exe-openmpi > dumpWSS.log 2>&1 & ''' ''run the code''

tail -f dumpWSS.log

=============================================================================================================================

[Incompressible Flow]

It is easy to do but I need some time to test it. Michele might have a working version.

'''You are strongly encouraged to build the code under your own directory, probably together with the phSolver toolsets using CMake.'''

--[[User:Ychen|Ychen]] 13:14, 26 July 2012 (MDT)

View VTK file in Matlab

2012-07-27T20:10:58Z

Ychen: Created page with "1. Convert vtk file to matlab data VTKreader Convert the vtk file of an all-triangle surface mesh to the Matlab data. The connectivity is stored in tri, coordinates are x, y, an..."

1. Convert vtk file to matlab data

VTKreader
Convert the vtk file of an all-triangle surface mesh to the Matlab data. The connectivity is stored in tri, coordinates are x, y, and z.

/users/ychen/NGC/surface.vtk

In Matlab command line type

cd /users/ychen/NGC

VTKreader('surface')

Run Inviscid Simulations

2012-07-27T19:52:09Z

Ychen: Created page with "If the inviscid case is 3-D, the meshing is actually easier. The tricky case is 2-D 1. Solver Choice For low-speed case we can use AcuSolve, to load it vglconnect -s viz001 s..."

If the inviscid case is 3-D, the meshing is actually easier. The tricky case is 2-D

1. Solver Choice

For low-speed case we can use AcuSolve, to load it

vglconnect -s viz001

soft add +acusolve

vglrun acuConsole

We can group the CAD model and meshing it. All the walls are slip-wall type

After we generate the input file for the solver we can run the solver using the command:

Category:Post-processing

2012-07-27T17:13:34Z

Ychen:

[[Collect Turbulence Statistics]]

[[Extract Wall Shear Stress]]

[[Paraview Automation with Python]]

[[Visualization using Matlab]]

[[Run Inviscid Simulations]]

Category:Post-processing

2012-07-27T17:13:27Z

Ychen:

[[Collect Turbulence Statistics]]

[[Extract Wall Shear Stress]]

[[Paraview Automation with Python]]

[[Visualization using Matlab]]

[[Run Inviscid Simulations]]

ParaView/Automation with Python

2012-07-26T21:20:17Z

Ychen:

1. Extract surface using paraview, the code will look like

from paraview.simple import *

<nowiki>#========== read in data file</nowiki>

print 'Read phasta file'

<nowiki>dataObj = PhastaReader( FileName= 'restart_one.pht' ) # you have to change the name of this pht file</nowiki>

SRF = ExtractSurface(Input=dataObj)

MGBLK1 = MergeBlocks(Input=SRF)

<nowiki>vtkWriter = DataSetWriter(FileName = 'Surface.vtk' , Input = MGBLK1, FileType = 1) # you have to change the name of this vtk file</nowiki>

vtkWriter.UpdatePipeline()

Get the python file from

/users/ychen/NGC/geometry2012/pvExtractSuface.py

ParaView/Automation with Python

2012-07-26T21:19:52Z

Ychen:

1. Extract surface using paraview, the code will look like

<nowiki>from paraview.simple import *

#========== read in data file

print 'Read phasta file'

dataObj = PhastaReader( FileName= 'restart_one.pht' ) # you have to change the name of this pht file

SRF = ExtractSurface(Input=dataObj)

MGBLK1 = MergeBlocks(Input=SRF)

vtkWriter = DataSetWriter(FileName = 'Surface.vtk' , Input = MGBLK1, FileType = 1) # you have to change the name of this vtk file

vtkWriter.UpdatePipeline()

Get the python file from

/users/ychen/NGC/geometry2012/pvExtractSuface.py</nowiki>

ParaView/Automation with Python

2012-07-26T21:09:47Z

Ychen:

1. Extract surface using paraview, the code will look like

from paraview.simple import *

<nowiki>#========== read in data file</nowiki>

print 'Read phasta file'

dataObj = PhastaReader( FileName= 'restart_one.pht' )

SRF = ExtractSurface(Input=dataObj)

MGBLK1 = MergeBlocks(Input=SRF)

vtkWriter = DataSetWriter(FileName = 'Surface.vtk' , Input = MGBLK1, FileType = 1)

vtkWriter.UpdatePipeline()

ParaView/Automation with Python

2012-07-26T21:00:02Z

Ychen:

1. Extract surface using paraview, the code will look like

from paraview.simple import *

\#========== read in data file

print 'Read phasta file'

dataObj = PhastaReader( FileName= 'restart_one.pht' )

SRF = ExtractSurface(Input=dataObj)

MGBLK1 = MergeBlocks(Input=SRF)

vtkWriter = DataSetWriter(FileName = 'Surface.vtk' , Input = MGBLK1, FileType = 1)

vtkWriter.UpdatePipeline()

ParaView/Automation with Python

2012-07-26T20:59:51Z

Ychen:

1. Extract surface using paraview, the code will look like

from paraview.simple import *

/#========== read in data file

print 'Read phasta file'

dataObj = PhastaReader( FileName= 'restart_one.pht' )

SRF = ExtractSurface(Input=dataObj)

MGBLK1 = MergeBlocks(Input=SRF)

vtkWriter = DataSetWriter(FileName = 'Surface.vtk' , Input = MGBLK1, FileType = 1)

vtkWriter.UpdatePipeline()

ParaView/Automation with Python

2012-07-26T20:59:33Z

Ychen:

1. Extract surface using paraview, the code will look like

''from paraview.simple import *

#========== read in data file

print 'Read phasta file'

dataObj = PhastaReader( FileName= 'restart_one.pht' )

SRF = ExtractSurface(Input=dataObj)

MGBLK1 = MergeBlocks(Input=SRF)

vtkWriter = DataSetWriter(FileName = 'Surface.vtk' , Input = MGBLK1, FileType = 1)

vtkWriter.UpdatePipeline()
''

ParaView/Automation with Python

2012-07-26T20:59:00Z

Ychen: Created page with "1. Extract surface using paraview from paraview.simple import * #========== read in data file print 'Read phasta file' dataObj = PhastaReader( FileName= 'restart_one.pht' ) ..."

1. Extract surface using paraview

from paraview.simple import *

#========== read in data file

print 'Read phasta file'

dataObj = PhastaReader( FileName= 'restart_one.pht' )

SRF = ExtractSurface(Input=dataObj)

MGBLK1 = MergeBlocks(Input=SRF)

vtkWriter = DataSetWriter(FileName = 'Surface.vtk' , Input = MGBLK1, FileType = 1)

vtkWriter.UpdatePipeline()

PHASTA/Extract Wall Shear Stress

2012-07-26T19:30:01Z

Ychen:

[Compressible Flow]
How to extract the wall shear stress (wss) vectors from the ybar file containing time-averaged statistics or the restart file containing instantaneous solutions?

1. Go to the directory containing the reduced one-part geombc.dat.1 file and the ybar file, for example, ybar.352000-372000.0

2. Run the following command one by one

'''mkdir wss''' ''create a new directory to put the wss vtk file''

'''cd wss'''

'''ln -s ../geombc.dat.1''' ''create a softlink of the geombc.dat.1 file''

'''cp ~/develop/phasta/wallshearC/wallshearC/input.config .'''

'''cp ~/develop/phasta/wallshearC/wallshearC/solver.inp .''' ''input.config and solver.inp here have no real meanings, the only reason is that the code need to read them''

'''ln -s ../ybar.352000-372000.0 ybar.372000.0 ''' ''create a softlink of the ybar file, notice that the link only takes the final time step in the file name''

or '''ln -s ../restart.372000.0''' ''if we are extracting wss from the solution field of the restart file''

'''echo 372000 > numstart.dat'''

'''echo 1 > itagybar.dat''' ''tell the code to extract wss from the ybar file ''

'''mpirun ~/develop/phasta/wallshearC/wallshearC/bin/x86_64_linux-icc/wallshearC.exe-openmpi > dumpWSS.log 2>&1 & ''' ''run the code''

tail -f dumpWSS.log

=============================================================================================================================

[Incompressible Flow]

It is easy to do but I need some time to test it. Michele might have a working version.

'''You are strongly encouraged to build the code under your own directory, probably together with the phSolver toolsets using CMake.'''

--[[User:Ychen|Ychen]] 13:14, 26 July 2012 (MDT)

PHASTA/Extract Wall Shear Stress

2012-07-26T19:15:22Z

Ychen:

[Compressible Flow]
How to extract the wall shear stress (wss) vectors from the ybar file containing time-averaged statistics or the restart file containing instantaneous solutions?

1. Go to the directory containing the reduced one-part geombc.dat.1 file and the ybar file, for example, ybar.352000-372000.0

2. Run the following command one by one

'''mkdir wss''' ''create a new directory to put the wss vtk file''

'''cd wss'''

'''ln -s ../geombc.dat.1''' ''create a softlink of the geombc.dat.1 file''

'''cp ~/develop/phasta/wallshearC/wallshearC/input.config .'''

'''cp ~/develop/phasta/wallshearC/wallshearC/solver.inp .''' ''input.config and solver.inp here have no real meanings, the only reason is that the code need to read them''

'''ln -s ../ybar.352000-372000.0 ybar.372000.0 ''' ''create a softlink of the ybar file, notice that the link only takes the final time step in the file name''

'''echo 372000 > numstart.dat'''

'''echo 1 > itagybar.dat''' ''tell the code to extract wss from the ybar file ''

'''mpirun ~/develop/phasta/wallshearC/wallshearC/bin/x86_64_linux-icc/wallshearC.exe-openmpi > dumpWSS.log 2>&1 & ''' ''run the code''

tail -f dumpWSS.log

=============================================================================================================================

[Incompressible Flow]

It is easy to do but I need some time to test it. Michele might have a working version.

'''You are strongly encouraged to build the code under your own directory, probably together with the phSolver toolsets using CMake.'''

--[[User:Ychen|Ychen]] 13:14, 26 July 2012 (MDT)

PHASTA/Extract Wall Shear Stress

2012-07-26T19:14:37Z

Ychen:

[Compressible Flow]
How to extract the wall shear stress (wss) vectors from the ybar file containing time-averaged statistics or the restart file containing instantaneous solutions?

1. Go to the directory containing the reduced one-part geombc.dat.1 file and the ybar file, for example, ybar.352000-372000.0

2. Run the following command one by one

'''mkdir wss''' ''create a new directory to put the wss vtk file''

'''cd wss'''

'''ln -s ../geombc.dat.1''' ''create a softlink of the geombc.dat.1 file''

'''cp ~/develop/phasta/wallshearC/wallshearC/input.config .'''

'''cp ~/develop/phasta/wallshearC/wallshearC/solver.inp .''' ''input.config and solver.inp here have no real meanings, the only reason is that the code need to read them''

'''ln -s ../ybar.352000-372000.0 ybar.372000.0 ''' ''create a softlink of the ybar file, notice that the link only takes the final time step in the file name''

'''echo 372000 > numstart.dat'''

'''echo 1 > itagybar.dat''' ''tell the code to extract wss from the ybar file ''

'''mpirun ~/develop/phasta/wallshearC/wallshearC/bin/x86_64_linux-icc/wallshearC.exe-openmpi > dumpWSS.log 2>&1 & ''' ''run the code''

tail -f dumpWSS.log

[Incompressible Flow]
It is easy to do but I need some time to test it. Michele might have a working version.

'''You are strongly encouraged to build the code under your own directory, probably together with the phSolver toolsets using CMake.'''

--[[User:Ychen|Ychen]] 13:14, 26 July 2012 (MDT)

PHASTA/Extract Wall Shear Stress

2012-07-26T19:11:42Z

Ychen:

How to extract the wall shear stress (wss) vectors from the ybar file containing time-averaged statistics or the restart file containing instantaneous solutions?

1. Go to the directory containing the reduced one-part geombc.dat.1 file and the ybar file, for example, ybar.352000-372000.0

2. Run the following command one by one

'''mkdir wss''' ''create a new directory to put the wss vtk file''

'''cd wss'''

'''ln -s ../geombc.dat.1''' ''create a softlink of the geombc.dat.1 file''

'''cp ~/develop/phasta/wallshearC/wallshearC/input.config .'''

'''cp ~/develop/phasta/wallshearC/wallshearC/solver.inp .''' ''input.config and solver.inp here have no real meanings, the only reason is that the code need to read them''

'''ln -s ../ybar.352000-372000.0 ybar.372000.0 ''' ''create a softlink of the ybar file, notice that the link only takes the final time step in the file name''

'''echo 372000 > numstart.dat'''

'''echo 1 > itagybar.dat''' ''tell the code to extract wss from the ybar file ''

'''mpirun ~/develop/phasta/wallshearC/wallshearC/bin/x86_64_linux-icc/wallshearC.exe-openmpi > dumpWSS.log 2>&1 & ''' ''run the code''

tail -f dumpWSS.log

PHASTA/Extract Wall Shear Stress

2012-07-26T19:05:52Z

Ychen:

How to extract the wall shear stress (wss) vectors from the ybar file containing time-averaged statistics or the restart file containing instantaneous solutions?

1. Go to the directory containing the reduced one-part geombc.dat.1 file and the ybar file, for example, ybar.352000-372000.0

2. Run the following command one by one

'''mkdir wss''' ''create a new directory to put the wss vtk file''

'''cd wss'''

'''ln -s ../geombc.dat.1''' ''create a softlink of the geombc.dat.1 file''

'''cp ~/develop/phasta/wallshearC/wallshearC/input.config .'''

'''cp ~/develop/phasta/wallshearC/wallshearC/solver.inp .''' ''input.config and solver.inp here have no real meanings, the only reason is that the code need to read them''

'''ln -s ../ybar.352000-372000.0 ybar.372000.0 ''' ''create a softlink of the ybar file, notice that the link only takes the final time step in the file name''

'''echo 372000 > numstart.dat'''

'''ehco 1 > itagybar.dat''' ''tell the code to extract wss from the ybar file ''

'''mpirun ~/develop/phasta/wallshearC/wallshearC/bin/x86_64_linux-icc/wallshearC.exe-openmpi > dumpWSS.log 2>&1 & ''' ''run the code''

tail -f dumpWSS.log

PHASTA/Extract Wall Shear Stress

2012-07-26T19:01:04Z

Ychen: Created page with "How to extract the wall shear stress vectors from the ybar file containing time-averaged statistics or the restart file containing instantaneous solutions? 1. Go to the director..."

How to extract the wall shear stress vectors from the ybar file containing time-averaged statistics or the restart file containing instantaneous solutions?

1. Go to the directory containing the reduced one-part geombc.dat.1 file and the ybar file, for example, ybar.352000-372000.0

2. Run the following command one by one

mkdir wss

cd wss

ln -s ../geombc.dat.1

cp ~/develop/phasta/wallshearC/wallshearC/input.config .

cp ~/develop/phasta/wallshearC/wallshearC/solver.inp .

ln -s ../ybar.352000-372000.0 ybar.372000.0

echo 372000 > numstart.dat

ehco 1 > itagybar.dat

mpirun ~/develop/phasta/wallshearC/wallshearC/bin/x86_64_linux-icc/wallshearC.exe-openmpi > dumpWSS.log 2>&1 &

tail -f dumpWSS.log

Category:Post-processing

2012-07-26T18:57:44Z

Ychen: Created page with "Collect Turbulence Statistics Extract Wall Shear Stress Paraview Automation with Python Visualization using Matlab"

[[Collect Turbulence Statistics]]

[[Extract Wall Shear Stress]]

[[Paraview Automation with Python]]

[[Visualization using Matlab]]

Main Page

2012-07-26T18:53:24Z

Ychen:

Welcome to the Wiki for the [https://redmine.scorec.rpi.edu/projects/phasta PHASTA] project at the University of Colorado, Boulder. The information contained within is mostly systems documentation and usage information for software which may or may not be public yet, but feel free to look around.

[http://fluid.colorado.edu/wiki/index.php/Special:AllPages List Of All Pages]

[[:Category:Software|Available Software]]

[[:Category:Software Engineering|Software Engineering]]

[[UNIX|Getting Started with UNIX]]

[[:Category:Post-processing|Post-processing and Visualization]]