extracting velocity field from 0/U file into matrix format

eaglemckenna · November 3, 2014, 11:17

Hi

Does anyone know a methodology or is there a command or scrip in OF to extract the velocity components from the final time step U file into column matrix format? It is not feasible to do this manually due to the number of rows as well as the parenthesis. I intend to perform numerical analysis on the data.

thanks
g

ssss · November 3, 2014, 14:01

Maybe this helps but it's not a decent way of doing it:

1) Go to paraView, last time step
2) Go to file save data
3) Save it in csv
4) Open the csv you will see some columns with {U:0} {U:1} {U:2} this are the velocity components.

Hope it helps.

eaglemckenna · November 4, 2014, 06:43

Thanks so much
Apologies for such an obvious question.

eaglemckenna · November 5, 2014, 11:29

Greeting fellow CFD users.

I have created an algorithm to determine POD velocity basis vector for x and y components at each node point for a 2d mesh for 9 observations with Matlab.
I have a column vector x denoting the x co-ordinates at each node and a column vector y for the y co-ordinates.
plot(x,y) gives me the 2d domain.
I have a column vector r denoting the resulting velocity of x and y at each node.
Now all I need to do is map the resultant velocity vector at each node onto the back onto the 2d grid.
What is the command to do this?
I have tried meshgrid, and ndgrid but have no success as of yet.
kind regards
gowain

November 3, 2014, 11:17	extracting velocity field from 0/U file into matrix format	#1
eaglemckenna Member Gowain Join Date: May 2014 Location: Ireland Posts: 86 Rep Power: 11	Hi Does anyone know a methodology or is there a command or scrip in OF to extract the velocity components from the final time step U file into column matrix format? It is not feasible to do this manually due to the number of rows as well as the parenthesis. I intend to perform numerical analysis on the data. thanks g

November 3, 2014, 14:01		#2
ssss Senior Member anonymous Join Date: Aug 2014 Posts: 205 Rep Power: 12	Maybe this helps but it's not a decent way of doing it: 1) Go to paraView, last time step 2) Go to file save data 3) Save it in csv 4) Open the csv you will see some columns with {U:0} {U:1} {U:2} this are the velocity components. Hope it helps. eaglemckenna likes this.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
SparceImage v1.7.x Issue on MAC OS X	rcarmi	OpenFOAM Installation	4	August 14, 2014 06:42
"parabolicVelocity" in OpenFoam 2.1.0 ?	sawyer86	OpenFOAM Running, Solving & CFD	21	February 7, 2012 11:44
Version 15 on Mac OS X	gschaider	OpenFOAM Installation	113	December 2, 2009 10:23
DecomposePar links against liblamso0 with OpenMPI	jens_klostermann	OpenFOAM Bugs	11	June 28, 2007 17:51
[blockMesh] Axisymmetrical mesh	Rasmus Gjesing (Gjesing)	OpenFOAM Meshing & Mesh Conversion	10	April 2, 2007 14:00

November 4, 2014, 06:43		#3
eaglemckenna Member Gowain Join Date: May 2014 Location: Ireland Posts: 86 Rep Power: 11	Thanks so much Apologies for such an obvious question.

November 5, 2014, 11:29		#4
eaglemckenna Member Gowain Join Date: May 2014 Location: Ireland Posts: 86 Rep Power: 11	Greeting fellow CFD users. I have created an algorithm to determine POD velocity basis vector for x and y components at each node point for a 2d mesh for 9 observations with Matlab. I have a column vector x denoting the x co-ordinates at each node and a column vector y for the y co-ordinates. plot(x,y) gives me the 2d domain. I have a column vector r denoting the resulting velocity of x and y at each node. Now all I need to do is map the resultant velocity vector at each node onto the back onto the 2d grid. What is the command to do this? I have tried meshgrid, and ndgrid but have no success as of yet. kind regards gowain