Drainage Time

By default, the drainage time is set to a constant value in the whole domain. In this case t_drain is simply set to an integer or float. We can also set t_drain to an array like of length Nx. In this case t_drain will vary accordingly with x.

Let’s take a look at a quick example. Let’s assume we want t_drain to decrease linearly with x. We could implement this as follows:

    import numpy as np
    from blobmodel import Model

    t_drain = np.linspace(2, 1, 100)

    tmp = Model(
        Nx=100,
        Ny=1,
        Lx=10,
        Ly=0,
        dt=1,
        T=1000,
        t_drain=t_drain,
        periodic_y=False,
        num_blobs=10,
    )
    tmp.make_realization()

The time averaged x-profile of n compared to a constant t_drain = 2 would then look like this: