GPUSPH / gpusph

The world's first CUDA implementation of Weakly-Compressible Smoothed Particle Hydrodynamics
164 stars 67 forks source link

about floating and boundary with addXYZFile #65

Open Dongxueyang opened 4 years ago

Dongxueyang commented 4 years ago

Hi, everyone: @Oblomov @Narcolessico @Azrael3000 I am making a simulation of a sphere (floating) fall into the water.

I generate the the boundary like this: GeometryID boundary = addXYZFile(GT_FIXED_BOUNDARY, FT_BORDER, “boundary.dat”,"boudnary.stl");

And generate the fluid like this: GeometryID fluid= addXYZFile(GT_FLUID, FT_SOLID, “fluid.dat”);

And generate the floating sphere like this: GeometryID floating_sphere= addSphere(GT_FLOATING_BODY, FT_BORDER, Point(50,50,30),5); setMassByDensity(floating_sphere,physparam->rho0[0]2.0); setParticleMassByDensity(floating_sphere,physparam->rho0[0]2.0);

After that I run make successfully. But when I run ./GPUSPH , the results is wrong. Only when time=0 I can see the sphere(floating). Next timestep, the sphere disapeared. I am in trouble with the problem, and could you give some advices? Maybe there is some wrong with the model. Hope for your reply. Dong

Dongxueyang commented 4 years ago

@Oblomov Do you see the message and the problem mentioned above?

Oblomov commented 4 years ago

Hello @Dongxueyang,

Is the sphere's initial placement in water? Particles and bodies disappearing in such a short time is usually a side-effect of initialization issues such as particles not being correctly deleted or overlapping in some way. Can you verify what the situation is like concerning the sphere placement? Does it happen if the sphere is placed outside of the water at the beginning?

(Also please be patient, we reply when we have the time.)