Shape inverse problem for Stokes–Brinkmann equations

In this paper, we propose an imaging technique for the detection of porous inclusions in a stationary flow governed by Stokes–Brinkmann equations. We introduce the velocity method to perform the shape deformation, and derive the structure of shape gradient for the cost functional based on the continuous adjoint method and the function space parametrization technique. Moreover, we present a gradient-type algorithm to the shape inverse problem. The numerical results demonstrate the proposed algorithm is feasible and effective for the quite high Reynolds numbers problems.     

一个新的边界层方程及其在形状控制中的应用

本文给出固壁边界上(即一个二维流形上) 的流体速度梯度和压力的二阶偏微分方程, 从而也给出边界上法向应力, 以及流体中运动物体所受的阻力和升力的计算公式. 本方法的创新在于边界上法向速度梯度不是通过在边界层内速度梯度的数值微分达到, 而是通过它与其他变量一起作为一组偏微分方程的解而得到, 证明边界层方程组的适定性问题, 并且给出解关于边界形状的Gâteaux 导数所满足的偏微分方程. 本文将本方法应用于飞机外形的形状最优控制, 给出阻力泛函关于形状第一变分的可计算形式. 数值例子表明, 用本方法得到的阻力精度比通用程序得到要高.     

The application of adjoint method for shape optimization in Stokes–Oseen flow

This paper presents a numerical method for shape optimization of a body immersed in an incompressible viscous flow governed by Stokes–Oseen equations. The purpose of this work is to optimize the shape that minimizes a given cost functional. Based on the continuous adjoint method, the shape gradient of the cost functional is derived by involving a Lagrangian functional with the function space parametrization technique. Then, a gradient‐type algorithm is applied to the shape optimization problem. The numerical examples indicate the proposed algorithm is feasible and effective in low Reynolds number flow. Copyright © 2016 John Wiley & Sons, Ltd.     

Cross-Diffusion Systems with Excluded-Volume Effects and Asymptotic Gradient Flow Structures

In this paper, we discuss the analysis of a cross-diffusion PDE system for a mixture of hard spheres, which was derived in Bruna and Chapman (J Chem Phys 137:204116-1–204116-16, 2012a) from a stochastic system of interacting Brownian particles using the method of matched asymptotic expansions. The resulting cross-diffusion system is valid in the limit of small volume fraction of particles. While the system has a gradient flow structure in the symmetric case of all particles having the same size and diffusivity, this is not valid in general. We discuss local stability and global existence for the symmetric case using the gradient flow structure and entropy variable techniques. For the general case, we introduce the concept of an asymptotic gradient flow structure and show how it can be used to study the behavior close to equilibrium. Finally, we illustrate the behavior of the model with various numerical simulations.     

Shape optimization of a body in the Oseen flow

This article is concerned with a numerical simulation of shape optimization of the Oseen flow around a solid body. The shape gradient for shape optimization problem in a viscous incompressible flow is computed by the velocity method. The flow is governed by the Oseen equations with mixed boundary conditions containing the pressure. The structure of continuous shape gradient of the cost functional is derived by using the differentiability of a minimax formulation involving a Lagrange functional with a function space parametrization technique. A gradient type algorithm is applied to the shape optimization problem. Numerical examples show that our theory is useful for practical purpose and the proposed algorithm is feasible. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Shape Optimization of a Nozzle with Specified Flow Field Including Viscosity Effect

In this article we intend to find the optimal shape of a nozzle respecting to some given target flow fields including viscosity effect. Via an approach based on measure theory which is not an iterative method and need not to any initial guess, each shape optimization problems are solved and consequently each geometry of the nozzle corresponding to prescribed flow fields is determined. Analyzing several case studies make us to confident on the use of the presented approach, because the obtained results give entirely the same as what we expect physically.     

Shape sensitivity analysis of an elastic contact problem: Convergence of the Nitsche based finite element approximation

In a recent work, we introduced a finite element approximation for the shape optimization of an elastic structure in sliding contact with a rigid foundation where the contact condition (Signorini’s condition) is approximated by Nitsche’s method and the shape gradient is obtained via the adjoint state method. The motivation of this work is to propose an a priori convergence analysis of the numerical approximation of the variables of the shape gradient (displacement and adjoint state) and to show some numerical results in agreement with the theoretical ones. The main difficulty comes from the non-differentiability of the contact condition in the classical sense which requires the notion of conical differentiability.     

On Distributed $H^1$ Shape Gradient Flows in Optimal Shape Design of Stokes Flows: Convergence Analysis and Numerical Applications

We consider optimal shape design in Stokes flow using $H^1$ shape gradient flows based on the distributed Eulerian derivatives. MINI element is used for discretizations of Stokes equation and Galerkin finite element is used for discretizations of distributed and boundary $H^1$ shape gradient flows. Convergence analysis with a priori error estimates is provided under general and different regularity assumptions. We investigate the performances of shape gradient descent algorithms for energy dissipation minimization and obstacle flow. Numerical comparisons in 2D and 3D show that the distributed $H^1$ shape gradient flow is more accurate than the popular boundary type. The corresponding distributed shape gradient algorithm is more effective.     

On a flow to untangle elastic knots

In this paper we present a method to untangle smooth knots by a gradient flow for a suitable energy. We show that the flow of smooth initial knots remains smooth for all time and approaches asymptotically an “optimal embedding” in its isotopy type. The method is to set up a gradient flow for the total energy of knots, which consists of bending energy and the Möbius energy of knots.     

The dynamics of elastic closed curves under uniform high pressure

We consider the dynamics of an inextensible elastic closed wire in the plane under uniform high pressure. In 1967, Tadjbakhsh and Odeh (J. Math. Anal. Appl. 18:59–74, 1967) posed a variational problem to determine the shape of a buckled elastic ring under uniform pressure. In order to comprehend a dynamics of the wire, we consider the following two mathematical questions: (i) can we construct a gradient flow for the Tadjbakhsh–Odeh functional under the inextensibility condition?; (ii) what is a behavior of the wire governed by the gradient flow near every critical point of the Tadjbakhsh–Odeh variational problem? For (i), first we derive a system of equations which governs the gradient flow, and then, give an affirmative answer to (i) by solving the system involving fourth order parabolic equations. For (ii), we first prove a stability and instability of each critical point by considering the second variation formula of the Tadjbakhsh–Odeh functional. Moreover, we give a lower bound of its Morse index. Finally we prove a dynamical aspects of the wire near each equilibrium state.     

A numerical method for mass conservative coupling between fluid flow and solute transport

We present a new coupled discretization approach for species transport in an incompressible fluid. The Navier-Stokes equations for the flow are discretized by the divergence-free Scott-Vogelius element on barycentrically refined meshes guaranteeing LBB stability. The convection-diffusion equation for species transport is discretized by the Voronoi finite volume method. In accordance to the continuous setting, due to the exact integration of the normal component of the flow through the Voronoi surfaces, the species concentration fulfills discrete global and local maximum principles. Besides of the numerical scheme itself, we present important aspects of its implementation. Further, for the case of homogeneous Dirichlet boundary conditions, we give a convergence proof for the coupled scheme. We report results of the application of the scheme to the interpretation of limiting current measurements in an electrochemical flow cell with cylindrical shape.     

Affine Invariant Detection: Edge Maps, Anisotropic Diffusion, and Active Contours

In this paper we undertake a systematic investigation of affine invariant object detection and image denoising. Edge detection is first presented from the point of view of the affine invariant scale-space obtained by curvature based motion of the image level-sets. In this case, affine invariant maps are derived as a weighted difference of images at different scales. We then introduce the affine gradient as an affine invariant differential function of lowest possible order with qualitative behavior similar to the Euclidean gradient magnitude. These edge detectors are the basis for the extension of the affine invariant scale-space to a complete affine flow for image denoising and simplification, and to define affine invariant active contours for object detection and edge integration. The active contours are obtained as a gradient flow in a conformally Euclidean space defined by the image on which the object is to be detected. That is, we show that objects can be segmented in an affine invariant manner by computing a path of minimal weighted affine distance, the weight being given by functions of the affine edge detectors. The gradient path is computed via an algorithm which allows to simultaneously detect any number of objects independently of the initial curve topology. Based on the same theory of affine invariant gradient flows we show that the affine geometric heat flow is minimizing, in an affine invariant form, the area enclosed by the curve.     

A note on minimization problems and multistep methods

Summary. We use the qualitative properties of the solution flow of the gradient equation to compute a local minimum of a real-valued function . Under the regularity assumption of all equilibria we show a convergence result for bounded trajectories of a consistent, strictly stable linear multistep method applied to the gradient equation. Moreover, we compare the asymptotic features of the numerical and the exact solutions as done by Humphries, Stuart (1994) and Schropp (1995) for one-step methods. In the case of -stable formulae this leads to an efficient solver for stiff minimization problems. Received July 10, 1995 / Revised version received June 27, 1996     

Analytical solutions to Stokes-type flows of inhomogeneous fluids

In this paper, we study the unsteady motion of an inhomogeneous incompressible viscous fluid, where the viscosity varies spatially according to various models. We study the Stokes-type flow for these types of fluids where in the first case the flow between two parallel plates is examined with one of the plates oscillating and in the second case when the flow is caused by a pulsatile pressure gradient. A general argument establishes the existence of oscillatory solutions to our problem. Exact solutions are obtained in terms of some special functions and comparisons are made with the cases of constant viscosity and the slow flow regimes.     

On the Voitkunskii Amfilokhiev Pavlovskii Model of Motion of Aqueous Polymer Solutions

We study the mathematical properties of the model of motion of aqueous polymer solutions (Voitkunskii, Amfilokhiev, Pavlovskii, 1970) and its modifications in the limiting case of small relaxation times (Pavlovskii, 1971). In both cases, we examine plane unsteady laminar flows. In the first case, the properties of the flows are similar to those of the flow of an ordinary viscous fluid. In the second case, there may exist weak discontinuities that are preserved during the motion. We also address the steady flow problem for a dilute aqueous polymer solution moving in a cylindrical tube under a longitudinal pressure gradient. In this case, a flow with rectilinear trajectories (an analog of the classical Poiseuille flow) is possible. However, in contrast to the latter, the pressure in this flow depends on all three spatial variables.     

Computational inversion of electron micrographs using L2‐gradient flows—convergence analysis

A gradient flow‐based explicit finite element method (L2GF) for reconstructing the 3D density function from a set of 2D electron micrographs has been proposed in recently published papers. The experimental results showed that the proposed method was superior to the other classical algorithms, especially for the highly noisy data. However, convergence analysis of the L2GF method has not been conducted. In this paper, we present a complete analysis on the convergence of L2GF method for the case of using a more general form regularization term, which includes the Tikhonov‐type regularizer and modified or smoothed total variation regularizer as two special cases. We further prove that the L²‐gradient flow method is stable and robust. These results demonstrate that the iterative variational reconstruction method derived from the L²‐gradient flow approach is mathematically sound and effective and has desirable properties. Copyright © 2013 John Wiley & Sons, Ltd.     

Dissociation effect on the hypersonic flow past a circular cylinder

The effects of dissociation of air on hypersonic flow past a circular cylinder at zero angle of incidence are considered under the assumptions that the shock wave is in the shape of a circular cylinder, the density ratio across the shock is constant, the flow behind the shock is at constant density and dissociation occurs only behind the shock wave. In the present paper, the velocity, pressure and drag coefficients, vorticity, shock detachment distance, stagnation point velocity gradient and sonic points on the shock and the surface have been obtained in the presence of dissociation. The results have been compared with the corresponding results obtained in the case when dissociation dose not occur and the corresponding results in the case of the sphere in the presence of dissociation.     

Identification of a free boundary in Norton–Hoff flows with thermal effects

This work deals with a free boundary identification problem in a steady viscoplastic flow. We provide a novel identification model based on a non-linear optimization. The fluid motion is governed by the incompressible Norton–Hoff model coupled with the heat equation. The viscosity of the fluid is modeled by the non-linear Arrhenius law. Our point of view is to treat the problem as a shape sensitivity of a cost functional formulated on the free boundary and governed by the normal component of the velocity of the flow. We analyze the mathematical statement of the forward problem. The equations related to the free boundary are simplified. Various properties of this optimization are proved. Since the state of Norton–Hoff model is not regular enough we introduce a parameter penalization. The shape gradient of the considered cost functional is given in the strong sense up to the parameter of penalization. We supply the expression of the shape gradient in a weak sense.     

Multiple mixed states of nodal solutions for nonlinear Schrödinger systems

In this paper we develop a general critical point theory to deal with existence and locations of multiple critical points produced by minimax methods in relation to multiple invariant sets of the associated gradient flow. The motivation is to study non-trivial nodal solutions with each component sign-changing for a class of nonlinear Schrödinger systems which arise from Bose–Einstein condensates theory. Our general method allows us to obtain infinitely many mixed states of nodal solutions for the repulsive case.     

A two-step extragradient method for variational inequalities

In this paper we consider an extragradient method for solving variational inequalities and related problems. On each iteration this method makes two trial steps along the gradient, and the value of the gradient at the second point is used at the first point as the iteration direction. We prove the convergence of this method in a general case. For problems with a bilinear functional we prove the geometric convergence rate.