A cut finite element method for a Stokes interface problem

We present a finite element method for the Stokes equations involving two immiscible incompressible fluids with different viscosities and with surface tension. The interface separating the two fluids does not need to align with the mesh. We propose a Nitsche formulation which allows for discontinuities along the interface with optimal a priori error estimates. A stabilization procedure is included which ensures that the method produces a well conditioned stiffness matrix independent of the location of the interface.     

The first eigenvalue of the p‐Laplace operator under powers of mean curvature flow

In this paper, we show the following main results. Let (Mⁿ,g(t)), t ∈ [0,T), be a solution of the unnormalized H^k ? flow on a closed manifold, and λ_1,p(t) be the first eigenvalue of the p‐Laplace operator. If there exists a nonnegative constant ε such that in M × [0,T) and in M × [0,T),then λ_1,p(t) is increasing and the differentiable almost everywhere along the unnormalized H^k ? flow on [0,T). At last, we discuss some useful monotonic quantities. Copyright © 2013 John Wiley & Sons, Ltd.     

Analysis of a combined mixed finite element and discontinuous Galerkin method for incompressible two‐phase flow in porous media

We analyze a combined method consisting of the mixed finite element method for pressure equation and the discontinuous Galerkin method for saturation equation for the coupled system of incompressible two‐phase flow in porous media. The existence and uniqueness of numerical solutions are established under proper conditions by using a constructive approach. Optimal error estimates in L²(H¹) for saturation and in L ^∞ (H(div)) for velocity are derived. Copyright © 2013 John Wiley & Sons, Ltd.     

Modeling,simulation and validation of material flow on conveyor belts

In this paper a model comparison approach based on material flow systems is investigated that is divided into a microscopic and a macroscopic model scale. On the microscopic model scale particles are simulated using a model based on Newton dynamics borrowed from the engineering literature. Phenomenological observations lead to a hyperbolic partial differential equation on the macroscopic model scale. Suitable numerical algorithms are presented and both models are compared numerically and validated against real-data test settings.     

A new kernel function for SPH with applications to free surface flows

Smoothed particle hydrodynamics (SPH) is a popular meshfree Lagrangian particle method, which uses a kernel function for numerical approximations. The kernel function is closely related to the computational accuracy and stability of the SPH method. In this paper, a new kernel function is proposed, which consists of two cosine functions and is referred to as double cosine kernel function. The newly proposed double cosine kernel function is sufficiently smooth, and is associated with an adjustable support domain. It also has smaller second order momentum, and therefore it can have better accuracy in terms of kernel approximation. SPH method with this double cosine kernel function is applied to simulate a dam-break flow and water entry of a horizontal circular cylinder. The obtained SPH results agree very well with the experimental results. The double cosine kernel function is also comparatively studied with two frequently used SPH kernel functions, Gaussian and cubic spline kernel functions.     

Scheduling the hospital-wide flow of elective patients

In this paper, we address the problem of planning the patient flow in hospitals subject to scarce medical resources with the objective of maximizing the contribution margin. We assume that we can classify a large enough percentage of elective patients according to their diagnosis-related group (DRG) and clinical pathway. The clinical pathway defines the procedures (such as different types of diagnostic activities and surgery) as well as the sequence in which they have to be applied to the patient. The decision is then on which day each procedure of each patient’s clinical pathway should be done, taking into account the sequence of procedures as well as scarce clinical resources, such that the contribution margin of all patients is maximized. We develop two mixed-integer programs (MIP) for this problem which are embedded in a static and a rolling horizon planning approach. Computational results on real-world data show that employing the MIPs leads to a significant improvement of the contribution margin compared to the contribution margin obtained by employing the planning approach currently practiced. Furthermore, we show that the time between admission and surgery is significantly reduced by applying our models.     

Asymmetric flow networks

This paper provides a new model of network formation that bridges the gap between the two benchmark game-theoretic models by Bala and Goyal (2000a) – the one-way flow model, and the two-way flow model – and includes both as limiting cases. As in both the said models, a link can be initiated unilaterally by any player with any other in what we call an “asymmetric flow” network, and the flow through a link towards the player who supports it is perfect. Unlike those models, there is friction or decay in the opposite direction. When this decay is complete there is no flow and this corresponds to the one-way flow model. The limit case when the decay in the opposite direction (and asymmetry) disappears corresponds to the two-way flow model. We characterize stable and strictly stable architectures for the whole range of parameters of this “intermediate” and more general model. A study of the efficiency of these architectures shows that in general stability and efficiency do not go together. We also prove the convergence of Bala and Goyal’s dynamic model in this context.     

Boundary of Anosov dynamics and evolution equations for surfaces

We show that a compact surface of genus greater than one, without focal points and a finite number of bubbles (“good” shaped regions of positive curvature) is in the closure of Anosov metrics. Compact surfaces of nonpositive curvature and genus greater than one are in the closure of Anosov metrics, by Hamilton's work about the Ricci flow. We generalize this fact to the above surfaces without focal points admitting regions of positive curvature using a “magnetic” version of the Ricci flow, the so‐called Ricci Yang‐Mills flow.