search

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

update_nonInv_Beams #18

Open cubayang opened 2 years ago

cubayang commented 2 years ago

Hi Nick,

I tried to simulate the swimmer case with pyIB2d, and I found that a code snippet in IBM_Driver.py:

if update_nonInv_Beams_Flag and nonInv_beams_Yes: from update_nonInv_Beams import update_nonInv_Beams

This function is application specific, located with main2d

    nonInv_beams_info = update_nonInv_Beams(dt,current_time,beams_info)

Should I use "nonInv_beams_info" to replace "beams_info" in "update_nonInv_Beams" function since it is the nonInv Beam?

Thanks in advance, Yang

nickabattista commented 2 years ago

Hi Yang,

Yes, I believe that should be changed, since the beam_Info corresponds to the torsional spring model for bending. Thank you for catching that. How is the implementation going otherwise?

cubayang commented 2 years ago

Hi Nick,

The implementation goes well except for the speed. The setup of the case uses a uniform Cartesian mesh and one processor, I am wondering if I can try with non-uniform Cartesian meshes and multiple processors to increase the computational speed?Thanks!