Numerical-Value Perturbation Method with Application of Solving Convective-Diffusion Integral Equation

The convective--diffusion equation is of primary importance in such fields as fluid dynamics and heat transfer hi the numerical methods solving the convective-diffusion equation, the finite volume method can use conveniently diversified grids (structured and unstructured grids) and is suitable for very complex geometry The disadvantage of FV methods compared to the finite difference method is that FV-methods of order higher than second are more difficult to develop in three-dimensional cases. The second-order central scheme (2cs) offers a good compromise among accuracy, simplicity and efficiency, however, it will produce oscillatory solutions when the grid Reynolds numbers are large and then very fine grids are required to obtain accurate solution. The simplest first-order upwind (IUW) scheme satisfies the convective boundedness criteria, however. Its numerical diffusion is large. The power-law scheme, QMCK and second-order upwind (2UW) schemes are also often used in some commercial codes. Their numerical accurate are roughly consistent with that of ZCS. Therefore, it is meaningful to offer higher-accurate three point FV scheme. In this paper, the numerical-value perturbational method suggested by Zhi Gao is used to develop an upwind and mixed FV scheme using any higher-order interpolation and second-order integration approximations, which is called perturbational finite volume (PFV) scheme. The PFV scheme uses the least nodes similar to the standard three-point schemes, namely, the number of the nodes needed equals to unity plus the face-number of the control volume. For instanc6, in the two-dimensional (2-D) case, only four nodes for the triangle grids and five nodes for the Cartesian grids are utilized, respectively. The PFV scheme is applied on a number of 1-D problems, 2~Dand 3-D flow model equations. Comparing with other standard three-point schemes, The PFV scheme has much smaller numerical diffusion than the first-order upwind (IUW) scheme, its numerical accuracy are also higher than the second-order central scheme (2CS), the power-law scheme (PLS), the QUICK scheme and the second-order upwind(ZUW) scheme.

Numerical Analysis of the Effects of Free-Stream Thermo-Chemical State on the Test Model Flow Field in High-Enthalpy Tunnel

The effects of the free-stream thermo-chemical state on the test model flow field in the high-enthalpy tunnel are studied numerically. The properties of the free-stream, which is in thermo-chemical non-equilibrium, are determined by calculating the nozzle flow field. A free-stream with total enthalpy equal to the real one in the tunnel while in thermo-chemical equilibrium is constructed artificially to simulate the natural atmosphere condition. The flow fields over the test models (blunt cone and Apollo command capsule model) under both the non-equilibrium and the virtual equilibrium free-stream conditions are calculated. By comparing the properties including pressure, temperature, species concentration and radiation distributions of these two types of flow fields, the effects of the non-equilibrium state of the free-stream in the high-enthalpy shock tunnel are analyzed.

Numerical simulation of a single fluid particle rising in non-Newtonian fluids

In this work, a level set method is developed for simulating the motion of a fluid particle rising in non-Newtonian fluids described by generalized Newtonian as well as viscoelastic model fluids. As the shear-thinning model we use a Carreau-Yasuda model, and the viscoelastic effect can be modeled with Oldroyd-B constitutive equations. The control volume formulation with the SIMPLEC algorithm incorporated is used to solve the governing equations on a staggered Eulerian grid. The level set method is implemented to compute the motion of a bubble in a Newtonian fluid as one of typical examples for validation, and the computational results are in good agreement with the reported experimental data．The level set method is also applied for simulating a Newtonian drop rising in Carreau-Yasuda and Oldroyd-B fluids．Numerical results including noticeably negative wake behind the drop and viscosity field are obtained, and compare satisfactorily with the known literature data.

Numerical Investigation On The Flowfield Of "Swallowtail" Cavity For Supersonic Mixing Enhancement

A "swallowtail" cavity for the supersonic combustor was proposed to serve as an efficient flame holder for scramjets by enhancing the mass exchange between the cavity and the main flow. A numerical study on the "swallowtail" cavity was conducted by solving the three-dimensional Reynolds-averaged Navier-Stokes equations implemented with a k-epsilon turbulence model in a multi-block mesh. Turbulence model and numerical algorithms were validated first, and then test cases were calculated to investigate into the mechanism of cavity flows. Numerical results demonstrated that the certain mass in the supersonic main flow was sucked into the cavity and moved spirally toward the combustor walls. After that, the flow went out of the cavity at its lateral end, and finally was efficiently mixed with the main flow. The comparison between the "swallowtail" cavity and the conventional one showed that the mass exchanged between the cavity and the main flow was enhanced by the lateral flow that was induced due to the pressure gradient inside the cavity and was driven by the three-dimensional vortex ring generated from the "swallowtail" cavity structure.

Direct numerical simulation of hypersonic boundary layer transition over a blunt cone with a small angle of attack

The direct numerical simulation of boundary layer transition over a 5° half-cone-angle blunt cone is performed. The free-stream Mach number is 6 and the angle of attack is 1°. Random wall blow-and-suction perturbations are used to trigger the transition. Different from the authors’ previous work [Li et al., AIAA J. 46, 2899(2008)], the whole boundary layer flow over the cone is simulated (while in the author’s previous work, only two 45° regions around the leeward and the windward sections are simulated). The transition location on the cone surface is determined through the rapid increase in skin fraction coefficient (Cf). The transition line on the cone surface shows a nonmonotonic curve and the transition is delayed in the range of 0° ≤ θ ≤ 30° (θ = 0° is the leeward section). The mechanism of the delayed transition is studied by using joint frequency spectrum analysis and linear stability theory (LST). It is shown that the growth rates of unstable waves of the second mode are suppressed in the range of 20° ≤ θ ≤ 30°, which leads to the delayed transition location. Very low frequency waves VLFWs� are found in the time series recorded just before the transition location, and the periodic times of VLFWs are about one order larger than those of ordinary Mack second mode waves. Band-pass filter is used to analyze the low frequency waves, and they are deemed as the effect of large scale nonlinear perturbations triggered by LST waves when they are strong enough.The direct numerical simulation of boundary layer transition over a 5° half-cone-angle blunt cone is performed. The free-stream Mach number is 6 and the angle of attack is 1°. Random wall blow-and-suction perturbations are used to trigger the transition. Different from the authors’ previous work [ Li et al., AIAA J. 46, 2899 (2008) ], the whole boundary layer flow over the cone is simulated (while in the author’s previous work, only two 45° regions around the leeward and the windward sections are simulated). The transition location on the cone surface is determined through the rapid increase in skin fraction coefficient (Cf). The transition line on the cone surface shows a nonmonotonic curve and the transition is delayed in the range of 20° ≤ θ ≤ 30° (θ = 0° is the leeward section). The mechanism of the delayed transition is studied by using joint frequency spectrum analysis and linear stability theory (LST). It is shown that the growth rates of unstable waves of the second mode are suppressed in the range of 20° ≤ θ ≤ 30°, which leads to the delayed transition location. Very low frequency waves (VLFWs) are found in the time series recorded just before the transition location, and the periodic times of VLFWs are about one order larger than those of ordinary Mack second mode waves. Band-pass filter is used to analyze the low frequency waves, and they are deemed as the effect of large scale nonlinear perturbations triggered by LST waves when they are strong enough.

A scaling law of photoelectron angular distributions in one-cycle laser pulses

The photoelectron angular distributions (PADs) from above-threshold ionization of atoms irradiated by one-cycle laser pulses satisfy a scaling law. The scaling law denotes that the main features of the PADs are determined by four dimensionless parameters: (1) the ponderomotive number u(p) = U-p/hw, the ponderomotive energy U-p in units of laser photon energy; (2) the binding number E-b = E-b/h(w), the atomic binding energy E-b in units of laser photon energy; (3) the number of absorbed photons q; (4) the carrier-envelope phase phi(0), the phase of the carrier wave with respect to the envelope. We verify the scaling law by theoretical analysis and numerical calculation, compared to that in long-pulse case. A possible experimental test to verify the scaling law is suggested.

NUMERICAL STUDY OF DAMMANN ARRAY ILLUMINATORS

The numerical solutions of binary-phase (0, tau) gratings for one-dimensional array illuminators up to 32 are presented. Some fabrication errors, which are due to position-quantization errors, phase errors, dilation (or erosion) errors, and the side-slope error, are calculated and show that even-number array illuminators are superior to odd-number array illuminators when these fabrication errors are considered. One (0, tau) binary-phase, 8 x 16 array illuminator made with the wet-chemical-etching method is given in this paper.

D77, one benzoic acid derivative, functions as a novel anti-HIV-1 inhibitor targeting the interaction between integrase and cellular LEDGF/p75

Integration of viral-DNA into host chromosome mediated by the viral protein HIV-1 integrase (IN) is an essential step in the HIV-1 life cycle. In this process, Lens epithelium-derived growth factor (LEDGF/p75) is discovered to function as a cellular co-fa

Site-directed gene integration in transgenic zebrafish mediated by cre recombinase using a combination of mutant Lox sites

With current gene-transfer techniques in fish, insertion of DNA into the genome occurs randomly and in many instances at multiple sites. Associated position effects, copy number differences, and multiple gene interactions make gene expression experiments difficult to interpret and fish phenotype less predictable. To meet different fish engineering needs, we describe here a gene targeting model in zebrafish. At first, four target zebrafish lines, each harboring a single genomic lox71 target site, were generated by zebrafish transgenesis. The zygotes of transgenic zebrafish lines were coinjected with capped Cre mRNA and a knockin vector pZklox66RFP. Site-specific integration event happened from one target zebrafish line. In this line two integrant zebrafish were obtained from more than 80,000 targeted embryos (integrating efficiency about 10(-4) to 10(-5)) and confirmed to have a sole copy of the integrating DNA at the target genome site. Genomic polymerase chain reaction analysis and DNA sequencing verified the correct gene target events where lox71 and lox66 have accurately recombined into double mutant lox72 and wild-type loxP. Each integrant zebrafish chosen for analysis harbored the transgene rfp at the designated egfp concatenates. Although the Cre-mediated recombination is site specific, it is dependent on a randomly placed target site. That is, a genomic target cannot be preselected for integration based solely on its sequence. Conclusively, an rfp reporter gene was successfully inserted into the egfp target locus of zebrafish genome by Cre-lox-mediated recombination. This site-directed knockin system using the lox71/lox66 combination should be a promising gene-targeting platform serving various purposes in fish genetic engineering.

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.

Numerical study of sub-wavelength plasmonic waveguide

The authors present an analysis of a plasmonic waveguide, simulated using a two-dimensional finite-difference time-domain technique. With the surface structures located on the surface of the metal, the device is able to confine and guide light waves in a sub-wavelength scale. And two waveguides can be placed within 150 nm (similar to 6% of the incident wavelength) that will helpful for the optoelectronic integration. Within the 20 mu m simulation region, it is found that the intensity of the guided light at the interface is roughly two to four times the peak intensity of the incident light, and the propagation length can reach approximately 40 Pm at the wavelength of 2.44 mu m. (c) 2007 Elsevier B.V. All rights reserved.

One-dimensional quantum waveguide theory of Rashba electrons

The ballistic spin transport in one-dimensional waveguides with the Rashba effect is studied. Due to the Rashba effect, there are two electron states with different wave vectors for the same energy. The wave functions of two Rashba electron states are derived, and it is found that their phase depend on the direction of the circuit and the spin directions of two states are perpendicular to the circuit, with the +pi/2 and -pi/2 angles, respectively. The boundary conditions of the wave functions and their derivatives at the intersection of circuits are given, which can be used to investigate the waveguide transport properties of Rashba spin electron in circuits of any shape and structure. The eigenstates of the closed circular and square loops are studied by using the transfer matrix method. The transfer matrix M(E) of a circular arc is obtained by dividing the circular arc into N segments and multiplying the transfer matrix of each straight segment. The energies of eigenstates in the closed loop are obtained by solving the equation det[M(E)-I]=0. For the circular ring, the eigenenergies obtained with this method are in agreement with those obtained by solving the Schrodinger equation. For the square loop, the analytic formula of the eigenenergies is obtained first The transport properties of the AB ring and AB square loop and double square loop are studied using the boundary conditions and the transfer matrix method In the case of no magnetic field, the zero points of the reflection coefficients are just the energies of eigenstates in closed loops. In the case of magnetic field, the transmission and reflection coefficients all oscillate with the magnetic field; the oscillating period is Phi(m)=hc/e, independent of the shape of the loop, and Phi(m) is the magnetic flux through the loop. For the double loop the oscillating period is Phi(m)=hc/2e, in agreement with the experimental result. At last, we compared our method with Koga's experiment. (C) 2009 American Institute of Physics. [doi: 10.1063/1.3253752]

Analysis and numerical simulation of eigenmode characteristics for semiconductor lasers with an equilateral triangle micro-resonator

The eigenmodes confined in the equilateral triangle resonator (ETR) are analyzed by deriving the eigenvalues and the mode field distributions and by the finite difference time domain (FDTD) technique. The analytical results show that the one-period-length for the mode light rays inside the ETR is the perimeter of the ETR, and the number of transverse modes is limited by the condition of total internal reflection. In addition, the sum of the longitudinal mode index and the transverse mode index should be an even number, which limits the number of confined modes again. Based on the FDTD technique and the Pade approximation, we calculate the mode resonant frequencies and the quality factors from the local maximum and the width of the spectral distribution of the intensity The numerical results of mode frequencies agree very well with the analytical results, and the quality factor of the fundamental mode is usually higher than that of the higher order transverse modes. The results show that the ETR is suitable to realize single-made operation as semiconductor microcavity lasers.

Numerical simulations of interrupted and conventional microchannel heat sinks

We provide three-dimensional numerical simulations of conjugate heat transfer in conventional and the newly proposed interrupted microchannel heat sinks. The new microchannel heat sink consists of a set of separated zones adjoining shortened parallel microchannels and transverse microchambers. Multi-channel effect, physical property variations, and axial thermal conduction are considered. It is found that flow rate variations in different channels can be neglected, while heat received by different channels accounts for 2% deviations from the averaged value when the heat flux at the back surface of the silicon chip reaches 100 W/cm(2). The computed hydraulic and thermal boundary layers are redeveloping in each separated zone due to shortened flow length for the interrupted microchannel heat sink. The periodic thermal developing flow is responsible for the significant heat transfer enhancement. Two effects influence pressure drops across the newly proposed microchannel heat sink. The first one is the pressure recovery effect in the microchamber, while the second one is the head loss when liquid leaves the microchamber and enters the next zone. The first effect compensates or suppresses the second one, leading to similar or decreased pressure drop than that for the conventional microchannel heat sink, with the fluid Prandtl number larger than unity.

OPTICAL-PHONON MODES AND FROHLICH POTENTIAL IN ONE-DIMENSIONAL QUANTUM-WELL WIRES

By extending our microscopic model on optical-phonon modes in quantum wells to one-dimensional (1D) quantum-well wires (QWW), the optical displacements and associated electrostatic potentials of optical-phonon modes in 1D QWW are calculated. The modes can be clearly divided into confined LO bulklike, TO bulklike modes, and extended interfacelike modes provided the bulk phonon dispersion is ignored. The character of each type of mode is illustrated with special attention to the interfacelike modes, which are hybrids of longitudinal- and transverse-optical waves from the corresponding bulk materials. Based on the numerical results, approximate analytical formulas for bulklike modes are presented. As in 2D wells, both the optical displacements and Frohlich potentials for the bulklike modes vanish at the interfaces. The finite dispersion of bulk phonons has a more pronounced effect on the 1D phonon modes because interfacelike modes show mixed characteristics of 2D interface and bulklike modes.