An implicit technique for solving 3D low Reynolds number moving free surface flows

This paper describes the development of an implicit finite difference method for solving transient three-dimensional incompressible free surface flows. To reduce the CPU time of explicit low-Reynolds number calculations, we have combined a projection method with an implicit technique for treating the pressure on the free surface. The projection method is employed to uncouple the velocity and the pressure fields, allowing each variable to be solved separately. We employ the normal stress condition on the free surface to derive an implicit technique for calculating the pressure at the free surface. Numerical results demonstrate that this modification is essential for the construction of methods that are more stable than those provided by discretizing the free surface explicitly. In addition, we show that the proposed method can be applied to viscoelastic fluids. Numerical results include the simulation of jet buckling and extrudate swell for Reynolds numbers in the range [0.01, 0.5]. (C) 2008 Elsevier Inc. All rights reserved.

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.

Implicit corporate political activity and elite formation

Nesta tese, analisamos como a elite empresarial progressista criou a organização da sociedade civil Rede Nossa São Paulo (RNSP), alcançando mudanças institucionais significativas, permitindo assim a consolidação da elite na esfera política. A pesquisa resultou em três artigos. O primeiro artigo discute como a RNSP se tornou um forte ator político na cidade de São Paulo e também no Brasil. Para abordar esta questão, mostramos como a RNSP usou a história retórica para se tornar um ator central na esfera política. No segundo artigo, propomos o conceito de atividade política corporativa implícita (ICPA), complementar a atividade política corporativa. Conceituamos ICPA como elites empresariais em conjunto com organizações da sociedade civil agindo para influenciar o governo. Com os limites entre o governo, as empresas e organizações da sociedade civil difusos; entendemos que este conceito é extremamente importante para chamar a atenção e criar novos caminhos para a pesquisa sobre a influência das empresas no governo. No último artigo, mostramos os micro fundamentos da ICPA. Especificamente, como as elites empresariais e corporações influenciam a RNSP e, indiretamente, o governo. Concluindo, contribuímos para a literatura sobre a influência das empresas no governo e na esfera pública indiretamente, por meio de organizações da sociedade civil. Teoricamente, estendemos a literatura de teoria institucional, história e poder

Optimal provision of implicit and explicit incentives in asset management contracts

This paper investigates the importance of the fiow of funds as an implicit incetive provided by investors to portfolio managers in a two-period relationship. We show that the fiow of funds is a powerful incentive in an asset management contract. We build a binomial moral hazard model to explain the main trade-ofIs in the relationship between fiow, fees and performance. The main assumption is that efIort depend" on the combination of implicit and explicit incentives while the probability distrioutioll function of returns depends on efIort. In the case of full commitment, the investor's relevant trade-ofI is to give up expected return in the second period vis-à-vis to induce efIort in the first período The more concerned the investor is with today's payoff. the more willing he will be to give up expected return in the following periods. That is. in the second period, the investor penalizes observed low returns by withdrawing resources from non-performing portfolio managers. Besides, he pays performance fee when the observed excess return is positive. When commitment is not a plausible hypothesis, we consider that the investor also learns some symmetríc and imperfect information about the ability of the manager to generate positive excess returno In this case, observed returns reveal ability as well as efIort choices exerted by the portfolio manager. We show that implicit incentives can explain the fiow-performance relationship and, conversely, endogenous expected return determines incentives provision and define their optimal leveIs. We provide a numerical solution in Matlab that characterize these results.

A combination of implicit and adaptative upwind tools for the numerical solution of incompressible free surface flows

This paper is concerned with the numerical solutions of time dependent two-dimensional incompressible flows. By using the primitive variables of velocity and pressure, the Navier-Stokes and mass conservation equations are solved by a semi-implicit finite difference projection method. A new bounded higher order upwind convection scheme is employed to deal with the non-linear (advective) terms. The procedure is an adaptation of the GENSMAC (J. Comput. Phys. 1994; 110: 171-186) methodology for calculating confined and free surface fluid flows at both low and high Reynolds numbers. The calculations were performed by using the 2D version of the Freeflow simulation system (J. Comp. Visual. Science 2000; 2:199-210). In order to demonstrate the capabilities of the numerical method, various test cases are presented. These are the fully developed flow in a channel, the flow over a backward facing step, the die-swell problem, the broken dam flow, and an impinging jet onto a flat plate. The numerical results compare favourably with the experimental data and the analytical solutions. Copyright (c) 2006 John Wiley & Sons, Ltd.

Dynamics of meso and thermo citrate synthases with implicit solvation

The dynamics of hydration of meso and thermo citrate synthases has been investigated using the EEF1 methodology implemented with the CHARNM program. The native enzymes are composed of two identical subunits, each divided into a small and large domain. The dynamics behavior of both enzymes at 30 degrees C and 60 degrees C has been compared. The results of simulations show that during the hydration process, each subunit follows a different pathway of hydration, in spite of the identical sequence. The hydrated structures were compared with the crystalline structure, and the root mean square deviation (RMSD) of each residue along the trajectory was calculated. The regions with larger and smaller mobility were identified. In particular, helices belonging to the small domain are more mobile than those of the large domain. In contrast, the residues that constitute the active site show a much lower displacement compared with the crystalline structure. Hydration free energy calculations point out that Thermoplasma acidophilum citrate synthase (TCS) is more stable than chicken citrate synthase (CCS), at high temperatures. Such result has been ascribed to the higher number of superficial charges in the thermophilic homologue, which stabilizes the enzyme, while the mesophilic homologue denatures. These results are in accord with the experimental found that TCS keeps activity at temperatures farther apart from the catalysis regular temperature than the CCS. (c) 2005 Wiley Periodicals, Inc.

Differential geometry of intersection curves in R(4) of three implicit surfaces

We present algorithms for computing the differential geometry properties of intersection Curves of three implicit surfaces in R(4), using the implicit function theorem and generalizing the method of X. Ye and T. Maekawa for 4-dimension. We derive t, n, b(1), b(2) vectors and curvatures (k(1), k(2), k(3)) for transversal intersections of the intersection problem. (C) 2008 Elsevier B.V. All rights reserved.

Cell enumeration and visualisation by transmission electron microscopy of Lactobacillus rhamnosus treated with cinnamon (Cinnamomum zeylanicum B.) essential oil

The use of essential oils (EOs) in functional foods containing probiotic microorganisms must consider the antimicrobial activity of these oils against beneficial bacteria such as Lactobacillus rhamnosus. This study aimed to evaluate the sensitivity of L. rhamnosus cultures treated with cinnamon EO through viable cell counts and visualisation by transmission electron microscopy. Cinnamon EO at a concentration of 0.04% had a bacteriostatic activity after 2 h of incubation. Although slight alterations were detected in the cell structure, this concentration was considered to be bactericidal, since it led to a significant reduction in cell numbers after 24 h. on the other hand, cinnamon EO at a 1.00% concentration decreased cell counts by 3 log units after 2 h incubation and no viable cell count was detected after 24 h. Transmission electron microscopy indicated that cells treated with 1.00% cinnamon EO were severely damaged and presented cell membrane disruption and cytoplasmic leakage.

Implicit differential equations with impasse singularities and singular perturbation problems

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

Extension of the Nambu-Jona-Lasinio model predictions at high densities and temperatures using an implicit regularization scheme

Traditional cutoff regularization schemes of the Nambu-Jona-Lasinio model limit the applicability of the model to energy-momentum scales much below the value of the regularizing cutoff. In particular, the model cannot be used to study quark matter with Fermi momenta larger than the cutoff. In the present work, an extension of the model to high temperatures and densities recently proposed by Casalbuoni, Gatto, Nardulli, and Ruggieri is used in connection with an implicit regularization scheme. This is done by making use of scaling relations of the divergent one-loop integrals that relate these integrals at different energy-momentum scales. Fixing the pion decay constant at the chiral symmetry breaking scale in the vacuum, the scaling relations predict a running coupling constant that decreases as the regularization scale increases, implementing in a schematic way the property of asymptotic freedom of quantum chromodynamics. If the regularization scale is allowed to increase with density and temperature, the coupling will decrease with density and temperature, extending in this way the applicability of the model to high densities and temperatures. These results are obtained without specifying an explicit regularization. As an illustration of the formalism, numerical results are obtained for the finite density and finite temperature quark condensate and applied to the problem of color superconductivity at high quark densities and finite temperature.

A generalized alternating-direction implicit scheme for incompressible magnetohydrodynamic viscous flows at low magnetic Reynolds number

This paper presents numerical simulations of incompressible fluid flows in the presence of a magnetic field at low magnetic Reynolds number. The equations governing the flow are the Navier-Stokes equations of fluid motion coupled with Maxwell's equations of electromagnetics. The study of fluid flows under the influence of a magnetic field and with no free electric charges or electric fields is known as magnetohydrodynamics. The magnetohydrodynamics approximation is considered for the formulation of the non-dimensional problem and for the characterization of similarity parameters. A finite-difference technique is used to discretize the equations. In particular, an extension of the generalized Peaceman and Rachford alternating-direction implicit (ADI) scheme for simulating two-dimensional fluid flows is presented. The discretized conservation equations are solved in stream function-vorticity formulation. We compare the ADI and generalized ADI schemes, and show that the latter is more efficient in simulating low Reynolds number and magnetic Reynolds number problems. Numerical results demonstrating the applicability of this technique are also presented. The simulation of incompressible magneto hydrodynamic fluid flows is illustrated by numerical solution for two-dimensional cases. (c) 2007 Elsevier B.V. All rights reserved.

An implicit BEM formulation to model strong discontinuities in solids

A boundary element method (BEM) formulation to predict the behavior of solids exhibiting displacement (strong) discontinuity is presented. In this formulation, the effects of the displacement jump of a discontinuity interface embedded in an internal cell are reproduced by an equivalent strain field over the cell. To compute the stresses, this equivalent strain field is assumed as the inelastic part of the total strain. As a consequence, the non-linear BEM integral equations that result from the proposed approach are similar to those of the implicit BEM based on initial strains. Since discontinuity interfaces can be introduced inside the cell independently on the cell boundaries, the proposed BEM formulation, combined with a tracking scheme to trace the discontinuity path during the analysis, allows for arbitrary discontinuity propagation using a fixed mesh. A simple technique to track the crack path is outlined. This technique is based on the construction of a polygonal line formed by segments inside the cells, in which the assumed failure criterion is reached. Two experimental concrete fracture tests were analyzed to assess the performance of the proposed formulation.

A homological version of the implicit function theorem

Our objective in this paper is to prove an Implicit Function Theorem for general topological spaces. As a consequence, we show that, under certain conditions, the set of the invertible elements of a topological monoid X is an open topological group in X and we use the classical topological group theory to conclude that this set is a Lie group.