Computation of nonautonomous invariant and inertial manifolds

We derive a numerical scheme to compute invariant manifolds for time-variant discrete dynamical systems, i.e., nonautonomous difference equations. Our universally applicable method is based on a truncated Lyapunov–Perron operator and computes invariant manifolds using a system of nonlinear algebraic equations which can be solved both locally using (nonsmooth) inexact Newton, and globally using continuation algorithms. Compared to other algorithms, our approach is quite flexible, since it captures time-dependent, nonsmooth, noninvertible or implicit equations and enables us to tackle the full hierarchy of strongly stable, stable and center-stable manifolds, as well as their unstable counterparts. Our results are illustrated using a test example and are applied to a population dynamical model and the Hénon map. Finally, we discuss a linearly implicit Euler–Bubnov–Galerkin discretization of a reaction diffusion equation in order to approximate its inertial manifold.     

Persistence of undercompressive phase boundaries for isothermal Euler equations including configurational forces and surface tension

The persistence of subsonic phase boundaries in a multidimensional Van der Waals fluid is analyzed. The phase boundary is considered as a sharp free boundary that connects liquid and vapor bulk phase dynamics given by the isothermal Euler equations. The evolution of the boundary is driven by effects of configurational forces as well as surface tension. To analyze this problem, the equations and trace conditions are linearized such that one obtains a general hyperbolic initial boundary value problem with higher‐order boundary conditions. A global existence theorem for the linearized system with constant coefficients is shown. The proof relies on the normal mode analysis and a linear form in suitable spaces that is defined using an associated adjoint problem. Especially, the associated adjoint problem satisfies the uniform backward in time Kreiss–Lopatinski? condition. A new energy‐like estimate that also includes surface energy terms leads finally to the uniqueness and regularity for the found solutions of the problem in weighted spaces. Copyright © 2016 John Wiley & Sons, Ltd.     

Bed of polydisperse viscous spherical drops under thermocapillary effects

Viscous flow past an ensemble of polydisperse spherical drops is investigated under thermocapillary effects. We assume that the collection of spherical drops behaves as a porous media and estimates the hydrodynamic interactions analytically via the so- called cell model that is defined around a specific representative particle. In this method, the hydrodynamic interactions are assumed to be accounted by suitable boundary conditions on a fictitious fluid envelope surrounding the representative particle. The force calculated on this representative particle will then be extended to a bed of spherical drops visualized as a Darcy porous bed. Thus, the “effective bed permeability” of such a porous bed will be computed as a function of various parameters and then will be compared with Carman–Kozeny relation. We use cell model approach to a packed bed of spherical drops of uniform size (monodisperse spherical drops) and then extend the work for a packed bed of polydisperse spherical drops, for a specific parameters. Our results show a good agreement with the Carman–Kozeny relation for the case of monodisperse spherical drops. The prediction of overall bed permeability using our present model agrees well with the Carman–Kozeny relation when the packing size distribution is narrow, whereas a small deviation can be noted when the size distribution becomes broader.     

Phase-field approach to fracture for finite-deformation contact problems

The present contribution deals with a variationally consistent Mortar contact algorithm applied to a phase-field fracture approach for finite deformations, see [4]. A phase-field approach to fracture allows for the numerical simulation of complex fracture patterns for three dimensional problems, extended recently to finite deformations (see [2] for more details). In a nutshell, the phase-field approach relies on a regularization of the sharp (fracture-) interface. In order to improve the accuracy, a fourth-order Cahn-Hilliard phase-field equation is considered, requiring global C¹ continuity (see [1]), which will be dealt with using an isogeometrical analysis (IGA) framework. Additionally, a newly developed hierarchical refinement scheme is applied to resolve for local physical phenomena e.g. the contact zone (see [3] for more details). The Mortar method is a modern and very accurate numerical method to implement contact boundaries. This approach can be extended in a straightforward manner to transient phase-field fracture problems. The performance of the proposed methods will be examined in a representative numerical example. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Minimization- and Saddle-Point-Based Modeling of Diffusion-Deformation-Processes in Hydrogels

Hydrogels have gained importance during the last years due to their wide range of synthetically fabricable elastic properties as well their increasing meaning in biomedical applications. Future exploitation of the vast prospects of hydrogels is however only feasible by establishing reliable material models that precisely capture their behavior in different environments. To this end, we propose a consistent variational framework for deformation-diffusion processes, offering a canonically compact approach to the chemo-mechanical coupling of hydrogels via a saddle-point as well as a new minimization formulation. The work depicts the construction of rate-type potentials for the chemo-mechanical evolution problem and their transformation into time-discrete incremental potentials. In terms of spatial discretization, the finite element method is employed, benefiting from the intrinsic symmetric structure of the variational foundation. While the saddle-point formulation yields the well-known LBB condition as a constraint for finite element interpolations, on the part of its minimizing counterpart H(Div, ℬ︁)-conforming elements have to be chosen. We illustrate appropriate solutions to both challenges, using mixed Taylor-Hood for the saddle-point and Raviart-Thomas elements for the minimization formulation and discuss advantages of the new approach. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

An anisotropic phase-field model for transversely isotropic barium titanate with bounded moduli

This contribution presents a phase-field model for transversely isotropic barium titanate, which allows for the adjustment of the full set of anisotropic material parameters. It is a direct extension of the work by Schrade et al. [1] who proposed a phase-field model for ferroelectrics in the framework of invariant theory. In the present contribution, the loss of positive definiteness is avoided by formulating energetic terms that provide upper and lower bounds for all material moduli involved. We show the characteristics of the formulation by a set of numerical examples in two and three dimensions. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Magneto-Mechanically Coupled Phase-Field Model at Large Deformation

The aggregate magneto-mechanical behavior of magneto-rheological elastomers (MREs) stems from the magnetic properties of the ferromagnetic inclusion and the mechanical properties of the matrix material. We propose a large deformation micro-magnetic theory, to predict the behavior and interaction of ferromagnetic particles inside an elastomeric matrix. A rate-type variational principle, with the magnetization as the order parameter is proposed. A large deformation Landau-Lifshitz-Gilbert equation for the time evolution of the magnetization, is obtained directly from the proposed variational principle. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On Kemnitz’ conjecture concerning lattice-points in the plane

In 1961, Erdős, Ginzburg and Ziv proved a remarkable theorem stating that each set of 2n−1 integers contains a subset of size n, the sum of whose elements is divisible by n. We will prove a similar result for pairs of integers, i.e. planar lattice-points, usually referred to as Kemnitz’ conjecture. Dedicated to Richard Askey on the occasion of his 70th birthday. 2000 Mathematics Subject Classification Primary—11B50.     

Renormalization as a functor on bialgebras

The Hopf algebra of renormalization in quantum field theory is described at a general level. The products of fields at a point are assumed to form a bialgebra B and renormalization endows T(T(B)⁺), the double tensor algebra of B, with the structure of a noncommutative bialgebra. When the bialgebra B is commutative, renormalization turns S(S(B)⁺), the double symmetric algebra of B, into a commutative bialgebra. The usual Hopf algebra of renormalization is recovered when the elements of S¹(B) are not renormalized, i.e., when Feynman diagrams containing one single vertex are not renormalized. When B is the Hopf algebra of a commutative group, a homomorphism is established between the bialgebra S(S(B)⁺) and the Faà di Bruno bialgebra of composition of series. The relation with the Connes-Moscovici Hopf algebra is given. Finally, the bialgebra S(S(B)⁺) is shown to give the same results as the standard renormalization procedure for the scalar field.