[achodankar]

Hello Mathias,
Thank you. I understand it well now.

Yours sincerely,

Abhijeet C.

[achodankar]

Hello Mathias,
This explanation is very helpful.

Thank you.

Yours sincerely,

Abhijeet C.

[achodankar]

Hello Mathias,
I will try to check if the compiler installation on the HPC is correct. Thank you for your prompt response.

Thank you.

Yours sincerely,

Abhijeet C.

[achodankar]

Hello Adrian,
Thank you so much for your help. I stored the scalar field of velocity separately using calculator tool in ParaView, and now I can plot the contours of velocity field. The problem is resolved.

Thank you so much!!! I really appreciate it.

Thank you.

Yours sincerely,

Abhijeet C.

[achodankar]

Hello Adrian,
The asymmetry issue persists even after increasing the resolution. I have ignored the issue temporarily, since for higher resolution it not that much difference. However, from the fundamental point of view, it is incorrect. So, I am pondering if anything can be tweaked in the stl reader.

Here is the link to the image of the geometry: https://ibb.co/zGbNp3s

Thank you.

Yours sincerely,

Abhijeet C.

[achodankar]

Hello Adrian,
Thank you for your response. My apologies for the late reply.

What is the use of the method parameter?
Choose indication method 0: fast, less stable 1: slow, more stable (for untight STLs)

Thank you.

Yours sincerely,

Abhijeet C

[achodankar]

Hello Adrian,
The arrangement of the nodes on the top, bottom, left and right side is not symmetric. Basically, it should be mirror images(Eg:-circular geometry). That’s not the case presently. I didn’t have that problem, when the geometries were created using indicators.

Thank you.

Yours sincerely,

Abhijeet C.

[achodankar]

Hello Adrian,
I understand the usage of N in general examples wherein the geometry is created using indicators. The way we define N will determine the deltax value. However, I fail to understand how it will work when the geometry is imported using an STL file. The imported geometry has predetermined deltax value based on the triangulations. SolidWorks or FreeCAD will just spit out the STL files. How is Openlb communicating with these STL files and setting or resetting the deltax values based on a given N value?

Thank you.

Yours sincerely,

Abhijeet C.

[achodankar]

Hello Developers,
The problem is seen when I import a 3d geometry consisting of curves. I tried with a circular and semicircular geometry, thereby facing the same issues for these geometries. Is there any particular part of openlb I should address to fix this issue? I would really appreciate your feedback.

Thank you.

Yours sincerely,

Abhijeet C.

[achodankar]

Hello Adrian,
Thank you very much for your prompt reply. Your suggestion did work out. I can import the stl file successfully in openlb. However, the problem is that the nodes at the wall are not symmetric at all. The top and bottom nodes are not symmetric. It applies the same to the left and right nodes. I have imported a simple cylinder geometry. I would be grateful if you could suggest me how do I solve this issue.

Thank you.

Yours sincerely,

Abhijeet C.

[achodankar]

Yes the standard example Poiseuille3d works. What can be the possible cause for the error? I have changed the dimensions of cylinder.

[achodankar]

Hello Julius,
I understand it well now.

Thank you.

Yours sincerely,

Abhijeet C.

[achodankar]

Hello Julius,
I followed your suggestion and the issue is fixed now.

Thank you very much!!! You have been of great help.

Yours sincerely,

Abhijeet C.

[achodankar]

Hello Julius,
I tried this out, but I obtain zero values for velocity in all files. Also, the points passed out as input are lattice coordinates, Correct?

{

int Ly = converter.getLatticeLength( ly );
int Lx = converter.getLatticeLength( lx );

for (int iX=0; iX<Lx; iX++)
{

ofstream mfile1;

mfile1.open(“v_nointerpolate_lat_”+ std::to_string(iX)+ “_”+ std::to_string(iT) + “.csv”);
mfile1 <<“x,” << “y_lattice,” << “U_lattice,”<< “y_phys,”<<“U_phys,”<<std::endl;

T num_phys[Ly],num_lat[Ly];
int y[Ly],x[1];

for ( int iY=0; iY<=Ly; ++iY )
{
int point[2];
point[0] = iX;
point[1] = iY;

//Store Lattice Velocity (Numerical)
T numerical_lat[2];

SuperLatticeVelocity2D<T, DESCRIPTOR> velocity_lat( sLattice);

velocity_lat( numerical_lat,point );

//Store Physical Velocity (Numerical and Analytical)
T numerical_phys[2];

SuperLatticePhysVelocity2D<T, DESCRIPTOR> velocity_phys( sLattice, converter );

velocity_phys( numerical_phys,point );

//storing data

mfile1 << point[0] <<“,”<< point[1]<<“,” <<numerical_lat[0]<<“,”<<point[1]*converter.getConversionFactorLength()<<“,”<<numerical_phys[0]<<endl;

}
mfile1.close();

}

[achodankar]

Hello Julius,
Thank you very much for your help. I will try it out.

Thank you.

Yours sincerely,

Abhijeet C.

[Author]

Posts