diff --git a/src/TestReactors.jl b/src/TestReactors.jl index d05e97f7..4b05c6d3 100644 --- a/src/TestReactors.jl +++ b/src/TestReactors.jl @@ -76,7 +76,7 @@ using DataFrames interfaces = [kLAkHCondensationEvaporationWithReservoir(domain, conds)] react = Reactor(domain, y0, (0.0, 140000.01), interfaces; p=p) - sol1 = solve(react.ode, react.recommendedsolver, abstol=1e-18, reltol=1e-6) + sol1 = solve(react.ode, react.recommendedsolver, abstol=1e-16, reltol=1e-6) phaseDict = readinput("../src/testing/TdependentkLAkH.rms") spcs = phaseDict["phase"]["Species"] @@ -90,7 +90,7 @@ using DataFrames interfaces = [kLAkHCondensationEvaporationWithReservoir(domain, conds)] react = Reactor(domain, y0, (0.0, 140000.01), interfaces; p=p) #Create the reactor object - sol2 = solve(react.ode, react.recommendedsolver, abstol=1e-18, reltol=1e-6) + sol2 = solve(react.ode, react.recommendedsolver, abstol=1e-16, reltol=1e-6) spcnames = getfield.(liq.species, :name) octaneind = findfirst(isequal("octane"), spcnames) @@ -136,7 +136,7 @@ using DataFrames domains = (domainliq, domaingas) interfaces = [vl, inletgas, outletgas] react, y0, p = Reactor(domains, (y0liq, y0gas), (0.0, tf), interfaces, (pliq, pgas, pinter)) - sol = solve(react.ode, react.recommendedsolver, abstol=1e-18, reltol=1e-6) + sol = solve(react.ode, react.recommendedsolver, abstol=1e-16, reltol=1e-6) name = "oxygen" ind = findfirst(x -> x == name, liqspcnames) @@ -581,7 +581,7 @@ using DataFrames react, y0, p = Reactor((domaingas, domaincat), (y0gas, y0cat), (0.0, 0.1), (inter,), (pgas, pcat, pinter)) - sol = solve(react.ode, CVODE_BDF(), abstol=1e-20, reltol=1e-6) + sol = solve(react.ode, CVODE_BDF(), abstol=1e-16, reltol=1e-6) ssys = SystemSimulation(sol, (domaingas, domaincat,), (inter,), p)