Superconvergence of Semidiscrete Splitting Positive Definite Mixed Finite Elements for Hyperbolic Optimal Control Problems

In this paper, we consider semidiscrete splitting positive definite mixed finite element methods for optimal control problems governed by hyperbolic equations with integral constraints. The state and costate are approximated by the lowest order Raviart-Thomas mixed rectangular finite element, and the control is approximated by piecewise constant functions. We derive some convergence and superconvergence results for the control, the state and the adjoint state. A numerical example is provided to demonstrate our theoretical results.

On the solvability of a semi-periodic boundary value problem for the nonlinear Goursat equation

In this paper, by means of a change of variables, a nonlinear semi-periodic boundary value problem for the Goursat equation is reduced to a linear gravity problem for hyperbolic equations. Reintroducing a new function, the obtained problem is reduced to a family of boundary value problems for ordinary differential equations and functional relations. When solving a family of boundary value problems for ordinary differential equations, the parameterization method is used. The application of this approach made it possible to establish the coefficients of the unique solvability of the semi-periodic problem for the Goursat equation and to propose constructive algorithms for finding an approximate solution.

Numerical simulation of thermal processes occurring in materials under the action of femtosecond laser pulses

In this work, a numerical study of the solutions of the parabolic and hyperbolic equations of heat conduction with the same physical parameters is carried out and a comparative analysis of the results obtained is carried out. The mathematical formulation of the problem is discussed. The action of the laser is taken into account through the source function, which was chosen as a double femtosecond laser pulse. In the hyperbolic equation, in contrast to the parabolic one, there is an additional parameter that characterizes the relaxation time of the heat flux. In addition, the source of the hyperbolic equation contains an additional term - the derivative of the power density of the source of the parabolic equation. This means that the temperature of the sample is influenced not only by the power density of the source, but also by the rate of its change. The profiles of the sample temperature at different times and its dynamics at different target depths are shown. The calculations were carried out for different time delays between pulses and for different relaxation parameters.

Age-structured delayed SIPCV epidemic model of HPV and cervical cancer cells dynamics I. Numerical method

The numerical method for simulation dynamics of nonlinear epidemic model of age-structured sub-populations of susceptible, infectious, precancerous and cancer cells and unstructured population of human papilloma virus (HPV) is developed (SIPCV model). Cell population dynamics is described by the initial-boundary value problem for the delayed semi-linear hyperbolic equations with age- and time-dependent coefficients and HPV dynamics is described by the initial problem for nonlinear delayed ODE. The model considers two time-delay parameters: the time between viral entry into a target susceptible cell and the production of new virus particles, and duration of the first stage of delayed immune response to HPV population growing. Using the method of characteristics and method of steps we obtain the exact solution of the SIPCV epidemic model in the form of explicit recurrent formulae. The numerical method designed for this solution and used the trapezoidal rule for integrals in recurrent formulae has a second order of accuracy. Numerical experiments with vanished mesh spacing illustrate the second order of accuracy of numerical solution with respect to the benchmark solution and show the dynamical regimes of cell-HPV population with the different phase portraits.