Computation of limit periodic continued fractions. A survey

Over the last 20 years a large number of algorithms has been published to improve the speed and domain of convergence of continued fractions. In this survey we show that these algorithms are strongly related. Actually, they essentially boil down to two main principles.We also prove some results on asymptotic expansions of tail values of limit periodic continued fractions.Dedicated to Luigi Gatteschi on his seventieth birthdayThis research was partially supported by The Norwegian Research Council and by the HMC project ROLLS, under contract CHRX-CT93-0416.     

Radiation bursts from particles in the field of compact, impenetrable, astrophysical objects

The radiation emitted by charged, scalar particles in a Schwarzschild field with maximal acceleration corrections is calculated classically and in the tree approximation of quantum field theory. In both instances the particles emit radiation that has characteristics similar to those of gamma-ray bursters.     

AGE method simulations of a turbulent far‐wake compared to spectral DNS

Turbulent flow simulation methods based on finite differences are attractive for their simplicity, flexibility and efficiency, but not always for accuracy or stability. This paper demonstrates that a good compromise is possible with the advected grid explicit (AGE) method. Starting from the same initial field as a previous spectral DNS, AGE method simulations of a planar turbulent wake were carried out as DNS, and then at three levels of reduced resolution. The latter cases were in a sense large‐eddy simulations (LES), although no specific sub‐grid‐scale model was used. Results for the two DNS methods, including variances and power spectra, were very similar, but the AGE simulation required much less computational effort. Small‐scale information was lost in the reduced resolution runs, but large‐scale mean and instantaneous properties were reproduced quite well, with further large reductions in computational effort. Quality of results becomes more sensitive to the value chosen for one of the AGE method parameters as resolution is reduced, from which it is inferred that the numerical stability procedure controlled by the parameter is acting in part as a sub‐grid‐scale model. Copyright © 2002 John Wiley & Sons, Ltd.     

Accelerating Effect in Kerr-Newman-Kasuya Field

We have obtained expressions of the accelerating effect in Kerr-Newman Kasuya field. These expressionsinclude four parameters: mass m, angular momentum a, electric charge q, and magnetic charge φ. Furthermore we studyits special case (vi = 0). We get the following conclusion. In the gravitation field of souse mass with electric charge qand magnetic charge b, the acceleration of test particle has not only radial component but also transverse component.When θ = 0, the acceleration is minimum, and when θ = π/2, the acceleration is maximum. Furthermore, we discussthe effects of electric charge q and magnetic charge φ respectively.     

A model study for the steady state error of numerical approximations of inviscid flow equations on Cartesian grid

The widely used locally adaptive Cartesian grid methods involve a series of abruptly refined interfaces. In this paper we consider the influence of the refined interfaces on the steady state errors for second‐order three‐point difference approximations of flow equations. Since the various characteristic components of the Euler equations should behave similarly on such grids with regard to refinement‐induced errors, it is sufficient enough to conduct the analysis on a scalar model problem. The error we consider is a global error, different to local truncation error, and reflects the interaction between multiple interfaces. The steady state error will be compared to the errors on smooth refinement grids and on uniform grids. The conclusion seems to support the numerical findings of Yamaleev and Carpenter (J. Comput. Phys. 2002; 181: 280–316) that refinement does not necessarily reduce the numerical error. Copyright © 2005 John Wiley & Sons, Ltd.     

Stability and Nonlinear Deformation of Cylindrical Grids

The general stability of single-layer cylindrical grids is studied in linear and nonlinear formulations. Dependence of the general buckling load on the geometry and stiffness parameters of a grid is established in an analytical form. Such grids are numerically analyzed for stability. It is established that the general buckling load is much less than the local buckling load. Typical general buckling modes are found. It is shown that such grids are weakly sensitive to imperfections     

Geometry of an Accelerated Rotating Disk

We analyze the geometry of a rotating disk with a tangential acceleration in the framework of the theory of Special Relativity, using the kinematic linear differential system that verifies the relative position vector of time-like curves in a Fermi reference. A numerical integration of these equations for a generic initial value problem is made up and the results are compared with those obtained in other works.     

Numerical Modeling of Interior Ballistics Processes in Light Gas Guns

An improved procedure for modeling the operation of a light-gas gun is proposed. The motion of working bodies in both the firing chamber and the light-gas chamber is studied within the framework of the mechanics of heterogeneous media. The problem of barrel heating taking into account its melting and removal of thermal ablation products into the medium inside the bore is solved in a coupled formulation. Heat and mass transfer and friction on the barrel surface are calculated using empirical dependences. The deformable piston is considered compressible and elastoviscoplastic. Allowance is made for the presence of a clearance between the lateral surface of the piston and the barrel bore walls and the associated gas flow between the firing and the light-gas chamber. Calculation results are given.     

Computations of two passing‐by high‐speed trains by a relaxation overset‐grid algorithm

This paper presents a relaxation algorithm, which is based on the overset grid technology, an unsteady three‐dimensional Navier–Stokes flow solver, and an inner‐ and outer‐relaxation method, for simulation of the unsteady flows of moving high‐speed trains. The flow solutions on the overlapped grids can be accurately updated by introducing a grid tracking technique and the inner‐ and outer‐relaxation method. To evaluate the capability and solution accuracy of the present algorithm, the computational static pressure distribution of a single stationary TGV high‐speed train inside a long tunnel is investigated numerically, and is compared with the experimental data from low‐speed wind tunnel test. Further, the unsteady flows of two TGV high‐speed trains passing by each other inside a long tunnel and at the tunnel entrance are simulated. A series of time histories of pressure distributions and aerodynamic loads acting on the train and tunnel surfaces are depicted for detailed discussions. Copyright © 2004 John Wiley & Sons, Ltd.     

A time‐accurate pseudo‐compressibility approach based on an unstructured hybrid finite volume technique applied to unsteady turbulent premixed flame propagation

A preconditioning approach based on the artificial compressibility formulation is extended to solve the governing equations for unsteady turbulent reactive flows with heat release, at low Mach numbers, on an unstructured hybrid grid context. Premixed reactants are considered and a flamelet approach for combustion modelling is adopted using a continuous quenched mean reaction rate. An overlapped cell‐vertex finite volume method is adopted as a discretisation scheme. Artificial dissipation terms for hybrid grids are explicitly added to ensure a stable, discretised set of equations. A second‐order, explicit, hybrid Runge–Kutta scheme is applied for the time marching in pseudo‐time. A time derivative of the dependent variable is added to recover the time accuracy of the preconditioned set of equations. This derivative is discretised by an implicit, second‐order scheme. The resulting scheme is applied to the calculation of an infinite planar (one‐dimensional) turbulent premixed flame propagating freely in reactants whose turbulence is supposed to be frozen, homogeneous and isotropic. The accuracy of the results obtained with the proposed method proves to be excellent when compared to the data available in the literature. Copyright © 2004 John Wiley & Sons, Ltd.