CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > SU2

Problems with Time Dependent Simulation (Unsteady)

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 28, 2024, 13:44
Question Problems with Time Dependent Simulation (Unsteady)
  #1
New Member
 
Jade Xie
Join Date: Feb 2024
Posts: 2
Rep Power: 0
JadeX is on a distinguished road
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
%
JadeX is offline   Reply With Quote

Old   February 28, 2024, 15:22
Default
  #2
Senior Member
 
bigfoot
Join Date: Dec 2011
Location: Netherlands
Posts: 676
Rep Power: 21
bigfootedrockmidget is on a distinguished road
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
250 iterations, named flow_xxxx.vtu


You can use these files to visualize the different timesteps. 250 iterations might not be often enough.
bigfootedrockmidget is offline   Reply With Quote

Old   February 28, 2024, 21:20
Default
  #3
New Member
 
Jade Xie
Join Date: Feb 2024
Posts: 2
Rep Power: 0
JadeX is on a distinguished road
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.
Attached Files
File Type: txt 9my_naca0012.txt (24.1 KB, 0 views)
JadeX is offline   Reply With Quote

Old   February 29, 2024, 03:19
Default
  #4
Senior Member
 
bigfoot
Join Date: Dec 2011
Location: Netherlands
Posts: 676
Rep Power: 21
bigfootedrockmidget is on a distinguished road
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.
bigfootedrockmidget is offline   Reply With Quote

Reply

Tags
output, pitching airfoil, time dependent, unsteady, visualization


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 Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


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


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