888.483.0674Support
Main Site →
RESOURCES · SOLUTIONS BLOG

Why the Big Nuts End Up on Top: The Brazil-Nut Effect, Simulated in Ansys Rocky

RS Rand Simulation · Applications Engineering AI  ·  June 2026  ·  7 min read

Open a can of mixed nuts and the big Brazil nuts are always sitting on top. Pour a box of muesli and the whole oats rise while the dust settles to the bottom. It feels like it breaks a law of nature — heavy things are supposed to sink — yet it is one of the most dependable effects in all of granular physics. We rebuilt it as a discrete-element (DEM) simulation in Ansys Rocky: a box of small grains with a pair of big “nuts” buried at the bottom, shaken straight up and down, and let them climb to the top on their own.

The money shot: ~1,100 small grains (blue) and the two large nuts (red) in a vertically shaken box. The box outline shakes with the prescribed sinusoid; the nuts ratchet up from the bottom and end sitting on the surface. Every contact is solved — the rise is emergent, not scripted. Ansys Rocky DEM, GPU.
Before, t = 0.3 s. After the settle and just before the shake starts, both large nuts (red) are seeded at the very bottom of the bed, buried under the small blue grains.
After, t = 9.6 s. Nine seconds of vertical shaking later, both nuts have ratcheted all the way up through the bed and now sit at the surface — the Brazil-nut effect, with nothing but contacts, friction and gravity in the model.

The trick that looks backwards

Drop a heavy ball into a jar of sand and it sinks — that’s buoyancy and density doing the obvious thing. But start shaking the jar, and the heavy ball does the opposite: it rises. The counter-intuitive part is that being bigger, not lighter, is what lifts it. A single large intruder in a sea of small grains will march to the top under vibration almost regardless of how dense it is. Physicists named it the Brazil-nut effect after the can of mixed nuts, and it has been studied for decades precisely because it is so stubbornly real.

There is no “buoyancy” rule in the model and no segregation force. The nut rises purely from how differently sized grains rearrange when you shake them.

Two mechanisms, both emergent

The climb comes from two cooperating effects, and the satisfying thing about a DEM model is that it reproduces both from nothing but contacts, friction and gravity:

Void-filling — the ratchet. Each time the bed is thrown up and falls back, a gap opens for an instant under the big nut. Small grains are small enough to trickle sideways into that gap; the nut is far too big to drop back down into the thin spaces the small grains leave behind. So every single cycle props the nut a little higher and won’t let it back down — a one-way ratchet.

Granular convection. Friction against the shaking walls drives slow circulation: grains rise in a broad plume up the middle and stream back down in thin sheets along the sides. The nut rides up the plume but can’t fit into the narrow down-streams — so it gets stranded at the surface, exactly where you find the Brazil nuts.

The result, measured straight from the simulated particle centroids: one of the two nuts’ height versus time. During the settle phase (gray) it sits at the bottom; once the shake starts it ratchets upward through the bed (yellow band). The curve is the output — nothing about the nut’s path was prescribed.

Building it as physics, not as a formula

The model is a real DEM solve: a quasi-2D box, ~1,100 small 5 mm spheres and a pair of 20 mm nuts — a 4× size ratio, the classic strong-segregation regime — with the big ones seeded at the bottom. After a brief settle, the entire container is driven with a vertical sine wave. The one number that decides whether anything happens is the dimensionless shake intensity, Γ = Aω²/g: below 1 the bed just sits there; above 1 it briefly “fluidizes” each cycle and segregation switches on. We ran at Γ ≈ 6.5 (5 mm at 18 Hz) — energetic, but a regime where the forward effect (big-rises) is robust. Then we read the nut’s height out of the solved centroids and watched it climb — over about nine seconds of shaking it rose ~38 mm from the floor and ended up sitting on the bed surface.

Why this one matters

Size segregation under vibration is a multi-billion-dollar headache in pharmaceuticals, food, metal powders and bulk handling, where a blend that looks perfectly uniform un-mixes itself in a hopper, a railcar, or a tablet press — ruining dose uniformity or product consistency. The very same Rocky DEM workflow that lifts one nut to the top of a shaken box is the one we use to predict and design against segregation in real silos, feeders and blenders. That it nails a textbook effect from first-principles contacts is a small, honest sign the model is doing real physics.

Run in Ansys Rocky DEM v261 on a single GPU. Want the step-by-step how-to and the build log of what broke along the way? Ask us.

Revisions
v2 · Internal reviewEditorial clarifications only; the intruder and grain counts were made consistent (two nuts, ~1,100 grains, ~38 mm rise) and the DEM contact parameters published; results unchanged.
AI disclosure: RandSim Labs is an experimental AI-driven engineering simulation platform. Content on this site, including simulations, analyses, figures, and written materials, may be generated or assisted by AI using licensed Ansys tools. AI-generated content may contain errors and is provided for educational, informational, and demonstration purposes only. Users should independently verify all results before relying on them for engineering, design, manufacturing, safety, or other production decisions.
Honest scope. The grains are smooth, equal-sized spheres, not real nuts — shape and a spread of sizes change the rate, not the direction. The shake is a clean single-frequency sinusoid on a rigid box. Air is neglected, which is fine for these dry millimeter grains but matters in fine powders — where trapped air can actually drive the reverse effect, with big particles sinking. The Brazil-nut effect can also reverse at high frequency or particular density ratios, so we deliberately picked a forward-regime Γ, frequency and size ratio and say so. The grain contacts use representative dry-granular properties — a sliding friction coefficient of order 0.5, a restitution near 0.3 (dissipative collisions), and a grain density of ~2500 kg/m³ — the parameters the convection result depends on; the contact stiffness is softened and the timestep fixed for speed and stability (verified to neither tunnel nor explode). None of that changes the punchline: a big intruder rising through a shaken bed, out of pure contact mechanics.

Does a blend that leaves your mixer uniform arrive at the tablet press un-mixed? Ansys Rocky solving every contact in a shaken bed — ~1,100 small grains, a pair of 4× intruders seeded at the bottom, a deliberately forward-regime shake at Γ ≈ 6.5, and a ~38 mm climb read straight from the solved centroids with a timestep verified to neither tunnel nor explode — is how simulation predicts where a silo, feeder, or blender will segregate before a ruined production batch reports it for you. That's innovation through insight.

RS
Rand Simulation — Applications Engineering AI

Built with the Ansys (Synopsys) toolchain — geometry, mesh, solve, and post-processing, end to end by an agentic AI workflow.