Embed_force for 3D

[stephen muwanguzi]

I am working on wave and cylinder structure interaction program in 3D. 

It seems that current embed_force function is built on assertion that dimension=2. 

Kindly wanted to ask if any body  has some advice on how to extend embed_force to calculate Pressure force, Fp and Viscous force Fmu for 3D. 

Thanks 
Stephen 

[Francesco Picella]

Hi,

Have a look at Arthur's implementation of embed_force
https://basilisk.fr/sandbox/ghigo/src/myembed.h

It is functional in 3D as well.

Besides, watch out on computing forces and torques on _moving embedded boundaries_.

Computing torques, there is an issue whether the body has a non zero angular velocity:

https://basilisk.fr/sandbox/fpicella/test/couette_torque.c

To avoid this issue, a more complex definition of embed_torque function is required.

At the moment, it is available only in 2D:

https://basilisk.fr/sandbox/fpicella/src/compute_embed_color_force_torque_RBM.h

Hope it helps

Francesco 

[kevin.a....@gmail.com]

I'm computing forces on 3D ships using embed.  I use an interpolation method that searches for cells near the ship panel of interest.  To reduce searching time, I ensure that the mesh is fine near the ship for all time steps, and I find a nearest cell location for each panel centroid at the beginning of the simulation.  I'm satisfied with accuracy and computational speed.  

Kevin