Skip to content
Snippets Groups Projects
Select Git revision
  • 3D-two-matrices
  • EffectiveDrop
  • Fix-mass-conservation
  • Reducing-input-parameters
  • Test_alpha_sensitivity
  • Troubleshooting-oscillations
  • flux_accuracy
  • master default
  • test_parallel
9 results
You can move around the graph by using the arrow keys.
Created with Raphaël 2.2.029Mar28252423222117161514832122Feb212015Jan24Nov26Oct2019141314Sep31Aug271091Apr30Mar231914922Feb201413108754330Jan24191816158Dec227Nov18171611943229Oct2623222120191413729Sep2825242221181716159875331Aug302726252418121065431Jul30292221201918171614131279Jun22May2015121110328Apr272220161412118762127Mar26252423201917151412109625Feb241716141312111031Jan2015146330Dec1918328Nov2726251918Implemented different Δx inside and outside in function "solve_D_adaptedGrid". Keeping solve_D separate for alpha analysis.setting α=1 by adapting Δt slightly improves mass conservation.Plotting m/m0 for different N, for 1D case and 3D case. Spherical mass conservation comparable to 1D case, but oscillating.Changed half-partitioning according to volume integral in spherical coordinates for 3D case.make_A, calc_b and solve_D now depend on geometry.Added percentage recovery time for the spherical case.spherical solution is working.Update example analytical 1D solution.NOW it works.calc_b_spherical addedfunction make_A_spherical added.1D analytical index issue solvedChanged second GP with D_in and D_out. Plot still to fixAdded half partitioning timeProblem in analytical (D_m,D_p inverted). Analytical outside fixed. Analytical inside wip.Introduced D_in D_out in solve_D. Matches analytical solution just if D_in=D_outIntroduce final time T as parameter.Change labels, timepoints, comment.Unify function argument order. Added doc strings.Comparison with Analytical solution works.Final interval now from 0 to 1.Smaller D gets rid of discontinuous boundary intensity increase.Comparison with Matlab numerics works.Different time steps work. alpha is now function.Creating matrix A is now a function make_A.Seems to work. Now try with larger time steps after initial equilibration.Small time step and setting one point right of boundary to 0 seems to help.Comparison with Matlab. Realised oscillation problem is in the initial conditions.Matrices for b and u are now output and plotted.Merge remote-tracking branch 'origin/EffectiveDrop' into EffectiveDropEffectiveDropEffectiveDropIntroduce solution arrayeff-Drop ghost points work above a certain number of cycleseff_Drop_DinDoutC.jl ghost points workingChange length to len.Remove factor 1/2 from alpha definition.eff_Drop_DinDoutC.jl added external GP, cut_ghost cuts external and boundary GP and returns X, U.Analytical solution starting inside and outside work.committing eff_Drop_DinDoutC.jl againRevert for Marta's code to work.Plotting optimized for interfacial resistance.