Eulerian method for multiphase interactions of soft solid bodies in fluids

We introduce an Eulerian approach for problems involving one or more soft solids immersed in a fluid, which permits mechanical interactions between all phases. The reference map variable is exploited to simulate finite-deformation constitutive relations in the solid(s) on the same fixed grid as the fluid phase, which greatly simplifies the coupling between phases. Our coupling procedure, a key contribution in the current work, is shown to be computationally faster and more stable than an earlier approach, and admits the ability to simulate both fluid--solid and solid--solid interaction between submerged bodies. The interface treatment is demonstrated with multiple examples involving a weakly compressible Navier--Stokes fluid interacting with a neo-Hookean solid, and we verify the method's convergence. The solid contact method, which exploits distance-measures already existing on the grid, is demonstrated with two examples. A new, general routine for cross-interface extrapolation is introduced and used as part of the new interfacial treatment

Efficient computation of two-dimensional steady free-surface flows

We consider a family of steady free-surface flow problems in two dimensions, concentrating on the effect of nonlinearity on the train of gravity waves that appear downstream of a disturbance. By exploiting standard complex variable techniques, these problems are formulated in terms of a coupled system of Bernoulli's equation and an integral equation. When applying a numerical collocation scheme, the Jacobian for the system is dense, as the integral equation forces each of the algebraic equations to depend on each of the unknowns. We present here a strategy for overcoming this challenge, which leads to a numerical scheme that is much more efficient than what is normally employed for these types of problems, allowing for many more grid points over the free surface. In particular, we provide a simple recipe for constructing a sparse approximation to the Jacobian that is used as a preconditioner in a Jacobian-free Newton-Krylov method for solving the nonlinear system. We use this approach to compute numerical results for a variety of prototype problems including flows past pressure distributions, a surface-piercing object and bottom topographies.Comment: 20 pages, 13 figures, under revie

A note on the Stokes phenomenon in flow under an elastic sheet: Stokes Phenomenon in flow under a sheet

The Stokes phenomenon is a class of asymptotic behaviour that was first discovered by Stokes in his study of the Airy function. It has since been shown that the Stokes phenomenon plays a significant role in the behaviour of surface waves on flows past submerged obstacles. A detailed review of recent research in this area is presented, which outlines the role that the Stokes phenomenon plays in a wide range of free surface flow geometries. The problem of inviscid, irrotational, incompressible flow past a submerged step under a thin elastic sheet is then considered. It is shown that the method for computing this wave behaviour is extremely similar to previous work on computing the behaviour of capillary waves. Exponential asymptotics are used to show that free-surface waves appear on the surface of the flow, caused by singular fluid behaviour in the neighbourhood of the base and top of the step. The amplitude of these waves is computed and compared to numerical simulations, showing excellent agreements between the asymptotic theory and computational solutions. This article is part of the theme issue 'Stokes at 200 (part 2)'

Eddy genesis and manipulation in plane laminar shear flow

Eddy formation and presence in a plane laminar shear flow configuration consisting of two infinitely long plates orientated parallel to each other is investigated theoretically. The upper plate, which is planar, drives the flow; the lower one has a sinusoidal profile and is fixed. The governing equations are solved via a full finite element formulation for the general case and semi-analytically at the Stokes flow limit. The effects of varying geometry (involving changes in the mean plate separation or the amplitude and wavelength of the lower plate) and inertia are explored separately. For Stokes flow and varying geometry, excellent agreement between the two methods of solution is found. Of particular interest with regard to the flow structure is the importance of the clearance that exists between the upper plate and the tops of the corrugations forming the lower one. When the clearance is large, an eddy is only present at sufficiently large amplitudes or small wavelengths. However, as the plate clearance is reduced, a critical value is found which triggers the formation of an eddy in an otherwise fully attached flow for any finite amplitude and arbitrarily large wavelength. This is a precursor to the primary eddy to be expected in the lid-driven cavity flow which is formed in the limit of zero clearance between the plates. The influence of the flow driving mechanism is assessed by comparison with corresponding solutions for the case of gravity-driven fluid films flowing over an undulating substrate. When inertia is present, the flow generally becomes asymmetrical. However, it is found that for large mean plate separations the flow local to the lower plate becomes effectively decoupled from the inertia dominated overlying flow if the wavelength of the lower plate is sufficiently small. In such cases the local flow retains its symmetry. A local Reynolds number based on the wavelength is shown to be useful in characterising these large-gap flows. As the mean plate separation is reduced, the form of the asymmetry caused by inertia changes, and becomes strongly dependent on the plate separation. For lower plate wavelengths which do not exhibit a cinematically induced secondary eddy, an inertially induced secondary eddy can be created if the mean plate separation is sufficiently small and the global Reynolds number sufficiently large

An Eulerian projection method for quasi-static elastoplasticity

A well-established numerical approach to solve the Navier--Stokes equations for incompressible fluids is Chorin's projection method, whereby the fluid velocity is explicitly updated, and then an elliptic problem for the pressure is solved, which is used to orthogonally project the velocity field to maintain the incompressibility constraint. In this paper, we develop a mathematical correspondence between Newtonian fluids in the incompressible limit and hypo-elastoplastic solids in the slow, quasi-static limit. Using this correspondence, we formulate a new fixed-grid, Eulerian numerical method for simulating quasi-static hypo-elastoplastic solids, whereby the stress is explicitly updated, and then an elliptic problem for the velocity is solved, which is used to orthogonally project the stress to maintain the quasi-staticity constraint. We develop a finite-difference implementation of the method and apply it to an elasto-viscoplastic model of a bulk metallic glass based on the shear transformation zone theory. We show that in a two-dimensional plane strain simple shear simulation, the method is in quantitative agreement with an explicit method. Like the fluid projection method, it is efficient and numerically robust, making it practical for a wide variety of applications. We also demonstrate that the method can be extended to simulate objects with evolving boundaries. We highlight a number of correspondences between incompressible fluid mechanics and quasi-static elastoplasticity, creating possibilities for translating other numerical methods between the two classes of physical problems.Comment: 49 pages, 20 figure