# Inverse Design Optimization

 Register Blogs Members List Search Today's Posts Mark Forums Read

 June 20, 2019, 03:37 Inverse Design Optimization #1 New Member   kevin havart Join Date: Jun 2019 Posts: 8 Rep Power: 6 Hello, I'm currently running an inverse design optimization on a nozzle. My goal pressure distribution come from said to be optimized nozzle. (Nozzle_Opti.png) I copy the pressure distribution at the outlet. Then I create a new nozzle which is not optimized (Nozzle_Base.png) and set the TargetCp.dat as following. Code: ```"Global_Index", "x_coord", "y_coord","Pressure","Pressure_Coefficient","Mach_Number" 3, 0.123, 0, 15000,-2.14844, 2.58 2, 0.123, 0.0143, 15000, -2.14844, 2.58 223, 0.123, 0.000595833, 15000, -2.14844, 2.58 224, 0.123, 0.00119167, 15000, -2.14844, 2.58 225, 0.123, 0.0017875, 15000, -2.14844, 2.58 226, 0.123, 0.00238333, 15000, -2.14844, 2.58 227, 0.123, 0.00297917, 15000, -2.14844, 2.58 228, 0.123, 0.003575, 15000, -2.14844, 2.58 229, 0.123, 0.00417083, 15000, -2.14844, 2.58 230, 0.123, 0.00476667, 15000, -2.14844, 2.58 231, 0.123, 0.0053625, 15000, -2.14844, 2.58 232, 0.123, 0.00595833, 15000, -2.14844, 2.58 233, 0.123, 0.00655417, 15000, -2.14844, 2.58 234, 0.123, 0.00715, 15000, -2.14844, 2.58 235, 0.123, 0.00774583, 15000, -2.14844, 2.58 236, 0.123, 0.00834167, 15000, -2.14844, 2.58 237, 0.123, 0.0089375, 15000, -2.14844, 2.58 238, 0.123, 0.00953333, 15000, -2.14844, 2.58 239, 0.123, 0.0101292, 15000, -2.14844, 2.58 240, 0.123, 0.010725, 15000, -2.14844, 2.58 241, 0.123, 0.0113208, 15000, -2.14844, 2.58 242, 0.123, 0.119167, 15000, -2.14844, 2.58 243, 0.123, 0.0125125, 15000, -2.14844, 2.58 244, 0.123, 0.0131083, 15000, -2.14844, 2.58 245, 0.123, 0.0137042, 15000, -2.14844, 2.58``` Then my config files is : Code: ```% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------% PHYSICAL_PROBLEM= EULER MATH_PROBLEM= DIRECT REGIME_TYPE= COMPRESSIBLE AXISYMMETRIC= NO RESTART_SOL= NO SYSTEM_MEASUREMENTS= SI % -------------------- COMPRESSIBLE FREE-STREAM DEFINITION --------------------% MACH_NUMBER= 0.8 AOA= 0 SIDESLIP_ANGLE= 0.0 INIT_OPTION= TD_CONDITIONS FREESTREAM_OPTION= TEMPERATURE_FS FREESTREAM_PRESSURE= 400000 FREESTREAM_TEMPERATURE= 288 REYNOLDS_LENGTH= 0.005 % % ---------------------- REFERENCE VALUE DEFINITION ---------------------------% REF_ORIGIN_MOMENT_X = 0.25 REF_ORIGIN_MOMENT_Y = 0.00 REF_ORIGIN_MOMENT_Z = 0.00 REF_LENGTH= 1.0 REF_AREA= 1.0 REF_DIMENSIONALIZATION= DIMENSIONAL % ---- IDEAL GAS, POLYTROPIC, VAN DER WAALS AND PENG ROBINSON CONSTANTS -------% FLUID_MODEL= IDEAL_GAS GAMMA_VALUE= 1.4 GAS_CONSTANT= 287.058 CRITICAL_TEMPERATURE= 131.00 CRITICAL_PRESSURE= 3588550.0 ACENTRIC_FACTOR= 0.035 % -------------------- BOUNDARY CONDITION DEFINITION --------------------------% MARKER_EULER= (Nozzle) MARKER_RIEMANN= (Inlet, TOTAL_CONDITIONS_PT, 400000.000000, 288, 1, 0.0, 0.0, Outlet, STATIC_PRESSURE, 20000.000000, 0.0, 0.0, 0.0, 0.0) MARKER_SYM= (Symmetry) % ------------------------ SURFACES IDENTIFICATION ----------------------------% MARKER_PLOTTING= ( Nozzle ) MARKER_MONITORING= ( Outlet ) % ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------% NUM_METHOD_GRAD= GREEN_GAUSS CFL_NUMBER= 2.0 CFL_ADAPT= YES CFL_ADAPT_PARAM= ( 1.5, 0.5, 0.5, 1000) MAX_DELTA_TIME= 1E6 RK_ALPHA_COEFF= ( 0.66667, 0.66667, 1.000000 ) % ----------- SLOPE LIMITER AND DISSIPATION SENSOR DEFINITION -----------------% MUSCL_FLOW= YES SLOPE_LIMITER_FLOW= NONE MUSCL_ADJFLOW= YES SLOPE_LIMITER_ADJFLOW= VENKATAKRISHNAN ADJ_SHARP_LIMITER_COEFF= 3.0 LIMITER_ITER= 999999 % --------------------- INVERSE DESIGN SIMULATION -----------------------------% INV_DESIGN_CP= YES INV_DESIGN_HEATFLUX= NO % ------------------------ LINEAR SOLVER DEFINITION ---------------------------% LINEAR_SOLVER= FGMRES LINEAR_SOLVER_PREC= LU_SGS LINEAR_SOLVER_ERROR= 1E-4 LINEAR_SOLVER_ITER= 4 % -------------------------- MULTIGRID PARAMETERS -----------------------------% MGLEVEL= 2 MGCYCLE= V_CYCLE MG_PRE_SMOOTH= ( 1, 2, 3, 3 ) MG_POST_SMOOTH= ( 0, 0, 0, 0 ) MG_CORRECTION_SMOOTH= ( 0, 0, 0, 0 ) MG_DAMP_RESTRICTION= 1.0 MG_DAMP_PROLONGATION= 1.0 % -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------% CONV_NUM_METHOD_FLOW= JST ENTROPY_FIX_COEFF= 0.0 TIME_DISCRE_FLOW= EULER_IMPLICIT RELAXATION_FACTOR_FLOW= 0.95 JST_SENSOR_COEFF= ( 0.5, 0.02 ) % ---------------- ADJOINT-FLOW NUMERICAL METHOD DEFINITION -------------------% FROZEN_LIMITER_DISC= NO FROZEN_VISC_DISC= NO INCONSISTENT_DISC= NO CONV_NUM_METHOD_ADJFLOW= JST TIME_DISCRE_ADJFLOW= EULER_IMPLICIT RELAXATION_FACTOR_ADJFLOW= 1.0 CFL_REDUCTION_ADJFLOW= 0.8 LIMIT_ADJFLOW= 1E6 MG_ADJFLOW= YES % --------------------------- CONVERGENCE PARAMETERS --------------------------% EXT_ITER= 10000 CONV_CRITERIA= RESIDUAL RESIDUAL_REDUCTION= 7 RESIDUAL_MINVAL= -20 STARTCONV_ITER= 10 % ------------------------- INPUT/OUTPUT INFORMATION --------------------------% MESH_FILENAME= Mesh_FFD.su2 MESH_FORMAT= SU2 MESH_OUT_FILENAME= mesh_out.su2 SOLUTION_FLOW_FILENAME= solution_flow.dat SOLUTION_ADJ_FILENAME= solution_adj.dat OUTPUT_FORMAT= PARAVIEW CONV_FILENAME= history BREAKDOWN_FILENAME= forces_breakdown.dat RESTART_FLOW_FILENAME= restart_flow.dat RESTART_ADJ_FILENAME= restart_adj.dat VOLUME_FLOW_FILENAME= flow VOLUME_ADJ_FILENAME= adjoint VALUE_OBJFUNC_FILENAME= of_eval.dat GRAD_OBJFUNC_FILENAME= of_grad.dat SURFACE_FLOW_FILENAME= surface_flow SURFACE_ADJ_FILENAME= surface_adjoint WRT_SOL_FREQ= 250 WRT_SOL_FREQ_DUALTIME= 1 WRT_CON_FREQ= 1 WRT_CON_FREQ_DUALTIME= 10 WRT_RESIDUALS= NO WRT_LIMITERS= NO WRT_SHARPEDGES= NO WRT_SURFACE= NO LOW_MEMORY_OUTPUT= NO CONSOLE_OUTPUT_VERBOSITY= HIGH WRT_BINARY_RESTART= YES READ_BINARY_RESTART= YES REORIENT_ELEMENTS= YES % ----------------------- DESIGN VARIABLE PARAMETERS --------------------------% DV_MARKER= (Nozzle) %DV_KIND= FFD_SETTING %DV_PARAM=(NOZZLE,1,1,0.0,1.0) %DV_VALUE= 0.1 DV_KIND=FFD_CONTROL_POINT_2D,FFD_CONTROL_POINT_2D,FFD_CONTROL_POINT_2D,FFD_CONTROL_POINT_2D,FFD_CONTROL_POINT_2D,FFD_CONTROL_POINT_2D,FFD_CONTROL_POINT_2D,FFD_CONTROL_POINT_2D,FFD_CONTROL_POINT_2D,FFD_CONTROL_POINT_2D DV_PARAM=(NOZZLE,1,1,0.0,1.0);(NOZZLE,2,1,0.0,1.0);(NOZZLE,3,1,0.0,1.0);(NOZZLE,4,1,0.0,1.0);(NOZZLE,5,1,0.0,1.0);(NOZZLE,6,1,0.0,1.0);(NOZZLE,7,1,0.0,1.0);(NOZZLE,8,1,0.0,1.0);(NOZZLE,9,1,0.0,1.0);(NOZZLE,10,1,0.0,1.0) DV_VALUE= 0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001 % -------------------- FREE-FORM DEFORMATION PARAMETERS -----------------------% FFD_TOLERANCE= 1E-10 FFD_ITERATIONS= 1000 FFD_DEFINITION= (NOZZLE, 0.0559, 0.0, 0.0, 0.125, 0.0, 0.0, 0.125, 0.05, 0.0, 0.0559, 0.05, 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_DEGREE= (10, 1, 0) FFD_CONTINUITY= NO_DERIVATIVE FFD_SYMMETRY_PLANE= YES % ------------------------ GRID DEFORMATION PARAMETERS ------------------------% DEFORM_LINEAR_ITER= 1000 DEFORM_NONLINEAR_ITER= 2 DEFORM_CONSOLE_OUTPUT= YES DEFORM_LINEAR_SOLVER_ERROR= 1E-14 DEFORM_STIFFNESS_TYPE= CONSTANT_STIFFNESS DEFORM_LINEAR_SOLVER= FGMRES DEFORM_LINEAR_SOLVER_PREC= LU_SGS DEFORM_TOL_FACTOR = 0.000001 DEFORM_COEFF = 1E6 VISUALIZE_DEFORMATION= YES % --------------------- OPTIMAL SHAPE DESIGN DEFINITION -----------------------% OPT_OBJECTIVE=INVERSE_DESIGN_PRESSURE OPT_CONSTRAINT= NONE OPT_GRADIENT_FACTOR= 1E-6 OPT_RELAX_FACTOR= 1E3 OPT_ITERATIONS= 40 OPT_ACCURACY= 1E-20 OPT_BOUND_UPPER= 0.1 OPT_BOUND_LOWER= 0 FIN_DIFF_STEP=0.0001 DEFINITION_DV= ( 15, 1 | Nozzle | NOZZLE, 1, 1, 0.0, 1.0); ( 15, 1 | Nozzle | NOZZLE, 2, 1, 0.0, 1.0); ( 15, 1 | Nozzle | NOZZLE, 3, 1, 0.0, 1.0 ); ( 15, 1 | Nozzle | NOZZLE, 4, 1, 0.0, 1.0 );( 15, 1 | Nozzle | NOZZLE, 5, 1, 0.0, 1.0); ( 15, 0.1 | Nozzle | NOZZLE, 6, 1, 0.0, 1.0); ( 15, 0.1 | Nozzle | NOZZLE, 7, 1, 0.0, 1.0 ); ( 15, 1 | Nozzle | NOZZLE, 8, 1, 0.0, 1.0 );( 15, 1 | Nozzle | NOZZLE, 9, 1, 0.0, 1.0 );( 15, 1 | Nozzle | NOZZLE, 10, 1, 0.0, 1.0 )``` I'm using FFD boxes for the deformation and I have already observed the deformations. The problem is that when I run the optimization I can't recover the optimized design and it stops after only 1 iteration, saying : Code: ```Found: Mesh_FFD.su2 Found: TargetCp.dat New Project: ./ Removing old designs in 10s. Done! Sequential Least SQuares Programming (SLSQP) parameters: Number of design variables: 10 ( 10 ) Objective function scaling factor: [1.0] Maximum number of iterations: 40 Requested accuracy: 1e-26 Initial guess for the independent variable(s): [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0] Lower and upper bound for each independent variable: [(0.0, 0.0001), (0.0, 0.0001), (0.0, 0.0001), (0.0, 0.0001), (0.0, 0.0001), (0.0, 0.0001), (0.0, 0.0001), (0.0, 0.0001), (0.0, 0.0001), (0.0, 0.0001)] NIT FC OBJFUN GNORM 5 1 2.557944E-07 7.675535E-04 Positive directional derivative for linesearch (Exit mode 8) Current function value: 2.55794353e-07 Iterations: 5 Function evaluations: 1``` Have one of you ever had this problem ? Thanks for your help, Kevin Nozzle_Opti.png Nozzle_Base.jpg

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are Off Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post ericclaro SU2 Shape Design 1 May 30, 2019 18:53 dominico SU2 Shape Design 1 September 19, 2017 08:56 Bartosz.Ziegler SU2 Shape Design 10 February 17, 2016 12:50 Simon Nava Main CFD Forum 8 May 22, 2007 14:55 Arnold Free Main CFD Forum 0 August 10, 1999 10:18

All times are GMT -4. The time now is 17:17.