The stability of the zero solution of a nonlinear nonautonomous Caputo fractional differential equation is studied using Lyapunov-like functions. The novelty of this paper is based on the new definition of the derivative of a Lyapunov-like function along the given fractional equation. Comparison results using this definition for scalar fractional differential equations are presented. Several sufficient conditions for stability, uniform stability and asymptotic uniform stability, based on the new definition of the derivative of Lyapunov functions and the new comparison result, are established.

I am presenting a survey of regularity results for both minima of variational integrals, and solutions to non-linear elliptic, and sometimes parabolic, systems of partial differential equations. I will try to take the reader to the Dark Side...

We propose a Halpern-type forward-backward splitting with inertial extrapolation step for finding a zero of the sum of accretive operators in Banach spaces. Strong convergence of the sequence of iterates generated by the method proposed is obtained under mild assumptions. We give some numerical results in compressed sensing to validate the theoretical analysis results. Our result is one of the lew available inertial-type methods for zeros of the sum of accretive operators in Banach spaces.

t=T − ψ(x)) 2 dx is applied and shown that this cost functional is Fréchet differentiable and its derivative can be formulated via the solution of an adjoint problem. In addition, Lipschitz continuity of the gradient is proved. These results help us to prove the monotonicity and convergence of the sequence {J′(f (n))}, where f (n) is the nth iteration of a gradient like method. At the end, the convexity of the Fréchet derivative is given.

Consider a second-order elliptic boundary value problem in three dimensions with locally smooth coefficients and solution. Discuss local superconvergence estimates for the tensor-product finite element approximation on a regular family of rectangular meshes. It will be shown that, by the estimates for the discrete Green’s function and discrete derivative Green’s function, and the relationship of norms in the finite element space such as L 2 -norms, W 1,∞ -norms, and negative-norms in locally smooth subsets of the domain Ω, locally pointwise superconvergence occurs in function values and derivatives.

The stress-strength model is proposed based on the m-generalized order statistics and the corresponding concomitant. For the dependency between m-generalized order statistics and its concomitant, a bivariate copula expansion is considered and the stressstrength model is obtained for two special cases of order statistics and upper record values. In the particular case of copula function, the generalized Farlie-Gumbel-Morgenstern bivariate distribution function is considered with proportional reversed hazard functions as marginal functions. Based on the order statistics and record values, two estimators of stress-strength are presented using a procedure similar to the inference functions for margins. Finally, a simulation study is carried out which shows the good performance of the proposed estimators for a finite sample.

Linear matrix approximation problems AX ≈ B are often solved by the total least squares minimization (TLS). Unfortunately, the TLS solution may not exist in general. The so-called core problem theory brought an insight into this effect. Moreover, it simplified the solvability analysis if B is of column rank one by extracting a core problem having always a unique TLS solution. However, if the rank of B is larger, the core problem may stay unsolvable in the TLS sense, as shown for the first time by Hnětynková, Plešinger, and Sima (2016). Full classification of core problems with respect to their solvability is still missing. Here we fill this gap. Then we concentrate on the so-called composed (or reducible) core problems that can be represented by a composition of several smaller core problems. We analyze how the solvability class of the components influences the solvability class of the composed problem. We also show on an example that the TLS solvability class of a core problem may be in some sense improved by its composition with a suitably chosen component. The existence of irreducible problems in various solvability classes is discussed.

We study an inverse eigenvalue problem (IEP) of reconstructing a special kind of symmetric acyclic matrices whose graph is a generalized star graph. The problem involves the reconstruction of a matrix by the minimum and maximum eigenvalues of each of its leading principal submatrices. To solve the problem, we use the recurrence relation of characteristic polynomials among leading principal minors. The necessary and sufficient conditions for the solvability of the problem are derived. Finally, a numerical algorithm and some examples are given.

Laplace interpolation is a popular approach in image inpainting using partial differential equations. The classic approach considers the Laplace equation with mixed boundary conditions. Recently a more general formulation has been proposed, where the differential operator consists of a point-wise convex combination of the Laplacian and the known image data. We provide the first detailed analysis on existence and uniqueness of solutions for the arising mixed boundary value problem. Our approach considers the corresponding weak formulation and aims at using the Theorem of Lax-Milgram to assert the existence of a solution. To this end we have to resort to weighted Sobolev spaces. Our analysis shows that solutions do not exist unconditionally. The weights need some regularity and must fulfil certain growth conditions. The results from this work complement findings which were previously only available for a discrete setup.

We study a method based on Balancing Domain Decomposition by Constraints (BDDC) for numerical solution of a single-phase flow in heterogeneous porous media. The method solves for both flux and pressure variables. The fluxes are resolved in three steps: the coarse solve is followed by subdomain solves and last we look for a divergence-free flux correction and pressures using conjugate gradients with the BDDC preconditioner. Our main contribution is an application of the adaptive algorithm for selection of flux constraints. Performance of the method is illustrated on the benchmark problem from the 10th SPE Comparative Solution Project (SPE 10). Numerical experiments in both 2D and 3D demonstrate that the first two steps of the method exhibit some numerical upscaling properties, and the adaptive preconditioner in the last step allows a significant decrease in the number of iterations of conjugate gradients at a small additional cost.

We consider the pricing of credit default swaps (CDSs) with the reference asset assumed to follow a geometric Brownian motion with a fast mean-reverting stochastic volatility, which is often observed in the financial market. To establish the pricing mechanics of the CDS, we set up a default model, under which the fair price of the CDS containing the unknown “no default” probability is derived first. It is shown that the “no default” probability is equivalent to the price of a down-and-out binary option written on the same reference asset. Based on the perturbation approach, we obtain an approximated but closed-form pricing formula for the spread of the CDS. It is also shown that the accuracy of our solution is in the order of $$\mathscr{O}(\varepsilon)$$ O ( ε ) .

Trust region methods are a class of effective iterative schemes in numerical optimization. In this paper, a new improved nonmonotone adaptive trust region method for solving unconstrained optimization problems is proposed. We construct an approximate model where the approximation to Hessian matrix is updated by the scaled memoryless BFGS update formula, and incorporate a nonmonotone technique with the new proposed adaptive trust region radius. The new ratio to adjusting the next trust region radius is different from the ratio in the traditional trust region methods. Under some suitable and standard assumptions, it is shown that the proposed algorithm possesses global convergence and superlinear convergence. Numerical results demonstrate that the proposed method is very promising.

This contribution presents the shape optimization problem of the plunger cooling cavity for the time dependent model of pressing the glass products. The system of the mould, the glass piece, the plunger and the plunger cavity is considered in four consecutive time intervals during which the plunger moves between 6 glass moulds.The state problem is represented by the steady-state Navier-Stokes equations in the cavity and the doubly periodic energy equation in the whole system, under the assumption of rotational symmetry, supplemented by suitable boundary conditions.The cost functional is defined as the squared weighted L 2 norm of the difference between a prescribed constant and the temperature of the plunger surface layer at the moment before separation of the plunger and the glass piece.The existence and uniqueness of the solution to the state problem and the existence of a solution to the optimization problem are proved.

We are interested in the numerical solution of a two-dimensional fluid-structure interaction problem. A special attention is paid to the choice of physically relevant inlet boundary conditions for the case of channel closing. Three types of the inlet boundary conditions are considered. Beside the classical Dirichlet and the do-nothing boundary conditions also a generalized boundary condition motivated by the penalization prescription of the Dirichlet boundary condition is applied. The fluid flow is described by the incompressible Navier-Stokes equations in the arbitrary Lagrangian-Eulerian (ALE) form and the elastic body creating a part of the channel wall is modelled with the aid of linear elasticity. Both models are coupled with the boundary conditions prescribed at the common interface.The elastic and the fluid flow problems are approximated by the finite element method. The detailed derivation of the weak formulation including the boundary conditions is presented. The pseudo-elastic approach for construction of the ALE mapping is used. Results of numerical simulations for three considered inlet boundary conditions are compared. The flutter velocity is determined for a specific model problem and it is shown that the boundary condition with the penalization approach is suitable for the case of the fluid flow in a channel with vibrating walls.

A pharmacodynamic model introduced earlier in the literature for in silico prediction of rifampicin-induced CYP3A4 enzyme production is described and some aspects of the involved curve-fitting based parameter estimation are discussed. Validation with our own laboratory data shows that the quality of the fit is particularly sensitive with respect to an unknown parameter representing the concentration of the nuclear receptor PXR (pregnane X receptor). A detailed analysis of the influence of that parameter on the solution of the model’s system of ordinary differential equations is given and it is pointed out that some ingredients of the analysis might be useful for more general pharmacodynamic models. Numerical experiments are presented to illustrate the performance of related parameter estimation procedures based on least-squares minimization.

We consider the implicit discretization of Nagumo equation on finite lattices and show that its variational formulation corresponds in various parameter settings to convex, mountain-pass or saddle-point geometries. Consequently, we are able to derive conditions under which the implicit discretization yields multiple solutions. Interestingly, for certain parameters we show nonuniqueness for arbitrarily small discretization steps. Finally, we provide a simple example showing that the nonuniqueness can lead to complex dynamics in which the number of bounded solutions grows exponentially in time iterations, which in turn implies infinite number of global trajectories.