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.020May829Apr282721151187654131Mar302928252423222117161514832122Feb212015Jan24Nov26Oct2019141314Sep31Aug271091Apr30Mar231914922Feb201413108754330Jan24191816158Dec227Nov18171611943229Oct2623222120191413729Sep2825242221181716159875331Aug302726252418121065431Jul30292221201918171614131279Jun22May2015121110328Apr272220161412118762127Mar26252423201917151412109625Feb241716141312111031Jan2015146330Dec1918328Nov2726251918Cosmetics.Adding plots for Dresden application (chemical reaction recovery).With k not working in InterfaceTry. NEXT: start fixing from line 369Not breaking for half-times.eff drop adating delta_x for D,1, adapting delta_D to N=51, alpha=1 for D>=1.Minor fix on adaptive delta_t.u_in_eq corrected based on geometry. Still not perfectly matching.Now MC without ghost points. problem with u_in_eq 3D.now u_in_eq correct based on geometry.k (interfacial resistance) is now parameter. Add plotting function.Add plotting script.Plotting now in different file.Try out interface resistance boundary conditions.NOW time stepping is working, indexing problem should be solved.Now final step always landing exactly at T. Fixed last value of t_vec.Steepest gradient set inside in case of Boundary="E".Assigning u0 if IC="Exp_data" and Boundary="E"Layout fixed.solved problem LinearInterpolation for last delta-t exceeding lat timestep.WIP on computing derivatives for the next time step explicitly within the while loop.Introduced function adapt_Δt, working for the case theta=0 (=keeping Δt constant).Now looping while t<=T. First step for adaptive time stepping.Merge remote-tracking branch 'origin/Troubleshooting-oscillations' into Troubleshooting-oscillationsLAST VERSION LOOPING: for i = 3:Nt.Now steepest gradient computed with central difference.NOW x=0 working for the 3D case. and solve_D computing steepest gradient at each timepoint.x=0 working for the 3D case.External ghost points not at x:L-Δxi,R+Δxo. Working for the 1D case. 3D case?make_A optimized with anonymous function. Running time ca. same.calc_b optimized with anonymous function. Runs faster!There was still a mistake in previous version. Working nowRevert "Now EB working"Now EB workingGiving Experimental Boundary from simulation is working!calc_b working for Neumann boundary condition ("N") and experimental ("E").make_A working for Neumann boundary condition ("N") and experimental ("E").make_A_inside_EB working! Tomorrow will integrate with make_A.AdaptedGrid incorporated in solve_D, function of parameter AG ("Y"|"N")make_A changed to set boundary L and R to zero fluxcleaned, all analysis in frap_theory_analysis/analysis.ipynb