Neumann Control of Unstable Parabolic Systems: Numerical Approach

The present article is concerned with the Neumann control of systems modeled by scalar or vector parabolic equations of reaction-advection-diffusion type with a particular emphasis on systems which are unstable if uncontrolled. To solve these problems, we use a combination of finite-difference methods for the time discretization, finite-element methods for the space discretization, and conjugate gradient algorithms for the iterative solution of the discrete control problems. We apply then the above methodology to the solution of test problems in two dimensions, including problems related to nonlinear models.     

External finite-element approximations of eigenfunctions in the case of multiple eigenvalues

The paper deals with the finite-element analysis of second-order elliptic eigenvalue problems when the approximate domains Ω_h are not subdomains of the original domain . The considerations are restricted to piecewise linear approximations. Special attention is devoted to the convergence of approximate eigenfunctions in the case of multiple exact eigenvalues. As yet the approximate solutions have been compared with linear combinations of exact eigenfunctions with coefficients depending on the mesh parameter h. We avoid this disadvantage.     

Free-Edge Stresses around Holes in Laminates by the Boundary Finite-Element Method

In the present paper the boundary finite-element method is presented as a highly efficient technique for the numerical investigation of the free-edge stresses around a circular hole in laminates. In this method, as in the boundary element method, only the boundary needs to be discretized, whereas the element formulation in essence is finite-element based. The surface discretization provides a high numerical efficiency and requires less computation time compared to finite-element analyses. Numerical results for the concentration of interlaminar stresses at holes in composite laminates show a very good agreement with comparative finite-element calculations.     

Numerical identification of the leading coefficient of a parabolic equation

For a multidimensional parabolic equation, we study the problem of finding the leading coefficient, which is assumed to depend only on time, on the basis of additional information about the solution at an interior point of the computational domain. For the approximate solution of the nonlinear inverse problem, we construct linearized approximations in time with the use of ordinary finite-element approximations with respect to space. The numerical algorithm is based on a special decomposition of the approximate solution for which the transition to the next time level is carried out by solving two standard elliptic problems. The capabilities of the suggested numerical algorithm are illustrated by the results of numerical solution of a model inverse two-dimensional problem.     

Error bounds for the finite-element approximation of the exterior Stokes equations in two dimensions

In this paper we design high-order (non)local artificial boundaryconditions (ABCs) which are different from those proposed byHan, H. & Bao, W. (1997 Numer. Math., 77, 347–363)for incompressible materials, and present error bounds for thefinite-element approximation of the exterior Stokes equationsin two dimensions. The finite-element approximation (especiallyits corresponding stiff matrix) becomes much simpler (sparser)when it is formulated in a bounded computational domain usingthe new (non)local approximate ABCs. Our error bounds indicatehow the errors of the finite-element approximations depend onthe mesh size, terms used in the approximate ABCs and the locationof the artificial boundary. Numerical examples of the exteriorStokes equations outside a circle in the plane are presented.Numerical results demonstrate the performance of our error bounds.     

Numerical schemes for the simulation of the two-dimensional Schrödinger equation using non-reflecting boundary conditions

This paper adresses the construction and study of a Crank-Nicolson-type discretization of the two-dimensional linear Schrödinger equation in a bounded domain with artificial boundary conditions set on the arbitrarily shaped boundary of . These conditions present the features of being differential in space and nonlocal in time since their definition involves some time fractional operators. After having proved the well-posedness of the continuous truncated initial boundary value problem, a semi-discrete Crank-Nicolson-type scheme for the bounded problem is introduced and its stability is provided. Next, the full discretization is realized by way of a standard finite-element method to preserve the stability of the scheme. Some numerical simulations are given to illustrate the effectiveness and flexibility of the method.

     

On exact and approximate boundary controllabilities for the heat equation: A numerical approach

The present article is concerned with the numerical implementation of the Hilbert uniqueness method for solving exact and approximate boundary controllability problems for the heat equation. Using convex duality, we reduce the solution of the boundary control problems to the solution of identification problems for the initial data of an adjoint heat equation. To solve these identification problems, we use a combination of finite difference methods for the time discretization, finite element methods for the space discretization, and of conjugate gradient and operator splitting methods for the iterative solution of the discrete control problems. We apply then the above methodology to the solution of exact and approximate boundary controllability test problems in two space dimensions. The numerical results validate the methods discussed in this article and clearly show the computational advantage of using second-order accurate time discretization methods to approximate the control problems.     

二维非线性对流扩散方程特征有限元的两重网络算法

针对二维非线性对流扩散方程,构造了特征有限元两重网格算法．该算法只需要在粗网格上进行非线性迭代运算,而在所需要求解的细网格上进行一次线性运算即可．对于非线性对流占优扩散方程,不仅可以消除因对流占优项引起的数值振荡现象,还可以加快收敛速度、提高计算效率．误差估计表明只要选取粗细网格步长满足一定的关系式,就可以使两重网格解与有限元解保持同样的计算精度．算例显示:两重网格算法比特征有限元算法的收敛速度明显加快．     

Projected equation methods for approximate solution of large linear systems

We consider linear systems of equations and solution approximations derived by projection on a low-dimensional subspace. We propose stochastic iterative algorithms, based on simulation, which converge to the approximate solution and are suitable for very large-dimensional problems. The algorithms are extensions of recent approximate dynamic programming methods, known as temporal difference methods, which solve a projected form of Bellman’s equation by using simulation-based approximations to this equation, or by using a projected value iteration method.     

Optimal control computation for parabolic systems with boundary conditions involving time delays

In this paper, we consider a class of optimal control problems involving a second-order, linear parabolic partial differential equation with Neumann boundary conditions. The time-delayed arguments are assumed to appear in the boundary conditions. A necessary and sufficient condition for optimality is derived, and an iterative method for solving this optimal control problem is proposed. The convergence property of this iterative method is also investigated.On the basis of a finite-element Galerkin's scheme, we convert the original distributed optimal control problem into a sequence of approximate problems involving only lumped-parameter systems. A computational algorithm is then developed for each of these approximate problems. For illustration, a one-dimensional example is solved.     

Accuracy of multivalue methods during change of time‐step size: Adams‐Moulton

Multivalue methods are a class of time‐stepping procedures for numerical solution of differential equations that progress to a new time level using the approximate solution for the function of interest and its derivatives at a single time level. The methods differ from multistep procedures, which make use of solutions to the differential equation at multiple time levels to advance to the new time level. Multistep methods are difficult to employ when a change in time‐step is desired, because the standard formulas (e.g., Adams‐Moulton or Gear) must be modified to accommodate the change. Multivalue methods seem to possess the desirable feature that the time‐step may be changed arbitrarily as one proceeds, since the solution proceeds from a single time level. However, in practice, changes in the time‐step introduce lower order errors or alter the coefficient in the truncation error term. Here, the multivalue Adams‐Moulton method is presented based on a general interpolation procedure. Modifications required to retain the high‐order accuracy of these methods during a change in time‐step are developed. Additionally, a formula for the unknown initial derivatives is presented. Finally, two examples are provided to illustrate the potential merit of the modification to the standard multivalue methods. © 2000 John Wiley & Sons, Inc. Numer Methods Partials Differential Eq 16: 312–326, 2000     

A mixed collocation–finite difference method for 3D microscopic heat transport problems

Three-dimensional time-dependent initial-boundary value problems of a novel microscopic heat equation are solved by the mixed collocation–finite difference method in and on the boundaries of a particle when the thickness is much smaller than both the length and width. The collocation method on fixed grid size is used to approximate the space operator, whereas the finite difference scheme is used for time discretization. This new mixed method is applied to a novel heat problem in a particle, in order to compute the temperature distribution in and on the particle's surface. The second derivatives of the basis functions for the spectral approximation are derived. Direct substitution of derivatives in the model transforms the differential equation into a linear system of equations that is solved by the specific preconditioned conjugate gradient method. The high-order accuracy and resolution achieved by the proposed method allows one to obtain engineering-accuracy solution on coarse meshes. The consistency, stability and convergence analysis are provided and numerical results are presented.     

Mathieu equation with application to analysis of dynamic characteristics of resonant inertial sensors

In this paper, Mathieu equation is applied to analyze the dynamic characteristics of resonant inertial sensors. Unlike previous work, Mathieu equation is not just a differential equation and analyzes the stability of the transition curves, but become an important method in analyzing parametric resonant characteristics and approximate output of resonant inertial sensors. It is demonstrated that the mathematical model of resonant inertial sensors is described by Mathieu equation. The relevant Mathieu equation theory and dynamic characteristics analysis methods were proposed, which include both stability and dynamic linear output. Finally, theoretical and experimental analysis show that the correlation of the theoretical curve and the experimental result coincide so perfectly, which means proposed analysis methods for Mathieu equation could be used to analyze the dynamic output characteristic of resonant inertial sensors. The theoretical analyzing approach of Mathieu equation and experimental results of resonant inertial sensors are obtained, which provide an application area for Mathieu equation and a reference for the robust design for resonant inertial sensors.     

B-splines methods with redefined basis functions for solving fourth order parabolic partial differential equations

In this work, we discuss two methods for solving a fourth order parabolic partial differential equation. In Method-I, we decompose the given equation into a system of second order equations and solve them by using cubic B-spline method with redefined basis functions. In Method-II, the equation is solved directly by applying quintic B-spline method with redefined basis functions. Stability of these methods have been discussed. Both methods are unconditionally stable. These methods are tested on four examples. The computed results are compared wherever possible with those already available in literature. We have developed Method-I for fourth order non homogeneous parabolic partial differential equation from which we can obtain displacement and bending moment both simultaneously, while Method-II gives only displacement. The results show that the derived methods are easily implemented and approximate the exact solution very well.     

A comparative study of fiber-matrix interactions by using micromechanical models

The purpose of this study is to describe the interfacial interactions in terms of stress distributions on short fibers in fiber-matrix unit-cell models. The fiber and matrix are subjected to tensile loading. The study consists of three main parts. First, fiber-matrix cell segments are modeled using a 3D finite-element analysis (FEA) with ANSYS. Three different finite-element geometrical unit-cell models are generated in order to simulate the Cox analytical model: a fiber-matrix combination, a single fiber, and a single matrix element. The second part contains the results of 3D FE analyses, which are applied to the Cox formulations by using a computer program developed. In the last part, the analytical solutions for distributions of normal and shear stresses are investigated. Cox 2D linear elasticity solutions, together with finite-element ones, are presented in detail in graphs. The interfacial interactions between the fibers and matrix are also discussed considering the relative changes in the distributions of normal and shear stresses. Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 44, No. 4, pp. 505–520, July–August, 2008.     

Adaptive multilevel methods in space and time for parabolic problems-the periodic case

The aim of this paper is to display numerical results that show the interest of some multilevel methods for problems of parabolic type. These schemes are based on multilevel spatial splittings and the use of different time steps for the various spatial components. The spatial discretization we investigate is of spectral Fourier type, so the approximate solution naturally splits into the sum of a low frequency component and a high frequency one. The time discretization is of implicit/explicit Euler type for each spatial component. Based on a posteriori estimates, we introduce adaptive one-level and multilevel algorithms. Two problems are considered: the heat equation and a nonlinear problem. Numerical experiments are conducted for both problems using the one-level and the multilevel algorithms. The multilevel method is up to 70% faster than the one-level method.

     

Nonlinear free dynamic response of laminated compressible cylindrical shell panels

The governing equations for a free dynamic response of a symmetrically laminated composite shell are used to analyze a nonlinear differential panel. The FEM and the Lindstedt–Poincare perturbation technique are invoked to construct a uniform asymptotic expansion of the solution to a nonlinear differential equation ofmotion. A comparison between numerical and finite-element methods for analyzing a symmetrically laminated graphite/epoxy shell panel is performed to show that the nonlinearities are of hardening type and are more repeated for smaller opening angles. It is also shown that large-amplitude motions are dominated by lower modes.     

Analytical solutions to the three-dimensional radio-nuclide transport equation for computer code verification

A series of analytical solutions to the three-dimensional radio-nuclide transport equation is derived and used to verify the numerical accuracy of the finite-element code NAMSOL. Good agreement is found for grid Peclet numbers up to 8. Above this value the solution becomes unstable as expected from the standard stability criterion.     

GLOBAL EXISTENCE OF WEAK SOLUTIONS TO THE NON-ISOTHERMAL NEMATIC LIQUID CRYSTALS IN 2D

In this article, we prove the global existence of weak solutions to the nonisothermal nematic liquid crystal system on T2, on the basis of a new approximate system which is different from the classical Ginzburg-Landau approximation. Local in space energy inequalities are employed to recover the estimates on the second order spatial derivatives of the director fields locally in time, which cannot be derived from the basic energy balance.It is shown that these weak solutions satisfy the temperature equation, and also the total energy equation but away from at most finite many "singular" times, at which the energy concentration occurs and the director field losses its second order derivatives.     

Multirate infinitesimal step methods for atmospheric flow simulation

The numerical solution of the Euler equations requires the treatment of processes in different temporal scales. Sound waves propagate fast compared to advective processes. Based on a spatial discretisation on staggered grids, a multirate time integration procedure is presented here generalising split-explicit Runge-Kutta methods. The advective terms are integrated by a Runge-Kutta method with a macro stepsize restricted by the CFL number. Sound wave terms are treated by small time steps respecting the CFL restriction dictated by the speed of sound.Split-explicit Runge-Kutta methods are generalised by the inclusion of fixed tendencies of previous stages. The stability barrier for the acoustics equation is relaxed by a factor of two.Asymptotic order conditions for the low Mach case are given. The relation to commutator-free exponential integrators is discussed. Stability is analysed for the linear acoustic equation. Numerical tests are executed for the linear acoustics and the nonlinear Euler equations.