Numerical assessment of stability of interface discontinuous finite element pressure spaces

The stability of two recently developed pressure spaces has been assessed numerically: The space proposed by Ausas et al. [R.F. Ausas, F.S. Sousa, G.C. Buscaglia, An improved finite element space for discontinuous pressures, Comput. Methods Appl. Mech. Engrg. 199 (2010) 1019-1031], which is capable of representing discontinuous pressures, and the space proposed by Coppola-Owen and Codina [A.H. Coppola-Owen, R. Codina, Improving Eulerian two-phase flow finite element approximation with discontinuous gradient pressure shape functions, Int. J. Numer. Methods Fluids, 49 (2005) 1287-1304], which can represent discontinuities in pressure gradients. We assess the stability of these spaces by numerically computing the inf-sup constants of several meshes. The inf-sup constant results as the solution of a generalized eigenvalue problems. Both spaces are in this way confirmed to be stable in their original form. An application of the same numerical assessment tool to the stabilized equal-order P-1/P-1 formulation is then reported. An interesting finding is that the stabilization coefficient can be safely set to zero in an arbitrary band of elements without compromising the formulation's stability. An analogous result is also reported for the mini-element P-1(+)/P-1 when the velocity bubbles are removed in an arbitrary band of elements. (C) 2012 Elsevier B.V. All rights reserved.

SCG postnatal remodelling - hypertrophy and neuron number stability - in Spix's Yellow-toothed Cavies (Galea spixii)

Whilst a fall in neuron numbers seems a common pattern during postnatal development, several authors have nonetheless reported an increase in neuron number, which may be associated with any one of a number of possible processes encapsulating either neurogenesis or late maturation and incomplete differentiation. Recent publications have thus added further fuel to the notion that a postnatal neurogenesis may indeed exist in sympathetic ganglia. In the light of these uncertainties surrounding the effects exerted by postnatal development on the number of superior cervical ganglion (SCG) neurons, we have used state-of-the-art design-based stereology to investigate the quantitative structure of SCG at four distinct timepoints after birth, viz., 1-3 days, 1 month, 12 months and 36 months. The main effects exerted by ageing on the SCG structure were: (i) a 77% increase in ganglion volume; (ii) stability in the total number of the whole population of SCG nerve cells (no change - either increase or decrease) during post-natal development; (iii) a higher proportion of uninucleate neurons to binucleate neurons only in newborn animals; (iv) a 130% increase in the volume of uninucleate cell bodies; and (v) the presence of BrdU positive neurons in animals at all ages. At the time of writing our results support the idea that neurogenesis takes place in the SCG of preas, albeit it warrants confirmation by further markers. We also hypothesise that a portfolio of other mechanisms: cell repair, maturation, differentiation and death may be equally intertwined and implicated in the numerical stability of SCG neurons during postnatal development. (C) 2011 ISDN. Published by Elsevier Ltd. All rights reserved.

GENERALIZED FINITE ELEMENT METHOD ON NONCONVENTIONAL HYBRID-MIXED FORMULATION

The generalized finite element method (GFEM) is applied to a nonconventional hybrid-mixed stress formulation (HMSF) for plane analysis. In the HMSF, three approximation fields are involved: stresses and displacements in the domain and displacement fields on the static boundary. The GFEM-HMSF shape functions are then generated by the product of a partition of unity associated to each field and the polynomials enrichment functions. In principle, the enrichment can be conducted independently over each of the HMSF approximation fields. However, stability and convergence features of the resulting numerical method can be affected mainly by spurious modes generated when enrichment is arbitrarily applied to the displacement fields. With the aim to efficiently explore the enrichment possibilities, an extension to GFEM-HMSF of the conventional Zienkiewicz-Patch-Test is proposed as a necessary condition to ensure numerical stability. Finally, once the extended Patch-Test is satisfied, some numerical analyses focusing on the selective enrichment over distorted meshes formed by bilinear quadrilateral finite elements are presented, thus showing the performance of the GFEM-HMSF combination.

A regional multirnodulus algorithm for blind equalization of QAM signals: Introduction and steady-state analysis

It is well known that constant-modulus-based algorithms present a large mean-square error for high-order quadrature amplitude modulation (QAM) signals, which may damage the switching to decision-directed-based algorithms. In this paper, we introduce a regional multimodulus algorithm for blind equalization of QAM signals that performs similar to the supervised normalized least-mean-squares (NLMS) algorithm, independently of the QAM order. We find a theoretical relation between the coefficient vector of the proposed algorithm and the Wiener solution and also provide theoretical models for the steady-state excess mean-square error in a nonstationary environment. The proposed algorithm in conjunction with strategies to speed up its convergence and to avoid divergence can bypass the switching mechanism between the blind mode and the decision-directed mode. (c) 2012 Elsevier B.V. All rights reserved.

MHD numerical simulations of colliding winds in massive binary systems - I. Thermal versus non-thermal radio emission

In the past few decades detailed observations of radio and X-ray emission from massive binary systems revealed a whole new physics present in such systems. Both thermal and non-thermal components of this emission indicate that most of the radiation at these bands originates in shocks. O and B-type stars and WolfRayet (WR) stars present supersonic and massive winds that, when colliding, emit largely due to the freefree radiation. The non-thermal radio and X-ray emissions are due to synchrotron and inverse Compton processes, respectively. In this case, magnetic fields are expected to play an important role in the emission distribution. In the past few years the modelling of the freefree and synchrotron emissions from massive binary systems have been based on purely hydrodynamical simulations, and ad hoc assumptions regarding the distribution of magnetic energy and the field geometry. In this work we provide the first full magnetohydrodynamic numerical simulations of windwind collision in massive binary systems. We study the freefree emission characterizing its dependence on the stellar and orbital parameters. We also study self-consistently the evolution of the magnetic field at the shock region, obtaining also the synchrotron energy distribution integrated along different lines of sight. We show that the magnetic field in the shocks is larger than that obtained when the proportionality between B and the plasma density is assumed. Also, we show that the role of the synchrotron emission relative to the total radio emission has been underestimated.

Architectures, stability and optimization for clock distribution networks

Synchronous telecommunication networks, distributed control systems and integrated circuits have its accuracy of operation dependent on the existence of a reliable time basis signal extracted from the line data stream and acquirable to each node. In this sense, the existence of a sub-network (inside the main network) dedicated to the distribution of the clock signals is crucially important. There are different solutions for the architecture of the time distribution sub-network and choosing one of them depends on cost, precision, reliability and operational security. In this work we expose: (i) the possible time distribution networks and their usual topologies and arrangements. (ii) How parameters of the network nodes can affect the reachability and stability of the synchronous state of a network. (iii) Optimizations methods for synchronous networks which can provide low cost architectures with operational precision, reliability and security. (C) 2011 Elsevier B. V. All rights reserved.

Modeling the emergence of HIV-1 drug resistance resulting from antiretroviral therapy: Insights from theoretical and numerical studies

The use of antiretroviral therapy has proven to be remarkably effective in controlling the progression of human immunodeficiency virus (HIV) infection and prolonging patient's survival. Therapy however may fail and therefore these benefits can be compromised by the emergence of HIV strains that are resistant to the therapy. In view of these facts, the question of finding the reason for which drug-resistant strains emerge during therapy has become a worldwide problem of great interest. This paper presents a deterministic HIV-1 model to examine the mechanisms underlying the emergence of drug-resistance during therapy. The aim of this study is to determine whether, and how fast, antiretroviral therapy may determine the emergence of drug resistance by calculating the basic reproductive numbers. The existence, feasibility and local stability of the equilibriums are also analyzed. By performing numerical simulations we show that Hopf bifurcation may occur. The model suggests that the individuals with drug-resistant infection may play an important role in the epidemic of HIV. (C) 2011 Elsevier Ireland Ltd. All rights reserved.

Stability analysis and area spectrum of three-dimensional Lifshitz black holes

In this work, we probe the stability of a z = 3 three-dimensional Lifshitz black hole by using scalar and spinorial perturbations. We found an analytical expression for the quasinormal frequencies of the scalar probe field, which perfectly agree with the behavior of the quasinormal modes obtained numerically. The results for the numerical analysis of the spinorial perturbations reinforce the conclusion of the scalar analysis, i.e., the model is stable under scalar and spinor perturbations. As an application we found the area spectrum of the Lifshitz black hole, which turns out to be equally spaced.

Structural stability of flexible lines in catenary configuration under torsion

Catenary risers can present during installation a very low tension close to seabed, which combined with torsion moment can lead to a structural instability, resulting in a loop. This is undesirable once it is possible that the loop turns into a kink, creating damage. This work presents a numerical methodology to analyze the conditions of loop formation in catenary risers. Stability criteria were applied to finite element models, including geometric nonlinearities and contact constraint due to riser-seabed interaction. The classical Greenhill's formula was used to predict the phenomenon and parametric analysis shows a “universal plot” able to predict instability in catenaries using a simple equation that can be applied for typical risers installation conditions and, generically, for catenary lines under torsion.