CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   CD-adapco (http://www.cfd-online.com/Forums/cd-adapco/)
-   -   velocity divergence as plot (http://www.cfd-online.com/Forums/cd-adapco/52704-velocity-divergence-plot.html)

Sabatini October 12, 2001 04:28

velocity divergence as plot
 
Hallo people Is there a way to extract the velocity divergence as plot Thanks

Paal October 16, 2001 06:56

Re: velocity divergence as plot
 
Hi,

I am not fully sure that I get your question right: Do you want to plot the diverged "results"?

These results are saved to a file - case.div, and can be loaded in the same manner as normal result files just replacing the case.pst with case.div


steve October 16, 2001 13:11

Re: velocity divergence as plot
 
I think to get what you want you have to load velocities (GETV,ALL) and then use OPERATE,CURL function. I think CURL is the same as divergence of a function?

Lars Ola Liavåg October 17, 2001 02:49

Re: velocity divergence as plot
 
The curl of a vector field is itself a vector (cross product of the "upside down" delta and the vector), whereas the divergence (the corresponding dot product) is a scalar. Hence, for the velocity vector V, div V = du/dx + dv/dy + dw/dz, whereas curl V = i(dw/dy - du/dz) + j(du/dz - dw/dx) + k(dv/dx - du/dy), all d's being partial differentials.

In PROSTAR, the OPERATE,CURL function takes the curl of the loaded vector field and returns the curl (vector) of it. You can also calculate the gradient vector, but I don't know how to calculate the divergence in PROSTAR. I guess contacting the support team is a good idea.


daniel January 21, 2002 10:30

Re: velocity divergence as plot
 
CURL is definitely NOT the same as the divergence of a function, but a DIV command should also be available the same way

cfd January 28, 2002 20:07

Re: velocity divergence as plot
 
Continuity equation:

div(rho u) = d(rho)/dt

For a converged steady state solution:

div(u) -> 0.

For a transient compressible time step solution:

div(u) = 1/rho [d(rho)/dt - u grad(rho)]

For incompressible time step solution: (rho = const.):

div(u) -> 0.


All times are GMT -4. The time now is 23:55.