[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[ESPResSo-users] constant velocity field in "object in fluid"

From: gianmarco
Subject: [ESPResSo-users] constant velocity field in "object in fluid"
Date: Mon, 22 Jul 2019 17:26:54 +0200
User-agent: Mozilla/5.0 (X11; Linux x86_64; rv:60.0) Gecko/20100101 Thunderbird/60.7.2

Hello everyone,

I'm writing this messsage because I'm trying to reproduce a fluid constant velocity field within the package "object in fluid". To do this, I created a fluid at rest, therefore by setting the parameter "ext_force_density" to a zero 3d vector. To impose a constant velocity field I inserted a wall on the top of the simulation  box, with the command


wall_shape = shapes.Wall(normal=[1, 0, 0], dist=Lz)

walls[0]=lbboundaries.LBBoundary(shape=wall_shape, velocity=np.array([vcell,0.0,0.0]))


where vcell is a parameter representing the boundary velocity. The resulting effect is that a mesh inserted in the fluid "feels" the effect of the velocity field, i. e. it moves at constant velocity. However, the configuration files produced with the command

lbf.print_vtk_velocity("fluid_%04i.vtk" %i)

show a uniform and constant velocity field, equal to zero. Does anyone know how to fix this problem?

Thank you in advance for you response.

Best regards


reply via email to

[Prev in Thread] Current Thread [Next in Thread]