CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Programming & Development (
-   -   Adjoint sensitivity on specified surface (

robyTKD April 29, 2013 05:05

Adjoint sensitivity on specified surface
Hi Foamers,

I would like to compute shape sensitivity of a specified surface after every time cycle of adjointShapeOptimizationFoam.
The formula is:

G = \frac{\vec{d} \cdot \nabla p}{C_{\infty}} - \nu \left[ \left( \vec{n} \cdot \nabla \right) \left( \vec{n} \cdot \nabla \right) \vec{v} \right] \cdot \vec{\Psi}^u


\vec{d} = (\cos (\alpha), \sin (\alpha))

C_{\infty} = \frac{1}{2} \rho U_{\infty}^2

\vec{v} primal vector velocity

p pressure

\vec{\Psi}^u vector adjoint velocity

Is it possible to create a post-processing script like forceCoeffs.C in order to store the sensitivity of every cell on the specified surface? Is it also possible to store its integral on the surface?

Another idea is to compute the sensitivity at each time step in the solver (adjointShapeOptimizationFoam.C), is it a simpler way to proceed?

Thank you,

All times are GMT -4. The time now is 04:58.