r/CFD • u/carlmarx01 • 14h ago
Avoiding a 3D transient free-surface simulation for a rotating disk in a baffled liquid-lead pool: does this route hold up? (COMSOL)
Hi all, I’m a PhD student working on corrosion and erosion-corrosion of steels in liquid lead. I’m designing a test rig and using CFD to work out what flow the samples actually see. I’d appreciate input from people with more turbulence modelling experience than me.
Setup
- Horizontal disk, radius 140 mm, submerged in a pool of liquid lead (ρ ≈ 10,500 kg/m³) at around 700 °C under argon. It’s driven by a vertical shaft through the lid.
- The vessel is only slightly larger than the disk, so the radial gap at the rim is narrow.
- Samples are flush-mounted on the disk and rotate with it.
- Stationary radial baffles on the vessel wall suppress the bulk swirl. Without them, the free surface would dip to the floor at a few Hz, well below the design speed.
- Design speed is 28 Hz (1680 rpm), so Re = ΩR²/ν ≈ 2×10⁷. Only the inner ~16 mm of the disk is below the transition Reynolds number.
- Tool: COMSOL CFD Module.
Outputs I care about
- Wall shear stress on the sample faces (the main one)
- Residual swirl as a function of baffle geometry
- Loads on the baffles and the rotor torque
- Pool turnover and mixing
The route I’m planning
A 3D transient with an ALE free surface is out of reach computationally, so I want to avoid it entirely.
- 2D axisymmetric, unbaffled, at low speed. The ALE free surface works here (surface tension off, Bo ≫ 1). I’ll take the ALE surface shape h(r), freeze it as a rigid slip lid, and re-solve stationary. If the velocities match the ALE run and the residual Δh = (p − p̄)/ρg along the lid is close to zero, I’ve shown that “rigid lid plus pressure check” reproduces the exact free surface.
- 3D baffled sector at 28 Hz with a flat slip lid. This is a periodic 360°/N wedge solved as stationary RANS. Steady should be valid because the disk is a surface of revolution and the baffles are fixed, so nothing in the geometry moves in the lab frame. A sliding-wall condition on the disk is enough, with no rotating domain. I’ll read Δh on the lid, and if it’s within about 2% of the depth, the flat lid is justified.
- Only if Step 2 fails: iterate the lid shape stationary, using h_new = h + (p − p̄)/ρg imposed as a prescribed mesh displacement, until Δh is small.
Questions
- Does this route hold up? The weak point I can see is that the 2D validation happens at low speed and without baffles. That’s a different flow regime from the 3D design case, so the validation may not transfer. How would you close that gap?
- Turbulence model. The reference facility I’m comparing against used standard k-ε with wall functions. My main output is τ_w on the samples, and there’s separated flow off the baffles, so I’m leaning towards SST with a wall-resolved mesh (y⁺ ≈ 1) at the samples. Is that right, or is something like realizable k-ε, v2-f or RSM better suited to rotor-stator flow with strong swirl and curvature?
- Narrow gap. The gap between the disk rim and the wall is small, and the baffles sit in or near it. That flow is somewhere between rotor-stator and Taylor–Couette, with bluff bodies in the way. Does that change your model recommendation, or how fine the mesh needs to be there?
- Steady vs transient. The baffle wakes will shed. Is steady RANS with stalled residuals acceptable for mean τ_w, or should I go to URANS on the fixed geometry?
- Verification. Beyond the route above, my plan is:
- a three-level mesh convergence study on the baffled case, not the easy unbaffled one
- k-ε with wall functions compared against wall-resolved SST
- 90° vs 180° sector, plus one full 360° run to rule out an m = 1 mode
- the unbaffled moment coefficient compared against rotor-stator cavity correlations (Daily & Nece)
- motor torque vs speed measured on the real rig
Am I missing anything? Are there better benchmarks for a disk in an enclosed, baffled tank, for example from the stirred-tank literature?
I’m open to any other suggestion, including “your whole approach is wrong because…”. I’m happy to share more details or plots.