Efficient solutions to robust, semi-implicit discretizations of the immersed boundary method

The immersed boundary method is a versatile tool for the investigation of flow-structure interaction. In a large number of applications, the immersed boundaries or structures are very stiff and strong tangential forces on these interfaces induce a well-known, severe time-step restriction for explicit discretizations. This excessive stability constraint can be removed with fully implicit or suitable semi-implicit schemes but at a seemingly prohibitive computational cost. While economical alternatives have been proposed recently for some special cases, there is a practical need for a computationally efficient approach that can be applied more broadly. In this context, we revisit a robust semi-implicit discretization introduced by Peskin in the late 1970s which has received renewed attention recently. This discretization, in which the spreading and interpolation operators are lagged. leads to a linear system of equations for the inter-face configuration at the future time, when the interfacial force is linear. However, this linear system is large and dense and thus it is challenging to streamline its solution. Moreover, while the same linear system or one of similar structure could potentially be used in Newton-type iterations, nonlinear and highly stiff immersed structures pose additional challenges to iterative methods. In this work, we address these problems and propose cost-effective computational strategies for solving Peskin`s lagged-operators type of discretization. We do this by first constructing a sufficiently accurate approximation to the system`s matrix and we obtain a rigorous estimate for this approximation. This matrix is expeditiously computed by using a combination of pre-calculated values and interpolation. The availability of a matrix allows for more efficient matrix-vector products and facilitates the design of effective iterative schemes. We propose efficient iterative approaches to deal with both linear and nonlinear interfacial forces and simple or complex immersed structures with tethered or untethered points. One of these iterative approaches employs a splitting in which we first solve a linear problem for the interfacial force and then we use a nonlinear iteration to find the interface configuration corresponding to this force. We demonstrate that the proposed approach is several orders of magnitude more efficient than the standard explicit method. In addition to considering the standard elliptical drop test case, we show both the robustness and efficacy of the proposed methodology with a 2D model of a heart valve. (C) 2009 Elsevier Inc. All rights reserved.

On state representations of nonlinear implicit systems

This work considers a semi-implicit system A, that is, a pair (S, y), where S is an explicit system described by a state representation (x)over dot(t) = f(t, x(t), u(t)), where x(t) is an element of R(n) and u(t) is an element of R(m), which is subject to a set of algebraic constraints y(t) = h(t, x(t), u(t)) = 0, where y(t) is an element of R(l). An input candidate is a set of functions v = (v(1),.... v(s)), which may depend on time t, on x, and on u and its derivatives up to a Finite order. The problem of finding a (local) proper state representation (z)over dot = g(t, z, v) with input v for the implicit system Delta is studied in this article. The main result shows necessary and sufficient conditions for the solution of this problem, under mild assumptions on the class of admissible state representations of Delta. These solvability conditions rely on an integrability test that is computed from the explicit system S. The approach of this article is the infinite-dimensional differential geometric setting of Fliess, Levine, Martin, and Rouchon (1999) (`A Lie-Backlund Approach to Equivalence and Flatness of Nonlinear Systems`, IEEE Transactions on Automatic Control, 44(5), (922-937)).

A Robust, Fully Adaptive Hybrid Level-Set/Front-Tracking Method for Two-Phase Flows with an Accurate Surface Tension Computation

We present a variable time step, fully adaptive in space, hybrid method for the accurate simulation of incompressible two-phase flows in the presence of surface tension in two dimensions. The method is based on the hybrid level set/front-tracking approach proposed in [H. D. Ceniceros and A. M. Roma, J. Comput. Phys., 205, 391400, 2005]. Geometric, interfacial quantities are computed from front-tracking via the immersed-boundary setting while the signed distance (level set) function, which is evaluated fast and to machine precision, is used as a fluid indicator. The surface tension force is obtained by employing the mixed Eulerian/Lagrangian representation introduced in [S. Shin, S. I. Abdel-Khalik, V. Daru and D. Juric, J. Comput. Phys., 203, 493-516, 2005] whose success for greatly reducing parasitic currents has been demonstrated. The use of our accurate fluid indicator together with effective Lagrangian marker control enhance this parasitic current reduction by several orders of magnitude. To resolve accurately and efficiently sharp gradients and salient flow features we employ dynamic, adaptive mesh refinements. This spatial adaption is used in concert with a dynamic control of the distribution of the Lagrangian nodes along the fluid interface and a variable time step, linearly implicit time integration scheme. We present numerical examples designed to test the capabilities and performance of the proposed approach as well as three applications: the long-time evolution of a fluid interface undergoing Rayleigh-Taylor instability, an example of bubble ascending dynamics, and a drop impacting on a free interface whose dynamics we compare with both existing numerical and experimental data.

A boundary element formulation for analysis of elastoplastic plates with geometrical nonlinearity

In this paper a new boundary element method formulation for elastoplastic analysis of plates with geometrical nonlinearities is presented. The von Mises criterion with linear isotropic hardening is considered to evaluate the plastic zone. Large deflections are assumed but within the context of small strain. To derive the boundary integral equations the von Karman`s hypothesis is taken into account. An initial stress field is applied to correct the true stresses according to the adopted criterion. Isoparametric linear elements are used to approximate the boundary unknown values while triangular internal cells with linear shape function are adopted to evaluate the domain value influences. The nonlinear system of equations is solved by using an implicit scheme together with the consistent tangent operator derived along the paper. Numerical examples are presented to demonstrate the accuracy and the validity of the proposed formulation.

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

Análise de um ciclone semi-estacionário na costa sul do Brasil associado a bloqueio atmosférico

Na primeira semana de maio de 2008, durante quatro dias, um ciclone em superfície permaneceu semi-estacionário na costa da região sul do Brasil. Este sistema foi responsável por chuvas e ventos fortes no Rio Grande do Sul e Santa Catarina, os quais causaram muitos danos (queda de árvores, enchentes e desabamentos). O objetivo deste trabalho é avaliar o processo de formação e entender os mecanismos responsáveis pelo lento deslocamento do ciclone, já que a maioria dos ciclones nesta região possui deslocamento mais rápido. A equação de desenvolvimento de Sutcliffe mostrou que a advecção de vorticidade absoluta ciclônica na média troposfera e a advecção de ar quente na camada entre 1000-500 hPa foram mecanismos importantes para a ciclogênese. Neste período, o intenso aquecimento diabático também contribuiu para a ciclogênese, à medida que se contrapôs ao intenso resfriamento adiabático devido aos movimentos verticais ascendentes. A advecção de vorticidade absoluta ciclônica que favoreceu a ciclogênese esteve associada a um Vórtice Ciclônico em Altos Níveis (VCAN), que se formou numa região de anomalia de vorticidade potencial. O VCAN se manteve semi-estacionário e compôs o setor norte de um bloqueio do tipo dipolo. Tal bloqueio intensificou um anticiclone em superfície, situado a sul/leste do ciclone, o que contribuiu para o ciclone se manter semi-estacionário. O movimento atípico e lento do ciclone para sul, e em alguns períodos para sudoeste, esteve associado com advecções de vorticidade absoluta ciclônica na média troposfera e de ar quente no seu setor sul. Somente quando o bloqueio em níveis médios e a anomalia de vorticidade potencial em níveis médios/altos se enfraqueceram, o ciclone em superfície se afastou da costa sul do Brasil.

Convecção linearmente organizada na área de Petrolina, semi-árido do Nordeste do Brasil: aspectos em meso e grande escala

É apresentado um estudo sobre sistemas convectivos linearmente organizados e observados por um radar meteorológico banda-C na região semi-árida do Nordeste do Brasil. São analisados três dias (27 a 29) de março de 1985, com ênfase na investigação do papel desempenhado por fatores locais e de grande escala no desenvolvimento dos sistemas. No cenário de grande escala, a área de cobertura do radar foi influenciada por um cavado de ar superior austral no dia 27 e por um vórtice ciclônico de altos níveis no dia 29. A convergência de umidade próxima à superfície favoreceu a atividade convectiva nos dias 27 e 29, enquanto que divergência de umidade próxima à superfície inibiu a atividade convectiva no dia 28. No cenário de mesoescala, foi observado que o aquecimento diurno é um fator importante para a formação de células convectivas, somando-se a ele o papel determinante da orografia na localização dos ecos. De maneira geral, as imagens de radar mostram os sistemas convectivos linearmente organizados em áreas elevadas e núcleos convectivos intensos envolvidos por uma área de precipitação estratiforme. Os resultados indicam que convergência do fluxo de umidade em grande escala e aquecimento radiativo, são fatores determinantes na evolução e desenvolvimento dos ecos na área de estudo.

Estudo da mobilização de metais e elementos traços em ambientes aquáticos do semi-árido brasileiro aplicando análises de componentes principais

The spatial and temporal retention of metals has been studied in water and sediments of the Gavião River, Anagé and Tremedal Reservoirs, located in the semi-arid region, Bahia - Brazil, in order to identify trends in the fluxes of metals from the sediments to the water column. The determination of metals was made by ICP OES and ET AAS. The application of statistical methods showed that this aquatic system presents suitable conditions to move Cd2+ and Pb2+ from the water column to the sediment.

Rapid detection of bovine coronavirus by a semi-nested RT-PCR

Bovine coronavirus (BCoV) is a member of the group 2 of the Coronavirus (Nidovirales: Coronaviridae) and the causative agent of enteritis in both calves and adult bovine, as well as respiratory disease in calves. The present study aimed to develop a semi-nested RT-PCR for the detection of BCoV based on representative up-to-date sequences of the nucleocapsid gene, a conserved region of coronavirus genome. Three primers were designed, the first round with a 463bp and the second (semi-nested) with a 306bp predicted fragment. The analytical sensitivity was determined by 10-fold serial dilutions of the BCoV Kakegawa strain (HA titre: 256) in DEPC treated ultra-pure water, in fetal bovine serum (FBS) and in a BCoV-free fecal suspension, when positive results were found up to the 10-2, 10-3 and 10-7 dilutions, respectively, which suggests that the total amount of RNA in the sample influence the precipitation of pellets by the method of extraction used. When fecal samples was used, a large quantity of total RNA serves as carrier of BCoV RNA, demonstrating a high analytical sensitivity and lack of possible substances inhibiting the PCR. The final semi-nested RT-PCR protocol was applied to 25 fecal samples from adult cows, previously tested by a nested RT-PCR RdRp used as a reference test, resulting in 20 and 17 positives for the first and second tests, respectively, and a substantial agreement was found by kappa statistics (0.694). The high sensitivity and specificity of the new proposed method and the fact that primers were designed based on current BCoV sequences give basis to a more accurate diagnosis of BCoV-caused diseases, as well as to further insights on protocols for the detection of other Coronavirus representatives of both Animal and Public Health importance.

Lysine and zinc chelate in diets for brown laying hens: effects on egg production and composition

It was evaluated the effects of levels of digestible lysine and chelate zinc combined in the diet for laying on the egg quality. It was used 720 birds, from 48 to 60 weeks of age, distributed in a completely randomized design in a 3 × 5 factorial scheme with three levels of zinc and five levels of lysine, applied into six replicates in the experimental units of eight birds per plot. The levels were: 137, 309 and 655 ppm zinc and 0.482, 0.527, 0.582, 0.644 and 0.732% digestible lysine. It was not observed any interaction among digestible lysine and zinc for the primary variables of fractions and egg composition. Levels of zinc reduced egg weight, suggesting the lowest efficiency in nutrient intake. At the highest dietary concentration of zinc, the addition of digestible lysine coincided with a linear increase in shell weight. However, zinc addition, regardless of lysine level in the diet, resulted in the reduction of egg weight and of the percentage of mineral matter in the yolk, limiting the efficiency of mineral deposition in this fraction of the egg. Concentration of zinc that produced the best results was 137 ppm inasmuch as higher quantities limit the use of digestible lysine, with effects harming composition and egg quality. The study indicates the following requirement for digestible lysine: 0.639% from the 48th to the 52nd week, 0.679% from the 52nd to 56th week, and 0.635% from the 56th to the 60th week. Considering the total period from 48th to the 60th week, the level 0.638% of lysine or the daily intake of 707 mg of the amino acid met the requirement for egg quality of semi-heavy layers used in this study.

Antileishmanial Activity of the Hydroalcoholic Extract of Miconia langsdorffii, Isolated Compounds, and Semi-Synthetic Derivatives

The in vitro activity of the crude hydroalcoholic extract of the aerial parts of Miconia langsdorffii Cogn. was evaluated against the promastigote forms of L. amazonensis, the causative agent of cutaneous leishmaniasis in humans. The bioassay-guided fractionation of this extract led to identification of the triterpenes ursolic acid and oleanolic acid as the major compounds in the fraction that displayed the highest activity. Several ursolic acid semi-synthetic derivatives were prepared, to find out whether more active compounds could be obtained. Among these ursolic acid-derived substances, the C-28 methyl ester derivative exhibited the best antileishmanial activity.

The merging cluster Abell 1758 revisited: multi-wavelength observations and numerical simulations

Context. Cluster properties can be more distinctly studied in pairs of clusters, where we expect the effects of interactions to be strong. Aims. We here discuss the properties of the double cluster Abell 1758 at a redshift z similar to 0.279. These clusters show strong evidence for merging. Methods. We analyse the optical properties of the North and South cluster of Abell 1758 based on deep imaging obtained with the Canada-France-Hawaii Telescope (CFHT) archive Megaprime/Megacam camera in the g' and r' bands, covering a total region of about 1.05 x 1.16 deg(2), or 16.1 x 17.6 Mpc(2). Our X-ray analysis is based on archive XMM-Newton images. Numerical simulations were performed using an N-body algorithm to treat the dark-matter component, a semi-analytical galaxy-formation model for the evolution of the galaxies and a grid-based hydrodynamic code with a parts per million (PPM) scheme for the dynamics of the intra-cluster medium. We computed galaxy luminosity functions (GLFs) and 2D temperature and metallicity maps of the X-ray gas, which we then compared to the results of our numerical simulations. Results. The GLFs of Abell 1758 North are well fit by Schechter functions in the g' and r' bands, but with a small excess of bright galaxies, particularly in the r' band; their faint-end slopes are similar in both bands. In contrast, the GLFs of Abell 1758 South are not well fit by Schechter functions: excesses of bright galaxies are seen in both bands; the faint-end of the GLF is not very well defined in g'. The GLF computed from our numerical simulations assuming a halo mass-luminosity relation agrees with those derived from the observations. From the X-ray analysis, the most striking features are structures in the metal distribution. We found two elongated regions of high metallicity in Abell 1758 North with two peaks towards the centre. In contrast, Abell 1758 South shows a deficit of metals in its central regions. Comparing observational results to those derived from numerical simulations, we could mimic the most prominent features present in the metallicity map and propose an explanation for the dynamical history of the cluster. We found in particular that in the metal-rich elongated regions of the North cluster, winds had been more efficient than ram-pressure stripping in transporting metal-enriched gas to the outskirts. Conclusions. We confirm the merging structure of the North and South clusters, both at optical and X-ray wavelengths.

Quantitative constraints on the transport properties of hot partonic matter from semi-inclusive single high transverse momentum pion suppression in Au plus Au collisions at root S(NN)=200 GeV

The PHENIX experiment has measured the suppression of semi-inclusive single high-transverse-momentum pi(0)'s in Au+Au collisions at root s(NN) = 200 GeV. The present understanding of this suppression is in terms of energy loss of the parent (fragmenting) parton in a dense color-charge medium. We have performed a quantitative comparison between various parton energy-loss models and our experimental data. The statistical point-to-point uncorrelated as well as correlated systematic uncertainties are taken into account in the comparison. We detail this methodology and the resulting constraint on the model parameters, such as the initial color-charge density dN(g)/dy, the medium transport coefficient <(q) over cap >, or the initial energy-loss parameter epsilon(0). We find that high-transverse-momentum pi(0) suppression in Au+Au collisions has sufficient precision to constrain these model-dependent parameters at the +/- 20-25% (one standard deviation) level. These constraints include only the experimental uncertainties, and further studies are needed to compute the corresponding theoretical uncertainties.

Amorphous HfO(2) and Hf(1-x)Si(x)O via a melt-and-quench scheme using ab initio molecular dynamics

We have performed ab initio molecular dynamics simulations to generate an atomic structure model of amorphous hafnium oxide (a-HfO(2)) via a melt-and-quench scheme. This structure is analyzed via bond-angle and partial pair distribution functions. These results give a Hf-O average nearest-neighbor distance of 2.2 angstrom, which should be compared to the bulk value, which ranges from 1.96 to 2.54 angstrom. We have also investigated the neutral O vacancy and a substitutional Si impurity for various sites, as well as the amorphous phase of Hf(1-x)Si(x)O(2) for x=0.25, 0375, and 0.5.

A Semi-Tethered Test for Power Assessment in Running

We analyzed the usefulness of a semi-tethered field running test (STR) and the relationships between indices of anaerobic power, anaerobic capacity and running performance in 9 trained male sprinters (22.2 +/- 2.9 yrs, 176 +/- 1 cm, 68.0 +/- 9.4 kg). STR involved an all out 120 m run attached to an apparatus that enabled power calculation from force and velocity measures. Subjects also carried out a cycloergometer Win-gate Anaerobic Test (WT), an all out 300 m run and had accessed their maximal accumulated oxygen deficit (MAOD) on a treadmill. Peak and mean powers attained in STR (1 720 +/- 221 and 1 391 +/- 201 W) were greater but significantly related (r=0.82; P<0.01) to those in the WT (808 +/- 130 and 603 +/- 87 W). In addition, power measures derived from the STR were stronger related to running performance compared to those from the WT (r=0.81-0.94 vs. 0.68-0.84; P<0.05). Relationships between MAOD and most power indices were only weak to moderate. These results support the usefulness of STR for specific power assessment in field running and suggest that anaerobic power and capacity are not related entities, irrespective of having been evaluated using similar or dissimilar exercise modes.