espressomd / espresso

The ESPResSo package
https://espressomd.org
GNU General Public License v3.0
229 stars 183 forks source link

Split LB tutorial into LB, polymer and Langevin simulation tutorials #3939

Closed jngrad closed 3 years ago

jngrad commented 4 years ago

The LB tutorial introduces the Flory theory of polymers, the Rouse regime (implicit solvent) and the Zimm regime (LB fluid). The only difference between Langevin Dynamics and lattice-Boltzmann is the absence resp. presence of the second term in the Kirkwood-Zimm equation of the diffusion coefficient of polymers:

 D=(D_0)/(N) + \frac{k_B T}{6 \pi \eta } \left\langle (1)/(R_h) \right\rangle

We could move part 2 and 3 to a self-contained polymer tutorial. The part 3 can be switched from LD to LB by changing a variable at the beginning of the file (both versions are tested in CI).

jngrad commented 4 years ago

The resulting LB tutorial could get an additional part where we simulate particles falling in a vacuum vs. in a LB fluid, to reproduce the video from the "ESPResSo features" lecture.

schlaicha commented 3 years ago

I've split the tutorial in #4052 and will briefly check the physics with Christian later. For the suggested LB tutorial is the script somewhere around? Is this at all realistic to have it interactively in a notebook?

schlaicha commented 3 years ago

Offline discussion with Christian:

schlaicha commented 3 years ago

I cannot change the issue title, can somebody either give me the rights (assign?) or change it correspondingly?

jngrad commented 3 years ago

@schlaicha You should now be able to change the title.

I don't know if there is a script for the sedimentation visualization (the video was probably made in the Tcl days).

schlaicha commented 3 years ago

@schlaicha You should now be able to change the title.

No, also I cannot add myself as assignee.

jngrad commented 3 years ago

ok, maybe GitHub prevents users from changing other people's ticket titles... I can take care of that, what title do you propose?

Unfortunately, we cannot give you more rights, because the next permission level is Triage, which would allow you to merge approved PRs by editing the labels.

schlaicha commented 3 years ago

I would just change according to the PR: "Split LB tutorial into LB, polymer and Langevin simulation tutorial"

biermanncarl commented 3 years ago
* add density-dependent diffusion ($\sim N^3/2$?)

According to Wikipedia, D should behave like 1/P, which is proportional to 1/rho. Unfortunately, I couldn't get this scaling out of a Langevin Dynamics setup. The problem is that kinetic gas theory expects the particles to move in a straight line, while Langevin Dynamics prevents exactly this. Lowering the friction coefficient will render the simulation numerically unstable, i.e. the particles will slowly gain kinetic energy. The only way to deal with this is to lower the time step, which increases computational effort to get useful statistics. Overall I have the impression that using this result in a Langevin Dynamics tutorial has proven to be unfeasible.

jngrad commented 3 years ago

We need 2 people: one to look at the Langevin tutorial, one to look at the LB tutorial.