How to set up the paddle motions for irregular waves?

[SurelyD]

Hi everyone,

A quick question about how to set the boundary conditions. I modified the code from WaveTank.cu to simulate waves in a basin and it worked well with regular waves. Now that I want to modify the file to generate irregular waves, I find a there is a problem I can not fix.

I figure out the way to calculate the time series of the paddle stroke s(t) following previous studies and the webpage of DualSPHysics (https://github.com/DualSPHysics/DualSPHysics/wiki/3.-SPH-formulation#39-wave-generation) and write the time series into a file. The problem is that I do not know how to read the data as paddle motions in the code... If you know how to do it, please feel free to drop a reply.

PS. the potentially related lines should be as below:

%%%
// The stroke value is given at free surface level H
float stroke = 0.2;
// m_mbamplitude is the maximal angular value for paddle angle
// Paddle angle is in [-m_mbamplitude, m_mbamplitude]
paddle_amplitude = atan(stroke/(2.0*(H - paddle_origin.z)));
cout << "\npaddle_amplitude (radians): " << paddle_amplitude << "\n";
paddle_omega = 2.0M_PI/0.8; // period T = 0.8 s
%%%
void
WaveTank::moving_bodies_callback(const uint index, Object object, const double t0, const double t1,
const float3& force, const float3& torque, const KinematicData& initial_kdata,
KinematicData& kdata, double3& dx, EulerParameters& dr)
{
dx= make_double3(0.0);
kdata.lvel=make_double3(0.0f, 0.0f, 0.0f);
if (t1> paddle_tstart && t1 < paddle_tend){
kdata.avel = make_double3(0.0, paddle_amplitudepaddle_omegasin(paddle_omega*(t1-paddle_tstart)),0.0);
EulerParameters dqdt = 0.5*EulerParameters(kdata.avel)*kdata.orientation;
dr = EulerParameters::Identity() + (t1-t0)dqdtkdata.orientation.Inverse();
dr.Normalize();
kdata.orientation = kdata.orientation + (t1 - t0)*dqdt;
kdata.orientation.Normalize();
}
else {
kdata.avel = make_double3(0.0,0.0,0.0);
kdata.orientation = kdata.orientation;
dr.Identity();
}
}

All comments and suggestions welcomed!

Thanks!