On the numerical solution of the 2d wave equation with compact fdtd schemes

This paper discusses compact-stencil finite difference time domain (FDTD) schemes for approximating the 2D wave equation in the context of digital audio. Stability, accuracy, and efficiency are investigated and new ways of viewing and interpreting the results are discussed. It is shown that if a tight accuracy constraint is applied, implicit schemes outperform explicit schemes. The paper also discusses the relevance to digital waveguide mesh modelling, and highlights the optimally efficient explicit scheme.

Doppler cooling of gallium atoms: 2. Simulation in complex multilevel systems

This paper derives a general procedure for the numerical solution of the Lindblad equations that govern the coherences arising from multicoloured light interacting with a multilevel system. A systematic approach to finding the conservative and dissipative terms is derived and applied to the laser cooling of p-block elements. An improved numerical method is developed to solve the time-dependent master equation and results are presented for transient cooling processes. The method is significantly more robust, efficient and accurate than the standard method and can be applied to a broad range of atomic and molecular systems. Radiation pressure forces and the formation of dynamic dark states are studied in the gallium isotope 66Ga.

Two-dimensional hygrothermal transfer in porous building materials

A two-dimensional mathematical model for evaluating the simultaneous heat and moisture migration in porous building materials was proposed. Vapor content and temperature were chosen as the principal driving potentials. The numerical solution was based on the control volume finite difference technique with fully implicit scheme in time. Two validation experiments were developed in this study. The evolution of transient moisture distributions in both one-dimensional and two-dimensional cases was measured. A comparison between experimental results and those obtained by the numerical model proves that they are fully consistent with each other. The model can be easily integrated into a whole building heat, air and moisture transfer model. Another main advantage of the present numerical method lies in the fact that the required moisture transport properties are comparatively simple and easy to determine.

Boson pairs in a one-dimensional split trap

We describe the properties of a pair of ultracold bosonic atoms in a one-dimensional harmonic trapping potential with a tunable zero-ranged barrier at the trap center. The full characterization of the ground state is done by calculating the reduced single-particle density, the momentum distribution, and the two-particle entanglement. We derive several analytical expressions in the limit of infinite repulsion (Tonks-Girardeau limit) and extend the treatment to finite interparticle interactions by numerical solution. As pair interactions in double wells form a fundamental building block for many-body systems in periodic potentials, our results have implications for a wide range of problems.

Improved conjugate gradient implementation for least squares support vector machines

As a promising method for pattern recognition and function estimation, least squares support vector machines (LS-SVM) express the training in terms of solving a linear system instead of a quadratic programming problem as for conventional support vector machines (SVM). In this paper, by using the information provided by the equality constraint, we transform the minimization problem with a single equality constraint in LS-SVM into an unconstrained minimization problem, then propose reduced formulations for LS-SVM. By introducing this transformation, the times of using conjugate gradient (CG) method, which is a greatly time-consuming step in obtaining the numerical solution, are reduced to one instead of two as proposed by Suykens et al. (1999). The comparison on computational speed of our method with the CG method proposed by Suykens et al. and the first order and second order SMO methods on several benchmark data sets shows a reduction of training time by up to 44%. (C) 2011 Elsevier B.V. All rights reserved.

The Reliability of the Arc-Length Method in the Analysis of Mode-Jumping Problems

The Arc-Length Method is a solution procedure that enables a generic non-linear problem to pass limit points. Some examples are provided of mode-jumping problems solutions using a commercial nite element package, and other investigations are carried out on a simple structure of which the numerical solution can be compared with an analytical one. It is shown that Arc-Length Method is not reliable when bifurcations are present in the primary equilibrium path; also the presence of very sharp snap-backs or special boundary conditions may cause convergence diÆculty at limit points. An improvement to the predictor used in the incremental procedure is suggested, together with a reliable criteria for selecting either solution of the quadratic arc-length constraint. The gap that is sometimes observed between the experimantal load level of mode-jumping and its arc-length prediction is explained through an example.

2DRMP-G: Migrating a large-scale numerical mathematical application to a grid environment

We report on the migration of a traditional, single architecture application to a grid application using heterogeneous resources. We focus on the use of the UK e-Science Level 2 grid (UKL2G) which provides a heterogeneous collection of resources distributed within the UK. We discuss the solution architecture, the performance of our application, its future development as a grid-based application and comment on the lessons we have learned in using a grid infrastructure for large-scale numerical problems.

Buckling and Postbuckling Structures: Experimental, Analytical and Numerical Studies

This book provides an in-depth treatment of the study of the stability of engineering structures. Contributions from internationally recognized leaders in the field ensure a wide coverage of engineering disciplines in which structural stability is of importance, in particular the analytical and numerical modelling of structural stability applied to aeronautical, civil, marine and offshore structures. The results from a number of comprehensive experimental test programs are also presented, thus enhancing our understanding of stability phenomena as well as validating the analytical and computational solution schemes presented. A variety of structural materials are investigated with special emphasis on carbon-fibre composites, which are being increasingly utilized in weight-critical structures. Instabilities at the meso- and micro-scales are also discussed. This book will be particularly relevant to professional engineers, graduate students and researchers interested in structural stability.

A numerical study of automotive turbocharger mixed flow turbine inlet geometry for off design performance

Mixed flow turbines represent a potential solution to the increasing requirement for high pressure, low velocity ratio operation in turbocharger applications. While literature exists for the use of these turbines at such operating conditions, there is a lack of detailed design guidance for defining the basic geometry of the turbine, in particular, the cone angle – the angle at which the inlet of the mixed flow turbine is inclined to the axis. This investigates the effect and interaction of such mixed flow turbine design parameters.
Computational Fluids Dynamics was initially used to investigate the performance of a modern radial turbine to create a baseline for subsequent mixed flow designs. Existing experimental data was used to validate this model.
Using the CFD model, a number of mixed flow turbine designs were investigated. These included studies varying the cone angle and the associated inlet blade angle.
The results of this analysis provide insight into the performance of a mixed flow turbine with respect to cone and inlet blade angle.

Temperature fields in a channel partially filled with a porous material under local thermal non-equilibrium condition-an exact solution

This work examines analytically the forced convection in a channel partially filled with a porous material and subjected to constant wall heat flux. The Darcy–Brinkman–Forchheimer model is used to represent the fluid transport through the porous material. The local thermal non-equilibrium, two-equation model is further employed as the solid and fluid heat transport equations. Two fundamental models (models A and B) represent the thermal boundary conditions at the interface between the porous medium and the clear region. The governing equations of the problem are manipulated, and for each interface model, exact solutions, for the solid and fluid temperature fields, are developed. These solutions incorporate the porous material thickness, Biot number, fluid to solid thermal conductivity ratio and Darcy number as parameters. The results can be readily used to validate numerical simulations. They are, further, applicable to the analysis of enhanced heat transfer, using porous materials, in heat exchangers.