|
|
Blitz Devel : |
From: Julian Cummings (cummings_at_[hidden])
Date: 2002-11-26 18:10:46
Hi Martin,
I looked at your "fix" for the cfd example, but I don't understand it.
You have added the divergence of the external force field to the
computation of the rhs of the pressure equation. But the external
force field is already being accounted for during the velocity update.
nextV = *V + delta_t * ( recip_rho * (
eta * Laplacian3DVec4(V,geom) - grad3D4(P, geom) + *force) - *advect);
So I don't think the force belongs in the pressure computation at all.
I am curious as to what the "bug" is that you were trying to fix.
Regards, Julian C.
Martin Verlaan en Hilde Poodt wrote:
> L.S.,
>
> After playing around with your nice library and cfd example, I think
> I may have found an answer to a bug in it. I enclosed the modified
> example with this e-mail.
>
> Best regards,
>
> Martin Verlaan
>
> -------------------------------------------------------
>
> --------------------------------------------------------------------------------
> Name: cfd.cpp
> cfd.cpp Type: text/x-c
> Encoding: base64
-- Dr. Julian C. Cummings E-mail: cummings_at_[hidden] California Institute of Technology Phone: 626-395-2543 1200 E. California Blvd., Mail Code 158-79 Fax: 626-584-5917 Pasadena, CA 91125