Comparison of habitat-based indices of abundance with fishery-independent biomass estimates from bottom trawl surveys

Rockfish species are notoriously difficult to sample with multispecies bottom trawl survey methods. Typically, biomass estimates have high coefficients of variation and can fluctuate outside the bounds of biological reality from year to year. This variation may be due in part to their patchy distribution related to very specific habitat preferences. We successfully modeled the distribution of five commercially important and abundant rockf ish species. A two-stage modeling method (modeling both presence-absence and abundance) and a collection of important habitat variables were used to predict bottom trawl survey catch per unit of effort. The resulting models explained between 22% and 66% of the variation in rockfish distribution. The models were largely driven by depth, local slope, bottom temperature, abundance of coral and sponge, and measures of water column productivity (i.e., phytoplankton and zooplankton). A year-effect in the models was back-transformed and used as an index of the time series of abundance. The abundance index trajectories of three of five species were similar to the existing estimates of their biomass. In the majority of cases the habitat-based indices exhibited less interannual variability and similar precision when compared with stratified survey-based biomass estimates. These indices may provide for stock assessment models a more stable alternative to current biomass estimates produced by the multispecies bottom trawl survey in the Gulf of Alaska.

Propagação de ondas em teorias alternativas da gravitação

Uma forma de generalizar a teoria de Einstein da gravitação é incorporar na lagrangiana termos que dependem de escalares formados com os tensores de Ricci e Riemann, tais como (Ricci)2, ou (Riemann)2. Estas teorias tem sido estudadas intensamente nos últimos anos, já que elas podem ser usadas para descrever a expansão acelerada do universo no modelo cosmológico standard. Entre os desfios de modificar a teoria de Einstein, se encontra o de limitar a ambiguidade na escolha da dependência da lagrangiana com os escalares antes mencionados. A proposta desta dissertação é a de colocar limites sobre as possíveis lagrangianas impondo que as ondas (isto é, perturbações lineares) se propaguem no vácuo sem que apareça, shocks.

Variability in supply and cross-shelf transport of pink shrimp (Farfantepenaeus duorarum) postlarvae into western Florida Bay

The variability in the supply of pink shrimp (Farfantepenaeus duorarum) postlarvae and the transport mechanisms of planktonic stages were investigated with field data and simulations of transport. Postlarvae entering the nursery grounds of Florida Bay were collected for three consecutive years at channels that connect the Bay with the Gulf of Mexico, and in channels of the Middle Florida Keys that connect the southeastern margin of the Bay with the Atlantic Ocean. The influx of postlarvae in the Middle Florida Keys was low in magnitude and varied seasonally and among years. In contrast, the greater postlarval influx occurred at the northwestern border of the Bay, where there was a strong seasonal pattern with peaks in influx from July through September each year. Planktonic stages need to travel up to 150 km eastward between spawning grounds (northeast of Dry Tortugas) and nursery grounds (western Florida Bay) in about 30 days, the estimated time of planktonic development for this species. A Lagrangian trajectory model was developed to estimate the drift of planktonic stages across the SW Florida shelf. The model simulated the maximal distance traveled by planktonic stages under various assumptions of behavior. Simulation results indicated that larvae traveling with the instantaneous current and exhibiting a diel behavior travel up to 65 km and 75% of the larvae travel only 30 km. However, the eastward distance traveled increased substantially when a larval response to tides was added to the behavioral variable (distance increased to 200 km and 85% of larvae traveled 150 km). The question is, when during larval development, and where on the shallow SW Florida shelf, does the tidal response become incorporated into the behavior of pink shrimp.

Problemas inversos em engenharia financeira: regularização com critério de entropia

Esta dissertação aplica a regularização por entropia máxima no problema inverso de apreçamento de opções, sugerido pelo trabalho de Neri e Schneider em 2012. Eles observaram que a densidade de probabilidade que resolve este problema, no caso de dados provenientes de opções de compra e opções digitais, pode ser descrito como exponenciais nos diferentes intervalos da semireta positiva. Estes intervalos são limitados pelos preços de exercício. O critério de entropia máxima é uma ferramenta poderosa para regularizar este problema mal posto. A família de exponencial do conjunto solução, é calculado usando o algoritmo de Newton-Raphson, com limites específicos para as opções digitais. Estes limites são resultados do princípio de ausência de arbitragem. A metodologia foi usada em dados do índice de ação da Bolsa de Valores de São Paulo com seus preços de opções de compra em diferentes preços de exercício. A análise paramétrica da entropia em função do preços de opções digitais sínteticas (construídas a partir de limites respeitando a ausência de arbitragem) mostraram valores onde as digitais maximizaram a entropia. O exemplo de extração de dados do IBOVESPA de 24 de janeiro de 2013, mostrou um desvio do princípio de ausência de arbitragem para as opções de compra in the money. Este princípio é uma condição necessária para aplicar a regularização por entropia máxima a fim de obter a densidade e os preços. Nossos resultados mostraram que, uma vez preenchida a condição de convexidade na ausência de arbitragem, é possível ter uma forma de smile na curva de volatilidade, com preços calculados a partir da densidade exponencial do modelo. Isto coloca o modelo consistente com os dados do mercado. Do ponto de vista computacional, esta dissertação permitiu de implementar, um modelo de apreçamento que utiliza o princípio de entropia máxima. Três algoritmos clássicos foram usados: primeiramente a bisseção padrão, e depois uma combinação de metodo de bisseção com Newton-Raphson para achar a volatilidade implícita proveniente dos dados de mercado. Depois, o metodo de Newton-Raphson unidimensional para o cálculo dos coeficientes das densidades exponenciais: este é objetivo do estudo. Enfim, o algoritmo de Simpson foi usado para o calculo integral das distribuições cumulativas bem como os preços do modelo obtido através da esperança matemática.

Double layers in gravity theories

Gravitational double layers, unlike their classical electromagnetic counterparts, are thought to be forbidden in gravity theories. It has been recently shown, however, that they are feasible in, for instance, gravity theories with a Lagrangian quadratic in the curvature. This is surprising with many potential consequences and the possibility of new physical behaviours. While a clear interpretation seems elusive, several lines of research are open. I present the field equations for double layers, the new physical quantities arising, and several explicit examples

Avaliação de ligações entre vigas de perfil I ou H e colunas em perfis tubulares circulares.

Considerando uma nova realidade com o incremento do uso de perfis tubulares, este trabalho apresenta uma análise de ligações soldadas tipo T com perfis tubulares circulares (CHS) para a coluna e perfil I ou H para a viga com flexão no plano, efetuado com base nas normas EC3, CIDECT e a NBR 16239 comparando com o critério de deformação limite, proposto por Lu et al., através de um modelo em elementos finitos desenvolvido no programa Ansys versão 12.0. A não-linearidade geométrica foi introduzida no modelo através da Formulação de Lagrange Atualizado. Com base nos resultados numéricos avaliados foram traçadas curvas momento-rotação para cada modelo com o objetivo de obter o momento resistente de cada ligação bem como a classificação da ligação quanto a capacidade de rotação, a influência dos parâmetros geométricos em cada modelo e o modo de falha que controlará o dimensionamento da ligação.

Avaliação de ligações entre perfis tubulares em aço submetidas a esforços de flexão.

Os perfis tubulares sem costura são largamente utilizados em diversos países, porém, no Brasil, o uso desses perfis na construção civil era bastante limitado, restringindo-se praticamente a coberturas espaciais. Considerando uma nova realidade para o uso de perfis tubulares, este trabalho apresenta uma análise de ligações tipo T com perfis tubulares quadrados (SHS) para o banzo e para o montante efetuada com base na norma europeia, Eurocode 3, no CIDECT, na NBR 16239:2013 e ISO 14346 através de um modelo em elementos finitos desenvolvido no programa Ansys. Verificou-se a influência do momento atuante no montante no comportamento global das ligações. As não-linearidades físicas e geométricas foram incorporadas aos modelos, a fim de se mobilizar totalmente a capacidade resistente desta ligação. A não-linearidade do material foi considerada através do critério de plastificação de von Mises através da lei constitutiva tensão versus deformação trilinear de forma a exibir um comportamento elasto-plástico com encruamento. A nãolinearidade geométrica foi introduzida no modelo através da Formulação de Lagrange Atualizado. Os resultados numéricos são avaliados para a ligação em estudo quanto aos modos de falha e a distribuição de tensões. A análise dos momentos resistentes obtidos em comparação com os resultados do modelo numérico, apresentou valores excessivamente a favor da segurança no cálculo utilizando as equações de dimensionamento. Um estudo para fatores de correção das equações de dimensionamento é proposto.

Implementação distribuída de auto-cura em redes inteligentes de distribuição de energia elétrica utilizando árvores de extensão mínima.

A propriedade de auto-cura, em redes inteligente de distribuição de energia elétrica, consiste em encontrar uma proposta de reconfiguração do sistema de distribuição com o objetivo de recuperar parcial ou totalmente o fornecimento de energia aos clientes da rede, na ocorrência de uma falha na rede que comprometa o fornecimento. A busca por uma solução satisfatória é um problema combinacional cuja complexidade está ligada ao tamanho da rede. Um método de busca exaustiva se torna um processo muito demorado e muitas vezes computacionalmente inviável. Para superar essa dificuldade, pode-se basear nas técnicas de geração de árvores de extensão mínima do grafo, representando a rede de distribuição. Porém, a maioria dos estudos encontrados nesta área são implementações centralizadas, onde proposta de reconfiguração é obtida por um sistema de supervisão central. Nesta dissertação, propõe-se uma implementação distribuída, onde cada chave da rede colabora na elaboração da proposta de reconfiguração. A solução descentralizada busca uma redução no tempo de reconfiguração da rede em caso de falhas simples ou múltiplas, aumentando assim a inteligência da rede. Para isso, o algoritmo distribuído GHS é utilizado como base na elaboração de uma solução de auto-cura a ser embarcada nos elementos processadores que compõem as chaves de comutação das linhas da rede inteligente de distribuição. A solução proposta é implementada utilizando robôs como unidades de processamento que se comunicam via uma mesma rede, constituindo assim um ambiente de processamento distribuído. Os diferentes estudos de casos testados mostram que, para redes inteligentes de distribuição compostas por um único alimentador, a solução proposta obteve sucesso na reconfiguração da rede, indiferentemente do número de falhas simultâneas. Na implementação proposta, o tempo de reconfiguração da rede não depende do número de linhas nela incluídas. A implementação apresentou resultados de custo de comunicação e tempo dentro dos limites teóricos estabelecidos pelo algoritmo GHS.

Jato transversal de gotas: simulações por ALE/FEM e efeitos interfaciais.

Um código computacional para escoamentos bifásicos incorporando metodologia híbrida entre oMétodo dos Elementos Finitos e a descrição Lagrangeana-Euleriana Arbitrária do movimento é usado para simular a dinâmica de um jato transversal de gotas na zona primária de quebra. Os corpos dispersos são descritos por meio de um método do tipo front-tracking que produz interfaces de espessura zero através de malhas formadas pela união de elementos adjacentes em ambas as fases e de técnicas de refinamento adaptativo. Condições de contorno periódicas são implementadas de modo variacionalmente consistente para todos os campos envolvidos nas simulações apresentadas e uma versão modificada do campo de pressão é adicionada à formulação do tipo um-fluido usada na equação da quantidade de movimento linear. Simulações numéricas diretas em três dimensões são executadas para diferentes configurações de líquidos imiscí veis compatíveis com resultados experimentais encontrados na literatura. Análises da hidrodinâmica do jato transversal de gotas nessas configurações considerando trajetórias, variação de formato de gota, espectro de pequenas perturbações, além de aspectos complementares relativos à qualidade de malha são apresentados e discutidos.

Uma Metodologia de Estudo de Simulação Tridimensional de Escoamento Turbulento Estratificado no Reservatório de Plantas Hidrelétricas.

Uma simulação numérica que leva em conta os efeitos de estratificação e mistura escalar (como a temperatura, salinidade ou substância solúvel em água) é necessária para estudar e prever os impactos ambientais que um reservatório de usina hidrelétrica pode produzir. Este trabalho sugere uma metodologia para o estudo de escoamentos ambientais, principalmente aqueles em que o conhecimento da interação entre a estratificação e mistura pode dar noções importantes dos fenômenos que ocorrem. Por esta razão, ferramentas de simulação numérica 3D de escoamento ambiental são desenvolvidas. Um gerador de malha de tetraedros do reservatório e o modelo de turbulência algébrico baseado no número de Richardson são as principais ferramentas desenvolvidas. A principal dificuldade na geração de uma malha de tetraedros de um reservatório é a distribuição não uniforme dos pontos relacionada com a relação desproporcional entre as escalas horizontais e verticais do reservatório. Neste tipo de distribuição de pontos, o algoritmo convencional de geração de malha de tetraedros pode tornar-se instável. Por esta razão, um gerador de malha não estruturada de tetraedros é desenvolvido e a metodologia utilizada para obter elementos conformes é descrita. A geração de malha superficial de triângulos utilizando a triangulação Delaunay e a construção do tetraedros a partir da malha triangular são os principais passos para o gerador de malha. A simulação hidrodinâmica com o modelo de turbulência fornece uma ferramenta útil e computacionalmente viável para fins de engenharia. Além disso, o modelo de turbulência baseado no número de Richardson leva em conta os efeitos da interação entre turbulência e estratificação. O modelo algébrico é o mais simples entre os diversos modelos de turbulência. Mas, fornece resultados realistas com o ajuste de uma pequena quantidade de parâmetros. São incorporados os modelos de viscosidade/difusividade turbulenta para escoamento estratificado. Na aproximação das equações médias de Reynolds e transporte de escalar é utilizando o Método dos Elementos Finitos. Os termos convectivos são aproximados utilizando o método semi-Lagrangeano, e a aproximação espacial é baseada no método de Galerkin. Os resultados computacionais são comparados com os resultados disponíveis na literatura. E, finalmente, a simulação de escoamento em um braço de reservatório é apresentada.

Elementos finitos em fluidos dominados pelo fenômeno de advecção: um método semi-Lagrangeano.

Os escoamentos altamente convectivos representam um desafio na simulação pelo método de elementos finitos. Com a solução de elementos finitos de Galerkin para escoamentos incompressíveis, a matriz associada ao termo convectivo é não simétrica, e portanto, a propiedade de aproximação ótima é perdida. Na prática as soluções apresentam oscilações espúrias. Muitos métodos foram desenvolvidos com o fim de resolver esse problema. Neste trabalho apresentamos um método semi- Lagrangeano, o qual é implicitamente um método do tipo upwind, que portanto resolve o problema anterior, e comparamos o desempenho do método na solução das equações de convecção-difusão e Navier-Stokes incompressível com o Streamline Upwind Petrov Galerkin (SUPG), um método estabilizador de reconhecido desempenho. No SUPG, as funções de forma e de teste são tomadas em espaços diferentes, criando um efeito tal que as oscilações espúrias são drasticamente atenuadas. O método semi-Lagrangeano é um método de fator de integração, no qual o fator é um operador de convecção que se desloca para um sistema de coordenadas móveis no fluido, mas restabelece o sistema de coordenadas Lagrangeanas depois de cada passo de tempo. Isto prevê estabilidade e a possibilidade de utilizar passos de tempo maiores.Existem muitos trabalhos na literatura analisando métodos estabilizadores, mas não assim com o método semi-Lagrangeano, o que representa a contribuição principal deste trabalho: reconhecer as virtudes e as fraquezas do método semi-Lagrangeano em escoamentos dominados pelo fenômeno de convecção.

Comportamento estrutural de vigas casteladas em aço inoxidável.

As vigas casteladas já são utilizadas em diversos tipos de estruturas para se vencer grandes vãos. Há uma nova realidade para o uso de perfis castelados, agora em aço inoxidável. Este trabalho apresenta uma análise comparativa entre vigas casteladas e vigas de alma cheia em perfis IPE, baseadas na norma europeia, no método da resistência contínua; e, também, através de um modelo em elementos finitos desenvolvido no programa Ansys. Este trabalho tem como objetivo verificar o comportamento e a resistência à flexão das vigas casteladas em aço inoxidável. As não-linearidades físicas e geométricas foram incorporadas aos modelos, a fim de se mobilizar totalmente a capacidade resistente desta estrutura. A não-linearidade do material foi considerada através do critério de plastificação de von Mises e da lei constitutiva tensão versus deformação multi-linear, de forma a exibir um comportamento elasto-plástico com encruamento. A não-linearidade geométrica foi considerada a partir de uma formulação de Lagrange atualizado. Os resultados numéricos das vigas em estudo foram avaliados quanto aos modos de falha e da distribuição de tensões. Os momentos resistentes analíticos foram comparados com os resultados do modelo numérico apresentando valores satisfatórios e a favor da segurança.

Avaliação de ligações tubulares tipo T entre CHS com reforço tipo chapa.

Os perfis tubulares sem costura são utilizados em diversos países, principalmente devido às vantagens associadas à estética a sua elevada resistência à torção, cargas axiais e efeitos combinados. Canadá, Inglaterra, Alemanha e Holanda utilizam de forma veemente estas estruturas e possuem produção contínua e industrializada com alto nível de desenvolvimento tecnológico. O Brasil, porém, se limitava praticamente ao uso destes perfis nas coberturas espaciais. Devido ao aumento da utilização desses tipos de estruturas, fez-se necessário o aprofundamento dos estudos com métodos de análise coerentes para utilização de perfis tubulares, principalmente em relação às ligações, pois são consideradas regiões vulneráveis neste tipo de estrutura. Para atender a necessidade de normatização deste procedimento desenvolveu-se uma norma brasileira específica para o dimensionamento de estruturas em perfis tubulares. Considerando esta perspectiva, este trabalho apresenta uma análise de ligações tipo T com reforço tipo chapa com perfis tubulares circulares (CHS) para o banzo e para o montante efetuada com base na norma europeia, Eurocode 3, no CIDECT, na NBR 16239:2013 e ISO 14346. Desenvolveu-se no programa Ansys um modelo de elementos finitos para cada tipo de ligação analisada, calibrado e validado com resultados experimentais e numéricos existentes na literatura. Verificou-se a influência da compressão atuante no montante no comportamento global das ligações. As não-linearidades físicas e geométricas foram incorporadas aos modelos a fim de se mobilizar totalmente a capacidade resistente desta ligação. A nãolinearidade do material foi considerada através do critério de plastificação de von Mises através da lei constitutiva tensão versus deformação bilinear de forma a exibir um comportamento elasto-plástico com encruamento. A não-linearidade geométrica foi introduzida no modelo através da Formulação de Lagrange Atualizada. A análise dos esforços resistentes obtidos em comparação com os resultados do modelo numérico, apresentaram valores a favor da segurança no cálculo utilizando as equações de dimensionamento. Por fim um estudo para fatores de correção das equações de dimensionamento foi também proposto.

DEVELOPMENT OF THE WAVE CONTOURING RAFT.

The wave contouring raft system is the outcome of ideas initiated and developed by Sir Christopher Cockerell from 1972 onwards. His objective was to develop a wave energy device which is within the bounds of current technology. It should consist of simple, relatively small units, amenable to quantity production, which would enable a power generating system to be built up and commissioned in stages according to needs and production capability. This thinking led to the investigation of chains of pontoons, hinged together so that the passage of a wave down the chain causes the pontoons to oscillate relative to one another. Energy is extracted from the sea by applying a torque about the hinges to damp the motion. The work has involved extensive model testing in wave tanks and the building and testing of a 3-unit 1/10 scale power generating installation in the Solent, as well as design studies for a full size installation for Atlantic conditions.

The effect of surface hardening on the elastic shakedown of elliptical contacts

The effect of varying both the aspect ratio and the coefficient of friction of contacts with elliptical geometry on their elastic shakedown performance has been examined theoretically for surfaces with two types of subsurface hardness or strength profiles. In stepwise hardening the hard layer is of uniform strength while in linear hardening its strength reduces from a maximum at the surface to that of the core at the base of the hardened layer. The shakedown load is expressed as the ratio of the maximum Hertzian pressure to the strength of the core material. As the depth of hardening, expressed as a multiple of the elliptical semi-axis, is increased so the potential shakedown load increases from a level that is appropriate to a uniform half-space of unhardened material to a value reflecting the hardness of the surface and near-surface material. In a step-hardened material, the shakedown limit for a surface 'pummelled' by the passage of a sequence of such loads reaches a cut-off or plateau value, which cannot be exceeded by further increases in hardening depth irrespective of the value of the friction coefficient. For a linear-hardened material the corresponding plateau is approached asymptotically. The work confirms earlier results on the upper bounds on shakedown of both point and line contacts and provides numerical values of shakedown loads for intermediate geometries. In general, the case depth required to achieve a given shakedown limit reduces in moving from a transversely moving nominal line load to an axisymmetric point load.