
[Sponsors] 
May 31, 2016, 15:34 
3D Windturbine simulation in SU2

#1 
Member
Kisorthman Vimalakanthan
Join Date: Apr 2011
Posts: 45
Rep Power: 13 
Dear All,
I'm trying to simulate a single blade 120deg axisymmetric model of the CX100 turbine using the incompressible RANS solver in SU2. The 3D grid was made in Pointwise using this webinar  https://www.youtube.com/watch?v=_o0KOJ7RJXc. The SU2 configuration file is also attached. You can download the grid using the link below or if this link expires, feel free to contact me. If anyone knows, please let me know if you know a place where I can share this indefinitely.. https://we.tl/Nz7L2hjxrv The attached images shows the boundary patches. With the current setup I'm not able to even reach one iteration. I've tried using the SU2_MSH processor and running the simulation with 'CX100_pre.su2' grid, which also fails. Could anyone please have a look at this case and suggest improvements that can actually enable me to run this problem. The .cfg file: Code:
%  DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION % % % Physical governing equations (EULER, NAVIER_STOKES, % WAVE_EQUATION, HEAT_EQUATION, LINEAR_ELASTICITY, % POISSON_EQUATION) PHYSICAL_PROBLEM= NAVIER_STOKES % % Specify turbulent model (NONE, SA, SA_NEG, SST) KIND_TURB_MODEL= SST % % Mathematical problem (DIRECT, CONTINUOUS_ADJOINT) MATH_PROBLEM= DIRECT % % Regime type (COMPRESSIBLE, INCOMPRESSIBLE, FREESURFACE) REGIME_TYPE= INCOMPRESSIBLE % % Restart solution (NO, YES) RESTART_SOL= NO %  INCOMPRESSIBLE FREESTREAM DEFINITION % % % Freestream density (1.2886 Kg/m^3 (air), 998.2 Kg/m^3 (water)) FREESTREAM_DENSITY= 1.2886 % % Freestream velocity (m/s) FREESTREAM_VELOCITY= ( 0.000, 0.000 , 4.000) % % Freestream viscosity (1.853E5 Ns/m^2 (air), 0.798E3 Ns/m^2 (water)) FREESTREAM_VISCOSITY= 1.853e05 %  REFERENCE VALUE DEFINITION % % % Reference origin for moment computation REF_ORIGIN_MOMENT_X = 0.00 REF_ORIGIN_MOMENT_Y = 0.00 REF_ORIGIN_MOMENT_Z = 0.00 % % Reference length for pitching, rolling, and yawing nondimensional moment REF_LENGTH_MOMENT= 1 % % Reference area for force coefficients (0 implies automatic calculation) REF_AREA= 0 %  BOUNDARY CONDITION DEFINITION % % % NavierStokes wall boundary marker(s) (NONE = no marker) MARKER_HEATFLUX= ( Blade, 0.0, Spinner, 0.0) % % Farfield boundary marker(s) (NONE = no marker) MARKER_FAR= ( Inlet, Outlet, TopWall ) % % Periodic boundary marker(s) (NONE = no marker) % Format: ( periodic marker, donor marker, rot_cen_x, rot_cen_y, rot_cen_z, rot_angle_xaxis, rot_angle_yaxis, rot_angle_zaxis, translation_x, translation_y, translation_z) MARKER_PERIODIC= ( Per1, Per2, 0.0, 0.0, 0.0, 0.0, 0.0, 120, 0.0, 0.0, 0.0) % % Marker(s) of the surface to be plotted or designed MARKER_PLOTTING= ( Blade ) % % Marker(s) of the surface where the functional (Cd, Cl, etc.) will be evaluated MARKER_MONITORING= ( Blade ) %  GRID ADAPTATION STRATEGY % % % Kind of grid adaptation (NONE, PERIODIC) KIND_ADAPT= PERIODIC % %  DYNAMIC MESH DEFINITION % % % Dynamic mesh simulation (NO, YES) GRID_MOVEMENT= YES % % Type of dynamic mesh (NONE, ROTATING_FRAME) GRID_MOVEMENT_KIND= ROTATING_FRAME % % Motion mach number (nondimensional). Used for intitializing a viscous flow % with the Reynolds number and for computing force coeffs. with dynamic meshes. % MACH_MOTION= 0.01 % % Coordinates of the motion origin MOTION_ORIGIN_X= 0.0 MOTION_ORIGIN_Y= 0.0 MOTION_ORIGIN_Z= 0.0 % % Angular velocity vector (rad/s) about the motion origi. Example 1250 RPM > 130.89969389957471 rad/s ROTATION_RATE_X = 0.0 ROTATION_RATE_Y = 0.0 ROTATION_RATE_Z = 5.811 % %  COMMON PARAMETERS DEFINING THE NUMERICAL METHOD % % % Numerical method for spatial gradients (GREEN_GAUSS, WEIGHTED_LEAST_SQUARES) NUM_METHOD_GRAD= WEIGHTED_LEAST_SQUARES % % CourantFriedrichsLewy condition of the finest grid CFL_NUMBER= 1 % % Adaptive CFL number (NO, YES) CFL_ADAPT= NO % % Parameters of the adaptive CFL number (factor down, factor up, CFL min value, % CFL max value ) CFL_ADAPT_PARAM= ( 2, 4, 500.0, 1000.0 ) %  SLOPE LIMITER DEFINITION % % % Reference element length for computing the slope and sharp edges limiters. REF_ELEM_LENGTH= 0.1 % % Coefficient for the limiter LIMITER_COEFF= 0.1 % % Coefficient for the sharp edges limiter SHARP_EDGES_COEFF= 3.0 % % Reference coefficient (sensitivity) for detecting sharp edges. REF_SHARP_EDGES= 3.0 % % Remove sharp edges from the sensitivity evaluation (NO, YES) SENS_REMOVE_SHARP= NO %  LINEAR SOLVER DEFINITION % % % Linear solver for implicit formulations (BCGSTAB, FGMRES) LINEAR_SOLVER= FGMRES % % Preconditioner of the Krylov linear solver (JACOBI, LINELET, LU_SGS) LINEAR_SOLVER_PREC= LU_SGS % % Minimum error of the linear solver for implicit formulations LINEAR_SOLVER_ERROR= 1E4 % % Max number of iterations of the linear solver for the implicit formulation LINEAR_SOLVER_ITER= 10 %  FLOW NUMERICAL METHOD DEFINITION % % % Convective numerical method (JST, LAXFRIEDRICH, CUSP, ROE, AUSM, HLLC, % TURKEL_PREC, MSW) CONV_NUM_METHOD_FLOW= ROE % % Spatial numerical order integration (1ST_ORDER, 2ND_ORDER, 2ND_ORDER_LIMITER) % SPATIAL_ORDER_FLOW= 2ND_ORDER_LIMITER % % Slope limiter (VENKATAKRISHNAN, MINMOD) SLOPE_LIMITER_FLOW= VENKATAKRISHNAN % % 1st, 2nd and 4th order artificial dissipation coefficients AD_COEFF_FLOW= ( 0.15, 0.5, 0.02 ) % % Time discretization (RUNGEKUTTA_EXPLICIT, EULER_IMPLICIT, EULER_EXPLICIT) TIME_DISCRE_FLOW= EULER_IMPLICIT %  TURBULENT NUMERICAL METHOD DEFINITION % % % Convective numerical method (SCALAR_UPWIND) CONV_NUM_METHOD_TURB= SCALAR_UPWIND % % Spatial numerical order integration (1ST_ORDER, 2ND_ORDER, 2ND_ORDER_LIMITER) % SPATIAL_ORDER_TURB= 1ST_ORDER % % Slope limiter (VENKATAKRISHNAN, MINMOD) SLOPE_LIMITER_TURB= VENKATAKRISHNAN % % Time discretization (EULER_IMPLICIT) TIME_DISCRE_TURB= EULER_IMPLICIT %  CONVERGENCE PARAMETERS % % % % Number of total iterations EXT_ITER= 4000 % % Convergence criteria (CAUCHY, RESIDUAL) % CONV_CRITERIA= RESIDUAL % % Residual reduction (order of magnitude with respect to the initial value) RESIDUAL_REDUCTION= 8 % % Min value of the residual (log10 of the residual) RESIDUAL_MINVAL= 9 % % Start convergence criteria at iteration number STARTCONV_ITER= 10 % % Number of elements to apply the criteria CAUCHY_ELEMS= 100 % % Epsilon to control the series convergence CAUCHY_EPS= 1E4 % % Function to apply the criteria (LIFT, DRAG, NEARFIELD_PRESS, SENS_GEOMETRY, % SENS_MACH, DELTA_LIFT, DELTA_DRAG) CAUCHY_FUNC_FLOW= LIFT %  INPUT/OUTPUT INFORMATION % % % Mesh input file MESH_FILENAME= CX100.su2 % % Mesh input file format (SU2, CGNS, NETCDF_ASCII) MESH_FORMAT= SU2 % % Mesh output file MESH_OUT_FILENAME= CX100_per.su2 % % Restart flow input file SOLUTION_FLOW_FILENAME= restart_flow.dat % % Output file format (PARAVIEW, TECPLOT, STL) OUTPUT_FORMAT= TECPLOT % % Output file convergence history (w/o extension) CONV_FILENAME= history % % Output file restart flow RESTART_FLOW_FILENAME= restart_flow.dat % % % Output file flow (w/o extension) variables VOLUME_FLOW_FILENAME= Results % % Output file surface flow coefficient (w/o extension) SURFACE_FLOW_FILENAME= surface_flow % % Writing solution file frequency WRT_SOL_FREQ= 500 % % Writing convergence history frequency WRT_CON_FREQ= 1 Any advice is greatly appreciated Thanks in advance, Kind regards, Kishore
__________________
Kisorthman Vimalakanthan Dept. of Power and Propulsion Cranfield University Email: k.vimalakanthan@gmail.com Last edited by k.vimalakanthan; June 1, 2016 at 04:52. 

May 31, 2016, 15:53 

#2 
New Member
Abhinandan Kumar Singh
Join Date: Apr 2016
Location: Germany
Posts: 4
Rep Power: 8 
hey,
Could you please elaborate what error you got. 

June 1, 2016, 05:27 

#3 
Member
Kisorthman Vimalakanthan
Join Date: Apr 2011
Posts: 45
Rep Power: 13 
This is the output when I run the case as it is:
Code:
  ___ _ _ ___   / __   _ ) Release 4.1.1 "Cardinal"   \__ \ _ / /   ___/\___//___ Suite (Computational Fluid Dynamics Code)      SU2 Lead Dev.: Dr. Francisco Palacios, Francisco.D.Palacios@boeing.com  Dr. Thomas D. Economon, economon@stanford.edu    SU2 Developers:    Prof. Juan J. Alonso's group at Stanford University.    Prof. Piero Colonna's group at Delft University of Technology.    Prof. Nicolas R. Gauger's group at Kaiserslautern U. of Technology.    Prof. Alberto Guardone's group at Polytechnic University of Milan.    Prof. Rafael Palacios' group at Imperial College London.    Copyright (C) 20122016 SU2, the opensource CFD code.     SU2 is free software; you can redistribute it and/or   modify it under the terms of the GNU Lesser General Public   License as published by the Free Software Foundation; either   version 2.1 of the License, or (at your option) any later version.     SU2 is distributed in the hope that it will be useful,   but WITHOUT ANY WARRANTY; without even the implied warranty of   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU   Lesser General Public License for more details.     You should have received a copy of the GNU Lesser General Public   License along with SU2. If not, see <http://www.gnu.org/licenses/>.    Physical Case Definition  Incompressible RANS equations. Turbulence model: Menter's SST Performing a dynamic mesh simulation: rotating reference frame. No restart solution, use the values at infinity (freestream). Dimensional simulation. The reference length/area will be computed using y(2D) or z(3D) projection. The reference length (moment computation) is 1. Reference origin (moment computation) is (0, 0, 0). Surface(s) where the force coefficients are evaluated: Blade. Surface(s) plotted in the output file: Blade. Surface(s) belonging to the FluidStructure Interaction problem: Input mesh file name: CX100.su2  Space Numerical Integration  Roe (with entropy fix) solver for the flow inviscid terms. Second order integration with slope limiter. Venkatakrishnan slopelimiting method, with constant: 0.1. The reference element size is: 0.1. Scalar upwind solver (first order) for the turbulence model. First order integration. Average of gradients with correction (viscous flow terms). Average of gradients with correction (viscous turbulence terms). Gradient Computation using weighted LeastSquares method. Artificial compressibility factor: 1.  Time Numerical Integration  Local time stepping (steady state simulation). Euler implicit method for the flow equations. No CFL adaptation. CourantFriedrichsLewy number: 1 Euler implicit time integration for the turbulence model.  Convergence Criteria  Maximum number of iterations: 4000. Reduce the density residual 8 orders of magnitude. The minimum bound for the density residual is 10^(9). Start convergence criteria at iteration 10.  Output Information  Writing a flow solution every 500 iterations. Writing the convergence history every 1 iterations. The output file format is Tecplot ASCII (.dat). Convergence history file name: history. Forces breakdown file name: forces_breakdown.dat. Surface flow coefficients file name: surface_flow. Flow variables file name: Results. Restart flow file name: restart_flow.dat.  Config File Boundary Information  Farfield boundary marker(s): Inlet, Outlet, TopWall. Periodic boundary marker(s): Per1, Per2. Constant heat flux wall boundary marker(s): Blade, Spinner.  Read Grid File Information  Three dimensional problem. 2759777 points. 7 surface markers. 34994 boundary elements in index 0 (Marker = Blade). 1318 boundary elements in index 1 (Marker = Inlet). 1322 boundary elements in index 2 (Marker = Outlet). 16969 boundary elements in index 3 (Marker = Per1). 16969 boundary elements in index 4 (Marker = Per2). 5388 boundary elements in index 5 (Marker = Spinner). 8056 boundary elements in index 6 (Marker = TopWall). 10565786 tetrahedra. 435120 hexahedra. 1035964 prisms. 63800 pyramids.  Geometry Preprocessing  Setting point connectivity. Renumbering points (Reverse Cuthill McKee Ordering). Recomputing point connectivity. Setting element connectivity. Checking the numerical grid orientation. Identifying edges and vertices. Computing centers of gravity. Setting the control volume structure. Volume of the computational grid: 2.01e+05. Searching for the closest normal neighbors to the surfaces. Compute the surface curvature. Max K: 1.23e+04. Mean K: 11.2. Standard deviation K: 195. Computing wall distances. Area projection in the zplane = 4.97.  Driver Preprocessing  Instantiating a single zone driver for the problem.  Iteration Preprocessing  Zone 1: Euler/NavierStokes/RANS flow iteration.  Solver Preprocessing  Force coefficients computed using MACH_MOTION. Viscous and Inviscid flow: rho_ref, and vel_ref are based on the freestream values, p_ref = rho_ref*vel_ref^2. The freestream value of the pressure is 0. Mach number: 0.0120497, computed using the Bulk modulus. Angle of attack (deg): 90, computed using the the freestream velocity. Side slip angle (deg): 0, computed using the the freestream velocity. Reynolds number: 278165, computed using freestream values. Only dimensional computation, the grid should be dimensional.  Input conditions: Bulk modulus: 142000 Pa. Artificial compressibility factor: 1 Pa. Freestream static pressure: 0 Pa. Freestream total pressure: 0 Pa. Freestream density: 1.2886 kg/m^3. Freestream velocity: (0, 0, 4) m/s. Magnitude: 4 m/s. Freestream viscosity: 1.853e05 N.s/m^2. Freestream turb. kinetic energy per unit mass: 0.06 m^2/s^2. Freestream specific dissipation: 417.248 1/s.  Reference values: Reference pressure: 20.6176 Pa. Reference density: 1.2886 kg/m^3. Reference velocity: 4 m/s. Reference length: 1 m. Reference viscosity: 5.1544 N.s/m^2.  Resulting nondimensional state: Mach number (nondim): 0.0120497 Reynolds number (nondim): 278165. Re length: 1 m. Freestream pressure (nondim): 0 Freestream density (nondim): 1 Freestream velocity (nondim): (0, 0, 1). Magnitude: 1 Freestream viscosity (nondim): 3.59499e06 Freestream turb. kinetic energy (nondim): 0.00375 Freestream specific dissipation (nondim): 104.312 Initialize Jacobian structure (NavierStokes). MG level: 0. Initialize Jacobian structure (SST model).  Integration and Numerics Preprocessing  Integration Preprocessing. Numerics Preprocessing. Setting dynamic mesh structure.  Begin Solver  Setting rotating frame grid velocities for zone 0. Rotational origin (x, y, z): ( 0, 0, 0 ) Angular velocity about x, y, z axes: ( 0, 0, 2.76552e+52 ) rad/s Segmentation fault If I use SU2_MSH to adapt the CX100.su2 grid file for periodic boundary conditions the output is: Code:
  ___ _ _ ___   / __   _ ) Release 4.1.1 "Cardinal"   \__ \ _ / /   ___/\___//___ Suite (Mesh Adaptation Code)      SU2 Lead Dev.: Dr. Francisco Palacios, Francisco.D.Palacios@boeing.com  Dr. Thomas D. Economon, economon@stanford.edu    SU2 Developers:    Prof. Juan J. Alonso's group at Stanford University.    Prof. Piero Colonna's group at Delft University of Technology.    Prof. Nicolas R. Gauger's group at Kaiserslautern U. of Technology.    Prof. Alberto Guardone's group at Polytechnic University of Milan.    Prof. Rafael Palacios' group at Imperial College London.    Copyright (C) 20122016 SU2, the opensource CFD code.     SU2 is free software; you can redistribute it and/or   modify it under the terms of the GNU Lesser General Public   License as published by the Free Software Foundation; either   version 2.1 of the License, or (at your option) any later version.     SU2 is distributed in the hope that it will be useful,   but WITHOUT ANY WARRANTY; without even the implied warranty of   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU   Lesser General Public License for more details.     You should have received a copy of the GNU Lesser General Public   License along with SU2. If not, see <http://www.gnu.org/licenses/>.    Physical Case Definition  Input mesh file name: CX100.su2  Grid adaptation strategy  Grid modification to run periodic bc problems.  Output Information  Output mesh file name: CX100_per.su2. Output mesh file name: CX100_per.su2. Restart flow file name: restart_flow.dat.  Config File Boundary Information  Farfield boundary marker(s): Inlet, Outlet, TopWall. Periodic boundary marker(s): Per1, Per2. Constant heat flux wall boundary marker(s): Blade, Spinner.  Read Grid File Information  Three dimensional problem. 2759777 points. 7 surface markers. 34994 boundary elements in index 0 (Marker = Blade). 1318 boundary elements in index 1 (Marker = Inlet). 1322 boundary elements in index 2 (Marker = Outlet). 16969 boundary elements in index 3 (Marker = Per1). 16969 boundary elements in index 4 (Marker = Per2). 5388 boundary elements in index 5 (Marker = Spinner). 8056 boundary elements in index 6 (Marker = TopWall). 10565786 tetrahedra. 435120 hexahedra. 1035964 prisms. 63800 pyramids.  Preprocessing computations  Setting local point and element connectivity. Check numerical grid orientation. Identify faces, edges and vertices. Set control volume structure. Volume of the computational grid: 201013.  Setting the periodic boundaries  Setting the periodic boundary conditions. Checking Per1 boundary against periodic donor, Per2. Bad match for point 1026073. Nearest donor distance: 8.5443461700e09. Bad match for point 1038733. Nearest donor distance: 8.5443461700e09. Bad match for point 1065889. Nearest donor distance: 1.7181944039e08. Bad match for point 1037906. Nearest donor distance: 1.7191530396e08. Bad match for point 1055474. Nearest donor distance: 1.7364121476e08. . . . Bad match for point 1189530. Nearest donor distance: 8.3956174894e09. Bad match for point 1190592. Nearest donor distance: 7.2492817214e09. Bad match for point 1189506. Nearest donor distance: 1.3692191254e08. Bad match for point 1189453. Nearest donor distance: 5.3120100478e09. !!! Warning !!! Bad matches found. Computation will continue, but be cautious. Number of new points: 11360. Number of new interior elements: 51634. Number of new boundary elements added to preexisting markers: 376. Completed in 165.0100000000 seconds on 1 core.  Exit Success (SU2_MSH)  Code:
  ___ _ _ ___   / __   _ ) Release 4.1.1 "Cardinal"   \__ \ _ / /   ___/\___//___ Suite (Computational Fluid Dynamics Code)      SU2 Lead Dev.: Dr. Francisco Palacios, Francisco.D.Palacios@boeing.com  Dr. Thomas D. Economon, economon@stanford.edu    SU2 Developers:    Prof. Juan J. Alonso's group at Stanford University.    Prof. Piero Colonna's group at Delft University of Technology.    Prof. Nicolas R. Gauger's group at Kaiserslautern U. of Technology.    Prof. Alberto Guardone's group at Polytechnic University of Milan.    Prof. Rafael Palacios' group at Imperial College London.    Copyright (C) 20122016 SU2, the opensource CFD code.     SU2 is free software; you can redistribute it and/or   modify it under the terms of the GNU Lesser General Public   License as published by the Free Software Foundation; either   version 2.1 of the License, or (at your option) any later version.     SU2 is distributed in the hope that it will be useful,   but WITHOUT ANY WARRANTY; without even the implied warranty of   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU   Lesser General Public License for more details.     You should have received a copy of the GNU Lesser General Public   License along with SU2. If not, see <http://www.gnu.org/licenses/>.    Physical Case Definition  Incompressible RANS equations. Turbulence model: Menter's SST Performing a dynamic mesh simulation: rotating reference frame. No restart solution, use the values at infinity (freestream). Dimensional simulation. The reference length/area will be computed using y(2D) or z(3D) projection. The reference length (moment computation) is 1. Reference origin (moment computation) is (0, 0, 0). Surface(s) where the force coefficients are evaluated: Blade. Surface(s) plotted in the output file: Blade. Surface(s) belonging to the FluidStructure Interaction problem: Input mesh file name: CX100_per.su2  Space Numerical Integration  Roe (with entropy fix) solver for the flow inviscid terms. Second order integration with slope limiter. Venkatakrishnan slopelimiting method, with constant: 0.1. The reference element size is: 0.1. Scalar upwind solver (first order) for the turbulence model. First order integration. Average of gradients with correction (viscous flow terms). Average of gradients with correction (viscous turbulence terms). Gradient Computation using weighted LeastSquares method. Artificial compressibility factor: 1.  Time Numerical Integration  Local time stepping (steady state simulation). Euler implicit method for the flow equations. No CFL adaptation. CourantFriedrichsLewy number: 1 Euler implicit time integration for the turbulence model.  Convergence Criteria  Maximum number of iterations: 4000. Reduce the density residual 8 orders of magnitude. The minimum bound for the density residual is 10^(9). Start convergence criteria at iteration 10.  Output Information  Writing a flow solution every 500 iterations. Writing the convergence history every 1 iterations. The output file format is Tecplot ASCII (.dat). Convergence history file name: history. Forces breakdown file name: forces_breakdown.dat. Surface flow coefficients file name: surface_flow. Flow variables file name: Results. Restart flow file name: restart_flow.dat.  Config File Boundary Information  Farfield boundary marker(s): Inlet, Outlet, TopWall. Periodic boundary marker(s): Per1, Per2. Constant heat flux wall boundary marker(s): Blade, Spinner.  Read Grid File Information  Three dimensional problem. 2749827 points and 21310 ghost points. 9 surface markers. 34994 boundary elements in index 0 (Marker = Blade). 1364 boundary elements in index 1 (Marker = Inlet). 1368 boundary elements in index 2 (Marker = Outlet). 16969 boundary elements in index 3 (Marker = Per1). 16969 boundary elements in index 4 (Marker = Per2). 5539 boundary elements in index 5 (Marker = Spinner). 8189 boundary elements in index 6 (Marker = TopWall). 2 periodic transformations. 10611732 tetrahedra. 435120 hexahedra. 1041592 prisms. 63860 pyramids.  Geometry Preprocessing  Setting point connectivity. Renumbering points (Reverse Cuthill McKee Ordering). Recomputing point connectivity. Setting element connectivity. Checking the numerical grid orientation. The surface element (1, 1325) doesn't have an associated volume element. Kind regards, Kishore
__________________
Kisorthman Vimalakanthan Dept. of Power and Propulsion Cranfield University Email: k.vimalakanthan@gmail.com 

June 1, 2016, 10:23 

#4 
Super Moderator
Tim Albring
Join Date: Sep 2015
Posts: 195
Rep Power: 9 
Hi Kishore,
as far as I know, the incompressible solver is currently not able to handle a rotating frame (at least the additional terms are not included in the flux computation). You could try to use the compressible solver. However, this will not solve your problem with the periodic boundaries. What exactly is the difference between the first run and the third one ? Tim 

June 1, 2016, 12:59 

#5  
Member
Kisorthman Vimalakanthan
Join Date: Apr 2011
Posts: 45
Rep Power: 13 
Quote:
Many thanks for your reply. The differences are that the first run was performed with the grid directly from pointwise 'CX100.su2' without any preprocessing for SU2 periodic simulation (halo cells, renumbering etc.), while the second run was performed with the grid outputted from the SU2_MSH 'CX100_pre.su2'  grid adapted for periodic boundaries. Have you or anyone successfully simulated a 3D periodic problem in SU2 using the grid from pointwise? Kind regards, Kishore
__________________
Kisorthman Vimalakanthan Dept. of Power and Propulsion Cranfield University Email: k.vimalakanthan@gmail.com 

June 1, 2016, 13:38 

#6  
New Member
Join Date: Dec 2013
Location: Italy
Posts: 29
Rep Power: 10 
Hi,
I tried to run your grid, but the file seems different, as there are only 5 markers and they are named: MARKER_TAG= Blade MARKER_ELEMS= 34994 MARKER_TAG= Farfield MARKER_ELEMS= 8944 MARKER_TAG= Per1 MARKER_ELEMS= 7393 MARKER_TAG= Per2 MARKER_ELEMS= 7393 MARKER_TAG= Spinner MARKER_ELEMS= 5388 I modified cfg file accordingly and the output of SU2_MSH gave log similar to yours with the following ending lines: Quote:
I'm using release 4.1.2 

June 1, 2016, 14:15 

#7  
Member
Kisorthman Vimalakanthan
Join Date: Apr 2011
Posts: 45
Rep Power: 13 
Quote:
I realised that I uploaded the wrong grid file and updated the link above this morning with the grid shown in the images, maybe you were using the old link. Regardless as you've already found out, they only vary in the boundary names. Its interesting that your are getting up to the angular velocity output from SU2_CFD, so the SU2_MSH in v4.1.2 is better at adapting periodic boundaries for this grid than the v4.1.1. Because my run failed on "The surface element (1, 1325) doesn't have an associated volume element.". Did you get bad matches too from SU2_MSH? I'm glad we got progress with handling the grid, maybe as Tim said we need to try the compressible solver to get the rotating frame and the periodic BCs to work? My next step will be to get v4.1.2 and setup a compressible case..maybe with very low mach numbers to skip any compressibility effects?  adjust the fluid models so it computes the same density for all pressure and temperatures perhaps. maybe I can do isothermal compressible simulation? where we change the fluid model to set same density for all pressures...I'm sure there is more to it than that... Has anyone tried to compute a incompressible solution using the compressible SU2 solver? Please let me know your thoughts and if any other results.... Thanks a lot, Kind regards, Kishore
__________________
Kisorthman Vimalakanthan Dept. of Power and Propulsion Cranfield University Email: k.vimalakanthan@gmail.com 

June 1, 2016, 23:13 

#8  
New Member
Yang Muchen
Join Date: Jan 2015
Posts: 7
Rep Power: 9 
Quote:
For the bad matches in SU2_MSH, since the nearest donor distance is about 1e8 level, which is close to the truncation error of single precision, I think you could try to output the mesh with doubleprecision in Pointwise. Besides,you could also adjust the tolerance in function CPhysicalGeometry::SetPeriodicBoundary,which is currently 1e10, and changing it to 1e8 may be helpful in this case. Muchen 

June 2, 2016, 12:59 

#9  
Member
Kisorthman Vimalakanthan
Join Date: Apr 2011
Posts: 45
Rep Power: 13 
Quote:
Many thanks for your advice. The grids were generated with double doubleprecision, so maybe that is the maximum matches I could get from Pointwise? So from what I understand if I have grids with 1e8 or more bad match I can set CPhysicalGeometry::SetPeriodicBoundary to higher values to not allow CFD_MSH to modify the location of the grid elements but still renumber them for periodic BC? Kind regards, Kishore
__________________
Kisorthman Vimalakanthan Dept. of Power and Propulsion Cranfield University Email: k.vimalakanthan@gmail.com 

June 3, 2016, 03:20 

#10  
New Member
Yang Muchen
Join Date: Jan 2015
Posts: 7
Rep Power: 9 
Quote:
Actually, no matter what the tolerance is, for each point on the periodic boundaries, the SU2_MSH command would still find its nearest point after the transformation defined in the config file and set it as the donor point . Since this is a pointtopoint approach, it is best to create identical meshes on the two periodic boundaries, which means that these two meshes could completely overlap with each other after certain translation or rotation. Besides, you could always check the "periodic_halo.dat" to see if the halo cells are created successfully. Muchen 

June 16, 2016, 08:52 

#11  
Member
Kisorthman Vimalakanthan
Join Date: Apr 2011
Posts: 45
Rep Power: 13 
Quote:
In Pointwise, I did make surface mesh for both periodic boundaries to be identical (1:1)  the surface mesh from one side was rotated and the volumes were generated between them to make it easy for SU2_MSH. Maybe there is some tolerance issue with the copy/rotation in pointwise... Still haven't tried the compressible solver with this grid, hopefully with the newer version V4.2 everything run smoothly.... Kind regards, Kishore
__________________
Kisorthman Vimalakanthan Dept. of Power and Propulsion Cranfield University Email: k.vimalakanthan@gmail.com 

June 16, 2016, 09:41 

#12 
New Member
Join Date: Jun 2016
Posts: 7
Rep Power: 8 
Hi Kishore,
I am also working on the 3d Periodic BCs for rotor simulation using SU2. the link to download your grid isn't available anymore, could you please update it again? And I want try your run to see what happen. Thanks! 

June 17, 2016, 04:45 

#13  
Member
Kisorthman Vimalakanthan
Join Date: Apr 2011
Posts: 45
Rep Power: 13 
Quote:
I've uploaded it again: https://we.tl/peNWItztZt Maybe you know a better place to permanently share this file.... Can you please post your configuration file that we already have for a 3D periodic rotor cases? and any improvements you suggest for this case. Any help is greatly appreciated, Kind regards, Kishore
__________________
Kisorthman Vimalakanthan Dept. of Power and Propulsion Cranfield University Email: k.vimalakanthan@gmail.com 

June 20, 2016, 05:27 

#14 
New Member
Join Date: Jun 2016
Posts: 7
Rep Power: 8 
Hi Kishore,
Now I am trying to run the Caradonna&Tung case using 3D rotating periodic BC, the procedures are: 1. Output the su2 mesh from pointwise. 2. Set up the configuration file, the periodic boundary and grid adaption parts are: % Periodic boundary marker(s) (NONE = no marker) % Format: ( periodic marker, donor marker, rot_cen_x, rot_cen_y, rot_cen_z, rot_angle_xaxis, rot_angle_yaxis, rot_angle_zaxis, translation_x, translation_y, translation_z) MARKER_PERIODIC= ( per1, per2, 0.0, 0.0, 0.0, 0.0, 0.0, 180.0, 0.0, 0.0, 0.0) %  GRID ADAPTATION STRATEGY % % % Kind of grid adaptation (NONE, PERIODIC) KIND_ADAPT= PERIODIC 3. Execute the SU2_MSH to output su2 grid with periodic halo element 4. Running the SU2_CFD using new su2 grid from step 3. However, there are two types of error are occurred at the last step. The first one is new halo element boundary surface cannot find associated element, the error information is "Checking the numerical grid orientation. The surface element (1, 10691) doesn't have an associated volume element. " This surface (1,10691) is a new boundary surface created in step 3. I have tested several grids for this case, and some grids could pass the step where the first error was occurred, but still cannot be start running. I think the second error is about MPI sending/receiving, the error information is segmentation fault when I executed the SU2_CFD in parallel. And the code can run several steps in serial execution, but the segmentation error will be occurred when the master node merge coordinates at solution writing out step. These are the errors and information I got using 3d rotation periodic BC, and maybe we can discuss about it. Thanks! 

February 8, 2019, 14:43 
Final report on modelling Wind turbine using SU2 and OpenFOAM

#15 
Member
Kisorthman Vimalakanthan
Join Date: Apr 2011
Posts: 45
Rep Power: 13 
Hi All,
I've managed to put together the steps and settings for wind turbine simulation using SU2 and OpenFOAM. It also includes some validation for the mexico rotor. http://publications.tno.nl/publicati...018R11648.pdf https://www.slideshare.net/slideshow...n9VMov6E11jeAd I hope you find it useful, Kind regards, Kishore
__________________
Kisorthman Vimalakanthan Dept. of Power and Propulsion Cranfield University Email: k.vimalakanthan@gmail.com 

Tags 
su2 3d periodicity, su2_cfd, su2_msh 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
CGNS Boundary conditions using SU2  denzell  SU2  3  July 9, 2018 05:58 
Some questions about flow boiling simulation in Fluent  beastieboys6  FLUENT  8  November 20, 2017 23:47 
Control simulation to apply different fields with chtMultiRegionFoam  jmdf  OpenFOAM Running, Solving & CFD  0  February 29, 2016 07:05 
restarting paused transient simulation using reactingFoam  JMDag2004  OpenFOAM Running, Solving & CFD  1  August 10, 2015 10:15 
Incompressible simulation  brugiere_olivier  SU2  2  April 15, 2014 10:12 