Analysis of super compact finite difference method and application to simulation of vortex-shock interaction

Turbulence and aeroacoustic noise high-order accurate schemes are required, and preferred, for solving complex flow fields with multi-scale structures. In this paper a super compact finite difference method (SCFDM) is presented, the accuracy is analysed and the method is compared with a sixth-order traditional and compact finite difference approximation. The comparison shows that the sixth-order accurate super compact method has higher resolving efficiency. The sixth-order super compact method, with a three-stage Runge-Kutta method for approximation of the compressible Navier-Stokes equations, is used to solve the complex flow structures induced by vortex-shock interactions. The basic nature of the near-field sound generated by interaction is studied.

Comparison Of The Quasi-Static Method And The Dynamic Method For Simulating Fracture Processes In Concrete

Concrete is heterogeneous and usually described as a three-phase material, where matrix, aggregate and interface are distinguished. To take this heterogeneity into consideration, the Generalized Beam (GB) lattice model is adopted. The GB lattice model is much more computationally efficient than the beam lattice model. Numerical procedures of both quasi-static method and dynamic method are developed to simulate fracture processes in uniaxial tensile tests conducted on a concrete panel. Cases of different loading rates are compared with the quasi-static case. It is found that the inertia effect due to load increasing becomes less important and can be ignored with the loading rate decreasing, but the inertia effect due to unstable crack propagation remains considerable no matter how low the loading rate is. Therefore, an unrealistic result will be obtained if a fracture process including unstable cracking is simulated by the quasi-static procedure.

MILU-CG method and the numerical study on the flow around a rotating circular cylinder

A hybrid finite difference method and vortex method (HDV), which is based on domain decomposition and proposed by the authors (1992), is improved by using a modified incomplete LU decomposition conjugate gradient method (MILU-CG), and a high order implicit difference algorithm. The flow around a rotating circular cylinder at Reynolds number R-e = 1000, 200 and the angular to rectilinear speed ratio alpha is an element of (0.5, 3.25) is studied numerically. The long-time full developed features about the variations of the vortex patterns in the wake, and drag, lift forces on the cylinder are given. The calculated streamline contours agreed well with the experimental visualized flow pictures. The existence of critical states and the vortex patterns at the states are given for the first time. The maximum lift to drag force ratio can be obtained nearby the critical states.

Domain decomposition hybrid method for numerical simulation of bluff body flows:theoretical model and application

The discrete vortex method is not capable of precisely predicting the bluff body flow separation and the fine structure of flow field in the vicinity of the body surface. In order to make a theoretical improvement over the method and to reduce the difficulty in finite-difference solution of N-S equations at high Reynolds number, in the present paper, we suggest a new numerical simulation model and a theoretical method for domain decomposition hybrid combination of finite-difference method and vortex method. Specifically, the full flow. field is decomposed into two domains. In the region of O(R) near the body surface (R is the characteristic dimension of body), we use the finite-difference method to solve the N-S equations and in the exterior domain, we take the Lagrange-Euler vortex method. The connection and coupling conditions for flow in the two domains are established. The specific numerical scheme of this theoretical model is given. As a preliminary application, some numerical simulations for flows at Re=100 and Re-1000 about a circular cylinder are made, and compared with the finite-difference solution of N-S equations for full flow field and experimental results, and the stability of the solution against the change of the interface between the two domains is examined. The results show that the method of the present paper has the advantage of finite-difference solution for N-S equations in precisely predicting the fine structure of flow field, as well as the advantage of vortex method in efficiently computing the global characteristics of the separated flow. It saves computer time and reduces the amount of computation, as compared with pure N-S equation solution. The present method can be used for numerical simulation of bluff body flow at high Reynolds number and would exhibit even greater merit in that case.

Investigation of the Mechanical Properties of Thin Films by Nanoindentation, Considering the Effects of Thickness and Different Coating-Substrate Combinations

In order to further investigate nanoindentation data of film-substrate systems and to learn more about the mechanical properties of nanometer film-substrate systems, two kinds of films on different substrate systems have been tested with a systematic variation in film thickness and substrate characteristics. The two kinds of films are aluminum and tungsten, which have been sputtered on to glass and silicon substrates, respectively. Indentation experiments were performed with a Nano Indent XP II with indenter displacements typically about two times the nominal film thicknesses. The resulting data are analyzed in terms of load-displacement curves and various comparative parameters, such as hardness, Young's modulus, unloading stiffness and elastic recovery. Hardness and Young's modulus are investigated when the substrate effects are considered. The results show how the composite hardness and Young's modulus are different for different substrates, different films and different film thicknesses. An assumption of constant Young's modulus is used for the film-substrate system, in which the film and substrate have similar Young's moduli. Composite hardness obtained by the Joslin and Oliver method is compared with the directly measured hardness obtained by the Oliver and Pharr method.

Nanoindentation of thin-film-substrate system: Determination of film hardness and Young's modulus

In the present paper, the hardness and Young's modulus of film-substrate systems are determined by means of nanoindentation experiments and modified models. Aluminum film and two kinds of substrates; i.e. glass and silicon, are studied. Nanoindentation XP II and continuous stiffness mode are used during the experiments. In order to avoid the influence of the Oliver and Pharr method used in the experiments, the experiment data are analyzed with the constant Young's modulus assumption and the equal hardness assumption. The volume fraction model (CZ model) proposed by Fabes et al. (1992) is used and modified to analyze the measured hardness. The method proposed by Doerner and Nix (DN formula) (1986) is modified to analyze the measured Young's modulus. Two kinds of modified empirical formula are used to predict the present experiment results and those in the literature, which include the results of two kinds of systems, i.e., a soft film on a hard substrate and a hard film on a soft substrate. In the modified CZ model, the indentation influence angle, phi, is considered as a relevant physical parameter, which embodies the effects of the indenter tip radius, pile-up or sink-in phenomena and deformation of film and substrate.

A BIEM optimization method for fracture dynamics inverse problem

In the present paper, based on the theory of dynamic boundary integral equation, an optimization method for crack identification is set up in the Laplace frequency space, where the direct problem is solved by the author's new type boundary integral equations and a method for choosing the high sensitive frequency region is proposed. The results show that the method proposed is successful in using the information of boundary elastic wave and overcoming the ill-posed difficulties on solution, and helpful to improve the identification precision.

A new expression of hardening coefficients for fcc-crystal and calibration of the material constants

In order to describe the effect of latent hardening on the macro-plastic behavior of foc-crystal, a new expression for hardening coefficient is proposed in which there are 12 material constants, each having clear physical meaning. And a method of material constant calibration is suggested and used to determine the material constants of copper and aluminum crystal. The simulated load-elongation curves along various crystallographic orientations are comparable with the experimental ones.

Ablation and ultrafast dynamics of zinc selenide under femtosecond laser irradiation

The ablation in zinc selenide (ZnSe) crystal is studied by using 150-fs, 800-nm laser system. The images of the ablation pit measured by scanning electronic microscope (SEM) show no thermal stress and melting dynamics. The threshold fluence is measured to be 0.7 J/cm2. The ultrafast ablation dynamics is studied by using pump and probe method. The result suggests that optical breakdown and ultrafast melting take place in ZnSe irradiated under femtosecond laser pulses.

Ultraviolet-infrared femtosecond laser-induced damage in fused silica and CaF2 crystals

The damage in fused silica and CaF2 crystals induced by wavelength tunable femtosecond lasers is studied. The threshold fluence is observed to increase rapidly with laser wavelength lambda in the region of 250-800 nm, while it is nearly a constant for 800 and probe method. The reflectivity increases rapidly in the latter half of pump pulse, which supports that impact ionization plays an important role in the generation of conduction band electrons (CBEs). We study the CBEs absorption via subconduction-band (sub-CB) transition, and develop a coupled avalanche model. Our results indicate that the CBEs absorption via sub-CB transition plays an important role in the damage in dielectrics irradiated by the visible and near ultraviolet femtosecond lasers. Our theory explains well the experiments.

Optical two-dimensional SW-banyan network: Optical implementation, routing control, and determination of permissible permutations

A 2-D SW-banyan network is introduced by properly folding the 1-D SW-banyan network, and its corresponding optical setup is proposed by means of polarizing beamsplitters and 2-D phase spatial light modulators. Then, based on the characteristics and the proposed optical setup, the control for the routing path between any source-destination pair is given, and the method to determine whether a given permutation is permissible or not is discussed. Because the proposed optical setup consists of only optical polarization elements, it is compact in structure, its corresponding energy loss and crosstalk are low, and its corresponding available number of channels is high. (C) 1996 Society of Photo-Optical Instrumentation Engineers.

Study of polymer thermal and mechanical properties using a fiber Bragg grating thermal sensor

A simple and practical method for the study of polymer thermal and mechanical properties using a fiber Bragg grating (FBG) sensor is presented for the first time, in which the FBG is embedded in a typical epoxy polymer. By measuring the sensitivity change of the FBG sensor, changes of the thermal-mechanical properties of the polymer with temperature and pressure can be measured. The experimental results show that this technique is capable of providing continuous in-line monitoring such properties with high sensitivity during transformation between the glassy state and the rubbery state of a polymer within the temperature and pressure range of 20 to 180 C and 0 to 15 MPa. (c) 2007 Society of Photo-Optical Instrumentation Engineers.

Study of polymer thermal and mechanical properties using a fiber Bragg grating thermal sensor

A simple and practical method for the study of polymer thermal and mechanical properties using a fiber Bragg grating (FBG) sensor is presented for the first time, in which the FBG is embedded in a typical epoxy polymer. By measuring the sensitivity change of the FBG sensor, changes of the thermal-mechanical properties of the polymer with temperature and pressure can be measured. The experimental results show that this technique is capable of providing continuous in-line monitoring such properties with high sensitivity during transformation between the glassy state and the rubbery state of a polymer within the temperature and pressure range of 20 to 180 C and 0 to 15 MPa. (c) 2007 Society of Photo-Optical Instrumentation Engineers.

Numerical analysis and comparison of three metal-oxide-type super-resolution near field structures

Based on the Fresnel-Kirchkoff diffraction theory, we build up a Gaussian diffraction model of metal-oxide-type super-resolution near field structure (super-RENS), which can describe far field optical properties. The spectral contrast induced by refractive index and the structural changes in AgOx, PtOx and PdOx thin films, which are the key functional layers in super-RENS, are studied by using this model. Comparison results indicate that the spectral contrast intensively on laser-induced distribution and change of the refractive index in the metal-oxide films. The readout mechanism of the metal-oxide-type super-RENS optical disc is further clarified. This Gaussian diffraction model can be used as a simple and effective method for choosing proper active materials in super-RENS.

Optical transitions and upconversion luminescence of Er3+/Yb3+-codoped halide modified tellurite glasses

Er3+-doped halide modified tellurite glasses were synthesized by conventional melting and quenching method. The Judd-Ofelt analysis was performed on the absorption spectra and the transition probabilities, excited state lifetimes, and the branching ratios were calculated and discussed. The intense infrared and visible fluorescence spectra under 980 nm excitation were obtained. Strong upconversion signal was observed at pumping power as low as 30 mW in the glasses with halide ions. The upconversion mechanisms and power dependent intensities were discussed, which showed two-photon process are involved for the green and red emissions. The decay times of the emitting states and the corresponding quantum efficiency were determined and explained. (C) 2004 American Institute of Physics.