update boundary layer tests for swap to head's method

byuflowlab · Nov 8, 2024 · 8d859c1 · 8d859c1
1 parent 5d4826a
commit 8d859c1
Showing 1 changed file with 86 additions and 81 deletions.
diff --git a/test/boundary_layer_tests.jl b/test/boundary_layer_tests.jl
@@ -1,21 +1,5 @@
 @testset "Boundary Layer Functions" begin
-    # - Stagnation Point - #
-    cp_duct = abs.(range(-10, 10; step=1)) .+ 1.1
-    panel_lengths = ones(length(cp_duct))
-    s, s_stagnation, lower_length, upper_length = dt.split_at_stagnation_point(
-        panel_lengths, cp_duct
-    )
-    @test s == cumsum([0.0; panel_lengths[1:(end - 1)]])
-    @test s_stagnation == 10.0
-    @test lower_length == 10.0
-    @test upper_length == 10.0
-
-    # @testset "Radius of Curvature Calculation" begin
-    s = range(1, 10)
-    controlpoint = [s'; s' .^ 2]
-    ss = 4.0
-    R = dt.calculate_radius_of_curvature(s, controlpoint, ss)
-    @test isapprox(R, 1.0 / 0.003816453371976)
+    ##### ----- Head's Method Tests ----- #####
 
     # @testset "Step Setter" begin
     N = 10
@@ -32,93 +16,114 @@
     s = dt.arc_lengths_from_panel_lengths(panel_lengths)
 
     @test s == 0:9
+end
 
-    # @testset "Boundary Layer Initialization" begin
-    x = 1.0
-    Rex = 1.0
-    @test 0.036 == dt.d2_init(x, Rex)
+# these are for the green's method which may get resurrected at some point...
+# @testset "Boundary Layer Functions" begin
+#     # - Stagnation Point - #
+#     cp_duct = abs.(range(-10, 10; step=1)) .+ 1.1
+#     panel_lengths = ones(length(cp_duct))
+#     s, s_stagnation, lower_length, upper_length = dt.split_at_stagnation_point(
+#         panel_lengths, cp_duct
+#     )
+#     @test s == cumsum([0.0; panel_lengths[1:(end - 1)]])
+#     @test s_stagnation == 10.0
+#     @test lower_length == 10.0
+#     @test upper_length == 10.0
 
-    @test dt.calculate_H12bar0(2.0, 0.0) == 1.0 / (1.0 - 6.55)
+#     # @testset "Radius of Curvature Calculation" begin
+#     s = range(1, 10)
+#     controlpoint = [s'; s' .^ 2]
+#     ss = 4.0
+#     R = dt.calculate_radius_of_curvature(s, controlpoint, ss)
+#     @test isapprox(R, 1.0 / 0.003816453371976)
 
-    @test dt.calculate_CEeq(1.0, 0.0, 1.0, 1.0) == sqrt((1.0 - 0.32) / 1.2 + 0.0001) - 0.01
+#     # @testset "Boundary Layer Initialization" begin
+#     x = 1.0
+#     Rex = 1.0
+#     @test 0.036 == dt.d2_init(x, Rex)
 
-    s_init = 0.1
-    r_init = 1.0
-    Ue = 1.0
-    M = 0.1
-    rhoe = 1.0
-    mue = 1e-7
+#     @test dt.calculate_H12bar0(2.0, 0.0) == 1.0 / (1.0 - 6.55)
 
-    states, Cf_init, H12_init = dt.initialize_turbulent_boundary_layer_states(
-        s_init, r_init, Ue, M, rhoe, mue
-    )
+#     @test dt.calculate_CEeq(1.0, 0.0, 1.0, 1.0) == sqrt((1.0 - 0.32) / 1.2 + 0.0001) - 0.01
 
-    @test isapprox(
-        states, [0.00022714464401286951, 1.3836007395668004, 0.10980454099439806]
-    )
-    @test isapprox(Cf_init, 0.003581879845103732)
-    @test isapprox(H12_init, 1.3883679410459338)
+#     s_init = 0.1
+#     r_init = 1.0
+#     Ue = 1.0
+#     M = 0.1
+#     rhoe = 1.0
+#     mue = 1e-7
 
-    # @testset "Schlichting" begin
-    @test dt.d2_init(1.0, 1.0) == 0.036
+#     states, Cf_init, H12_init = dt.initialize_turbulent_boundary_layer_states(
+#         s_init, r_init, Ue, M, rhoe, mue
+#     )
 
-    # - H12bar_init - #
-    @test dt.H12bar_init(1.0, 1.0) == dt.calculate_H12bar0(1.0, 1.0)
+#     @test isapprox(
+#         states, [0.00022714464401286951, 1.3836007395668004, 0.10980454099439806]
+#     )
+#     @test isapprox(Cf_init, 0.003581879845103732)
+#     @test isapprox(H12_init, 1.3883679410459338)
 
-    # - CE_init - #
-    @test dt.CE_init(1.0, 0.1, 0.1) == dt.calculate_CEeq(1.0, 0.1, 1, 0.1)
+#     # @testset "Schlichting" begin
+#     @test dt.d2_init(1.0, 1.0) == 0.036
 
-    # - Fc - #
-    @test dt.Fc(0) == 1.0
+#     # - H12bar_init - #
+#     @test dt.H12bar_init(1.0, 1.0) == dt.calculate_H12bar0(1.0, 1.0)
 
-    # - FR - #
-    @test dt.FR(0) == 1.0
+#     # - CE_init - #
+#     @test dt.CE_init(1.0, 0.1, 0.1) == dt.calculate_CEeq(1.0, 0.1, 1, 0.1)
 
-    # - Reynolds - #
-    @test dt.calculate_Re(rhoe, Ue, 1.0, mue) == rhoe * Ue / mue
+#     # - Fc - #
+#     @test dt.Fc(0) == 1.0
 
-    # - H12bar0 - #
-    @test dt.calculate_H12bar0(0, 1.0) == 1.0
+#     # - FR - #
+#     @test dt.FR(0) == 1.0
 
-    # - Cf - #
-    @test isapprox(dt.calculate_Cf(2.2, 1.0, 1.0), 0.0, atol=eps())
+#     # - Reynolds - #
+#     @test dt.calculate_Re(rhoe, Ue, 1.0, mue) == rhoe * Ue / mue
 
-    # - H12 - #
-    @test dt.calculate_H12(0.0, sqrt(5)) == 1.0
+#     # - H12bar0 - #
+#     @test dt.calculate_H12bar0(0, 1.0) == 1.0
 
-    # - Ctau - #
-    @test dt.calculate_Ctau(1.0, 0.0, sqrt(10)) == 2.0 * (0.024 + 1.2)
+#     # - Cf - #
+#     @test isapprox(dt.calculate_Cf(2.2, 1.0, 1.0), 0.0, atol=eps())
 
-    # - F - #
-    @test dt.calculate_F(1.0, 1.0) == (0.02 + 1.0 + 0.8 * 1.0 / 3) / 1.01
+#     # - H12 - #
+#     @test dt.calculate_H12(0.0, sqrt(5)) == 1.0
 
-    # - H1 - #
-    @test dt.calculate_H1(2.0) == 3.15 + 1.72 - 0.01
+#     # - Ctau - #
+#     @test dt.calculate_Ctau(1.0, 0.0, sqrt(10)) == 2.0 * (0.024 + 1.2)
 
-    # - dH12bardH1 - #
-    @test dt.calculate_dH12bardH1(2.0) == -1.0 / (1.72 + 0.02)
+#     # - F - #
+#     @test dt.calculate_F(1.0, 1.0) == (0.02 + 1.0 + 0.8 * 1.0 / 3) / 1.01
 
-    # - Ri - #
-    @test dt.calculate_richardson_number(1.0, 1.0, 0.0, 1.0, 1.0) == 2.0 / 3.0 * (1.0 + 0.3)
+#     # - H1 - #
+#     @test dt.calculate_H1(2.0) == 3.15 + 1.72 - 0.01
 
-    # - Secondary Influences - #
-    Ri = 1.0
-    @test dt.longitudinal_curvature_influence(M, Ri) == 8.014
-    @test dt.lateral_strain_influence(ones(6)...) == -8.0 - 1 / 3
-    @test dt.dilation_influence(ones(7)...) == 10 + 1 / 3
+#     # - dH12bardH1 - #
+#     @test dt.calculate_dH12bardH1(2.0) == -1.0 / (1.72 + 0.02)
 
-    # - d2dUedsUeeq0 - #
-    @test dt.calculate_d2dUedsUeeq0(ones(4)...) == 0.625
+#     # - Ri - #
+#     @test dt.calculate_richardson_number(1.0, 1.0, 0.0, 1.0, 1.0) == 2.0 / 3.0 * (1.0 + 0.3)
 
-    # - CEeq0 - #
-    @test dt.calculate_CEeq0(ones(4)...) == 0.0
+#     # - Secondary Influences - #
+#     Ri = 1.0
+#     @test dt.longitudinal_curvature_influence(M, Ri) == 8.014
+#     @test dt.lateral_strain_influence(ones(6)...) == -8.0 - 1 / 3
+#     @test dt.dilation_influence(ones(7)...) == 10 + 1 / 3
 
-    # - Ctaueq0 - #
-    @test isapprox(dt.calculate_Ctaueq0(ones(3)...), 1.936)
+#     # - d2dUedsUeeq0 - #
+#     @test dt.calculate_d2dUedsUeeq0(ones(4)...) == 0.625
 
-    # - CEeq - #
-    @test isapprox(dt.calculate_CEeq(ones(4)...), 0.6907204085147591)
+#     # - CEeq0 - #
+#     @test dt.calculate_CEeq0(ones(4)...) == 0.0
 
-    # - d2dUedsUeeq - #
-    @test dt.calculate_d2dUedsUeeq(ones(4)...) == -0.25
-end
+#     # - Ctaueq0 - #
+#     @test isapprox(dt.calculate_Ctaueq0(ones(3)...), 1.936)
+
+#     # - CEeq - #
+#     @test isapprox(dt.calculate_CEeq(ones(4)...), 0.6907204085147591)
+
+#     # - d2dUedsUeeq - #
+#     @test dt.calculate_d2dUedsUeeq(ones(4)...) == -0.25
+# end