Search

We consider a construction proposed in \cite{acharyaQAM} that builds on the notion of weak solutions for incompressible fluids to provide a scheme that generates variationally a certain type of dual solutions. If these dual solutions are regular enough one can use them to recover standard solutions. The scheme provides a generalisation of a construction of Y. Brenier for the Euler equations. We rigorously analyze the scheme, extending the work of Y.Brenier for Euler, and also provide an extension of it to the case of the Navier-Stokes equations. Furthermore we obtain the inviscid limit of Navier-Stokes to Euler as a $\Gamma$-limit.

We consider the motion of a large number of heavy particles in a Newtonian fluid occupying a bounded spatial domain. When we say "heavy", we mean a particle with a mass density that approaches infinity at an appropriate rate as its radius vanishes. We show that the collective effect of heavy particles on the fluid motion is similar to the Brinkman perturbation of the Navier-Stokes system identified in the homogenization process.

We show that the collective effect of $N$ rigid bodies $(\mathcal{S}_{n,N})_{n=1}^N$ of diameters $(r_{n,N})_{n=1}^N$ immersed in an incompressible non--Newtonian fluid is negligible in the asymptotic limit $N \to \infty$ as long as their total packing volume $\sum_{n=1}^N r_{n,N}^d$, $d=2,3$ tends to zero exponentially -- $\sum_{n=1}^N r_{n,N}^d \approx A^{-N}$ -- for a certain constant $A > 1$. The result is rather surprising and in a sharp contrast with the associated homogenization problem, where the same number of obstacles can completely stop the fluid motion in the case of shear thickening viscosity. A large class of non--Newtonian fluids is included, for which the viscous stress is a subdifferential of a convex potential.

This paper analytically explores a simplified model for the hydrodynamics of nematic liquid crystal colloids. We integrate a Stokes equation for the velocity field with a Ginzburg-Landau transported heat flow for the director field. The study focuses on a bounded spatial domain containing periodically distributed colloidal particles, which impose no-anchoring conditions on the nematic liquid crystal. By progressively reducing the particle size to zero and simultaneously increasing the number of particles, we delve into the associated homogenisation problem. Our analysis uncovers a form of decoupling where the velocity field asymptotically satisfies a Darcy equation, independent of the director, while the director follows a gradient flow, unaffected by the velocity field. One of the most intricate aspects of the homogenisation process is the absence of an extension operator for the director field that preserves the uniform estimates related to the system's energy. We address this challenge with a novel variation of the Aubin-Lions lemma, specifically adapted for homogenisation problems.

Rational computer-aided design of multiphase polymer materials is vital for rapid progress in many important applications, such as: diagnostic tests, drug delivery, coatings, additives for constructing materials, cosmetics, etc. Several property predictive models, including the prospective Population Balance Model for Latex Particles Morphology Formation (LPMF PBM), have already been developed for such materials. However, they lack computational efficiency, and the accurate prediction of materials' properties still remains a great challenge. To enhance performance of the LPMF PBM, we explore the feasibility of reducing its complexity through disregard of the aggregation terms of the model. The introduced nondimensionalization approach, which we call Optimal Scaling with Constraints, suggests a quantitative criterion for locating regions of slow and fast aggregation and helps to derive a family of dimensionless LPMF PBM of reduced complexity. The mathematical analysis of this new family is also provided. When compared with the original LPMF PBM, the resulting models demonstrate several orders of magnitude better computational efficiency.

We consider a system of colloidal particles embedded in a paranematic -- an isotropic phase of a nematogenic medium above the temperature of the nematic-to-isotropic transition. In this state, the nematic order is induced by the boundary conditions in a narrow band around each particle and it decays exponentially in the bulk. We develop rigorous asymptotics of the linearization of the appropriate variational model that allow us to describe weak far-field interactions between the colloidal particles in two dimensional paranematic suspensions. We demonstrate analytically that decay rates of solutions to the full nonlinear and linear problems are similar and verify numerically that the interactions between the particles in these problems have similar dependence on the distance between the particles. Finally, we perform Monte-Carlo simulations for a system of colloidal particles in a paranematic and describe the statistical properties of this system., Comment: 50 pages, 13 figures. Submitted

The paper concerns the analysis of global minimizers of a Dirichlet-type energy functional defined on the space of vector fields $H^1(S,T)$, where $S$ and $T$ are surfaces of revolution. The energy functional we consider is closely related to a reduced model in the variational theory of micromagnetism for the analysis of observable magnetization states in curved thin films. We show that axially symmetric minimizers always exist, and if the target surface $T$ is never flat, then any coexisting minimizer must have line symmetry. Thus, the minimization problem reduces to the computation of an optimal one-dimensional profile. We also provide a necessary and sufficient condition for energy minimizers to be axially symmetric.

We study the scaled anisotropic co-rotational Beris-Edwards system modeling the hydrodynamic motion of nematic liquid crystals in dimension two. We prove the global well-posedness with small analytic data in a thin strip domain. Moreover, we justify the limit to a system involving the hydrostatic Navier-Stokes system with analytic data and prove the convergence., Comment: second, revised version

We consider the motion of $N$ rigid bodies -- compact sets $(\mathcal{S}^1_\varepsilon, \cdots, \mathcal{S}^N_\varepsilon )_{\varepsilon > 0}$ -- immersed in a viscous incompressible fluid contained in a domain in the Euclidean space $\mathbb{R}^d$, $d=2,3$. We show the fluid flow is not influenced by the presence of the infinitely many bodies in the asymptotic limit $\varepsilon \to 0$ and $N=N(\varepsilon)\rightarrow\infty$ as soon as \[ {\rm diam}[\mathcal{S}^i_\varepsilon ] \to 0 \ \mbox{as}\ \varepsilon \to 0 ,\ i=1,\cdots, N(\varepsilon). \] The result depends solely on the geometry of the bodies and is independent of their mass densities. Collisions are allowed and the initial data are arbitrary with finite energy.

In this paper we study probabilistic and neural network approximations for solutions to Poisson equation subject to Holder data in general bounded domains of $\mathbb{R}^d$. We aim at two fundamental goals. The first, and the most important, we show that the solution to Poisson equation can be numerically approximated in the sup-norm by Monte Carlo methods, and that this can be done highly efficiently if we use a modified version of the walk on spheres algorithm as an acceleration method. This provides estimates which are efficient with respect to the prescribed approximation error and with polynomial complexity in the dimension and the reciprocal of the error. A crucial feature is that the overall number of samples does not not depend on the point at which the approximation is performed. As a second goal, we show that the obtained Monte Carlo solver renders in a constructive way ReLU deep neural network (DNN) solutions to Poisson problem, whose sizes depend at most polynomialy in the dimension $d$ and in the desired error. In fact we show that the random DNN provides with high probability a small approximation error and low polynomial complexity in the dimension.

We consider the motion of a small rigid object immersed in a viscous compressible fluid in the 3-dimensional Eucleidean space. Assuming the object is a ball of a small radius $\varepsilon$ we show that the behavior of the fluid is not influenced by the object in the asymptotic limit $\varepsilon \to 0$. The result holds for the isentropic pressure law $p(\varrho) = a \varrho^\gamma$ for any $\gamma > \frac{3}{2}$ under mild assumptions concerning the rigid body density. In particular, the latter may be bounded as soon as $\gamma > 3$. The proof uses a new method of construction of the test functions in the weak formulation of the problem, and, in particular, a new form of the so-called Bogovskii operator.

For the Landau-de Gennes functional on 3D domains, \begin{equation*} I_{\varepsilon}(Q,\Omega):=\int_{\Omega}\left\{\frac{1}{2}|\nabla Q|^2+\frac{1}{\varepsilon^2}\left( -\frac{a^2}{2}\mathrm{tr}(Q^2)-\frac{b^2}{3}\mathrm{tr}(Q^3)+\frac{c^2}{4}[\mathrm{tr}(Q^2)]^2 \right) \right\}\,dx, \end{equation*} it is well-known that under suitable boundary conditions, the global minimizer $Q_\varepsilon$ converges strongly in $H^1(\Omega)$ to a uniaxial minimizer $Q_*=s_+(n_*\otimes n_*-\frac{1}{3}\mathrm{Id})$ up to some subsequence $\varepsilon_n\rightarrow\infty$ , where $n_*\in H^1(\Omega,\mathbb{S}^2)$ is a minimizing harmonic map. In this paper we further investigate the structure of $Q_{\varepsilon}$ near the core of a point defect $x_0$ which is a singular point of the map $n_*$. The main strategy is to study the blow-up profile of $Q_{\varepsilon_n}(x_n+\varepsilon_n y)$ where $\{x_n\}$ are carefully chosen and converge to $x_0$. We prove that $Q_{\varepsilon_n}(x_n+\varepsilon_n y)$ converges in $C^2_{loc}(\mathbb{R}^n)$ to a tangent map $Q(x)$ which at infinity behaves like a "hedgehog" solution that coincides with the asymptotic profile of $n_*$ near $x_0$. Moreover, such convergence result implies that the minimizer $Q_{\varepsilon_n}$ can be well approximated by the Oseen-Frank minimizer $n_*$ outside the $O(\varepsilon_n)$ neighborhood of the point defect., Comment: 26 pages

We consider the motion of a compressible viscous fluid containing a moving rigid body confined to a planar domain $\Omega \subset R^2$. The main result states that the influence of the body on the fluid is negligible if (i) the diameter of the body is small and (ii) the fluid is nearly incompressible (the low Mach number regime). The specific shape of the body as well as the boundary conditions on the fluid--body interface are irrelevant and collisions with the boundary $\partial \Omega$ are allowed. The rigid body motion may be enforced externally or governed solely by its interaction with the fluid.

We study globally bounded entire minimizers $u:\mathbb{R}^n\rightarrow\mathbb{R}^m$ of Allen-Cahn systems for potentials $W\geq 0$ with $\{W=0\}=\{a_1,...,a_N\}$ and $W(u)\sim |u-a_i|^\alpha$ near $u=a_i$, $0<\alpha<2$. Such solutions are, over large regions, identically equal to some zeroes of the potential $a_i$'s. We establish the estimates \begin{equation*} \mathcal{L}^n(I_0\cap B_r(x_0))\leq c_1r^{n-1},\quad \mathcal{H}^{n-1}(\partial^* I_0\cap B_r(x_0))\geq c_2r^{n-1}, \quad r\geq r_0(x_0) \end{equation*} for the diffuse interface $I_0:=\{x\in\mathbb{R}^n: \min_{1\leq i\leq N}|u(x)-a_i|>0\}$ and the free boundary $\partial I_0$. Furthermore, if $\alpha=1$ we establish the upper bound \begin{equation*} \mathcal{H}^{n-1}(\partial^* I_0\cap B_r(x_0))\leq c_3r^{n-1}, \quad r\geq r_0(x_0). \end{equation*}, Comment: 25 pages

We study the effect of boundary rugosity in nematic liquid crystalline systems. We consider a highly general formulation of the problem, able to simultaneously deal with several liquid crystal theories. We use techniques of Gamma convergence and demonstrate that the effect of fine-scale surface oscillations may be replaced by an effective homogenised surface energy on a simpler domain. The homogenisation limit is then quantitatively studied in a simplified setting, obtaining convergence rates.

We study entire minimizers of the Allen-Cahn systems. The specific feature of our systems are potentials having a finite number of global minima, with sub-quadratic behaviour locally near their minima. The corresponding formal Euler-Lagrange equations are supplemented with free boundaries. We do not study regularity issues but focus on qualitative aspects. We show the existence of entire solutions in an equivariant setting connecting the minima of $ W $ at infinity, thus modeling many coexisting phases, possessing free boundaries and minimizing energy in the symmetry class. We also present a very modest result of existence of free boundaries under no symmetry hypotheses. The existence of a free boundary can be related to the existence of a specific sub-quadratic feature, a dead core, whose size is also quantified.

We consider a Landau-de Gennes model for a polydisperse, inhomogeneous suspension of colloidal inclusions in a nematic host, in the dilute regime. We study the homogenised limit and compute the effective free energy of the composite material. By suitably choosing the shape of the inclusions and imposing a quadratic, Rapini-Papoular type surface anchoring energy density, we obtain an effective free energy functional with an additional linear term, which may be interpreted as an "effective field" induced by the inclusions. Moreover, we compute the effective free energy in a regime of "very strong anchoring", that is, when the surface energy effects dominate over the volume free energy., Comment: 24 pages, 1 figure

In modelling of chemical, physical or biological systems it may occur that the coefficients, multiplying various terms in the equation of interest, differ greatly in magnitude, if a particular system of units is used. Such is, for instance, the case of the Population Balance Equations (PBE) proposed to model the Latex Particles Morphology formation. The obvious way out of this difficulty is the use of dimensionless scaled quantities, although often the scaling procedure is not unique. In this paper, we introduce a conceptually new general approach, called Optimal Scaling (OS). The method is tested on the known examples from classical and quantum mechanics, and applied to the Latex Particles Morphology model, where it allows us to reduce the variation of the relevant coefficients from 49 to just 4 orders of magnitudes. The PBE are then solved by a novel Generalised Method Of Characteristics, and the OS is shown to help reduce numerical error, and avoid unphysical behaviour of the solution. Although inspired by a particular application, the proposed scaling algorithm is expected find application in a wide range of chemical, physical and biological problems., Comment: 29 pages, 7 figures, 5 tables, 34 references, 3 appendices. Other author's papers can be downloaded at http://www.denys-dutykh.com/

We consider a variational two-dimensional Landau-de Gennes model in the theory of nematic liquid crystals in a disk of radius $R$. We prove that under a symmetric boundary condition carrying a topological defect of degree $\frac{k}{2}$ for some given {\bf even} non-zero integer $k$, there are exactly two minimizers for all large enough $R$. We show that the minimizers do not inherit the full symmetry structure of the energy functional and the boundary data. We further show that there are at least five symmetric critical points.

We study the asymptotic behavior of the minimisers of the Landau-de Gennes model for nematic liquid crystals in a two-dimensional domain in the regime of small elastic constant. At leading order in the elasticity constant, the minimum-energy configurations can be described by the simpler Oseen-Frank theory. Using a refined notion of $\Gamma$-development we recover Landau-de Gennes corrections to the Oseen-Frank energy. We provide an explicit characterisation of minimizing $Q$-tensors at this order in terms of optimal Oseen-Frank directors and observe the emerging biaxiality. We apply our results to distinguish between optimal configurations in the class of conformal director fields of fixed topological degree saturating the lower bound for the Oseen-Frank energy.

The main aim of this note is to prove a sharp Poincar\'e-type inequality for vector-valued functions on $\mathbb{S}^2$, that naturally emerges in the context of micromagnetics of spherical thin films., Comment: 2 figures

We consider a Landau-de Gennes model for a suspension of small colloidal inclusions in a nematic host. We impose suitable anchoring conditions at the boundary of the inclusions, and we work in the dilute regime - i.e., the size of the inclusions is much smaller than the typical separation distance between them, so that the total volume occupied by the inclusions is small. By studying the homogenised limit, and proving rigorous convergence results for local minimisers, we compute the effective free energy for the doped material. In particular, we show that not only the phase transition temperature, but any coefficient of the quartic Landau-de Gennes bulk potential can be tuned, by suitably choosing the surface anchoring energy density., Comment: 31 pages

For $\epsilon>0$, we consider the Ginzburg-Landau functional for $\mathbb R^N$-valued maps defined in the unit ball $B^N\subset \mathbb R^N$ with the vortex boundary data $x$ on $\partial B^N$. In dimensions $N\geq 7$, we prove that for every $\epsilon>0$, there exists a unique global minimizer $u_\epsilon$ of this problem; moreover, $u_\epsilon$ is symmetric and of the form $u_\epsilon(x)=f_\epsilon(|x|)\frac{x}{|x|}$ for $x\in B^N$.

This paper considers the network flow stabilization problem in power systems and adopts an output regulation viewpoint. Building upon the structure of a heterogeneous port-Hamiltonian model, we integrate network aspects and develop a systematic control design procedure. First, the passive output is selected to encode two objectives: consensus in angular velocity and constant excitation current. Second, the non-Euclidean nature of the angle variable reveals the geometry of a suitable target set, which is compact and attractive for the zero dynamics. On this set, circuit-theoretic aspects come into play, giving rise to a network potential function which relates the electrical circuit variables to the machine rotor angles. As it turns out, this energy function is convex in the edge variables, concave in the node variables and, most importantly, can be optimized via an intrinsic gradient flow, with its global minimum corresponding to angle synchronization. The third step consists of explicitly deriving the steady-state-inducing control action by further refining this sequence of control-invariant sets. Analogously to solving the so called regulator equations, we obtain an impedance-based network flow map leading to novel error coordinates and a shifted energy function. The final step amounts to decoupling the rotor current dynamics via feedback-linearziation resulting in a cascade which is used to construct an energy-based controller hierarchically., Comment: In preparation for MTNS 2018

In this paper we study a limit connecting a scaled wave map with the heat flow into the unit sphere $\mathbb{S}^2$. We show quantitatively how that the two equations are connected by means of an initial layer correction. This limit is motivated as a first step into understanding the limit of zero inertia for the hyperbolic-parabolic Ericksen-Leslie's liquid crystal model., Comment: 21 pages