Weak solutions to problems involving inviscid fluids

We consider an abstract functional-differential equation derived from the pressure-less Euler system with variable coefficients that includes several systems of partial differential equations arising in the fluid mechanics. Using the method of convex integration we show the existence of infinitely many weak solutions for prescribed initial data and kinetic energy

Critical dynamics of an isothermal compressible non-ideal fluid

A pure fluid at its critical point shows a dramatic slow-down in its dynamics, due to a divergence of the order-parameter susceptibility and the coefficient of heat transport. Under isothermal conditions, however, sound waves provide the only possible relaxation mechanism for order-parameter fluctuations. Here we study the critical dynamics of an isothermal, compressible non-ideal fluid via scaling arguments and computer simulations of the corresponding fluctuating hydrodynamics equations. We show that, below a critical dimension of 4, the order-parameter dynamics of an isothermal fluid effectively reduces to "model A," characterized by overdamped sound waves and a divergent bulk viscosity. In contrast, the shear viscosity remains finite above two dimensions. Possible applications of the model are discussed.Comment: 19 pages, 7 figures; v3: minor corrections and clarifications; as published in Phys. Rev.

Quantum dot formation on a strain-patterned epitaxial thin film

We model the effect of substrate strain patterning on the self-assembly of quantum dots (QDs). When the surface energy is isotropic, we demonstrate that strain patterning via embedded substrate inclusions may result in ordered, self-organized QD arrays. However, for systems with strong cubic surface energy anisotropy, the same patterning does not readily lead to an ordered array of pyramids at long times. We conclude that the form of the surface energy anisotropy strongly influences the manner in which QDs self-assemble into regular arrays.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/87827/2/133102_1.pd

On a Cahn--Hilliard--Darcy system for tumour growth with solution dependent source terms

We study the existence of weak solutions to a mixture model for tumour growth that consists of a Cahn--Hilliard--Darcy system coupled with an elliptic reaction-diffusion equation. The Darcy law gives rise to an elliptic equation for the pressure that is coupled to the convective Cahn--Hilliard equation through convective and source terms. Both Dirichlet and Robin boundary conditions are considered for the pressure variable, which allows for the source terms to be dependent on the solution variables.Comment: 18 pages, changed proof from fixed point argument to Galerkin approximatio

Thermodynamics of non-local materials: extra fluxes and internal powers

The most usual formulation of the Laws of Thermodynamics turns out to be suitable for local or simple materials, while for non-local systems there are two different ways: either modify this usual formulation by introducing suitable extra fluxes or express the Laws of Thermodynamics in terms of internal powers directly, as we propose in this paper. The first choice is subject to the criticism that the vector fluxes must be introduced a posteriori in order to obtain the compatibility with the Laws of Thermodynamics. On the contrary, the formulation in terms of internal powers is more general, because it is a priori defined on the basis of the constitutive equations. Besides it allows to highlight, without ambiguity, the contribution of the internal powers in the variation of the thermodynamic potentials. Finally, in this paper, we consider some examples of non-local materials and derive the proper expressions of their internal powers from the power balance laws.Comment: 16 pages, in press on Continuum Mechanics and Thermodynamic

Mass conservative and energy stable finite difference methods for the quasi-incompressible Navier–Stokes–Cahn–Hilliard system:Primitive variable and projection-type schemes

In this paper we describe two fully mass conservative, energy stable, finite difference methods on a staggered grid for the quasi-incompressible Navier-Stokes-Cahn-Hilliard (q-NSCH) system governing a binary incompressible fluid flow with variable density and viscosity. Both methods, namely the primitive method (finite difference method in the primitive variable formulation) and the projection method (finite difference method in a projection-type formulation), are so designed that the mass of the binary fluid is preserved, and the energy of the system equations is always non-increasing in time at the fully discrete level. We also present an efficient, practical nonlinear multigrid method - comprised of a standard FAS method for the Cahn-Hilliard equation, and a method based on the Vanka-type smoothing strategy for the Navier-Stokes equation - for solving these equations. We test the scheme in the context of Capillary Waves, rising droplets and Rayleigh-Taylor instability. Quantitative comparisons are made with existing analytical solutions or previous numerical results that validate the accuracy of our numerical schemes. Moreover, in all cases, mass of the single component and the binary fluid was conserved up to 10 to -8 and energy decreases in time

Extended Smoothed Boundary Method for Solving Partial Differential Equations with General Boundary Conditions on Complex Boundaries

In this article, we describe an approach for solving partial differential equations with general boundary conditions imposed on arbitrarily shaped boundaries. A continuous function, the domain parameter, is used to modify the original differential equations such that the equations are solved in the region where a domain parameter takes a specified value while boundary conditions are imposed on the region where the value of the domain parameter varies smoothly across a short distance. The mathematical derivations are straightforward and generically applicable to a wide variety of partial differential equations. To demonstrate the general applicability of the approach, we provide four examples herein: (1) the diffusion equation with both Neumann and Dirichlet boundary conditions; (2) the diffusion equation with both surface diffusion and reaction; (3) the mechanical equilibrium equation; and (4) the equation for phase transformation with the presence of additional boundaries. The solutions for several of these cases are validated against corresponding analytical and semi-analytical solutions. The potential of the approach is demonstrated with five applications: surface-reaction-diffusion kinetics with a complex geometry, Kirkendall-effect-induced deformation, thermal stress in a complex geometry, phase transformations affected by substrate surfaces, and a self-propelled droplet.Comment: This document is the revised version of arXiv:0912.1288v

A numerical method for the quasi-incompressible Cahn-Hilliard-Navier-Stokes equations for variable density flows with a discrete energy law

In this paper, we investigate numerically a diffuse interface model for the Navier-Stokes equation with fluid-fluid interface when the fluids have different densities \cite{Lowengrub1998}. Under minor reformulation of the system, we show that there is a continuous energy law underlying the system, assuming that all variables have reasonable regularities. It is shown in the literature that an energy law preserving method will perform better for multiphase problems. Thus for the reformulated system, we design a $C^0$ finite element method and a special temporal scheme where the energy law is preserved at the discrete level. Such a discrete energy law (almost the same as the continuous energy law) for this variable density two-phase flow model has never been established before with $C^0$ finite element. A Newton's method is introduced to linearise the highly non-linear system of our discretization scheme. Some numerical experiments are carried out using the adaptive mesh to investigate the scenario of coalescing and rising drops with differing density ratio. The snapshots for the evolution of the interface together with the adaptive mesh at different times are presented to show that the evolution, including the break-up/pinch-off of the drop, can be handled smoothly by our numerical scheme. The discrete energy functional for the system is examined to show that the energy law at the discrete level is preserved by our scheme

A new ghost cell/level set method for moving boundary problems:application to tumor growth

In this paper, we present a ghost cell/level set method for the evolution of interfaces whose normal velocity depend upon the solutions of linear and nonlinear quasi-steady reaction-diffusion equations with curvature-dependent boundary conditions. Our technique includes a ghost cell method that accurately discretizes normal derivative jump boundary conditions without smearing jumps in the tangential derivative; a new iterative method for solving linear and nonlinear quasi-steady reaction-diffusion equations; an adaptive discretization to compute the curvature and normal vectors; and a new discrete approximation to the Heaviside function. We present numerical examples that demonstrate better than 1.5-order convergence for problems where traditional ghost cell methods either fail to converge or attain at best sub-linear accuracy. We apply our techniques to a model of tumor growth in complex, heterogeneous tissues that consists of a nonlinear nutrient equation and a pressure equation with geometry-dependent jump boundary conditions. We simulate the growth of glioblastoma (an aggressive brain tumor) into a large, 1 cm square of brain tissue that includes heterogeneous nutrient delivery and varied biomechanical characteristics (white matter, gray matter, cerebrospinal fluid, and bone), and we observe growth morphologies that are highly dependent upon the variations of the tissue characteristics—an effect observed in real tumor growth

Polymer-induced tubulation in lipid vesicles

A mechanism of extraction of tubular membranes from a lipid vesicle is presented. A concentration gradient of anchoring amphiphilic polymers generates tubes from bud-like vesicle protrusions. We explain this mechanism in the framework of the Canham-Helfrich model. The energy profile is analytically calculated and a tube with a fixed length, corresponding to an energy minimum, is obtained in a certain regime of parameters. Further, using a phase-field model, we corroborate these results numerically. We obtain the growth of tubes when a polymer source is added, and the bud-like shape after removal of the polymer source, in accordance with recent experimental results