A Multiscale-Domain Decomposition Method for High Reynolds Number Flows

The high Reynolds number flow contains a wide range of length and time scales, and the flow domain can be divided into several sub-domains with different characteristic scales. In some sub-domains, the viscosity dissipation scale can only be considered in a certain direction; in some sub-domains, the viscosity dissipation scales need to be considered in all directions; in some sub-domains, the viscosity dissipation scales are unnecessary to be considered at all. For laminar boundary layer region, the characteristic length scales in the streamwise and normal directions are L and L Re-1/ 2 , respectively. The characteristic length scale and the velocity scale in the outer region of the boundary layer are L and U, respectively. In the neighborhood region of the separated point, the length scale l<method is developed for the high Reynolds number flow. First, the whole domain is decomposed to different sub-domains with the different characteristic scales. Then the different dominant equation of all sub-domains is defined according to the diffusion parabolized (DP) theory of viscous flow. Finally these different equations are solved simultaneously in whole computational region. For numerical tests of high Reynolds numerical flows, two-dimensional supersonic flows over rearward and frontward steps as well as an interaction flow between shock wave and boundary layer were solved numerically. The pressure distributions and local coefficients of skin friction on the wall are given. The numerical results obtained by the multiscale-domain decomposition algorithm are well agreement with those by NS equations. Comparing with the usual method of solving the Navier-Stokes equations in the whole flow, under the same numerical accuracy, the present multiscale domain decomposition method decreases CPU consuming about 20% and reflects the physical mechanism of practical flow more accurately.

Integration property of a wedge-elliptic cone waverider with scramjet engine

The hypersonic waverider forebody is designed in this paper. For the present waverider, the undersurface is carved out as a stream surface of a hypersonic inviscid flow field around wedge-elliptic cone, and the upper surface is assumed to be a freestream surface. A finite-volume code is used to generate the three-dimensional flow field. The leading edge is determined by satisfying the condition that the lip is situated at the intersection line of shocks.

A smoothing technique for discrete delta functions with application to immersed boundary method in moving boundary simulations.

The effects of complex boundary conditions on flows are represented by a volume force in the immersed boundary methods. The problem with this representation is that the volume force exhibits non-physical oscillations in moving boundary simulations. A smoothing technique for discrete delta functions has been developed in this paper to suppress the non-physical oscillations in the volume forces. We have found that the non-physical oscillations are mainly due to the fact that the derivatives of the regular discrete delta functions do not satisfy certain moment conditions. It has been shown that the smoothed discrete delta functions constructed in this paper have one-order higher derivative than the regular ones. Moreover, not only the smoothed discrete delta functions satisfy the first two discrete moment conditions, but also their derivatives satisfy one-order higher moment condition than the regular ones. The smoothed discrete delta functions are tested by three test cases: a one-dimensional heat equation with a moving singular force, a two-dimensional flow past an oscillating cylinder, and the vortex-induced vibration of a cylinder. The numerical examples in these cases demonstrate that the smoothed discrete delta functions can effectively suppress the non-physical oscillations in the volume forces and improve the accuracy of the immersed boundary method with direct forcing in moving boundary simulations.

An Instrumented Indentation Method for Young's Modulus Measurement with Accuracy Estimation

The relationships between indentation responses and Young's modulus of an indented material were investigated by employing dimensional analysis and finite element method. Three representative tip bluntness geometries were introduced to describe the shape of a real Berkovich indenter. It was demonstrated that for each of these bluntness geometries, a set of approximate indentation relationships correlating the ratio of nominal hardness/reduced Young's modulus H (n) /E (r) and the ratio of elastic work/total work W (e)/W can be derived. Consequently, a method for Young's modulus measurement combined with its accuracy estimation was established on basis of these relationships. The effectiveness of this approach was verified by performing nanoindentation tests on S45C carbon steel and 6061 aluminum alloy and microindentation tests on aluminum single crystal, GCr15 bearing steel and fused silica.

On optimal message vector length for block single parallel partition algorithm in a three-dimensional ADI solver

It has long been recognized that many direct parallel tridiagonal solvers are only efficient for solving a single tridiagonal equation of large sizes, and they become inefficient when naively used in a three-dimensional ADI solver. In order to improve the parallel efficiency of an ADI solver using a direct parallel solver, we implement the single parallel partition (SPP) algorithm in conjunction with message vectorization, which aggregates several communication messages into one to reduce the communication costs. The measured performances show that the longest allowable message vector length (MVL) is not necessarily the best choice. To understand this observation and optimize the performance, we propose an improved model that takes the cache effect into consideration. The optimal MVL for achieving the best performance is shown to depend on number of processors and grid sizes. Similar dependence of the optimal MVL is also found for the popular block pipelined method.

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.

MODIFIED ENGAGEMENT METHOD FOR MATRIX OPERATION

We describe a modified engagement method for matrix operation based on a two-dimensional crossed-ring interconnection network, Our method incorporates fewer steps than that reported by Bocker et al. [Appl. Opt. 22, 804 (1983)], and its performance is found to be the most efficient (minimum steps) in comparison with other systolic and/or engagement methods for matrix operation. Thus, it may be helpful for other optical and electronic implementations of matrix operations. One compact optoelectronic integrity approach for implementing the modified engagement method is briefly described. (C) 1995 Optical Society of America

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.

Absolute left-handed behaviors in a triangular elliptical-rod photonic crystal

In this paper we theoretically study the left-handed behaviors in a two-dimensional triangular photonic crystal made of elliptical rods in air. An absolute left-handed region is found in the second photonic band by using the plane wave expansion method to analyze the photonic band structure and equifrequency contours. Typical left-handed behaviors such as negative refraction, flat superlensing and plano-concave lensing are demonstrated by the finite-difference time-domain simulations. These behaviors are also compared with the quasi-negative refraction and the resulted focusing effects in a square-lattice two-dimensional photonic crystal. (c) 2005 Optical Society of America

Negative refraction in a honeycomb-lattice photonic crystal

In this paper we theoretically investigate a photonic crystal with dielectric rods in a honeycomb lattice. Two left-handed frequency regions are found in the second and third photonic band by using the plane wave expansion method to analyze the photonic band structure and equifrequency contours. Subwavelength imaging by the photonic crystal flat lens are systematically studied by numerical simulations using the multiple scattering method. Different from the photonic crystals with noncircular dielectric rods in air, this structure is almost isotropic at the optimal frequency for superlensing. As a comparison, flat slab focusing is also demonstrated at other frequencies in the two left-handed regions. (c) 2006 Elsevier Ltd. All rights reserved.

Influence of negative ion element impurities on laser induced damage threshold of HfO2 thin film

Negative ion element impurities breakdown model in HfO2 thin film was reported in this paper. The content of negative ion elements were detected by glow discharge mass spectrum analysis (GDMS); HfO2 thin films were deposited by the electron-beam evaporation method. The weak absorption and laser induced damage threshold (LIDT) of HfO2 thin films were measured to testify the negative ion element impurity breakdown model. It was found that the LIDT would decrease and the absorption would increase with increasing the content of negative ion element. These results indicated that negative ion elements were harmful impurities and would speed up the damage of thin film. (c) 2006 Elsevier B.V. All rights reserved.

Monitoring bioaccumulation and toxic effects of hexachlorobenzene using the polyurethane foam unit method in the microbial communities of the Fuhe River, Wuhan

Hexachlorobenzene (HCB) is a chlorinated aromatic hydrocarbon that was widely used for seed dressing in prevention of fungal growth on crops, and also as a component of fireworks, ammunition, and synthetic rubbers. Because of its resistance to degradation and mobility, HCB is widely distributed throughout the environment and is accumulated through food chains in different ecosystems. In this study, a preliminary investigation was carried out on the bioaccumulation and the toxic effects of HCB in the microbial (protozoan in particular) communities in the Fuhe River, Wuhan, a water body receiving industrial wastewaters containing HCB and other pollutants, using the standardized polyurethane foam units (PFU) method. Field samples were taken from eight stations established along the Fuhe River in January and August 2006. The concentration ratios of HCB in microbial communities and in water were 9.66-18.64, and the microbial communities accumulated 13.29-56.88 mu g/L of HCB in January and 0.82-10.25 mu g/L HCB in August. Correlation analysis showed a negative correlation between the HCB contents in the microbial assemblage, and the number of species and the diversity index of the protozoan communities. This study demonstrated the applicability of the PFU method in monitoring the effects of HCB on the level of microbial communities.

Improvement and observation of immunoelectron microscopic method for the localization of frog Rana grylio virus (RGV) in infected fish cells

In this paper, to understand the roles of amorphous structures which were observed within the viromatrix of Rana grylio virus (RGV), an improved immunoelectron microscopy (IEM) method was developed to detect the localization of RGV in carp Epithelipma papulosum cyprinid (EPC) cells. Infected EPC cells were fixed with 4% paraformaldehyde-0.25% glutaraldehyde mixture, dehydrated completely, and embedded in LR White resin. This method allowed good ultrastructural preservation and specific labeling with anti-RGV antibodies. The results of IEM showed that colloidal gold mainly bound to the capsids of viral particles at the stage of viral assembly, while during the viral maturation colloidal gold bound to the envelop of virions. In addition, within the viromatrix, the amorphous structures, including dense floccules, membranous materials and tubules, also had strong colloidal gold signals, revealing that those amorphous structures were participated in RGV assembly. In contrast, no significant gold labeling signals were obtained in negative controls. The present study not only provided further evidence that amorphous structures within the viromatrix were involved in the process of RGV assembly, but also developed an improved IEM method for studying the interaction between iridovirus and host cells. (C) 2006 Elsevier Ltd. All rights reserved.

Integration of double-fluorescence expression vectors into zebrafish genome for the selection of site-directed knockout/knockin

Production of zebrafish by modifying endogenous growth hormone (GH) gene through homologous recombination is described here. We first constructed the targeting vectors pGHT1.7k and pGHT2.8k, which were used for the knockout/knockin of the endogenous GH gene of zebrafish, and injected these two vectors into the embryos of zebrafish. Overall, the rate of targeted integration with the characteristic of germ line transmission in zebrafish was 1.7x10(-6). In one experimental patch, the integrating efficiency of pGHT2.8k was higher than that of pGHT1.7k, but the lethal effect of pGHT2.8k was stronger than that of pGHT1.7k. The clones with the correct integration of target genes were identified by a simple screening procedure based on green fluorescent protein (GFP) and RFP dual selection, which corresponded to homologous recombination and random insertion, respectively. The potential homologous recombination zebrafish was further bred to produce a heterozygous F-1 generation, selected based on the presence of GFP. The potential targeted integration of exogenous GH genes into a zebrafish genome at the P-0 generation was further verified by polymerase chain reaction and Southern blot analysis. Approximately 2.5% of potential founder knockout and knockin zebrafish had the characteristic of germ line transmission. In this study, we developed an efficient method for producing the targeted gene modification in zebrafish for future studies on genetic modifications and gene functions using this model organism.

Characterization of transgene integration pattern in F4 hGH-transgenic common carp (Cyprinus carpio L.)

The integration pattern and adjacent host sequences of the inserted pMThGH-transgene in the F4 hGH-transgenic common carp were extensively studied. Here we show that each F4 transgenic fish contained about 200 copies of the pMThGH-transgene and the transgenes were integrated into the host genome generally with concatemers in a head-to-tail arrangement at 4-5 insertion sites. By using a method of plasmid rescue, four hundred copies of transgenes from two individuals of F4 transgenic fish, A and B, were recovered and clarified into 6 classes. All classes of recovered transgenes contained either complete or partial pMThGH sequences. The class I, which comprised 83% and 84.5% respectively of the recovered transgene copies from fish A and B, had maintained the original configuration, indicating that most transgenes were faithfully inherited during the four generations of reproduction. The other five classes were different from the original configuration in both molecular weight and restriction map, indicating that a few transgenes had undergone mutation, rearrangement or deletion during integration and germline transmission. In the five types of aberrant transgenes, three flanking sequences of the host genome were analyzed. These sequences were common carp beta-actin gene, common carp DNA sequences homologous to mouse phosphoglycerate kinase-1 and human epidermal keratin 14, respectively.