Exponentially accurate solution tracking for nonlinear ODEs, the higher order Stokes phenomenon and double transseries resummation

We demonstrate the conjunction of new exponential-asymptotic effects in the context of a second order nonlinear ordinary differential equation with a small parameter. First, we show how to use a hyperasymptotic, beyond-all-orders approach to seed a numerical solver of a nonlinear ordinary differential equation with sufficiently accurate initial data so as to track a specific solution in the presence of an attractor. Second, we demonstrate the necessary role of a higher order Stokes phenomenon in analytically tracking the transition between asymptotic behaviours in a heteroclinic solution. Third, we carry out a double resummation involving both subdominant and sub-subdominant transseries to achieve the two-dimensional (in terms of the arbitrary constants) uniform approximation that allows the exploration of the behaviour of a two parameter set of solutions across wide regions of the independent variable. This is the first time all three effects have been studied jointly in the context of an asymptotic treatment of a nonlinear ordinary differential equation with a parameter. This paper provides an exponential asymptotic algorithm for attacking such problems when they occur. The availability of explicit results would depend on the individual equation under study

Semiconvection: numerical simulations

A grid of numerical simulations of double-diffusive convection is presented for the astrophysical case where viscosity (Prandtl number Pr) and solute diffusivity (Lewis number Le) are much smaller than the thermal diffusivity. As in laboratory and geophysical cases convection takes place in a layered form. The proper translation between subsonic flows in a stellar interior and an incompressible (Boussinesq) fluid is given, and the validity of the Boussinesq approximation for the semiconvection problem is checked by comparison with fully compressible simulations. The predictions of a simplified theory of mixing in semiconvection given in a companion paper are tested against the numerical results, and used to extrapolate these to astrophysical conditions. The predicted effective He-diffusion coefficient is nearly independent of the double-diffusive layering thickness $d$. For a fiducial main sequence model (15 $M_\odot$) the inferred mixing time scale is of the order $10^{10}$ yr. An estimate for the secular increase of $d$ during the semiconvective phase is given. It can potentially reach a significant fraction of a pressure scale height.Comment: arXiv admin note: substantial text overlap with arXiv:1012.585

A Matrix Element for Chaotic Tunnelling Rates and Scarring Intensities

It is shown that tunnelling splittings in ergodic double wells and resonant widths in ergodic metastable wells can be approximated as easily-calculated matrix elements involving the wavefunction in the neighbourhood of a certain real orbit. This orbit is a continuation of the complex orbit which crosses the barrier with minimum imaginary action. The matrix element is computed by integrating across the orbit in a surface of section representation, and uses only the wavefunction in the allowed region and the stability properties of the orbit. When the real orbit is periodic, the matrix element is a natural measure of the degree of scarring of the wavefunction. This scarring measure is canonically invariant and independent of the choice of surface of section, within semiclassical error. The result can alternatively be interpretated as the autocorrelation function of the state with respect to a transfer operator which quantises a certain complex surface of section mapping. The formula provides an efficient numerical method to compute tunnelling rates while avoiding the need for the exceedingly precise diagonalisation endemic to numerical tunnelling calculations.Comment: Submitted to Annals of Physics. This work has been submitted to Academic Press for possible publicatio

On Hirth Ring Couplings: Design Principles Including the Effect of Friction

Rings with Hirth couplings are primarily used for the accurate positioning of axial-symmetric components in the machine tool industry and, generally, in mechanical components. It is also possible to use Hirth rings as connection tools. Specific industries with special milling and grinding machines are able to manufacture both tailor made and standard Hirth rings available on stock. Unfortunately, no international standard (for instance ISO, DIN or AGMA) is available for the production and the design of such components. In the best-case scenario, it is possible to find simplified design formulae in the catalogue of the suppliers. The aim of this work is to provide some accurate formulae and computational methods for design to provide better awareness on the limitations and the potential of this type of connection. The work consists of five parts: (i) a review of the base calculation derived mainly from the catalogues of manufacturers; (ii) an improved calculation based on a new analytical method including the friction phenomenon; (iii) an experimentation run for validating the method; (iv) a case study applied to a machine tool; and, (v) a closed form formulation to determine an upper threshold for friction, thus ensuring the Hirth coupling regular performance

Validation of Simulation: Patterns in the Social and Natural Sciences

In most cases, the meaning of computer simulation is strongly connected to the idea numerical calculations. A computer simulation is a numerical solution of a complex mathematical problem. Therefore, the problem of validation of its results should be only a problem of judging the underlying computational methods. However, it will be argued, that this is not the case. It is consensus in literature that validation constitutes one of the central epistemological problems of computer simulation methods. Especially in the case of simulations in the social sciences the answers given by many authors are not satisfactory. The following paper attempts to show how the characteristics of simulation, i.e. the imitation of a dynamic, constitute the problem of validation even in the case of the natural sciences and what consequences arise. Differences as well as common grounds between social and natural sciences will be discussed.Generative Mechanism, Imitation, Patterns, Simulation, Validation

A study of jet impingment on curved surfaces followed by oblique introduction into a freestream flow

Technology of thrust reversers with particular application to STOL aircraf

Numerical calculation of Bessel, Hankel and Airy functions

The numerical evaluation of an individual Bessel or Hankel function of large order and large argument is a notoriously problematic issue in physics. Recurrence relations are inefficient when an individual function of high order and argument is to be evaluated. The coefficients in the well-known uniform asymptotic expansions have a complex mathematical structure which involves Airy functions. For Bessel and Hankel functions, we present an adapted algorithm which relies on a combination of three methods: (i) numerical evaluation of Debye polynomials, (ii) calculation of Airy functions with special emphasis on their Stokes lines, and (iii) resummation of the entire uniform asymptotic expansion of the Bessel and Hankel functions by nonlinear sequence transformations. In general, for an evaluation of a special function, we advocate the use of nonlinear sequence transformations in order to bridge the gap between the asymptotic expansion for large argument and the Taylor expansion for small argument ("principle of asymptotic overlap"). This general principle needs to be strongly adapted to the current case, taking into account the complex phase of the argument. Combining the indicated techniques, we observe that it possible to extend the range of applicability of existing algorithms. Numerical examples and reference values are given.Comment: 18 pages; 7 figures; RevTe

Trajectory generation for road vehicle obstacle avoidance using convex optimization

This paper presents a method for trajectory generation using convex optimization to find a feasible, obstacle-free path for a road vehicle. Consideration of vehicle rotation is shown to be necessary if the trajectory is to avoid obstacles specified in a fixed Earth axis system. The paper establishes that, despite the presence of significant non-linearities, it is possible to articulate the obstacle avoidance problem in a tractable convex form using multiple optimization passes. Finally, it is shown by simulation that an optimal trajectory that accounts for the vehicle’s changing velocity throughout the manoeuvre is superior to a previous analytical method that assumes constant speed