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.

Surface Activity of the Triflate Ion at the Air/Water Interface and Properties of N,N,N-Trimethyl-N-Dodecylammonium Triflate Aqueous Solutions

The surface activity of salts added to water is Air orders of magnitude lower than that of surfactants. Sodium trifluoromethanesulfonate (NaTf) produced a change in surface tension. with concentration, Delta gamma/Delta c, of -13.2 mN.L/m.mol. This value is ca. 4-fold larger than those of simple salts and that of methanesulfonate. This unexpected surface effect suggested that positively charged micelles containing Tf could exhibit interesting properties. Dodecyltrimethylammonium triflate (DTATf) had a higher Kraft temperature (37 degrees C) and a lower cmc (5 x 10(-3)M) and degree of dissociation (0.11) than the chloride and bromide salts of DTA. Above the Kraft temperature, at a characteristic temperature t(1), the addition of NaTf above 0.05 M. to a DTATf solution induced phase separation. By increasing the temperature of the two-phase system to above t(1), a homogeneous, transparent solution was obtained at a characteristic temperature t(2). These results, together with well-known triflate properties, led us to suggest that the Tf ion pairs With DTA and that the -CF(3) group may be dehydrated in the interfacial region, resulting in new and interesting self-aggregated structures.

Analysis and Evaluation of EndpointSecurity Solutions

The main objective for this degree project was to analyze the Endpoint Security Solutions developed by Cisco, Microsoft and a third minor company solution represented by InfoExpress. The different solutions proposed are Cisco Network Admission Control, Microsoft Network Access Protection and InfoExpress CyberGatekeeper. An explanation of each solution functioning is proposed as well as an analysis of the differences between those solutions. This thesis work also proposes a tutorial for the installation of Cisco Network Admission Control for an easier implementation. The research was done by reading articles on the internet and by experimenting the Cisco Network Admission Control solution. My background knowledge about Cisco routing and ACL was also used. Based on the actual analysis done in this thesis, a conclusion was drawn that all existing solutions are not yet ready for large-scale use in corporate networks. Moreover all solutions are proprietary and incompatible. The future possible standard for Endpoint solution might be driven by Cisco and Microsoft and a rude competition begins between those two giants.

From XML to relational view updates: applying old solutions to solve a new problem

XML has become an important medium for data exchange, and is frequently used as an interface to - i.e. a view of - a relational database. Although lots of work have been done on querying relational databases through XML views, the problem of updating relational databases through XML views has not received much attention. In this work, we give the rst steps towards solving this problem. Using query trees to capture the notions of selection, projection, nesting, grouping, and heterogeneous sets found throughout most XML query languages, we show how XML views expressed using query trees can be mapped to a set of corresponding relational views. Thus, we transform the problem of updating relational databases through XML views into a classical problem of updating relational databases through relational views. We then show how updates on the XML view are mapped to updates on the corresponding relational views. Existing work on updating relational views can then be leveraged to determine whether or not the relational views are updatable with respect to the relational updates, and if so, to translate the updates to the underlying relational database. Since query trees are a formal characterization of view de nition queries, they are not well suited for end-users. We then investigate how a subset of XQuery can be used as a top level language, and show how query trees can be used as an intermediate representation of view de nitions expressed in this subset.

Contractual solutions to the holdup problem: a survey

In this survey, we presented the general idea and main results from what we understand that are the most important contributions to contractual solutions to the holdup problem literature. The aim of this paper is to push the previous analysis, uniform the notation and provide a snapshot on the most recent literature, as well as bring topics for future inquires on this issue.

Two field BPS solutions for generalized Lorentz breaking models

In this work we present nonlinear models in two-dimensional space-time of two interacting scalar fields in the Lorentz and CPT violating scenarios. We discuss the soliton solutions for these models as well as the question of stability for them. This is done by generalizing a model recently published by Barreto and collaborators and also by getting new solutions for the model introduced by them.

Bounded solutions for nonconserving-parity pseudoscalar potentials

The Dirac equation is analyzed for nonconserving-parity pseudoscalar radial potentials in 3+1 dimensions. It is shown that despite the nonconservation of parity this general problem can be reduced to a Sturm-Liouville problem of nonrelativistic fermions in spherically symmetric effective potentials. The searching for bounded solutions is done for the power-law and Yukawa potentials. The use of the methodology of effective potentials allow us to conclude that the existence of bound-state solutions depends whether the potential leads to a definite effective potential-well structure or to an effective potential less singular than -1/4r(2).

Bounded solutions of neutral fermions with a screened Coulomb potential

The intrinsically relativistic problem of a fermion subject to a pseudoscalar screened Coulomb plus a uniform background potential in two-dimensional space-time is mapped into a Sturm-Liouville. This mapping gives rise to an effective Morse-like potential and exact bounded solutions are found. It is shown that the uniform background potential determinates the number of bound-state solutions. The behaviour of the eigenenergies as well as of the upper and lower components of the Dirac spinor corresponding to bounded solutions is discussed in detail and some unusual results are revealed. An apparent paradox concerning the uncertainty principle is solved by recurring to the concepts of effective mass and effective Compton wavelength. (c) 2005 Elsevier B.V. All rights reserved.

Finite-well potential in the 3D nonlinear Schrodinger equation: application to Bose-Einstein condensation

Using variational and numerical solutions we show that stationary negative-energy localized (normalizable) bound states can appear in the three-dimensional nonlinear Schrodinger equation with a finite square-well potential for a range of nonlinearity parameters. Below a critical attractive nonlinearity, the system becomes unstable and experiences collapse. Above a limiting repulsive nonlinearity, the system becomes highly repulsive and cannot be bound. The system also allows nonnormalizable states of infinite norm at positive energies in the continuum. The normalizable negative-energy bound states could be created in BECs and studied in the laboratory with present knowhow.

Freezing of the QCD coupling constant and solutions of Schwinger-Dyson equations

We compare phenomenological values of the frozen QCD running coupling constant (alpha(s)) with two classes of infrared finite solutions obtained through nonperturbative Schwinger-Dyson equations. We use these same solutions with frozen coupling constants as well as their respective nonperturbative gluon propagators to compute the QCD prediction for the asymptotic pion form factor. Agreement between theory and experiment on alpha(s)(0) and F (pi)(Q(2)) is found only for one of the Schwinger-Dyson equation solutions.

A class of soliton solutions for the N=2 super mKdV/Sinh-Gordon hierarchy

Employing Hirota's method, a class of soliton solutions for the N = 2 super mKdV equations is proposed in terms of a single Grassmann parameter. Such solutions are shown to satisfy two copies of N = 1 supersymmetric mKdV equations connected by nontrivial algebraic identities. Using the super Miura transformation, we obtain solutions of the N = 2 super KdV equations. These are shown to generalize solutions derived previously. By using them KdV/sinh-Gordon hierarchy properties we generate the solutions of the N = 2 super sinh-Gordon as well.

Numerical and variational solutions of the dipolar Gross-Pitaevskii equation in reduced dimensions

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

Solutions of the bound-state Faddeev-Yakubovsky equations in three dimensions by using NN and 3N potential models

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

Lower bounds on blow up solutions of the three-dimensional Navier-Stokes equations in homogeneous Sobolev spaces

Suppose that u(t) is a solution of the three-dimensional Navier-Stokes equations, either on the whole space or with periodic boundary conditions, that has a singularity at time T. In this paper we show that the norm of u(T - t) in the homogeneous Sobolev space (H)over dot(s) must be bounded below by c(s)t(-(2s-1)/4) for 1/2 < s < 5/2 (s not equal 3/2), where c(s) is an absolute constant depending only on s; and by c(s)parallel to u(0)parallel to((5-2s)/5)(L2)t(-2s/5) for s > 5/2. (The result for 1/2 < s < 3/2 follows from well-known lower bounds on blowup in Lp spaces.) We show in particular that the local existence time in (H)over dot(s)(R-3) depends only on the (H)over dot(s)-norm for 1/2 < s < 5/2, s not equal 3/2. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4762841]

Photo-electrochemical investigation of anodic oxide films on cast Ti-Mo alloys I. Anodic behaviour and effect of alloy composition

The anodic behaviour of cast Ti-Mo alloys, having different Mo contents (6-20 wt.%), was investigated in acidic and neutral aerated aqueous solutions. All sample showed a valve-metal behaviour, owing to formation and thickening of barrier-type anodic oxides displaying interference colours Growth kinetics. of passive films is influenced by both anodizing electrolyte and composition of the starting alloy. This last parameter was found to change also the solid-state properties of the films, explored by photoelectrochemical and impedance spectroscopy experiments. Thicker films (U(f) = 8 V/MSE) grown on alloys richer in Mo showed more resistive character and a photocurrent sign inversion under negative bias, that revealed an insulating character, whereas corresponding films grown on alloys with lower Mo content, as well as thinner films, behaved as n-type semiconductors. Results are discussed in terms of formation of a mixed Ti-Mo oxide phase. (C) 2008 Elsevier Ltd. All rights reserved