[liugx212]

Dear all:

I want to modify the x, y and z value of U. What do I do? Could you please give help?

Thank you very much,
xiang

[alberto]

If you mean the U components of the initial condition, open the

/case_name/0

directory and edit the U file according to your needs.

Alberto

[alberto]

Sorry, read "initial condition" as "initial or baoundary conditions".

[liugx212]

Dear Alberto Passalacqua:

Thanks a lot.
I have read the initial value. I mean I want to change them after read in icoFoam.C. And I have add the code. But some errors is here when I run it. Could you please give more help.

code:
// modifying the value of the U
volScalarField Ux = U.component(vector::X);
volScalarField Uy = U.component(vector::Y);
volScalarField Uz = U.component(vector::Z);

const vectorField& centers = mesh.C().internalField();
forAll (centers, celli)
{
double angle = 180;
Ux = Ux*Foam::sin(angle);
Uy = Uy*Foam::cos(angle);
Uz = 0;
}

error:

--> FOAM FATAL ERROR : Different dimensions for =
dimensions : [0 1 -1 0 0 0 0] = [0 0 0 0 0 0 0]


From function dimensionSet::operator=(const dimensionSet& ds) const
in file dimensionSet/dimensionSet.C at line 158.

FOAM aborting

Aborted


Thank you very much.
xiang

[niklas]

There are two major mistakes in the above code.
First, you are working on a copy of the U field.
what you probably want is to create a reference.
You should then use volScalarField& instead.

Second, you are looping over the cells which probably means that you want to get the cell values,i.e
Ux[celli] = U[celli].component(vector::X) etc...

third, Xian you have posted an identical question on another topic (under a really useles topic since it applies on every post in here), this is REALLY bad. Now we have 2 forums to talk about the same thing and you have done this before: I am begging you on my knees.
Please dont do it again or you will go on my ignore list.

N

[liugx212]

Dear Niklas Nordin:

Thanks a lot.

You are right. But I only worry nobody to see. Since somebody were noticed, I surely do as you said.

Thanks again.