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.

Instability of the vertically forced surface wave in a circular cylindrical container

The nonlinear amplitude equation, which was derived by Jian Yongjun employing expansion of two-time scales in inviscid fluids in a vertically oscillating circular cylindrical vessel, is modified by introducing a damping term due to the viscous dissipation of this system. Instability of the surface wave is analysed and properties of the solutions of the modified equation are determined together with phase-plane trajectories. A necessary condition of forming a stable surface wave is obtained and unstable regions are illustrated. Research results show that the stable pattern of surface wave will not lose its stability to an infinitesimal disturbance.

Vertically forced surface wave in weakly viscous fluids bounded in a circular cylindrical vessel

Two-time scale perturbation expansions were developed in weakly viscous fluids to investigate surface wave motions by linearizing the Navier-Stokes equation in a circular cylindrical vessel which is subject to a vertical oscillation. The fluid field was divided into an outer potential flow region and an inner boundary layer region. A linear amplitude equation of slowly varying complex amplitude, which incorporates a damping term and external excitation, was derived for the weakly viscid fluids. The condition for the appearance of stable surface waves was obtained and the critical curve was determined. In addition, an analytical expression for the damping coefficient was determined and the relationship between damping and other related parameters (such as viscosity, forced amplitude, forced frequency and the depth of fluid, etc.) was presented. Finally, the influence both of the surface tension and the weak viscosity on the mode formation was described by comparing theoretical and experimental results. The results show that when the forcing frequency is low, the viscosity of the fluid is prominent for the mode selection. However, when the forcing frequency is high, the surface tension of the fluid is prominent.

Internal resonance of an L-shaped beam with a limit stop .2. Forced vibration

A limit stop is placed at the elbow of an L-shaped beam whose linear natural frequencies are nearly commensurable. As a result of this hardening device the non-linear system exhibits multiple internal resonances, which involve various degree of coupling between the first five modes of the beam in free vibration. A point load is so placed as to excite several modes and the resulting forced vibration is examined. In the undamped case, three in-phase and two out-of-phase solution branches have been found. The resonance curve is extremely complicated, with multiple branches and interactions between the first four modes. The amplitudes of the higher harmonics are highly influenced by damping, the presence of which can effectively attenuate internal resonances. Consequently parts of the resonance curve may be eliminated, with the resulting response comprising different distinctive branches. (C) 1996 Academic Press Limited

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.

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.

2D Kinetics and Forced Dissociation of Selectin-ligand Bindings

Cell adhesion is crucial to many pathophysiological processes, such as inflammatory reaction and tumor metastasis. It is mediated by specific interactions between receptors and ligands, and provides the physical linkages among cells. For example, interactions between selectins and glycoconjugate ligands mediate leukocyte initially tethering to and subsequently rolling on vascular surfaces in sites of inflammation or injury, which is determined by their fast kinetic rates. To mediate cell adhesion, the interacting receptors and ligands must anchor to apposing surfaces of two cells or a cell and the substratum, i.e. , the so-called two-dimensional (2D) binding, which differs from interactions in the fluid phase, i.e. , the three-dimensional (3D) binding. How structural variations and surface environments of interacting molecules affect their 2D kinetics, and how external forces manipulate their dissociation has little been known quantitatively, and nowadays attracts more and more attentions.

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.

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.