[Top][All Lists]

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

Re: [ESPResSo-users] Lattice Boltzmann - Force Exerted on Object by Flui

From: Kai Szuttor
Subject: Re: [ESPResSo-users] Lattice Boltzmann - Force Exerted on Object by Fluid
Date: Mon, 25 Jun 2018 16:16:03 +0200
User-agent: NeoMutt/20170714-126-deb55f (1.8.3)

Hi Esteban,

1. I can not give a exper answer here because I'm not familiar with the
   OIF implementation, but in principal you can get
   the total force on the object by summing up the forces you get with
   F_particle = -fric * (v_particle - v_fluid)
   where v_particle is the particle velocit and v_fluid is the fluid
   velocity at the position of the particle.
2. I can not reproduce your issue, I tested with the following minimal

import espressomd
system = espressomd.System(box_l = [10.0, 10.0, 10.0])
system.time_step = 0.01 = 0.4
lb_fluid =, dens=1.0, visc=1.0, fric=1.0, 

I hope that helps.


print lb_fluid.get_interpolated_velocity([1.0, 2.3, 5.2])

On Mon, Jun 25, 2018 at 12:54:57PM +0000, Ongini  Esteban wrote:
>    Dear Espresso Developers,
>    I am currently implementing  a simulation using the LB in espresso along
>    with the OIF extension . I have two questions regarding the Lattice
>    Boltzmann code:
>     1. I was wondering if it is possible to output the force exerted on an
>        object by the fluid using the code ?
>     2. The second question is regarding how use
>        get_interpolated_velocity(pos) function? I have tried to implement
>        this however, I receive the following error message:  
>                '' object has no attribute
>    'get_interpolated_velocity'   
>    Thank You,
>    Esteban

Attachment: signature.asc
Description: PGP signature

reply via email to

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