D. I. Meiron

Analysis of iterative methods for the steady and unsteady Stokes problem: application to spectral element discretizations

A new and detailed analysis of the basic Uzawa algorithm for decoupling of the pressure and the velocity in the steady and unsteady Stokes operator is presented. The paper focuses on the following new aspects: explicit construction of the Uzawa pressure-operator spectrum for a semiperiodic model problem; general relationship of the convergence rate of the Uzawa procedure to classical inf-sup discretization analysis; and application of the method to high-order variational discretization.

A virtual test facility for the efficient simulation of solid material response under strong shock and detonation wave loading

A virtual test facility (VTF) for studying the three-dimensional dynamic response of solid materials subject to strong shock and detonation waves has been constructed as part of the research program of the Center for Simulating the Dynamic Response of Materials at the California Institute of Technology. The compressible fluid flow is simulated with a Cartesian finite volume method and treating the solid as an embedded moving body, while a Lagrangian finite element scheme is employed to describe the structural response to the hydrodynamic pressure loading. A temporal splitting method is applied to update the position and velocity of the boundary between time steps. The boundary is represented implicitly in the fluid solver with a level set function that is constructed on-the-fly from the unstructured solid surface mesh. Block-structured mesh adaptation with time step refinement in the fluid allows for the efficient consideration of disparate fluid and solid time scales. We detail the design of the employed object-oriented mesh refinement framework AMROC and outline its effective extension for fluid–structure interaction problems. Further, we describe the parallelization of the most important algorithmic components for distributed memory machines and discuss the applied partitioning strategies. As computational examples for typical VTF applications, we present the dynamic deformation of a tantalum cylinder due to the detonation of an interior solid explosive and the impact of an explosion-induced shock wave on a multi-material soft tissue body.

Dynamical aspects of vortex reconnection of perturbed anti-parallel vortex tubes

The phenomenon of vortex reconnection is analysed numerically and the results are compared qualitatively with the predictions of a model of reconnection recently proposed by Saffman. Using spectral methods over both uniform and strained meshes, numerical simulations are performed of two nearly parallel, counter-rotating vortex tubes, over the range of Reynolds numbers Re = 1000–3500. The calculations utilizing a uniform mesh are performed for Re ≤ 1500 with a resolution of 128 points in each direction. The calculations utilizing a stretched mesh are performed for 1500 < Re ≤ 3500 with a resolution of up to 160 points in each direction and with a fourfold stretching about the region of reconnection. We present results for the variation of the maximum of vorticity, the time to reconnection, and other diagnostics of this flow as functions of the Reynolds number. From numerical simulation of the model equations, we infer and demonstrate the existence of exact solutions to the model to which its solutions arising from more general initial conditions are attracted at late times. In the limit of infinite Reynolds number, the model predicts eventual saturation of the axial strain, a feature observed in the recent work of Pumir & Siggia and also observed in our full numerical simulations. In this respect the model captures the observed local dynamics of vortex stretching. However, because the global effects of external flows are not included in the model, the model predicts that the axial strain eventually decays and the maximum vorticity grows linearly at late times. In contrast, from the full simulations, we see the possible emergence of the behaviour of the axial strain at infinite Reynolds number. As our simulations are affected by non-local effects, we do observe saturation of the strain but no subsequent decay. It is also shown analytically that the model predicts a reconnection time which varies logarithmically with increasing Reynolds number. Comparison with the full numerical simulations shows a much slower variation of the reconnection time with increasing Reynolds number than predicted by the model. Other points of agreement and disagreement between the Saffman model and the simulations are discussed, Reconnection is also discussed from the point of view of its relation to the possible onset of nearly singular behaviour of the Euler equation. In agreement with the recent numerical results of Pumir & Siggia, our results suggest that no singularity in the vorticity will form in a finite time for this initial condition.

Numerical simulation of morphological development during ostwald ripening

Abstract A boundary integral technique is employed to determine the morphological evolution of a small number of particles during Ostwald ripening in two dimensions. The approach allows the bodies to change shape consistent with interparticle diffusional interactions and the interfacial concentrations as given by the Gibbs-Thomson equation. It is shown that the strong interparticle diffusional interactions which occur at small interparticle separations can induce significant motions of the centers of mass of the particles. Such motion is shown to be a strong function of the spatial distribution of particles. The generality of the mechanism responsible for the particle migration suggests that particle motion is a generic aspect of the ripening process at high volume fractions of coarsening phase. It was found that significant shape distortions of particles during ripening requires particle arrangements which induce significant diffusional screening of regions of interface. Through particle arrangements similar to those found in solid-liquid systems during liquid phase sintering, it is shown that the formation of regions of flat interface between particles is completely consistent with an Ostwald ripening mechanism.

A Virtual Test Facility for the Simulation of Dynamic Response in Materials

The Center for Simulating Dynamic Response of Materials at the California Institute of Technology is constructing a virtual shock physics facility for studying the response of various target materials to very strong shocks. The Virtual Test Facility (VTF) is an end-to-end, fully three-dimensional simulation of the detonation of high explosives (HE), shock wave propagation, solid material response to pressure loading, and compressible turbulence. The VTF largely consists of a parallel fluid solver and a parallel solid mechanics package that are coupled together by the exchange of boundary data. The Eulerian fluid code and Lagrangian solid mechanics model interact via a novel approach based on level sets. The two main computational packages are integrated through the use of Pyre, a problem solving environment written in the Python scripting language. Pyre allows application developers to interchange various computational models and solver packages without recompiling code, and it provides standardized access to several data visualization engines and data input mechanisms. In this paper, we outline the main components of the VTF, discuss their integration via Pyre, and describe some recent accomplishments in large-scale simulation using the VTF.

The linear two-dimensional stability of inviscid vortex streets of finite-cored vortices

The stability of two-dimensional infinitesimal disturbances of the inviscid Karman vortex street of finite-area vortices is reexamined. Numerical results are obtained for the growth rate and oscillation frequencies of disturbances of arbitrary subharmonic wavenumber and the stability boundaries are calculated. The stabilization of the pairing instability by finite area demonstrated by Saffman & Schatzman (1982) is confirmed, and also Kida’s (1982) result that this is not the most unstable disturbance when the area is finite. But, contrary to Kida’s quantitative predictions, it is now found that finite area does not stabilize the street to infinitesimal two-dimensional disturbances of arbitrary wavelength and that it is always unstable except for one isolated value of the aspect ratio which depends upon the size of the vortices. This result does agree, however, with those of a modified version of Kida’s analysis.

Overhanging interfacial gravity waves of large amplitude

Methods to investigate the existence of overhanging gravity waves of permanent form at the interface between two uniform fluids of different density are discussed. Numerical results which demonstrate their existence are presented.

Atwood ratio dependence of Richtmyer-Meshkov flows under reshock conditions using large-eddy simulations

We study the shock-driven turbulent mixing that occurs when a perturbed planar density interface is impacted by a planar shock wave of moderate strength and subsequently reshocked. The present work is a systematic study of the influence of the relative molecular weights of the gases in the form of the initial Atwood ratio A. We investigate the cases A = ± 0.21, ±0.67 and ±0.87 that correspond to the realistic gas combinations air–CO_2, air–SF_6 and H_2–air. A canonical, three-dimensional numerical experiment, using the large-eddy simulation technique with an explicit subgrid model, reproduces the interaction within a shock tube with an endwall where the incident shock Mach number is ~1.5 and the initial interface perturbation has a fixed dominant wavelength and a fixed amplitude-to-wavelength ratio ~0.1. For positive Atwood configurations, the reshock is followed by secondary waves in the form of alternate expansion and compression waves travelling between the endwall and the mixing zone. These reverberations are shown to intensify turbulent kinetic energy and dissipation across the mixing zone. In contrast, negative Atwood number configurations produce multiple secondary reshocks following the primary reshock, and their effect on the mixing region is less pronounced. As the magnitude of A is increased, the mixing zone tends to evolve less symmetrically. The mixing zone growth rate following the primary reshock approaches a linear evolution prior to the secondary wave interactions. When considering the full range of examined Atwood numbers, measurements of this growth rate do not agree well with predictions of existing analytic reshock models such as the model by Mikaelian (Physica D, vol. 36, 1989, p. 343). Accordingly, we propose an empirical formula and also a semi-analytical, impulsive model based on a diffuse-interface approach to describe the A-dependence of the post-reshock growth rate.

Calculation of steady three-dimensional deep-water waves

Steady three-dimensional symmetric wave patterns for finite-amplitude gravity waves on deep water are calculated from the full unapproximated water-wave equations as well as from an approximate equation due to Zakharov. These solutions are obtained as bifurcations from plane Stokes waves. The results are in good agreement with the experimental observations of Su.

On the stability of gas bubbles rising in an inviscid fluid

The stability to three dimensional disturbances of bubbles rising rectilinearly in an inviscid fluid is studied numerically. It is found, in contrast with earlier work, that the interaction of hydrodynamic pressure forces and surface tension does not lead to linear instability of the bubble path.