Skip to content

Commit ba923bd

Browse files
committed
docs: add EL bubble translational-motion equations; fix IB patch-types table
- equations.md 6.2: document the moving-bubble equation of motion -- tracer (dx/dt = u_l) and Newton's 2nd law (m dv/dt = F_D + F_p + F_g) -- with the non-dimensional drag/pressure/gravity force laws and their citations. This is the canonical @ref equations home; case.md keeps the concise inline version. - case.md: fix the Immersed Boundary Patch Types table -- the 3D Cylinder row carried a 4th cell with no matching header column, which broke the table render. Added a Notes column.
1 parent 9767418 commit ba923bd

2 files changed

Lines changed: 26 additions & 7 deletions

File tree

docs/documentation/case.md

Lines changed: 7 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -1320,14 +1320,14 @@ Boundary is at polar angle \f$\theta = \mathrm{atan2}(y - y_{\mathrm{centroid}},
13201320

13211321
### Immersed Boundary Patch Types {#immersed-boundary-patch-types}
13221322

1323-
| # | Name | Dim. |
1324-
| ---: | :----: | :--- |
1325-
| 2 | 2D Circle | 2 |
1326-
| 3 | 2D Rectangle | 2 |
1327-
| 4 | 2D Airfoil | 2 |
1328-
| 8 | 3D Sphere | 3 |
1323+
| # | Name | Dim. | Notes |
1324+
| ---: | :----: | :---: | :--- |
1325+
| 2 | 2D Circle | 2 | |
1326+
| 3 | 2D Rectangle | 2 | |
1327+
| 4 | 2D Airfoil | 2 | |
1328+
| 8 | 3D Sphere | 3 | |
13291329
| 10 | 3D Cylinder | 3 | `length_x` sets the axial length of the cylinder. |
1330-
| 11 | 3D Airfoil | 3 |
1330+
| 11 | 3D Airfoil | 3 | |
13311331

13321332
### Acoustic Supports {#acoustic-supports}
13331333

docs/documentation/equations.md

Lines changed: 19 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -522,6 +522,25 @@ with \f$\sigma = \varepsilon_b \max(\Delta x^{1/3}_\text{cell},\;R_\text{bubble}
522522

523523
Each bubble is tracked individually with Keller-Miksis dynamics and 4th-order adaptive Runge-Kutta time integration.
524524

525+
**Translational motion (`vel_model > 0`):**
526+
527+
Bubbles may also translate through the carrier flow. Let \f$\mathbf{x}_b\f$ be the bubble position, \f$\mathbf{u}_b\f$ its velocity, \f$a\f$ its radius, and \f$\mathbf{u}_l(\mathbf{x}_b)\f$ the carrier velocity interpolated to the bubble location (a Lagrange polynomial of order `fd_order`, which must be set when `vel_model > 0`).
528+
529+
- **Tracer bubbles (`vel_model = 1`)** follow the local carrier velocity:
530+
\f[\frac{d\mathbf{x}_b}{dt} = \mathbf{u}_l(\mathbf{x}_b).\f]
531+
532+
- **Newton's second law (`vel_model = 2`)** integrates the bubble momentum:
533+
\f[m_b\,\frac{d\mathbf{u}_b}{dt} = \mathbf{F}_D + \mathbf{F}_p + \mathbf{F}_g, \qquad \frac{d\mathbf{x}_b}{dt} = \mathbf{u}_b,\f]
534+
with bubble mass \f$m_b\f$ and forces acting on the slip velocity \f$\mathbf{u}_\text{rel} = \mathbf{u}_b - \mathbf{u}_l\f$:
535+
536+
| Force | `case.py` control | Non-dimensional form |
537+
|---|---|---|
538+
| Drag \f$\mathbf{F}_D\f$ | `drag_model` | \f$-\,c_D\,\pi\,a\,\mathbf{u}_\text{rel} / \text{Re}\f$, with \f$c_D = 4\f$ free-slip (\cite Hadamard1911; \cite Rybczynski1911), \f$c_D = 6\f$ no-slip Stokes (\cite Stokes1851), \f$c_D = 12\f$ Levich (\cite Levich1962) |
539+
| Pressure \f$\mathbf{F}_p\f$ | `pressure_force` | \f$-V_b\,\nabla p\f$, with bubble volume \f$V_b = \frac{4}{3}\pi a^3\f$ |
540+
| Gravity \f$\mathbf{F}_g\f$ | `gravity_force` | \f$m_b\,\mathbf{g}\f$, with \f$\mathbf{g}\f$ the body-force acceleration |
541+
542+
Here \f$\text{Re}\f$ is the mixture Reynolds number — the same `fluid_pp%%Re(1)` that scales the viscous stress tensor (@ref sec-two-viscosities) — so the drag scales with the liquid viscosity. The three drag models increase in magnitude free-slip \f$<\f$ no-slip \f$<\f$ Levich; see \cite Magnaudet2000 for a review of these bubble-drag regimes.
543+
525544
---
526545

527546
## 7. Fluid-Structure Interaction

0 commit comments

Comments
 (0)