Re: [ff3d-users] Navier-Stokes problem

[DELPINO Stephane]

Hello. I just reminded that I made a mistake.

Stephane Del Pino a écrit :
> Yes. The force on the object can be approximated by the integral over the 
> object of 1/epsilon*U :
>         double fx = int[M](1/epsilon*ux);
>         double fy = int[M](1/epsilon*uy);
>         double fz = int[M](1/epsilon*uz);
> where U = (ux,uy,uz) and the epsilon is used for the penalty term P in the 
> example for Dominique, 1/epsilon = 10^3.
>   
One should have read

double fx = int[M](ki/epsilon*ux);

where ki is a function that is 0 outside the obstacle and 1 inside. For 
instance:
   function ki = one(<1,0,0>);
where <1,0,0> is the reference of the object.

Regards,
Stephane.