Convergent discrete numerical solutions of strongly coupled mixed parabolic systems

This work has been partially supported by the Spanish D.G.I.C.Y.T. grant BMF 2000-0206-C04-04Jódar Sánchez, LA.; Casabán, M. (2003). Convergent discrete numerical solutions of strongly coupled mixed parabolic systems. UTILITAS MATHEMATICA. 63:151-172. http://hdl.handle.net/10251/161860S1511726

Exact Solutions and Continuous Numerical Approximations of Coupled Systems of Diffusion Equations with Delay

In this work, we obtain exact solutions and continuous numerical approximations for mixed problems of coupled systems of diffusion equations with delay. Using the method of separation of variables, and based on an explicit expression for the solution of the separated vector initial-value delay problem, we obtain exact infinite series solutions that can be truncated to provide analytical–numerical solutions with prescribed accuracy in bounded domains. Although usually implicit in particular applications, the method of separation of variables is deeply correlated with symmetry ideas.This research was funded by Ministerio de Economía y Competitividad grant number CGL2017-89804-R

Mixed Hyperbolic - Second-Order Parabolic Formulations of General Relativity

Two new formulations of general relativity are introduced. The first one is a parabolization of the Arnowitt, Deser, Misner (ADM) formulation and is derived by addition of combinations of the constraints and their derivatives to the right-hand-side of the ADM evolution equations. The desirable property of this modification is that it turns the surface of constraints into a local attractor because the constraint propagation equations become second-order parabolic independently of the gauge conditions employed. This system may be classified as mixed hyperbolic - second-order parabolic. The second formulation is a parabolization of the Kidder, Scheel, Teukolsky formulation and is a manifestly mixed strongly hyperbolic - second-order parabolic set of equations, bearing thus resemblance to the compressible Navier-Stokes equations. As a first test, a stability analysis of flat space is carried out and it is shown that the first modification exponentially damps and smoothes all constraint violating modes. These systems provide a new basis for constructing schemes for long-term and stable numerical integration of the Einstein field equations.Comment: 19 pages, two column, references added, two proofs of well-posedness added, content changed to agree with submitted version to PR

Numerical approximation of corotational dumbbell models for dilute polymers

We construct a general family of Galerkin methods for the numerical approximation of weak solutions to a coupled microscopic-macroscopic bead-spring model that arises from the kinetic theory of dilute solutions of polymeric liquids with noninteracting polymer chains. The model consists of the unsteady incompressible Navier-Stokes equations in a bounded domain Ω in R d, d=2 or 3, for the velocity and the pressure of the fluid, with an elastic extra-stress tensor as right-hand side in the momentum equation. The extra-stress tensor stems from the random movement of the polymer chains and is defined through the associated probability density function which satisfies a Fokker-Planck type parabolic equation, a crucial feature of which is the presence of a centre-of-mass diffusion term. We focus on finitely-extensible nonlinear elastic, FENE-type, dumbbell models. In the case of a corotational drag term we perform a rigorous passage to the limit as the spatial and temporal discretization parameters tend to zero, and show that a (sub)sequence of numerical solutions converges to a weak solution of this coupled Navier-Stokes-Fokker-Planck system

The LifeV library: engineering mathematics beyond the proof of concept

LifeV is a library for the finite element (FE) solution of partial differential equations in one, two, and three dimensions. It is written in C++ and designed to run on diverse parallel architectures, including cloud and high performance computing facilities. In spite of its academic research nature, meaning a library for the development and testing of new methods, one distinguishing feature of LifeV is its use on real world problems and it is intended to provide a tool for many engineering applications. It has been actually used in computational hemodynamics, including cardiac mechanics and fluid-structure interaction problems, in porous media, ice sheets dynamics for both forward and inverse problems. In this paper we give a short overview of the features of LifeV and its coding paradigms on simple problems. The main focus is on the parallel environment which is mainly driven by domain decomposition methods and based on external libraries such as MPI, the Trilinos project, HDF5 and ParMetis. Dedicated to the memory of Fausto Saleri.Comment: Review of the LifeV Finite Element librar

Fast derivatives of likelihood functionals for ODE based models using adjoint-state method

We consider time series data modeled by ordinary differential equations (ODEs), widespread models in physics, chemistry, biology and science in general. The sensitivity analysis of such dynamical systems usually requires calculation of various derivatives with respect to the model parameters. We employ the adjoint state method (ASM) for efficient computation of the first and the second derivatives of likelihood functionals constrained by ODEs with respect to the parameters of the underlying ODE model. Essentially, the gradient can be computed with a cost (measured by model evaluations) that is independent of the number of the ODE model parameters and the Hessian with a linear cost in the number of the parameters instead of the quadratic one. The sensitivity analysis becomes feasible even if the parametric space is high-dimensional. The main contributions are derivation and rigorous analysis of the ASM in the statistical context, when the discrete data are coupled with the continuous ODE model. Further, we present a highly optimized implementation of the results and its benchmarks on a number of problems. The results are directly applicable in (e.g.) maximum-likelihood estimation or Bayesian sampling of ODE based statistical models, allowing for faster, more stable estimation of parameters of the underlying ODE model.Comment: 5 figure

A thermodynamically consistent model of magneto-elastic materials under diffusion at large strains and its analysis

The theory of elastic magnets is formulated under possible diffusion and heat flow governed by Fick's and Fourier's laws in the deformed (Eulerian) configuration, respectively. The concepts of nonlocal nonsimple materials and viscous Cahn-Hilliard equations are used. The formulation of the problem uses Lagrangian (reference) configuration while the transport processes are pulled back. Except the static problem, the demagnetizing energy is ignored and only local non-selfpenetration is considered. The analysis as far as existence of weak solutions of the (thermo)dynamical problem is performed by a careful regularization and approximation by a Galerkin method, suggesting also a numerical strategy. Either ignoring or combining particular aspects, the model has numerous applications as ferro-to-paramagnetic transformation in elastic ferromagnets, diffusion of solvents in polymers possibly accompanied by magnetic effects (magnetic gels), or metal-hydride phase transformation in some intermetalics under diffusion of hydrogen accompanied possibly by magnetic effects (and in particular ferro-to-antiferromagnetic phase transformation), all in the full thermodynamical context under large strains

Resonance and frequency-locking phenomena in spatially extended phytoplankton-zooplankton system with additive noise and periodic forces

In this paper, we present a spatial version of phytoplankton-zooplankton model that includes some important factors such as external periodic forces, noise, and diffusion processes. The spatially extended phytoplankton-zooplankton system is from the original study by Scheffer [M Scheffer, Fish and nutrients interplay determines algal biomass: a minimal model, Oikos \textbf{62} (1991) 271-282]. Our results show that the spatially extended system exhibit a resonant patterns and frequency-locking phenomena. The system also shows that the noise and the external periodic forces play a constructive role in the Scheffer's model: first, the noise can enhance the oscillation of phytoplankton species' density and format a large clusters in the space when the noise intensity is within certain interval. Second, the external periodic forces can induce 4:1 and 1:1 frequency-locking and spatially homogeneous oscillation phenomena to appear. Finally, the resonant patterns are observed in the system when the spatial noises and external periodic forces are both turned on. Moreover, we found that the 4:1 frequency-locking transform into 1:1 frequency-locking when the noise intensity increased. In addition to elucidating our results outside the domain of Turing instability, we provide further analysis of Turing linear stability with the help of the numerical calculation by using the Maple software. Significantly, oscillations are enhanced in the system when the noise term presents. These results indicate that the oceanic plankton bloom may partly due to interplay between the stochastic factors and external forces instead of deterministic factors. These results also may help us to understand the effects arising from undeniable subject to random fluctuations in oceanic plankton bloom.Comment: Some typos errors are proof, and some strong relate references are adde