rework boundary layer

byuflowlab · Jan 12, 2025 · 1d59939 · 1d59939
1 parent e5781ef
commit 1d59939
Show file tree

Hide file tree

Showing 20 changed files with 1,658 additions and 434 deletions.
diff --git a/Project.toml b/Project.toml
@@ -17,6 +17,7 @@ LineSearches = "d3d80556-e9d4-5f37-9878-2ab0fcc64255"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 MINPACK = "4854310b-de5a-5eb6-a2a5-c1dee2bd17f9"
 NLsolve = "2774e3e8-f4cf-5e23-947b-6d7e65073b56"
+OrdinaryDiffEq = "1dea7af3-3e70-54e6-95c3-0bf5283fa5ed"
 PreallocationTools = "d236fae5-4411-538c-8e31-a6e3d9e00b46"
 Printf = "de0858da-6303-5e67-8744-51eddeeeb8d7"
 ProgressMeter = "92933f4c-e287-5a05-a399-4b506db050ca"

diff --git a/src/DuctAPE.jl b/src/DuctAPE.jl
@@ -44,6 +44,9 @@ using ForwardDiff # used for jacobian for newton solver
 
 # For boundary layer stuff
 using Roots
+using OrdinaryDiffEq:
+    ODEProblem, terminate!, ContinuousCallback, RadauIIA5
+using OrdinaryDiffEq: solve as diffeq_solve
 
 # - Utility Packages - #
 using FLOWMath # used for various items, mostly interpolation
@@ -85,6 +88,8 @@ export ChainSolverOptions,
     FixedPointOptions,
     CSORSolverOptions
 export BoundaryLayerOptions
+export RK2, RK4
+export RadauIIA5, RadauIIA3
 
 # - Preprocess - #
 export setup_analysis
@@ -190,6 +195,7 @@ include("postprocess/viscous_drag.jl")
 # include("visualization/plot_recipe_defaults.jl")
 include("visualization/plot_recipes.jl")
 include("visualization/calculate_streamlines.jl")
+    include("visualization/calculate_velocity_field.jl")
 include("visualization/convenience_plots.jl")
 
 ##### ----- DEBUGGING ----- #####

diff --git a/src/analysis/analyses.jl b/src/analysis/analyses.jl
@@ -66,7 +66,7 @@ function analyze(
         #TODO: probably just call  the post-process function directly and return a reset_container! of the output
         if return_inputs
             return [],#zero_outputs(),
-            (; solve_parameter_tuple..., airfoils, idmaps, panels, problem_dimensions),
+            (; solve_parameter_tuple..., airfoils, idmaps, panels, problem_dimensions, reference_parameters, multi_point = [operating_point]),
             false
         else
             return [],#zero_outputs(),
@@ -410,6 +410,8 @@ function analyze(
             linsys=(; solve_parameter_tuple.linsys..., A_bb_LU),
             idmaps,
             problem_dimensions,
+            reference_parameters,
+            multi_point=operating_point,
         ),
         options.solver_options.converged
     else