Skip to content

Commit

Permalink
Changed variable naming in TestReactors
Browse files Browse the repository at this point in the history
`cat` is a default Julia function for concatenating
arrays. As a precaution these are renamed to
`catalyst` to avoid conflicts
  • Loading branch information
ssun30 committed Nov 7, 2024
1 parent f9a0509 commit b582e25
Showing 1 changed file with 6 additions and 6 deletions.
12 changes: 6 additions & 6 deletions src/TestReactors.jl
Original file line number Diff line number Diff line change
Expand Up @@ -536,15 +536,15 @@ using DataFrames
interfacerxns = phaseDict[Set(["gas", "surface"])]["Reactions"]

ig = IdealGas(gasspcs, gasrxns; name="gas")
cat = IdealSurface(surfacespcs, surfacerxns, 2.486e-5; name="surface")
catalyst = IdealSurface(surfacespcs, surfacerxns, 2.486e-5; name="surface")

initialconds = Dict(["T" => 800.0, "P" => 1.0e5, "O2" => 0.2, "N2" => 0.7, "CH4" => 0.1])
domaingas, y0gas, pgas = ConstantTPDomain(phase=ig, initialconds=initialconds,)

V = 8.314 * 800.0 / 1.0e5
A = 1.0e5 * V
initialconds = Dict(["T" => 800.0, "A" => A, "vacantX" => cat.sitedensity * A])
domaincat, y0cat, pcat = ConstantTAPhiDomain(phase=cat, initialconds=initialconds,)
initialconds = Dict(["T" => 800.0, "A" => A, "vacantX" => catalyst.sitedensity * A])
domaincat, y0cat, pcat = ConstantTAPhiDomain(phase=catalyst, initialconds=initialconds,)

inter, pinter = ReactiveInternalInterfaceConstantTPhi(domaingas, domaincat, interfacerxns, 800.0, A)

Expand All @@ -567,15 +567,15 @@ using DataFrames
interfacerxns = phaseDict[Set(["gas", "surface"])]["Reactions"]

ig = IdealGas(gasspcs, gasrxns; name="gas")
cat = IdealSurface(surfacespcs, surfacerxns, 2.486e-5; name="surface")
catalyst = IdealSurface(surfacespcs, surfacerxns, 2.486e-5; name="surface")

initialconds = Dict(["T" => 800.0, "P" => 1.0e5, "O2" => 0.2, "N2" => 0.7, "CH4" => 0.1])
domaingas, y0gas, pgas = ConstantVDomain(phase=ig, initialconds=initialconds,)

V = 8.314 * 800.0 / 1.0e5
A = 1.0e5 * V
initialconds = Dict(["T" => 800.0, "A" => A, "vacantX" => cat.sitedensity * A])
domaincat, y0cat, pcat = ConstantTAPhiDomain(phase=cat, initialconds=initialconds,)
initialconds = Dict(["T" => 800.0, "A" => A, "vacantX" => catalyst.sitedensity * A])
domaincat, y0cat, pcat = ConstantTAPhiDomain(phase=catalyst, initialconds=initialconds,)

inter, pinter = ReactiveInternalInterface(domaingas, domaincat, interfacerxns, A)

Expand Down

0 comments on commit b582e25

Please sign in to comment.