A diagonally inverted LU implicit multigrid scheme for the 3-D Navier-Stokes equations and a two equation model of turbulence

An LU implicit multigrid algorithm is developed to calculate 3-D compressible viscous flows. This scheme solves the full 3-D Reynolds-Averaged Navier-Stokes equation with a two-equation kappa-epsilon model of turbulence. The flow equations are integrated by an efficient, diagonally inverted, LU implicit multigrid scheme while the kappa-epsilon equations are solved, uncoupled from the flow equations, by a block LU implicit algorithm. The flow equations are solved within the framework of the multigrid method using a four-grid level W-cycle, while the kappa-epsilon equations are iterated only on the finest grid. This treatment of the Reynolds-Averaged Navier-Stokes equations proves to be an efficient method for calculating 3-D compressible viscous flows

Numerical analysis of three-dimensional viscous internal flows

A 3-D Navier-Stokes code has been developed for analysis of turbomachinery blade rows and other internal flows. The Navier-Stokes equations are written in a Cartesian coordinate system rotating about the x-axis, and then mapped to a general body-fitted coordinate system. Streamwise viscous terms are neglected using the thin-layer assumption, and turbulence effects are modeled using the Baldwin-Lomax turbulence model. The equations are discretized using finite differences on stacked C-type grids and are solved using a multistage Runge-Kutta algorithm with a spatially-varying time step and implicit residual smoothing. Calculations have been made of a horseshoe vortex formed in front of a flat plate with a round leading edge standing in a turbulent endwall boundary layer. Comparisons are made with experimental data taken by Eckerle and Langston for a circular cylinder under similar conditions. Computer and measured results are compared in terms of endwall flow visualization pictures and total pressure loss contours and vector plots on the symmetry plane. Calculated details of the primary vortex show excellent agreement with the experimental data. The calculations also show a small secondary vortex not seen experimentally

A nonoscillatory, characteristically convected, finite volume scheme for multidimensional convection problems

A new, nonoscillatory upwind scheme is developed for the multidimensional convection equation. The scheme consists of an upwind, nonoscillatory interpolation of data to the surfaces of an intermediate finite volume; a characteristic convection of surface data to a midpoint time level; and a conservative time integration based on the midpoint rule. This procedure results in a convection scheme capable of resolving discontinuities neither aligned with, nor convected along, grid lines

A diagonally inverted LU implicit multigrid scheme

A new Diagonally Inverted LU Implicit scheme is developed within the framework of the multigrid method for the 3-D unsteady Euler equations. The matrix systems that are to be inverted in the LU scheme are treated by local diagonalizing transformations that decouple them into systems of scalar equations. Unlike the Diagonalized ADI method, the time accuracy of the LU scheme is not reduced since the diagonalization procedure does not destroy time conservation. Even more importantly, this diagonalization significantly reduces the computational effort required to solve the LU approximation and therefore transforms it into a more efficient method of numerically solving the 3-D Euler equations

Unsteady Shear Disturbances Within a Two Dimensional Stratified Flow

The origin and evolution of shear disturbances within a stratified, inviscid, incompressible flow are investigated numerically by a Clebsch/Weber decomposition based scheme. In contrast to homogeneous flows, within which vorticity can be redistributed but not generated, the presence of a density stratification can render an otherwise irrotational flow vortical. In this work, a kinematic decomposition of the unsteady Euler equations separates the unsteady velocity field into rotational and irrotational components. The subsequent evolution of these components is used to study the influence various velocity disturbances have on both stratified and homogeneous flows. In particular, the flow within a two-dimensional channel is used to investigate the evolution of rotational disturbances, generated or convected, downstream from an unsteady inflow condition. Contrasting simulations of both stratified and homogeneous flows are used to distinguish between redistributed inflow vorticity and that which is generated by a density stratification

Status of ECR ion sources at JAERI

At the Takasaki site of Japan Atomic Energy Research Institute, four ECR ion sources were purchased or developed so far. This paper will report their performance, modification and status. The outlines for each source are as follows; 1. OCTOPUS purchased from IBA s.a. has been in use with a cyclotron since 1990. The gas feed system was modified to change gas species within 10 minutes to avoid impurity ions in the cocktail beam acceleration technique of the cyclotron. 2. ECR-18 with 18-GHz microwave has a solenoid coil between a pair of mirror coils to change mirror ratio in a wide range. A bump between mirror peaks in the original axial field distribution was removed by halving the solenoid length. The performance in generating high charge state ions was significantly improved as a result. 3. HYPERNANOGAN was purchased from PANTECHNIK s.a. and installed in the cyclotron system this year. Test operation was successfully made with generation of Ar, Pb and Ta ions. 4. MINI ECR is a full permanent magnet source with 10 GHz microwave. As a result of modifications for vacuum and the application of a bias probe, 40 microA of Ar4+ beam was obtained, for example. Using this source, a new technique was developed to vaporize very refractory materials into plasma by fluorinating them with SF6 plasma. The source is in use for a 400 kV ion implantor since last year

Estimation of Spin-Spin Interaction by Weak Measurement Scheme

Precisely knowing an interaction Hamiltonian is crucial to realize quantum information tasks, especially to experimentally demonstrate a quantum computer and a quantum memory. We propose a scheme to experimentally evaluate the spin-spin interaction for a two-qubit system by the weak measurement technique initiated by Yakir Aharonov and his colleagues. Furthermore, we numerically confirm our proposed scheme in a specific system of a nitrogen vacancy center in diamond. This means that the weak measurement can also be taken as a concrete example of the quantum process tomography.Comment: 4 pages, 1 table, 2 figures, to appear in Europhysics Letter

Evolution of superconductivity in isovalent Te-substituted KxFe2-ySe2 crystals

We report the evolution of superconductivity and the phase diagram of the KxFe2-ySe2-zTez (z=0-0.6) crystals grown by a simple one-step synthesis. No structural transition is observed in any crystals, while lattice parameters exhibit a systematic expansion with Te content. The Tc exhibits a gradual decrease with increasing Te content from Tconset = 32.9 K at z = 0 to Tconset = 27.9 K at z = 0.5, followed by a sudden suppression of superconductivity at z = 0.6. Upon approaching a Te concentration of 0.6, the shielding volume fraction decreases and eventually drops to zero. Simultaneously, hump positions in r-T curve shift to lower temperatures. These results suggest that isovalent substitution of Te for Se in KxFe2-ySe2 crystals suppresses the superconductivity in this system.Comment: 10 pages, 1 table, 8 figure

X-ray Temperature and Mass Measurements to the Virial Radius of Abell 1413 with Suzaku

We present X-ray observations of the northern outskirts of the relaxed galaxy cluster A1413 with Suzaku, whose XIS instrument has the low intrinsic background needed to make measurements of these low surface brightness regions. We excise 15 point sources superimposed on the image above a flux of $1\times 10^{-14}$ \fluxunit (2--10keV) using XMM-Newton and Suzaku images of the cluster. We quantify all known systematic errors as part of our analysis, and show our statistical errors encompasses them for the most part. Our results extend previous measurements with Chandra and XMM-Newton, and show a significant temperature drop to about 3keV at the virial radius, $r_{200}$. Our entropy profile in the outer region ($> 0.5 r_{200}$) joins smoothly onto that of XMM-Newton, and shows a flatter slope compared with simple models, similar to a few other clusters observed at the virial radius. The integrated mass of the cluster at the virial radius is approximately $7.5\times10^{14}M_{\odot}$ and varies by about 30% depending on the particular method used to measure it.Comment: 32pages, 9 figures, accepted for publication in PAS

Irradiation-induced Ag nanocluster nucleation in silicate glasses: analogy with photography

The synthesis of Ag nanoclusters in sodalime silicate glasses and silica was studied by optical absorption (OA) and electron spin resonance (ESR) experiments under both low (gamma-ray) and high (MeV ion) deposited energy density irradiation conditions. Both types of irradiation create electrons and holes whose density and thermal evolution - notably via their interaction with defects - are shown to determine the clustering and growth rates of Ag nanocrystals. We thus establish the influence of redox interactions of defects and silver (poly)ions. The mechanisms are similar to the latent image formation in photography: irradiation-induced photoelectrons are trapped within the glass matrix, notably on dissolved noble metal ions and defects, which are thus neutralized (reverse oxidation reactions are also shown to exist). Annealing promotes metal atom diffusion, which in turn leads to cluster nuclei formation. The cluster density depends not only on the irradiation fluence, but also - and primarily - on the density of deposited energy and the redox properties of the glass. Ion irradiation (i.e., large deposited energy density) is far more effective in cluster formation, despite its lower neutralization efficiency (from Ag+ to Ag0) as compared to gamma photon irradiation.Comment: 48 pages, 18 figures, revised version publ. in Phys. Rev. B, pdf fil