Skip to content
Snippets Groups Projects
Select Git revision
  • 11f0c22b82ae81ced6133afbe2ea23cd75b4ed13
  • 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

test_analytical.py

Blame
    • Markus Holzer's avatar
    Adapted Parameters
    Markus Holzer authored
    4813be0c
    History
    Markus Holzer's avatar 4813be0c
    History
    test_analytical.py 2.58 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 
    

    

    

    import numpy as np

from lbmpy.phasefield_allen_cahn.parameter_calculation import calculate_dimensionless_rising_bubble, \
    calculate_parameters_rti, calculate_parameters_taylor_bubble

from lbmpy.phasefield_allen_cahn.analytical import analytic_rising_speed


def test_analytical():
    parameters = calculate_dimensionless_rising_bubble(reference_time=18000,
                                                       density_heavy=1.0,
                                                       bubble_radius=16,
                                                       bond_number=30,
                                                       reynolds_number=420,
                                                       density_ratio=1000,
                                                       viscosity_ratio=100)

    assert np.isclose(parameters.density_light, 0.001)
    assert np.isclose(parameters.gravitational_acceleration, -9.876543209876543e-08)

    parameters = calculate_parameters_rti(reference_length=128,
                                          reference_time=18000,
                                          density_heavy=1.0,
                                          capillary_number=9.1,
                                          reynolds_number=128,
                                          atwood_number=1.0,
                                          peclet_number=744,
                                          density_ratio=3,
                                          viscosity_ratio=3)

    assert np.isclose(parameters.density_light, 1/3)
    assert np.isclose(parameters.gravitational_acceleration, -3.9506172839506174e-07)
    assert np.isclose(parameters.mobility, 0.0012234169653524492)

    rs = analytic_rising_speed(1 - 6, 20, 0.01)
    assert np.isclose(rs, 16666.666666666668)

    parameters = calculate_parameters_taylor_bubble(reference_length=192,
                                                    reference_time=36000,
                                                    density_heavy=1.0,
                                                    diameter_outer=0.0254,
                                                    diameter_inner=0.0127)

    assert np.isclose(parameters.density_heavy, 1.0)
    assert np.isclose(parameters.density_light, 0.001207114228456914)
    assert np.isclose(parameters.dynamic_viscosity_heavy, 5.733727652152216e-05)
    assert np.isclose(parameters.dynamic_viscosity_light, 1.0417416358027054e-06)
    assert np.isclose(parameters.gravitational_acceleration, -7.407407407407407e-08)
    assert np.isclose(parameters.surface_tension, 3.149857262258028e-05, rtol=1e-05)