147 lines
3.2 KiB
C++
147 lines
3.2 KiB
C++
/*------------------------------- phasicFlow ---------------------------------
|
|
O C enter of
|
|
O O E ngineering and
|
|
O O M ultiscale modeling of
|
|
OOOOOOO F luid flow
|
|
------------------------------------------------------------------------------
|
|
Copyright (C): www.cemf.ir
|
|
email: hamid.r.norouzi AT gmail.com
|
|
------------------------------------------------------------------------------
|
|
Licence:
|
|
This file is part of phasicFlow code. It is a free software for simulating
|
|
granular and multiphase flows. You can redistribute it and/or modify it under
|
|
the terms of GNU General Public License v3 or any other later versions.
|
|
|
|
phasicFlow is distributed to help others in their research in the field of
|
|
granular and multiphase flows, but WITHOUT ANY WARRANTY; without even the
|
|
implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
|
|
|
|
-----------------------------------------------------------------------------*/
|
|
|
|
#ifndef __AdamsBashforth5_hpp__
|
|
#define __AdamsBashforth5_hpp__
|
|
|
|
|
|
#include "AdamsBashforth4.hpp"
|
|
#include "pointFields.hpp"
|
|
|
|
namespace pFlow
|
|
{
|
|
|
|
|
|
/**
|
|
* Fifth order Adams-Bashforth integration method for solving ODE
|
|
*
|
|
* This is a one-step integration method and does not have prediction step.
|
|
*/
|
|
class AdamsBashforth5
|
|
:
|
|
public AdamsBashforth4
|
|
{
|
|
private:
|
|
|
|
friend class processorAB4BoundaryIntegration;
|
|
|
|
realx3PointField_D dy4_;
|
|
|
|
protected:
|
|
|
|
const auto& dy4()const
|
|
{
|
|
return dy4_;
|
|
}
|
|
|
|
auto& dy4()
|
|
{
|
|
return dy4_;
|
|
}
|
|
|
|
public:
|
|
|
|
/// Type info
|
|
TypeInfo("AdamsBashforth5");
|
|
|
|
// - Constructors
|
|
|
|
/// Construct from components
|
|
AdamsBashforth5(
|
|
const word& baseName,
|
|
pointStructure& pStruct,
|
|
const word& method,
|
|
const realx3Field_D& initialValField,
|
|
bool keepHistory);
|
|
|
|
|
|
|
|
/// Destructor
|
|
~AdamsBashforth5() override =default;
|
|
|
|
/// Add a this to the virtual constructor table
|
|
add_vCtor(
|
|
integration,
|
|
AdamsBashforth5,
|
|
word);
|
|
|
|
|
|
// - Methods
|
|
|
|
void updateBoundariesSlaveToMasterIfRequested()override;
|
|
|
|
/// return integration method
|
|
word method()const override
|
|
{
|
|
return "AdamsBashforth5";
|
|
}
|
|
|
|
bool correct(
|
|
real dt,
|
|
realx3PointField_D& y,
|
|
realx3PointField_D& dy,
|
|
real damping = 1.0) override;
|
|
|
|
bool correctPStruct(
|
|
real dt,
|
|
pointStructure& pStruct,
|
|
realx3PointField_D& vel) override;
|
|
|
|
};
|
|
|
|
|
|
/*template<typename activeFunctor>
|
|
bool pFlow::AdamsBashforth5::intRange(
|
|
real dt,
|
|
realx3Vector_D& y,
|
|
realx3Vector_D& dy,
|
|
activeFunctor activeP )
|
|
{
|
|
auto d_dy = dy.deviceViewAll();
|
|
auto d_y = y.deviceViewAll();
|
|
auto d_history = history_.deviceViewAll();
|
|
auto activeRng = activeP.activeRange();
|
|
|
|
Kokkos::parallel_for(
|
|
"AdamsBashforth5::correct",
|
|
rpIntegration (activeRng.first, activeRng.second),
|
|
LAMBDA_HD(int32 i){
|
|
if( activeP(i))
|
|
{
|
|
|
|
d_y[i] += dt*(
|
|
static_cast<real>(1901.0 / 720.0) * d_dy[i]
|
|
- static_cast<real>(2774.0 / 720.0) * d_history[i].dy1_
|
|
+ static_cast<real>(2616.0 / 720.0) * d_history[i].dy2_
|
|
- static_cast<real>(1274.0 / 720.0) * d_history[i].dy3_
|
|
+ static_cast<real>( 251.0 / 720.0) * d_history[i].dy4_
|
|
);
|
|
d_history[i] = {d_dy[i] ,d_history[i].dy1_, d_history[i].dy2_, d_history[i].dy3_};
|
|
}
|
|
});
|
|
Kokkos::fence();
|
|
|
|
return true;
|
|
}*/
|
|
|
|
} // pFlow
|
|
|
|
#endif //
|