[anja]

Is it possible with these lines:

scalar dx = x - Cp*vector(1,0,0);
scalar dy = y - Cp*vector(0,1,0);
scalar dz = z - Cp*vector(0,0,1);
scalar d = sqrt(dx*dx+dy*dy+dz*dz);

to find the distance of each each cell to the centre of the inlet?

Or is this basically not the way to do that?

Thanks
Anja

[eugene]

There is almost always a better way to do things in OF.

scalarField d = mag(mesh.C() - Cp);

[anja]

Hi,

Thanks, it seems to work, or at least there was no error reported while compiling.

The problem is that I can't use my application yet. When I try the way Fabian was describing:
Usage: ssimpleFoam <root> <case> [-parallel]

ERROR: wrong number of arguments, expected 2 found 3

So what do I have to do that I can use my velocity profil for an inlet?

Thanks
Anja

[mirsaeedi]

Hi fabian

when i used your ssimpleFoam and wmake it i recieved below error:

Foam::error::printStack(Foam:stream&)
Foam::error::abort()
Foam::operator+(Foam::dimensionSet const&, Foam::dimensionSet const&)
Foam::tmp<foam::geometricfield<double,> > Foam::operator+<double,>(Foam::tmp<foam::geometric field<double,> > const&, Foam::tmp<foam::geometricfield<double,> > const&)
mySpecialFoam [0x80c551b]
mySpecialFoam [0x805b769]
__libc_start_main
__gxx_personality_v0
Aborted


i think it is from line " scalarField y=centre.component(vector::Y);"

Do you have any suggestion?

Thanx
Misaeedi.

[mirsaeedi]

I have another question?

I can not register in workshop of openfoam and i am new with openfoam and c++.
My MSC thesis is about solving Supercavitation phenomena with LES method.

is openFoam an apropriate tool to do this thesis?if it is, can i learn openfoam and it's programming on myself. ofcourse i have read the Userguide and programmer's Guide but i think it is not suffiecient for me and my thesis.

thanks
Morteza

[aldo.iannetti]

Hi
I need to use (in my OF 1.6) a velocity profile coming from an external file (reading and interpolate points)
Is there a tool already implemented in my version?
(I'm a new user)
Thanks

Originally Posted by francois

Thank you Hrvoje for this precision !

I will try to make a utility that reads a file with an inlet velocity profile in it and interpolate it on the mesh.

In your previous mail you mentioned a class that can do interpolation for profiles:

src/sampling/interpolateXY/interpolateXY.H

I will try to use it for my utility.

Thank's for taking some of your time to answer questions from newbies !

Francois

[maddalena]

Originally Posted by aldo.iannetti

Hi
I need to use (in my OF 1.6) a velocity profile coming from an external file (reading and interpolate points)
Is there a tool already implemented in my version?
(I'm a new user)
Thanks

Hello Aldo and welcome!
Search for the timeVaryingMappedFixedValue boundary condition, it can do what you are looking for.
Enjoy,

mad

[aldo.iannetti]

Hi Mad and thanks,
Correct me if I'm wrong,timeVaryingMappedFixedValue is for transient bc, I need a spatial varing bc in a way that depends on a spatial interpolation of a number of points given by the user trough an external .dat file.
Do you mean I can use timeVaryingMappedFixedValue modifying it? How?

Aldo

Originally Posted by maddalena

Hello Aldo and welcome!
Search for the timeVaryingMappedFixedValue boundary condition, it can do what you are looking for.
Enjoy,

mad

[maddalena]

Well,
you can use timeVaryingMappedFixedValue with only one sample time, thus it is not time varying! What I usually do is sample data from a OF converged case on an OF format file, adjust their position manually and apply timeVaryingMappedFixedValue bc for velocity in a different case (a subcase). However, I guess it can be used with tests sample data in the same way...
Hope it helps,
cheers,

mad

[maysmech]

Dear Foamers,
I need to define Driver and Seegmiller velocity inlet profile for backward facing step.
How can i write it in 0/U flie?

[xiyuqiu]

Hello, is it possible to generate a velocity profile by simply running a development flow case and use the result for the new case's initial condition?

[maddalena]

Originally Posted by xiyuqiu

Hello, is it possible to generate a velocity profile by simply running a development flow case and use the result for the new case's initial condition?

Of course it is!

mad