Issue with Drainage Simulation in OpenPNM (example in tutorials)

[purunyr]

Hello,

I am currently following the tutorials for drainage simulations, but I'm encountering a problem during the execution of the simulation.

Problem Description:
I have written code to create a network and set the parameters for a drainage simulation. However, when I run the simulation, it hangs with the output message:

Performing drainage simulation: 0%| | 0/25 [00:00<?, ?it/s]

This message persists without any further progress.

Here is my code (just same as the code provided in the tutorial):

---

import openpnm as op
import numpy as np
import matplotlib.pyplot as plt
op.visualization.set_mpl_style()

np.random.seed(5)
pn = op.network.Demo(shape=[25, 25, 1], spacing=1e-4)
air = op.phase.Air(network=pn)

air['pore.contact_angle'] = 120
air['pore.surface_tension'] = 0.072
f = op.models.physics.capillary_pressure.washburn
air.add_model(propname='throat.entry_pressure',
model=f,
surface_tension='throat.surface_tension',
contact_angle='throat.contact_angle',
diameter='throat.diameter',)

#Using the Drainage Algorithm
drn = op.algorithms.Drainage(network=pn, phase=air)
drn.set_inlet_BC(pores=pn.pores('left'))
drn.run()

---

I have tried various solutions, but it is not working :(

Questions:

1. What could be causing the simulation to hang at 0%?
2. Are there any additional checks or parameter settings that I might be missing?

Any guidance or suggestions would be greatly appreciated!

Thank you in advance for your assistance!