Full and Model-Reduced Structure-Preserving Simulation of Incompressible Fluids

This thesis outlines the construction of several types of structured integrators for incompressible fluids. We first present a vorticity integrator, which is the Hamiltonian counterpart of the existing Lagrangian-based fluid integrator. We next present a model-reduced variational Eulerian integrator for incompressible fluids, which combines the efficiency gains of dimension reduction, the qualitative robustness to coarse spatial and temporal resolutions of geometric integrators, and the simplicity of homogenized boundary conditions on regular grids to deal with arbitrarily-shaped domains with sub-grid accuracy.

Both these numerical methods involve approximating the Lie group of volume-preserving diffeomorphisms by a finite-dimensional Lie-group and then restricting the resulting variational principle by means of a non-holonomic constraint. Advantages and limitations of this discretization method will be outlined. It will be seen that these derivation techniques are unable to yield symplectic integrators, but that energy conservation is easily obtained, as is a discretized version of Kelvin's circulation theorem.

Finally, we outline the basis of a spectral discrete exterior calculus, which may be a useful element in producing structured numerical methods for fluids in the future.

Tratamento Homeopático: buscando estratégias para avaliação de seus resultados e efetividade

Face à crescente procura por outras modalidades terapêuticas que abordam o ser humano de forma holística e a introdução das mesmas no SUS, torna-se muito importante a avaliação da efetividade e segurança dessas formas de cuidado. A Homeopatia faz parte deste conjunto de terapêuticas e, para sua avaliação, pode existir a necessidade de se valer de múltiplos instrumentos para abarcar os vários aspectos de uma resposta integral ao tratamento. Este trabalho teve como objetivo identificar e elaborar categorias de análises e instrumentos que permitam avaliar e mensurar a efetividade deste tratamento, bem como testá-los, considerando-se as características desta racionalidade. Foram levantados, na literatura nacional e internacional, trabalhos sobre efetividade do tratamento homeopático, em busca da definição do estado da arte mas também dos principais problemas, limitações e possibilidades dessas avaliações tendo em vista seu resultado integral. Finda esta etapa, a pesquisa destinou-se a elaboração, proposição e testagem de uma metodologia considerada mais adequada a avaliar o tratamento homeopático nesta perspectiva. Um estudo observacional foi realizado em serviço público homeopático no município de Juiz de Fora, com tratamento individualizado, no qual foi utilizada uma estratégia de avaliação composta por três componentes: (1) avaliação de qualidade de vida pelo instrumento SF-36; (2) análises em busca de objetivar e quantificar queixas clínicas e outros atributos de natureza subjetiva (sensação de bem-estar, sono, estado cognitivo e memória, vida sexual, sensação de felicidade) por meio da utilização de uma escala visual analógica (EVA), na mensuração da intensidade e de opções fechadas, a exemplo do SF-36, na estimativa da frequência desses aspectos e (3) entrevistas qualitativas por intermédio de questionário semiestruturado, com a finalidade de abordar questões relacionadas a biopatografia e mudança da atitude vital (como pacientes enfrentam os problemas do cotidiano, fatores deflagradores das queixas, como se sentem e como reagem, além de indagar seus projetos de vida e felicidade). A aplicação do questionário SF-36 apresentou algumas dificuldades de compreensão pelos participantes, talvez devido à baixa escolaridade dos entrevistados, mas mostrou-se útil à pesquisa, embora demonstre limitações na avaliação do aspecto integral do resultado da terapêutica analisada. O acompanhamento das queixas clínicas, sensação de bem-estar, sono e estado cognitivo e memória foram captados e mensurados de forma satisfatória tanto pela EVA (intensidade dos sintomas) quanto pelas respostas fechadas para medir a frequência. Situações como as avaliações da biopatografia e da sexualidade foram insuficientes para serem adequadamente avaliadas pelo pesquisador e o paciente somente. A participação do médico assistente poderia contribuir nestes casos. Questões mais abrangentes na avaliação da mudança na atitude vital, como reação diante de fatores desequilibrantes e projeto de vida e felicidade, necessitam de metodologia qualitativa até que se possa avançar nas pesquisas à espera de soluções futuras. A combinação dessas estratégias em estudos controlados, randomizados, com amostras de magnitude satisfatória, preferencialmente em rede e que explicitem as condições nas quais o atendimento homeopático ocorreu e como se chegou a cada prescrição, podem ter utilidade para a avaliação da efetividade da dimensão integral do tratamento homeopático.

Microscopic Origin of Macroscopic Strength in Granular Media: A Numerical and Analytical Approach

Constitutive modeling in granular materials has historically been based on macroscopic experimental observations that, while being usually effective at predicting the bulk behavior of these type of materials, suffer important limitations when it comes to understanding the physics behind grain-to-grain interactions that induce the material to macroscopically behave in a given way when subjected to certain boundary conditions.

The advent of the discrete element method (DEM) in the late 1970s helped scientists and engineers to gain a deeper insight into some of the most fundamental mechanisms furnishing the grain scale. However, one of the most critical limitations of classical DEM schemes has been their inability to account for complex grain morphologies. Instead, simplified geometries such as discs, spheres, and polyhedra have typically been used. Fortunately, in the last fifteen years, there has been an increasing development of new computational as well as experimental techniques, such as non-uniform rational basis splines (NURBS) and 3D X-ray Computed Tomography (3DXRCT), which are contributing to create new tools that enable the inclusion of complex grain morphologies into DEM schemes.

Yet, as the scientific community is still developing these new tools, there is still a gap in thoroughly understanding the physical relations connecting grain and continuum scales as well as in the development of discrete techniques that can predict the emergent behavior of granular materials without resorting to phenomenology, but rather can directly unravel the micro-mechanical origin of macroscopic behavior.

In order to contribute towards closing the aforementioned gap, we have developed a micro-mechanical analysis of macroscopic peak strength, critical state, and residual strength in two-dimensional non-cohesive granular media, where typical continuum constitutive quantities such as frictional strength and dilation angle are explicitly related to their corresponding grain-scale counterparts (e.g., inter-particle contact forces, fabric, particle displacements, and velocities), providing an across-the-scale basis for better understanding and modeling granular media.

In the same way, we utilize a new DEM scheme (LS-DEM) that takes advantage of a mathematical technique called level set (LS) to enable the inclusion of real grain shapes into a classical discrete element method. After calibrating LS-DEM with respect to real experimental results, we exploit part of its potential to study the dependency of critical state (CS) parameters such as the critical state line (CSL) slope, CSL intercept, and CS friction angle on the grain's morphology, i.e., sphericity, roundness, and regularity.

Finally, we introduce a first computational algorithm to ``clone'' the grain morphologies of a sample of real digital grains. This cloning algorithm allows us to generate an arbitrary number of cloned grains that satisfy the same morphological features (e.g., roundness and aspect ratio) displayed by their real parents and can be included into a DEM simulation of a given mechanical phenomenon. In turn, this will help with the development of discrete techniques that can directly predict the engineering scale behavior of granular media without resorting to phenomenology.

Kires: a data-centric telerehabilitation system based on kinect

185 p.

A working guide to the assessment, implementation and post project monitoring of fisheries habitat improvement schemes in South Cumbria

Due to changes in land use over the last century, the physical nature of many streams and rivers in the British Isles has probably changed. In some cases this change may be large for example as a result of flood defence schemes and is easily observed, whilst in other cases altered land use, farming, forestry or urbanization may have resulted in more subtle changes to river features. This working guide draws together a way of assessing habitat in any stream or river and determine sites or reaches on the assessed watercourse that may benefit from habitat improvement schemes. It will determine a method of measuring existing habitat in a broad sense, whilst referring to R and D studies currently being undertaken in this area. A method of prioritising any proposed habitat restoration work will be suggested. The limitations of fisheries improvement schemes in terms of cross functional acceptance (flood defence and conservation) will be examined along with suggested proposals for some example watercourses. The need for pre and post enhancement monitoring will be discussed as will the requirement for maintenance programs on schemes. Finally methods for determining the cost benefits of small schemes will be examined, compared to other currently used enhancement strategies. This will allow small scale revenue schemes to be used to back up pre project cost benefit analysis as required in future capital submissions.

Pacific Salmon, Oncorhynchus spp., and the Definition of "Species" Under the Endangered Species Act

For purposes ofthe Endangered Species Act (ESA), a "species" is defined to include "any distinct population segment of any species of vertebrate fish or wildlife which interbreeds when mature. "Federal agencies charged with carrying out the provisions of the ESA have struggled for over a decade to develop a consistent approach for interpreting the term "distinct population segment." This paper outlines such an approach and explains in some detail how it can be applied to ESA evaluations of anadromous Pacific salmonids. The following definition is proposed: A population (or group of populations) will be considered "distinct" (and hence a "species ")for purposes of the ESA if it represents an evolutionarily significant unit (ESU) of the biological species. A population must satisfy two criteria to be considered an ESU: 1) It must be substantially reproductively isolated from other conspecific population units, and 2) It must represent an important component in the evolutionary legacy of the species. Isolation does not have to be absolute, but it must be strong enough to permit evolutionarily important differences to accrue in different population units. The second criterion would be met if the population contributes substantially to the ecological/genetic diversity of the species as a whole. Insights into the extent of reproductive isolation can be provided by movements of tagged fish, natural recolonization rates observed in other populations, measurements of genetic differences between populations, and evaluations of the efficacy of natural barriers. Each of these methods has its limitations. Identification of physical barriers to genetic exchange can help define the geographic extent of distinct populations, but reliance on physical features alone can be misleading in the absence of supporting biological information. Physical tags provide information about the movements of individual fish but not the genetic consequences of migration. Furthermore, measurements ofc urrent straying or recolonization rates provide no direct information about the magnitude or consistency of such rates in the past. In this respect, data from protein electrophoresis or DNA analyses can be very useful because they reflect levels of gene flow that have occurred over evolutionary time scales. The best strategy is to use all available lines of evidence for or against reproductive isolation, recognizing the limitations of each and taking advantage of the often complementary nature of the different types of information. If available evidence indicates significant reproductive isolation, the next step is to determine whether the population in question is of substantial ecological/genetic importance to the species as a whole. In other words, if the population became extinct, would this event represent a significant loss to the ecological/genetic diversity of thes pecies? In making this determination, the following questions are relevant: 1) Is the population genetically distinct from other conspecific populations? 2) Does the population occupy unusual or distinctive habitat? 3) Does the population show evidence of unusual or distinctive adaptation to its environment? Several types of information are useful in addressing these questions. Again, the strengths and limitations of each should be kept in mind in making the evaluation. Phenotypic/life-history traits such as size, fecundity, and age and time of spawning may reflect local adaptations of evolutionary importance, but interpretation of these traits is complicated by their sensitivity to environmental conditions. Data from protein electrophoresis or DNA analyses provide valuable insight into theprocessofgenetic differentiation among populations but little direct information regarding the extent of adaptive genetic differences. Habitat differences suggest the possibility for local adaptations but do not prove that such adaptations exist. The framework suggested here provides a focal point for accomplishing the majorgoal of the Act-to conserve the genetic diversity of species and the ecosystems they inhabit. At the same time, it allows discretion in the listing of populations by requiring that they represent units of real evolutionary significance to the species. Further, this framework provides a means of addressing several issues of particular concern for Pacific salmon, including anadromous/nonanadromous population segments, differences in run-timing, groups of populations, introduced populations, and the role of hatchery fish.

Caracterização de Bacilos Gram-Negativos Não Fermentadores não usuais em bacteremias pelas técnicas de Matrix-Assisted Laser Desorption IonizationTime of Flight Mass Spectrometry, sequenciamento de DNA e método fenotípico convencional

Alguns Bastonetes Gram-negativos não fermentadores (BGNNF) costumam ser considerados clinicamente pouco significantes e a sua implicação em infecções é subestimada. Devido à similaridade fenotípica, mudanças taxonômicas, baixa reatividade bioquímica e limitações nos bancos de dados em sistemas comerciais, a identificação de BGNNF é frequentemente equivocada, culminando com a denominação de diferentes micro-organismos apenas como BGNNF, por falta de melhor diferenciação. O objetivo desse estudo foi avaliar, por métodos fenotípico convencional, proteômico e molecular, a identificação de BGNNF incomuns isolados em hemoculturas de pacientes atendidos em um hospital universitário no Rio de Janeiro. Foram selecionadas 78 amostras isoladas de hemoculturas caracterizadas no laboratório clinico como BGNNF para a identificação por sequenciamento dos genes 16S RNA e recA, por um conjunto amplo de testes fenotípicos manuais e por MALDI-TOF MS. Os micro-organismos predominantes na amostragem foram genotipados pela técnica de eletroforese em gel de campo pulsado (PFGE). Pelo sequenciamento do gene 16S rRNA, a maioria das amostras (n=31; 40%) foi incluída no gênero Burkholderia, seguido de Pseudomonas stutzeri (10%) e Delftia acidovorans (4%). Os demais isolados foram agrupados em 27 diferentes espécies. O sequencimento do gene recA identificou a maioria das espécies de Burkholderia como Burkholderia contaminans (n=19; 24%). Os testes fenotípicos incluíram as 31 amostras apenas no CBc e para as outras 47 amostras, a concordância com o sequenciamento do gene 16S rRNA em nível de espécie foi de 64% (n=30) e apenas em gênero a concordância foi de 17% (n=8). A análise comparativa geral da identificação por MALDI-TOF MS com o sequenciamento do gene16S rRNA mostrou que 42% (n=33) das 78 amostras foram concordantes em nível de espécie e 45% (n=35) apenas em gênero. Excluindo as amostras do CBc, houve um aumento da concordância em nível de espécie para 60%. As discordâncias parecem ser devido às diferenças nos perfis proteicos das amostras em relação às amostras-referência do banco de dados do equipamento e podem ser aprimorados com a atualização de perfis no sistema. A análise do polimorfismo genético de B. contaminans mostrou a ausência de um clone disseminado causando surto, além da provável origem ambiental das infecções. Os setores de nefrologia e hemodiálise contribuíram com maior número de pacientes com amostras positivas (5 pacientes e 9 amostras). Os grupos clonais BcoD e BcoE foram encontrados em pacientes assistidos no mesmo setor com diferença de quatro meses (BcoD, nefrologia) e 1,5 ano (BcoE, hemodilálise), entre as culturas, respectivamente. As discordâncias entre as técnicas ocorreram principalmente devido a dificuldade de identificação das espécies do CBc. Os BGNNF incomuns são de difícil caracterização independente da metodologia usada e nenhum método por si só foi capaz de identificar todas as amostras.

A stochastic model of temporal variations in monthly temperature, precipitation, snowfall, and resulting snowpack

We report a Monte Carlo representation of the long-term inter-annual variability of monthly snowfall on a detailed (1 km) grid of points throughout the southwest. An extension of the local climate model of the southwestern United States (Stamm and Craig 1992) provides spatially based estimates of mean and variance of monthly temperature and precipitation. The mean is the expected value from a canonical regression using independent variables that represent controls on climate in this area, including orography. Variance is computed as the standard error of the prediction and provides site-specific measures of (1) natural sources of variation and (2) errors due to limitations of the data and poor distribution of climate stations. Simulation of monthly temperature and precipitation over a sequence of years is achieved by drawing from a bivariate normal distribution. The conditional expectation of precipitation. given temperature in each month, is the basis of a numerical integration of the normal probability distribution of log precipitation below a threshold temperature (3°C) to determine snowfall as a percent of total precipitation. Snowfall predictions are tested at stations for which long-term records are available. At Donner Memorial State Park (elevation 1811 meters) a 34-year simulation - matching the length of instrumental record - is within 15 percent of observed for mean annual snowfall. We also compute resulting snowpack using a variation of the model of Martinec et al. (1983). This allows additional tests by examining spatial patterns of predicted snowfall and snowpack and their hydrologic implications.

Ecological and economic consequences of hypoxia: Topic 2 Report for the Integrated Assessment on Hypoxia in the Gulf of Mexico

In this report we have attempted to evaluate the ecological and economic consequences of hypoxia in the northern Gulf of Mexico. Although our initial approach was to rely on published accounts, we quickly realized that the body of published literature deahng with hypoxia was limited, and we would have to conduct our own exploratory analysis of existing Gulf data, or rely on published accounts from other systems to infer possible or potential effects of hypoxia. For the economic analysis, we developed a conceptual model of how hypoxia-related impacts could affect fisheries. Our model included both supply and demand components. The supply model had two components: (1) a physical production function for fish or shrimp, and (2) the cost of fishing. If hypoxia causes the cost of a unit of fishing effort to change, then this will result in a shift in supply. The demand model considered how hypoxia might affect the quality of landed fish or shrimp. In particular, the market value per pound is lower for small shrimp than for large shrimp. Given the limitations of the ecological assessment, the shallow continental shelf area affected by hypoxia does show signs of hypoxia-related stress. While current ecological conditions are a response to a variety of stressors, the effects of hypoxia are most obvious in the benthos that experience mortality, elimination of larger long-lived species, and a shifting of productivity to nonhypoxic periods (energy pulsing). What is not known is whether hypoxia leads to higher productivity during productive periods, or simply to a reduction of productivity during oxygen-stressed periods. The economic assessment based on fisheries data, however, failed to detect effects attributable to hypoxia. Overall, fisheries landings statistics for at least the last few decades have been relatively constant. The failure to identify clear hypoxic effects in the fisheries statistics does not necessarily mean that they are absent. There are several possibilities: (1) hypoxic effects are small relative to the overall variability in the data sets evaluated; (2) the data and the power of the analyses are not adequate; and (3) currently there are no hypoxic effects on fisheries. Lack of identified hypoxic effects in available fisheries data does not imply that effects would not occur should conditions worsen. Experience with other hypoxic zones around the globe shows that both ecological and fisheries effects become progressively more severe as hypoxia increases. Several large systems around the globe have suffered serious ecological and economic consequences from seasonal summertime hypoxia; most notable are the Kattegat and Black Sea. The consequences range from localized loss of catch and recruitment failure to complete system-wide loss of fishery species. If experiences in other systems are applicable to the Gulf of Mexico, then in the face of worsening hypoxic conditions, at some point fisheries and other species will decline, perhaps precipitously.

Spiny lobsters of Ceylon

Spiny lobsters of the family Palinuridae are known to exist in parts of Asia such as India, Japan, Indonesia and Malaya. In Ceylon spiny lobsters are caught chiefly by a primitive type of gear designed to catch these species and incidentally by nets set to catch fish. The limitations of these types of gear and the results of experimental fishing in Ceylon with different designs of lobster traps have been discussed earlier. As a result of these experiments a very effective trap for capturing spiny lobsters was developed and it was shown that spiny lobsters were present in large quantities (De Bruin 1960). This created great interest among skin-divers who found lobster fishing very remunerative.

Modeling of unsteady transitional boundary layers

In turbomachinery, a considerable proportion of the blade surface area can be covered by transitional boundary layers. This means that accurate prediction of the profile loss and boundary layer behavior in general depends on the accurate modeling of the transitional boundary layers, especially at low Reynolds numbers. This paper presents a model for determining the intermittency resulting from the unsteady transition caused by the passage of wakes over a blade surface. The model is founded on work by Emmons (1951) who showed that the intermittency could be calculated from a knowledge of the behavior of randomly formed turbulent spots. The model is used to calculate the development of the boundary layer on the rotor of a low Reynolds number single-stage turbine. The predictions are compared with experimental results obtained using surface-mounted hot-film anemometers and hot-wire traverses of the rotor midspan boundary layer at two different rotor-stator gaps. The validity and limitations of the model are discussed.

Optimization of semiconductor optical amplifiers for all optical wavelength conversion

In recent years a variety of experimental and theoretical work has been reported on the use of semiconductor optical amplifiers for high speed wavelength conversion. However little work has addressed the dynamic limitations of this conversion process in detail with a view to device optimization. In this paper, a detailed study of the conversion process is carried out in order to optimize device parameters and drive conditions for increased conversion speed and improved modulation index.

Growth of Skipjack Katsuwonus pelamis from Mozambique

Growth of skipjack, Katsuwonus pelamis in Mozambican waters was studied, using fork length data collected in 1983-1985, during an experimental pole-and-line/live-bait fishing program. Tentative estimates of growth parameters are FL∞=80 cm, K=0.6 yearˉ¹. Some evidence for seasonal growth is presented, along with a discussion of the limitations of the data set used in the investigation.

An overview of Sequential Monte Carlo methods for parameter estimation in general state-space models

Nonlinear non-Gaussian state-space models arise in numerous applications in control and signal processing. Sequential Monte Carlo (SMC) methods, also known as Particle Filters, provide very good numerical approximations to the associated optimal state estimation problems. However, in many scenarios, the state-space model of interest also depends on unknown static parameters that need to be estimated from the data. In this context, standard SMC methods fail and it is necessary to rely on more sophisticated algorithms. The aim of this paper is to present a comprehensive overview of SMC methods that have been proposed to perform static parameter estimation in general state-space models. We discuss the advantages and limitations of these methods. © 2009 IFAC.

Vibration of fluid loaded conical shells.

An analytical model is presented to describe the vibration of a truncated conical shell with fluid loading in the low frequency range. The solution for the dynamic response of the shell is presented in the form of a power series. Fluid loading is taken into account by dividing the shell into narrow strips which are considered to be locally cylindrical. Analytical results are presented for different boundary conditions and have been compared with the computational results from a boundary element model. Limitations of the model to the low frequency range are discussed.