diff --git a/.pre-commit-config.yaml b/.pre-commit-config.yaml
index b3a59792..80231586 100644
--- a/.pre-commit-config.yaml
+++ b/.pre-commit-config.yaml
@@ -28,7 +28,7 @@ repos:
   - id: isort
     name: isort
 - repo: https://github.com/psf/black
-  rev: 23.12.1
+  rev: 24.1.1
   hooks:
   - id: black
 - repo: https://github.com/pycqa/flake8
diff --git a/src/planingfsi/__init__.py b/src/planingfsi/__init__.py
index f7b2a14f..b54327b2 100644
--- a/src/planingfsi/__init__.py
+++ b/src/planingfsi/__init__.py
@@ -8,6 +8,7 @@
 and the structural solver considers a large-deformation simple beam element.
 
 """
+
 import logging
 
 logging.basicConfig()
diff --git a/src/planingfsi/__main__.py b/src/planingfsi/__main__.py
index b93348c9..240134c1 100644
--- a/src/planingfsi/__main__.py
+++ b/src/planingfsi/__main__.py
@@ -1,4 +1,5 @@
 """Package-level accessor to run the main CLI program with python -m planingfsi."""
+
 from planingfsi.cli import cli
 
 if __name__ == "__main__":
diff --git a/src/planingfsi/cli.py b/src/planingfsi/cli.py
index 13e09e8c..c63b4f4e 100644
--- a/src/planingfsi/cli.py
+++ b/src/planingfsi/cli.py
@@ -5,6 +5,7 @@
 files, assembles the problem, and runs it.
 
 """
+
 from __future__ import annotations
 
 from pathlib import Path
diff --git a/src/planingfsi/config.py b/src/planingfsi/config.py
index cac9a721..2c02d9f5 100644
--- a/src/planingfsi/config.py
+++ b/src/planingfsi/config.py
@@ -7,6 +7,7 @@
 instance, which then serves as a reference point elsewhere in the code.
 
 """
+
 from __future__ import annotations
 
 import math
diff --git a/src/planingfsi/fe/femesh.py b/src/planingfsi/fe/femesh.py
index 13225024..b963a7fc 100644
--- a/src/planingfsi/fe/femesh.py
+++ b/src/planingfsi/fe/femesh.py
@@ -1,4 +1,5 @@
 """The `femesh` module handles mesh generation for planingFSI cases."""
+
 from __future__ import annotations
 
 import abc
diff --git a/src/planingfsi/math_helpers.py b/src/planingfsi/math_helpers.py
index 6d1bb9ac..03fe204c 100644
--- a/src/planingfsi/math_helpers.py
+++ b/src/planingfsi/math_helpers.py
@@ -1,4 +1,5 @@
 """General math helpers."""
+
 from __future__ import annotations
 
 from collections.abc import Callable
diff --git a/src/planingfsi/potentialflow/pressureelement.py b/src/planingfsi/potentialflow/pressureelement.py
index d0cb600b..c12949dd 100644
--- a/src/planingfsi/potentialflow/pressureelement.py
+++ b/src/planingfsi/potentialflow/pressureelement.py
@@ -1,4 +1,5 @@
 """Module containing definitions of different types of pressure element."""
+
 from __future__ import annotations
 
 import abc
diff --git a/src/planingfsi/potentialflow/pressurepatch.py b/src/planingfsi/potentialflow/pressurepatch.py
index bb9e77cd..7a2f6990 100644
--- a/src/planingfsi/potentialflow/pressurepatch.py
+++ b/src/planingfsi/potentialflow/pressurepatch.py
@@ -1,4 +1,5 @@
 """Classes representing a pressure patch on the free surface."""
+
 from __future__ import annotations
 
 import abc
@@ -653,9 +654,7 @@ def _calculate_shear_stress(self) -> None:
         def get_shear_stress(xx: float) -> float:
             """Calculate the shear stress at a given location."""
             re_x = self.config.flow.flow_speed * xx / self.config.flow.kinematic_viscosity
-            return (
-                0.332 * self.config.flow.density * self.config.flow.flow_speed**2 * re_x**-0.5
-            )
+            return 0.332 * self.config.flow.density * self.config.flow.flow_speed**2 * re_x**-0.5
 
         x = self.element_coords[:-1]
         assert self.interpolator is not None
diff --git a/src/planingfsi/potentialflow/solver.py b/src/planingfsi/potentialflow/solver.py
index 1ab16c52..bae18d8d 100644
--- a/src/planingfsi/potentialflow/solver.py
+++ b/src/planingfsi/potentialflow/solver.py
@@ -1,4 +1,5 @@
 """Fundamental module for constructing and solving planing potential flow problems."""
+
 from __future__ import annotations
 
 from pathlib import Path
diff --git a/src/planingfsi/simulation.py b/src/planingfsi/simulation.py
index 9cfd1e0f..5cf0f5f1 100644
--- a/src/planingfsi/simulation.py
+++ b/src/planingfsi/simulation.py
@@ -1,4 +1,5 @@
 """High-level control of a `planingfsi` simulation."""
+
 from __future__ import annotations
 
 from pathlib import Path
diff --git a/src/planingfsi/solver.py b/src/planingfsi/solver.py
index 6450484d..7c8e51ad 100644
--- a/src/planingfsi/solver.py
+++ b/src/planingfsi/solver.py
@@ -1,4 +1,5 @@
 """Solvers for single and multi-dimensional nonlinear problems."""
+
 from __future__ import annotations
 
 from collections.abc import Callable
diff --git a/src/planingfsi/trig.py b/src/planingfsi/trig.py
index d4a29695..151fe821 100644
--- a/src/planingfsi/trig.py
+++ b/src/planingfsi/trig.py
@@ -1,4 +1,5 @@
 """Convenient trigonometric functions."""
+
 import math
 from typing import Union
 
diff --git a/src/planingfsi/unit.py b/src/planingfsi/unit.py
index b50b8f0a..9e455dc5 100644
--- a/src/planingfsi/unit.py
+++ b/src/planingfsi/unit.py
@@ -4,6 +4,7 @@
     which will return 0.3048 meters.
 
 """
+
 import math
 
 # Acceleration due to gravity
diff --git a/src/planingfsi/writers.py b/src/planingfsi/writers.py
index dc80ab39..4c757583 100644
--- a/src/planingfsi/writers.py
+++ b/src/planingfsi/writers.py
@@ -1,4 +1,5 @@
 """Functions for writing results to a file on disk."""
+
 from __future__ import annotations
 
 from pathlib import Path