Three-dimensional instability of an oscillating viscous flow past a circular cylinder

A systematically numerical study of the sinusoidally oscillating viscous flow around a circular cylinder was performed to investigate vortical instability by solving the three-dimensional incompressible Navier-Stokes equations. The transition from two- to three-dimensional flow structures along the axial direction due to the vortical instability appears, and the three-dimensional structures lie alternatively on the two sides of the cylinder. Numerical study has been taken for the Keulegan-Carpenter( KC) numbers from 1 to 3.2 and frequency parameters from 100 to 600. The force behaviors are also studied by solving the Morison equation. Calculated results agree well with experimental data and theoretical prediction.

Material parameter extraction in THz-TDS using a converging beam transfer function

We demonstrate a parameter extraction algorithm based on a theoretical transfer function, which takes into account a converging THz beam. Using this, we successfully extract material parameters from data obtained for a quartz sample with a THz time domain spectrometer. © 2010 IEEE.

A bulk acoustic mode single-crystal silicon microresonator with a high-quality factor

This paper details a bulk acoustic mode resonator fabricated in single-crystal silicon with a quality factor of 15 000 in air, and over a million below 10 mTorr at a resonant frequency of 2.18 MHz. The resonator is a square plate that is excited in the square-extensional mode and has been fabricated in a commercial foundry silicon-on-insulator (SOI) MEMS process through MEMSCAP. This paper also presents a simple method of extracting resonator parameters from raw measurements heavily buried in electrical feedthrough. Its accuracy has been demonstrated through a comparison between extracted motional resistance values measured at different voltage biases and those predicted from an analytical model. Finally, a method of substantially cancelling electrical feedthrough through system-level electronic implementation is also introduced. © 2008 IOP Publishing Ltd.

Computation of stress intensity factors by the sub-region mixed finite element method of lines

Based on the sub-region generalized variational principle, a sub-region mixed version of the newly-developed semi-analytical 'finite element method of lines' (FEMOL) is proposed in this paper for accurate and efficient computation of stress intensity factors (SIFs) of two-dimensional notches/cracks. The circular regions surrounding notch/crack tips are taken as the complementary energy region in which a number of leading terms of singular solutions for stresses are used, with the sought SIFs being among the unknown coefficients. The rest of the arbitrary domain is taken as the potential energy region in which FEMOL is applied to obtain approximate displacements. A mixed system of ordinary differential equations (ODEs) and algebraic equations is derived via the sub-region generalized variational principle. A singularity removal technique that eliminates the stress parameters from the mixed equation system eventually yields a standard FEMOL ODE system, the solution of which is no longer singular and is simply and efficiently obtained using a standard general-purpose ODE solver. A number of numerical examples, including bi-material notches/cracks in anti-plane and plane elasticity, are given to show the generally excellent performance of the proposed method.

Monte Carlo Valuation of natural gas investments

This paper deals with the valuation of energy assets related to natural gas. In particular, we evaluate a baseload Natural Gas Combined Cycle (NGCC) power plant and an ancillary instalation, namely a Liquefied Natural Gas (LNG) facility, in a realistic setting; specifically, these investments enjoy a long useful life but require some non-negligible time to build. Then we focus on the valuation of several investment options again in a realistic setting. These include the option to invest in the power plant when there is uncertainty concerning the initial outlay, or the option's time to maturity, or the cost of CO2 emission permits, or when there is a chance to double the plant size in the future. Our model comprises three sources of risk. We consider uncertain gas prices with regard to both the current level and the long-run equilibrium level; the current electricity price is also uncertain. They all are assumed to show mean reversion. The two-factor model for natural gas price is calibrated using data from NYMEX NG futures contracts. Also, we calibrate the one-factor model for electricity price using data from the Spanish wholesale electricity market, respectively. Then we use the estimated parameter values alongside actual physical parameters from a case study to value natural gas plants. Finally, the calibrated parameters are also used in a Monte Carlo simulation framework to evaluate several American-type options to invest in these energy assets. We accomplish this by following the least squares MC approach.

Population characteristics of gulf menhaden, Brevoortia patronus

The status of the Gulf menhaden, Brevoortia patronus, fishery was assessed with purse-seine landings data from 1946 to 1997 and port sampling data from 1964 to 1997. These data were analyzed to determine growth rates, biological reference points for fi shing mortality from yield per recruit and maximum spawning potential analyses, spawner-recruit relationships, and maximum sustainable yield (MSY). The separable virtual population approach was used for the period 1976–97 (augmented by earlier analyses for 1964–75) to obtain point estimates of stock size, recruits to age 1, spawning stock size, and fishing mortality rates. Exploitation rates for age-1 fi sh ranged between 11% and 45%, for age-2 fi sh between 32% and 72%, and for age-3 fi sh between 32% and 76%. Biological reference points from yield per recruit (F0.1: 1.5–2.5/yr) and spawning potential ratio (F20: 1.3–1.9/yr and F30: 0.8–1.2/yr) were obtained for comparison with recent estimates of F (0.6–0.8/yr). Recent spawning stock estimates (as biomass or eggs) are above the long-term average, while recent recruits to age 1 are comparable to the long-term average. Parameters from Ricker-type spawner-recruit relations were estimated, although considerable unexplained variability remained. Recent survival to age-1 recruitment has generally been below that expected based on the Ricker spawner-recruit relation. Estimates of long-term MSY from PRODFIT and ASPIC estimation of production model ranged between 717,000 t and 753,000 t, respectively. Declines in landings between 1988 and 1992 raised concerns about the status of the Gulf menhaden stock. Landings have fl uctuated without trend since 1992, averaging about 571,000 t. However, Gulf menhaden are short lived and highly fecund. Thus, variation in recruitment to age 1, largely mediated by environmental conditions, infl uences fi shing success over the next two years (as age-1 and age-2 fi sh). Comparisons of recent estimates of fi shing mortality to biological reference points do not suggest overfishing. (PDF file contains 22 pages.)

Population Characteristics of Gulf Menhaden, Brevoortia patronus

The status of the gulf menhaden, Brevoortia patronus, fishery was assessed with purseseine landing data from 1946 to 1992 and port sampling data from 1964 to 1992. These data were analyzed to determine growth rates, biological reference points for fishing mortality from yield per recruit and maximum spawning potential analyses, spawner-recruit relationships, and maximum sustainable yield (MSY). Virtual population approaches were used to obtain point estimates of stock size, recruits to age I, spawning stock size, and fishing mortality rates. Exploitation rates ranged between 14% and 45% for age-1 fish, between 30% and 72% for age-2 fish, and between 36% and 71% for age-3 fish. Biological reference points from yield per recruit (FO. I: 0.7-0.9 yr-1) and maximum spawning potential (F20: 1.62.9 yr-l and F30: 1.0-2.1 yr-1) were obtained for comparison with recent estimates of F (0.4-0.8 yr-l). Parameters from Ricker-type spawner-recruit relations were estimated, although considerable unexplained variability remained. Estimates of long-term MSY from fits of the generalized production model ranged between 664,000 metric tons (t) and 897,000 t. Declines in landings since 1988 have raised concerns about the status of the gulf menhaden stock. However, gulf menhaden are short lived and highly fecund. Thus, variation in recruitment to age 1 largely mediated by environmental conditions influences fishing success over the next two years (as age-1 and age-2 fish). Comparisons of recent estimates of fishing mortality to biological reference points do not suggest overfishing. (PDF file contains 26 pages.)

Sensor networks for geospatial event detection - theory and applications

This thesis presents theories, analyses, and algorithms for detecting and estimating parameters of geospatial events with today's large, noisy sensor networks. A geospatial event is initiated by a significant change in the state of points in a region in a 3-D space over an interval of time. After the event is initiated it may change the state of points over larger regions and longer periods of time. Networked sensing is a typical approach for geospatial event detection. In contrast to traditional sensor networks comprised of a small number of high quality (and expensive) sensors, trends in personal computing devices and consumer electronics have made it possible to build large, dense networks at a low cost. The changes in sensor capability, network composition, and system constraints call for new models and algorithms suited to the opportunities and challenges of the new generation of sensor networks. This thesis offers a single unifying model and a Bayesian framework for analyzing different types of geospatial events in such noisy sensor networks. It presents algorithms and theories for estimating the speed and accuracy of detecting geospatial events as a function of parameters from both the underlying geospatial system and the sensor network. Furthermore, the thesis addresses network scalability issues by presenting rigorous scalable algorithms for data aggregation for detection. These studies provide insights to the design of networked sensing systems for detecting geospatial events. In addition to providing an overarching framework, this thesis presents theories and experimental results for two very different geospatial problems: detecting earthquakes and hazardous radiation. The general framework is applied to these specific problems, and predictions based on the theories are validated against measurements of systems in the laboratory and in the field.

Fermionic quantum systems. Part I: Phase transitions in quantum dots. Part II: Nuclear matter on a lattice

In the first part I perform Hartree-Fock calculations to show that quantum dots (i.e., two-dimensional systems of up to twenty interacting electrons in an external parabolic potential) undergo a gradual transition to a spin-polarized Wigner crystal with increasing magnetic field strength. The phase diagram and ground state energies have been determined. I tried to improve the ground state of the Wigner crystal by introducing a Jastrow ansatz for the wave function and performing a variational Monte Carlo calculation. The existence of so called magic numbers was also investigated. Finally, I also calculated the heat capacity associated with the rotational degree of freedom of deformed many-body states and suggest an experimental method to detect Wigner crystals.

The second part of the thesis investigates infinite nuclear matter on a cubic lattice. The exact thermal formalism describes nucleons with a Hamiltonian that accommodates on-site and next-neighbor parts of the central, spin-exchange and isospin-exchange interaction. Using auxiliary field Monte Carlo methods, I show that energy and basic saturation properties of nuclear matter can be reproduced. A first order phase transition from an uncorrelated Fermi gas to a clustered system is observed by computing mechanical and thermodynamical quantities such as compressibility, heat capacity, entropy and grand potential. The structure of the clusters is investigated with the help two-body correlations. I compare symmetry energy and first sound velocities with literature and find reasonable agreement. I also calculate the energy of pure neutron matter and search for a similar phase transition, but the survey is restricted by the infamous Monte Carlo sign problem. Also, a regularization scheme to extract potential parameters from scattering lengths and effective ranges is investigated.

Modelagem do Equilíbrio de Fases entre Hidrocarbonetos Leves e Pesados

O comportamento de fases para sistemas binários com um hidrocarboneto leve e um pesado é muito importante tanto para o projeto real de um processo quanto para o desenvolvimento de modelos teóricos. Para atender a crescente demanda por informação experimental de equilíbrio de fases a altas pressões, o objetivo deste estudo é obter uma metodologia que substitua parcialmente ou maximize a pouca informação experimental disponível. Para isto propõe-se a modelagem do equilíbrio de fases em misturas de hidrocarboneto leve com um pesado, sem o conhecimento da estrutura molecular do pesado, inferindo-se os parâmetros do modelo a partir da modelagem de dados de ponto de bolha obtidos na literatura. Esta metodologia implica não só na descrição do equilíbrio de fases de um sistema como na estimação das propriedades críticas do pesado, de difícil obtenção devido ao craqueamento destes a altas temperaturas. Neste contexto, este estudo apresenta uma estratégia que estima indiretamente as propriedades críticas dos compostos pesados. Para isto, foram correlacionados dados experimentais de ponto de bolha de misturas binárias contendo um hidrocarboneto leve e um pesado, usando-se dois modelos: o de Peng-Robinson e o TPT1M (Teoria da Polimerização Termodinâmica de primeira ordem de Wertheim modificada). Os parâmetros ajustados com o modelo de Peng-Robinson correspondem diretamente às propriedades críticas do composto pesado, enquanto os ajustados com o modelo TPT1M foram usados para obtê-las. Esta estratégia fornece parâmetros dependentes do modelo, porém permite o cálculo de outras propriedades termodinâmicas, como a extrapolação da temperatura dos dados estudados. Além disso, acredita-se que a correlação dos parâmetros obtidos com as propriedades críticas disponíveis ajudará na caracterização de frações pesadas de composição desconhecida

A singular perturbation method of calculating the behavior of supercavitating hydrofoils with rounded noses

A simple, direct and accurate method to predict the pressure distribution on supercavitating hydrofoils with rounded noses is presented. The thickness of body and cavity is assumed to be small. The method adopted in the present work is that of singular perturbation theory. Far from the leading edge linearized free streamline theory is applied. Near the leading edge, however, where singularities of the linearized theory occur, a non-linear local solution is employed. The two unknown parameters which characterize this local solution are determined by a matching procedure. A uniformly valid solution is then constructed with the aid of the singular perturbation approach.

The present work is divided into two parts. In Part I isolated supercavitating hydrofoils of arbitrary profile shape with parabolic noses are investigated by the present method and its results are compared with the new computational results made with Wu and Wang's exact "functional iterative" method. The agreement is very good. In Part II this method is applied to a linear cascade of such hydrofoils with elliptic noses. A number of cases are worked out over a range of cascade parameters from which a good idea of the behavior of this type of important flow configuration is obtained.

Some of the computational aspects of Wu and Wang's functional iterative method heretofore not successfully applied to this type of problem are described in an appendix.

Synthesis and magnetic properties of Ni-doped zinc oxide powders

Polycrystalline Zn1-xNixO diluted magnetic semiconductors have been successfully synthesized by an auto-combustion method. X-ray diffraction measurements indicated that the 5 at% Ni-cloped ZnO had the pure wurtzite structure. Refinements of cell parameters from powder diffraction data revealed that the cell parameters of Zn0.95Ni0.05O were a little bit larger than ZnO. Transmission electron microscopy observation showed that the as-synthesized powders were of the size similar to 60 nm. Magnetic investigations showed that the nanocystalline Zn0.95Ni0.05O possessed room temperature ferromagnetisin with the saturation magnetic moment of 0.1 emu/g (0.29 mu(B)/Ni2+). (c) 2005 Elsevier B.V. All rights reserved.

Demographic characteristics of exploited tropical lutjanids: a comparative analysis

Demographic parameters from seven exploited coral reef lutjanid species were compared as a case study of the implications of intrafamily variation in life histories for multispecies harvest management. Modal lengths varied by 4 cm among four species (Lutjanus fulviflamma, L. vitta, L. carponotatus, L. adetii), which were at least 6 cm smaller than the modal lengths of the largest species (L. gibbus, Symphorus nematophorus, Aprion virescens). Modal ages, indicating ages of full selection to fishing gear, were 10 years or less for all species, but maximum ages ranged from 12 (L. gibbus) to 36 years (S. nematophorus). Each species had a unique growth pattern, with differences in length-at-age and mean asymptotic fork length (L∞), but smaller species generally grew fast during the first 1–2 years of life and larger species grew more slowly over a longer period. Total mortality rates varied among species; L. gibbus had the highest mortality and L. fulviflamma, the lowest mortality. The variability in life history strategies of these tropical lutjanids makes generalizations about lutjanid life histories difficult, but the fact that all seven had characteristics that would make them particularly vulnerable to fishing indicates that harvest of tropical lutjanids should be managed with caution.

Critério estatístico para obtenção de valores de NSPT para previsão da capacidade de carga de estacas por métodos semi empíricos.

Uma das tarefas mais desafiadoras do engenheiro na área da Geotecnia é a escolha dos valores de parâmetros geotécnicos obtidos de ensaios de campo ou laboratório e que serão utilizados nos modelos analíticos ou numéricos na fase de projeto de fundações. Diante das incertezas inerentes aos ensaios de SPT e da heterogeneidade de abordagens para a utilização dos valores de NSPT, é proposta neste estudo, a aplicação de um critério estatístico para obtenção de valores de NSPT, a partir da construção de intervalos de confiança de 95% de probabilidade em torno da reta ajustada de regressão linear simples entre a variável aleatória NSPT e a profundidade. Os valores obtidos de NSPT pelo critério aplicado foram utilizados na previsão da capacidade de carga de 19 estacas isoladas a partir da utilização de três métodos semi-empíricos: Aoki-Velloso (1975) com coeficientes alterados por Monteiro (1997), Décourt & Quaresma (1978) alterado pelo método de Décourt (1996) e Método de Alonso (1996). As cargas de ruptura dessas 19 estacas ensaiadas através de Provas de Carga Estática foram obtidas pelos métodos de extrapolação de Van Der Veen (1953) e Décourt (1996) e serviram para comparação e consequente validação do critério estatístico. Adicionalmente, com fulcro no item 6.2.1.2.1 da ABNT NBR 6122:2010 Resistência calculada por método semi-empírico, foram avaliados os fatores de segurança em relação às cargas de projeto, inclusive, também se utilizando da premissa de reconhecimento de regiões representativas, levando em conta o número de ensaios de SPT executados, fato que promove uma diminuição da incerteza dos parâmetros, apontando a um menor fator de segurança. A dissertação enfatiza as vantagens de um adequado tratamento estatístico dos parâmetros geotécnicos, a exemplo da recomendação já existente nas normas internacionais como Eurocódigo e outras. O critério construído permite e encoraja análises e decisões racionais no universo das partes interessadas consumidores, projetistas, fiscais de obras, contratantes e comunidade científica promovendo as discussões de forma mais objetiva e harmoniosa sobre o tema.

Análise das deformações dos aterros que compõem a obra do Arco Metropolitano do Rio de Janeiro.

O presente trabalho tem como objetivo a análise numérica do comportamento dos aterros instrumentados que compõem a obra do Arco Metropolitano do Rio de Janeiro. Os resultados da instrumentação de campo serão retroanalisados, juntamente com dados disponíveis na literatura, buscando-se a definição de parâmetros confiáveis, representativos do comportamento da argila compressível da região em estudo. O Arco Metropolitano do Rio de Janeiro é uma rodovia projetada que servirá como ligação entre as principais rodovias que cortam o município do Rio de Janeiro. Devido a presença de grandes espessuras de solo compressível em alguns trechos da região, cortados pelo traçado da rodovia, instrumentos de monitoramento, como placas de recalque e inclinômetros, foram utilizados para avaliar os deslocamentos verticais e horizontais dos aterros durante o processo construtivo. Para este trabalho foram selecionados trechos de aterros representativos, devido à magnitude dos recalques, qualidade do resultado da instrumentação e diferentes métodos construtivos. A partir da análise dos parâmetros de projeto e dos parâmetros encontrados na literatura, procede-se à simulação numérica do processo construtivo dos aterros selecionados com o programa PLAXIS, de elementos finitos, através de modelagem bidimensional. Os resultados numéricos são confrontados com a instrumentação de campo (placas de recalque e inclinômetros) e com os resultados de previsões teóricas (teoria de adensamento unidimensional). Os resultados comprovaram que a modelagem numérica mostrou-se uma ferramenta adequada para a previsão dos recalques totais, tempos de adensamento e ganho de resistência ao longo do tempo. A retroanálise do comportamento de aterros sobre solos moles permite a reavaliação das premissas de projeto, uma vez que as limitações das teorias de análise e a dificuldade na seleção de parâmetros, muitas vezes acarretam em estimativas de recalque incoerentes com as observações de campo.