Efficient fully nonlinear data assimilation for geophysical fluid dynamics

A potential problem with Ensemble Kalman Filter is the implicit Gaussian assumption at analysis times. Here we explore the performance of a recently proposed fully nonlinear particle filter on a high-dimensional but simplified ocean model, in which the Gaussian assumption is not made. The model simulates the evolution of the vorticity field in time, described by the barotropic vorticity equation, in a highly nonlinear flow regime. While common knowledge is that particle filters are inefficient and need large numbers of model runs to avoid degeneracy, the newly developed particle filter needs only of the order of 10-100 particles on large scale problems. The crucial new ingredient is that the proposal density cannot only be used to ensure all particles end up in high-probability regions of state space as defined by the observations, but also to ensure that most of the particles have similar weights. Using identical twin experiments we found that the ensemble mean follows the truth reliably, and the difference from the truth is captured by the ensemble spread. A rank histogram is used to show that the truth run is indistinguishable from any of the particles, showing statistical consistency of the method.

On the validity of single-parcel energetics to assess the importance of internal energy and compressibility effects in stratified fluids

It is often assumed on the basis of single-parcel energetics that compressible effects and conversions with internal energy are negligible whenever typical displacements of fluid parcels are small relative to the scale height of the fluid (defined as the ratio of the squared speed of sound over gravitational acceleration). This paper shows that the above approach is flawed, however, and that a correct assessment of compressible effects and internal energy conversions requires considering the energetics of at least two parcels, or more generally, of mass conserving parcel re-arrangements. As a consequence, it is shown that it is the adiabatic lapse rate and its derivative with respect to pressure, rather than the scale height, which controls the relative importance of compressible effects and internal energy conversions when considering the global energy budget of a stratied fluid. Only when mass conservation is properly accounted for is it possible to explain why available internal energy can account for up to 40 percent of the total available potential energy in the oceans. This is considerably larger than the prediction of single-parcel energetics, according to which this number should be no more than about 2 percent.

Osteocyte Density in the Peri-Implant Bone of Immediately Loaded and Submerged Dental Implants

Background: The role of osteocytes in bone structure and function remains partially unresolved. Their participation in mechanotransduction, i.e., the conversion of a physical stimulus into a cellular response, has been hypothesized. The present study was an evaluation of the osteocyte density in the peri-implant bone of immediately loaded and submerged dental implants. Methods: Fourteen male patients were included in the study; all of them were partially edentulous and needed a posterior mandibular restoration. Implants were inserted in these areas; half of the sample was loaded immediately (included in a fixed provisional prosthesis on the same day as implant surgery), whereas the other half was left to heal submerged. Fourteen implants (seven immediately loaded and seven unloaded) were retrieved with a trephine after a healing period of 8 weeks. The specimens were treated to obtain thin ground sections, and histomorphometry was used to evaluate the osteocyte index in the peri-implant bone. Results: A higher and statistically significant number of osteocytes was found in the peri-implant bone around immediately loaded implants (P=0.0081). A correlation between the percentage of bone-implant contact and osteocyte density was found for immediately loaded implants (P=0.0480) but not for submerged implants (P=0.2667). Conclusion: The higher number of osteocytes in the peri-implant bone around immediately loaded implants could be related to the functional adaptation required by the loading stimulus, which also explains the hypothesized involvement of the osteocytes in the maintenance of the bone matrix. J Periodontol 2009;80:499-504.

P-T-Fluid evolution and graphite deposition during retrograde metamorphism in Ribeira Fold Belt, SE Brazil: Oxygen fugacity, fluid inclusions and C-O-H isotopic evidence

Combined fluid inclusion (FI) microthermometry, Raman spectroscopy, X-ray diffraction, C-O-H isotopes and oxygen fugacities of granulites from central Ribeira Fold Belt, SE Brazil, provided the following results: i) Magnetite-Hematite fO(2) estimates range from 10(-11.5) bar (QFM + 1) to 10(-18.3) bar (QFM - 1) for the temperature range of 896 degrees C-656 degrees C, implying fO(2) decrease from metamorphic peak temperatures to retrograde conditions; ii) 5 main types of fluid inclusions were observed: a) CO(2) and CO(2)-N(2) (0-11 mol%) high to medium density (1.01-0.59 g/cm(3)) FI; b) CO(2) and CO(2)-N(2) (0-36 mol%) low density (0.19-0.29 g/cm(3)) FI; c) CO(2) (94-95 mol%)-N(2) (3 mol%)-CH(4) (2-3 mol%)-H(2)O (water phi(v) (25 degrees C) = 0.1) FI; d) low-salinity H(2)O-CO(2) FI; and e) late low-salinity H(2)O FI; iii) Raman analyses evidence two graphite types in khondalites: an early highly ordered graphite (T similar to 450 degrees C) overgrown by a disordered kind (T similar to 330 degrees C); iv) delta(18)O quartz results of 10.3-10.7 parts per thousand, imply high-temperature CO(2) delta(18)O values of 14.4-14.8 parts per thousand, suggesting the involvement of a metamorphic fluid, whereas lower temperature biotite delta(18)O and delta D results of 7.5-8.5 parts per thousand and -54 to -67 parts per thousand respectively imply H(2)O delta(18)O values of 10-11 parts per thousand and delta D(H2O) of -23 to -36 parts per thousand suggesting delta(18)O depletion and increasing fluid/rock ratio from metamorphic peak to retrograde conditions. Isotopic results are compatible with low-temperature H(2)O influx and fO(2) decrease that promoted graphite deposition in retrograde granulites, simultaneous with low density CO(2), CO(2)-N(2) and CO(2)-N(2)-CH(4)-H(2)O fluid inclusions at T = 450-330 degrees C. Graphite delta(13)C results of -10.9 to -11.4 parts per thousand imply CO(2) delta(13)C values of -0.8 to -1.3 parts per thousand suggesting decarbonation of Cambrian marine carbonates with small admixture of lighter biogenic or mantle derived fluids. Based on these results, it is suggested that metamorphic fluids from the central segment of Ribeira Fold Belt evolved to CO(2)-N(2) fluids during granulitic metamorphism at high fO(2), followed by rapid pressure drop at T similar to 400-450 degrees C during late exhumation that caused fO(2) reduction induced by temperature decrease and water influx, turning carbonic fluids into CO(2)-H(2)O (depleting biotite delta(18)O and delta D values), and progressively into H(2)O. When fO(2) decreased substantially by mixture of carbonic and aqueous fluids, graphite deposited forming khondalites. (C) 2010 Elsevier Ltd. All rights reserved.

Three-dimensional, fully adaptive simulations of phase-field fluid models

We present an efficient numerical methodology for the 31) computation of incompressible multi-phase flows described by conservative phase-field models We focus here on the case of density matched fluids with different viscosity (Model H) The numerical method employs adaptive mesh refinements (AMR) in concert with an efficient semi-implicit time discretization strategy and a linear, multi-level multigrid to relax high order stability constraints and to capture the flow`s disparate scales at optimal cost. Only five linear solvers are needed per time-step. Moreover, all the adaptive methodology is constructed from scratch to allow a systematic investigation of the key aspects of AMR in a conservative, phase-field setting. We validate the method and demonstrate its capabilities and efficacy with important examples of drop deformation, Kelvin-Helmholtz instability, and flow-induced drop coalescence (C) 2010 Elsevier Inc. All rights reserved

Desenvolvimento de um programa de simulação computacional de sistemas de aquecimento solar para água

Esta Tese apresenta uma análise do comportamento térmico de um sistema de aquecimento solar operando por termossifão. Neste tipo de sistema o fluido no coletor solar é circulado por convecção natural, que acontece devido à diferença de massa específica da água ao longo circuito. Nestes sistemas a vazão mássica varia ao longo do dia e do ano, dependendo, dentre outros fatores, da irradiância solar absorvida, do perfil de temperaturas da água no sistema, da geometria, do volume e do perfil de demanda de água quente. Para uma avaliação detalhada do comportamento térmico de aquecedores solares operando por termossifão foram realizados ensaios experimentais e cálculos teóricos. Os resultados dos experimentos concordaram com aqueles apresentados na literatura e sua análise fundamentou o desenvolvimento do aplicativo TermoSim, um programa de simulação computacional do comportamento térmico de sistemas de aquecimento de água com energia solar. O tratamento matemático adotado no TermoSim compreende a modelagem dos coletores solares de acordo com a teoria de Hottel-Bliss-Whillier. O reservatório térmico é modelado com estratificação térmica, convecção e condução entre as camadas. A vazão mássica é obtida a partir do balanço da quantidade de movimento no circuito. Os modelos matemáticos empregados na construção do aplicativo TermoSim foram validados através do confronto dos resultados simulados com medidas experimentais. Foi demonstrado que a utilização destes modelos é adequada e permite reproduzir com precisão o comportamento térmico dos coletores solares e do reservatório térmico. Além do programa TermoSim, foi também desenvolvido o programa TermoDim, que é uma ferramenta para o dimensionamento de sistemas de aquecimento solar, que requer apenas o conhecimento dos parâmetros geométricos do sistema, dados meteorológicos em média mensal e informação a respeito do volume de demanda. O TermoDim é apropriado para estimar o desempenho de aquecedores solares operando por termossifão com tanques verticais e horizontais. O método de dimensionamento do TermoDim é baseado na correlação para a eficiência média mensal obtida neste trabalho a partir de um grande número de simulações.

Physical-chemical analysis of the cerebrospinal fluid of healthy dogs submitted to different storage periods and temperatures

Disfunções envolvendo o sistema nervoso são de grande importância na Medicina Veterinária, pois tratam-se de enfermidades de elevada incidência e com poucos subsídios auxiliares no seu diagnóstico, prognóstico e na avaliação de terapias empregadas. Ainda hoje, o diagnóstico baseia-se, em grande parte, no histórico e no exame clínico neurológico. Dessa forma, a análise dos constituintes do fluido cefalorraquidiano torna-se uma das poucas alternativas de acesso clínico ao sistema nervoso central (SNC). Mesmo com a grande utilidade do exame físico-químico e citoscópico do liquor na neurologia veterinária, poucos são os estudos sobre a estabilidade dos seus constituintes sob estocagem. Dessa forma, o presente trabalho teve como finalidade verificar a influência da temperatura e do tempo de conservação nas características físico-químicas do liquor de cães hígidos. Para tanto, foram coletadas amostras de LCR, através da punção da cisterna cerebelo-medular de cães clinicamente sadios, as quais foram submetidas à análise da densidade específica, do pH, da glicorraquia, das proteínas totais e das atividades das enzimas creatina quinase (CK) e aspartato aminotransferase (AST), após conservá-las em diferentes temperaturas (25°C, 4°C e -4°C) e por diferentes períodos de tempo (logo após a colheita, 24 horas, 48 horas, uma semana e um mês). Dentre os resultados obtidos, foi possível verificar, principalmente, que houve estabilidade dos parâmetros estudados por até um mês de estocagem nas amostras mantidas sob a temperaturas de congelamento de -4°C.

Análise da deformação pós-Rifte na Bacia do Araripe, Nordeste do Brasil

The Araripe Basin is located over Precambrian terrains of the Borborema Province, being part of Northeast Brazil inner basins. Its origin is related to the fragmentation of the Gondwana supercontinent and consequently opening of South Atlantic during early Cretaceous. The basin has a sedimentary infill encompassing four distinct evolution stages, comprising Paleozoic syneclisis, pre-rift, rift and post-rift. The target of this study comprises the post-rift section of the basin focusing deformational styles which affect evaporates from Ipubi Member of the Santana Formation, which is composed by gypsum and anidrite layers interbedded with shales. These units occur widespread across the basin. In the central part of the basin, near Nova Olinda-Santana do Cariri, evaporites are affected by an essentialy brittle deformation tipified by fibrous gypsum filled fractures, cutting massive layers of gypsum and anidrite. Veins with variable orientations and dips are observed in the region distributed over three main populations: i) a dominant NWSE with shallow to moderate NE dipping population, consisting of gypsum filled veins in which fibers are normal to vein walls; i) NE-SW veins with moderate SE dips containing subhorizontal growth fibers; and iii) N-S veins with shallow E-W dips with fibers oblique to vein walls. In the west portion of the basin, near Trindade-Ipubi-Araripina towns, evaporate layers are dominantly constituted by gypsum/anidrite finely stratified, showing a minor density of veins. These layers are affected by a unique style of deformation, more ductile, typified by gentle to open horizontal normal folding with several tens of meters length and with double plunging NW-SE or NE-SW hinges, configuring domic features. In detail, gypsum/anidrite laminae are affected by metre to decimeter scale close to tight folding, usually kinked, with broken hinges, locally turning into box folds. Veins show NE-SW main directions with shallow NE dips, growth fibers are parallel to vein walls, constituting slickenfibers. This region is marked by faults that affect Araripina Formation with NW-SE, NE-SW and E-W directions. The main structural styles and general orientations of structures which affected the post-rift section of Araripe Basin yielded important kinematic information analysis which led us to infer a E-W to NE-SW extension direction to the northeastern part of the Basin, whereas in the southeastern part, extension occurred in N-S direction. Thus, it was possible to determine a regional kinematic setting, through this analysis, characterizing a NE-SW to ENE-WSW system for the post-rift section, which is compatible with the tension settings for the Sout American Plate since Albian. Local variations at the fluid pressure linked (or not) to sedimentary overload variation define local tension settings. This way, at the northeastern portion of the basin, the post-rift deformation was governed by a setting which σ 1 is sub-horizontal trending NE-SW and, σ 3 is sub-vertical, emphasizing a reverse fault situation. At the southwestern portion however there was characterized a strike slip fault setting, featuring σ 1 trending ENEWSW and σ3 trending NNW-SSE

Glial reactivity in dogs with visceral leishmaniasis: correlation with T lymphocyte infiltration and with cerebrospinal fluid anti-Leishmania antibody titres

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

LEAST ACTION PRINCIPLE and THE INCOMPRESSIBLE EULER EQUATIONS WITH VARIABLE DENSITY

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Limits of validity for a semiclassical mean-field two-fluid model for Bose-Einstein condensation thermodynamics

We reinvestigate the Bose-Einstein condensation (BEC) thermodynamics of a weakly interacting dilute Bose gas under the action of a trap using a semi-classical two-fluid mean-field model in order to find the domain of applicability of the model. Such a model is expected to break down once the condition of diluteness and weak interaction is violated. We find that this breakdown happens for values of coupling and density near the present experimental scenario of BEG. With the increase of the interaction coupling and density the model may lead to unphysical results for thermodynamic observables. (C) 2000 Published by Elsevier B.V. B.V, All rights reserved.

Universal scaling in a trapped Fermi super-fluid in the BCS-unitarity crossover

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

The mineralizing fluid in the Piaba gold deposit, São Luís cratonic fragment (NW-Maranhão, Brazil) based on fluid inclusion studies in quartz veins

O depósito aurífero de Piaba tornou-se a primeira mina em operação no fragmento cratônico São Luís, noroeste do Maranhão. Seu ambiente geológico compreende rochas metavulcanossedimentares do Grupo Aurizona e granitoides da Suíte Tromaí, entre outras unidades menores, formadas em ambiente de arcos de ilhas entre 2240 e 2150 Ma, juntamente com outras unidades menores. A mineralização é caracterizada por uma trama stockwork de veios e vênulas de quartzo com seus halos de alteração (clorita + muscovita + carbonato + pirita + calcopirita e ouro) hospedada em um granodiorito granofírico fino (Granófiro Piaba) e em rocha subvulcânica andesítica do Grupo Aurizona. O corpo mineralizado é espacialmente limitado à zona de cisalhamento ENE-WSW rúptil-dúctil (Falha Piaba). Estudos petrográficos, microtermométricos e por espectroscopia microRaman no quartzo definiram inclusões aquo-carbônicas bifásicas e trifásicas, produzidas por aprisionamento heterogêneo durante separação de fases, e fluidos aquosos tardios. A solução mineralizadora corresponde a um fluido aquo-carbônico composto por CO₂ (5 - 24 mol%, densidade de 0,96 - 0,99 g/cm³), H₂O (74 - 93 mol%), N₂ (< 1 mol%), CH₄ (<1mol%) e 5,5 % em peso NaCl equivalente. O minério depositou a 267 - 302ºC e 1,25 - 2,08 kbar, correspondendo a profundidades de 4 a 7 km, em consonância com o regime estrutural. A composição e o intervalo de P-T do fluido mineralizador, combinados com o caráter redutor (log ƒ_O2 -31,3 a -34,3) e a sulfetação da rocha hospedeira, sugerem que o ouro foi transportado como um complexo sulfetado. O minério foi depositado em consequência da separação de fase, redução da atividade de enxofre e da ƒO₂ pela interação fluido-rocha.

Sampling density and proportion for the characterization of the variability of Oxisol attributes on different materials

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

A snapshot of suspended sediment and fluid mud occurrence in a mixed-energy embayment, Tijucas Bay, Brazil

Along the southern Brazilian coast, Tijucas Bay is known for its unique muddy tidal flats associated with chenier plains. Previous field observations pointed to very high suspended sediment concentrations (SSCs) in the inner parts of the bay, and in the estuary of the Tijucas River, suggesting the presence of fluid mud. In this study, the occurrences of suspended sediments and fluid mud were examined during a larger-scale, high-resolution 2-day field campaign on 1-2 May 2007, encompassing survey lines spanning nearly 80 km, 75 water sampling stations for near-bottom density estimates, and ten sediment sampling stations. Wave refraction modeling provided qualitative wave energy estimates as a function of different incidence directions. The results show that SSC increases toward the inner bay near the water surface, but seaward near the bottom. This suggests that suspended sediment is supplied by the local rivers, in particular the Tijucas. Near-surface SSCs were of the order of 50 mg l(-1) close to the shore, but exceeded 100 mg l(-1) near the bottom in the deeper parts of the bay. Fluid mud thickness and location given by densimetry and echo-sounding agreed in some places, although being mostly discordant. The best agreement was observed where wave energy was high during the campaign. The discrepancy between the two methods may be an indication for the existence of fluid mud, which is recorded by one method but not the other. Agreement is considered to be an indication of fluidization, whereas disagreement indicates more consolidation. Wave modeling suggests that waves from the ENE and SE are the most effective in supplying energy to the inner bay, which may induce the liquefaction of mud deposits to form fluid mud. Nearshore mud resuspension and weak horizontal currents result in sediment-laden offshore flow, which explains the higher SSCs measured in the deeper parts of the bay, besides providing a mechanism for fine-sediment export to the inner shelf.