how to find the maximum vorticity and streamfunction in driven caivty
 

Dear friends:
i solved the 2D lid driven cavity problem using SIMPLE algorithm. could you please tell me how to obtain the maximum vorticity and stream function in the flow field.
thanks

did you find the answer?

what is the meaning?

At steady state, compute the vorticity from your velocity field and solve the elliptic equation Lap (Psi) = vorticity. This way you get both functions and all the values you need

thank you . i can obtain the whole field of them. but my problem is to find the maximum values so as to compare it with the benchmark solutions in the public paper. Since we only have the field on grid points. then how to find the maximum value?

Let's imagine your are in 2D and vorticity is stored in an array VOR(NI,NJ).

VORMAX=VOR(1,1)

DO J=1,NJ
DO I=1,NI

IF ( VOR(I,J).GT.VORMAX ) VORMAX = VOR(I,J)

END DO
END DO

assume fortran coding

thank you , but i can not fully agree with you. in this way, the maximum value will depend on the mesh density.

Establish a grid-independent solution. Otherwise any approach you follow will depend on the grid, whether you use an interpolation method or simply find the max over the grid.

But even that will be grid dependent if the underlying solution is grid dependent.

leflix is right, just add the i,j values where the maximum is obtained.
What you ask for is simply wrong. In a numerical solution nothing else exists but the grid value. This is established by the Nyquist theorem.
When you want to compute some extra-grid value, is common to use some interpolation to reconstruct a continous function. But the values you obtain are somehow arbitrarily, depending on the order of the polynomial reconstruction.
Therefore, just use the grid value, they are meaningfull of your solution

Hi Ztdep
why you say that it will be grid dependent?if the main solution be grid independent,then vorticity values will be independent and you can use the code like a postProcessing utility for obtaining maximum vortivity(at least by theory).

open foam
 

hi everybody
can you help me about openFoam ?please
I have run a program but the software shows some error .if you can help me I will send errors for you.:)

hi
whats your problem?which university are you?

openFoam
 

hi I am in razi university from kermanshah.
I am working with openfoam I am going to analyse ac cavity .my solver is boussinesqsimple .when I run the program there is some error in fvslution(include solver for pressuer & tempreture & velocity and SIMPLE solver) file .I have copied this file from another exaxmple without any changing.
I should be say that the software reads constant & 0 files.