ConstantReservoirDiffusion Interface

src/Domain.jl

@@ -1852,6 +1852,8 @@ export calcthermo
             cond = kLAs .* inter.molefractions .* inter.P ./ kHs * V
 
             dydt[d.indexes[1]:d.indexes[2]] .-= (evap .- cond)
+        elseif isa(inter, ConstantReservoirDiffusion) && d == inter.domain 
+            dydt[d.indexes[1]:d.indexes[2]] = inter.A .* d.diffusivity .* (inter.c .- ns./V) / inter.layer_thickness
         elseif isa(inter, VolumetricFlowRateInlet) && d == inter.domain
             dydt[d.indexes[1]:d.indexes[2]] .+= inter.Vin(t) * inter.cs
         elseif isa(inter, VolumetricFlowRateOutlet) && d == inter.domain
@@ -3084,6 +3086,12 @@ end
             @simd for i in domain.indexes[1]:domain.indexes[2]
                 @inbounds @fastmath jac[i, i] -= inter.Vout(t) / V
             end
+        elseif isa(inter, ConstantReservoirDiffusion) && domain == inter.domain 
+            # dn/dt .+= inter.A .* domain.diffusivity .* (inter.c .- ns./V) / inter.layer_thickness
+            # d/dni(dni/dt) .-= (inter.A .* d.diffusivity / (V * inter.layer_thickness) 
+            @simd for i in domain.indexes[1]:domain.indexes[2]
+                @inbounds @fastmath jac[i, i] -= inter.A * domain.diffusivity / (V * inter.layer_thickness)
+            end
         end
     end
 

src/Interface.jl

@@ -271,11 +271,23 @@ struct Outlet{V,FF<:Function} <: AbstractBoundaryInterface
 end
 export Outlet
 
+struct ConstantReservoirDiffusion{S,V<:AbstractArray,U<:Real}
+    domain::S
+    c::V
+    A::U
+    layer_thickness::U
+end
+
+function ConstantReservoirDiffusion(domain::V, conddict::Dict{X1,X}, A::B, layer_thickness::B) where {V,X1,X,B<:Real}
+    cs = makespcsvector(domain.phase, conddict)
+    return ConstantReservoirDiffusion(domain,cs,A,layer_thickness)
+end
+
 """
 kLAkHCondensationEvaporationWithReservoir adds evaporation and condensation to
 (1) a liquid phase domain with a constant composition vapor resevoir, where number of moles, P, and T need to be specified, or
 (2) a gas phase domain with a constant composition liquid resevoir, where number of moles, V, and T need to be specified.
-kLA and kH are used to model cond/evap. 
+kLA and kH are used to model cond/evap.
 kLA is liquid volumetric mass transfer coefficient with unit 1/s , and kH is Henry's law constant defined as gas phase partial pressure of solute over liquid phase concentration of solute.
 """
 
@@ -363,7 +375,7 @@ end
 export VolumetricFlowRateInlet
 
 """
-VolumeMaintainingOutlet is designed for gas phase domain such that the flow rate of this outlet will adjust to maintain the volume of the 
+VolumeMaintainingOutlet is designed for gas phase domain such that the flow rate of this outlet will adjust to maintain the volume of the
     domain to be constant. This is particularly useful to simulate any vapor-liquid phase system where the gas phase outlet
     is determined by the amount of evaporation.
 """
@@ -538,4 +550,4 @@ function evaluate(ri::FragmentBasedReactiveFilmGrowthInterfaceConstantT, dydt, V
     end
 end
 
-export evaluate
+export evaluate