[abbas]

Thanks Adrian, I tried other 3d model like 3d cylinder, it runs without problem. the 2d falling droplet also works without problem, but it gives error for 3d falling droplet. I did not change anything in the code I used the Ubuntu windows shell installed and run it.

[abbas]

Dear Adrian,
Thanks again for your reply, Can you introduce a paper that use the same method which is used in the free surface simulation in open LB. I would like to simulate droplet impact on the solid surface, I believe the interaction between the solid and liquid is not considered in the free surface model.
What is the solution for it can I combine the free energy model and free surface model for droplet impact on the solid surface with different hydrophobicity?

one more question when I tried to run the 3D falling droplet i got the following error, can you help?
Segmentation fault (core dumped)\

Thanks
Abbas

[abbas]

Hi all,
Can you explain how this command change only the flow cells into the material number 1 and kepps the boundary cells as 2?

superGeometry.rename( 2,1,{1,1} );

I found this function
void rename(int fromM, int toM);
/// replace one material that fulfills an indicator functor condition with another

so the materials are change from 2 to 1, then what is {1,1};

Thanks

[abbas]

ok thanks, got it.
can you explain how the following code generate a circle?
I think only rectangular and circle and triangle can be defined using the library of Openlb and if I want to create for example a hexagonal I have to write my own code.

template <typename S>
class IndicatorCircle2D : public IndicatorF2D<S> {
private:
const Vector<S,2> _center;
const S _radius;
const S _radius2;
public:
IndicatorCircle2D(Vector<S,2> center, S radius);
Vector<S,2> const& getCenter() const;
S const getRadius() const;
bool operator() (bool output[], const S input[]) override;
bool distance(S& distance, const Vector<S,2>& origin, const Vector<S,2>& direction, int iC=-1) override;
bool normal(Vector<S,2>& normal, const Vector<S,2>& origin, const Vector<S,2>& direction, int iC=-1) override;
S signedDistance(const Vector<S,2>& input) override;
using IndicatorF2D<S>::distance;
};

[abbas]

or for example if I want to put a hexagonal obstacle what part of the code should be changed?

[abbas]

Thanks Adrian, The thing I don’t understand is there is no sign of the defining the geometry for example the circular obstacle. I only found the center of the obstacle information in cylinder2d.cpp file.
In the userguide it is mentioned that the geometry is read through an .stl file which I also cannot find it in the directory of the cylinder2D example. for instance if I wanted to change the circular obstacle into a rectangular one what should I do?

[Author]

Posts