[BramT]

Hi all,

My problem is as follows. After defining an FFD box, no matter what value I assign to DV_VALUE, the FFD_CONTROL_POINTS don't budge (checked using visual inspection in ParaView through ffd_boxes.dat). However, I'm using the exact same approach as for an older version of SU2, which worked. Below I will try to explain my process in creating the FFD box and attempting to move one of its control points. Hopefully someone can help me along in the right direction in getting this to work in SU2 v4.0.

First, I've defined the FFD box as follows:

FFD_TOLERANCE = 1E-10
FFD_ITERATIONS = 500
FFD_DEFINITION = (WING, 5.28885, 0.00000,-0.40000, 5.99403, 0.00000,-0.40000, 6.48809, 6.34030, 0.60000, 6.11228, 6.34030, 0.60000, 5.28885, 0.00000, 0.80000, 5.99403, 0.00000, 0.80000, 6.48809, 6.34030, 1.00000, 6.11228, 6.34030, 1.00000)
FFD_DEGREE = (2, 10, 1)
FFD_CONTINUITY = 2ND_DERIVATIVE
DV_KIND = FFD_SETTING
DV MARKER = ( WING_UPPER )

Then I run SU2_DEF for this .cfg-file to add the FFD box to the output mesh. This results in the following addition to my mesh output file (seems to be in order):

FFD_NBOX= 1
FFD_NLEVEL= 1
FFD_TAG= WING
FFD_LEVEL= 0
FFD_DEGREE_I= 2
FFD_DEGREE_J= 10
FFD_DEGREE_K= 1
FFD_PARENTS= 0
FFD_CHILDREN= 0
FFD_CORNER_POINTS= 8
5.28885	0	-0.4
5.99403	0	-0.4
6.48809	6.3403	0.6
6.11228	6.3403	0.6
5.28885	0	0.8
5.99403	0	0.8
6.48809	6.3403	1
6.11228	6.3403	1
FFD_CONTROL_POINTS= 66
0	0	0	5.28885	0	-0.4
0	0	1	5.28885	0	0.8
0	1	0	5.371193	0.63403	-0.3
0	1	1	5.371193	0.63403	0.82
0	2	0	5.453536	1.26806	-0.2
0	2	1	5.453536	1.26806	0.84
0	3	0	5.535879	1.90209	-0.1
0	3	1	5.535879	1.90209	0.86
0	4	0	5.618222	2.53612	0
0	4	1	5.618222	2.53612	0.88
0	5	0	5.700565	3.17015	0.1
0	5	1	5.700565	3.17015	0.9
0	6	0	5.782908	3.80418	0.2
0	6	1	5.782908	3.80418	0.92
0	7	0	5.865251	4.43821	0.3
0	7	1	5.865251	4.43821	0.94
0	8	0	5.947594	5.07224	0.4
0	8	1	5.947594	5.07224	0.96
0	9	0	6.029937	5.70627	0.5
0	9	1	6.029937	5.70627	0.98
0	10	0	6.11228	6.3403	0.6
0	10	1	6.11228	6.3403	1
1	0	0	5.64144	0	-0.4
1	0	1	5.64144	0	0.8
1	1	0	5.7073145	0.63403	-0.3
1	1	1	5.7073145	0.63403	0.82
1	2	0	5.773189	1.26806	-0.2
1	2	1	5.773189	1.26806	0.84
1	3	0	5.8390635	1.90209	-0.1
1	3	1	5.8390635	1.90209	0.86
1	4	0	5.904938	2.53612	0
1	4	1	5.904938	2.53612	0.88
1	5	0	5.9708125	3.17015	0.1
1	5	1	5.9708125	3.17015	0.9
1	6	0	6.036687	3.80418	0.2
1	6	1	6.036687	3.80418	0.92
1	7	0	6.1025615	4.43821	0.3
1	7	1	6.1025615	4.43821	0.94
1	8	0	6.168436	5.07224	0.4
1	8	1	6.168436	5.07224	0.96
1	9	0	6.2343105	5.70627	0.5
1	9	1	6.2343105	5.70627	0.98
1	10	0	6.300185	6.3403	0.6
1	10	1	6.300185	6.3403	1
2	0	0	5.99403	0	-0.4
2	0	1	5.99403	0	0.8
2	1	0	6.043436	0.63403	-0.3
2	1	1	6.043436	0.63403	0.82
2	2	0	6.092842	1.26806	-0.2
2	2	1	6.092842	1.26806	0.84
2	3	0	6.142248	1.90209	-0.1
2	3	1	6.142248	1.90209	0.86
2	4	0	6.191654	2.53612	0
2	4	1	6.191654	2.53612	0.88
2	5	0	6.24106	3.17015	0.1
2	5	1	6.24106	3.17015	0.9
2	6	0	6.290466	3.80418	0.2
2	6	1	6.290466	3.80418	0.92
2	7	0	6.339872	4.43821	0.3
2	7	1	6.339872	4.43821	0.94
2	8	0	6.389278	5.07224	0.4
2	8	1	6.389278	5.07224	0.96
2	9	0	6.438684	5.70627	0.5
2	9	1	6.438684	5.70627	0.98
2	10	0	6.48809	6.3403	0.6
2	10	1	6.48809	6.3403	1
FFD_SURFACE_POINTS= 123
WING_UPPER	127853	8.226298993063333e-01	8.985001057312524e-01	6.493972666542595e-01
WING_UPPER	27364	9.690594492835374e-01	9.384025615966917e-01	6.468360537974274e-01
WING_UPPER	25792	6.196250163451079e-01	9.254079532276076e-01	6.755166121471692e-01
WING_UPPER	27297	9.637467276860744e-01	5.975994043142797e-01	5.820256224404885e-01
WING_UPPER	27857	2.063432598605249e-01	5.789738905519113e-01	6.486300793194301e-01
WING_UPPER	127835	2.393827774123823e-01	5.557402534255793e-01	6.436903783756670e-01
WING_UPPER	27300	5.941398387919182e-01	5.727936514089166e-01	6.150234253695677e-01
WING_UPPER	26993	5.124724642910951e-01	2.443376983386993e-01	5.966089009849213e-01
WING_UPPER	26991	7.337216689852011e-01	2.658247953333422e-01	5.767706414147422e-01
WING_UPPER	26990	4.854531403688018e-01	2.719351857164670e-01	6.007841746332410e-01
WING_UPPER	27947	9.662934438295363e-01	2.439137921130930e-01	5.515967224269148e-01
WING_UPPER	127818	7.640984376281444e-01	2.407883200221859e-01	5.722341902480614e-01
WING_UPPER	26992	9.278967790827599e-01	2.173479804085976e-01	5.539687394971078e-01
WING_UPPER	28276	8.316318431601394e-01	2.942185404699176e-01	5.687176784824913e-01
WING_UPPER	127814	1.309782264169977e-02	4.924015086135676e-01	6.553120581815435e-01
WING_UPPER	27233	2.726044697457866e-01	4.776612427355620e-01	6.340054724802396e-01
WING_UPPER	27915	2.732683252839532e-01	5.023503114381348e-01	6.360932065328483e-01
WING_UPPER	24799	2.490321531348474e-02	4.716076902757700e-01	6.527812490528665e-01
WING_UPPER	26223	8.477980011341163e-01	1.363152764211193e-01	5.578108753785620e-01
WING_UPPER	27945	3.349004626080829e-01	1.783337731376947e-01	6.089690792287388e-01
WING_UPPER	27017	3.672108208701332e-01	2.094287217108024e-01	6.078753172361538e-01
WING_UPPER	27023	5.505424728200542e-02	2.018028093663182e-01	6.337460570093585e-01
WING_UPPER	127608	7.900361986796682e-01	4.395365609126674e-01	5.837296651470666e-01
WING_UPPER	24750	9.577100289328770e-01	4.017846606095486e-01	5.632972986859264e-01
WING_UPPER	27248	8.851102630774988e-01	4.638326041523937e-01	5.763294315288588e-01
WING_UPPER	27250	5.181023883731359e-01	4.568903998820189e-01	6.106886502905092e-01
WING_UPPER	126590	4.072521496339458e-01	9.979560741363684e-01	7.187434813777137e-01
WING_UPPER	24990	5.777660456544114e-01	9.799450708696307e-01	6.976898342430636e-01
WING_UPPER	1965	8.390620766585954e-01	9.980205872721121e-01	6.806554121708749e-01
WING_UPPER	24797	8.976212290792989e-01	3.223647218143263e-01	5.638887291016695e-01
WING_UPPER	27103	6.138530012081629e-01	3.180500009569069e-01	5.916958899640808e-01
WING_UPPER	26999	5.630680549048220e-01	2.920652788136938e-01	5.947915530724819e-01
WING_UPPER	28233	8.869947620912314e-01	9.731783648104146e-01	6.666588276553664e-01
WING_UPPER	24988	5.552557381903700e-01	9.565470115365695e-01	6.911589939907652e-01
WING_UPPER	25128	7.555847094250663e-01	3.686335250535980e-02	5.625585687301246e-01
WING_UPPER	25137	5.984091814451001e-01	1.878737646938907e-02	5.773629400490495e-01
WING_UPPER	1919	8.321677620451017e-01	0.000000000000000e+00	5.532086590925852e-01
WING_UPPER	26250	9.554314192618468e-01	5.953685475071115e-02	5.428052210155739e-01
WING_UPPER	28114	3.308324887136690e-02	9.026592910364326e-01	7.126876503137387e-01
WING_UPPER	27361	1.973070335934081e-01	8.800317309939586e-01	6.954266341056098e-01
WING_UPPER	27368	3.846542509983683e-01	9.015050838454756e-01	6.876072740160947e-01
WING_UPPER	25791	9.950614083877612e-02	9.221146847884096e-01	7.141884672235105e-01
WING_UPPER	26976	6.782819892236971e-01	2.148851277295009e-01	5.792281477816601e-01
WING_UPPER	27946	3.257679262458614e-01	2.896008467984497e-01	6.160146649815754e-01
WING_UPPER	26998	1.464975046961341e-01	3.131330948666556e-01	6.324162378439361e-01
WING_UPPER	26995	7.258802900708010e-02	2.844677947927773e-01	6.368316486641689e-01
WING_UPPER	26997	3.955924117116066e-01	3.111584506974868e-01	6.112857911990101e-01
WING_UPPER	26996	2.893260530774340e-01	2.648852111020792e-01	6.177427953644981e-01
WING_UPPER	26245	6.146780148121669e-01	1.586553187578022e-01	5.821850910685294e-01
WING_UPPER	27018	6.315768542683767e-01	1.861337528445658e-01	5.820254119555408e-01
WING_UPPER	27774	3.902224466115193e-01	1.544678770034665e-01	6.029387751054692e-01
WING_UPPER	27775	1.503888567491181e-01	1.526995829095619e-01	6.235081122746838e-01
WING_UPPER	27772	4.855862785168744e-02	1.279842010659773e-01	6.307029740276898e-01
WING_UPPER	27941	2.630966613838909e-01	2.408300651604563e-01	6.186106958241288e-01
WING_UPPER	27002	6.682528426276886e-02	2.561783810518326e-01	6.357001178683028e-01
WING_UPPER	27025	7.778206271688509e-02	2.293533296135193e-01	6.333716696395461e-01
WING_UPPER	27917	9.306792295467401e-01	6.788942623335471e-01	5.958179520246558e-01
WING_UPPER	24761	9.242056033297922e-01	7.048589980987890e-01	6.002544569992467e-01
WING_UPPER	27212	6.110842574816362e-01	6.898685939705294e-01	6.280001833480566e-01
WING_UPPER	27189	7.028880518997911e-01	6.678089998427332e-01	6.165658982380363e-01
WING_UPPER	27293	9.463132734093330e-01	6.561620136583054e-01	5.911008089200359e-01
WING_UPPER	27231	5.612718159883857e-01	5.109695249445669e-01	6.118891415901143e-01
WING_UPPER	27232	2.621565049354669e-01	5.288827498686750e-01	6.393519779732723e-01
WING_UPPER	27234	5.564431300563668e-01	4.846319111337790e-01	6.098488315366124e-01
WING_UPPER	27912	9.511249780893161e-01	3.514872160569729e-01	5.601759708840301e-01
WING_UPPER	27607	6.707336783246510e-01	3.469944895195508e-01	5.881906018208770e-01
WING_UPPER	27610	7.773023735794446e-01	3.753491268359815e-01	5.799094121851485e-01
WING_UPPER	27904	9.600424249716272e-01	6.307219619670356e-01	5.864442301114613e-01
WING_UPPER	24768	6.664510679894105e-01	6.474229948418759e-01	6.171868765335121e-01
WING_UPPER	27294	6.490009469353895e-01	6.228366969912262e-01	6.158516344875473e-01
WING_UPPER	27299	9.205964115204106e-01	5.449743776116998e-01	5.803296998371589e-01
WING_UPPER	27360	2.188202077825722e-01	8.167557739605544e-01	6.805177712289505e-01
WING_UPPER	24733	4.321971989075017e-01	8.538563219018627e-01	6.719678816612007e-01
WING_UPPER	27355	6.347391223566178e-01	8.779592420630560e-01	6.608608735058676e-01
WING_UPPER	24754	5.763589189385444e-01	5.385362148518055e-01	6.132907332687056e-01
WING_UPPER	27629	2.903041962468326e-01	5.975552575892125e-01	6.439490462162361e-01
WING_UPPER	27773	3.016252738032389e-01	1.286265658401208e-01	6.094753144872452e-01
WING_UPPER	27771	1.905968469574787e-01	1.031903309218393e-01	6.179196592444833e-01
WING_UPPER	27608	2.239123533261586e-01	4.068027957990684e-01	6.321683201167207e-01
WING_UPPER	27611	1.207617552025279e-01	3.730155179752936e-01	6.382707734776633e-01
WING_UPPER	27609	4.546425313921721e-01	3.746359355404213e-01	6.100650212616245e-01
WING_UPPER	27606	3.378244478630481e-01	3.401862314385984e-01	6.180551660192708e-01
WING_UPPER	27348	6.816861671730818e-01	7.936928176126240e-01	6.387186017365248e-01
WING_UPPER	24737	7.326968453689370e-01	8.307170512474930e-01	6.412894004894611e-01
WING_UPPER	27247	1.472419453226504e-01	4.349346689306059e-01	6.405289926617233e-01
WING_UPPER	27241	2.624998219882771e-01	4.534184938097648e-01	6.327475431616625e-01
WING_UPPER	27230	8.114925205545813e-01	4.934304941622386e-01	5.861461405119541e-01
WING_UPPER	27227	5.614376394540139e-01	6.019498486476335e-01	6.214075841840121e-01
WING_UPPER	27630	3.445325112602489e-01	6.199630241658172e-01	6.421208516288941e-01
WING_UPPER	27211	5.288530589213575e-01	7.608072793455378e-01	6.458488672928481e-01
WING_UPPER	27206	6.330219444043915e-02	7.373744444158379e-01	6.788797103732457e-01
WING_UPPER	27207	7.222894566626384e-02	7.672463524103362e-01	6.828470783009838e-01
WING_UPPER	27213	5.451491502135047e-01	7.214304863228650e-01	6.382361456098045e-01
WING_UPPER	27200	2.171462590541912e-01	7.888840279717637e-01	6.754865474521289e-01
WING_UPPER	27628	7.590960995077814e-02	6.421128295226506e-01	6.656795783013454e-01
WING_UPPER	27193	2.002781795809631e-01	7.000269258917088e-01	6.634083771869337e-01
WING_UPPER	27188	3.662874686827030e-01	6.660772754019416e-01	6.459032488486635e-01
WING_UPPER	27187	4.131100465070147e-01	6.396569335147092e-01	6.386998891350368e-01
WING_UPPER	27032	8.724693790757823e-01	1.896641174442033e-01	5.582555014097054e-01
WING_UPPER	26251	4.647291876873049e-01	1.026870715440646e-01	5.936357244710582e-01
WING_UPPER	26249	6.549614385941567e-01	7.257955772948976e-02	5.740989428347034e-01
WING_UPPER	26253	8.958921697363189e-01	8.511311635289240e-02	5.504070292180907e-01
WING_UPPER	26244	8.050398627898158e-01	1.663452702444919e-01	5.638049605993616e-01
WING_UPPER	26241	5.610489473168904e-01	1.309078115925124e-01	5.860010182790396e-01
WING_UPPER	26224	7.415359193443715e-01	1.054547960133274e-01	5.670357984911217e-01
WING_UPPER	25790	2.666223113838946e-01	9.419940870254324e-01	7.089317334776143e-01
WING_UPPER	25788	2.281858832999392e-01	9.745414980205830e-01	7.230786236140189e-01
WING_UPPER	1922	4.156055602408112e-02	5.528408851449792e-11	6.261210093049641e-01
WING_UPPER	25136	1.501346375569656e-01	2.483934708830480e-02	6.180442820436578e-01
WING_UPPER	1921	2.199434178533936e-01	8.575072025401097e-23	6.111627022425333e-01
WING_UPPER	25135	3.951405578642770e-01	1.435861889970770e-02	5.962280029772223e-01
WING_UPPER	1920	3.928910474545399e-01	2.949550924621363e-27	5.958608488241831e-01
WING_UPPER	25129	4.322955072317019e-01	4.276054401065890e-02	5.939505746266024e-01
WING_UPPER	25126	3.302138126523320e-01	7.625969306525000e-02	6.045558459372192e-01
WING_UPPER	25123	9.059439112512011e-02	8.162886325507031e-02	6.253391734173848e-01
WING_UPPER	25121	1.339118189073332e-01	5.626638794919999e-02	6.206277163106709e-01
WING_UPPER	1917	5.742146688289456e-01	1.311481254929932e-29	5.789277652899424e-01
WING_UPPER	24760	9.569089183025173e-01	7.362615260641979e-01	6.019092484329639e-01
WING_UPPER	24755	8.785217979817210e-01	5.169201936582367e-01	5.820101178500799e-01
WING_UPPER	24748	6.018222441378711e-01	4.047787403935025e-01	5.990179120462047e-01
WING_UPPER	24744	4.500272686876063e-01	4.335562633992351e-01	6.149471819851093e-01
WING_UPPER	24736	9.075697231695496e-01	8.561036854612449e-01	6.304782359516546e-01
WING_UPPER	24699	7.311464050736993e-02	6.157214800574932e-01	6.629383264067141e-01

Once the box is added, I change/add the following to the .cfg-file:

DV_KIND = FFD_CONTROL_POINT; FFD_CONTROL_POINT % ALSO TRIED WITH FFD_SETTING; FFD SETTING %
DV_MARKER = ( WING_UPPER )
DV_PARAM = ( WING, 1, 1, 1, 0, 0, 1 ); ( WING, 1, 2, 1, 0, 0, 1)
DV_VALUE = 0.1, 0.1

Again, I run SU2_DEF with this updated .cfg-file as input. Once it is done I compare the old ffd_boxes.dat with the new ffd_boxes.dat, in paraview, and they overlap exactly.. no changes in any of the control point locations.

If anybody knows what I'm doing wrong, please let me know! Also, if additional information is required to help solve this issue, feel free to let me know so I can add it to this post.

Regards,
BramT

[BramT]

Maybe the following is off help in identifying the issue to someone with better understanding of the code than me:

After spending too much time on the above issue, I figured I might find a clue by trying to adjust the FFD box from one of the TestCases (TestCases/optimization_euler/steady_oneram6 to be specific).

Initially, I messed around with the top layer of control points of the FFD box and I started noticing something weird: it seemed like the inboard 3 rows of control points, in spanwise direction, were ignoring any kind of input, while the other control points all were behaving as expected. In order to confirm my hunch I translated ALL the upper surface control points upwards by 0.2 and inspected the results in paraview, of which the results can be seen in the (linked) image below.

[demanosalvas]

BramT,

In version 4.0 we added the functionality to be able to maintain 1st and 2nd derivative continuity of the surface when being deformed. This translates into the code not allowing you to move the first plane of the FFD after intercepting the surface for 1st derivative continuity, and the first two planes after intercepting the surface for 2nd derivative continuity. In your case you have the FFD_CONTINUITY set to 2ND_DERIVATIVE and this might be the reason why you can't move some of the point.

Hope this helps,

David

[BramT]

Originally Posted by demanosalvas

BramT,

In version 4.0 we added the functionality to be able to maintain 1st and 2nd derivative continuity of the surface when being deformed. This translates into the code not allowing you to move the first plane of the FFD after intercepting the surface for 1st derivative continuity, and the first two planes after intercepting the surface for 2nd derivative continuity. In your case you have the FFD_CONTINUITY set to 2ND_DERIVATIVE and this might be the reason why you can't move some of the point.

Hope this helps,

David

Thanks for your reply David!

I'll upgrade to SU2 v4.0 soon and verify if this is indeed the case.

Regards,
BramT

[Xiaosong]

Hi,

I am working on the same test with you. What I try to do is to change the sweep angle of wing, using SU2_DEF, to a specific value. However, it seems no FFD method works for this. Could you please offer me some advices? Or, do you know how to do translation for FFD box, please?

Thanks!

XIaosong