Optimization with the Source term in NS equation

milind543 · November 3, 2014, 00:57

I am interested in shape optimization of a 3D wing shape in Incompressible flow conditions with a source term in the NS equation. can I get help in doing so? However I already have my own subroutine which calculates the flux and I want to put as a source term in the governing equation and finally use it in the optimization procedure this will also include to modify the adjoint equation. Please let me know if this is possible.

Thanks

economon · January 5, 2015, 14:53

Hi,

Sorry for a slow reply. Indeed, it is straightforward to include source terms in the flow and adjoint problems.

You can add an additional source term in the direct problem in the routine CEulerSolver::Source_Residual(...) in the file solver_direct_mean.cpp. You might follow the rotating frame source term at the top of the routine as an example.

Similarly, the source term can be included for the adjoint equations: after the proper derivation, the required adjoint source term should be the transpose of the Jacobian of the direct source term with respect to the conservative variables multiplied by the adjoint state. This can then be included in the CAdjNSSolver::Source_Residual(...) routine in the solver_adjoint_mean.cpp file.

Hope this helps,

Thomas D. Economon
SU2 lead developer

November 3, 2014, 00:57	Optimization with the Source term in NS equation	#1
milind543 New Member milind Join Date: Sep 2011 Posts: 6 Rep Power: 14	I am interested in shape optimization of a 3D wing shape in Incompressible flow conditions with a source term in the NS equation. can I get help in doing so? However I already have my own subroutine which calculates the flux and I want to put as a source term in the governing equation and finally use it in the optimization procedure this will also include to modify the adjoint equation. Please let me know if this is possible. Thanks

Thread Tools	Search this Thread
Show Printable Version Email this Page	Search this Thread: Advanced Search
Display Modes
Linear Mode Switch to Hybrid Mode Switch to Threaded Mode

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
SparceImage v1.7.x Issue on MAC OS X	rcarmi	OpenFOAM Installation	4	August 14, 2014 06:42
centOS 5.6 : paraFoam not working	yossi	OpenFOAM Installation	2	October 9, 2013 01:41
add source term to transport equation	m003020002	Fluent UDF and Scheme Programming	0	September 10, 2013 03:09
friction forces icoFoam	ofslcm	OpenFOAM	3	April 7, 2012 10:57
"parabolicVelocity" in OpenFoam 2.1.0 ?	sawyer86	OpenFOAM Running, Solving & CFD	21	February 7, 2012 11:44

January 5, 2015, 14:53		#2
economon Super Moderator Thomas D. Economon Join Date: Jan 2013 Location: Stanford, CA Posts: 271 Rep Power: 14	Hi, Sorry for a slow reply. Indeed, it is straightforward to include source terms in the flow and adjoint problems. You can add an additional source term in the direct problem in the routine CEulerSolver::Source_Residual(...) in the file solver_direct_mean.cpp. You might follow the rotating frame source term at the top of the routine as an example. Similarly, the source term can be included for the adjoint equations: after the proper derivation, the required adjoint source term should be the transpose of the Jacobian of the direct source term with respect to the conservative variables multiplied by the adjoint state. This can then be included in the CAdjNSSolver::Source_Residual(...) routine in the solver_adjoint_mean.cpp file. Hope this helps, Thomas D. Economon SU2 lead developer