[WIP] Simultaneous network solve: injection and production networks, autochoke, gas-lift response - #7348

Draft

hnil wants to merge 80 commits into

hnil:simultaneous-network

Commits on Aug 24, 2026

Update injection network leaf node rates.

Atgeirr Flø Rasmussen
authored and
hnil
committed
470dd52
Add anyNetworkActive() and activeNetworks() helpers.

Atgeirr Flø Rasmussen
authored and
hnil
committed
93805ab
Make updateActiveState() honor injection networks.

Atgeirr Flø Rasmussen
authored and
hnil
committed
6ad4125
Check active status of all networks.

Atgeirr Flø Rasmussen
authored and
hnil
committed
f449efb
WIP stage 1

Atgeirr Flø Rasmussen
authored and
hnil
committed
80c1381
Check well VFP table before setting dynamic THP limit.

Atgeirr Flø Rasmussen
authored and
hnil
committed
bdd2402
Check correct rates storage for empty leaf node guard.

Atgeirr Flø Rasmussen
authored and
hnil
committed
2d45479
Suggestions from copilot

Atgeirr Flø Rasmussen
authored and
hnil
committed
cb0046d
More copilot.

Atgeirr Flø Rasmussen
authored and
hnil
committed
1dca228
Do not extrapolate VFP tables for network branches

hnil
and
claude
committed
700b05a
Bracketing update of network node pressures; fix injector coupling

hnil
and
claude
committed
b88a7a3
Do not shut an injector that cannot operate at the network THP

hnil
and
claude
committed
e1e08c8
Store injection network leaf rates per phase

hnil
and
claude
committed
23dd772
Export injection network node/branch values for GPRG / GPRW

hnil
and
claude
committed
fe70531
Restrict the bracketing node-pressure update to injection networks by default

hnil
and
claude
committed
b1dc9f7
Add GNETINJE_GAS-01 / GNETINJE_WAT-01 regression tests; GNETINJE is supported

hnil
and
claude
committed
7d5283f
Re-solve the production network before the injection networks

hnil
and
claude
committed
c554b4d
Optional Anderson acceleration of the network pressures (off; does not help yet)

hnil
and
claude
committed
8d98377
Optional well-index proxy for the network balance (off; does not help yet)

hnil
and
claude
committed
8aa9ea6
Add a standalone bench for the injection-network solve

hnil
and
claude
committed
f9132d7
Test globalisation in the network bench

hnil
and
claude
committed
18988cb
Add the full network formulation to the bench, next to the eliminated one

hnil
and
claude
committed
4c32918
Bench: table limits as bounds, a group multiplier, and a case you can build

hnil
and
claude
committed
2838f48
Bench: build a case from a deck, with the operating point set or read

hnil
and
claude
committed
e4490eb
Bench: run the full formulation across both test cases

hnil
and
claude
committed
2e0a1c2
Move the full network formulation into the library

hnil
and
claude
committed
3646b29
Solve the injection networks simultaneously, behind --network-solver=newton

hnil
and
claude
committed
a70f852
Analytic network Jacobian and network group control, both opt in

hnil
and
claude
committed
779b814
Take network guide rates from well potential, not from the current rate

hnil
and
claude
committed
9e00f3d
Build the network solve from data every rank has

hnil
and
claude
committed
65ff13b
Say why the network solve handed back to the relaxed update

hnil
and
claude
committed
6544cfa
Report why a network solve stopped, and record the active set

hnil
and
claude
committed
4e7bac7
Cover the guide-refresh path in the bench, and rule it out on its own

hnil
and
claude
committed
1305c6f
Replay failed network solves in the bench, and fix what it found

hnil
and
claude
committed
bb4085e
Count every well the group allocated against its target

hnil
and
claude
committed
f89dc64
Prototype the same formulation for a production network

hnil
and
claude
committed
c0a767c
Fix the group row's derivative, and give the test a case that can see it

hnil
and
claude
committed
0f281cc
Check the group equations against the rule they replace

hnil
and
claude
committed
093853f
Decide each control from the rate it allows, and test the group logic harder

hnil
and
claude
committed
ab7accb
Resolve the group split where it is chosen, not across Newton iterations

hnil
and
claude
committed
33c0554
Give the production network the same control rule, and group control

hnil
and
claude
committed
f245860
Measure the flipping, and say where the new split does not help

hnil
and
claude
committed
fec5f81
Enforce rate limits from a converged point, and stop falling back

hnil
and
claude
committed
4cd511a
Solve production networks simultaneously too

hnil
and
claude
committed
821aced
Solve multi-root networks as forests, and refuse what the system cannot model

hnil
and
claude
committed
576fe56
Model efficiencies, lift gas and satellites instead of refusing them

hnil
and
claude
committed
bff9d42
Autochoke: take the target from the control the block decided on

hnil
and
claude
committed
29a1cd8
Autochoke inside the simultaneous solve, opt-in

hnil
and
claude
committed
5a4207c
Solve a frozen production network once per sub-loop

hnil
and
claude
committed
379d0bd
Answer gas lift trials from the network solve (experimental, off)

hnil
and
claude
committed
ee6bed8
Make the legacy autochoke's bracket sampling a parameter

hnil
and
claude
committed
d4d1323
Cache the tubing allowance on a pressure grid, exact across jumps

hnil
and
claude
committed
eb07333
Analytic Jacobian for the production system

hnil
and
claude
committed
564be5d
Production dump and replay, a tied control, and what the bench said

hnil
and
claude
committed
22fa09c
Complementarity rows for well limits, opt-in and not yet better

hnil
and
claude
committed
23fec38
Choke as a complementarity row, Shut control, scan-decided branch
hnil
committed
1100666
Complementarity: dead wells shut, death sticky, componentwise step
hnil
committed
4f298de
Complementarity: shut wells without inflow performance; dead-well tests
hnil
committed
efcaf39
Pin the active set's dead-well hole in a test
hnil
committed
fdb495f
Drop the two node-pressure experiments that did not pay off

hnil
and
claude
committed
7c25558
Review: one leaf-node entry per node, not per phase

hnil
and
claude
committed
b5a50cb
Review: carry the IPR in the simulator's own convention

hnil
and
claude
committed
1b32904
Review: let Dune do the dense solve

hnil
and
claude
committed
ace7891
Review: name the injection system, and say what solve() does

hnil
and
claude
committed
d37d21e
Review: no default arguments on solve()

hnil
and
claude
committed
4bad252

Commits on Aug 25, 2026

Say where the explicit group share could go next

hnil
and
claude
committed
ca143b4
Review: a single solve(), with the settings stated at every call

hnil
and
claude
committed
76dcc46
Review: one file per system, and the solver in its own

hnil
and
claude
committed
04a77a3

Commits on Aug 26, 2026

Commits on Aug 27, 2026

Refresh the branch pressure drops the simultaneous solve invalidates

hnil
and
claude
committed
7a8b7fa
Deck tests for what the network solve costs
hnil
committed
d0d32e2
Say what the complementarity rows mean, for a reader who was not there
hnil
committed
4cab575
Reach the production solve on plain network decks
hnil
committed
bfb8f0c
Agreement tests: the simultaneous solve against the relaxed update
hnil
committed
b10321a
Review: put the network solver headers in a network/ subdir

hnil
and
claude
committed
a2f4989

Commits on Aug 28, 2026

Copilot review: the domain mapping, the debug trace, and an indent

hnil
and
claude
committed
71efe14

Commits on Aug 31, 2026

Drop the last of the unused injection phase
hnil
committed
e60d49d