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.

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

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.

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)).

New methods to quantify NH3 volatilization from fertilized surface soil with urea

Gaseous N losses from soil are considerable, resulting mostly from ammonia volatilization linked to agricultural activities such as pasture fertilization. The use of simple and accessible measurement methods of such losses is fundamental in the evaluation of the N cycle in agricultural systems. The purpose of this study was to evaluate quantification methods of NH3 volatilization from fertilized surface soil with urea, with minimal influence on the volatilization processes. The greenhouse experiment was arranged in a completely randomized design with 13 treatments and five replications, with the following treatments: (1) Polyurethane foam (density 20 kg m-3) with phosphoric acid solution absorber (foam absorber), installed 1, 5, 10 and 20 cm above the soil surface; (2) Paper filter with sulfuric acid solution absorber (paper absorber, 1, 5, 10 and 20 cm above the soil surface); (3) Sulfuric acid solution absorber (1, 5 and 10 cm above the soil surface); (4) Semi-open static collector; (5) 15N balance (control). The foam absorber placed 1 cm above the soil surface estimated the real daily rate of loss and accumulated loss of NH3N and proved efficient in capturing NH3 volatized from urea-treated soil. The estimates based on acid absorbers 1, 5 and 10 cm above the soil surface and paper absorbers 1 and 5 cm above the soil surface were only realistic for accumulated N-NH3 losses. Foam absorbers can be indicated to quantify accumulated and daily rates of NH3 volatilization losses similarly to an open static chamber, making calibration equations or correction factors unnecessary.

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.

Contraceptive methods with male participation: a perspective of brazilian couples

OBJECTIVE: To assess the perspectives of couples who requested vasectomy in a public health service on the use of male participation contraceptive methods available in Brazil: male condoms, natural family planning/calendar, coitus interruptus and vasectomy. METHODS: A qualitative study with semi-structured interviews was held with 20 couples who had requested vasectomy at the Human Reproduction Unit of the Universidade Estadual de Campinas, Brazil. Data analysis was carried out through thematic content analysis. FINDINGS: The couples did not, in general, know any effective contraceptive options for use by men and/or participating in their use, except for vasectomy. The few methods with male participation that they knew of were perceived to interfere in spontaneity and in pleasure of intercourse. Men accepted that condom use in extra-conjugal relations offered them protection from sexually transmitted diseases; that their wives might also participate in extra-marital relationships was not considered. DISCUSSION: The few contraceptive options with male participation lead to difficulty in sharing responsibilities between men and women. On the basis of perceived gender roles, women took the responsibility for contraception until the moment when the situation became untenable, and they faced the unavoidable necessity of sterilization. CONCLUSIONS: Specific actions are necessary for men to achieve integral participation in relation to reproductive sexual health. These include education and discussions on gender roles, leading to greater awareness in men of the realities of sexual and reproductive health

A Study on Unconditionally Stable FDTD Methods for the Modeling of Metamaterials

We assess the performance of three unconditionally stable finite-difference time-domain (FDTD) methods for the modeling of doubly dispersive metamaterials: 1) locally one-dimensional FDTD; 2) locally one-dimensional FDTD with Strang splitting; and (3) alternating direction implicit FDTD. We use both double-negative media and zero-index media as benchmarks.

Analysis of a semi-rigid connection for precast concrete

This paper presents a study of a specific type of beam-to-column connection for precast concrete structures. Furthermore, an analytical model to determine the strength and the stiffness of the connection, based on test results of two prototypes, is proposed. To evaluate the influence of the strength and stiffness of the connection on the behaviour of the structure, the results of numerical simulations of a typical multi-storey building with semi-rigid connections are also presented and compared with the results using pinned and rigid connections. The main conclusions are: (a) the proposed design model can reasonably evaluate the studied connection strength; (b) the evaluation of strength is more accurate than that of stiffness; (c) for a typical structure, it is possible to increase the number of storeys of the structure from two to four with lower horizontal displacement at the top, and only a small increase of the column base bending moment by replacing the pinned connections with semi-rigid ones; and (d) although there is significant uncertainty in the connection stiffness, the results show that the displacements at the top of the structure, and the column base moments present low susceptibility deviations to this parameter.

Repeated disturbances and canopy disturbance regime in a tropical semi-deciduous forest

The canopy disturbance regime and the influence of gap methods on the interpretation of forest structure and dynamics were evaluated in a tropical semi-deciduous forest in south-eastern Brazil. We encountered a gap density of 11.2 gaps ha(-1) and an average size which varied from 121 to 333 m(2) depending on the gap delimitation method considered (minimum gap size was 10 m(2)). Although average size was slightly higher, the median value obtained (78 m(2)) was comparable to other tropical forest sites and the gap size-class distribution found supported the pattern described for such forest sites. Among 297 gap makers, snapping and uprooting were the most common modes of disturbance. The number and basal area of gap makers were good predictors of gap size. Almost 25% of all gaps suffered from repeated disturbance events that brought about larger gap sizes. Such processes, along with delimitation methods, strongly influenced the estimation of turnover rate and therefore the interpretation of forest dynamics. These results demonstrated the importance of further studies on repeated disturbances, which is often neglected in forest studies.

Tree-ring formation, radial increment periodicity, and phenology of tree species from a seasonal semi-deciduous forest in southeast Brazil

Many tropical tree species produce growth rings in response to seasonal environmental factors that influence the activity of the vascular cambium. We applied the following methods to analyze the annual nature of treering formation of 24 tree species from a seasonal semi-deciduous forest of southeast Brazil: describing wood anatomy and phenology, counting tree rings after cambium markings, and using permanent dendrometer bands. After 7 years of systematic observations and measurements, we found the following: the trees lost their leaves during the dry season and grew new leaves at the end of the same season; trunk increment dynamics corresponded to seasonal changes in precipitation, with higher increment (active period) during the rainy season (October-April) and lower increment (dormant period) during the dry season (May-September); the number of tree rings formed after injuries to the cambium coincided with the number of years since the extraction of the wood samples. As a result of these observations, it was concluded that most study trees formed one growth ring per year. This suggests that tree species from the seasonal semi-deciduous forests of Brazil have an annual cycle of wood formation. Therefore, these trees have potential for use in future studies of tree age and radial growth rates, as well as to infer ecological and regional climatic conditions. These future studies can provide important information for the management and conservation of these endangered forests.

Distinct conformational properties determined by implicit and explicit representation of protein-solvent interactions. An analytical and computer simulation study

In the protein folding problem, solvent-mediated forces are commonly represented by intra-chain pairwise contact energy. Although this approximation has proven to be useful in several circumstances, it is limited in some other aspects of the problem. Here we show that it is possible to achieve two models to represent the chain-solvent system. one of them with implicit and other with explicit solvent, such that both reproduce the same thermodynamic results. Firstly, lattice models treated by analytical methods, were used to show that the implicit and explicitly representation of solvent effects can be energetically equivalent only if local solvent properties are time and spatially invariant. Following, applying the same reasoning Used for the lattice models, two inter-consistent Monte Carlo off-lattice models for implicit and explicit solvent are constructed, being that now in the latter the solvent properties are allowed to fluctuate. Then, it is shown that the chain configurational evolution as well as the globule equilibrium conformation are significantly distinct for implicit and explicit solvent systems. Actually, strongly contrasting with the implicit solvent version, the explicit solvent model predicts: (i) a malleable globule, in agreement with the estimated large protein-volume fluctuations; (ii) thermal conformational stability, resembling the conformational hear resistance of globular proteins, in which radii of gyration are practically insensitive to thermal effects over a relatively wide range of temperatures; and (iii) smaller radii of gyration at higher temperatures, indicating that the chain conformational entropy in the unfolded state is significantly smaller than that estimated from random coil configurations. Finally, we comment on the meaning of these results with respect to the understanding of the folding process. (C) 2009 Elsevier B.V. All rights reserved.

Evaluation of flow regimes in a semi-cylindrical spouted bed through statistical, mutual information, spectral and Hurst`s analysis

Hydrodynamic studies were conducted in a semi-cylindrical spouted bed column of diameter 150 mm, height 1000 mm, conical base included angle of 60 degrees and inlet orifice diameter 25 mm. Pressure transducers at several axial positions were used to obtain pressure fluctuation time series with 1.2 and 2.4 mm glass beads at U/U-ms from 0.3 to 1.6, and static bed depths from 150 to 600 mm. The conditions covered several flow regimes (fixed bed, incipient spouting, stable spouting, pulsating spouting, slugging, bubble spouting and fluidization). Images of the system dynamics were also acquired through the transparent walls with a digital camera. The data were analyzed via statistical, mutual information theory, spectral and Hurst`s Rescaled Range methods to assess the potential of these methods to characterize the spouting quality. The results indicate that these methods have potential for monitoring spouted bed operation.

Approximate controllability of second-order distributed implicit functional systems

In this paper we study the approximate controllability of control systems with states and controls in Hilbert spaces, and described by a second-order semilinear abstract functional differential equation with infinite delay. Initially we establish a characterization for the approximate controllability of a second-order abstract linear system and, in the last section, we compare the approximate controllability of a semilinear abstract functional system with the approximate controllability of the associated linear system. (C) 2008 Elsevier Ltd. All rights reserved.

The influence of oscillations on global existence for a class of semi-linear wave equations

The goal of this paper is to study the global existence of small data solutions to the Cauchy problem for the nonlinear wave equation u(tt) - a(t)(2) Delta u = u(t)(2) - a(t)(2)vertical bar del u vertical bar(2). In particular we are interested in statements for the 1D case. We will explain how the interplay between the increasing and oscillating behavior of the coefficient will influence global existence of small data solutions. Copyright c 2011 John Wiley & Sons, Ltd.