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

Implementation of continuous adjoint SA model

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 23, 2018, 05:34
Default Implementation of continuous adjoint SA model
  #1
New Member
 
c3k1uI5gzT
Join Date: Jul 2018
Posts: 2
Rep Power: 0
c3k1uI5gzT is on a distinguished road
Hi,

I am trying to implement a continuous adjoint formulation of the Spalart-Allmaras turbulence model. This comes down to programming three additional (adjoint) equations for the adjoint velocity, adjoint pressure, and the adjoint SA-variable. The first two are similar to the velocity and pressure equation with some minor modifications and work fine. The third is similar to the SA model transport equation.

My question is: What would be the most convenient way to implement this third equation, provided that
- it has a form similar to the SA model transport equation
- it uses both the velocity and the adjoint velocity
- it uses the SA coefficients/quantities

Proposed solution 1: implement it in a way similar to UEqn, pEqn, (as I did for UaEqn, paEqn for the other adjoint variables)
- difficulty: how to access the SA-quantities (e.g. fv1, fv2, coefficients, etc.) from the solver

Proposed solution 2: implement it in a way similar to the SA turbulence model, i.e. as a turbulence model
- difficulty: would require the use of two turbulence models/adjoint equation is not strictly a turbulence model

Thanks!
c3k1uI5gzT is offline   Reply With Quote

Old   August 13, 2018, 05:16
Default
  #2
New Member
 
c3k1uI5gzT
Join Date: Jul 2018
Posts: 2
Rep Power: 0
c3k1uI5gzT is on a distinguished road
For anyone wondering, I went for solution 2. I implemented the adjoint nuTilda equation in the same file as the nuTilda equation (spalartAllmaras.C). For the variables I needed which were not available in the turbulence model I used for example

Code:
this->db().objectRegistry::lookupObject<volVectorField> ("Ua");
c3k1uI5gzT is offline   Reply With Quote

Reply

Tags
additional equation, adjoint

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Issues in FGM combustion model implementation Lisandro Maders OpenFOAM Programming & Development 15 April 22, 2020 20:18
Turbulence model implementation tilasoldo OpenFOAM Programming & Development 2 July 29, 2018 16:09
BL k-ε-v2/k model implementation - OpenFOAM 4 tomaszkura OpenFOAM Programming & Development 1 April 27, 2018 04:20
Status of continuous RANS adjoint praveen SU2 1 April 11, 2013 14:12
UDF for linear PTT model implementation dpitz Fluent UDF and Scheme Programming 1 February 4, 2011 11:08


All times are GMT -4. The time now is 15:48.