From 5a7db94c268425eb228ca1a9353b72d8f5d01421 Mon Sep 17 00:00:00 2001 From: MaxVZIKM <147824975+MaxVZIKM@users.noreply.github.com> Date: Wed, 3 Apr 2024 15:58:27 +0200 Subject: [PATCH 1/2] Update topology_optimization.jl Multiplication of the number of iterations needed for stability (n_j) with the viscosity (\eta) within the function update_density. If n_j is bigger than 1, behavior similar to checkerboarding could be observed otherwise. --- docs/src/literate-gallery/topology_optimization.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/docs/src/literate-gallery/topology_optimization.jl b/docs/src/literate-gallery/topology_optimization.jl index a1f26d730f..8ad7356c09 100644 --- a/docs/src/literate-gallery/topology_optimization.jl +++ b/docs/src/literate-gallery/topology_optimization.jl @@ -287,7 +287,7 @@ function update_density(dh, states, mp, ρ, topology, Δh) Δχ = pΨ/p_Ω + mp.β*∇²χ - χn1 = compute_χn1(χn, Δχ, ρ, mp.η, mp.χ_min) + χn1 = compute_χn1(χn, Δχ, ρ, mp.η*n_j, mp.χ_min) if(j Date: Thu, 18 Apr 2024 13:58:37 +0200 Subject: [PATCH 2/2] Added reference for density update --- docs/src/literate-gallery/topology_optimization.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/docs/src/literate-gallery/topology_optimization.jl b/docs/src/literate-gallery/topology_optimization.jl index 8ad7356c09..9fc9ab3bce 100644 --- a/docs/src/literate-gallery/topology_optimization.jl +++ b/docs/src/literate-gallery/topology_optimization.jl @@ -273,7 +273,7 @@ end #md nothing # hide # Finally, we put everything together to update the density. The loop ensures the stability of the -# updated solution. +# updated solution, see [JanHacJun2019regularizedthermotopopt](@cite) equation (42) for reference. function update_density(dh, states, mp, ρ, topology, Δh) n_j = Int(ceil(6*mp.β/(mp.η*Δh^2))) # iterations needed for stability