コード例 #1
0
ファイル: test_unit.py プロジェクト: McBs/lettuce
def test_consistency(lattice):
    units = UnitConversion(lattice, reynolds_number=1000, mach_number=0.05, characteristic_length_lu=100,
                           characteristic_length_pu=2*np.pi, characteristic_velocity_pu=2, characteristic_density_pu=0.7)
    rho = 0.9
    u = 0.1
    assert (
        units.convert_density_to_pu(rho)*units.convert_velocity_to_pu(u)**2
        == pytest.approx(units.convert_energy_to_pu(rho*u*u))
    )
    assert units.convert_velocity_to_pu(u)**2 == pytest.approx(units.convert_incompressible_energy_to_pu(u*u))
コード例 #2
0
ファイル: test_unit.py プロジェクト: McBs/lettuce
def test_conversion_reversible(lattice):
    units = UnitConversion(lattice, reynolds_number=1000, mach_number=0.05, characteristic_length_lu=100,
                           characteristic_length_pu=2*np.pi, characteristic_velocity_pu=2, characteristic_density_pu=0.7)
    assert units.convert_velocity_to_lu(units.convert_velocity_to_pu(2.0)) == pytest.approx(2.0)
    assert units.convert_time_to_lu(units.convert_time_to_pu(2.0)) == pytest.approx(2.0)
    assert units.convert_length_to_lu(units.convert_length_to_pu(2.0)) == pytest.approx(2.0)
    assert units.convert_density_to_lu(units.convert_density_to_pu(2.0)) == pytest.approx(2.0)
    assert units.convert_pressure_to_lu(units.convert_pressure_to_pu(2.0)) == pytest.approx(2.0)
    assert units.convert_density_lu_to_pressure_pu(units.convert_pressure_pu_to_density_lu(2.0)) == pytest.approx(2.0)
    assert units.convert_energy_to_lu(units.convert_energy_to_pu(2.0)) == pytest.approx(2.0)
    assert units.convert_incompressible_energy_to_lu(units.convert_incompressible_energy_to_pu(2.0)) == pytest.approx(2.0)