CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

Absence of a term in adjointShapeOptimizationFoam outlet boundary conditions

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

Like Tree3Likes
  • 3 Post By robyTKD

Reply
 
LinkBack Thread Tools Display Modes
Old   April 8, 2013, 08:02
Default Absence of a term in adjointShapeOptimizationFoam outlet boundary conditions
  #1
Member
 
Roberto Pieri
Join Date: Feb 2012
Location: Milan
Posts: 57
Rep Power: 5
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,
Roberto
samiam1000, Dome and Garry974 like this.
robyTKD is offline   Reply With Quote

Old   April 10, 2013, 06:05
Default
  #2
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:

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

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

Let's hope for some help

Domenico
Dome is offline   Reply With Quote

Reply

Tags
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
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
Domain Imbalance HMR CFX 3 March 6, 2011 21:10
Help with boundary conditions Dan CFX 0 April 3, 2006 11:32


All times are GMT -4. The time now is 19:51.