Skip to content
This repository has been archived by the owner on Feb 14, 2024. It is now read-only.

Commit

Permalink
Update optimize_simple_bl.jl
Browse files Browse the repository at this point in the history
Catch up to changes in default pressure setup
moyner committed Feb 5, 2023
1 parent c293da4 commit 8018656
Showing 1 changed file with 2 additions and 0 deletions.
2 changes: 2 additions & 0 deletions darcy/optimization/optimize_simple_bl.jl
Original file line number Diff line number Diff line change
@@ -18,6 +18,8 @@ function setup_bl(;nc = 100, time = 1.0, nstep = 100, poro = 0.1, perm = 9.8692e
# Define system and realize on grid
sys = ImmiscibleSystem((LiquidPhase(), VaporPhase()))
model = SimulationModel(G, sys)
# Avoid dampening and limits on pressure - simple synthetic case can have negative or zero pressures.
model.primary_variables[:Pressure] = Pressure(minimum = -Inf, max_rel = nothing)
kr = BrooksCoreyRelPerm(sys, [2.0, 2.0])
replace_variables!(model, RelativePermeabilities = kr)
tot_time = sum(tstep)

0 comments on commit 8018656

Please sign in to comment.