Smoothed Particle Hydrodynamics (SPH): an Overview and Recent Developments

Smoothed particle hydrodynamics (SPH) is a meshfree particle method based on Lagrangian formulation, and has been widely applied to different areas in engineering and science. This paper presents an overview on the SPH method and its recent developments, including (1) the need for meshfree particle methods, and advantages of SPH, (2) approximation schemes of the conventional SPH method and numerical techniques for deriving SPH formulations for partial differential equations such as the Navier-Stokes (N-S) equations, (3) the role of the smoothing kernel functions and a general approach to construct smoothing kernel functions, (4) kernel and particle consistency for the SPH method, and approaches for restoring particle consistency, (5) several important numerical aspects, and (6) some recent applications of SPH. The paper ends with some concluding remarks.

Parallel and Convergent Evolution of the Dim-Light Vision Gene RH1 in Bats (Order: Chiroptera)

Rhodopsin, encoded by the gene Rhodopsin (RH1), is extremely sensitive to light, and is responsible for dim-light vision. Bats are nocturnal mammals that inhabit poor light environments. Megabats (Old-World fruit bats) generally have well-developed eyes, while microbats (insectivorous bats) have developed echolocation and in general their eyes were degraded, however, dramatic differences in the eyes, and their reliance on vision, exist in this group. In this study, we examined the rod opsin gene (RH1), and compared its evolution to that of two cone opsin genes (SWS1 and M/LWS). While phylogenetic reconstruction with the cone opsin genes SWS1 and M/LWS generated a species tree in accord with expectations, the RH1 gene tree united Pteropodidae (Old-World fruit bats) and Yangochiroptera, with very high bootstrap values, suggesting the possibility of convergent evolution. The hypothesis of convergent evolution was further supported when nonsynonymous sites or amino acid sequences were used to construct phylogenies. Reconstructed RH1 sequences at internal nodes of the bat species phylogeny showed that: (1) Old-World fruit bats share an amino acid change (S270G) with the tomb bat; (2) Miniopterus share two amino acid changes (V104I, M183L) with Rhinolophoidea; (3) the amino acid replacement I123V occurred independently on four branches, and the replacements L99M, L266V and I286V occurred each on two branches. The multiple parallel amino acid replacements that occurred in the evolution of bat RH1 suggest the possibility of multiple convergences of their ecological specialization (i.e., various photic environments) during adaptation for the nocturnal lifestyle, and suggest that further attention is needed on the study of the ecology and behavior of bats.

Integration of a photonic crystal polarization beam splitter and waveguide bend

In this work, we present the design of an integrated photonic-crystal polarization beam splitter (PC-PBS) and a low-loss photonic-crystal 60 waveguide bend. Firstly, the modal properties of the PC-PBS and the mechanism of the low-loss waveguide bend are investigated by the two-dimensional finite-difference time-domain (FDTD) method, and then the integration of the two devices is studied. It shows that, although the individual devices perform well separately, the performance of the integrated circuit is poor due to the multi-mode property of the PC-PBS. By introducing deformed airhole structures, a single-mode PC-PBS is proposed, which significantly enhance the performance of the circuit with the extinction ratios remaining above 20dB for both transverse-electric (TE) and transverse-magnetic (TM) polarizations. Both the specific result and the general idea of integration design are promising in the photonic crystal integrated circuits in the future. (C) 2009 Optical Society of America

Multiplicity factor and diffraction geometry factor for single crystal X-ray diffraction analysis and measurement of phase content in cubic GaN/GaAs(001) epilayers

On the basis of integrated intensity of rocking curves, the multiplicity factor and the diffraction geometry factor for single crystal X-ray diffraction (XRD) analysis were proposed and a general formula for calculating the content of mixed phases was obtained. With a multifunction four-circle X-ray double-crystal diffractometer, pole figures of cubic (002), {111} and hexagonal {1010} and reciprocal space mapping were measured to investigate the distributive character of mixed phases and to obtain their multiplicity factors and diffraction geometry factors. The contents of cubic twins and hexagonal inclusions were calculated by the integrated intensities of rocking curves of cubic (002), cubic twin {111}, hexagonal {1010} and {1011}.

A global shallow water model using high order multi-moment constrained finite volume method and icosahedral grid

A novel accurate numerical model for shallow water equations on sphere have been developed by implementing the high order multi-moment constrained finite volume (MCV) method on the icosahedral geodesic grid. High order reconstructions are conducted cell-wisely by making use of the point values as the unknowns distributed within each triangular cell element. The time evolution equations to update the unknowns are derived from a set of constrained conditions for two types of moments, i.e. the point values on the cell boundary edges and the cell-integrated average. The numerical conservation is rigorously guaranteed. in the present model, all unknowns or computational variables are point values and no numerical quadrature is involved, which particularly benefits the computational accuracy and efficiency in handling the spherical geometry, such as coordinate transformation and curved surface. Numerical formulations of third and fourth order accuracy are presented in detail. The proposed numerical model has been validated by widely used benchmark tests and competitive results are obtained. The present numerical framework provides a promising and practical base for further development of atmospheric and oceanic general circulation models. (C) 2009 Elsevier Inc. All rights reserved.

Directing single-walled carbon nanotubes to self-assemble at water/oil interfaces and facilitate electron transfer

Both the behavior and the general key factors for assembling flexible SWNT films at the water/oil interface were investigated; the electron transfer, one of the most fundamental chemical processes, at the SWNT-sandwiched water/oil interface was also firstly illustrated using scanning electrochemical microscopy.

Effective dielectric response of nonlinear composites with external ac and dc electric field

The perturbation method is developed to investigate the effective nonlinear dielectric response of Kerr composites when the external ac and dc electric field is applied. Under the external ac and dc electric field E-app=E-a(1+sin omegat), the effective coupling nonlinear response can be induced by the cubic nonlinearity of Kerr nonlinear materials at the zero frequency, the finite basic frequency omega, the second and the third harmonics, 2omega and 3omega, and so on. As an example, we have investigated the cylindrical inclusions randomly embedded in a host and derived the formulas of the effective nonlinear dielectric response at harmonics in dilute limit. For a higher concentration of inclusions, we have proposed a nonlinear effective-medium approximation by introducing the general effective nonlinear response. With the relationships between the effective nonlinear response at harmonics and the general effective nonlinear response, we have derived a set of formulas of the effective nonlinear dielectric responses at harmonics for a larger volume fraction. (C) 2004 American Institute of Physics.

The positive floating potential of high-temperature particles in plasma

It has been predicted that the floating potential of particles in plasma may become positive when the particle surface temperature is high enough, but, to our knowledge, no positive floating potential has been obtained yet. In the present paper the floating potential theory of high-temperature particles in plasma is developed to cover the positive potential range for the first time, and a general approximate analytical formula for the positive floating potential with a thin plasma sheath and subsonic plasma flow is derived from the new model recently proposed by the authors. The results show that when the floating potential is positive, the net flux of charge incident on the particle approaches a constant similar to the 'electron saturation' phenomena in the case of the electric probes.

A consistent model of isolated force-free magnetic arches

The magnetic flux tube concentrating strong magnetic field is the basic configuration of magneticfield in the solar atmosphere. In the present paper, the equilibrium of isolated magnetic flux tube inthe solar atmosphere is discussed. In the viewpoint of mathematics, the boundary condition is nonlinearand the position of boundary needs to be determined by the physical condition although the equation ofmagnetic potential is linear for the linear force-free field. Analytical solutions to the arches of bothuniform circular cross-section and non-uniform cross section have been obtained. The results show thatthe nonlinear problem may have or not have any solution according to different azimuthal components of the magnetic field; the number of solutions to the nonlinear problem is four at most, and two in some cases. In the present paper, the analytical solutions to the approximations of both fat and slender arches are given in detail, and the general features of magnetic arch structure are shown.

Choice of calibration equations in calibrating microwave test fixtures

The error theory of linear equation system has been applied to the calibration procedure of microwave network analyser in this article. A new explanation for the choice of the linear calibration equations is proposed and a general principle for choosing calibration equations is presented. The method can also be used to predict the occurrence of the problem of frequency limitation at some periodic frequencies. This principle is employed to the thru-short-delay (TSD) method and the solution using the chosen equations gives the most accurate results. A good agreement between the theory and the experiment has been obtained.

Experimental Investigation into the Production Behavior of Methane Hydrate in Porous Sediment with Hot Brine Stimulation

The gas production behavior from methane hydrate in porous sediment by injecting the brine with the salinity of 0−24 wt % and the temperature of −1 to 130 °C was investigated in a one-dimensional experimental apparatus. The results show that the gas production process consists of three periods: the free gas production, the hydrate dissociation, and the general gas reservoir production. The hydrate dissociation accompanies the temperature decrease with the injection of the brine (NaCl solution), and the dissociation duration is shortened with the increase of the salinity. With the injection of hot brine, instantaneous hydrate dissociation rate also increases with the increase of the salinity. However, while the NaCl concentration is beyond a certain value, the rate has no longer continued increasing. Thermal efficiency and energy ratio for the hydrate production can be enhanced by injecting hot brine, and the enhanced effectiveness is quite good with the injection of high salinity at lower temperature.

Numerical investigation of particle velocity distributions in aeolian sand transport

Particle velocity distribution in a blowing sand cloud is a reflection of saltation movement of many particles. Numerical analysis is performed for particle velocity distribution with a discrete particle model. The probability distributions of resultant particle velocity in the impact-entrainment process, particle horizontal and vertical velocities at different heights and the vertical velocity of ascending particles are analyzed. The probability distributions of resultant impact and lift-off velocities of saltating particles can be expressed by a log-normal function, and that of impact angle comply with an exponential function. The probability distribution of particle horizontal and vertical velocities at different heights shows a typical single-peak pattern. In the lower part of saltation layer, the particle horizontal velocity distribution is positively skewed. Further analysis shows that the probability density function of the vertical velocity of ascending particles is similar to the right-hand part of a normal distribution function, and a general equation is acquired for the probability density function of non-dimensional vertical velocity of ascending particles which is independent of diameter of saltating particles, wind strength and height. These distributions in the present numerical analysis are consistent with reported experimental results. The present investigation is important for understanding the saltation state in wind-blown sand movement. (C) 2009 Elsevier B.V. All rights reserved.

Statistical distribution of nonlinear random wave height

A statistical model of random wave is developed using Stokes wave theory of water wave dynamics. A new nonlinear probability distribution function of wave height is presented. The results indicate that wave steepness not only could be a parameter of the distribution function of wave height but also could reflect the degree of wave height distribution deviation from the Rayleigh distribution. The new wave height distribution overcomes the problem of Rayleigh distribution that the prediction of big wave is overestimated and the general wave is underestimated. The prediction of small probability wave height value of new distribution is also smaller than that of Rayleigh distribution. Wave height data taken from East China Normal University are used to verify the new distribution. The results indicate that the new distribution fits the measurements much better than the Rayleigh distribution.