CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

Absence of a term in adjointShapeOptimizationFoam outlet boundary conditions

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree4Likes
  • 3 Post By robyTKD
  • 1 Post By fumiya

LinkBack Thread Tools Display Modes
Old   April 8, 2013, 08:02
Default Absence of a term in adjointShapeOptimizationFoam outlet boundary conditions
Roberto Pieri
Join Date: Feb 2012
Location: Milan
Posts: 57
Rep Power: 8
robyTKD is on a distinguished road
Hi Foamers,

I have a question about boundary conditions of adjointShapeOptimizationFoam. In particular, following Othmer's procedure in A continuous adjoint formulation for the computation of topological and surface sensitivities of ducted flows, boundary conditions for the Outlet are:

q = \textbf{u} \cdot \textbf{v} + u_n v_n + \nu (\textbf{n} \cdot \nabla) u_n + \frac{\partial J_{\Gamma}}{v_n}

0 = v_n \textbf{u}_t + \nu (\textbf{n} \cdot \nabla) \textbf{u}_t + \frac{\partial J_{\Gamma}}{\textbf{v}_t}

but I can't find \nu (\textbf{n} \cdot \nabla) \textbf{u}_t and \nu (\textbf{n} \cdot \nabla) u_n in adjointOutletPressureFvPatchScalarField.C and adjointOutletVelocityFvPatchScalarField.C.

Are they neglected because they are irrelevant or is simply the normal gradient of adjoint velocity forced to be null?

Best regards,
samiam1000, Dome and Garry974 like this.
robyTKD is offline   Reply With Quote

Old   April 10, 2013, 06:05
New Member
Domenico Sgrò
Join Date: Dec 2012
Location: Milan, IT
Posts: 2
Rep Power: 0
Dome is on a distinguished road
Hi Roberto!
I have exactly the same doubt, since in the code:

operator==((phiap/patch().magSf() - 1.0)*phip/patch().magSf() + (Up & Uap));

vectorField::operator=(phiap*patch().Sf()/sqr(patch().magSf()) + UtHat);
there isn't any "normal gradient" term.

Let's hope for some help

Dome is offline   Reply With Quote

Old   January 3, 2016, 01:13
Senior Member
Fumiya Nozaki
Join Date: Jun 2010
Location: Yokohama, Japan
Posts: 217
Blog Entries: 1
Rep Power: 11
fumiya is on a distinguished road

You can find how to implement these terms in the following document:

Description of adjointShapeOptimizationFoam and how to implement new objective functions

Best regards,
sylvester likes this.
fumiya is online now   Reply With Quote


adjointshapeoptimization, boundaries condition

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Domain Imbalance HMR CFX 5 October 10, 2016 05:57
Radiation interface hinca CFX 15 January 26, 2014 18:11
ribbed channel / simpleFoam / boundary conditions beeo OpenFOAM Pre-Processing 20 July 17, 2013 08:39
domain imbalance for enrgy equation happy CFX 14 September 6, 2012 01:54
Help with boundary conditions Dan CFX 0 April 3, 2006 11:32

All times are GMT -4. The time now is 10:57.