Issue with acceleration of colloidal gel poiseuille duct flow

34 views
Skip to first unread message

clara bae

unread,
Apr 25, 2023, 8:30:47 PM4/25/23
to hoomd-users
Dear all,

Hello, I am running a simulation with hoomd v3.1 and I am facing issue with simulation accelerating.

I am currently running a poiseuille flow between two walls. The box size is 300x30x30 (x y z). Initially I charge the particle and use the hoomd.md.external.field.Electric() to impose body force. 

Simulation typically do not accelerate for F=0.001 and 0.002 but things starts to accelerate at F=0.003. I am imposing different attraction strength for colloid-colloid and colloid-wall. Typically, for colloid-colloid of D0=6kT it shows acceleration no matter what the colloid-wall interaction is. However for colloid-colloid of D0=12kT colloid-wall D0=12kT, it was not accelerating. (but with colloid-colloid of D=12kT and colloid-wall D0=6kT, the simulation plateau for some period of time and start accelerating again.)

I cannot seem to solve this issue.

Below I have attached the plots of <average velocity in flow direction >vs <time> and python simulation code. 

Any suggestion to resolve this issue is greatly appreciate it. 
Thank you
sim_poi_p0.1.py
Screen Shot 2023-04-25 at 8.30.05 PM.png
Screen Shot 2023-04-25 at 8.15.21 PM.png
Screen Shot 2023-04-25 at 8.30.28 PM.png

Joshua Anderson

unread,
May 1, 2023, 11:00:26 AM5/1/23
to hoomd...@googlegroups.com
Clara,

This is a question for your research advisor. I do not find it surprising that a applying a constant force leads to an acceleration.
------
Joshua A. Anderson, Ph.D.
Research Area Specialist, Chemical Engineering, University of Michigan

--
You received this message because you are subscribed to the Google Groups "hoomd-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to hoomd-users...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/hoomd-users/5eab3efc-eb4d-4a8c-8821-32f6f25f2a30n%40googlegroups.com.
<sim_poi_p0.1.py><Screen Shot 2023-04-25 at 8.30.05 PM.png><Screen Shot 2023-04-25 at 8.15.21 PM.png><Screen Shot 2023-04-25 at 8.30.28 PM.png>

clara bae

unread,
May 4, 2023, 12:40:16 AM5/4/23
to hoomd-users
Oh, ok. Thank you!
Reply all
Reply to author
Forward
0 new messages