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.

周期性极化铌酸锂晶体中半非共线型宽带光学参变放大理论研究

基于周期性极化铌酸锂晶体(PPLN)的准相位匹配光参变放大过程,通过倾斜周期极化铌酸锂晶体中极化域(极化光栅)一定角度,实现了介于共线匹配方式和非共线匹配方式之间的一种半非共线型准相位匹配方式,并以该匹配方式下的各光矢量几何关系得出相位匹配曲线,找到在特定抽运光和信号光波长下能获得宽带增益放大的周期极化长度。并研究其极化倾斜角度与温度特性。模拟计算表明,在合适的角度与温度条件下,该方式可以532 nm抽运光抽运的信号光在800 nm和1064 nm处均获得宽带光参变放大。

A fully-nonlocal energy-based formulation and high-performance realization of the quasicontinuum method

The quasicontinuum (QC) method was introduced to coarse-grain crystalline atomic ensembles in order to bridge the scales from individual atoms to the micro- and mesoscales. Though many QC formulations have been proposed with varying characteristics and capabilities, a crucial cornerstone of all QC techniques is the concept of summation rules, which attempt to efficiently approximate the total Hamiltonian of a crystalline atomic ensemble by a weighted sum over a small subset of atoms. In this work we propose a novel, fully-nonlocal, energy-based formulation of the QC method with support for legacy and new summation rules through a general energy-sampling scheme. Our formulation does not conceptually differentiate between atomistic and coarse-grained regions and thus allows for seamless bridging without domain-coupling interfaces. Within this structure, we introduce a new class of summation rules which leverage the affine kinematics of this QC formulation to most accurately integrate thermodynamic quantities of interest. By comparing this new class of summation rules to commonly-employed rules through analysis of energy and spurious force errors, we find that the new rules produce no residual or spurious force artifacts in the large-element limit under arbitrary affine deformation, while allowing us to seamlessly bridge to full atomistics. We verify that the new summation rules exhibit significantly smaller force artifacts and energy approximation errors than all comparable previous summation rules through a comprehensive suite of examples with spatially non-uniform QC discretizations in two and three dimensions. Due to the unique structure of these summation rules, we also use the new formulation to study scenarios with large regions of free surface, a class of problems previously out of reach of the QC method. Lastly, we present the key components of a high-performance, distributed-memory realization of the new method, including a novel algorithm for supporting unparalleled levels of deformation. Overall, this new formulation and implementation allows us to efficiently perform simulations containing an unprecedented number of degrees of freedom with low approximation error.

Tm:Y_2SiO_5晶体的生长和光谱性质研究

采用中频感应提拉法生长了高质量的Tm:Y2SiO5(Tm:YSO)晶体,测定了晶体的晶格常数和分凝系数.运用劳厄照相法确定了单斜晶系Tm:YSO晶体的三个偏振轴〈010〉,D1和D2,在室温下测量了三个偏振轴方向的吸收光谱、荧光光谱和荧光寿命,计算了晶体吸收峰的吸收线宽和吸收截面.研究发现,相对于其他两个偏振轴方向,D1方向在790 nm处出现较强的吸收峰,同时在2μm附近出现了一定强度的发射峰,D1方向的吸收截面较大,荧光寿命较长.Tm:YSO晶体适用于AlGaAs二极管抽运固体激光器,在2μm波段固体激光器的应用上将有很大的发展潜力.

Dividing population genetic distance data with the software Partitioning Optimization with Restricted Growth Strings (PORGS): an application for Chinook salmon (Oncorhynchus tshawytscha), Vancouver Island, British Columbia

A new method of finding the optimal group membership and number of groupings to partition population genetic distance data is presented. The software program Partitioning Optimization with Restricted Growth Strings (PORGS), visits all possible set partitions and deems acceptable partitions to be those that reduce mean intracluster distance. The optimal number of groups is determined with the gap statistic which compares PORGS results with a reference distribution. The PORGS method was validated by a simulated data set with a known distribution. For efficiency, where values of n were larger, restricted growth strings (RGS) were used to bipartition populations during a nested search (bi-PORGS). Bi-PORGS was applied to a set of genetic data from 18 Chinook salmon (Oncorhynchus tshawytscha) populations from the west coast of Vancouver Island. The optimal grouping of these populations corresponded to four geographic locations: 1) Quatsino Sound, 2) Nootka Sound, 3) Clayoquot +Barkley sounds, and 4) southwest Vancouver Island. However, assignment of populations to groups did not strictly reflect the geographical divisions; fish of Barkley Sound origin that had strayed into the Gold River and close genetic similarity between transferred and donor populations meant groupings crossed geographic boundaries. Overall, stock structure determined by this partitioning method was similar to that determined by the unweighted pair-group method with arithmetic averages (UPGMA), an agglomerative clustering algorithm.

Compressão de dados sísmicos com perda controlada

Neste trabalho apresentamos um novo método de compressão, com perda controlada de dados, que tem a vantagem de ter uma taxa significativa de compressão sem introduzir nenhuma perda superior a um parâmetro escolhido pelo usuário. Esta abordagem é uma abordagem mista, pois usa técnicas de compactação de dados tanto com perda quanto sem perda. Isto quer dizer que conseguimos um método que alia as vantagens da alta compressão, sem introduzir distorções indesejáveis nos dados. Mostramos como a massa de dados utilizada nos nossos estudos é obtida e a sua importância na prospecção de depósitos de hidrocarbonetos. É apresentado um levantamento bibliográfico com técnicas de compressão aplicadas a dados sísmicos tipicamente utilizadas em aplicações comerciais. Por fim, apresentamos os resultados da compressão utilizando o método em conjuntos de dados sísmicos reais. Para 1% de erro, os arquivos de dados sísmicos compactados passaram a ter algo próximo a 25% de seus tamanhos originais, o que representa um fator de compressão de aproximadamente 4

Propriedades luminescentes do SrGa2O4 dopado com íons de Ni2+

O objetivo deste trabalho é a síntese e investigação estrutural e óptica de amostras SrGa2O4 dopados com 1% de íons Ni2+. Estas amostras foram sintetizados por reação do estado sólido convencional, utilizando como materiais de partida de alta pureza Ga2O3, SrCO3 e NiO em quantidades estequiométricas. As amostras foram caracterizadas estruturalmente pelo método de difração de raios - X( XRD ) e as medições de difração mostraram que as amostras têm uma única fase monoclínica. Os padrões de XRD também foram refinados pelo método de Rietveld, que permitiu a determinação dos parâmetros de célula unitária. A Caracterização óptica das amostras puras e dopadas SrGa2O4 foram realizadas as medições a partir de fotoluminescência, de excitação e de absorção fotoacústica, à temperatura ambiente. Os espectros de emissão mostraram três bandas de emissão localizadas em 557 nm, 661 nm e 844 nm e foram identificadas essas bandas, respectivamente, com as seguintes transições eletrônicas :1T2 (1D) → 3A2 (3F), 3T1 (3F)→ 3A2 (3F) e 1T2 (1D) → 3T2 (3F). Os espectros de excitação mostraram seis bandas de absorção associadas às transições electrônicas do nível 3A2 (3F) para o 3T1 (3P) , T1 (3P), 1A1 (1G), 1T2 (1D), 3T1 (3F), 1E (1D) e 1T2, 1E (1G). Medidas de absorção fotoacústica também foram realizados com o fim de verificar as transições ópticas observadas nos espectros de excitação e de identificar novas bandas de absorção óptica. Os resultados demonstraram que os íons de Ni2+ ocupam dois locais octaédricos diferentes na amostra SrGa2O4 dopado. A partir das transições ópticas observadas nos espectros de excitação e fotoacústica, determinou-se o parâmetro de cristal de campo, dq, e parâmetros Racah, B e C. A proporção Dq / B ≈ 1.2 para ambos os locais são típicos para Ni2+ íons inseridos em redes de óxido e em coordenação octaédrica.

Implementação paralelizada de métodos de resolução de sistemas algébricos na simulação de reservatórios de gás

Desde a década de 1960, devido à pertinência para a indústria petrolífera, a simulação numérica de reservatórios de petróleo tornou-se uma ferramenta usual e uma intensa área de pesquisa. O principal objetivo da modelagem computacional e do uso de métodos numéricos, para a simulação de reservatórios de petróleo, é o de possibilitar um melhor gerenciamento do campo produtor, de maneira que haja uma maximização na recuperação de hidrocarbonetos. Este trabalho tem como objetivo principal paralelizar, empregando a interface de programação de aplicativo OpenMP (Open Multi-Processing), o método numérico utilizado na resolução do sistema algébrico resultante da discretização da equação que descreve o escoamento monofásico em um reservatório de gás, em termos da variável pressão. O conjunto de equações governantes é formado pela equação da continuidade, por uma expressão para o balanço da quantidade de movimento e por uma equação de estado. A Equação da Difusividade Hidráulica (EDH), para a variável pressão, é obtida a partir deste conjunto de equações fundamentais, sendo então discretizada pela utilização do Método de Diferenças Finitas, com a escolha por uma formulação implícita. Diferentes testes numéricos são realizados a fim de estudar a eficiência computacional das versões paralelizadas dos métodos iterativos de Jacobi, Gauss-Seidel, Sobre-relaxação Sucessiva, Gradientes Conjugados (CG), Gradiente Biconjugado (BiCG) e Gradiente Biconjugado Estabilizado (BiCGStab), visando a uma futura aplicação dos mesmos na simulação de reservatórios de gás. Ressalta-se que a presença de heterogeneidades na rocha reservatório e/ou às não-linearidades presentes na EDH para o escoamento de gás aumentam a necessidade de métodos eficientes do ponto de vista de custo computacional, como é o caso de estratégias usando OpenMP.

Plankton community composition in the Three Gorges Reservoir Region revealed by PCR-DGGE and its relationships with environmental factors

To explore the relationships between community composition and the environment in a reservoir ecosystem, plankton communities from the Three Gorges Reservoir Region were studied by PCR-denaturing gradient gel electrophoresis fingerprinting. Bacterial and eukaryotic operational taxonomic units (OTUs), generated by DGGE analysis of the PCR-amplified 16S and 18S rRNA genes, were used as surrogates for the dominant "biodiversity units". OTU composition among the sites was heterogeneous; 46.7% of the total bacteria] OTUs (45) and 64.1% of the eukaryotic OTUs (39) were identified in less than half of the sampling sites. Unweighted pair group method with arithmetic averages (UPGMA) clustering of the OTUs suggested that the plankton communities in the Xiangxi Rive sites were not always significantly different from those from the Yangtze River sites, despite clear differences in their environmental characterizations. Canonical correspondence analysis (CCA) was applied to further investigate the relationships between OTU composition and the environmental factors. The first two CCA ordination axes suggested that the bacterial community composition was primarily correlated with the variables of NO3--N, dissolved oxygen (DO), and SiO32--Si, whereas, the eukaryotic community was mainly correlated with the concentrations of DO, PO43--P, and SiO32--Si.

Highly efficient photoluminescence of Er2SiO5 films grown by reactive magnetron sputtering method

E2SiO5 thin films were fabricated on Si substrate by reactive magnetron sputtering method with subsequent annealing treatment. The morphology properties of as-deposited films have been studied by scanning electron microscope. The fraction of erbium is estimated to be 23.5 at% based on Rutherford backscattering measurement in as-deposited Er-Si-O film. X-ray diffraction measurement revealed that Er2SiO5 crystalline structure was formed as sample treated at 1100 degrees C for 1 h in O-2 atmosphere. Through proper thermal treatment, the 1.53 mu m Er3+-related emission intensity can be enhanced by a factor of 50 with respect to the sample annealed at 800 degrees C. Analysis of pump-power dependence of Er3+ PL intensity indicated that the upconversion phenomenon could be neglected even under a high photon flux of 10(21) (photons/cm(2)/sec). Temperature-dependent photoluminescence (PL) of Er2SiO5 was studied and showed a weak thermal quenching factor of 2. Highly efficienct photoluminescence of Er2SiO5 films has been demonstrated with Er3+ concentration of 10(22)/cm(3), and it opens a promising way towards future Si-based light source for Si photonics. (C) 2009 Elsevier B.V. All rights reserved.

Mode Characteristics for Square Resonators With a Metal Confinement Layer

Microsquare resonators laterally confined by SiO2/Au/air multilayer structure are investigated by light ray method with reflection phase-shift of the multiple layers and two-dimensional (2-D) finite-difference time-domain (FDTD) technique. The reflectivity and phase shift of the mode light ray on the sides of the square resonator with the semiconductor/SiO2/Au/air multilayer structure are calculated for TE and TM modes by transfer matrix method. Based on the reflection phase shift and the reflectivity, the mode wavelength and factor are calculated by the resonant condition and the mirror loss, which are in agreement well with that obtained by the FDTD simulation. We find that the mode factor increases greatly with the increase of the SiO2 layer thickness, especially as d < 0.3 mu m. For the square resonator with side length 2 mu m and refractive index 3.2, anticrossing mode couplings are found for confined TE modes at wavelength about 1.6 mu m at d = 0.11 mu m, and confined TM modes at d = 0.71 mu m, respectively.

Single-mode condition for silicon rib waveguides with trapezoidal cross-section

Single-mode condition for silicon rib waveguides with trapezoidal cross-section was obtained using a numerical method based on imaginary-distance beam propagation method with non-uniform discretization. Both quasi-transverse-electric and quasi-transverse-magnetic modes were investigated. Simulated single-mode condition is given by a modified equation. Comparison with reported results shows that the Marcatili's method is in a better agreement with our results. (C) 2003 Elsevier B.V. All rights reserved.

1.55 mu m Ge islands resonant-cavity-enhanced narrowband detector

The high quality Ge islands material with 1.55 mu m photo-response grown on Sol substrate is reported. Due to the modulation of the cavity formed by the mirrors at the surface and the buried SiO2 interface, seven sharp and strong peaks with narrow linewidth are found. And a 1.55 mu m Ge islands resonant-cavity-enhanced (RCE) detector with narrowband was fabricated by a simple method. The bottom mirror was deposited in the hole formed by anisotropically etching, in a basic solution from the backside of the sample with the buried SiO2 layer in silicon-on-insulator substrate as the etch-stop layer. Reflectivity spectrum indicates that the mirror deposited in the hole has a reflectivity as high as 99% in the range of 1.2-1.65 mu m. The peak responsivity of the RCE detector at 1543.8 nm is 0.028 mA/W and a full width at half maximum of 5 nm is obtained. Compared with the conventional p-i-n photodetector, the responsivity of RCE detector has a nearly threefold enhancement.

Free radial oscillation of a sphere in relativistic elasticity theory

The spherically symmetric free radial oscillation in the first post-Newtonian approximation for a homogeneous and isotropic elastic sphere with a constant density is studied. Based on the Xu, Wu, and Soffel formalism, the relation of the oscillation frequency of the sphere with the radius, mass density, and elastic constants of the sphere is derived by using the successive approximation method.

Fluorescence Hydrogen Peroxide Probe Based on a Microstructured Polymer Optical Fiber Modified with a Titanium Dioxide Film

A highly sensitive microstructured polymer optical fiber (MPOF) probe for hydrogen peroxide was made by forming a rhodamine 6G-doped titanium dioxide film on the side walls of array holes in an MPOF. It was found that hydrogen peroxide only has a response to the MPOF probe in a certain concentration of potassium iodide in sulfuric acid solution. The calibration graph of fluorescence intensity versus hydrogen peroxide concentration is linear in the range of 1.6 x 10(-7) mol/L to 9.6 x 10(-5) mol/L. The method, with high sensitivity and a wide linear range, has been applied to the determination of trace amounts of hydrogen peroxide in a few real samples, such as rain water and contact lens disinfectant, with satisfactory results.