Boundary treatment `pow`

A boundary treatment based on the work of Ebadi et. al. (2020).

from mnflow.mfda.cad.dld.theme.block import DLD

dld = DLD(
    Np=10,
    Nw=8,
    gap_w=1.,
    boundary_treatment='pow',
    num_unit=2,
    rot_last=90,
    # constraints OFF
    turn_off_constraints=True,
    # image
    opt_save_image=True,
    img_dpu=40,
    # for compatibility; added in v.0.0.1
    opt_acc_balance_pressure=False,
)

----------------------------------------
core.DLD___Np:10_Nw:8_gap_w:1.000_pitch_w:2.000_gap_a:1.000_pitch_a:2.000_height:4.000_boundary_treatment:pow
block.DLD___num_unit:2_opt_mirror:False_array_counts:[1, 1]_opt_mirror_before_array:[False, False]
----------------------------------------
{'Np': 10,
 'Nw': 8,
 'area': 847.7199999999999,
 'bb': [(-39.0, -5.193), (1.0, 16.0)],
 'count of 1D arrays of core.DLD': 1,
 'd_c': 0.4635835700756275,
 'lx': 40.0,
 'ly': 21.192999999999998,
 'performance': {'Flow rate @ 1 bar/area (micro-liter/min/mm-sq)': 674.3980393789292,
                 'die area (mm-sq)': 0.0008477199999999999,
                 'gap over crit. dia.': 2.157108371715726,
                 'volumetric flow rate at 1 bar (micro-liter/min)': 0.5717007059423058},
 'resistance (Pa.sec/m^3)': 1.0495001908578192e+16,
 'volumetric flow rate at 1 bar (m^3/sec)': 9.528345099038431e-12,
 'volumetric flow rate at 1 bar (milli-liter/hr)': 0.03430204235653835}

Output layout:

Layout