3D design optimization in a rotating frame?

ITA91 · May 19, 2015, 11:49

Hi there,

I'm trying to perform an optimization case in SU2 using a relatively simple 3D mesh in a rotating frame. I'm using a Caradonna & Tung rotor and I've based my cfg file on the cfg file provided in the tutorials for the optimization case of the NACA0012 (Euler). First I'd like to get it working for the Euler equations before moving on, in order to save time. When running shape_optimization.py the flow solution converges relatively quickly (as it should), but afterwards the adjoint solution is giving me trouble. I've tried several things to fix this, like using a coarser/finer mesh, adjusting the reduction factor of the CFL number of the adjoint problem (CFL_REDUCTION_ADJFLOW), using or not using multigrid, lowering the rotational speed (default is 2500 rpm, meaning that the used MACH_MOTION=0.877) but nothing seems to work. Every time the residuals of the adjoint problem increase (sometimes slowly, sometimes quickly, depending on the configuration) and after a while an increasing number of non-physical points appear. My main question is: before I waste too much time trying to fix my case, has anyone yet succeeded in running an optimization case for a rotating frame in 3D? For 2D it works, a normal flow analysis with a rotating frame in 3D works and a normal optimization case in 3D with a rigid frame works, but I haven't found a successful combination of them yet.

Thanks in advance for your help, I'm an inexperienced CFD user, so any help is welcome!

economon · May 19, 2015, 18:01

Hi,

Yes, we have had some success running design in these scenarios, and you might be interested in the following references:

http://adl.stanford.edu/papers/AIAA-2012-3018.pdf
http://adl.stanford.edu/papers/AIAA-2013-2580.pdf

However, the adjoint can be notoriously difficult to converge, and in our experience, these problems often had issues near the leading edges of the blades (when solved in the rotating frame). Many of the ideas that you tried above are good, and we regularly also try limiting the absolute value of the adjoint density variable to avoid divergence. You can do this with the following option:

%
% Limit value for the adjoint variable
LIMIT_ADJFLOW= 1E6

You might print out the adjoint solution and identify where in the mesh you are having issues, and then set the maximum value to something reasonable to keep this region from diverging.

Hope this helps,
Tom

ITA91 · May 20, 2015, 10:01

Thanks a lot for your quick and helpful reply, Tom! I'm definitely interested in the papers you referred to. I hope that I can find a solution soon, the largest residuals indeed occur at the leading edge (at the tip of the blade, to be exact).

Just out of curiosity: in the first paper you refer to you also use the Caradonna & Tung rotor and you've used a periodic BC there. Have you also used multigrid? So far I haven't been able to use multigrid in combination with a periodic BC, I always get an error when it tries to set the multigrid structure in that case.

May 19, 2015, 11:49	3D design optimization in a rotating frame?	#1
ITA91 New Member Join Date: Mar 2014 Location: Enschede, The Netherlands Posts: 8 Rep Power: 12	Hi there, I'm trying to perform an optimization case in SU2 using a relatively simple 3D mesh in a rotating frame. I'm using a Caradonna & Tung rotor and I've based my cfg file on the cfg file provided in the tutorials for the optimization case of the NACA0012 (Euler). First I'd like to get it working for the Euler equations before moving on, in order to save time. When running shape_optimization.py the flow solution converges relatively quickly (as it should), but afterwards the adjoint solution is giving me trouble. I've tried several things to fix this, like using a coarser/finer mesh, adjusting the reduction factor of the CFL number of the adjoint problem (CFL_REDUCTION_ADJFLOW), using or not using multigrid, lowering the rotational speed (default is 2500 rpm, meaning that the used MACH_MOTION=0.877) but nothing seems to work. Every time the residuals of the adjoint problem increase (sometimes slowly, sometimes quickly, depending on the configuration) and after a while an increasing number of non-physical points appear. My main question is: before I waste too much time trying to fix my case, has anyone yet succeeded in running an optimization case for a rotating frame in 3D? For 2D it works, a normal flow analysis with a rotating frame in 3D works and a normal optimization case in 3D with a rigid frame works, but I haven't found a successful combination of them yet. Thanks in advance for your help, I'm an inexperienced CFD user, so any help is welcome!

May 19, 2015, 18:01		#2
economon Super Moderator Thomas D. Economon Join Date: Jan 2013 Location: Stanford, CA Posts: 271 Rep Power: 14	Hi, Yes, we have had some success running design in these scenarios, and you might be interested in the following references: http://adl.stanford.edu/papers/AIAA-2012-3018.pdf http://adl.stanford.edu/papers/AIAA-2013-2580.pdf However, the adjoint can be notoriously difficult to converge, and in our experience, these problems often had issues near the leading edges of the blades (when solved in the rotating frame). Many of the ideas that you tried above are good, and we regularly also try limiting the absolute value of the adjoint density variable to avoid divergence. You can do this with the following option: % % Limit value for the adjoint variable LIMIT_ADJFLOW= 1E6 You might print out the adjoint solution and identify where in the mesh you are having issues, and then set the maximum value to something reasonable to keep this region from diverging. Hope this helps, Tom ITA91 likes this.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Convert DFT of pressure in rotating frame to stationary frame	CFD-123	FLUENT	0	February 12, 2015 17:51
Rotating rotor inside a frame	wllmk1	OpenFOAM	38	April 8, 2014 06:58
Different solutions for Rotating Reference Frame	fafou09	FLUENT	1	July 24, 2012 10:06
question about governing equation in CFX using rotating/non rotating reference frame	rystokes	CFX	0	January 12, 2010 06:14
Info: Short Course On Thermal Design of Electronic Equipment	Arnold Free	Main CFD Forum	0	August 10, 1999 10:18

May 20, 2015, 10:01		#3
ITA91 New Member Join Date: Mar 2014 Location: Enschede, The Netherlands Posts: 8 Rep Power: 12	Thanks a lot for your quick and helpful reply, Tom! I'm definitely interested in the papers you referred to. I hope that I can find a solution soon, the largest residuals indeed occur at the leading edge (at the tip of the blade, to be exact). Just out of curiosity: in the first paper you refer to you also use the Caradonna & Tung rotor and you've used a periodic BC there. Have you also used multigrid? So far I haven't been able to use multigrid in combination with a periodic BC, I always get an error when it tries to set the multigrid structure in that case.