Hi All,
I'm using the FSI module (Fluid Structure Interaction) to simulate a sphere in a fluid. My problem is that I want to simulate an elastic sphere rotating in a viscous fluid. So, I've created my geometry as a region of fluid, with a solid sphere in it. Normally I would use the 'sliding wall' boundary condition on the sphere, but it is causing the 'fluid solid interface' node to be overwritten. This is an undesirable effect, so I went in to the equation view and modified the VxWall, VyWall, and VzWall nodes to reflect their individual components of velocity at the surface of the sphere. I defined a spherical coordinate system for the sphere and solved for the individual components of velocity as:
u = -sysU.r*Omega*sin(sysU.phi*sin(sysU.theta)
v = sysU.r*Omega*sin(sysU.theta)*cos(sysU.phi)
w = 0
I had success with a similar approach when I was only trying to get the elastic sphere to translate about one axis (z), but cannot seem to get this current simulation to work. The error that I am getting is "Undefined value found in the equation residual vector. There are 476 degrees of freedom giving NaN/Inf in the vector for the variable x". Is there a better way to do this, or do you have any suggestions?
Thank you,
Adam
I'm using the FSI module (Fluid Structure Interaction) to simulate a sphere in a fluid. My problem is that I want to simulate an elastic sphere rotating in a viscous fluid. So, I've created my geometry as a region of fluid, with a solid sphere in it. Normally I would use the 'sliding wall' boundary condition on the sphere, but it is causing the 'fluid solid interface' node to be overwritten. This is an undesirable effect, so I went in to the equation view and modified the VxWall, VyWall, and VzWall nodes to reflect their individual components of velocity at the surface of the sphere. I defined a spherical coordinate system for the sphere and solved for the individual components of velocity as:
u = -sysU.r*Omega*sin(sysU.phi*sin(sysU.theta)
v = sysU.r*Omega*sin(sysU.theta)*cos(sysU.phi)
w = 0
I had success with a similar approach when I was only trying to get the elastic sphere to translate about one axis (z), but cannot seem to get this current simulation to work. The error that I am getting is "Undefined value found in the equation residual vector. There are 476 degrees of freedom giving NaN/Inf in the vector for the variable x". Is there a better way to do this, or do you have any suggestions?
Thank you,
Adam