mo8it/ReCo.jl
1
0
Fork 0
mirror of https://gitlab.rlp.net/mobitar/ReCo.jl.git synced 2024-12-21 00:51:21 +00:00
Code Issues Releases Activity

Release hard fix of sigma

Browse source
This commit is contained in:
Mo8it 2022-01-20 17:33:17 +01:00
parent 658c33f5a2
commit 383a77e2bf
2 changed files with 5 additions and 5 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

2
src/run.jl
View file

@ -81,7 +81,7 @@ function run_sim(
n_snapshots=n_snapshots, n_snapshots=n_snapshots,
c₁=4 * sim_consts.ϵ * 6 * sim_consts.σ^6 * sim_consts.δt * sim_consts.μ, c₁=4 * sim_consts.ϵ * 6 * sim_consts.σ^6 * sim_consts.δt * sim_consts.μ,
c₂=2 * sim_consts.σ^6, c₂=2 * sim_consts.σ^6,
c₃=sqrt(2 * sim_consts.D₀ * sim_consts.δt), c₃=sqrt(2 * sim_consts.Dₜ * sim_consts.δt),
c₄=sqrt(2 * sim_consts.Dᵣ * sim_consts.δt), c₄=sqrt(2 * sim_consts.Dᵣ * sim_consts.δt),
v₀δt=sim_consts.v₀ * sim_consts.δt, v₀δt=sim_consts.v₀ * sim_consts.δt,
μ=sim_consts.μ, μ=sim_consts.μ,

8
src/setup.jl
View file

@ -55,10 +55,10 @@ function gen_sim_consts(
@assert packing_ratio > 0 @assert packing_ratio > 0
μ = 1.0 μ = 1.0
D₀ = 1.0 Dₜ = 1.0
particle_radius = 0.5 particle_radius = 0.5
Dᵣ = 3 * D₀ / ((2 * particle_radius)^2) Dᵣ = 3 * Dₜ / ((2 * particle_radius)^2)
σ = 1.0 σ = 2 * particle_radius
ϵ = 100.0 ϵ = 100.0
interaction_r = 2^(1 / 6) * σ interaction_r = 2^(1 / 6) * σ
@ -92,7 +92,7 @@ function gen_sim_consts(
packing_ratio, packing_ratio,
# Internal # Internal
μ, μ,
D₀, Dₜ,
particle_radius, particle_radius,
Dᵣ, Dᵣ,
σ, σ,