A flexible and generic approach to dynamic modelling of supply chains

In this paper, we present a new modelling approach for realistic supply chain simulation. The model provides an experimental environment for informed comparison between different supply chain policies. A basic simulation model for a generic node, from which a supply chain network can be built, has been developed using an object-oriented approach. This generic model allows the incorporation of the information and physical systems and decision-making policies used by each node. The object-oriented approach gives the flexibility in specifying the supply chain configuration and operation decisions, and policies. Stochastic simulations are achieved by applying Latin Supercube Sampling to the uncertain variables in descending order of importance, which reduces the number of simulations required. We also present a case study to show that the model is applicable to a real-life situation for dynamic stochastic studies.     

On a class of preconditioners for solving the Helmholtz equation

In 1983, a preconditioner was proposed [J. Comput. Phys. 49 (1983) 443] based on the Laplace operator for solving the discrete Helmholtz equation efficiently with CGNR. The preconditioner is especially effective for low wavenumber cases where the linear system is slightly indefinite. Laird [Preconditioned iterative solution of the 2D Helmholtz equation, First Year's Report, St. Hugh's College, Oxford, 2001] proposed a preconditioner where an extra term is added to the Laplace operator. This term is similar to the zeroth order term in the Helmholtz equation but with reversed sign. In this paper, both approaches are further generalized to a new class of preconditioners, the so-called “shifted Laplace” preconditioners of the form Δφ−k²φ with . Numerical experiments for various wavenumbers indicate the effectiveness of the preconditioner. The preconditioner is evaluated in combination with GMRES, Bi-CGSTAB, and CGNR.     

Fermions in Odd Space-Time Dimensions: Back to Basics

It is a well-known feature of odd space-time dimensions d that there exist two inequivalent fundamental representations A and B of the Dirac gamma matrices. Moreover, the parity transformation swaps the fermion fields living in A and B. As a consequence, a parity-invariant Lagrangian can only be constructed by incorporating both the representations. Based upon these ideas and contrary to long-held belief, we show that in addition to a discrete exchange symmetry for the massless case, we can also define chiral symmetry provided the Lagrangian contains fields corresponding to both the inequivalent representations. We also study the transformation properties of the corresponding chiral currents under parity and charge-conjugation operations. We work explicitly in 2 + 1 dimensions and later show how some of these ideas generalize to an arbitrary number of odd dimensions.     

Heavy Majorana neutrinos in e + e − → W + W−

In this paper we investigate the effects of heavy Majorana neutrinos in the reaction e ⁺ e ^? → W ⁺ W ^?. We consider neutrino masses in the 1–10 TeV region. We show that at LEP II and NLC energies it is possible to use this processes to verify indirect evidence of heavy neutral particles with mixing angles of the order sin² α = 0.01. We discuss the unitarity restrictions that can be obtained for vector singlet and fermion-mirror-fermion models.     

Vibrational spectrum of uridine obtained by means of nonresonant-background-free CARS spectroscopy

Received: 18 June 1996/Revised version: 3 January 1997     

The multiple-scale wave equation

An analytic and numerical study of the behavior of the linear nonhomogeneous wave equation of the form ε²u_tt = Δu + tf with high wave speed (ε 1) is carried out. This study was initially motivated by meteorological observations which have indicated the presence of large spatial scale gravity waves in the neighborhood of a number of summer and winter storms, mainly from visible images of ripples in clouds in satellite photos. There is a question as to whether the presence of these waves is caused by the nearby storms. Since the linear wave equation is an approximation to the full system describing pressure waves in the atmosphere, yet is considerably more tractable, we have chosen to analyze the behavior of the linear nonhomogeneous wave equation with high wave speed. The analysis is shown to be valid in one, two, and three space dimensions. Partly because of the high wave speed, the solution is known to consist of behavior which changes on two different time scales, one rapid and one slow. Additionally, because of the presence of the nonhomogeneous forcing term tf, we show that there is a component of the solution which will vary only on a very large spatial scale. Since even the linearized wave equation can give rise to persistent large spatial scale waves under the right conditions, the implication is that certain storms could be responsible for the generation of large-scale waves. Numerical simulations in one and two dimensions confirm analytic results.     

Double-diffusive convection for a heated cylinder submerged in a salt-stratified fluid layer

Double-diffusive convection due to a cylindrical source submerged in a salt-stratified solution is numerically investigated in this study. For proper simulation of the vortex generated around the cylinder, a computational domain with irregular shape is employed. Flow conditions depend strongly on the thermal Rayleigh number, Ra _T , and the buoyancy ratio, R _ρ. There are two types of onset of instability existing in the flow field. Both types are due to either the interaction of the upward temperature gradient and downward salinity gradient or the interaction of the lateral temperature gradient and downward salinity gradient. The onset of layer instability due to plume convection is due to the former, whereas, the onset of layer instability of layers around the cylinder is due to the latter. Both types can be found in the flow field. The transport mechanism of layers at the top of the basic plume belongs to former while that due to basic plume and layer around the cylinder are the latter. The increase in Ra _T reinforces the plume convection and reduces the layer numbers generated around the cylinder for the same buoyancy ratio. For the same Ra _T , the increase of R _ρ suppresses the plume convection but reinforces the layers generated around the cylinder. The profiles of local Nusselt number reflects the heat transfer characteristics of plume convection and layered structure. The profiles of averaged Nusselt number are between the pure conduction and natural convection modes and the variation is due to the evolution of layers. Received on 13 September 1996     

Time-harmonic response of a vertical crack in plates

A technique for the characterization of an interior crack perpendicular to the surface of the plates is presented. A time-harmonic line load is applied on the upper surface of the plate. The scattered displacement field is calculated by using the strip element method. It is found that there are additional superimposed oscillations in the curve of the absolute values of the displacement in the region between the load and the crack. The mean value of the peaks of the additional oscillations increases with the crack length, and is less dependent on the distance from the load to the crack. The crack position is determined from the displacement distribution curve and the crack length is determined from the mean value of the peaks of the additional oscillations. Approximate polynomial formula of degree four for predicting the length of interior vertical cracks in isotropic homogeneous plates is also proposed. Numerical examples are given to demonstrate the present technique.