|
[Sponsors] |
May 3, 2017, 22:20 |
Mesh not modified while shape optimization
|
#1 |
New Member
LiuYang
Join Date: Mar 2015
Posts: 4
Rep Power: 11 |
I tried to use SU2 to modify a part of a solid surface, for vehicle aerodynamic optimization. The FFD box is constructed around a part of the surface. But the iteration stops after one step and the mesh seems to keep undeformed. SU2_CFD runs OK, so I think the settings in the config file seems OK. Could anyone please tell me what is the source of this problem? Thank you in advance.
The following is my parameters. % -------------------- FREE-FORM DEFORMATION PARAMETERS -----------------------% % % Tolerance of the Free-Form Deformation point inversion FFD_TOLERANCE= 1E-8 % % Maximum number of iterations in the Free-Form Deformation point inversion FFD_ITERATIONS= 500 % % FFD box definition: 3D case (FFD_BoxTag, X1, Y1, Z1, X2, Y2, Z2, X3, Y3, Z3, X4, Y4, Z4, % X5, Y5, Z5, X6, Y6, Z6, X7, Y7, Z7, X8, Y8, Z8) % 2D case (FFD_BoxTag, X1, Y1, 0.0, X2, Y2, 0.0, X3, Y3, 0.0, X4, Y4, 0.0, % 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0) FFD_DEFINITION= (EngineCoverTip, -0.84, -0.36, 0.59, -0.8, -0.36, 0.59, -0.8, 0.36, 0.59, -0.84, 0.36, 0.59, -0.84, -0.36, 0.62, -0.8, -0.36, 0.62, -0.8, 0.36, 0.62, -0.84, 0.36, 0.62) % % FFD box degree: 3D case (x_degree, y_degree, z_degree) % 2D case (x_degree, y_degree, 0) FFD_DEGREE= (6, 8, 6) % % Surface continuity at the intersection with the FFD (1ST_DERIVATIVE, 2ND_DERIVATIVE) FFD_CONTINUITY= 1ST_DERIVATIVE % ----------------------- DESIGN VARIABLE PARAMETERS --------------------------% % % Kind of deformation (FFD_SETTING, FFD_CONTROL_POINT_2D, FFD_CAMBER_2D, FFD_THICKNESS_2D, % HICKS_HENNE, COSINE_BUMP, PARABOLIC, % NACA_4DIGITS, DISPLACEMENT, ROTATION, FFD_CONTROL_POINT, % FFD_DIHEDRAL_ANGLE, FFD_TWIST_ANGLE, FFD_ROTATION, % FFD_CAMBER, FFD_THICKNESS, FFD_CONTROL_SURFACE, SURFACE_FILE, AIRFOIL) DV_KIND= FFD_CONTROL_POINT % % Marker of the surface in which we are going apply the shape deformation DV_MARKER= ( car_ground ) % % Parameters of the shape deformation % - FFD_CONTROL_POINT ( FFD_BoxTag, i_Ind, j_Ind, k_Ind, x_Disp, y_Disp, z_Disp ) % - FFD_DIHEDRAL_ANGLE ( FFD_BoxTag, x_Orig, y_Orig, z_Orig, x_End, y_End, z_End ) % - FFD_TWIST_ANGLE ( FFD_BoxTag, x_Orig, y_Orig, z_Orig, x_End, y_End, z_End ) % - FFD_ROTATION ( FFD_BoxTag, x_Orig, y_Orig, z_Orig, x_End, y_End, z_End ) % - FFD_CAMBER ( FFD_BoxTag, i_Ind, j_Ind ) % - FFD_THICKNESS ( FFD_BoxTag, i_Ind, j_Ind ) % - FFD_VOLUME ( FFD_BoxTag, i_Ind, j_Ind ) DV_PARAM= ( EngineCoverTip, 1, 0, 0, 1.0, 0.0, 0.0 ) % % New value of the shape deformation DV_VALUE= 0.0 % ------------------------ GRID DEFORMATION PARAMETERS ------------------------% % % Linear solver or smoother for implicit formulations (FGMRES, RESTARTED_FGMRES, BCGSTAB) DEFORM_LINEAR_SOLVER= FGMRES % % Number of smoothing iterations for mesh deformation DEFORM_LINEAR_ITER= 800 % % Number of nonlinear deformation iterations (surface deformation increments) DEFORM_NONLINEAR_ITER= 2 % % Print the residuals during mesh deformation to the console (YES, NO) DEFORM_CONSOLE_OUTPUT= NO % % Factor to multiply smallest cell volume for deform tolerance (0.001 default) DEFORM_TOL_FACTOR = 1E-6 % % Type of element stiffness imposed for FEA mesh deformation (INVERSE_VOLUME, % WALL_DISTANCE, CONSTANT_STIFFNESS) DEFORM_STIFFNESS_TYPE= WALL_DISTANCE % % Visualize the deformation (NO, YES) VISUALIZE_DEFORMATION= NO % --------------------- OPTIMAL SHAPE DESIGN DEFINITION -----------------------% % Available flow based objective functions or constraint functions % DRAG, LIFT, SIDEFORCE, EFFICIENCY, % FORCE_X, FORCE_Y, FORCE_Z, % MOMENT_X, MOMENT_Y, MOMENT_Z, % THRUST, TORQUE, FIGURE_OF_MERIT, % EQUIVALENT_AREA, NEARFIELD_PRESSURE, % FREE_SURFACE % % Available geometrical based objective functions or constraint functions % MAX_THICKNESS, 1/4_THICKNESS, 1/2_THICKNESS, 3/4_THICKNESS, AREA, AOA, CHORD, % MAX_THICKNESS_SEC1, MAX_THICKNESS_SEC2, MAX_THICKNESS_SEC3, MAX_THICKNESS_SEC4, MAX_THICKNESS_SEC5, % 1/4_THICKNESS_SEC1, 1/4_THICKNESS_SEC2, 1/4_THICKNESS_SEC3, 1/4_THICKNESS_SEC4, 1/4_THICKNESS_SEC5, % 1/2_THICKNESS_SEC1, 1/2_THICKNESS_SEC2, 1/2_THICKNESS_SEC3, 1/2_THICKNESS_SEC4, 1/2_THICKNESS_SEC5, % 3/4_THICKNESS_SEC1, 3/4_THICKNESS_SEC2, 3/4_THICKNESS_SEC3, 3/4_THICKNESS_SEC4, 3/4_THICKNESS_SEC5, % AREA_SEC1, AREA_SEC2, AREA_SEC3, AREA_SEC4, AREA_SEC5, % AOA_SEC1, AOA_SEC2, AOA_SEC3, AOA_SEC4, AOA_SEC5, % CHORD_SEC1, CHORD_SEC2, CHORD_SEC3, CHORD_SEC4, CHORD_SEC5 % % Available design variables % HICKS_HENNE ( 1, Scale | Mark. List | Lower(0)/Upper(1) side, x_Loc ) % COSINE_BUMP ( 2, Scale | Mark. List | Lower(0)/Upper(1) side, x_Loc, x_Size ) % SPHERICAL ( 3, Scale | Mark. List | ControlPoint_Index, Theta_Disp, R_Disp ) % NACA_4DIGITS ( 4, Scale | Mark. List | 1st digit, 2nd digit, 3rd and 4th digit ) % DISPLACEMENT ( 5, Scale | Mark. List | x_Disp, y_Disp, z_Disp ) % ROTATION ( 6, Scale | Mark. List | x_Axis, y_Axis, z_Axis, x_Turn, y_Turn, z_Turn ) % FFD_CONTROL_POINT ( 7, Scale | Mark. List | FFD_BoxTag, i_Ind, j_Ind, k_Ind, x_Mov, y_Mov, z_Mov ) % FFD_DIHEDRAL_ANGLE ( 8, Scale | Mark. List | FFD_BoxTag, x_Orig, y_Orig, z_Orig, x_End, y_End, z_End ) % FFD_TWIST_ANGLE ( 9, Scale | Mark. List | FFD_BoxTag, x_Orig, y_Orig, z_Orig, x_End, y_End, z_End ) % FFD_ROTATION ( 10, Scale | Mark. List | FFD_BoxTag, x_Orig, y_Orig, z_Orig, x_End, y_End, z_End ) % FFD_CAMBER ( 11, Scale | Mark. List | FFD_BoxTag, i_Ind, j_Ind ) % FFD_THICKNESS ( 12, Scale | Mark. List | FFD_BoxTag, i_Ind, j_Ind ) % FFD_VOLUME ( 13, Scale | Mark. List | FFD_BoxTag, i_Ind, j_Ind ) % FOURIER ( 14, Scale | Mark. List | Lower(0)/Upper(1) side, index, cos(0)/sin(1) ) % % Optimization objective function with scaling factor % ex= Objective * Scale OPT_OBJECTIVE= DRAG * 0.1 % % Optimization constraint functions with scaling factors, separated by semicolons % ex= (Objective = Value ) * Scale, use '>','<','=' OPT_CONSTRAINT = NONE % Maximum number of iterations OPT_ITERATIONS= 100 % % Requested accuracy OPT_ACCURACY= 1E-9 % % Upper bound for each design variable OPT_BOUND_UPPER= 0.1 % % Lower bound for each design variable OPT_BOUND_LOWER= -0.1 % % Optimization design variables, separated by semicolons DEFINITION_DV= ( 7, 1.0 | car_ground | EngineCoverTip, 0, 0, 0, 1.0, 0.0, 0.0 );...... Last edited by dolfyn; May 9, 2017 at 02:25. Reason: Not stated clearly |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
decomposePar problem: Cell 0contains face labels out of range | vaina74 | OpenFOAM Pre-Processing | 37 | July 20, 2020 06:38 |
Shape Optimization with PBC | jubroll | SU2 Shape Design | 1 | June 12, 2015 01:38 |
Another shape optimization tutorial error | grjmpower | SU2 Shape Design | 6 | February 7, 2013 04:12 |
fluent add additional zones for the mesh file | SSL | FLUENT | 2 | January 26, 2008 12:55 |
Icemcfd 11: Loss of mesh from surface mesh option? | Joe | CFX | 2 | March 26, 2007 19:10 |