A high order moving boundary treatment for compressible inviscid flows

We develop a high order numerical boundary condition for compressible inviscid flows involving complex moving geometries. It is based on finite difference methods on fixed Cartesian meshes which pose a challenge that the moving boundaries intersect the grid lines in an arbitrary fashion. Our method is an extension of the so-called inverse Lax–Wendroff procedure proposed in [17] for conservation laws in static geometries. This procedure helps us obtain normal spatial derivatives at inflow boundaries from Lagrangian time derivatives and tangential derivatives by repeated use of the Euler equations. Together with high order extrapolation at outflow boundaries, we can impose accurate values of ghost points near the boundaries by a Taylor expansion. To maintain high order accuracy in time, we need some special time matching technique at the two intermediate Runge–Kutta stages. Numerical examples in one and two dimensions show that our boundary treatment is high order accurate for problems with smooth solutions. Our method also performs well for problems involving interactions between shocks and moving rigid bodies.     

Physalis method for heterogeneous mixtures of dielectrics and conductors: Accurately simulating one million particles using a PC

Prosperetti’s seminal Physalis method, an Immersed Boundary/spectral method, had been used extensively to investigate fluid flows with suspended solid particles. Its underlying idea of creating a cage and using a spectral general analytical solution around a discontinuity in a surrounding field as a computational mechanism to enable the accommodation of physical and geometric discontinuities is a general concept, and can be applied to other problems of importance to physics, mechanics, and chemistry. In this paper we provide a foundation for the application of this approach to the determination of the distribution of electric charge in heterogeneous mixtures of dielectrics and conductors. The proposed Physalis method is remarkably accurate and efficient. In the method, a spectral analytical solution is used to tackle the discontinuity and thus the discontinuous boundary conditions at the interface of two media are satisfied exactly. Owing to the hybrid finite difference and spectral schemes, the method is spectrally accurate if the modes are not sufficiently resolved, while higher than second-order accurate if the modes are sufficiently resolved, for the solved potential field. Because of the features of the analytical solutions, the derivative quantities of importance, such as electric field, charge distribution, and force, have the same order of accuracy as the solved potential field during postprocessing. This is an important advantage of the Physalis method over other numerical methods involving interpolation, differentiation, and integration during postprocessing, which may significantly degrade the accuracy of the derivative quantities of importance. The analytical solutions enable the user to use relatively few mesh points to accurately represent the regions of discontinuity. In addition, the spectral convergence and a linear relationship between the cost of computer memory/computation and particle numbers results in a very efficient method. In the present paper, the accuracy of the method is numerically investigated by example computations using one dielectric particle, one isolated conductor particle, one conductor particle connected to an external source with imposed voltage, and two conductor/dielectric particles with strong interactions. The efficiency of the method is demonstrated with one million particles, which suggests that the method can be used for many important engineering applications of broad interest.     

A Cahn-Hilliard model in a domain with non-permeable walls

We consider in this article a Cahn-Hilliard model in a bounded domain with non-permeable walls, characterized by dynamic-type boundary conditions. Dynamic boundary conditions for the Cahn-Hilliard system have recently been proposed by physicists in order to account for the interactions with the walls in confined systems and are obtained by writing that the total bulk mass is conserved and that there is a relaxation dynamics on the boundary. However, in the case of non-permeable walls, one should also expect some mass on the boundary. It thus seems more realistic to assume that the total mass, in the bulk and on the boundary, is conserved, which leads to boundary conditions of a different type. For the resulting mathematical model, we prove the existence and uniqueness of weak solutions and study their asymptotic behavior as time goes to infinity.     

Instability of two-phase flows: A lower bound on the dimension of the global attractor of the Cahn-Hilliard-Navier-Stokes system

We consider a model for the flow of a mixture of two viscous and incompressible fluids in a two or three dimensional channel-like domain. The model consists of the Navier-Stokes equations governing the fluid velocity coupled with a convective Cahn-Hilliard equation for the relative density of atoms of one of the fluids. We prove the instability of certain stationary solutions for such a system endowed with periodic boundary conditions on elongated domains (0,2π/α₀)×(0,2π) or (0,2π/α₀)×(0,2π)×(0,2π/β₀) for a special class of periodic body forces, provided that α₀ and β₀ are small enough. As a consequence, we deduce a lower bound for the Hausdorff dimension of the global attractor.     

Finite temperature Casimir effect of massive fermionic fields in the presence of compact dimensions

We consider the finite temperature Casimir effect of a massive fermionic field confined between two parallel plates, with MIT bag boundary conditions on the plates. The background spacetime is Mp+1×Tq which has q dimensions compactified to a torus. On the compact dimensions, the field is assumed to satisfy periodicity boundary conditions with arbitrary phases. Both the high temperature and the low temperature expansions of the Casimir free energy and the force are derived explicitly. It is found that the Casimir force acting on the plates is always attractive at any temperature regardless of the boundary conditions assumed on the compact torus. The asymptotic limits of the Casimir force in the small plate separation limit are also obtained.     

Initial conditions and global event properties from color glass condensate

Perturbative unitarization from non-linear effects is thought to deplete the gluon density for transverse momenta below the saturation scale. Such effects also modify the distribution of gluons produced in heavy-ion collisions in transverse impact parameter space. I discuss some of the consequences for the initial conditions for hydrodynamic models of heavy-ion collisions and for hard “tomographic” probes. Also, I stress the importance of realistic modelling of the fluctuations of the valence sources for the small-x fields in the impact parameter plane. Such models can now be combined with solutions of running–coupling Balitsky–Kovchegov evolution to obtain controlled predictions for initial conditions at the LHC.     

Convergence results of Landweber iterations for linear systems

The Landweber scheme is a method for algebraic image reconstructions. The convergence behavior of the Landweber scheme is of both theoretical and practical importance. Using the diagonalization of matrix, we derive a neat iterative representation formula for the Landweber schemes and consequently establish the convergence conditions of Landweber iteration. This work refines our previous convergence results on the Landweber scheme.     

Spectral asymptotics and basis properties of fourth order differential operators with regular boundary conditions

In this paper, we consider the problem where λ is a spectral parameter; q(x) ∈ L₁(0,1) is complex‐valued function; α_s, s = 1,2,3, are arbitrary complex constants that satisfy α₂ = α₁ + α₃ and σ = 0,1. The boundary conditions of this problem are regular, but not strongly regular. Asymptotic formulae for eigenvalues and eigenfunctions of the considered boundary value problem are established. It is proved that all the eigenvalues, except for finite number, are simple and the system of root functions of this spectral problem forms a basis in the space L_p(0,1), 1 < p < ∞ , when ; moreover, this basis is unconditional for p = 2. We note that the considered problem was previously investigated in the condition of α₂ ≠ α₁ + α₃. Copyright © 2013 John Wiley & Sons, Ltd.     

On a constructive approach to optimality conditions for convex SIP problems with polyhedral index sets

In the paper, we consider a problem of convex Semi-Infinite Programming with an infinite index set in the form of a convex polyhedron. In study of this problem, we apply the approach suggested in our recent paper [Kostyukova OI, Tchemisova TV. Sufficient optimality conditions for convex Semi Infinite Programming. Optim. Methods Softw. 2010;25:279–297], and based on the notions of immobile indices and their immobility orders. The main result of the paper consists in explicit optimality conditions that do not use constraint qualifications and have the form of criterion. The comparison of the new optimality conditions with other known results is provided.     

Equivalent conditions on periodic feedback stabilization for linear periodic evolution equations

This paper studies the periodic feedback stabilization for a class of linear T  -periodic evolution equations. Several equivalent conditions on the linear periodic feedback stabilization are obtained. These conditions are related to the following subjects: the attainable subspace of the controlled evolution equation under consideration; the unstable subspace (of the evolution equation with the null control) provided by the Kato projection; the Poincaré map associated with the evolution equation with the null control; and two unique continuation properties for the dual equations on different time horizons [0,T]$[0,T]$ and [0,_n0T]$[0,n_{0}T]$ (where _n0$n_0$ is the sum of algebraic multiplicities of distinct unstable eigenvalues of the Poincaré map). It is also proved that a T-periodic controlled evolution equation is linear T-periodic feedback stabilizable if and only if it is linear T-periodic feedback stabilizable with respect to a finite-dimensional subspace. Some applications to heat equations with time-periodic potentials are presented.