A research into the numerical method of Dirichlet's problem of complex Monge-Ampère equation on Cartan-Hartogs domain of the third type

Monge-Ampère equation is a nonlinear equation with high degree, therefore its numerical solution is very important and very difficult. In present paper the numerical method of Dirichlet's problem of Monge-Ampère equation on Cartan-Hartogs domain of the third type is discussed by using the analytic method. Firstly, the Monge-Ampère equation is reduced to the nonlinear ordinary differential equation, then the numerical method of the Dirichlet problem of Monge-Ampère equation becomes the numerical method of two point boundary value problem of the nonlinear ordinary differential equation. Secondly, the solution of the Dirichlet problem is given in explicit formula under the special case, which can be used to check the numerical solution of the Dirichlet problem.     

The construction of numerical integration rules of degree three for product regions

Numerical integration formulas in n-dimensional Euclidean space of degree three are discussed. In this paper, for the product regions a method is presented to construct numerical integration formulas of degree three with 2n real points and positive weights. The presented problem is a little different from those dealt with by other authors. All the corresponding one-dimensional integrals can be different from each other and they are also nonsymmetrical. In this paper an n-dimensional numerical integration problem is turned into n one-dimensional moment problems, which simplifies the construction process. Some explicit numerical formulas are given. Furthermore, a more generalized numerical integration problem is considered, which will shed light on the final solution to the third degree numerical integration problem.     

Numerical solutions of the m-membranes problem

In this paper numerical approximation for the m-membrane problem is considered. We make a change of variables that leads to a different expression of the quadratic functional that allows after discretizing the problem to reformulate it as finite dimensional bound constrained quadratic problem. To our knowledge this is the first paper on numerical approximation of the m-membrane problem. We reformulate the m-membrane problem as a bound constraint quadratic minimization problem. The bound constraint quadratic form is solved with the gradient projection method.     

Numerical method for determining the inhomogeneity boundary in the Dirichlet problem for Laplace’s equation in a piecewise homogeneous medium

The Dirichlet problem for Laplace’s equation in a two-dimensional domain filled with a piecewise homogeneous medium is considered. The boundary of the inhomogeneity is assumed to be unknown. The inverse problem of determining the inhomogeneity boundary from additional information on the solution of the Dirichlet problem is considered. A numerical method based on the linearization of the nonlinear operator equation for the unknown boundary is proposed for solving the inverse problem. The results of numerical experiments are presented.     

The transmutation operator method for efficient solution of the inverse Sturm‐Liouville problem on a half‐line

The inverse Sturm‐Liouville problem on a half‐line is considered. With the aid of a Fourier‐Legendre series representation of the transmutation integral kernel and the Gel'fand‐Levitan equation, the numerical solution of the problem is reduced to a system of linear algebraic equations. The potential q is recovered from the first coefficient of the Fourier‐Legendre series. The resulting numerical method is direct and simple. The results of the numerical experiments are presented.     

Parameter identification in multistage population dynamics model

We present a numerical analysis to solve a parameter identification problem. We identify the demographical parameters of a multistage population dynamics model (Ainseba et al., 2011 [12]). Our nonlinear optimization problem with constraints is solved by a Quasi-Newton method. The convergence proof of this numerical method is performed here. Some numerical applications of it are also given at the end of the paper.     

Isospectral flow method for nonnegative inverse eigenvalue problem with prescribed structure

The nonnegative inverse eigenvalue problem is that given a family of complex numbers λ={λ₁,…,λ_n}, find a nonnegative matrix of order n with spectrum λ. This problem is difficult and remains unsolved partially. In this paper, we focus on its generalization that the reconstructed nonnegative matrices should have some prescribed entries. It is easy to see that this new problem will come back to the common nonnegative inverse eigenvalue problem if there is no constraint of the locations of entries. A numerical isospectral flow method which is developed by hybridizing the optimization theory and steepest descent method is used to study the reconstruction. Moreover, an error estimate of the numerical iteration for ordinary differential equations on the matrix manifold is presented. After that, a numerical method for the nonnegative symmetric inverse eigenvalue problem with prescribed entries and its error estimate are considered. Finally, the approaches are verified by the numerical test results.     

A semidefinite programming approach for solving fractional optimal control problems

This paper presents a numerical scheme for solving fractional optimal control. The fractional derivative in this problem is in the Riemann–Liouville sense. The proposed method, based upon the method of moments, converts the fractional optimal control problem to a semidefinite optimization problem; namely, the nonlinear optimal control problem is converted to a convex optimization problem. The Grunwald–Letnikov formula is also used as an approximation for fractional derivative. The solution of fractional optimal control problem is found by solving the semidefinite optimization problem. Finally, numerical examples are presented to show the performance of the method.     

Implicit algorithms for solving the Cauchy problem for the equations of the dynamics of mechanical systems

It is shown that in the numerical solution of the Cauchy problem for systems of second-order ordinary differential equations, when solved for the highest-order derivative, it is possible to construct simple and economical implicit computational algorithms for step-by-step integration without using laborious iterative procedures based on processes of the Newton-Raphson iterative type. The initial problem must first be transformed to a new argument — the length of its integral curve. Such a transformation is carried out using an equation relating the initial parameter of the problem to the length of the integral curve. The linear acceleration method is used as an example to demonstrate the procedure of constructing an implicit algorithm using simple iterations for the numerical solution of the transformed Cauchy problem. Propositions concerning the computational properties of the iterative process are formulated and proved. Explicit estimates are given for an integration stepsize that guarantees the convergence of the simple iterations. The efficacy of the proposed procedure is demonstrated by the numerical solution of three problems. A comparative analysis is carried out of the numerical solutions obtained with and without parametrization of the initial problems in these three settings. As a qualitative test the problem of the celestial mechanics of the “Pleiades” is considered. The second example is devoted to modelling the non-linear dynamics of an elastic flexible rod fixed at one end as a cantilever and coiled in its initial (static) state into a ring by a bending moment. The third example demonstrates the numerical solution of the problem of the “unfolding” of a mechanical system consisting of three flexible rods with given control input.     

Multi-objective Pareto-optimal control: an application to wastewater management

This work treats, within a multi-objective framework, of an economical-ecological problem related to the optimal management of a wastewater treatment system consisting of several purifying plants. The problem is formulated as a multi-objective parabolic optimal control problem and it is studied from a cooperative point of view, looking for Pareto-optimal solutions. The weighting method is used here to characterize the Pareto solutions of our problem. To obtain them, a numerical algorithm—based in a characteristics-Galerkin discretization—is proposed and numerical results for a real world situation in the estuary of Vigo (NW Spain) are also presented.     

Exact controllability of a spherical cap: Numerical implementation of HUM

This paper deals with the numerical implementation of the exact boundary controllability of the Reissner model for shallow spherical shells (Ref. 1). The problem is attacked by the Hilbert uniqueness method (HUM, Refs. 2–4), and we propose a semidiscrete method for the numerical approximation of the minimization problem associated to the exact controllability problem. The numerical results compare well with the results obtained by a finite difference and conjugate gradient method in Ref. 5.This work was done when the first two authors were at CNR-IAC, Rome, Italy as Graduate Students.     

An efficient method for solving non-linear singularly perturbed two points boundary-value problems of fractional order

In this paper, we discuss a numerical solution of a class of non-linear fractional singularly perturbed two points boundary-value problem. The method of solution consists of solving reduced problem and boundary layer correction problem. A series method is used to solve the boundary layer correction problem, and then the series solutions is approximated by the Pade’ approximant of order [m, m]. Some theoretical results are established and proved. Two numerical examples are discussed to illustrate the efficiency of the present scheme.     

The random projection method for a model problem of combustion with stiff chemical reactions

In this paper we extend the random projection method, recently proposed by the author and S. Jin [J. Comput. Phys. 163 (2000) 216] for under resolved numerical simulations of a qualitative model problem for combustion with stiff chemical reactions:

In this problem, the reaction time is small, making the problem numerically stiff. A classic spurious numerical phenomenon – the incorrect shock speed – occurs when the reaction time scale is not properly resolved numerically. The random projection method is introduced recently to handle this kind of numerical difficulty. The key idea in this method is to randomize the ignition temperature in a suitable domain. Several numerical experiments demonstrate the reliability and robustness of this method.     

带有不连续系数的线性输运方程差分格式的收敛性

提出了一个基于三角形网格的显式差分格式逼近带有不连续系数的线性输运方程. 通过对数值解的有界性、TVD(total variation decreasing)和空间、时间方向的平移估计, 利用Kolmogorov紧性原理证明了数值解在L¹_loc模下收敛于初值问题的唯一弱解.从而得到了初值问题解的存在唯一性和关于初值的稳定性. 数值算例表明本文提出的格式计算方便而且比 Lax-Friedrichs格式更有效.       

Numerical Computation of Multivariate Normal Probabilities

Abstract

The numerical computation of a multivariate normal probability is often a difficult problem. This article describes a transformation that simplifies the problem and places it into a form that allows efficient calculation using standard numerical multiple integration algorithms. Test results are presented that compare implementations of two algorithms that use the transformation with currently available software.     

A nonlinear optimal control problem for the nonsteady temperature of a layer with constraints on the thermal stress

We propose a numerical method of constructing the optimal heating regime for a thermally stressed unbounded layer with constraints on the control and thermal stresses. Solving the nonlinear optimization problem for rapidity is reduced to solving the inverse problem of thermoelasticity. The results of numerical studies are presented. Translated fromMatematichni Metody i Fiziko-Mekhanichni Polya, Vol. 38, 1995.     

Interface problem in holonomic elastoplasticity

The three-dimensional interface problem with the homogeneous Lamé system in an unbounded exterior domain and holonomic material behaviour in a bounded interior Lipschitz domain is considered. Existence and uniqueness of solutions of the interface problem are obtained rewriting the exterior problem in terms of boundary integral operators following the symmetric coupling procedure. The numerical approximation of the solutions consists in coupling of the boundary element method (BEM) and the finite element method (FEM). A Céa-like error estimate is presented for the discrete solutions of the numerical procedure proving its convergence.     

Direct numerical method for an inverse problem of a parabolic partial differential equation

A coefficient inverse problem of the one-dimensional parabolic equation is solved by a high-order compact finite difference method in this paper. The problem of recovering a time-dependent coefficient in a parabolic partial differential equation has attracted considerable attention recently. While many theoretical results regarding the existence and uniqueness of the solution are obtained, the development of efficient and accurate numerical methods is still far from satisfactory. In this paper a fourth-order efficient numerical method is proposed to calculate the function u(x,t) and the unknown coefficient a(t) in a parabolic partial differential equation. Several numerical examples are presented to demonstrate the efficiency and accuracy of the numerical method.     

Simple layer potential method for domains having external corners

This note investigates the simple layer potential method for domains having external corners. The singular behaviour of simple layer density at the corners is studied for the Neumann problem of Helmholtz's equation. A numerical technique of solving the integral equation for this problem is proposed. This technique takes the singularity of the solution into consideration. Some numerical examples are given to show the applicability of the present method.