Velocity-based moving mesh methods for nonlinear partial differential equations

This article describes a number of velocity-based moving mesh numerical methods formultidimensional nonlinear time-dependent partial differential equations (PDEs). It consists of a short historical review followed by a detailed description of a recently developed multidimensional moving mesh finite element method based on conservation. Finite element algorithms are derived for both mass-conserving and non mass-conserving problems, and results shown for a number of multidimensional nonlinear test problems, including the second order porous medium equation and the fourth order thin film equation as well as a two-phase problem. Further applications and extensions are referenced.

Coalescence of particles by differential sedimentation

We consider a three dimensional system consisting of a large number of small spherical particles, distributed in a range of sizes and heights (with uniform distribution in the horizontal direction). Particles move vertically at a size-dependent terminal velocity. They are either allowed to merge whenever they cross or there is a size ratio criterion enforced to account for collision efficiency. Such a system may be described, in mean field approximation, by the Smoluchowski kinetic equation with a differential sedimentation kernel. We obtain self-similar steady-state and time-dependent solutions to the kinetic equation, using methods borrowed from weak turbulence theory. Analytical results are compared with direct numerical simulations (DNS) of moving and merging particles, and a good agreement is found.

A non-autonomous strongly damped wave equation: Existence and continuity of the pullback attractor

In this paper we consider the strongly damped wave equation with time-dependent terms u(tt) - Delta u - gamma(t)Delta u(t) + beta(epsilon)(t)u(t) = f(u), in a bounded domain Omega subset of R(n), under some restrictions on beta(epsilon)(t), gamma(t) and growth restrictions on the nonlinear term f. The function beta(epsilon)(t) depends on a parameter epsilon, beta(epsilon)(t) -> 0. We will prove, under suitable assumptions, local and global well-posedness (using the uniform sectorial operators theory), the existence and regularity of pullback attractors {A(epsilon)(t) : t is an element of R}, uniform bounds for these pullback attractors, characterization of these pullback attractors and their upper and lower semicontinuity at epsilon = 0. (C) 2010 Elsevier Ltd. All rights reserved.

Numerical prediction of three-dimensional time-dependent viscoelastic extrudate swell using differential and algebraic models

This study investigates the numerical simulation of three-dimensional time-dependent viscoelastic free surface flows using the Upper-Convected Maxwell (UCM) constitutive equation and an algebraic explicit model. This investigation was carried out to develop a simplified approach that can be applied to the extrudate swell problem. The relevant physics of this flow phenomenon is discussed in the paper and an algebraic model to predict the extrudate swell problem is presented. It is based on an explicit algebraic representation of the non-Newtonian extra-stress through a kinematic tensor formed with the scaled dyadic product of the velocity field. The elasticity of the fluid is governed by a single transport equation for a scalar quantity which has dimension of strain rate. Mass and momentum conservations, and the constitutive equation (UCM and algebraic model) were solved by a three-dimensional time-dependent finite difference method. The free surface of the fluid was modeled using a marker-and-cell approach. The algebraic model was validated by comparing the numerical predictions with analytic solutions for pipe flow. In comparison with the classical UCM model, one advantage of this approach is that computational workload is substantially reduced: the UCM model employs six differential equations while the algebraic model uses only one. The results showed stable flows with very large extrudate growths beyond those usually obtained with standard differential viscoelastic models. (C) 2010 Elsevier Ltd. All rights reserved.

On the dominance of roots of characteristic equations for neutral functional differential equations

We present a sufficient condition for a zero of a function that arises typically as the characteristic equation of a linear functional differential equations of neutral type, to be simple and dominant. This knowledge is useful in order to derive the asymptotic behaviour of solutions of such equations. A simple characteristic equation, arisen from the study of delay equations with small delay, is analyzed in greater detail. (C) 2009 Elsevier Inc. All rights reserved.

Invariants of binary differential equations

In this paper, we study binary differential equations a(x, y)dy (2) + 2b(x, y) dx dy + c(x, y)dx (2) = 0, where a, b, and c are real analytic functions. Following the geometric approach of Bruce and Tari in their work on multiplicity of implicit differential equations, we introduce a definition of the index for this class of equations that coincides with the classical Hopf`s definition for positive binary differential equations. Our results also apply to implicit differential equations F(x, y, p) = 0, where F is an analytic function, p = dy/dx, F (p) = 0, and F (pp) not equal aEuro parts per thousand 0 at the singular point. For these equations, we relate the index of the equation at the singular point with the index of the gradient of F and index of the 1-form omega = dy -aEuro parts per thousand pdx defined on the singular surface F = 0.

A DISCRETIZATION SCHEME FOR AN ONE-DIMENSIONAL REACTION-DIFFUSION EQUATION WITH DELAY AND ITS DYNAMICS

The goal of this paper is to present an approximation scheme for a reaction-diffusion equation with finite delay, which has been used as a model to study the evolution of a population with density distribution u, in such a way that the resulting finite dimensional ordinary differential system contains the same asymptotic dynamics as the reaction-diffusion equation.

On the Hamilton-Jacobi equation in the framework of generalized functions

In this work we study, in the framework of Colombeau`s generalized functions, the Hamilton-Jacobi equation with a given initial condition. We have obtained theorems on existence of solutions and in some cases uniqueness. Our technique is adapted from the classical method of characteristics with a wide use of generalized functions. We were led also to obtain some general results on invertibility and also on ordinary differential equations of such generalized functions. (C) 2011 Elsevier Inc. All rights reserved.

Evaluation of the light polymerization efficiency of copolymers used in dental formulations by differential scanning calorimetry

Different compositions of visible-light-curable triethylene glycol dimethacrylate/bisglycidyl methacrylate copolymers used in dental resin formulations were prepared through copolymerization photoinitiated by a camphorquinone/ethyl 4-dimethylaminobenzoate system irradiated with an Ultrablue IS light-emitting diode. The obtained copolymers were evaluated with differential scanning calorimetry. From the data for the heat of polymerization, before and after light exposure, obtained from exothermic differential scanning calorimetry curves, the light polymerization efficiency or degree of conversion of double bonds was calculated. The glass-transition temperature also was determined before and after photopolymerization. After the photopolymerization, the glass-transi-tion temperature was not well defined because of the breadth of the transition region associated with the properties of the photocured dimethacrylate. The glass-transition temperature after photopolymerization was determined experimentally and compared with the values determined with the Fox equation. In all mixtures, the experimental value was lower than the calculated value. Scanning electron microscopy was used to analyze the morphological differences in the prepared copolymer structures. (C) 2007 Wiley Periodicals, Inc.

Dynamic hedging with stocastic differential utility

In this paper we study the dynamic hedging problem using three different utility specifications: stochastic differential utility, terminal wealth utility, and we propose a particular utility transformation connecting both previous approaches. In all cases, we assume Markovian prices. Stochastic differential utility, SDU, impacts the pure hedging demand ambiguously, but decreases the pure speculative demand, because risk aversion increases. We also show that consumption decision is, in some sense, independent of hedging decision. With terminal wealth utility, we derive a general and compact hedging formula, which nests as special all cases studied in Duffie and Jackson (1990). We then show how to obtain their formulas. With the third approach we find a compact formula for hedging, which makes the second-type utility framework a particular case, and show that the pure hedging demand is not impacted by this specification. In addition, with CRRA- and CARA-type utilities, the risk aversion increases and, consequently the pure speculative demand decreases. If futures price are martingales, then the transformation plays no role in determining the hedging allocation. We also derive the relevant Bellman equation for each case, using semigroup techniques.

Computer simulation of the stochastic transport equation

Trabalho apresentado no 37th Conference on Stochastic Processes and their Applications - July 28 - August 01, 2014 -Universidad de Buenos Aires

A Higher order and stable method for the numerical integration of Random Differential Equations

Trabalho apresentado no Congresso Nacional de Matemática Aplicada à Indústria, 18 a 21 de novembro de 2014, Caldas Novas - Goiás

Stability of multistep second derivative methods for integro-differential equations

The stability of multistep second derivative methods for integro-differential equations is examined through a test equation which allows for the construction of the associated characteristic polynomial and its region of stability (roots in the unit circle) at a proper parameter space. (c) 2004 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Exact closed-form solutions of the Dirac equation for a class of effective Morse potentials

Exact bounded solutions for a fermion subject to exponential scalar potential in 1 + 1 dimensions are found in closed form. We discuss the existence of zero modes which are related to the ultrarelativistic limit of the Dirac equation and are responsible for the induction of a fractional fermion number on the vacuum.

Low-energy direct muon transfer from H to Ne10+, S16+ and Ar18+ using the two-state close-coupling approximation to the Faddeev-Hahn-type equation

We perform a three-body calculation of direct muon-transfer rates from thermalized muonic hydrogen isotopes to bare nuclei Ne10+, S16+ and Ar18+ employing integro-differential Faddeev-Hahn-type equations in configuration space with a two-state close-coupling approximation scheme. All Coulomb potentials including the strong final-state Coulomb repulsion are treated exactly. A long-range polarization potential is included in the elastic channel to take into account the high polarizability of the muonic hydrogen. The transfer rates so-calculated are in good agreement with recent experiments. We find that the muon is captured predominantly in the n = 6, 9 and 10 states of muonic Ne10+, S16+ and Ar18+, respectively.