Notes on diffusion-reaction in a dermograph: steady-state 1D model, boundary conditions, and Thiele modulus

Overview

• The instructor begins with a casual analogy: if you don’t like your Quiz score, you have other attempts; the numbers change between attempts, highlighting that parameters can vary between runs. This sets up the idea that model parameters (boundary conditions, diffusion constants, etc.) are not fixed constants across all problems.
• Then the lecture moves to diffusion in the context of a diffusion-reaction problem: adding a reaction term to the diffusion equation to model consumption of oxygen by cells inside a scaffold.
• Core idea: diffusion is the process by which molecules move from regions of high concentration to low concentration. Diffusion is an ensemble (average) effect, not necessarily the path of any single molecule.
• Fick’s Law: the diffusive flux is proportional to the concentration gradient. The diffusion constant D links flux to the gradient via flux J = -D ∂C/∂x (in one dimension).
• The diffusion constant D depends on temperature and viscosity, via the Stokes–Einstein relation: D is larger at higher temperature and smaller in more viscous solutions. Temperature changes are commonly reported with diffusion constants; viscosity slows diffusion.
• Intuition from the diffusion constant: heating the medium increases molecular motion, hence faster diffusion; higher viscosity slows diffusion.
• In a practical model (the “thermograph” or dermograph), we consider a slab of hydrogel with cells embedded, surrounded by a liquid bath containing nutrients, salts, sugars, and oxygen.
• Geometry: the problem is treated as a 3D slab, but diffusion is taken as the dominant direction (through the thickness) so the problem is reduced to 1D along x (through the thickness from the media into the scaffold).
• Boundary conditions motivation:
  • The diffusion from the media into the scaffold dominates; diffusion from the other edges is neglected.
  • Thickness L is the relevant length scale; x runs from 0 at one boundary to L at the boundary in contact with the media.

Governing equations and physical setup

• Mass balance approach yields a diffusion–reaction equation. In one dimension, with a reaction term R, the governing equation is: $\frac{\partial C}{\partial t} = D \frac{\partial^2 C}{\partial x^2} + R.$
  • D is the diffusion coefficient; C is the oxygen concentration; R is the reaction term (negative for consumption).
• In this context, cells within the scaffold consume oxygen via aerobic respiration (sugar + oxygen --> CO₂ + water). This links to a redox chemistry interpretation where a net consumption rate R is introduced in the mass balance.
• We can also view the CO₂ production as a parallel transport consideration (diffusing out of the scaffold) but the primary transport problem here is oxygen diffusion into the scaffold and consumption by cells.
• Units consistency note: the reaction term R must have the same units as ∂C/∂t, i.e. concentration per time, so that the equation is dimensionally consistent. If D has units of
  $\mathrm{m^2\,s^{-1}}$ and C has units of $\mathrm{mol\,m^{-3}}$, then R has units of $\mathrm{mol\,m^{-3}\,s^{-1}}$.
• For this lecture, we adopt a common simplifying assumption for the kinetics: the reaction rate R is treated as a constant (zero-order in C) for the purpose of solving the equation. Other kinetics (first-order: R ∝ C, or more complex dependencies) are mentioned as possibilities, but the constant-R assumption is used for the derivation here.

Steady-state simplification

• Steady-state assumption (often a good engineering approximation when you are interested in the long-time concentration profile after transients settle):
  $\frac{\partial C}{\partial t} = 0$
• Under steady state, the diffusion–reaction equation reduces to: $D \frac{d^2 C}{d x^2} + R = 0.$
  • This is a second-order ordinary differential equation in x.
  • The solution will be a quadratic in x once integrated twice.

Boundary conditions and physical reasoning

• Boundary condition 1 (edge in contact with media): the concentration at the outer boundary is the media concentration, i.e.,
  $C(L) = C<em>0,$ where $C0$ is the concentration of oxygen in the surrounding media.
• Boundary condition 2 (no flux at the inner “wall” of the scaffold): the wall is assumed not to be permeable to oxygen, so there is no flux through that boundary: \left.\frac{dC}{dx}\right|_{x=0} = 0.$n- Rationale for boundary conditions:
  • The left boundary (x = 0) is a wall with no oxygen permeation, hence zero flux.
  • The right boundary (x = L) is in contact with the well-mixed media with fixed oxygen concentration C₀.
  • The media is assumed to be well stirred so that the oxygen concentration in the bulk media is uniform, ensuring C(L) = C0 is a good boundary condition.
• Alternative geometry comment (symmetry): there is a discussion about modeling the problem with a geometry where the scaffold sits on a boundary and diffusion is symmetric about the midplane. In that symmetric setup, one boundary condition could be C(0) = 0 with symmetry implying dC/dx|_{x=0} = 0 in the mirrored formulation, yielding the same solution shape when properly transformed. The key point is that the boundary conditions and geometry can be framed in multiple equivalent ways, but the resulting concentration profile (and the role of the Thiele modulus) remains the same.

Solve for the concentration profile

• General integration of the steady-state equation: D \frac{d^2 C}{d x^2} + R = 0 \quad\Rightarrow\quad \frac{d^2 C}{d x^2} = -\frac{R}{D}.
  • Integrate once: \frac{dC}{dx} = -\frac{R}{D} x + C_1.
  • Integrate again: C(x) = -\frac{R}{2D} x^2 + C1 x + C2.
• Apply boundary conditions:
  • From no-flux at x = 0: \left.\frac{dC}{dx}\right|{x=0} = -\frac{R}{D} \cdot 0 + C1 = 0 \Rightarrow C_1 = 0.
  • From C(L) = C0: C0 = -\frac{R}{2D} L^2 + 0 \cdot L + C2 \Rightarrow C2 = C0 + \frac{R}{2D} L^2.
• Substituting back: the concentration profile in dimensional form is
  \boxed{\;C(x) = C_0 + \frac{R}{2D}\left( L^2 - x^2 \right)\;}.
• Interpretation: since R is negative for consumption, the term (\frac{R}{2D}(L^2 - x^2)) reduces C(x) below C0 as x moves away from the boundary in contact with the media.

Dimensionless form and the Thiele modulus

• Define dimensionless quantities:
  • Dimensionless concentration: \theta(\xi) = \frac{C(x)}{C_0}$where$\xi = \frac{x}{L}.
• Define the dimensionless parameter (Thiele modulus, also referred to here as the Heu–/modulus in the lecture; commonly called the Thiele modulus): \phi = \frac{R L^2}{2 D C_0}.
  • Note: R is negative for consumption, so φ is typically negative; often magnitude |φ| is used to discuss the strength of consumption relative to diffusion.
• Rewriting the solution in dimensionless form:
  • Start from the dimensional solution and divide by C0:
    \theta(\xi) = \frac{C(x)}{C0} = 1 + \frac{R L^2}{2 D C0} \left( 1 - \xi^2 \right) = 1 + \phi \left( 1 - \xi^2 \right).
  • Since (\xi = x/L) and (0 \le \xi \le 1), this is
    \boxed{\;\theta(\xi) = 1 + \phi(1 - \xi^2)\;}.
• Physical interpretation of φ (the Thiele modulus):
  • If φ is large in magnitude (strong consumption relative to diffusion), the concentration drops more quickly with x, and the interior of the scaffold experiences lower oxygen levels.
  • If φ is small in magnitude (diffusion dominates or reaction is weak), the concentration stays closer to the boundary value C0 throughout most of the scaffold.
• Relationships to design knobs:
  • To reduce depletion (keep O2 higher inside the scaffold), you can increase D (improve diffusion), increase C0 (raise media oxygen concentration), or decrease L (make the scaffold thinner).
  • The Thiele modulus provides a single dimensionless criterion to compare these competing effects; it helps determine whether diffusion or reaction dominates the transport process.
• Additional non-dimensionalization note (planning ahead without solving): one can non-dimensionalize the equation before integrating by setting
  • Dimensionless concentration: \Theta = \frac{C}{C_0}
  • Dimensionless distance: \xi = \frac{x}{L}
  • Then the dimensionless equation becomes a form that reveals that the transport is controlled by the ratio of reaction strength to diffusion and the geometry (L). This approach often suggests which terms can be neglected in certain limits (e.g., small φ or large φ).

Physical interpretation and practical implications

• Effect of a large Thiele modulus (|φ| large):
  • Strong consumption relative to diffusion.
  • The concentration declines more steeply into the scaffold; interior cells see much lower oxygen levels.
• Effect of a small Thiele modulus (|φ| small):
  • Diffusion dominates; the O₂ concentration remains closer to C0 throughout the scaffold.
• Design criteria discussion:
  • In lab/engineering design, you might set a performance criterion (e.g., ensure C(x) remains above a threshold Cmin at a certain depth within the scaffold).
  • You can then decide which levers to tune: increase diffusion constant D (choice of materials, porosity, temperature), reduce thickness L, or increase the initial media concentration C0. Since D is a physical property, changing it may be limited by material choice; the more controllable knob is often the external oxygenation (C0) or geometry (L).
• Conceptual insight: the dimensionless Thiele modulus captures the balance between two competing processes (diffusion vs consumption) and provides a compact way to reason about system behavior without solving the full equation each time.

Alternative modeling perspective and symmetry note

• The same 1D solution can arise from a different geometric setup that is symmetric about a mid-plane: if you model diffusion with a left boundary that appears as a wall (no flux) and a right boundary in contact with media, you can flip the geometry and obtain an equivalent concentration profile. In that symmetric view, one boundary condition could be C(0) = 0 with a zero-flux condition at the symmetry plane; the exact solution in terms of C(x) vs x/L looks the same after appropriate reflection. The key point remains the boundary conditions and the dimensionless modulus govern the shape of the profile, not the specific orientation of the slab.

Take-home messages

• A steady-state diffusion with a constant consumption rate yields a quadratic C(x) profile with a simple closed-form solution given the BCs.
• Dimensional form: \boxed{C(x) = C_0 + \frac{R}{2D}\left(L^2 - x^2\right)} with R < 0 for consumption.
• Dimensionless form highlights the Thiele modulus: \boxed{\theta(\xi) = 1 + \phi(1 - \xi^2),\quad \phi = \frac{R L^2}{2 D C0}}$where$\xi = \frac{x}{L}$and$\theta = \frac{C}{C0}$$.
• The thickness of the scaffold, the diffusion coefficient, the boundary oxygen concentration, and the consumption rate all combine to determine whether oxygen is sufficiently available throughout the scaffold; this is succinctly captured by the Thiele modulus.
• Non-dimensionalization can be used proactively to identify dominant terms and to guide experimental design or material choices without performing a full solution.