Modeling Multivariate Interest Rates Using Time-Varying Copulas and Reducible Nonlinear Stochastic Differential Equations

We propose a new approach for modeling nonlinear multivariate interest rate processes based on time-varying copulas and reducible stochastic differential equations (SDEs). In the modeling of the marginal processes, we consider a class of nonlinear SDEs that are reducible to Ornstein--Uhlenbeck (OU) process or Cox, Ingersoll, and Ross (1985) (CIR) process. The reducibility is achieved via a nonlinear transformation function. The main advantage of this approach is that these SDEs can account for nonlinear features, observed in short-term interest rate series, while at the same time leading to exact discretization and closed-form likelihood functions. Although a rich set of specifications may be entertained, our exposition focuses on a couple of nonlinear constant elasticity volatility (CEV) processes, denoted as OU-CEV and CIR-CEV, respectively. These two processes encompass a number of existing models that have closed-form likelihood functions. The transition density, the conditional distribution function, and the steady-state density function are derived in closed form as well as the conditional and unconditional moments for both processes. In order to obtain a more flexible functional form over time, we allow the transformation function to be time varying. Results from our study of U.S. and UK short-term interest rates suggest that the new models outperform existing parametric models with closed-form likelihood functions. We also find the time-varying effects in the transformation functions statistically significant. To examine the joint behavior of interest rate series, we propose flexible nonlinear multivariate models by joining univariate nonlinear processes via appropriate copulas. We study the conditional dependence structure of the two rates using Patton (2006a) time-varying symmetrized Joe--Clayton copula. We find evidence of asymmetric dependence between the two rates, and that the level of dependence is positively related to the level of the two rates. (JEL: C13, C32, G12) Copyright The Author 2010. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oxfordjournals.org, Oxford University Press.

Discretization and solution of the inhomogeneous Bogoliubov-de Gennes equations

In the presence of inhomogeneities, defects and currents, the equations describing a Bose-condensed ensemble of alkali atoms have to be solved numerically. By combining both linear and nonlinear equations within a Discrete Variable Representation framework, we describe a computational scheme for the solution of the coupled Bogoliubov-de Gennes (BdG) and nonlinear Schrodinger (NLS) equations for fields in a 3D spheroidal potential. We use the method to calculate the collective excitation spectrum and quasiparticle mode densities for excitations of a Bose condensed gas in a spheroidal trap. The method is compared against finite-difference and spectral methods, and we find the DVR computational scheme to be superior in accuracy and efficiency for the cases we consider. (C) 2004 Elsevier B.V. All rights reserved.

A review of the equations used to predict the velocity distribution within a ship’s propeller jet

Predicting the velocity within the ship’s propeller jet is the initial step to investigate the scouring made by the propeller jet. Albertson et al. (1950) suggested the investigation of a submerged jet can be undertaken through observation of the plain water jet from an orifice. The plain water jet investigation of Albertson et al. (1950) was based on the axial momentum theory. This has been the basis of all subsequent work with propeller jets. In reality, the velocity characteristic of a ship’s propeller jet is more complicated than a plain water jet. Fuehrer and Römisch (1977), Blaauw and van de Kaa (1978), Berger et al. (1981), Verhey (1983) and Hamill (1987) have carried out investigations using physical model. This paper reviews the state-of-art of the equations used to predict the time-averaged axial, tangential and radial components of velocity within the zone of flow establishment and the zone of established flow of a ship’s propeller jet.

Critical assessment of two-qubit post-Markovian master equations

A post-Markovian master equation has been recently proposed as a tool to describe the evolution of a system coupled to a memory-keeping environment [A. Shabani and D. A. Lidar, Phys. Rev. A 71, 020101 ( R) ( 2005)]. For a single qubit affected by appropriately chosen environmental conditions, the corresponding dynamics is always legitimate and physical. Here we extend such a situation to the case of two qubits, only one of which experiences the environmental effects. We show how, despite the innocence of such an extension, the introduction of the second qubit should be done cum grano salis to avoid consequences such as the breaking of the positivity of the associated dynamical map. This hints at the necessity of using care when adopting phenomenologically derived models for evolutions occurring outside the Markovian framework.

Cold-formed steel sections with web openings subjected to web crippling under two-flange loading conditions - Part II: Parametric study and proposed design equations

A parametric study of cold-formed steel sections with web openings subjected to web crippling was undertaken using finite element analysis, to investigate the effects of web holes and cross-section sizes on the web crippling strengths of channel sections subjected to web crippling under both interior-two-flange (ITF) and end-two-flange (ETF) loading conditions. In both loading conditions, the hole was centred beneath the bearing plate. It was demonstrated that the main factors influencing the web crippling strength are the ratio of the hole depth to the flat depth of the web, and the ratio of the length of bearing plates to the flat depth of the web. In this paper, design recommendations in the form of web crippling strength reduction factors are proposed, that are conservative to both the experimental and finite element results.

Optimal proportional-integral-derivative control of shape memory alloy actuators using genetic algorithm and the Preisach model

Shape memory alloy (SMA) actuators, which have the ability to return to a predetermined shape when heated, have many potential applications in aeronautics, surgical tools, robotics, and so on. Although the number of applications is increasing, there has been limited success in precise motion control owing to the hysteresis effect of these smart actuators. The present paper proposes an optimization of the proportional-integral-derivative (PID) control method for SMA actuators by using genetic algorithm and the Preisach hysteresis model.

R-matrix-Floquet theory of multiphoton processes:VIII. A linear equations method

A new linear equations method for calculating the R-matrix, which arises in the R-matrix-Floquet theory of multiphoton processes, is introduced. This method replaces the diagonalization of the Floquet Hamiltonian matrix by the solution of a set of linear simultaneous equations which are solved, in the present work, by the conjugate gradient method. This approach uses considerably less computer memory and can be readily ported onto parallel computers. It will thus enable much larger problems of current interest to be treated. This new method is tested by applying it to three-photon ionization of helium at frequencies where double resonances with a bound state and autoionizing states are important. Finally, an alternative linear equations method, which avoids the explicit calculation of the R-matrix by incorporating the boundary conditions directly, is described in an appendix.

A comparison of stress intensity factors through crack closure integral and other approaches using eXtended element-free Galerkin method

In this paper, a new approach for extracting stress intensity factors (SIFs) by the extended element-free Galerkin method, through a crack closure integral (CCI) scheme, is proposed. The CCI calculation is used in conjunction with a local smoothing technique to improve the accuracy of the computed SIFs in a number of case studies of linear elastic fracture mechanics. The cases involve problems of mixed-mode, curved crack and thermo-mechanical loading. The SIFs by CCI, displacement and stress methods are compared with those based on the M-integral technique reported in the literature. The proposed CCI method involves very simple relations, and still gives good accuracy. The convergence of the results is also examined.

Modified crack closure integral technique for extraction of SIFs in meshfee methods

A new approach for extracting stress intensity factors (SIFs) by the element-free Galerkin (EFG) class of methods through a modified crack closure integral (MCCI) scheme is proposed. Its primary feature is that it allows accurate calculation of mode I and mode II SIFs with a relatively simple and straightforward analysis even when a coarser nodal density is employed. The details of the adoption of the MCCI technique in the EFG method are described. Its performance is demonstrated through a number of case studies including mixed-mode and thermal problems in linear elastic fracture mechanics (LEFM). The results are compared with published theoretical solutions and those based on the displacement method, stress method, crack closure integral in conjunction with local smoothing (CCI–LS) technique, as well as the M-integral method. Its advantages are discussed.

Sound Synthesis for Contact-Driven Musical Instruments via Discretisation of Hamilton's Equations

Physical modelling of musical instruments involves studying nonlinear interactions between parts of the instrument. These can pose several difficulties concerning the accuracy and stability of numerical algorithms. In particular, when the underlying forces are non-analytic functions of the phase-space variables, a stability proof can only be obtained in limited cases. An approach has been recently presented by the authors, leading to unconditionally stable simulations for lumped collision models. In that study, discretisation of Hamilton’s equations instead of the usual Newton’s equation of motion yields a numerical scheme that can be proven to be energy conserving. In this paper, the above approach is extended to collisions of distributed objects. Namely, the interaction of an ideal string with a flat barrier is considered. The problem is formulated within the Hamiltonian framework and subsequently discretised. The resulting nonlinearmatrix equation can be shown to possess a unique solution, that enables the update of the algorithm. Energy conservation and thus numerical stability follows in a way similar to the lumped collision model. The existence of an analytic description of this interaction allows the validation of the model’s accuracy. The proposed methodology can be used in sound synthesis applications involving musical instruments where collisions occur either in a confined (e.g. hammer-string interaction, mallet impact) or in a distributed region (e.g. string-bridge or reed-mouthpiece interaction).