A mathematical model for phase separation: A generalized Cahn–Hilliard equation

In this paper we present a mathematical model to describe the phenomenon of phase separation, which is modelled as space regions where an order parameter changes smoothly. The model proposed, including thermal and mixing effects, is deduced for an incompressible fluid, so the resulting differential system couples a generalized Cahn–Hilliard equation with the Navier–Stokes equation. Its consistency with the second law of thermodynamics in the classical Clausius–Duhem form is finally proved. Copyright © 2011 John Wiley & Sons, Ltd.     

A New Unified Approach for the Simulation of a Wide Class of Directional Distributions

The need for effective simulation methods for directional distributions has grown as they have become components in more sophisticated statistical models. A new acceptance–rejection method is proposed and investigated for the Bingham distribution on the sphere using the angular central Gaussian distribution as an envelope. It is shown that the proposed method has high efficiency and is also straightforward to use. Next, the simulation method is extended to the Fisher and Fisher–Bingham distributions on spheres and related manifolds. Together, these results provide a widely applicable and efficient methodology to simulate many of the standard models in directional data analysis. An R package simdd, available in the online supplementary material, implements these simulation methods.     

A note on the Cauchy problem of the generalized Davey–Stewartson equations

In this note, we establish some local and global existence results for the Cauchy problem of a class of nonlinear dispersive equations which generalize the nonlinear Schrödinger equations and the Davey–Stewartson equations. These results improve some previously obtained results by some other authors when they are restricted to certain special equations.     

A sparse grid method for the Navier–Stokes equations based on hyperbolic cross

A sparse grid method for the time‐dependent Navier–Stokes equations based on hyperbolic cross approximation is considered in this article. Subsequent truncation of the associated series expansion results in a sparse grid discretization. Stability and convergence of the fully discrete sparse grid method are established. Finally, the numerical experiment is presented to show the effectiveness of this sparse grid method. Copyright © 2013 John Wiley & Sons, Ltd.     

A multicriteria model for risk sorting of natural gas pipelines based on ELECTRE TRI integrating Utility Theory

This paper proposes a multicriteria model for assessing risk in natural gas pipelines, and for classifying sections of pipeline into risk categories. The model integrates Utility Theory and the ELECTRE TRI method. It aims to help transmission and distribution companies, when engaged in risk management and decision-making, to consider the multiple dimensions of risk that may arise from pipeline accidents. Pipeline hazard scenarios are presented, and it is argued that the assessment of risk in natural gas pipelines should not be based solely on probabilities of human fatalities, but should involve a wider perspective that simultaneously takes into consideration the human, environmental and financial dimensions of impacts of pipeline accidents. Finally, in order to verify the effectiveness of the model set out, a numerical application based on a real case study is presented.     

A continuum theory for first‐order phase transitions based on the balance of structure order

First‐order phase transitions are modelled by a non‐homogeneous, time‐dependent scalar‐valued order parameter or phase field. The time dependence of the order parameter is viewed as arising from a balance law of the structure order. The gross motion is disregarded and hence the body is regarded merely as a heat conductor. Compatibility of the constitutive functions with thermodynamics is exploited by expressing the second law through the classical Clausius–Duhem inequality. First, a model for conductors without memory is set up and the order parameter is shown to satisfy a maximum theorem. Next, heat conductors with memory are considered. Different evolution problems are established through a system of differential equations whose form is related to the manner in which the memory property is represented. Copyright © 2007 John Wiley & Sons, Ltd.     

A monotone iterative scheme for a nonlinear second order equation based on a generalized anti–maximum principle

In this paper, a generalized anti–maximum principle for the second order differential operator with potentials is proved. As an application, we will give a monotone iterative scheme for periodic solutions of nonlinear second order equations. Such a scheme involves the L^p norms of the growth, 1 ≤ p ≤ ∞, while the usual one is just the case p = ∞.     

A series of high‐order quasi‐compact schemes for space fractional diffusion equations based on the superconvergent approximations for fractional derivatives

Based on the superconvergent approximation at some point (depending on the fractional order α, but not belonging to the mesh points) for Grünwald discretization to fractional derivative, we develop a series of high‐order quasi‐compact schemes for space fractional diffusion equations. Because of the quasi‐compactness of the derived schemes, no points beyond the domain are used for all the high‐order schemes including second‐order, third‐order, fourth‐order, and even higher‐order schemes; moreover, the algebraic equations for all the high‐order schemes have the completely same matrix structure. The stability and convergence analysis for some typical schemes are made; the techniques of treating the fractional derivatives with nonhomogeneous boundaries are introduced; and extensive numerical experiments are performed to confirm the theoretical analysis or verify the convergence orders. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1345–1381, 2015     

A numerically efficient closed-form representation of mean-variance hedging for exponential additive processes based on Malliavin calculus

ABSTRACT

We focus on mean-variance hedging problem for models whose asset price follows an exponential additive process. Some representations of mean-variance hedging strategies for jump-type models have already been suggested, but none is suited to develop numerical methods of the values of strategies for any given time up to the maturity. In this paper, we aim to derive a new explicit closed-form representation, which enables us to develop an efficient numerical method using the fast Fourier transforms. Note that our representation is described in terms of Malliavin derivatives. In addition, we illustrate numerical results for exponential Lévy models.     

A collocation method based on Bernoulli operational matrix for numerical solution of generalized pantograph equation

This paper presents a direct solution technique for solving the generalized pantograph equation with variable coefficients subject to initial conditions, using a collocation method based on Bernoulli operational matrix of derivatives. Only small dimension of Bernoulli operational matrix is needed to obtain a satisfactory result. Numerical results with comparisons are given to confirm the reliability of the proposed method for generalized pantograph equations.     

A finite element variational multiscale method for steady‐state natural convection problem based on two local gauss integrations

In this article, supposing that the velocity, pressure, and temperature are approximated by the elements , and applying the orthogonal projection technique, we introduce two Gauss integrations as a stabilizing term in the common variational multiscale (VMS) method and derive a new VMS (Two Gauss VMS) method for steady‐state natural convection problem. Comparing with the common VMS method, the Two Gauss VMS method does not need to introduce any extra variable and reduces the degrees of freedom of the discrete system a lot, but gets the same stabilized result. The effectiveness and stability of the Two Gauss VMS method are further demonstrated through two numerical examples. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 361–375, 2014     

A gray model for increasing sequences with nonhomogeneous index trends based on fractional‐order accumulation

By comparing the class ratio deviation and restoring error of first‐order accumulation with that of fractional‐order accumulation, a gray model for monotonically increasing sequences can obtain optimal simulation accuracy via selecting a proper cumulative order. In this study, a gray model for increasing sequences with nonhomogeneous index trends based on fractional‐order accumulation is proposed. To reduce the modeling error caused by the background value and to improve the prediction accuracy of the model, an optimized model using the 3/8 Simpson formula is constructed. Finally, the 2 proposed models are used to predict the total energy consumption in China and the monthly sales of new products in an enterprise. Compared with the GM(1,1) model based on fractional‐order accumulation, the proposed model exhibits better simulation and prediction accuracy.     

A numerical approach based on ln‐shifted Legendre polynomials for solving a fractional model of pollution

The model of pollution for a system of 3 lakes interconnected by channels is extended using Caputo‐Hadamard fractional derivatives of different orders α_i∈(0,1), i=1,2,3. A numerical approach based on ln‐shifted Legendre polynomials is proposed to solve the considered fractional model. No discretization is needed in our approach. Some numerical experiments are provided to illustrate the presented method.     

Modeling autocorrelation functions of long-range dependent teletraffic series based on optimal approximation in Hilbert space—A further study

This paper discusses optimal approximation of autocorrelation functions of teletraffic series by introducing a generalization of autocorrelation function form of fractional Gaussian noise (FGN). The demonstrations with real-traffic series are given.     

A reduced‐order extrapolated natural boundary element method based on POD for the 2D hyperbolic equation in unbounded domain

In this article, we primarily focuses to study the order‐reduction for the classical natural boundary element (NBE) method for the two‐dimensional (2D) hyperbolic equation in unbounded domain. To this end, we first build a semi‐discretized format about time for the hyperbolic equation and discuss the existence, stability, and convergence of the time semi‐discretized solutions. We then establish the classical fully discretized NBE format from the time semi‐discretized one and analyze the existence, stability, and convergence of the classical NBE solutions. Next, using proper orthogonal decomposition method, we build a reduced‐order extrapolated NBE (ROENBE) format containing very few unknowns but having adequately high accuracy, and we also discuss the existence, stability, and convergence of the ROENBE solutions. Finally, we use some numerical examples to show that the ROENBE method is far superior to the classical NBE one. It shows that the ROENBE method is reliable and effective for solving the 2D hyperbolic equation with the unbounded domain.     

A positive cell vertex Godunov scheme for a Beeler–Reuter based model of cardiac electrical activity

The monodomain model is a widely used model in electrocardiology to simulate the propagation of electrical potential in the myocardium. In this paper, we investigate a positive nonlinear control volume finite element scheme, based on Godunov's flux approximation of the diffusion term, for the monodomain model coupled to a physiological ionic model (the Beeler–Reuter model) and using an anisotropic diffusion tensor. In this scheme, degrees of freedom are assigned to vertices of a primal triangular mesh, as in conforming finite element methods. The diffusion term which involves an anisotropic tensor is discretized on a dual mesh using the diffusion fluxes provided by the conforming finite element reconstruction on the primal mesh and the other terms are discretized by means of an upwind finite volume method on the dual mesh. The scheme ensures the validity of the discrete maximum principle without any restriction on the transmissibility coefficients. By using a compactness argument, we obtain the convergence of the discrete solution and as a consequence, we get the existence of a weak solution of the original model. Finally, we illustrate by numerical simulations that the proposed scheme successfully removes nonphysical oscillations in the propagation of the wavefront and maintains conduction velocity close to physiological values.     

A series of exact solutions for (2 + 1)-dimensional Wick-type stochastic generalized Broer-Kaup system via a modified variable-coefficient projective Riccati equation mapping method

A modified variable-coefficient projective Riccati equation mapping method is applied to (2 + 1)-dimensional Wick-type stochastic generalized Broer-Kaup system. With the help of Hermit transformation, we obtain a series of new exact stochastic solutions to the stochastic Broer-Kaup system in the white noise environment.