Turek-hron case Eulerian grid

nickabattista / IB2d

An easy to use immersed boundary method in 2D, with full implementations in MATLAB and Python that contains over 75 built-in examples, including multiple options for fiber-structure models and advection-diffusion, Boussinesq approximations, and/or artificial forcing.

GNU General Public License v3.0

164 stars 90 forks source link

Turek-hron case Eulerian grid #11

Open yanggthomas opened 5 years ago

yanggthomas commented 5 years ago

Hi, the Eulerian grid description in channel_channel.m file is different from that in input2d

nickabattista commented 5 years ago

Thanks - if you noticed the resolution remains the same Lx/Nx = 1.0/512 (input2d) = 2.0/1024 (channel_channel.m) and similar for the Ly/Ny = 0.25/128 (input2d) = 1.0/512 (channel_channel.m). Since the Channel_Channel.m file is creating all the necessary input files (.vertex, .spring, .target, etc.), the fact that the grid descriptions are different don't matter as long as the resolutions are consistent. It probably should be remade to avoid confusion; thanks for pointing that out!