A new domain ConstantTVGasDomain

[donerancl]

I made a new domain for constant temperature, constant volume gas. The pressure is calculated from the number of moles as simulation progresses.

I have used it in a simulation with a vapor liquid interface with no issues.

[codecov]

Codecov Report

Attention: Patch coverage is 0% with 162 lines in your changes missing coverage. Please review.

Project coverage is 47.78%. Comparing base (1f37bfd) to head (a628c05).

❗ Current head a628c05 differs from pull request most recent head 512174c

Please upload reports for the commit 512174c to get more accurate results.

Files	Patch %	Lines
src/Domain.jl	0.00%	162 Missing ⚠️

Additional details and impacted files

@@            Coverage Diff             @@
##             main     #258      +/-   ##
==========================================
- Coverage   48.71%   47.78%   -0.94%     
==========================================
  Files          31       31              
  Lines        8313     8475     +162     
==========================================
  Hits         4050     4050              
- Misses       4263     4425     +162     

☔ View full report in Codecov by Sentry.
📢 Have feedback on the report? Share it here.

[donerancl]

I found a small bug in the rops function that we could include in this PR. In Simulation.jl line 254 raises an error if the array in the second term is empty (when using interfaces without reactions, e.g., VaporLiquidMassTransferInternalInterfaceConstantT)
Nrxns = sum([length(sim.domain.phase.reactions) for sim in ssys.sims]) + sum([length(inter.reactions) for inter in ssys.interfaces if hasproperty(inter, :reactions)])

[hwpang]

Thanks for the PR! I see that the test is failing with

Test vapor-liquid phase multi-domain reactor simulation with VaporLiquidMassTransferInternalInterfaceConstantT and VolumeMaintainingOutlet interface: Test Failed at /home/runner/work/ReactionMechanismSimulator.jl/ReactionMechanismSimulator.jl/src/TestReactors.jl:143
  Expression: ≈((sol(sol.t[end]))[ind], 0.11758959354431776, rtol = 1.0e-5)
   Evaluated: 1.8655128553425897 ≈ 0.11758959354431776 (rtol=1.0e-5)

Could you take a look into it?

[donerancl]

@hwpang Just saw this and took a look. The temperature (298 K) for the test was below the range used to extrapolate kH and kLA, so my branch calculated the values of kH(T) and kLA(T) at the lowest temperature in the given range (350 K), making the results different. We could allow extrapolation if the kLA and kH values are not negative. It is probably alright in this case I would think.