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.

A fluid dynamicist studies how fluids move and how their motion exchanges mass, momentum and energy. The work can combine mathematical models, controlled measurements and computation. Choosing what to resolve, what to approximate and which observations test the result is as central as solving an equation. A fluid dynamicist may work on applications ranging from aerodynamic drag force to transport in a blood.

For an original example, consider steady water flow through a duct with negligible density changes. Mass conservation gives $\rho A_1{U}_1=\rho A_2{U}_2$, with $U$ the cross-sectionally averaged speed.^[ref] If a circular duct doubles its radius while maintaining the same volume flow rate, its area becomes four times larger and its mean speed falls to one quarter. A researcher can compare that prediction with measured flow and velocity. Predicting the required pressure difference is another question, requiring information about losses and the velocity field.

In turbulence research, direct numerical simulation and large eddy simulation resolve different ranges of motion. A useful validation compares the same observable and averaging convention rather than simply comparing two pictures of complicated flow. Parviz Moin's institutional research profile describes this combination of computational methods with questions about turbulent-flow physics.^[ref]

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.

A large eddy simulation computes the time-dependent large-scale motion of a turbulent flow while representing the influence of unresolved motion through a subgrid-scale model. A spatial averaging filter specifies which information is retained. The distinction concerns spatial resolution, rather than simply replacing a fluctuating flow by its time average.^[ref]

For an elementary illustration, take one component of a periodic velocity field,An ideal filter retaining Fourier modes with integer index at most four keeps the first term and removes the second, whose index is ten. The removed term still contributes to products such as $u_x^2$. Filtering the nonlinear Navier-Stokes equations consequently introduces a subgrid-scale stress that cannot generally be reconstructed from the retained velocity field alone.

This gives a computational compromise. Direct numerical simulation must resolve the dynamically relevant dissipative motion, whereas a large eddy simulation resolves a chosen coarser field and models its missing influence. Doubling grid counts in all three directions multiplies the cell count by eight, before any additional time-step cost. The filter, numerical discretization and closure must work together; a coarse grid by itself does not certify a useful approximation.

An early concrete demonstration was James W. Deardorff's 1970 calculation of turbulent channel flow using 6720 uniform grid intervals and strain-dependent coefficients for subgrid effects. The paper compared computed statistics with experimental measurements and reported agreement ranging from good to marginal. That comparison tests a particular simulation, rather than establishing that every resolved structure or model is exact.^[ref]

Implementations can supply subgrid dissipation in different ways. The Cardinal documentation, for example, describes a spectral filter used to drain energy from the smallest retained wavelengths; that is a particular numerical treatment, not the definition of all large eddy simulation.^[ref]

Parviz Moin is a researcher in computational fluid dynamics and turbulence. Stanford's institutional profile identifies him as the founding director of the Center for Turbulence Research and describes his use of direct numerical simulation and large eddy simulation to investigate turbulent shear flows.^[ref]

With Massimo Germano, Ugo Piomelli and William H. Cabot, he coauthored work on the dynamic subgrid-scale model. The approach uses information at different resolved scales to determine a model coefficient during the calculation.^[ref] The participant account places the collaboration at the 1990 summer program and distinguishes that report from the subsequent 1991 journal paper.^[ref]

A useful scientific distinction is between a numerical method's resolution and its physical accuracy. A finely detailed flow visualization can still give an inaccurate mean drag force. Conversely, a carefully tested coarse model can predict a specified average while omitting individual small eddies. Moin's computational research makes the choice of observable, reference calculation and physical question part of evaluating a simulation.

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.

Raffaele Marino is a physicist whose research includes turbulence in space plasmas and fluids with density stratification. A first-authored study investigated intense vertical drafts using direct numerical simulation of equations with the Boussinesq approximation. In some simulated regimes, roughly $10\%$ of the volume accounted for as much as $50\%$ of the total dissipation.^[ref]

These figures concern the modeled flows, rather than a universal measured fraction of the ocean. They illustrate why the average can hide strong spatial concentration. Dividing the two fractions gives $0.50/0.10=5$: the active region's mean dissipation density is five times the whole-volume mean. The remaining $90\%$ supplies the other half, so its mean is only $0.50/0.90$ times the whole-volume mean. The contrast between the two regions is consequently ninefold.

The result gives an illuminating connection between intermittency and physical mixing: a relatively small region can dominate energy conversion. The simulations test a proposed mechanism for variability; identifying the same mechanism in a particular environment also requires observations of that environment.

A shell model of turbulence replaces a spatially resolved fluid by a sequence of amplitudes associated with progressively smaller length scales. A typical choice assigns a complex velocity amplitude $u_n$ to the wave-number shell $k_n=k_0{\lambda }^n$, with $\lambda >1$. Neighboring shells exchange energy through nonlinear couplings, while forcing supplies energy and viscosity damps high wave numbers. The model retains a scale cascade at far lower cost than direct numerical simulation, while discarding most spatial geometry.^[ref]

The basic energy bookkeeping is transparent. WriteIf the nonlinear coupling and end conditions satisfy ${\sum }_n\mathrm{Re}(u_n^{\ast }{N}_n)=0$, differentiating givesThe nonlinearity transfers energy between shells without creating its total. Viscous damping becomes faster at larger $k_n$: an isolated unforced linear mode decays as $e^{-\nu k_n^{2}t}$.

For an illustrative constant-flux estimate, a transfer time ${\tau }_n\sim (k_n|u_n|{)}^{-1}$ gives $\epsilon \sim |u_n{|}^2/{\tau }_n\sim k_n|u_n{|}^3$, hence $|u_n|\sim (\epsilon /k_n{)}^{1/3}$. Deviations of higher moments from a single scaling rule can then diagnose intermittency. The estimate is dimensional reasoning, not an exact solution. The published improved coupling was designed to reduce spurious shell correlations while preserving useful conservation and scaling properties.^[ref]

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.

William H. Cabot coauthored the dynamic subgrid-scale model developed with Massimo Germano, Ugo Piomelli and Parviz Moin. Their 1990 report compared modeled unresolved stresses with data from direct numerical simulation and tested large eddy simulation of transitional and turbulent channel flow.^[ref] A participant retrospective identifies Cabot as a Stanford research associate during the collaboration.^[ref]

The numerical comparison illustrates an important issue in turbulence research: matching the objects being compared. A coarse simulation describes resolved motion, while a finely resolved reference contains additional fluctuations. The report filtered the reference data before some comparisons, so the comparison did not automatically count intentionally omitted small scales as errors in a resolved field.^[ref]

For an original arithmetic illustration, fine-grid velocities $0$ and $2$ have mean $1$ and mean square $2$. The square of their mean is only $1$. Replacing both samples by the mean preserves that mean but loses a fluctuation contribution of $1$. This elementary difference is why an averaged velocity and an averaged momentum flux contain different information, and why subgrid-scale stress requires explicit treatment.