A relaxation scheme for conservation laws with a discontinuous coefficient

We study a relaxation scheme of the Jin and Xin type for conservation laws with a flux function that depends discontinuously on the spatial location through a coefficient . If , we show that the relaxation scheme produces a sequence of approximate solutions that converge to a weak solution. The Murat-Tartar compensated compactness method is used to establish convergence. We present numerical experiments with the relaxation scheme, and comparisons are made with a front tracking scheme based on an exact Riemann solver.     

Computational Modeling of Organizations Comes of Age

As they are maturing—i.e., as they are becoming validated, calibrated and refined—computational emulation models of organizations are evolving into: powerful new kinds of organizational design tools for predicting and mitigating organizational risks; and flexible new kinds of organizational theorem-provers for validating extant organization theory and developing new theory. Over the past 50 years, computational modeling and simulation have had enormous impacts on the rate of advancement of knowledge in fields like physics, chemistry and, more recently, biology; and their subsequent application has enabled whole new areas of engineering practice. In the same way, as our young discipline comes of age, computational organizational models are beginning to impact behavioral, organizational and economic science, and management consulting practice. This paper attempts to draw parallels between computational modeling in natural sciences and computational modeling of organizations as a contributor to both social science and management practice.To illustrate the lifecycle of a computational organizational model that is now relatively mature, this paper traces the evolution of the Virtual Design Team (VDT) computational modeling and simulation research project at Stanford University from its origins in 1988 to the present. It lays out the steps in the process of validating VDT as a computational emulation model of organizations to the point that VDT began to influence management practice and, subsequently, to advance organizational science. We discuss alternate research trajectories that can be taken by computational and mathematical modelers who prefer the typical natural science validation trajectory—i.e., who attempt to impact organizational science first and, perhaps subsequently, to impact management practice.The paper concludes with a discussion of the current state-of-the-art of computational modeling of organizations and some thoughts about where, and how rapidly, the field is headed.     

Non-destructive testing of tubes using a time reverse numerical simulation (TRNS) method

A method for the detection of defects in cylindrical structures and the determination of their positions and orientations is presented in this paper. The scattered field, which is generated by the interaction of excited guided waves with a defect, is evaluated with an approach named time reverse numerical simulation method (TRNS). Since the excited waves and the scattered field propagate along the sample, the time-consuming scanning of the whole tube can be eliminated. The scattered displacement field is measured in three dimensions over time with a laser vibrometer at different locations distributed equally around the circumference at a fixed axial coordinate far away from the defect. Instead of analyzing the complicated time signals directly, they are played back in time. If the recorded displacement histories of the scattered field are reversed in time and played back in an identical structure, the waves travel back the same path and interfere to a maximum at their origin. The result is an amplitude increase at the position of the defect where the scattered field was generated. Instead of playing back the recorded time signals in an experiment, this step is replaced by a numerical simulation. Only this enables the visualization and detection of the amplitude increase. As long as the simulation is of high accuracy, the position of the maximum interference corresponds exactly to the location of the defect in the experiment, although no defect is implemented in the simulation.     

Calculation of the Electron Efficiency for a Two-Cavity Gyroklystron Amplifier by Using a Self-Consistent Code

The beam-wave interaction in a Ka-band, two-cavity fundamental gyroklystron amplifier is studied by using a self-consistent nonlinear simulation code. The electron efficiency for this gyroklystron amplifier is calculated, and the effect of various parameters, such as beam voltage, beam current, electron guiding center radius, velocity pitch ratio and drift tube length on the electron efficiency is discussed in detail.     

Simulation of the Field Emission Arrays with Parallel Wedge Emitter, Gates and Co-Planar Lens

Using the particle-in-cell (PIC) simulation code MAGIC, the wedge emitter, parallel gate and anode with and without co-planar lens have been simulated. The electron beam streamlines and the X-V_x phase space have been evaluated, which show that the co-planar lens focuses the electron. Simulation results are presented for a single FEA's cell with the anode biased at 500V and spaced 100m from the cathode, which indicated that the electron beam width at anode is about 60.2m for without lens FEAs, after add the lens with –70V, the electron beam width will be decrease to 8.1m, meanwhile the maximum transverse velocity will decrease from ±1×10⁶m/s to ±1.46×10⁵m/s.     

Generalization of the lineshape useful in magnetic resonance spectroscopy

A new lineshape function is derived from the Tsallis distribution to describe electron paramagnetic resonance (EPR) spectra, and possibly nuclear magnetic resonance (NMR) spectra as well. This lineshape generalizes the Gaussian and Lorentzian lineshapes that are widely used in simulations. The main features of this lineshape function are presented: the normalization, moments, and first derivative. A number of experimental EPR spectra are compared with the results of simulations employing the new lineshape function. The results show that the new lineshape often provides a better approximation of the experimental spectrum. It is also shown that the new parameter of the lineshape function can be used to quantify the intermolecular spin-spin interactions.     

High order accurate,one-sided finite-difference approximations to concentration gradients at the boundaries,for the simulation of electrochemical reaction-diffusion problems in one-dimensional space geometry

Accurate calculation of concentration gradients at the boundaries is crucial in electrochemical kinetic simulations, owing to the frequent occurrence of gradient-dependent boundary conditions, and the importance of the gradient-dependent electric current. By using the information about higher spatial derivatives of the concentrations, contained in the time-dependent, kinetic reaction-diffusion partial differential equation(s) in one-dimensional space geometry, under appropriate assumptions it is possible to increase the accuracy orders of the conventional, one-sided n-point finite-difference formulae for the concentration gradients at the boundaries, without increasing n. In this way a new class of high order accurate gradient approximations is derived, and tested in simulations of potential-step chronoamperometric and current-step chronopotentiometric transients for the Reinert-Berg system. The new formulae possess advantages over the conventional gradient approximations. For example, they allow one to obtain a third order accuracy by using two space points only, or fourth order accuracy by using three points, and yet they yield smaller errors than the conventional four-point, or five-point formulae, respectively. Needing fewer points, for approximating the gradients with a given accuracy, simplifies also the solution of the linear algebraic equations arising from the application of implicit time integration schemes.     

Weighted Discrepancy and High-Dimensional Numerical Integration

The concept of weighted discrepancy of sequences was introduced by Sloan and Woniakowski when they proved a general form of a Koksma–Hlawka inequality for the numerical integration of functions. This version takes imbalances in the importance of the projections of the integrand into account.In this paper we give estimates for the weighted discrepancy of several important point sets. Further we carry out various (high-dimensional) numerical integration experiments and we compare the results with the error bounds provided by the generalized Koksma–Hlawka inequality and by the estimates for the weighted discrepancy. Finally we discuss various consequences.     

Exposure Simulation Model for Chemically Amplified Resists

In this study, after Dill’s model is discussed for transmittance and refractive indices of the non-chemically amplified resists, G- and I-line novolak resists, and the chemically amplified resists, a modification of Dill’s model as a new exposure model is introduced. The simulation results obtained using this new model with the multi-thin film interface method and the Berning theory have shown a good matching to the experimental data. Also, the simulated transmittance change due to the exposure parameters are used to analyze the influence of the coefficients on the transmittance.     

Modelling of magnetic effects in near-field optics

Green's dyadic technique represents a powerful tool for calculations in electrodynamics, especially in modelling optical properties of nanoscopic objects. The method does not only provide field distributions, but also maps of susceptibilities and densities of states. Whereas the formalism is well established for dielectrics and electric fields, I present here a straight forward extension to tensors of both electric and magnetic type as well as mixed ones and furthermore to the situation where objects with dielectric and magnetic permeabilities are present together. As examples, characteristic field patterns are compared for elementary dielectric and magnetic perturbations. Green's tensors calculated for a coral structure reveal that mixed susceptibilities can exhibit other symmetries than pure electric or magnetic ones. Maps of all tensor components can thus give essential clues to the interpretation of near-field images. Received 15 December 2002 Published online 20 June 2003 RID="a" ID="a"Files “maths.ps” and “tensors.ps” are only available in electronic form at http://www.edpsciences.org RID="b" ID="b"e-mail: Ursula.Schroeter@uni-konstanz.de