Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

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

Closed
jngrad opened this issue Oct 8, 2020 · 10 comments · Fixed by #4052
Closed

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

jngrad opened this issue Oct 8, 2020 · 10 comments · Fixed by #4052

Comments

@jngrad
Copy link
Member

jngrad commented Oct 8, 2020

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
Copy link
Member Author

jngrad commented Oct 8, 2020

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
Copy link
Contributor

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
Copy link
Contributor

Offline discussion with Christian:

  • Further split into a Langevin dynamics tutorial (now polymer part 1)
  • add density-dependent diffusion ($\sim N^3/2$?)
  • use more appropriate names for LB tutorial (theory and Poiseuille flow instead of part 1 and part 2)
  • Check is visualization of the sedimentation can be done (see my note above). @jngrad can you comment on this?
    Updating issue title and closing PR until this is added.

@schlaicha
Copy link
Contributor

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

@jngrad
Copy link
Member Author

jngrad commented Dec 18, 2020

@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
Copy link
Contributor

schlaicha commented Dec 18, 2020

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

No, also I cannot add myself as assignee.

@jngrad
Copy link
Member Author

jngrad commented Dec 18, 2020

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
Copy link
Contributor

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

@jngrad jngrad changed the title Move polymer part of LB tutorial to a polymer tutorial Split LB tutorial into LB, polymer and Langevin simulation tutorials Dec 18, 2020
@biermanncarl
Copy link
Contributor

* 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
Copy link
Member Author

jngrad commented Aug 9, 2021

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

@kodiakhq kodiakhq bot closed this as completed in #4052 Aug 17, 2021
kodiakhq bot added a commit that referenced this issue Aug 17, 2021
#4052)

Fixes #3939

Description of changes:
- Create new tutorial 'polymers'
- Move LB tutorials part 2+3 to polymers
- Update numbering for consistency
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

Successfully merging a pull request may close this issue.

6 participants