include("eff_Drop.jl") P = 10; Δt = 0.001; D_i = 1; D_o = 50; BP=1; x0=0.5; T=20; perc=0.5; L=0; R=10; k="none"; ## Solve diffusion equation in phase-separated media A, u, x, b, Nb, t_perc, αi, αo= @time solve_D("3D","FRAP","Y","N",zeros(9), P,D_i,D_o,0,Δt,L,R,BP,x0,T,perc,zeros(9),zeros(9), k) u = u[:, vcat(1:Nb-1, Nb+2:size(u, 2))] x = vcat(collect(x[1:Nb-1]), collect(x[Nb+2:end])) writedlm(string("interface_ref_none.csv"), hcat(x, u'), ',') ## k = 0.05 A, u1, x, b, Nb, t_perc, αi, αo= @time solve_D("3D","FRAP","Y","N",zeros(9), P,D_i,D_o,0,Δt,L,R,BP,x0,T,perc,zeros(9),zeros(9), k) u1 = u1[:, vcat(1:Nb-1, Nb+2:size(u1, 2))] x = vcat(collect(x[1:Nb-1]), collect(x[Nb+2:end])) writedlm(string("interface_ref_k_", k, ".csv"),hcat(x, u1'), ',') plot_model(x, u, Nb, T, Δt) plot_model(x, u1, Nb, T, Δt)