An implicit immersed boundary method for three-dimensional fluid–membrane interactions

We present an implicit immersed boundary method for the incompressible Navier–Stokes equations capable of handling three-dimensional membrane–fluid flow interactions. The goal of our approach is to greatly improve the time step by using the Jacobian-free Newton–Krylov method (JFNK) to advance the location of the elastic membrane implicitly. The most attractive feature of this Jacobian-free approach is Newton-like nonlinear convergence without the cost of forming and storing the true Jacobian. The Generalized Minimal Residual method (GMRES), which is a widely used Krylov-subspace iterative method, is used to update the search direction required for each Newton iteration. Each GMRES iteration only requires the action of the Jacobian in the form of matrix–vector products and therefore avoids the need of forming and storing the Jacobian matrix explicitly. Once the location of the boundary is obtained, the elastic forces acting at the discrete nodes of the membrane are computed using a finite element model. We then use the immersed boundary method to calculate the hydrodynamic effects and fluid–structure interaction effects such as membrane deformation. The present scheme has been validated by several examples including an oscillatory membrane initially placed in a still fluid, capsule membranes in shear flows and large deformation of red blood cells subjected to stretching force.     

Numerical simulations of two-dimensional foam by the immersed boundary method

In this paper, we present an immersed boundary (IB) method to simulate a dry foam, i.e., a foam in which most of the volume is attributed to its gas phase. Dry foam dynamics involves the interaction between a gas and a collection of thin liquid-film internal boundaries that partition the gas into discrete cells or bubbles. The liquid-film boundaries are flexible, contract under the influence of surface tension, and are permeable to the gas, which moves across them by diffusion at a rate proportional to the local pressure difference across the boundary. Such problems are conventionally studied by assuming that the pressure is uniform within each bubble. Here, we introduce instead an IB method that takes into account the non-equilibrium fluid mechanics of the gas. To model gas diffusion across the internal liquid-film boundaries, we allow normal slip between the boundary and the gas at a velocity proportional to the (normal) force generated by the boundary surface tension. We implement this method in the two-dimensional case, and test it by verifying the von Neumann relation, which governs the coarsening of a two-dimensional dry foam. The method is further validated by a convergence study, which confirms its first-order accuracy.     

An accurate and efficient method for the incompressible Navier–Stokes equations using the projection method as a preconditioner

The projection method is a widely used fractional-step algorithm for solving the incompressible Navier–Stokes equations. Despite numerous improvements to the methodology, however, imposing physical boundary conditions with projection-based fluid solvers remains difficult, and obtaining high-order accuracy may not be possible for some choices of boundary conditions. In this work, we present an unsplit, linearly-implicit discretization of the incompressible Navier–Stokes equations on a staggered grid along with an efficient solution method for the resulting system of linear equations. Since our scheme is not a fractional-step algorithm, it is straightforward to specify general physical boundary conditions accurately; however, this capability comes at the price of having to solve the time-dependent incompressible Stokes equations at each timestep. To solve this linear system efficiently, we employ a Krylov subspace method preconditioned by the projection method. In our implementation, the subdomain solvers required by the projection preconditioner employ the conjugate gradient method with geometric multigrid preconditioning. The accuracy of the scheme is demonstrated for several problems, including forced and unforced analytic test cases and lid-driven cavity flows. These tests consider a variety of physical boundary conditions with Reynolds numbers ranging from 1 to 30000. The effectiveness of the projection preconditioner is compared to an alternative preconditioning strategy based on an approximation to the Schur complement for the time-dependent incompressible Stokes operator. The projection method is found to be a more efficient preconditioner in most cases considered in the present work.     

A high-order cell-centered Lagrangian scheme for compressible fluid flows in two-dimensional cylindrical geometry

The goal of this paper is to present high-order cell-centered schemes for solving the equations of Lagrangian gas dynamics written in cylindrical geometry. A node-based discretization of the numerical fluxes is obtained through the computation of the time rate of change of the cell volume. It allows to derive finite volume numerical schemes that are compatible with the geometric conservation law (GCL). Two discretizations of the momentum equations are proposed depending on the form of the discrete gradient operator. The first one corresponds to the control volume scheme while the second one corresponds to the so-called area-weighted scheme. Both formulations share the same discretization for the total energy equation. In both schemes, fluxes are computed using the same nodal solver which can be viewed as a two-dimensional extension of an approximate Riemann solver. The control volume scheme is conservative for momentum, total energy and satisfies a local entropy inequality in its first-order semi-discrete form. However, it does not preserve spherical symmetry. On the other hand, the area-weighted scheme is conservative for total energy and preserves spherical symmetry for one-dimensional spherical flow on equi-angular polar grid. The two-dimensional high-order extensions of these two schemes are constructed employing the generalized Riemann problem (GRP) in the acoustic approximation. Many numerical tests are presented in order to assess these new schemes. The results obtained for various representative configurations of one and two-dimensional compressible fluid flows show the robustness and the accuracy of our new schemes.     

An accurate adaptive solver for surface-tension-driven interfacial flows

A method combining an adaptive quad/octree spatial discretisation, geometrical Volume-Of-Fluid interface representation, balanced-force continuum-surface-force surface-tension formulation and height-function curvature estimation is presented. The extension of these methods to the quad/octree discretisation allows adaptive variable resolution along the interface and is described in detail. The method is shown to recover exact equilibrium (to machine accuracy) between surface-tension and pressure gradient in the case of a stationary droplet, irrespective of viscosity and spatial resolution. Accurate solutions are obtained for the classical test case of capillary wave oscillations. An application to the capillary breakup of a jet of water in air further illustrates the accuracy and efficiency of the method. The source code of the implementation is freely available as part of the Gerris flow solver.     

Conservative numerical simulation of multi-component transport in two-dimensional unsteady shallow water flow

An explicit finite volume model to simulate two-dimensional shallow water flow with multi-component transport is presented. The governing system of coupled conservation laws demands numerical techniques to avoid unrealistic values of the transported scalars that cannot be avoided by decreasing the size of the time step. The presence of non conservative products such as bed slope and friction terms, and other source terms like diffusion and reaction, can make necessary the reduction of the time step given by the Courant number. A suitable flux difference redistribution that prevents instability and ensures conservation at all times is used to deal with the non-conservative terms and becomes necessary in cases of transient boundaries over dry bed. The resulting method belongs to the category of well-balanced Roe schemes and is able to handle steady cases with flow in motion. Test cases with exact solution, including transient boundaries, bed slope, friction, and reaction terms are used to validate the numerical scheme. Laboratory experiments are used to validate the techniques when dealing with complex systems as the κ–?$\kappa ''–''?$ model. The results of the proposed numerical schemes are compared with the ones obtained when using uncoupled formulations.     

A conservative discrete compatibility-constraint low-Mach pressure-correction algorithm for time-accurate simulations of variable density flows

In this paper, we develop the discrete compatibility-constraint pressure-correction algorithm for transient simulations of variable density flows at low-Mach numbers. The constraint for the velocity field is constructed from a combination of the discrete equations of continuity and scalar (e.g. energy) transport, imposing that the newly predicted state must be compatible, in agreement with the equation of state. This way, mass and scalar conservation are guaranteed and the equation of state is exactly fulfilled at every time step. For comparison reasons, two other types of well-known pressure-correction algorithms are also used. The first class, denoted as continuity-constraint pressure-correction, is based on a constraint for the velocity field that is derived solely from the continuity equation. The second class, denoted as analytical compatibility-constraint pressure-correction, constructs the constraint from an analytical combination of the material derivative of the equation of state and the continuity and scalar equations. The algorithms are tested for three example fluid configurations: a single-fluid ideal gas, a two-fluid inert mixture and a two-fluid reacting mixture. The latter is special in the sense that the equation of state is non-linear and not everywhere differentiable. The continuity-constraint pressure-correction algorithm yields unstable solutions if density ratios are high. The analytical compatibility-constraint pressure-correction algorithm yields stable results, but the predicted states do not correspond to the equation of state. The discrete compatibility-constraint pressure-correction algorithm performs well on all test cases: the simulation results are stable and exactly match the equation of state.     

Uniqueness of balanced metrics on holomorphic vector bundles

Let E→M$E\to M$ be a holomorphic vector bundle over a compact Kähler manifold (M,ω)$(M,\omega )$. We prove that if E$E$ admits a ω$\omega$-balanced metric (in X. Wang’s terminology (Wang, 2005 [3])) then it is unique. This result together with Biliotti and Ghigi (2008) [14] implies the existence and uniqueness of ω$\omega$-balanced metrics of certain direct sums of irreducible homogeneous vector bundles over rational homogeneous varieties. We finally apply our result to show the rigidity of ω$\omega$-balanced Kähler maps into Grassmannians.     

Constructive inverse function theorems

An explicit method for finding every coefficient of the reversed series of a power series in one variable is presented. We also show how implicit functions on the plane may be solved for one of the variables. Our approach to the Lagrange inversion formula is based on factorization properties of partitions of integers and distributions of distinguishable objects. These results and techniques may have many applications in numerical analysis, combinatorics, and soliton theory.     

Bounds for the chaotic region in the Lorenz model

In a previous paper, the authors made an extensive numerical study of the Lorenz model, changing all three parameters of the system. We conjectured that the region of parameters where the Lorenz model is chaotic is bounded for fixed r. In this paper, we give a theoretical proof of the conjecture by obtaining theoretical bounds for the chaotic region and by using Fenichel theory. The theoretical bounds are complemented with numerical studies performed using the Maximum Lyapunov Exponent and OFLI2 techniques, and a comparison of both sets of results is shown. Finally, we provide a complete three-dimensional model of the chaotic regime depending on the three parameters.     

Separation of acoustic waves in isentropic flow perturbations

The present contribution investigates the mechanisms of sound generation and propagation in the case of highly-unsteady flows. Based on the linearisation of the isentropic Navier-Stokes equation around a new pathline-averaged base flow, it is demonstrated for the first time that flow perturbations of a non-uniform flow can be split into acoustic and vorticity modes, with the acoustic modes being independent of the vorticity modes. Therefore, we can propose this acoustic perturbation as a general definition of sound.     

The cohomology of the classical and quantum Weyl CCR in curved spaces

An analysis is presented of the cohomological underpinnings for the Weyl group of the canonical commutation relations on manifolds of constant negative curvature. Several uniqueness results are obtained leading from purely classical considerations to the group associated with the systems of imprimitivity of the orthodox approach to quantum mechanics.     

C *-algebraic quantization and the origin of topological quantum effects

Concepts from the theory of abstract operator algebras are used to solve the problem of quantizing a particle moving on an arbitrary locally compact homogeneous space. Inequivalent quantizations are identified with inequivalent irreducible representations of the corresponding C ^*-algebra. Topological terms in the action (or Hamiltonian) are found to be representation-dependent, and are automatically induced by the quantization procedure. Known charge quantization conditions turn out to be identically satisfied. Several examples are considered, among them the Dirac monopole and the Aharonov-Bohm effect.     

Synchronization of switched system and application in communication

In this Letter, synchronization of switched system is investigated based on Lyapunov method. A sufficient condition is derived to ensure the synchronization between two switched systems, and a new communication scheme is also proposed based on this. Furthermore, some secure analysis works, such as return map attack and moving average error attack, are also given to show the security of the proposed scheme. Finally, simulation examples are constructed to verify the theoretical analysis and its application for communication.     

On the universal covering group of the real symplectic group

A model for the universal covering group of the symplectic group as a Lie group, and some calculations based on the model, as well as defining a similar model for the Lagrangian Grassmannian and relating our construction to the Maslov Index.     

Geometric quantization of an OSp(1/2) coadjoint orbit

In this Letter, we show how the complete geometric quantization extends to specific supersymplectic supermanifolds. More precisely, we extend this procedure to OSp(1/2)-coadjoint orbits, which are graded extensions of elliptic Sp(2, )-coadjoint orbits. Our approach exploits results obtained in a previous work, where the notion of a super-Kähler supermanifold was defined, and the former orbits were shown to be nontrivial examples of such a notion. As their underlying Kähler manifolds, these supermanifolds carry a natural (super-Kähler) polarization, a crucial notion that was so far lacking. Geometric quantization leads here to a nontrivial representation of osp(1/2), which is realized in a space of square integrable holomorphic sections of a super-Hermitian complex line bundle sheaf-with-connection over the homogenous space OSp(1/2)/U(1).