The dynamical process of a coronal transient associated with an eruptive prominence .2. Analytical solutions in finite regions

piston

Elastic-plastic finite element analyses of short cracks

Stress and strain distributions and crack opening displacement characteristics of short cracks have been studied in single edge notch bend and centre cracked panel specimens using elastic–plastic finite element analyses incorporating both a non strain hardening and a power law hardening behaviour. J contour integral solutions to describe stress strain conditions at crack tips for short cracks differ from those for long cracks. The analyses show that (i) short cracks can propagate at stress levels lower than those required for long cracks and (ii) a two-parameter description of crack tip fields is necessary for crack propagation.

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.

A Finite Element Analysis of a Crack Penetrating or Deflecting into an Interface in a Thin Laminate

The crack tip driving force of a crack growing from a pre-crack that is perpendicular to and terminating at an interface between two materials is investigated using a linear fracture mechanics theory. The analysis is performed both for a crack penetrating the interface, growing straight ahead, and for a crack deflecting into the interface. The results from finite element calculations are compared with asymptotic solutions for infinitesimally small crack extensions. The solution is found to be accurate even for fairly large amounts of crack growth. Further, by comparing the crack tip driving force of the deflected crack with that of the penetrating crack, it is shown how to control the path of the crack by choosing the adhesion of the interface relative to the material toughness.

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.

Elastic deformation of soft membrane with finite thickness induced by a sessile liquid droplet

In this paper, the role of vertical component of Surface tension of a droplet on the elastic deformation of a finite-thickness flexible membrane was theoretically analyzed using Hankel transformation. The vertical displacement at the Surface was derived and can be reduced to Lester's or Rusanov's solutions when the thickness is infinite. Moreover, some Simulations of the effect of a liquid droplet on a membrane with a finite thickness were made. The numerical results showed that there exists a saturated membrane thickness of the order of millimeter, when the thickness of a membrane is larger than such a value, the membrane can be regarded as a half-infinite body. Further numerical calculations for soft membrane whose thickness is far below the saturated thickness were made. By comparison between the maximum vertical displacement of an ultrathin soft membrane and a half-infinite body, we found that Lester's or Rusanov's solutions for a half-infinite body cannot correctly describe Such cases. In other words, the thickness of a soft membrane has great effect on the surface deformation of the ultrathin membrane induced by a liquid droplet. (C) 2009 Elsevier Inc. All rights reserved.

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.

Analysis of an ultrashort pulsed finite beam diffracted by volume gratings

We obtain analytical solutions of the coupled wave equations that describe the Bragg diffraction of ultrashort pulsed finite beams by a thick planar grating, using two-dimensional coupled wave theory. The diffraction properties for the case of an ultrashort pulsed finite beam with Gaussian profiles in both the time and spatial domains are investigated. The spectral bandwidth of the diffracted beam, the Bragg selectivity bandwidth and the diffraction efficiency of the volume grating are influenced by the geometry parameter and the input bandwidth. Therefore extra attention should be paid to designing optical elements based on volume gratings for use with ultrashort pulsed waves in applications of pulse shaping and processing.

Three-dimensional coupled wave study for finite-sized anisotropic volume holographic gratings under ultrashort pulsed beam readout

The three-dimensional coupled wave theory is extended to systematically investigate the diffraction properties of finite-sized anisotropic volume holographic gratings (VHGs) under ultrashort pulsed beam (UPB) readout. The effects of the grating geometrical size and the polarizations of the recording and readout beams on the diffraction properties are presented, in particular under the influence of grating material dispersion. The wavelength selectivity of the finite-sized VHG is analyzed. The wavelength selectivity determines the intensity distributions of the transmitted and diffracted pulsed beams along the output face of the VHG. The distortion and widening of the diffracted pulsed beams are different for different points on the output face, as is numerically shown for a VHG recorded in a LiNbO3 crystal. The beam quality is analyzed, and the variations of the total diffraction efficiency are shown in relation to the geometrical size of the grating and the temporal width of the readout UPB. In addition, the diffraction properties of the finite-sized and one-dimensional VHG for pulsed and continuous-wave readout are compared. The study shows the potential application of VHGs in controlling spatial and temporal features of UPBs simultaneously. (C) 2007 Optical Society of America

Anisotropic Bragg diffraction of finite-sized volume holographic grating in photorefractive crystals

Anisotropic diffraction of uniform plane wave by finite-sized volume holographic grating in photorefractive crystals is considered. It is found that the anisotropic diffraction can take place when some special conditions are satisfied. The diffracted image is obtained in experiment for the anisotropic Bragg diffraction in Fe:LiNbO3 crystals. A coupled wave analysis is presented to study the properties of anisotropic diffraction. An analytical integral solution for the amplitudes of the diffracted beams is submitted. A trade off between high diffraction efficiency and the deterioration of reconstruction fidelity is analyzed. Numerical evaluations also show that the finite-sized anisotropic volume grating exhibits strong angular and wavelength selectivity. All the results are useful for the optimizing design of VHOE based on finite-sized volume grating structures. (c) 2006 Elsevier GmbH. All rights reserved.

Finite element analysis on a circular wedge prism with 400 mm diameter

The paper comprehensively analyzes the distortions of a circular wedge prism with 400 mm diameter in a scanner by method of optical-mechanical-thermal integrating analysis. The structure and intensity of the prism assembly is verified and checked, and the surface deformations of the prism under gravity load, as well as the thermo-elastic distortions of the prism, are analyzed in detail and evaluated, which is finally contrasted with the measured values of Zygo Mark interferometer. The results show: the maximal distortion of the prism assembly is 10 nm magnitude and the maximal stress is 0.441 Mpa, which has much tolerance to the precision requirement of structure and the admissible stress of material; the influence of heat effect on the surface deformations of prism is proved to be far greater than the influence of gravity load, so some strict temperature-controlled measures are to be considered when the scanner is used. (c) 2006 Elsevier GmbH. All rights reserved.

Finite element beam propagation method for analysis of plasmonic waveguide

The basic idea of the finite element beam propagation method (FE-BPM) is described. It is applied to calculate the fundamental mode of a channel plasmonic polariton (CPP) waveguide to confirm its validity. Both the field distribution and the effective index of the, fundamental mode are given by the method. The convergence speed shows the advantage and stability of this method. Then a plasmonic waveguide with a dielectric strip deposited on a metal substrate is investigated, and the group velocity is negative for the fundamental mode of this kind of waveguide. The numerical result shows that the power flow direction is reverse to that of phase velocity.