In questo articolo si investigano le proprietà di stabilità asintotica dei metodi numerici per equazioni differenziali con ritardo, prendendo in esame l'equazione test: $$U^{\prime }\left(t\right)=aU\left(t\right)+bU\left(t-\tau \right)$$ dove $a$, $b\in \mathbb{R}$, $\tau >0$ e $g\left(t\right)$ è una funzione a valori reali e continua. In particolare, viene analizzata la dipendenza dal ritardo della stabilità numerica dei metodi di collocazione Gaussiana. Nel recente lavoro [GH99], la stabilità di questi metodi è stata dimostrata facendo uso di un approccio geometrico, basato sul legame tra la proprietà...

We present two defect correction schemes to accelerate the Petrov-Galerkin finite element methods [19] for nonlinear Volterra integro-differential equations. Using asymptotic expansions of the errors, we show that the defect correction schemes can yield higher order approximations to either the exact solution or its derivative. One of these schemes even does not impose any extra regularity requirement on the exact solution. As by-products, all of these higher order numerical methods can also be...

In this paper, a new numerical method for solving the nonlinear constrained optimal control with quadratic performance index is presented. The method is based upon B-spline functions. The properties of B-spline functions are presented. The operational matrix of derivative (${𝐃}_{\phi }$) and integration matrix ($𝐏$) are introduced. These matrices are utilized to reduce the solution of nonlinear constrained quadratic optimal control to the solution of nonlinear programming one to which existing well-developed...

A class of linear elliptic operators has an important qualitative property, the so-called maximum principle. In this paper we investigate how this property can be preserved on the discrete level when an interior penalty discontinuous Galerkin method is applied for the discretization of a 1D elliptic operator. We give mesh conditions for the symmetric and for the incomplete method that establish some connection between the mesh size and the penalty parameter. We then investigate the sharpness of...

This paper discusses analytical and numerical issues related to elliptic equations with random coefficients which are generally nonlinear functions of white noise. Singularity issues are avoided by using the Itô-Skorohod calculus to interpret the interactions between the coefficients and the solution. The solution is constructed by means of the Wiener Chaos (Cameron-Martin) expansions. The existence and uniqueness of the solutions are established under rather weak assumptions, the main of which...

We present a numerical algorithm to solve the micromagnetic equations based on tangential-plane minimization for the magnetization update and a homothethic-layer decomposition of outer space for the computation of the demagnetization field. As a first application, detailed results on the flower-vortex transition in the cube of Micromagnetic Standard Problem number 3 are obtained, which confirm, with a different method, those already present in the literature, and validate our method and code. We...

We present a numerical algorithm to solve the micromagnetic equations based on tangential-plane minimization for the magnetization update and a homothethic-layer decomposition of outer space for the computation of the demagnetization field. As a first application, detailed results on the flower-vortex transition in the cube of Micromagnetic Standard Problem number 3 are obtained, which confirm, with a different method, those already present in the literature, and validate our method and...

We prove the existence of solutions to nonlinear parabolic problems of the following type: $$\left\{\begin{array}{cc}{\displaystyle \frac{\partial b\left(u\right)}{\partial t}}+A\left(u\right)=f+\mathrm{div}\left(\Theta (x;t;u)\right)\hfill & \text{in}Q,\hfill \\ u(x;t)=0\hfill & \text{on}\partial \Omega \times [0;T],\hfill \\ b\left(u\right)(t=0)=b\left(u_0\right)\hfill & \text{on}\Omega ,\hfill \end{array}\right.$$ where $b:ℝ\to ℝ$ is a strictly increasing function of class ${𝒞}^1$, the term $$A\left(u\right)=-\mathrm{div}\left(a\right(x,t,u,\nabla u\left)\right)$$ is an operator of Leray-Lions type which satisfies the classical Leray-Lions assumptions of Musielak type, $\Theta :\Omega \times [0;T]\times ℝ\to ℝ$ is a Carathéodory, noncoercive function which satisfies the following condition: ${sup}_{\left|s\right|\le k}\left|\Theta (·,·,s)\right|\in E_{\psi }\left(Q\right)$ for all $k>0$, where $\psi$ is the Musielak complementary function of $\Theta$, and the second term $f$ belongs to $L^1\left(Q\right)$.