[ycui]

Dear Foamer:

The integration schemes in the Lagrangian frame only contain Euler & Analytical, so I'm trying to implement the Runge Kutta method inside.

The member functions in IntegrationScheme.H is:

//- Perform the Integration
virtual integrationResult integrate
(
const Type& phi,
const scalar dt,
const Type& alphaBeta,
const scalar beta
) const;

And from KinematicParcel.C it is used in this way:

IntegrationScheme<vector>::integrationResult Ures =
td.cloud().UIntegrator().integrate(U_, dt, abp, bp);

with:

const vector abp = (Feff.Sp()*Uc_ + (Feff.Su() + Su))/massEff;
const scalar bp = Feff.Sp()/massEff;

My question is: what's the meaning of alphaBeta & beta in integrate()?

Best regards,

ycui

[ycui]

I find the answer: there are user defined type, different integration scheme has a different types.

Originally Posted by ycui

Dear Foamer:

The integration schemes in the Lagrangian frame only contain Euler & Analytical, so I'm trying to implement the Runge Kutta method inside.

The member functions in IntegrationScheme.H is:


And from KinematicParcel.C it is used in this way:

with:


My question is: what's the meaning of alphaBeta & beta in integrate()?

Best regards,

ycui

[cdunn6754]

They are just constants for your integration. You can see how they are defined, at least for your case, about 4 lines above the call to integrate in KinematicParcel.C

[rqli_npu]

hello ,did you successfully use RK nethod , i also wang to do this thing, can you give me some help