Skip to content

Commit 21cc358

Browse files
mjohnson541ssun30
authored andcommitted
fix FragmentBasedReactiveFilmGrowthInterfaceConstantT interface
1 parent 774a8e0 commit 21cc358

File tree

1 file changed

+6
-3
lines changed

1 file changed

+6
-3
lines changed

src/Interface.jl

Lines changed: 6 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -464,9 +464,12 @@ function FragmentBasedReactiveFilmGrowthInterfaceConstantT(domainfilm, domain2,
464464

465465
rxnarray, fragmentbasedrxnarray = getfragmentbasedinterfacereactioninds(domainfilm, domain2, reactions)
466466

467-
kfs = getkf.(reactions, nothing, T, 0.0, 0.0, Ref([]), 0.0, 0.0)
468-
Kc = getKc.(reactions, domainfilm.phase, domain2.phase, Ref(domainfilm.Gs), Ref(domain2.Gs), T, 0.0)
469-
krevs = kfs ./ Kc
467+
M,Nrp1,Nrp2 = getstoichmatrix(domainfilm,domain2,reactions)
468+
Gpart = ArrayPartition(domainfilm.Gs,domain2.Gs)
469+
dGrxns = -M*Gpart
470+
kfs = getkf.(reactions,nothing,T,0.0,0.0,Ref([]),0.0,0.0,dGrxns,0.0)
471+
Kc = getKcs(domain1.phase,domain2.phase,T,Nrp1,Nrp2,dGrxns)
472+
krevs = kfs./Kc
470473

471474
M, Nrp1, Nrp2 = getstoichmatrix(domainfilm, domain2, reactions)
472475
reversibility = Array{Bool,1}(getfield.(reactions, :reversible))

0 commit comments

Comments
 (0)