Version 1.7: Example Applications
There is no mesh advection, except the interface that is moved in a Lagrangian fashion. For the structure a classic Lagrangian point-of-view is adopted, for the fluid the typical Eulerian description.
The SURC solvers are global-timestep, consistent mass-matrix CFD solvers whose formulation is based on the Charactersitic-Based split method, which is mostly the work of O.C. Zienkievicz and P. Nithiarasu.
The SURC solvers are not yet adapted to problems with inhomogenous pressure boundary conditions. The work is still heavily in progress, so please do not consider this production code.
The local time-step/compressibility experiments failed, likely due to the stronger coupling between the elements with higher spatial accuracy. The local time-stepping that is often used with CBS and linear interpolation, is hence no option.
The SURC4 and SURC6 solvers use an implicit pressure solution, and both converge slowly, but very steadily after a short start-up phase. The SURC5 solver uses pseudo-compressibility and is due to the absence of the pressure split matrix or a simultaneous pressure solution, apt to be the most performant choice for FSI. The down-side of this method is less steady convergence. The most recent version, SURC 7, adds adaptive timesteps to the SURC 6 scheme.
All methods produce the same velocity profiles.
Please note that higher velocities require higher mesh resolutions for stability. E.g. for 1000 a resolution of 8 or higher is recommendable.
Flow-direction-indicating particles released from the driving wall. Twice rendered without color and below with coloring and blending turned on.
It is configured for use with a version of Animath that uses Sparse Lib/IML++ as iterative matrix solver.
The approach is that both compressibility and the size of the timesteps are made dependent on the relative difference. Such that great changes of velocity between two consecutive timesteps, increase medium compressibility and decrease the size of the next step. The results are comparative to reference data of Chorin already at low mesh resolutions.