fix minor bug in drag

byuflowlab · Feb 14, 2025 · e53d0d3 · e53d0d3
1 parent 30505c2
commit e53d0d3
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 2 deletions.
diff --git a/src/analysis/analyses.jl b/src/analysis/analyses.jl
@@ -213,6 +213,7 @@ function analyze(
             linsys=(; solve_parameter_tuple.linsys..., A_bb_LU),
             idmaps,
             problem_dimensions,
+            options = deepcopy(options),
         ),
         any(options.solver_options.converged)
     else

diff --git a/src/postprocess/viscous_drag.jl b/src/postprocess/viscous_drag.jl
@@ -156,11 +156,17 @@ function compute_single_side_drag_coefficient(
     # printdebug("Hsep:", Hsep)
     # printdebug("δ_2:", usep[1])
 
-    cdsq = squire_young(
+    cdsqy = squire_young(
         usep[1], boundary_layer_functions.edge_velocity(s_sep), Vref[], Hsep
     )
 
-    cd = FLOWMath.ksmin([single_side_boundary_layer_options.separation_penalty; cdsq])
+    if single_side_boundary_layer_options.separation_penalty < eps()
+        cd = cdsqy
+    else
+        cd = FLOWMath.ksmin(
+            [single_side_boundary_layer_options.separation_penalty; cdsqy], 100
+        )
+    end
 
     cdadd = FLOWMath.ksmax(
         [