---
title: 'Cassie: Wetting State & Bipedal Robot'
url: https://www.emergentmind.com/topics/cassie
type: topic
---

# Cassie: Wetting State & Bipedal Robot

“Cassie” denotes two distinct technical objects in contemporary research. In interfacial science, the term refers to the Cassie or Cassie–Baxter wetting state, a composite state in which a liquid sits on top of a textured surface while vapor or air remains trapped in the underlying corrugations; this state underlies much of the modern literature on superhydrophobicity, Cassie–Wenzel transitions, and transport over textured interfaces [1810.02104]. In robotics, Cassie is a 20-DOF, underactuated, human-scale biped developed by Agility Robotics, which has become a common experimental platform for research on hybrid locomotion, trajectory optimization, feedback control, reinforcement learning, wheeled riding, and expressive motion [1907.11353].

## 1. Cassie as a wetting state

In the wetting literature, the Cassie state is the composite-interface state in which the liquid rests on top of pillars, ridges, or other asperities, with vapor trapped below; the Wenzel state is the fully intruded state in which the liquid completely wets the corrugations [1810.02104]. On a rough composite surface with solid fraction \(f_s\) and intrinsic Young angle \(\theta_Y\), the classical Cassie–Baxter relation is

\[
\cos\theta^* = f_s(\cos\theta_Y + 1) - 1,
\]

whereas the Wenzel relation is

\[
\cos\theta_W = r\cos\theta_Y,
\]

with \(r\) the roughness ratio [2111.08499, 2605.07688]. In this formulation, the Cassie state is associated with large apparent contact angles, low adhesion, and enhanced slip, while the Wenzel state corresponds to the loss of the composite interface and the disappearance of most superhydrophobic properties [1810.02104].

The classical Cassie–Baxter equation is explicitly static. It assumes equilibrium, a well-defined solid fraction, and no penetration into the asperities, and it does not describe spreading, recoil, infiltration, or contact-line dissipation [2111.08499]. This limitation becomes central once wetting is treated as a transition problem rather than a static geometry problem.

A rigorous two-scale homogenization treatment modifies the classical picture. For purely chemically inhomogeneous surfaces, the modified Cassie law becomes

\[
\cos\theta_a = \int_0^1 \cos\bigl(\theta_s(Y)\bigr)\,dY,
\]

which is a line-average of the local \(\cos\theta_s\) along the contact line rather than an area-average over the substrate [1609.07628]. In that framework, the modified equations correspond to local minimizers of the total interface energy and thereby account for hysteresis, whereas the classical area-averaged laws correspond to a global-minimizer idealization [1609.07628].

## 2. Free-energy landscapes and the Cassie–Wenzel transition

A central result of the modern Cassie literature is that the Cassie–Wenzel transition is not determined solely by a single Laplace-pressure threshold. For intrusion and extrusion in a single truncated or inverted truncated conical pore, the macroscopic capillary-thermodynamics description uses the grand potential

\[
G = F - \Delta p\,V,
\qquad
F = \gamma_{lv}S_{lv} - \gamma_{lv}\cos\theta_Y\,S_{sl},
\]

with \(\Delta p = p_{\text{liquid}} - p_{\text{vapor}}\) and filling fraction \(\zeta=z/H\in[0,1]\) [2010.11532]. For the truncated-cone geometry, the excess grand potential relative to the dry Cassie state has the cubic form

\[
\Delta G(\zeta)=\pi R_0^2\gamma_{lv}\,\Delta g(\zeta),
\qquad
\Delta g(\zeta)=g_1\zeta+g_2\zeta^2+g_3\zeta^3,
\]

so the transition is governed by the full one-dimensional free-energy landscape rather than by a local criterion alone [2010.11532].

The usual local-stability condition at \(\zeta=0\) yields a critical pressure

\[
\Delta p_c = -\frac{2\gamma_{lv}}{R_0}\left[a+\cos\theta_Y\sqrt{1+a^2}\right],
\]

which reduces for a cylinder \((a=0)\) to the familiar \(\Delta p_c=-2\gamma_{lv}\cos\theta_Y/R_0\) [2010.11532]. However, this criterion identifies only the loss of local stability of the dry state. It does not determine the location of a new interior minimum, the remaining barrier to full filling, or the reversibility of the pathway [2010.11532].

The distinction between gradual and abrupt collapse follows directly from the geometry of the free-energy profile. For a narrowing pore \((\alpha>0)\) with a hydrophobic wall, when \(\tilde p\) rises above \(\tilde p_c\) the meniscus slides gradually from \(\zeta=0\) to a finite-depth minimum and only at a larger pressure \(\tilde p_s\) does the pore fill completely. For an inverted widening pore \((\alpha<0)\), the Cassie state loses local stability abruptly at \(\tilde p=\tilde p_c\) and the meniscus jumps to the bottom at a spinodal [2010.11532]. The same model also predicts metastable minima under both compression and stretching, together with irreversibility driven by bottom-wall adsorption [2010.11532].

A related capillary result for spherical cavities shows that the free-energy barrier from complete drying to complete wetting is attained when the meniscus becomes flat; the free-energy maximum is given by an analytic formula, and the effect of line tension is largest at that maximum [1608.08286]. Positive line tension raises the barrier and stabilizes the Cassie state, whereas negative line tension lowers the barrier and destabilizes it [1608.08286].

## 3. Microscopic, dynamic, and non-equilibrium extensions

At nanometric scales, atomistic and mesoscale calculations replace the single filling coordinate by field-based descriptions. In a submerged hydrophobic cavity, the atomistic string method in the coarse-grained density field finds a most probable wetting path in which a meniscus initially pinned at the top corners bends and slides downward, then breaks symmetry near the bottom by nucleating a liquid finger on one side, leaving a vapor bubble in the opposite corner before final collapse to Wenzel [1411.4683]. In the Lennard-Jones cases reported there, the free-energy barrier is \(\Delta F^\ddagger_{\text{square}}\simeq 34\,k_BT\) for a square groove and \(\Delta F^\ddagger_{\text{rect}}\simeq 44\,k_BT\) for a rectangular groove, with the larger barrier associated with the larger aspect ratio [1411.4683].

The choice of collective variables materially changes the computed mechanism. On pillared substrates, using only the overall liquid density in the interpillar volume produces discontinuous, non-physical meniscus evolution, spurious kinks in the free-energy profile, and substantial errors in both state stability and barrier height; in a \(3\times 3\) pillared system at \(\Delta P\approx 0\), fine-grid CVs on an 864-cell grid give \(\Delta F_{CW}\approx 940\,k_BT\) and \(\Delta F^\ddagger\approx 420\,k_BT\), whereas a single-CV description shifts these values to \(\Delta F_{CW}\approx 1\,340\,k_BT\) and \(\Delta F^\ddagger\approx 380\,k_BT\) [1810.02104]. Finite-size effects are similarly strong: a \(1\times 1\) periodic-pillar system overestimates the barrier by \(\sim 180\)–\(200\,k_BT\) and suppresses the translational-symmetry breaking of the meniscus [1810.02104].

Microscopic DFT studies of nanoscale corrugation show that macroscopic force-balance remains qualitatively useful, but only if the Young angle is replaced by an effective wall angle \(\theta_{\text{eff}}\) that encodes finite interface thickness, finite-range solid–liquid forces, and line-tension corrections [1607.01636]. For nanostripes and square nanopillars, \(\theta_{\text{eff}}\) converges to a plateau about \(5^\circ\)–\(20^\circ\) below the corresponding Young angle even for corrugation periods of order \(30\,\text{nm}\), implying impalement pressures smaller than predicted by a purely macroscopic theory [1607.01636]. For binary liquid mixtures on nano-pitted surfaces, the transition cannot be predicted from a single contact angle parameter; composition, \(A\)–\(B\) interaction strength, wall affinities, and pit dimensions all shift both the intrusion threshold and the Cassie/Wenzel free-energy difference [1910.13766].

The static Cassie–Baxter average also fails for dynamic imbibition. Starting from the Lucas–Washburn description with spatially varying wettability, the crossing time depends on \(1/\cos\theta_i\) segment by segment, leading in the many-segment limit to the effective contact-angle relation

\[
\cos\theta_{\rm eff}
=
\left(
r_A\frac{1}{\cos\theta_A}
+
r_B\frac{1}{\cos\theta_B}
\right)^{-1},
\]

rather than the arithmetic average of the classical Cassie–Baxter law [2501.10255]. Ordered arrangements are likewise nontrivial: the “more hydrophobic-first” configuration minimizes the crossing time [2501.10255].

Non-equilibrium thermal effects provide an additional route to Cassie stabilization. On supercooled square-pillar silicon surfaces with \(w=s=100\,\mu\text{m}\), \(h=500\,\mu\text{m}\), and \(\theta_Y=120^\circ\), localized freezing at the droplet bottom can arrest penetration and produce a final Wenzel-to-Cassie transition after impact [2605.07688]. The regime map reported there places the Cassie regime at low \(We\) and high \(Ste\), specifically low \(We\lesssim 7\) and high \(Ste\gtrsim 0.27\), with partial Wenzel and full Wenzel occupying the complementary regions [2605.07688]. This suggests that freezing kinetics can act as a control parameter alongside geometry and wettability.

## 4. Hydrodynamics and heat transfer in the Cassie state

The Cassie state is not only a wetting morphology but also a transport regime. For shear flow over periodic rectangular grooves filled with an immiscible secondary fluid, eliminating the subphase dynamics yields a locally varying Navier slip condition on the meniscus,

\[
u_{\text{primary}}=\ell(x)\,\frac{\partial u_{\text{primary}}}{\partial n},
\]

with \(\ell(x)\) determined by groove geometry and the viscosity ratio \(N=\mu_1/\mu_2\) [1311.7006]. Closed-form expressions for the effective slip length show that the surface response is anisotropic: longitudinal flow experiences larger local and effective slip than transverse flow because the trapped-fluid vortex structure differs between the two orientations [1311.7006]. The analysis further shows that for water–air systems \((N\approx 55)\), groove aspect ratios \(A\approx 0.2\)–\(0.3\) already recover more than \(90\%\) of the maximal slip, so very deep grooves are unnecessary [1311.7006].

For laminar, fully developed internal flow through Cassie-state microchannels textured with isoflux ridges, matched asymptotics yields both hydrodynamic slip lengths and Nusselt numbers in terms of the small parameter \(\epsilon=d/H\ll 1\) and solid fraction \(\phi_s\) [2211.13751]. With ridges parallel to the flow, the dimensionless slip length satisfies

\[
\tilde b = \epsilon\,\lambda(\phi_s) + O(e^{-\pi/\epsilon}),
\qquad
\lambda(\phi_s)=\frac{2}{\pi}\ln\!\left[\sec\!\left(\frac{\pi\delta}{2}\right)\right],
\]

and the new closed-form Nusselt result is accurate to \(O(\epsilon^2)\) for any solid fraction, correcting earlier \(O(\epsilon^1)\) expressions that break down as \(\phi_s\to 0\) [2211.13751]. With ridges transverse to the flow, the leading inertial correction to the slip length is \(O(\epsilon^3 Re)\), and the Nusselt-number accuracy depends additionally on \(Re\), \(Pe\), and \(Pr\) [2211.13751]. The inner thermal problem can be summed in terms of polylogarithms, linking Cassie-state microchannel heat transfer to classical spreading-resistance theory [2211.13751].

## 5. Cassie as a bipedal robot platform

In robotics, Cassie is a human-scale biped developed by Agility Robotics. In the rigid-body model used for much of the locomotion literature, it is a 20-DOF, underactuated robot with a 6-DOF unactuated pelvis and two 7-DOF legs, each leg having five actuated joints and two passive leaf-spring joints [1809.07279]. In the hovershoe study, the leg coordinates are identified as hip roll, hip yaw, hip pitch, knee pitch, shin pitch, tarsus pitch, and toe pitch, with the pelvis coordinates collected in the generalized vector

\[
q=[q_x,q_y,q_z,q_{\text{yaw}},q_{\text{pitch}},q_{\text{roll}},q_{1L}\ldots q_{7L},q_{1R}\ldots q_{7R}]^T
\]

[1907.11353]. A full-body compliant model retains spring dynamics and uses \(n=22\) DOF, with constrained equations of motion of the form

\[
D(q)\ddot q + C(q,\dot q)\dot q + G(q)
=
B\,u + J_v(q)^T\lambda + k_s q + k_b\dot q
\]

within each contact domain [1904.11104].

A dominant model-based control line for Cassie combines virtual constraints, hybrid gait libraries, and direct-collocation optimization. One implementation regulates nine universal outputs
\[
y=h(q,\tau,\alpha):=h_0(q)-h_d(\tau,\alpha),
\]
where \(h_d\) is parameterized by spline coefficients and \(\tau\) is a phase variable running from 0 to 1 over a nominal step time \(T=0.4\,\text{s}\) [1809.07279]. Seven periodic gaits were optimized offline for average sagittal speeds from \(-0.5\) to \(1.0\,\text{m/s}\), then interpolated online as a function of filtered forward velocity [1809.07279]. On hardware, this controller enabled standing in place and walking over sidewalks, grass, snow, sand, and burning brush; the standing controller also enabled Segway riding at speeds up to \(\approx 4\,\text{m/s}\) [1809.07279].

Trajectory generation at full order has been accelerated by FROST and C-FROST. For a 20-DOF floating-base Cassie-series model, the hybrid direct-collocation problem is transcribed over right- and left-stance phases with impact resets and periodicity constraints [1807.06614]. C-FROST’s C++ backend and multi-threaded sparse Jacobian evaluation reduce a single walk-in-place optimization from \(90.12\,\text{s}\) in MATLAB/FROST to \(24.51\,\text{s}\) in single-threaded C-FROST, and large batches of gaits can be generated at average times of \(\approx 4.1\,\text{s}\) per gait on a 72-thread AWS EC2 configuration [1807.06614]. That framework was used to compute 1,331 gaits over a three-parameter grid of forward speed, lateral speed, and step height [1807.06614].

A complementary line incorporates passive compliance directly into the dynamics and controller design. “Dynamic Walking with Compliance” uses full-body optimization together with task-space PD and approximate gravity compensation, and reports stable stepping in place over \(10\,\text{s}\) with close agreement between optimized and experimental center-of-mass, joint, and spring-deflection trajectories [1904.11104]. The same controller was reported to work over indoor tile floors, outdoor grass, gravel, and mild slopes [1904.11104].

## 6. Learning, multimodal locomotion, and expressive behavior

Deep reinforcement learning has been used to synthesize Cassie feedback controllers directly from a realistic model. In one PPO-based formulation, the state has dimension 80, the action is a 10-dimensional offset to reference joint angles, and the reward is an exponential imitation objective combining joint tracking, pelvis position, pelvis orientation, and spring tracking [1803.05580]. Relative to a hand-designed reference controller, the learned policy remains stable on sinusoidal terrain up to amplitude \(h=0.15\) rather than \(h=0.07\), and withstands pelvis pushes of \(140\,\text{N}\) forward, \(90\,\text{N}\) backward, and \(50\,\text{N}\) lateral [1803.05580]. Interpolating between speed-specific policies further improves terrain robustness, to \(h=0.22\) in the reported tests [1803.05580].

An iterative-design variant introduces Deterministic Action Stochastic State (DASS) tuples,
\[
\mathcal D=\{(s_i,m_e(s_i))\}_{i=1}^N,
\]
as a mechanism for reward-function redefinition, policy regularization, and distillation across architectures [1903.09537]. The mixed objective combines policy-gradient updates with supervised imitation on DASS data, allowing successive redesign of gait style while constraining deviation from a previously working controller [1903.09537]. Hardware results include stable walking from \(-0.4\) to \(1.2\,\text{m/s}\), transfer from simulation without dynamics randomization, and compact variable-speed policies represented by 5–10k tuples, with a final small-network compression reported at 6k tuples [1903.09537].

Cassie has also served as a platform for multimodal locomotion. In autonomous riding of hovershoes, the robot is modeled jointly with two independent wheeled platforms, and the feedback architecture augments a nominal balancing controller with PD loops for foot alignment, forward velocity, and turning [1907.11353]. The perception stack combines Intel RealSense T265 stereo VIO, D435i depth sensing, RTAB-Map occupancy mapping, Dijkstra global planning, and Timed Elastic Band local planning [1907.11353]. Reported experimental capabilities include step velocity commands \(0.5\to 1.0\to 1.5\,\text{m/s}\) with RMSE \(\approx 0.147\,\text{m/s}\), yaw-rate steps \(0\to 0.5\,\text{rad/s}\to 0\) with RMSE \(\approx 0.145\,\text{rad/s}\), slopes up to \(+7^\circ\) and down to \(-13^\circ\), and stair descent over single steps at \(79\%\) and \(171\%\) of the hovershoe wheel radius [1907.11353].

The platform has also been used for affective robotics. “Animated Cassie” constructs an animation-to-hardware pipeline in which keyframed motions created in Blender are converted to dynamically feasible standing and walking motions by direct-collocation optimization and then executed through a finite-state machine and high-rate controllers [2009.02846]. The robot was used to perform “tired,” “happy,” and “curious” motions, as well as a stand-to-walk-to-stand “laser chase” story, with the optimized standing motion time-stretched by about \(2\times\) relative to the raw animation and with reported \(100\%\) success in mode transitions and less than \(5\%\) velocity error in walking tests [2009.02846]. In the robotics literature, Cassie therefore functions both as a benchmark locomotion platform and as an experimental substrate for dynamic character animation.

Source: https://www.emergentmind.com/topics/cassie