On the stability of a finite-difference scheme for nonlocal parabolic boundary-value problems

We deal with the stability analysis of difference schemes for a one-dimensional parabolic equation subject to integral conditions. It is based on the spectral structure of the transition matrix of a difference scheme. The stability domain is defined by using the hyperbola which is the locus of points where the transition matrix has trivial eigenvalues. The stability conditions obtained are much more general compared with those known in the literature. We analyze three separate cases of nonlocal integral conditions and solve an example illustrating the efficiency of the technique.     

Dynamic stability of a curved rail under a moving load

The stability of a curved rail under a constant moving load has been investigated using a linear theory; critical speeds of the moving load, and the dynamic rail deflections and rotation were calculated. The effect of the foundation was included through distributed linear springs. It was assumed that the moving load remains in constant contact with the rail and travels along a fixed path on the rail head.     

On the use of AGE algorithm with a high accuracy Numerov type variable mesh discretization for 1D non-linear parabolic equations

In this paper, the use of N-AGE and Newton-N-AGE iterative methods on a variable mesh for the solution of one dimensional parabolic initial boundary value problems is considered. Using three spatial grid points, a two level implicit formula based on Numerov type discretization is discussed. The local truncation error of the method is of O(k²h_l^-1 +kh_l +h_l³)O({k^2h_l^{-1} +kh_l +h_l^3}), where h _l  > 0 and k > 0 are the step lengths in space and time directions, respectively. We use a special technique to handle singular parabolic equations. The advantage of using these algorithms is highlighted computationally.     

Numerical stability study and error estimation for two implicit schemes in a moving boundary problem

Two algorithms are described [Ferris D. H. (fixed time‐step method) and Gupta and Kumar (variable time‐step method)] that solve a mathematical model for the study of the one‐dimensional moving boundary problem with implicit boundary conditions. Landau's transformation is used, in order to work with a fixed number of nodes at each time‐step. The p.d.e. is discretized using an implicit finite difference scheme. The mathematical model describes the oxygen diffusion in absorbing tissues. An important application is the estimation of time‐variant radiation treatments of cancerous tumors. © 2000 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16: 42–61, 2000     

Predictive control of uncertain nonlinear parabolic PDE systems using a Galerkin/neural-network-based model

In this paper, a model predictive control (MPC) scheme for a class of parabolic partial differential equation (PDE) systems with unknown nonlinearities, arising in the context of transport-reaction processes, is proposed. A spatial operator of a parabolic PDE system is characterized by a spectrum that can be partitioned into a finite slow and an infinite fast complement. In this view, first, Galerkin method is used to derive a set of finite dimensional slow ordinary differential equation (ODE) system that captures the dominant dynamics of the initial PDE system. Then, a Multilayer Neural Network (MNN) is employed to parameterize the unknown nonlinearities in the resulting finite dimensional ODE model. Finally, a Galerkin/neural-network-based ODE model is used to predict future states in the MPC algorithm. The proposed controller is applied to stabilize an unstable steady-state of the temperature profile of a catalytic rod subject to input and state constraints.     

Convergence analysis of a covolume scheme for Maxwell's equations in three dimensions

This paper contains error estimates for covolume discretizations of Maxwell's equations in three space dimensions. Several estimates are proved. First, an estimate for a semi-discrete scheme is given. Second, the estimate is extended to cover the classical interlaced time marching technique. Third, some of our unstructured mesh results are specialized to rectangular meshes, both uniform and nonuniform. By means of some additional analysis it is shown that the spatial convergence rate is one order higher than for the unstructured case.

     

An efficient boundary element method for a class of parabolic differential equations using discretization in time

A new boundary element method using discretization in time is proposed to solve a class of parabolic differential equations. The method treats the term containing the time derivative as a forcing term. This necessitates the introduction of additional unknowns in the interior of the domain. At the same time, however, values for the dependent variable are determined directly in the interior. The boundary element formulation is reduced to essentially solving a Poisson equation. The accuracy and efficiency of the method are demonstrated with several examples.     

Non-linear parametric vibration and stability analysis for two dynamic models of axially moving Timoshenko beams

Non-linear parametric vibration and stability of an axially moving Timoshenko beam are considered for two dynamic models; the first one, with considering only the transverse displacement and the second one, with considering both longitudinal and transverse displacements. The set of non-linear partial-differential equations of both models are derived using an energy approach. The method of multiple scales is applied directly to both models, and using the equation order one, the mode shape equations and natural frequencies are obtained. Then, for the equation order epsilon, the solvability conditions are considered for the resonance case and the stability boundaries are formulated analytically via Routh–Hurwitz criterion. Eventually, some numerical examples are provided to show the differences in the behavior of the above-mentioned non-linear models.     

Uniform stability for solutions of a structural acoustics PDE model with no added dissipative feedback

A rate of rational decay is obtained for smooth solutions of a PDE model, which has been used in the literature to describe structural acoustic flows. This structural acoustics model is composed of two distinct PDE systems: (i) a wave equation, to model the interior acoustic flow within the given cavity Ω and (ii) a structurally damped elastic equation, to describe time‐evolving displacements along the flexible portion Γ₀ of the cavity walls. Moreover, the extent of damping in this elastic component is quantified by parameter η∈[0,1]. The coupling between these two distinct dynamics occurs across the boundary interface Γ₀. Our main result is the derivation of uniform decay rates for classical solutions of this particular structural acoustic PDE, decay rates that are obtained without incorporating any additional boundary dissipative feedback mechanisms. In particular, in the case that full Kelvin–Voight damping is present in fourth‐order elastic dynamics, that is, the structural acoustics system as it appears in the literature, solutions that correspond to smooth initial data decay at a rate of . By way of deriving these stability results, necessary a priori inequalities for a certain static structural acoustics PDE model are generated here; these inequalities ultimately allow for an application of a recently derived resolvent criterion for rational decay. Copyright © 2016 John Wiley & Sons, Ltd.     

Adaptive mesh refinement for the control of cost and quality in finite element analysis

Adaptive strategies are a necessary tool to make finite element analysis applicable to engineering practice. In this paper, attention is restricted to mesh adaptivity. Traditionally, the most common mesh adaptive strategies for linear problems are used to reach a prescribed accuracy. This goal is best met with an h-adaptive scheme in combination with an error estimator. In an industrial context, the aim of the mechanical simulations in engineering design is not only to obtain greatest quality but more often a compromise between the desired quality and the computation cost (CPU time, storage, software, competence, human cost, computer used). In this paper, we propose the use of alternative mesh refinement criteria with an h-adaptive procedure for 3D elastic problems. The alternative mesh refinement criteria (MR) are based on: prescribed number of elements with maximum accuracy, prescribed CPU time with maximum accuracy and prescribed memory size with maximum accuracy. These adaptive strategies are based on a technique of error in constitutive relation (the process could be used with other error estimators) and an efficient adaptive technique which automatically takes into account the steep gradient areas. This work proposes a 3D method of adaptivity with the latest version of the INRIA automatic mesh generator GAMHIC3D.     

An implicit finite volume method for the solution of 3D low Mach number viscous flows using a local preconditioning technique

This paper presents a cell-centered high order finite volume scheme for the solution of the three-dimensional (3D) Navier–Stokes equations with low Mach number. The system of non-linear equations is solved by means of a fully implicit pseudo-transient scheme. Each pseudo-time step is solved by a Newton-GMRes procedure. A local preconditioning technique is used to scale the speed of sound and to improve the system condition number for low Mach number and low cell Reynolds number. This preconditioning is applied to the AUSM+up flux vector splitting function. The method is tested on 2D and 3D low Mach number laminar flows.     

Mathematical analysis and stability of a chemotaxis model with logistic term

In this paper we study a non‐linear system of differential equations arising in chemotaxis. The system consists of a PDE that describes the evolution of a population and an ODE which models the concentration of a chemical substance. We study the number of steady states under suitable assumptions, the existence of one global solution to the evolution problem in terms of weak solutions and the stability of the steady states. Copyright © 2004 John Wiley & Sons, Ltd.     

Modeling and stability analysis for a cholera model with vaccination

In this paper, we consider a cholera model with vaccination. The disease‐free equilibrium of the system is globally asymptotically stable when the basic reproduction number . If , the disease persists and the unique endemic equilibrium is globally asymptotically stable in the interior of the feasible region under some conditions, which is obtained by compound matrices and geometric approaches. We perform sensitivity analysis of on the parameters in order to determine their relative importance to disease transmission and prevalence. Numerical simulations are presented to illustrate the results. Copyright © 2011 John Wiley & Sons, Ltd.     

An Arbitrary Lagrangian Eulerian formulation for a 3D eutrophication model in a moving domain

This work analyzes a realistic mathematical model that governs eutrophication (an ecological process involving nutrients, phytoplankton, zooplankton, organic detritus and dissolved oxygen) into a moving aquatic domain. As a main result, we obtain existence-uniqueness results for the solution of the system within the general framework of non-cylindrical domains (based on studying the properties of a generic parabolic problem). The fact of dealing with moving domains, and the lack of regularity, preclude the use of standard semigroup approach, forcing us towards the utilization of Arbitrary Lagrangian Eulerian techniques.     

A fast and unconditionally positive finite-difference discretization of a multidimensional equation in nonlinear population dynamics

Departing from a nonlinear diffusion–reaction equation which describes the growth of biological films, we derive a finite-difference discretization which preserves unconditionally the positivity and the boundedness of approximations. The design of this method follows a non-traditional approach in the estimation of first-order partial derivatives, and the technique is a variable step-size and exact methodology for which the properties of existence and uniqueness of non-negative and bounded solutions hold for any initial profile which is likewise non-negative and bounded. As a consequence of the exactness, the computer implementation requires less resources and yields faster results than some linear schemes available in the literature. Qualitative comparisons against known linear and nonlinear techniques show that our method produces similar computer results in the two-dimensional case. Moreover, our scheme is also able to simulate the development of microbial films in the three-dimensional scenario. This is a feature that is not inherent to the linear methodologies considered, in view of the large amount of computational resources that such approaches require.     

An adaptive algorithm in time domain for dynamic analysis of a simply supported beam subjected to a moving vehicle

This paper presents an adaptive algorithm in the time domain for the dynamic analysis of a simply supported beam subjected to the moving load and moving vehicle with/without varying surface roughness. By expanding variables at a discretized time interval, a coupled spatial‐temporal problem can be converted into a series of recursive space problems that are solved by finite element method (FEM), and a piecewised adaptive computing procedure can be carried out for different sizes of time steps. The proposed approach is numerically verified via the comparison with analytical and the Runge–Kutta method‐based solutions, and satisfactory results have been achieved. Copyright © 2011 John Wiley & Sons, Ltd.     

Uniform l1 behaviour for time discretization of a Volterra equation with completely monotonic kernel: I. stability

This paper is the first of two papers on the time discretizationof the equation u_t + ^t ₀ ß (t — s) Au (s) ds= 0, t > 0, u (0) = u₀, where A is a self-adjoint denselydefined linear operator on a Hilbert space H with a completeeigensystem {_m, _m}_{m = 1}, and ß (t) is completely monotonicand locally integrable, but not constant. The equation is discretizedin time using first-order differences in combination with order-oneconvolution quadrature. The stability properties of the timediscretization are derived in the l¹_t (0, ; H) norm.     

Existence, uniqueness and stability analysis of allelopathic stimulatory phytoplankton model

In this paper we consider the two species competitive delay plankton allelopathy stimulatory model system. We show the existence and uniqueness of the solution of the deterministic model. Moreover, we study the persistence of the model and the stability properties of its equilibrium points. We illustrate the theoretical results by some numerical simulations.     

The stability and convergence of two linearized finite difference schemes for the nonlinear epitaxial growth model

The numerical simulation of the dynamics of the molecular beam epitaxy (MBE) growth is considered in this article. The governing equation is a nonlinear evolutionary equation that is of linear fourth order derivative term and nonlinear second order derivative term in space. The main purpose of this work is to construct and analyze two linearized finite difference schemes for solving the MBE model. The linearized backward Euler difference scheme and the linearized Crank‐Nicolson difference scheme are derived. The unique solvability, unconditional stability and convergence are proved. The linearized Euler scheme is convergent with the convergence order of O(τ + h²) and linearized Crank‐Nicolson scheme is convergent with the convergence order of O(τ² + h²) in discrete L²‐norm, respectively. Numerical stability with respect to the initial conditions is also obtained for both schemes. Numerical experiments are carried out to demonstrate the theoretical analysis. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011