A Cip/Multi-Moment Finite Volume Method For Shallow Water Equations With Source Terms

A novel finite volume method has been presented to solve the shallow water equations. In addition to the volume-integrated average (VIA) for each mesh cell, the surface-integrated average (SIA) is also treated as the model variable and is independently predicted. The numerical reconstruction is conducted based on both the VIA and the SIA. Different approaches are used to update VIA and SIA separately. The SIA is updated by a semi-Lagrangian scheme in terms of the Riemann invariants of the shallow water equations, while the VIA is computed by a flux-based finite volume formulation and is thus exactly conserved. Numerical oscillation can be effectively avoided through the use of a non-oscillatory interpolation function. The numerical formulations for both SIA and VIA moments maintain exactly the balance between the fluxes and the source terms. 1D and 2D numerical formulations are validated with numerical experiments. Copyright (c) 2007 John Wiley & Sons, Ltd.

A Modified Free Volume Model For Characterizing Of Rate Effect In Bulk Metallic Glasses

We investigate the plastic deformation and constitutive behaviour of bulk metallic glasses (BMGs). A dimensionless Deborah number De(ID) = t(r)/t(i) is proposed to characterize the rate effect in BMGs, where t(r) is the structural relaxing characteristic time of BMGs under shear load, t(i) is the macroscopic imposed characteristic time of applied stress or the characteristic time of macroscopic deformation. The results demonstrate that the modified free volume model can characterize the strain rate effect in BMGs effectively.

Comparative observation of Ar, Ar-H-2 and Ar-N-2 DC arc plasma jets and their arc root behaviour at reduced pressure

Results observed experimentally are presented, about the DC arc plasma jets and their arc-root behaviour generated at reduced gas pressure without or with an applied magnetic field. Pure argon, argon-hydrogen or argon- nitrogen mixture was used as the plasma-forming gas. A specially designed copper mirror was used for a better observation of the arc-root behaviour on the anode surface of the DC non-transferred arc plasma torch. It was found that in the cases without an applied magnetic field, the laminar plasma jets were stable and approximately axisymmetrical. The arc-root attachment on the anode surface was completely diffusive when argon was used as the plasma-forming gas, while the arc-root attachment often became constrictive when hydrogen or nitrogen was added into the argon. As an external magnetic field was applied, the arc root tended to rotate along the anode surface of the non-transferred arc plasma torch.

Influence of liquid bridge volume on the onset of oscillation in floating-zone convection .1. Experiments

The onset of oscillation in the floating zone convection driven by the gradient of surface tension was experimentally studied, and discussions were concentrated on the influence of liquid bridge volume on the onset of oscillation. Distributions of critical applied temperature difference and frequency depending on the volume of the liquid bridge were obtained, and there was a gap range of liquid volume which separated the curve of marginal stability into two parts for fixed rod diameter and aspect ratio. The results imply that the volume of the liquid bridge is a sensitive critical parameter for the onset of oscillation. The implication on the instability is also discussed in the present paper.

Influence of liquid bridge volume on the onset of oscillation in floating-zone convection .2. Numerical-simulation

The onset of oscillation in the floating zone convection driven by the gradient of surface tension was studied numerically for an unsteady and two-dimensional model, and studies were concentrated on the influence of liquid bridge volume on the onset of oscillation in comparison with the experimental results in the Paper I. The numerical results agree with the experimental ones presented in the previous paper, in which the distributions of critical applied temperature difference depending on the volume of liquid bridge and a gap range of liquid volume in marginal stability curve were obtained.

Comparative Observation of DC Arc Plasma Jets and Their Arc Root Behaviors at Reduced Pressure

Experimentally observed, results are presented for the DCarcplasmajets and theirarc-rootbehaviors generated atreduced gas pressure and without or with an' applied magnetic field. Pure argon, argon -hydrogen or argon-nitrogen mixture is used as the plasma-forming gas. A specially designed copper mirror is constructed and used for better observing the arc-root behavior on the anode surface of the DC non-transferred arcplasma torch. It is shown that for the cases without applied magnetic field, the laminar plasmajets are stable and approximately axisymmetrical. The arc-root attachment on the anode surface is completely diffusive when argon is used as the plasma-forming gas, while the arc-root attachment often becomes constrictive when hydrogen or nitrogen is added into the argon. When an external magnetic field is applied, the arcroot tends to rotate along the anode surface of the non-transferred arcplasma torch.

3D Finite Volume Scheme for Czochralski Crystal Growth

A general three-dimensional model is developed for simulation of the growth process of silicon single crystals by Czochralski technique. The numerical scheme is based on the curvilinear non-orthogonal finite volume discretization. Numerical solutions show that the flow and temperature fields in the melt are asymmetric and unsteady for 8’’ silicon growth. The effects of rotation of crystal on the flow structure are studied. The rotation of crystal forms the Ekman layer in which the temperature gradient along solid/melt surface is small.

3D Finite Volume Scheme for Czochralski Single Crystal Growth

Czochralski (Cz) technique, which is used for growing single crystals, has dominated the production of single crystals for electronic applications. The Cz growth process involves multiple phases, moving interface and three-dimensional behavior. Much has been done to study these phenomena by means of numerical methods as well as experimental observations. A three-dimensional curvilinear finite volume based algorithm has been developed to model the Cz process. A body-fitted transformation based approach is adopted in conjunction with a multizone adaptive grid generation (MAGG) technique to accurately handle the three-dimensional problems of phase-change in irregular geometries with free and moving surfaces. The multizone adaptive model is used to perform a three-dimensional simulation of the Cz growth of silicon single crystals.Since the phase change interface are irregular in shape and they move in response to the solution, accurate treatment of these interfaces is important from numerical accuracy point of view. The multizone adaptive grid generation (MAGG) is the appropriate scheme for this purpose. Another challenge encountered is the moving and periodic boundary conditions, which is essential to the numerical solution of the governing equations. Special treatments are implemented to impose the periodic boundary condition in a particular direction and to determine the internal boundary position and shape varying with the combination of ambient physicochemical transport process and interfacial dynamics. As indicated above that the applications and processes characterized by multi-phase, moving interfaces and irregular shape render the associated physical phenomena three-dimensional and unsteady. Therefore a generalized 3D model rather than a 2D simulation, in which the governing equations are solved in a general non-orthogonal coordinate system, is constructed to describe and capture the features of the growth process. All this has been implemented and validated by using it to model the low pressure Cz growth of silicon. Accuracy of this scheme is demonstrated by agreement of simulation data with available experimental data. Using the quasi-steady state approximation, it is shown that the flow and temperature fields in the melt under certain operating conditions become asymmetric and unsteady even in the absence of extrinsic sources of asymmetry. Asymmetry in the flow and temperature fields, caused by high shear initiated phenomena, affects the interface shape in the azimuthal direction thus results in the thermal stress distribution in the vicinity, which has serious implications from crystal quality point of view.

Perturbational finite volume scheme for the one-dimensional Navier-Stokes equations

Starting from the second-order finite volume scheme,though numerical value perturbation of the cell facial fluxes, the perturbational finite volume (PFV) scheme of the Navier-Stokes (NS) equations for compressible flow is developed in this paper. The central PFV scheme is used to compute the one-dimensional NS equations with shock wave.Numerical results show that the PFV scheme can obtain essentially non-oscillatory solution.

A Hybrid Scheme Based on Finite Element/Volume Methods for Two Immiscible Fluid Flows

We have successfully extended our implicit hybrid finite element/volume (FE/FV) solver to flows involving two immiscible fluids. The solver is based on the segregated pressure correction or projection method on staggered unstructured hybrid meshes. An intermediate velocity field is first obtained by solving the momentum equations with the matrix-free implicit cell-centered FV method. The pressure Poisson equation is solved by the node-based Galerkin FE method for an auxiliary variable. The auxiliary variable is used to update the velocity field and the pressure field. The pressure field is carefully updated by taking into account the velocity divergence field. This updating strategy can be rigorously proven to be able to eliminate the unphysical pressure boundary layer and is crucial for the correct temporal convergence rate. Our current staggered-mesh scheme is distinct from other conventional ones in that we store the velocity components at cell centers and the auxiliary variable at vertices. The fluid interface is captured by solving an advection equation for the volume fraction of one of the fluids. The same matrix-free FV method, as the one used for momentum equations, is used to solve the advection equation. We will focus on the interface sharpening strategy to minimize the smearing of the interface over time. We have developed and implemented a global mass conservation algorithm that enforces the conservation of the mass for each fluid.

High Order Multi-Moment Constrained Finite Volume Method. Part I: Basic Formulation

A new high-order finite volume method based on local reconstruction is presented in this paper. The method, so-called the multi-moment constrained finite volume (MCV) method, uses the point values defined within single cell at equally spaced points as the model variables (or unknowns). The time evolution equations used to update the unknowns are derived from a set of constraint conditions imposed on multi kinds of moments, i.e. the cell-averaged value and the point-wise value of the state variable and its derivatives. The finite volume constraint on the cell-average guarantees the numerical conservativeness of the method. Most constraint conditions are imposed on the cell boundaries, where the numerical flux and its derivatives are solved as general Riemann problems. A multi-moment constrained Lagrange interpolation reconstruction for the demanded order of accuracy is constructed over single cell and converts the evolution equations of the moments to those of the unknowns. The presented method provides a general framework to construct efficient schemes of high orders. The basic formulations for hyperbolic conservation laws in 1- and 2D structured grids are detailed with the numerical results of widely used benchmark tests. (C) 2009 Elsevier Inc. All rights reserved.

Development of arc root attachment in the nozzle of 1 kW N2 and H2-N2 arcjet thrusters

Arc root behavior affects the energy transfer and nozzle erosion in an arcjet thruster. To investigate the development of arc root attachment in 1 kW class N2 and H2-N2 arcjet thrusters from the time of ignition to the stably working condition, a kinetic series of end-on view images of the nozzle obtained by a high-speed video camera was analyzed. The addition of hydrogen leads to higher arc voltage levels and the determining factor for the mode of arc root attachment was found to be the nozzle temperature. At lower nozzle temperatures, constricted type attachment with unstable motions of the arc root was observed, while a fully diffused and stable arc root was observed at elevated nozzle temperatures.

COLOR-HOLOGRAPHY USING THE ANGULAR SELECTIVITY OF VOLUME RECORDING MEDIA

A display hologram of an object can be recorded and reconstructed in three primary colors if the angular selectivity of volume recording media is exploited. Three holograms are recorded in the same medium, each at a different primary color. These three holograms are reconstructed by simultaneous illumination of the hologram with the original reference beams. By proper choice of the angles that the reference beams make to the hologram, it is possible to suppress strongly cross talk between the different reconstructions (e.g., the red object reconstruction in green light). The technique exhibits high resolution, high diffraction efficiency, and vivid colors. Through the addition of three holographically recorded volume gratings it is possible to reconstruct the hologram with a beam of white light. The saturation and brightness of each primary color in the reconstruction can be adjusted by selection of an appropriate thickness for the corresponding grating.