|
[Sponsors] |
Problems with Time Dependent Simulation (Unsteady) |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
February 28, 2024, 13:44 |
Problems with Time Dependent Simulation (Unsteady)
|
#1 |
New Member
Jade Xie
Join Date: Feb 2024
Posts: 2
Rep Power: 0 |
Hi. I'm working on the fluid field around NACA0012 whose pitching angle changed periodically, and I hope to get the information on every single grid as time goes, but I failed.
Actually I get only one file named 'flow.vtu', it does include information for every grid but only shows a state that does not change by time, which means it has no time information. I don't know which time does it show, or it's displaying the average situation? I've set the TIME_DOMAIN = YES, and also TIME_STEP, MAX_TIME, but it doesn't work. I've pasted anything I think is related to time and output down there, so if anyone may find out what I'm missing or doing wrong in my code. If needed, I can paste all my code. And by the way, I'm visualizing the .vku by ParaView, and I'm not sure whether this is related to my problem. If you have any softare recommended, please tell me. Thanks a lot. Thank you in advance for any responses. -----------here is part of my code-------------------- -----------apologies for its untidiness---------------- Code:
% ------------------------------- SOLVER CONTROL ------------------------------% % % Maximum number of inner iterations INNER_ITER= 9999 % % Maximum number of outer iterations (only for multizone problems) OUTER_ITER= 1 % % Maximum number of time iterations TIME_ITER= 1 % % Convergence field CONV_FIELD= DRAG % % Min value of the residual (log10 of the residual) CONV_RESIDUAL_MINVAL= -10 % % Start convergence criteria at iteration number CONV_STARTITER= 10 % % Number of elements to apply the criteria CONV_CAUCHY_ELEMS= 100 % % Epsilon to control the series convergence CONV_CAUCHY_EPS= 1E-6 % % Iteration number to begin unsteady restarts RESTART_ITER= 0 % % ------------------------- TIME-DEPENDENT SIMULATION -------------------------------% % % Time domain simulation TIME_DOMAIN= YES % % Unsteady simulation (NO, TIME_STEPPING, DUAL_TIME_STEPPING-1ST_ORDER, % DUAL_TIME_STEPPING-2ND_ORDER, HARMONIC_BALANCE) TIME_MARCHING= DUAL_TIME_STEPPING-2ND_ORDER % % Time Step for dual time steng simulations (s) -- Only used when UNST_CFL_NUMBER = 0.0 % For the DG-FEM solver it is used as a synchronization time when UNST_CFL_NUMBER != 0.0 TIME_STEP= 0.1111 % % Total Physical Time for dual time stepping simulations (s) MAX_TIME= 200.0 % % Unsteady Courant-Friedrichs-Lewy number of the finest grid UNST_CFL_NUMBER= 0.0 % TIME_DISCRE_RADIATION = EULER_IMPLICIT % ------------------------- SCREEN/HISTORY VOLUME OUTPUT --------------------------% % % Screen output fields (use 'SU2_CFD -d <config_file>' to view list of available fields) SCREEN_OUTPUT= (INNER_ITER, RMS_DENSITY, RMS_MOMENTUM-X, RMS_MOMENTUM-Y, RMS_ENERGY) % % History output groups (use 'SU2_CFD -d <config_file>' to view list of available fields) HISTORY_OUTPUT= (ITER, TIME_DOMAIN, RMS_RES) % % User defined functions available on screen and history output. See TestCases/user_defined_functions/. CUSTOM_OUTPUTS= '' % % Volume output fields/groups (use 'SU2_CFD -d <config_file>' to view list of available fields) VOLUME_OUTPUT= (COORDINATES, SOLUTION, PRIMITIVE) % % Writing frequency for screen output SCREEN_WRT_FREQ_INNER= 1 % SCREEN_WRT_FREQ_OUTER= 1 % SCREEN_WRT_FREQ_TIME= 1 % % Writing frequency for history output HISTORY_WRT_FREQ_INNER= 1 % HISTORY_WRT_FREQ_OUTER= 1 % HISTORY_WRT_FREQ_TIME= 1 % % list of writing frequencies corresponding to the list in OUTPUT_FILES OUTPUT_WRT_FREQ= 10, 250, 42 % % Output the performance summary to the console at the end of SU2_CFD WRT_PERFORMANCE= NO % % Overwrite or append iteration number to the restart files when saving WRT_RESTART_OVERWRITE= YES % % Overwrite or append iteration number to the surface files when saving WRT_SURFACE_OVERWRITE= YES % % Overwrite or append iteration number to the volume files when saving WRT_VOLUME_OVERWRITE= YES % % % ------------------------- INPUT/OUTPUT FILE INFORMATION --------------------------% % % Mesh input file MESH_FILENAME= naca0012.su2 % % Mesh input file format (SU2, CGNS) MESH_FORMAT= SU2 % % Mesh output file MESH_OUT_FILENAME= mesh_out.su2 % % Restart flow input file SOLUTION_FILENAME= solution_flow.dat % % Restart adjoint input file SOLUTION_ADJ_FILENAME= solution_adj.dat % % Output tabular file format (TECPLOT, CSV) TABULAR_FORMAT= CSV % % Files to output % Possible formats : (TECPLOT_ASCII, TECPLOT, SURFACE_TECPLOT_ASCII, % SURFACE_TECPLOT, CSV, SURFACE_CSV, PARAVIEW_ASCII, PARAVIEW_LEGACY, SURFACE_PARAVIEW_ASCII, % SURFACE_PARAVIEW_LEGACY, PARAVIEW, SURFACE_PARAVIEW, RESTART_ASCII, RESTART, CGNS, SURFACE_CGNS, STL_ASCII, STL_BINARY) % default : (RESTART, PARAVIEW, SURFACE_PARAVIEW) OUTPUT_FILES= (RESTART, PARAVIEW, SURFACE_PARAVIEW) % % Output file convergence history (w/o extension) CONV_FILENAME= history % % Output file with the forces breakdown BREAKDOWN_FILENAME= forces_breakdown.dat % % Output file restart flow RESTART_FILENAME= restart_flow.dat % % Output file restart adjoint RESTART_ADJ_FILENAME= restart_adj.dat % % Output file flow (w/o extension) variables VOLUME_FILENAME= flow % % Output file adjoint (w/o extension) variables VOLUME_ADJ_FILENAME= adjoint % % Output Objective function VALUE_OBJFUNC_FILENAME= of_eval.dat % % Output objective function gradient (using continuous adjoint) GRAD_OBJFUNC_FILENAME= of_grad.dat % % Output file surface flow coefficient (w/o extension) SURFACE_FILENAME= surface_flow % % Output file surface adjoint coefficient (w/o extension) SURFACE_ADJ_FILENAME= surface_adjoint % % Read binary restart files (YES, NO) READ_BINARY_RESTART= YES % % Reorient elements based on potential negative volumes (YES/NO) REORIENT_ELEMENTS= YES % |
|
February 28, 2024, 15:22 |
|
#2 |
Senior Member
bigfoot
Join Date: Dec 2011
Location: Netherlands
Posts: 676
Rep Power: 21 |
Hi,
Code:
WRITE_VOLUME_OVERWRITE= YES If you change this to NO, paraview files will not be overwritten but instead you will get a paraview file every Code:
OUTPUT_WRT_FREQ= 10, 250, 42 You can use these files to visualize the different timesteps. 250 iterations might not be often enough. |
|
February 28, 2024, 21:20 |
|
#3 |
New Member
Jade Xie
Join Date: Feb 2024
Posts: 2
Rep Power: 0 |
Hi,
Thank you for for response! I revised my code as you mentioned above, but I got error 'Appending iterations to the filename (WRT_VOLUME_OVERWRITE=NO) is incompatible with transient problems.' At first, I thought it's the problem with SOLVER= INC_RANS, but after I tried INC_NAVIER_STOKES, I still got the same error. What might be the problem? I've upload my .cfg, thanks a lot. |
|
February 29, 2024, 03:19 |
|
#4 |
Senior Member
bigfoot
Join Date: Dec 2011
Location: Netherlands
Posts: 676
Rep Power: 21 |
Oops, my mistake. In unsteady simulations, the timestep should already be written to the file every time you save the file.So if you write
Code:
OUTPUT_WRT_FREQ= 100, 1, 100 Do you then get flow_00001.vtu, flow_00002.vtu, etc after every timestep? This should be the case in unsteady flows. |
|
Tags |
output, pitching airfoil, time dependent, unsteady, visualization |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
laplacianFoam - solver settings for steady state diffusion problems | Adam_K | OpenFOAM Running, Solving & CFD | 5 | June 8, 2021 15:33 |
Unexpected deltaT decrease in pimpleFoam simulation | robyTKD | OpenFOAM Running, Solving & CFD | 9 | June 27, 2014 07:52 |
plot over time | fferroni | OpenFOAM Post-Processing | 7 | June 8, 2012 08:56 |
IcoFoam parallel woes | msrinath80 | OpenFOAM Running, Solving & CFD | 9 | July 22, 2007 03:58 |
Could anybody help me see this error and give help | liugx212 | OpenFOAM Running, Solving & CFD | 3 | January 4, 2006 19:07 |