how to probe a single velocity component using probeLocations

gasagna · August 15, 2014, 08:05

Hi,

I am using with success the utility probeLocations to probe the pressure field at selected points. Probing the velocity field also works in general, but I need only one of the velocity components. I could parse the generated file for post processing, but it would be more simple if i only get one velocity component only, as the resulting file would be simpler.

Is this possible?

Thanks,

Davide

samieh · August 19, 2014, 11:53

hi
i am trying to do the same job.would you please tell me how you did that? i want to calculate the pressure on some points on my column but i do not know the file which should be added in contdic
thanks alot

Pravin Kadu · January 29, 2018, 01:49

You can write your own solver by adding the fields . I did it for me.

To create fields (createFields.H):
Info<< "Reading field Uxvel\n" <<endl;
volScalarField Uxvel
(
IOobject
(
"Uxvel",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
);

Info<< "Reading field Uyvel\n" <<endl;
volScalarField Uyvel
(
IOobject
(
"Uyvel",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
);

Info<< "Reading field Uzvel\n" <<endl;
volScalarField Uzvel
(
IOobject
(
"Uzvel",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
);

For solver (solver is renamed - pimpleFoamComponentsCalc.C)-
Uxvel = U.component(0);
Uyvel = U.component(1);
Uzvel = U.component(2);

August 15, 2014, 08:05	how to probe a single velocity component using probeLocations	#1
gasagna New Member Davide Lasagna Join Date: Jul 2014 Posts: 3 Rep Power: 11	Hi, I am using with success the utility probeLocations to probe the pressure field at selected points. Probing the velocity field also works in general, but I need only one of the velocity components. I could parse the generated file for post processing, but it would be more simple if i only get one velocity component only, as the resulting file would be simpler. Is this possible? Thanks, Davide

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Axisymmetric Case 3rd Velocity Component	Mentalo	OpenFOAM Running, Solving & CFD	2	June 8, 2014 12:45
Suspected division by zero on velocity reported by point probe	dahicke	STAR-CCM+	1	July 4, 2013 02:16
Stn frame and radial velocity component in NREL	mohammad	CFX	1	October 25, 2011 10:13
sampling a strain-rate tensor single component	vkrastev	OpenFOAM Post-Processing	0	March 2, 2010 07:21
Terrible Mistake In Fluid Dynamics History	Abhi	Main CFD Forum	12	July 8, 2002 09:11

August 19, 2014, 11:53		#2
samieh New Member Ontario Join Date: Aug 2014 Posts: 2 Rep Power: 0	hi i am trying to do the same job.would you please tell me how you did that? i want to calculate the pressure on some points on my column but i do not know the file which should be added in contdic thanks alot

January 29, 2018, 01:49		#3
Pravin Kadu New Member Kadu Join Date: Jun 2015 Location: India Posts: 6 Rep Power: 10	You can write your own solver by adding the fields . I did it for me. To create fields (createFields.H): Info<< "Reading field Uxvel\n" <<endl; volScalarField Uxvel ( IOobject ( "Uxvel", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh ); Info<< "Reading field Uyvel\n" <<endl; volScalarField Uyvel ( IOobject ( "Uyvel", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh ); Info<< "Reading field Uzvel\n" <<endl; volScalarField Uzvel ( IOobject ( "Uzvel", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh ); For solver (solver is renamed - pimpleFoamComponentsCalc.C)- Uxvel = U.component(0); Uyvel = U.component(1); Uzvel = U.component(2);