Closed ChristianBauerEng closed 1 year ago
What is the output of SU2_DEF? That usually means the linear solver did not converge.
Good point! Indeed the solver has a hard time converging and the resudual barely drops below 4E-8 (far away from the limit specified in the tutorial).
--------------------- Surface grid deformation (ZONE 0) -----------------
Performing the deformation of the surface grid.
1 Free Form Deformation boxes.
1 Free Form Deformation nested levels.
FFD box tag: BOX. FFD box level: 0. Degrees: 4, 1.
FFD Blending using Bezier Curves.
Number of parent boxes: 0. Number of child boxes: 0.
Corner points: 4. Control points: 20. Surface points: 3404.
Writing a Paraview file of the FFD boxes.
----------------- FFD technique (parametric -> cartesian) ---------------
Checking FFD box dimension.
Checking FFD box intersections with the solid surfaces.
SU2 is fixing user's input planes.
Update cartesian coord | FFD box: BOX. Max Diff: 0.00246933.
Writing a Paraview file of the FFD boxes.
------------------- Volumetric grid deformation (ZONE 0) ----------------
Performing the deformation of the volumetric grid.
Computing volumes of the grid elements.
Min. volume: 1.15412e-13, max. volume: 7.81947e-08.
Min. distance: 5e-07, max. distance: 0.0061806.
# FGMRES residual history
# Residual tolerance target = 1e-14
# Initial residual norm = 9.32748e+06
0 1
10 4.70151e-07
20 1.0396e-07
30 9.99835e-08
40 9.99358e-08
50 9.94837e-08
60 9.79571e-08
70 9.48291e-08
80 9.26963e-08
90 9.25628e-08
100 9.2495e-08
110 9.13192e-08
120 8.82461e-08
130 8.5977e-08
140 8.55573e-08
150 8.55348e-08
160 8.52077e-08
170 8.32895e-08
180 8.13964e-08
190 8.10219e-08
...
Any suggestions on which settings may help the solver to converge?
You can try CONJUGATE_GRADIENT instead of FGMRES, or running on a single core, or CONSTANT_STIFFNESS, or setting a large value (e.g. 1E6) for DEFORM_COEFF
Using CG did the trick it seems, and the solver now converges < 1e-18. Thanks a lot!
However, one small issue remains: near the tip, the deformation looks quite strange:
White box is undeformed FFD-box, translucent green is the deformed one.
BC on the bottom is axis, right and top boundaries are walls. Top and right wall are included in DV_MARKER and DV-direction is (0., 1., 0.).
Any idea what could be behind this behaviour?
The ffd doesn't cover the entire right surface, so some points on it are fixed while other are dragged to a far location. The deformation you are trying to apply is discontinuous. I would not move the bottom points of the ffd.
Describe the bug I'd like to use the FFD shape optimization feature of SU2. From a geometrical point of view the domain is very similar to the "2 way mixing channel optimization" tutorial example. I've therefore followed the steps of this tutorial and tested the mesh deformation with some dummy values.
The result of the deformation can be seen here:![grafik](https://user-images.githubusercontent.com/42602495/211366671-1344820b-87f6-429f-b1d4-c396a6467440.png)
The lower halve of the picture shows the original, undeformed mesh. The upper halve contains in green the original FFD box (ffd_box_0.vtk), in white the deformed FFD box and the deformed mesh (volume_deformed_0000.vtu) in blue/white.
From the definition of the deformation I'd expect that the cells outside the FFD box to remain largely unchanged. Instead, here the largest deformations occur.
Following the mixer-example the FFD information was added to the mesh by running SU2_DEF with a config file containing this information:
This added the following info to the mesh:
The deformation was then tested by running SU2_DEF with the two leftmost FFD control points being fixed:
Desktop (please complete the following information):