Researcher profile

Steven A. Orszag

Steven A. Orszag contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

cond-mat.mtrl-sci physics.flu-dyn astro-ph.IM chao-dyn nlin.CD physics.class-ph physics.comp-ph

Trust snapshot

Quick read

Trust 17 - Baseline

4works

0followers

7topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Open graph Browse works

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2012arXiv

Nonlinear mushy-layer convection with chimneys: stability and optimal solute fluxes

We model buoyancy-driven convection with chimneys -- channels of zero solid fraction -- in a mushy layer formed during directional solidification of a binary alloy in two-dimensions. A large suite of numerical simulations is combined with scaling analysis in order to study the parametric dependence of the flow. Stability boundaries are calculated for states of finite-amplitude convection with chimneys, which for a narrow domain can be interpreted in terms of a modified Rayleigh number criterion based on the domain width and mushy-layer permeability. For solidification in a wide domain with multiple chimneys, it has previously been hypothesised that the chimney spacing will adjust to optimise the rate of removal of potential energy from the system. For a wide variety of initial liquid concentration conditions, we consider the detailed flow structure in this optimal state and derive scaling laws for how the flow evolves as the strength of convection increases. For moderate mushy-layer Rayleigh numbers these flow properties support a solute flux that increases linearly with Rayleigh number. This behaviour does not persist indefinitely, however, with porosity-dependent flow saturation resulting in sub-linear growth of the solute flux for sufficiently large Rayleigh numbers. Finally, we consider the influence of the porosity dependence of permeability, with a cubic function and a Carmen-Kozeny permeability yielding qualitatively similar system dynamics and flow profiles for the optimal states.

preprint2011arXiv

Lattice Boltzmann Method for Electromagnetic Wave Propagation

We present a new Lattice Boltzmann (LB) formulation to solve the Maxwell equations for electromagnetic (EM) waves propagating in a heterogeneous medium. By using a pseudo-vector discrete Boltzmann distribution, the scheme is shown to reproduce the continuum Maxwell equations. The technique compares well with a pseudo-spectral method at solving for two-dimensional wave propagation in a heterogeneous medium, which by design contains substantial contrasts in the refractive index. The extension to three dimensions follows naturally and, owing to the recognized efficiency of LB schemes for parallel computation in irregular geometries, it gives a powerful method to numerically simulate a wide range of problems involving EM wave propagation in complex media.

preprint2010arXiv

Optimal and hysteretic fluxes in alloy solidification: Variational principles and chimney spacing

We take a numerical approach to analyze the mechanisms controlling the spacing of chimneys -- channels devoid of solid -- in two-dimensional mushy layers formed by solidifying a binary alloy. Chimneys are the principal conduits through which buoyancy effects transport material out of the mushy layer and into the liquid from which it formed. Experiments show a coarsening of chimney spacing and we pursue the hypothesis that this observation is a consequence of a variational principle: the chimney spacing adjusts to optimize material transport and hence maximize the rate of removal of potential energy stored in the mushy layer. The optimal solute flux increases approximately linearly with the mushy layer Rayleigh number. However, for spacings below a critical value the chimneys collapse and solute fluxes cease, revealing a hysteresis between chimney convection and no flow.

preprint1996arXiv

The Effect of Small-Scale Forcing on Large-Scale Structures in Two-Dimensional Flows

The effect of small scale forcing on large scale structures in $β$-plane two-dimensional (2D) turbulence is studied using long-term direct numerical simulations (DNS). We find that nonlinear effects remain strong at all times and for all scales and establish an inverse energy cascade that extends to the largest scales available in the system. The large scale flow develops strong spectral anisotropy: $k^{-5/3}$ Kolmogorov scaling holds for almost all $ϕ$, $ϕ=\arctan (k_y/k_x)$, except in the small vicinity of $k_x=0$, where Rhines's $k^{-5}$ scaling prevails. Due to the $k^{-5}$ scaling, the spectral evolution of $β$-plane turbulence becomes extremely slow which, perhaps, explains why this scaling law has never before been observed in DNS. Simulations with different values of $β$ indicate that the $β$-effect diminishes at small scales where the flow is nearly isotropic. Thus, for simulations of $\b$-plane turbulence forced at small scales sufficiently removed from the scales where $\b$-effect is strong, large eddy simulation (LES) can be used. A subgrid scale (SGS) parameterization for such LES must account for the small scale forcing that is not explicitly resolved and correctly accommodate two inviscid conservation laws, viz. energy and enstrophy. This requirement gives rise to a new anisotropic stabilized negative viscosity (SNV) SGS representation which is discussed in the context of LES of isotropic 2D turbulence.