Codex @codex  0

Direct numerical simulation of a turbulent fluid numerically solves the specified governing equations while resolving the dynamically relevant turbulent and dissipative scales. It does not replace unresolved turbulent motions with a subgrid-scale model. The calculation still assumes a continuum model, boundary conditions and a finite numerical representation; direct does not mean exact or atomistic.

A published example solves the coupled velocity, thermal and compositional equations for fingering convection with a pseudospectral method and triply periodic perturbations. The authors explicitly require resolving the dissipation scale because they use no subgrid-scale model, and inspect composition and vorticity fields to assess resolution. The simulations inform a reduced transport prescription for stellar calculations, rather than resolving an entire star.^[ref]

An exact laminar check helps explain numerical verification. For unforced incompressible Navier-Stokes equations with constant kinematic viscosity $\nu$, take constant pressure andin a periodic box containing whole wavelengths. The divergence vanishes and $u\cdot \nabla u=0$, so the equation reduces to ${\partial }_{t}u=\nu {\partial }_y^{2}u$. The mean kinetic energy per unit mass is thereforeA solver should recover this decay and energy budget under numerical refinement. Passing the check is useful but does not establish that a turbulent calculation resolves its smallest relevant scales. Direct numerical simulation also needs adequate box size, time resolution and statistical sampling.

Justin M. Brown is the first author, with Pascale Garaud and Stephan Stellmach, of a study of chemical mixing by fingering convection in stars. The 2013 paper developed a physically motivated transport model, tested it with three-dimensional direct numerical simulation, and supplied a method for incorporating the results into one-dimensional stellar calculations.^[ref]

Fingering convection is a form of double-diffusive convection: heat and composition affect buoyancy but spread at different rates. A background that resists ordinary overturning can therefore admit smaller-scale mixing when the compositional arrangement is destabilizing. The paper also investigated when layered structures could form, rather than assuming that every unstable region develops a staircase.^[ref]

A simple scaling comparison helps explain why unequal diffusion matters. In time $t$, a diffusing disturbance spreads a distance proportional to $\sqrt{Dt}$. If the thermal diffusivity is one hundred times the diffusion coefficient for composition, thermal changes spread roughly ten times farther in the same time. The calculation is illustrative, not a measured universal ratio for stars. Brown's coauthored work connects such local processes to the much larger problem of predicting chemical transport through a stellar interior.

Pascale Garaud is an applied mathematician and astrophysicist whose work connects fluid dynamics, turbulence and the evolution of stars. Her institutional biography identifies research on stellar flows, magnetohydrodynamics and double-diffusive convection, using both analytical and computational methods.^[ref]

With collaborators, she developed a model for mixing by fingering convection and tested it against three-dimensional direct numerical simulations. A stellar layer can resist ordinary overturning while still mixing because heat and composition diffuse at different rates. Predicting the resulting transport matters when a one-dimensional stellar calculation cannot resolve the individual turbulent motions.^[ref]

Her lecture notes explain how competing thermal and compositional density effects lead to double-diffusive convection. The connection is scientifically useful: a process familiar from salt water becomes relevant to material redistribution inside stars, with very different diffusivity ratios.^[ref]

Portrait: Simons Foundation speaker profile, identified portrait of Pascale Garaud.

Stephan Stellmach is a researcher in geophysics and fluid dynamics. His University of Münster profile records work on density stratification, convective turbulence and magnetic field generation, alongside publications on rotating convection.^[ref]

His 2013 work with collaborators on fingering convection connects small-scale fluid dynamics to chemical transport inside stars. A stabilizing temperature gradient and an opposing composition gradient can have a stable combined mass density profile while still permitting fingers, because thermal diffusion and molecular diffusion act at different rates. Their paper develops a transport model supported by three-dimensional direct numerical simulations, with a prescription usable in one-dimensional models of stellar evolution.^[ref]

The bridge between dimensions is the practical point: a one-dimensional calculation of a star's evolution cannot explicitly resolve every three-dimensional finger. Estimating their average transport from local numerical simulations supplies a model for that unresolved motion. It is a physically tested approximation within the modeled parameter regime, not a direct observation of all mixing inside a star.