Recent higher-order explicit Runge–Kutta methods are compared with the classic fourth-order (RK4) method in long-term integration of both energy-conserving and lossy systems. By comparing quantity of function evaluations against accuracy for systems with and without known solutions, optimal methods are proposed. For a conservative system, we consider positional accuracy for Newtonian systems of two or three bodies and total angular momentum for a simplified Solar System model, over moderate astronomical timescales (tens of millions of years). For a nonconservative system, we investigate a relativistic two-body problem with gravitational wave emission. We find that methods of tenth and twelfth order consistently outperform lower-order methods for the systems considered here.

We consider the pricing of discretely sampled volatility swaps under a modified Heston model, whose risk-neutralized volatility process contains a stochastic long-run variance level. We derive an analytical forward characteristic function under this model, which has never been presented in the literature before. Based on this, we further obtain an analytical pricing formula for volatility swaps which can guarantee the computational accuracy and efficiency. We also demonstrate the significant impact of the introduced stochastic long-run variance level on volatility swap prices with synthetic as well as calibrated parameters.

When faced with the task of solving hyperbolic partial differential equations (PDEs), high order, strong stability-preserving (SSP) time integration methods are often needed to ensure preservation of the nonlinear strong stability properties of spatial discretizations. Among such methods, SSP second derivative time-stepping schemes have been recently introduced and used for evolving hyperbolic PDEs. In previous works, coupling of forward Euler and a second derivative formulation led to sufficient conditions for a second derivative general linear method (SGLM), which preserve the strong stability properties of spatial discretizations. However, for such methods, the types of spatial discretizations that can be used are limited. In this paper, we use a formulation based on forward Euler and Taylor series conditions to extend the SSP SGLM framework. We investigate the construction of SSP second derivative diagonally implicit multistage integration methods (SDIMSIMs) as a subclass of SGLMs with order $p=r=s$ and stage order $q=p,p-1$ up to order eight, where r is the number of external stages and s is the number of internal stages of the method. Proposed methods are examined on some one-dimensional linear and nonlinear systems to verify their theoretical order, and show potential of these schemes in preserving some nonlinear stability properties such as positivity and total variation.

This paper investigates spatial data on the unit sphere. Traditionally, isotropic Gaussian random fields are considered as the underlying mathematical model of the cosmic microwave background (CMB) data. We discuss the generalized multifractional Brownian motion and its pointwise Hölder exponent on the sphere. The multifractional approach is used to investigate the CMB data from the Planck mission. These data consist of CMB radiation measurements at narrow angles of the sky sphere. The results obtained suggest that the estimated Hölder exponents for different CMB regions do change from location to location. Therefore, the CMB temperature intensities are multifractional. The methodology developed is used to suggest two approaches for detecting regions with anomalies in the cleaned CMB maps.

We discuss a bi-objective two-stage assignment problem (BiTSAP) that aims at minimizing two objective functions: one comprising a nonlinear cost function defined explicitly in terms of assignment variables and the other a total completion time. A two-stage assignment problem deals with the optimal allocation of n jobs to n agents in two stages, where $n_1$ out of n jobs are primary jobs which constitute Stage-1 and the rest of the jobs are secondary jobs constituting Stage-2. The paper proposes an algorithm that seeks an optimal solution for a BiTSAP in terms of various efficient time-cost pairs. An algorithm for ranking all feasible assignments of a two-stage assignment problem in order of increasing total completion time is also presented. Theoretical justification and numerical illustrations are included to support the proposed algorithms.

We consider fully three-dimensional time-dependent outflow from a source into a surrounding fluid of different density. The source is distributed over a sphere of finite radius. The nonlinear problem is formulated using a spectral approach in which two streamfunctions and the density are represented as a Fourier-type series with time-dependent coefficients that must be calculated. Linearized theories are also discussed and an approximate stability condition for early stages in the outflow is derived. Nonlinear solutions are presented and different outflow shapes adopted by the fluid interface are investigated.

We explain some key challenges when dealing with a single- or multi-objective optimization problem in practice. To overcome these challenges, we present a mathematical program that optimizes the Nash social welfare function. We refer to this mathematical program as the Nash social welfare program (NSWP). An interesting property of the NSWP is that it can be constructed for any single- or multi-objective optimization problem. We show that solving the NSWP could result in more desirable solutions in practice than its single- or multi-objective counterpart. We also discuss several promising approaches that could be employed to solve the NSWP in practice.

This paper looks at adapting the method of Medvedev and Scaillet for pricing short-term American options to evaluate short-term convertible bonds. However unlike their method, we provide explicit formulae for the coefficients of our series solution. This means that we do not need to solve complicated recursive systems, and can efficiently provide fast solutions. We also compare the method with numerical solutions, and find that it performs extremely well, giving accurate bond prices as well as accurate optimal conversion prices.