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

Renamed mobility

Browse source
This commit is contained in:
Mo8it 2022-01-23 16:14:37 +01:00
parent 03f8968ba0
commit 7aabbcb309
3 changed files with 6 additions and 6 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

2
src/data.jl
View file

@ -3,7 +3,7 @@ struct SimConsts
v₀::Float64 v₀::Float64
δt::Float64 δt::Float64
packing_ratio::Float64 packing_ratio::Float64
μ::Float64 μₜ::Float64
Dₜ::Float64 Dₜ::Float64
particle_radius::Float64 particle_radius::Float64
Dᵣ::Float64 Dᵣ::Float64

4
src/run.jl
View file

@ -84,12 +84,12 @@ function run_sim(
skin_r=sim_consts.skin_r, skin_r=sim_consts.skin_r,
skin_r²=sim_consts.skin_r^2, skin_r²=sim_consts.skin_r^2,
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.μₜ,
interaction_r=sim_consts.interaction_r, interaction_r=sim_consts.interaction_r,
interaction_r²=sim_consts.interaction_r^2, interaction_r²=sim_consts.interaction_r^2,
n_particles=sim_consts.n_particles, n_particles=sim_consts.n_particles,

6
src/setup.jl
View file

@ -6,7 +6,7 @@ const DEFAULT_PARENT_DIR = ""
const DEFAULT_COMMENT = "" const DEFAULT_COMMENT = ""
const DEFAULT_PARTICLE_RADIUS = 0.5 const DEFAULT_PARTICLE_RADIUS = 0.5
const DEFAULT_Dₜ = 1.0 const DEFAULT_Dₜ = 1.0
const DEFAULT_μ = 1.0 const DEFAULT_μₜ = 1.0
const DEFAULT_ϵ = 100.0 const DEFAULT_ϵ = 100.0
function initial_particle_grid_pos( function initial_particle_grid_pos(
@ -59,7 +59,7 @@ function gen_sim_consts(
@assert δt in 1e-7:1e-7:1e-4 @assert δt in 1e-7:1e-7:1e-4
@assert packing_ratio > 0 @assert packing_ratio > 0
μ = DEFAULT_μ μₜ = DEFAULT_μₜ
Dₜ = DEFAULT_Dₜ Dₜ = DEFAULT_Dₜ
particle_radius = DEFAULT_PARTICLE_RADIUS particle_radius = DEFAULT_PARTICLE_RADIUS
Dᵣ = 3 * Dₜ / ((2 * particle_radius)^2) Dᵣ = 3 * Dₜ / ((2 * particle_radius)^2)
@ -101,7 +101,7 @@ function gen_sim_consts(
δt, δt,
packing_ratio, packing_ratio,
# Internal # Internal
μ, μₜ,
Dₜ, Dₜ,
particle_radius, particle_radius,
Dᵣ, Dᵣ,