I would like to modify damBreak3d.cpp in such a way that there is a constant inlet fluid.
I found into freeSurfaceBoundaryCondition3D.h the following code:
template< typename T,template class Descriptor>
class PouringLiquid3D : public BoxProcessingFunctional3D {
public:
PouringLiquid3D(Dynamics<T,Descriptor>* dynamicsTemplate_, Array<T,3> injectionVelocity_)
: dynamicsTemplate(dynamicsTemplate_), injectionVelocity(injectionVelocity_)
{ }
PouringLiquid3D(PouringLiquid3D<T,Descriptor> const& rhs)
: dynamicsTemplate(rhs.dynamicsTemplate->clone()),
injectionVelocity(rhs.injectionVelocity)
{ }
PouringLiquid3D<T,Descriptor>* operator=(PouringLiquid3D<T,Descriptor> const& rhs)
{
PouringLiquid3D<T,Descriptor>(rhs).swap(this);
return this;
}
void swap(PouringLiquid3D<T,Descriptor>& rhs) {
std::swap(dynamicsTemplate, rhs.dynamicsTemplate);
std::swap(injectionVelocity, rhs.injectionVelocity);
}
virtual ~PouringLiquid3D() {
delete dynamicsTemplate;
}
virtual void processGenericBlocks(Box3D domain, std::vector<AtomicBlock3D> atomicBlocks);
virtual PouringLiquid3D<T,Descriptor> clone() const;
virtual void getTypeOfModification(std::vectormodif::ModifT& modified) const {
modified[0] = modif::dataStructure; // Fluid
modified[1] = modif::staticVariables; // rhoBar.
modified[2] = modif::staticVariables; // j.
modified[3] = modif::staticVariables; // Mass
modified[4] = modif::staticVariables; // Mass-fraction
modified[5] = modif::staticVariables; // Flag-status
modified[6] = modif::nothing; // Normal
modified[7] = modif::nothing; // Interface lists
modified[8] = modif::nothing; // Curvature.
modified[9] = modif::nothing; // Outside density.
}
private:
Dynamics<T,Descriptor>* dynamicsTemplate;
Array<T,3> injectionVelocity;
};
…but i don’t know how to define:
std::vector<AtomicBlock3D*> atomicBlocks
since that in the original code are defined only MultiBlockLattice3D!
Every idea is welcome!!!
Best Regards!!