Release Notes OVERFLOW 2.4
2.4d 28 April 2025
Improvements:
o Reworked namelist parsing so that repeating each namelist is no longer
needed. See change 130.
o Modified FOMOCO routines to reduce MPI overhead. NFOMO will no longer be
reset by dynamic cases to output every step. See change 133.
o CFL ramping reworked to accommodate grid sequencing.
o Added a ramping factor for rotor disk simulations
(variable IRAMP in &RDISK_TRIM for MODEL=0).
o Added tri_to_vtk utility to visualize all tri files (including triq_dp)
by converting to VTK format.
Bug fixes:
o Fixed 2D grid adaption bug.
o Fixed missing q.save file when using the SAVE_SPLIT option.
o Fixed potential segfault in INTOUT_WRITS.
o Fixed bug in trip line ITTYP=103 and enabled the simultaneous specification
of ITTYP 102 and 103.
o Fixed indexing error for multi-equation turbulence models when using GLS.
o Disabled temperature dependent gamma due to incorrect treatment.
Removed features:
o DEBUG_SPLIT preprocessor option.
o CONSERVE_MEM.
o Vortex generator source term BC (IBTYP=601).
o Non-point matched fold-over BCs (IBTYP=54, 55, 56).
o Steger-Warming flux-split LHS option (ILHS=1).
o Yee symmetric TVD and AUSM+ schemes (IRHS=2, 3).
o Smoothing options (ISPEC=-1,1,3).
o Baldwin-Barth turbulence model (NQT=100).
o Simple wall distance calculator (WALLDIST=-1,1).
o SCE central scheme (IUPC=0).
2.4c 8 March 2024
Improvements:
o Improved force and moment calculation with non-zero turbulent eddy viscosity at boundaries, to accommodate surface roughness models.
o Enabled Knopp wall roughness model for Wilcox k-omega turbulence model.
o Allowed force and moment calculation on data surface grids.
o Added DEBUG=6 output to q.trans file for AFT and LM transition models.
o Added skin friction coefficient as variable 6 in DEBUG=1 q.turb file.
o Created split-grid restart capability (when using the same number of MPI groups) by using SAVE_SPLIT=.T. in namelist GLOBAL. This avoids residual spikes on restart due to auto-split grid boundaries.
o Changed from Spalart Allmaras SA-neg-noft2 to SA-neg (i.e., turned on the ft2 term). This affects NQT=102,104,302, and should improve initial convergence with little effect on final result.
o Added a CFL ramping restart capability.
o Added option to write (non-standard) double-precision triq files for improved acoustic calculations.
o Moved source files into modules to support GPU porting and updated Makefile system.
Rotorcraft interface:
o Fixed rotor interface force and moment reporting to use (each) rotor center as the moment reference center for inertial, rotor, and blade reference frames. Included PSIOFFSET in the rotor frame definition. (From David Palmer, Archer Aviation.)
o Added additional rotor force and moment history output files as .rotorframe.txt and .bladeframe.txt to simplify post-processing. Output is in the form “blade# step# forces(1:3) moments(1:3)” (From David Palmer, Archer Aviation.)
Bugs fixed:
o Fixed DDES “tiny” numbers to better represent cell size/face area tiny values. This can improve convergence. Tiny numbers adjusted elsewhere as well.
o Fixed bugs when using Global Linear Solve with NQT=301 (SA model with AFT transition model).
o Improved DCF overset grid donor search near sharp trailing edges.
o Adjusted IBLANK usage in central differencing stencils near physical and hole boundaries.
o Allowed use of BC#201 for data surface grids.
2.4b 4 November 2022
CMake:
o Moved “makeall_cmake” script to “cmakeall”.
o Can specify FLAG_VERSION in cmakeall for different Intel compiler flags.
o Improved cmakeall and cmake.sys documentation.
Compiler options:
o Added -mcmodel=medium to gfortran flags to avoid link problems.
o Added -fp-model precise to Intel compiler flags to avoid 32-bit floating-point errors.
o Updated Makelist (for compile testing).
Improvements:
o Added DELTA>1 eigenvalue smoothing for HLLE++ left-hand side per suggestion of Josh Elliott (Boeing).
o Added dual number utility code and example.
o Allowed ILHS=16/17/18 when using multigrid.
o Improved use of “fold-over” grid subsets for DCF inverse maps to eliminate orphans.
o Updated test/airfoil_drop_2d Mach=0.2 case to use NQT=104.
Bugs fixed:
o Fixed misplaced parentheses in NQT=103/104 SA diffusion terms.
o Fixed some XML parser problems (for reading Config.xml, Scenario.xml).
o Fixed rotor disk model problem with multigrid and grid splitting.
o Fixed two problems with SPLITM subsets and XINTOUT mixed mode, causing extra fringe points from Pegasus5 to be shown in both PLOT3D-style output files and TRIQ files.
2.4a 8 September 2022
Bugs fixed:
o Fixed GLOBAL_LINEAR_SOLVE to work with data surface grids.
o Fixed some $OMP PRIVATE lists, especially for low-Mach preconditioning and/or C-grids.
o Fixed default GMP Constraint Center to be the body center of mass (CM).
o Improved some srb_axi and airfoil_drop_2d test case input files.
o Removed Baldwin-Barth from flat_plate* test cases.
2.4 20 July 2022
The big change for OVERFLOW Version 2.4 is the addition of the Parallel SSOR
linear solver, essentially a restricted additive Schwarz option, as described
in AIAA Paper # 2021-2747, which can be accessed through a new logical
variable GLOBAL_LINEAR_SOLVE in &GLOBAL. The basic idea is to make the
overset boundaries more implicit by communicating delta Q updates during the
linear solver process, doing a better job of solving the linear problem and
hopefully resulting in better nonlinear convergence, at the cost of more MPI
communication and higher memory usage. More details on how to use this new
option are discussed in doc/readmes/README.gls. The 2.4 release should
maintain backwards compatibility for most users, with slight changes to wall
functions, but for those users who have their own modifications to the code,
be aware that there were fairly significant changes at the OVERFL level to
accommodate inter-processor communication during the linear solve.
Further notes:
1: For startup and for higher speed flows, the new GLOBAL_LINEAR_SOLVE option
can be used in conjunction with ILHS = 6/7/8 to use more dissipative
Steger-Warming flux Jacobians as well as implicit physical and overset
boundary conditions; this should help users who have had difficulty getting
the Improved SSOR (ILHS = 26/27/28) to start. Be aware that starting
simulations using a constant DTPHYS may result in very high CFL numbers near
the wall and might destabilize the simulation.
2: The Improved and Parallel SSOR schemes have been adjusted to allow for
two-point boundary condition linearizations, improving convergence when using
symmetry planes.