Hi Kiny, On 14/11/2019 06:11, Kiny Wan wrote: > I found the gradient of velocity in fluidforce.py, expressed as > | > # Gradient of velocity > | > gradu = (gradrhou - gradrho[:, None]*u[None, 1:-1]/rho) / rho > > Is this form right? I think that the "rho" should be deleted as > | > gradu =(gradrhou -gradrho[:,None]*u[None,1:-1])/rho
In 3D:
u[None, 1:-1] = [rho*u, rho*v, rho*w] (the momenta)
and thus we divide through by rho here to get the actual velocities
which are required as per the chain rule.
Regards, Freddie.
--
You received this message because you are subscribed to the Google Groups "PyFR
Mailing List" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To view this discussion on the web, visit
https://groups.google.com/d/msgid/pyfrmailinglist/162ec720-37eb-5dbc-53af-7e2bf2554ecf%40witherden.org.
signature.asc
Description: OpenPGP digital signature
