Skip to content
Snippets Groups Projects
Select Git revision
  • c93b5cbaa1682f53702c30b585e15e47f60bdbbc
  • master default protected
  • suffa/cumulantfourth_order_correction_with_psm
  • mr_refactor_wfb
  • Sparse
  • WallLaw
  • improved_comm
  • release/1.3.7
  • release/1.3.6
  • release/1.3.5
  • release/1.3.4
  • release/1.3.3
  • release/1.3.2
  • release/1.3.1
  • release/1.3
  • release/1.2
  • release/1.1.1
  • release/1.1
  • release/1.0.1
  • release/1.0
  • release/0.4.4
  • release/0.4.3
  • release/0.4.2
  • release/0.4.1
  • release/0.4.0
  • release/0.3.4
  • release/0.3.3
27 results

stencils.py

Blame
    • stencils.py 4.28 KiB
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    14
    15
    16
    17
    18
    19
    20
    21
    22
    23
    24
    25
    26
    27
    28
    29
    30
    31
    32
    33
    34
    35
    36
    37
    38
    39
    40
    41
    42
    43
    44
    45
    46
    47
    48
    49
    50
    51
    52
    53
    54
    55
    56
    57
    58
    59
    60
    61
    62
    63
    64
    65
    66
    67
    68
    69
    70
    71
    72
    73
    74
    75
    76
    77
    78
    79
    80
    81
    82
    83
    84
    85
    86
    87
    88 
    

    

    

    def get_stencil(name, ordering='walberla'):
    """
    Stencils are tuples of discrete velocities. They are commonly labeled in the 'DxQy' notation, where d is the
    dimension (length of the velocity tuples) and y is number of discrete velocities.

    :param name: DxQy notation
    :param ordering: the LBM literature does not use a common order of the discrete velocities, therefore here
                     different common orderings are available. All orderings lead to the same method, it just has
                     to be used consistently. Here more orderings are available to compare intermediate results with
                     the literature.
    """
    try:
        return get_stencil.data[name.upper()][ordering.lower()]
    except KeyError:
        err_msg = ""
        for stencil, ordering_names in get_stencil.data.items():
            err_msg += "  %s: %s\n" % (stencil, ", ".join(ordering_names.keys()))

        raise ValueError("No such stencil available. "
                         "Available stencils: <stencil_name>( <ordering_names> )\n" + err_msg)


get_stencil.data = {
    'D2Q9': {
        'walberla': ((0, 0),
                     (0, 1), (0, -1), (-1, 0), (1, 0),
                     (-1, 1), (1, 1), (-1, -1), (1, -1),),
        'counterclockwise': ((0, 0),
                             (1, 0), (0, 1), (-1, 0), (0, -1),
                             (1, 1), (-1, 1), (-1, -1), (1, -1)),
        'braunschweig': ((0, 0),
                         (-1, 1), (-1, 0), (-1, -1), (0, -1),
                         (1, -1), (1, 0), (1, 1), (0, 1)),
        'uk': ((0, 0),
               (1, 0), (-1, 0), (0, 1), (0, -1),
               (1, 1), (-1, -1), (-1, 1), (1, -1),
               )
    },
    'D3Q15': {
        'walberla':
            ((0, 0, 0),
             (0, 1, 0), (0, -1, 0), (-1, 0, 0), (1, 0, 0), (0, 0, 1), (0, 0, -1),
             (1, 1, 1), (-1, 1, 1), (1, -1, 1), (-1, -1, 1), (1, 1, -1), (-1, 1, -1), (1, -1, -1), (-1, -1, -1)),
        'premnath': ((0, 0, 0),
                     (1, 0, 0), (-1, 0, 0), (0, 1, 0), (0, -1, 0), (0, 0, 1), (0, 0, -1),
                     (1, 1, 1), (-1, 1, 1), (1, -1, 1), (-1, -1, 1),
                     (1, 1, -1), (-1, 1, -1), (1, -1, -1), (-1, -1, -1)),
    },
    'D3Q19': {
        'walberla': ((0, 0, 0),
                     (0, 1, 0), (0, -1, 0), (-1, 0, 0), (1, 0, 0), (0, 0, 1), (0, 0, -1),
                     (-1, 1, 0), (1, 1, 0), (-1, -1, 0), (1, -1, 0),
                     (0, 1, 1), (0, -1, 1), (-1, 0, 1), (1, 0, 1),
                     (0, 1, -1), (0, -1, -1), (-1, 0, -1), (1, 0, -1)),
        'braunschweig': ((0, 0, 0),
                         (1, 0, 0), (-1, 0, 0),
                         (0, 1, 0), (0, -1, 0),
                         (0, 0, 1), (0, 0, -1),
                         (1, 1, 0), (-1, -1, 0),
                         (1, -1, 0), (-1, 1, 0),
                         (1, 0, 1), (-1, 0, -1),
                         (1, 0, -1), (-1, 0, 1),
                         (0, 1, 1), (0, -1, -1),
                         (0, 1, -1), (0, -1, 1)),
        'premnath': ((0, 0, 0),
                     (1, 0, 0), (-1, 0, 0), (0, 1, 0), (0, -1, 0), (0, 0, 1), (0, 0, -1),
                     (1, 1, 0), (-1, 1, 0), (1, -1, 0), (-1, -1, 0),
                     (1, 0, 1), (-1, 0, 1), (1, 0, -1), (-1, 0, -1),
                     (0, 1, 1), (0, -1, 1), (0, 1, -1), (0, -1, -1)),
    },
    'D3Q27': {
        'walberla': ((0, 0, 0),
                     (0, 1, 0), (0, -1, 0), (-1, 0, 0), (1, 0, 0), (0, 0, 1), (0, 0, -1),
                     (-1, 1, 0), (1, 1, 0), (-1, -1, 0), (1, -1, 0),
                     (0, 1, 1), (0, -1, 1), (-1, 0, 1), (1, 0, 1),
                     (0, 1, -1), (0, -1, -1), (-1, 0, -1), (1, 0, -1),
                     (1, 1, 1), (-1, 1, 1), (1, -1, 1), (-1, -1, 1), (1, 1, -1), (-1, 1, -1), (1, -1, -1),
                     (-1, -1, -1)),
        'premnath': ((0, 0, 0),
                     (1, 0, 0), (-1, 0, 0), (0, 1, 0), (0, -1, 0), (0, 0, 1), (0, 0, -1),
                     (1, 1, 0), (-1, 1, 0), (1, -1, 0), (-1, -1, 0),
                     (1, 0, 1), (-1, 0, 1), (1, 0, -1), (-1, 0, -1),
                     (0, 1, 1), (0, -1, 1), (0, 1, -1), (0, -1, -1),
                     (1, 1, 1), (-1, 1, 1), (1, -1, 1), (-1, -1, 1),
                     (1, 1, -1), (-1, 1, -1), (1, -1, -1), (-1, -1, -1))
    }
}