Black-box decomposition approach for computational hemodynamics: One-dimensional models

Increasing efforts exist in integrating different levels of detail in models of the cardiovascular system. For instance, one-dimensional representations are employed to model the systemic circulation. In this context, effective and black-box-type decomposition strategies for one-dimensional networks are needed, so as to: (i) employ domain decomposition strategies for large systemic models (1D-1D coupling) and (ii) provide the conceptual basis for dimensionally-heterogeneous representations (1D-3D coupling, among various possibilities). The strategy proposed in this article works for both of these two scenarios, though the several applications shown to illustrate its performance focus on the 1D-1D coupling case. A one-dimensional network is decomposed in such a way that each coupling point connects two (and not more) of the sub-networks. At each of the M connection points two unknowns are defined: the flow rate and pressure. These 2M unknowns are determined by 2M equations, since each sub-network provides one (non-linear) equation per coupling point. It is shown how to build the 2M x 2M non-linear system with arbitrary and independent choice of boundary conditions for each of the sub-networks. The idea is then to solve this non-linear system until convergence, which guarantees strong coupling of the complete network. In other words, if the non-linear solver converges at each time step, the solution coincides with what would be obtained by monolithically modeling the whole network. The decomposition thus imposes no stability restriction on the choice of the time step size. Effective iterative strategies for the non-linear system that preserve the black-box character of the decomposition are then explored. Several variants of matrix-free Broyden`s and Newton-GMRES algorithms are assessed as numerical solvers by comparing their performance on sub-critical wave propagation problems which range from academic test cases to realistic cardiovascular applications. A specific variant of Broyden`s algorithm is identified and recommended on the basis of its computer cost and reliability. (C) 2010 Elsevier B.V. All rights reserved.

An improved finite element space for discontinuous pressures

We consider incompressible Stokes flow with an internal interface at which the pressure is discontinuous, as happens for example in problems involving surface tension. We assume that the mesh does not follow the interface, which makes classical interpolation spaces to yield suboptimal convergence rates (typically, the interpolation error in the L(2)(Omega)-norm is of order h(1/2)). We propose a modification of the P(1)-conforming space that accommodates discontinuities at the interface without introducing additional degrees of freedom or modifying the sparsity pattern of the linear system. The unknowns are the pressure values at the vertices of the mesh and the basis functions are computed locally at each element, so that the implementation of the proposed space into existing codes is straightforward. With this modification, numerical tests show that the interpolation order improves to O(h(3/2)). The new pressure space is implemented for the stable P(1)(+)/P(1) mini-element discretization, and for the stabilized equal-order P(1)/P(1) discretization. Assessment is carried out for Poiseuille flow with a forcing surface and for a static bubble. In all cases the proposed pressure space leads to improved convergence orders and to more accurate results than the standard P(1) space. In addition, two Navier-Stokes simulations with moving interfaces (Rayleigh-Taylor instability and merging bubbles) are reported to show that the proposed space is robust enough to carry out realistic simulations. (c) 2009 Elsevier B.V. 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.

Continuously deformable topological structure

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Control aspects in nonlinear Hill's equation

The Hill's equations-even in the linear original version are a describer of phenomenon having chaotic flavor, giving sometimes very unusual situations. The theory of the so called intervals of instability in the equation provides the precise description for most of these phenomena. Considerations on nonlinearities into the Hill's equation is a quite recent task. The linearized version for almost of these systems it reduces to the Hill's classical linear one. In this paper, some indicative facts are pointed out on the possibility of having the linear system stabilizable and/or exactly controllable. As consequence of such an approach we get results having strong classical aspects, like the one talking about location of parameters in intervals of stability. A result for nonlinear proper periodic controls, is considered too. (C) 2010 Elsevier B.V. All rights reserved.

O sobre-sinal e a resposta inversa em sistemas de controle lineares contínuos no tempo

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

Short term hydroelectric scheduling combining network flow and interior point approaches

In this paper, short term hydroelectric scheduling is formulated as a network flow optimization model and solved by interior point methods. The primal-dual and predictor-corrector versions of such interior point methods are developed and the resulting matrix structure is explored. This structure leads to very fast iterations since it avoids computation and factorization of impedance matrices. For each time interval, the linear algebra reduces to the solution of two linear systems, either to the number of buses or to the number of independent loops. Either matrix is invariant and can be factored off-line. As a consequence of such matrix manipulations, a linear system which changes at each iteration has to be solved, although its size is reduced to the number of generating units and is not a function of time intervals. These methods were applied to IEEE and Brazilian power systems, and numerical results were obtained using a MATLAB implementation. Both interior point methods proved to be robust and achieved fast convergence for all instances tested. (C) 2004 Elsevier Ltd. All rights reserved.

Eventually minimal curves

A curve defined over a finite field is maximal or minimal according to whether the number of rational points attains the upper or the lower bound in Hasse-Weil's theorem, respectively. In the study of maximal curves a fundamental role is played by an invariant linear system introduced by Ruck and Stichtenoth in [6]. In this paper we define an analogous invariant system for minimal curves, and we compute its orders and its Weierstrass points. In the last section we treat the case of curves having genus three in characteristic two.

Analyses of the modifications in the pi circuits for inclusion of frequency influence in transmission line representation

In this article, it is represented by state variables phase a transmission line which parameters are considered frequency independently and frequency dependent. It is analyzed what is the reasonable number of pi circuits and the number of blocks composed by parallel resistor and inductor in parallel for reduction of numerical oscillations. It is simulated the numerical routine with and without the effect of frequency in the longitudinal parameters. Initially, it is used state variables and pi circuits representing the transmission line composing a linear system which is solved by numerical routines based on the trapezoidal rule. The effect of frequency on the line is synthesized by resistors and inductors in parallel and this representation is analyzed in details. It is described transmission lines and the frequency influence in these lines through the state variables.

A radar reflectivity-runoff model for use in flood warning

This model connects directly the radar reflectivity data and hydrological variable runoff. The catchment is discretized in pixels (4 Km × 4 Km) with the same resolution of the CAPPI. Careful discretization is made so that every grid catchment pixel corresponds precisely to CAPPI grid cell. The basin is assumed a linear system and also time invariant. The forecast technique takes advantage of spatial and temporal resolutions obtained by the radar. The method uses only the measurements of the factor reflectivity distribution observed over the catchment area without using the reflectivity - rainfall rate transformation by the conventional Z-R relationships. The reflectivity values in each catchment pixel are translated to a gauging station by using a transfer function. This transfer function represents the travel time of the superficial water flowing through pixels in the drainage direction ending at the gauging station. The parameters used to compute the transfer function are concentration time and the physiographic catchment characteristics. -from Authors

Nonmodal energetics of electromagnetic drift waves

The linear properties of an electromagnetic drift-wave model are examined. The linear system is non-normal in that its eigenvectors are not orthogonal with respect to the energy inner product. The non-normality of the linear evolution operator can lead to enhanced finite-time growth rates compared to modal growth rates. Previous work with an electrostatic drift-wave model found that nonmodal behavior is important in the hydrodynamic limit. Here, similar behavior is seen in the hydrodynamic regime even with the addition of magnetic fluctuations. However, unlike the results for the electrostatic drift-wave model, nonmodal behavior is also important in the adiabatic regime with moderate to strong magnetic fluctuations. © 2000 American Institute of Physics.

Constraining the free parameter of the high parton density effects

The high parton density effects are strongly dependent of the spatial gluon distribution within the proton, with radius R, which cannot be derived from perturbative QCD. In this paper we assume that the unitarity corrections are present in the HERA kinematical region and constrain the value of R using the data for the proton structure function and its slope. We obtain that the gluons are not distributed uniformly in the whole proton disc, but behave as concentrated in smaller regions. (C) 2000 Elsevier Science B.V.

Analytical/numerical hybrid solution for one-dimensional ablation problem

Ablation is a thermal protection process with several applications in engineering, mainly in the field of airspace industry. The use of conventional materials must be quite restricted, because they would suffer catastrophic flaws due to thermal degradation of their structures. However, the same materials can be quite suitable once being protected by well-known ablative materials. The process that involves the ablative phenomena is complex, could involve the whole or partial loss of material that is sacrificed for absorption of energy. The analysis of the ablative process in a blunt body with revolution geometry will be made on the stagnation point area that can be simplified as a one-dimensional plane plate problem, hi this work the Generalized Integral Transform Technique (GITT) is employed for the solution of the non-linear system of coupled partial differential equations that model the phenomena. The solution of the problem is obtained by transforming the non-linear partial differential equation system to a system of coupled first order ordinary differential equations and then solving it by using well-established numerical routines. The results of interest such as the temperature field, the depth and the rate of removal of the ablative material are presented and compared with those ones available in the open literature.

Voltammetric determination of dipyrone using a N,N′- ethylenebis(salicylideneaminato)oxovanadium(IV) modified carbon-paste electrode

The preparation and electrochemical characterization of a carbon paste electrode modified with N.N′-ethylenebis(salicylideneiminato) oxovanadium(IV) complex ([VO(Salen)]) as well as its behavior as electrocatalyst toward the oxidation of dipyrone were investigated. The electrochemical behavior of the modified electrode and the electrooxidation of dipyrone were explored using cyclic voltammetry. The voltammetric response of the modified electrode is based on two reactions. One electrochemical related to the oxidation of the metallic center of the [VO(Salen)] and the other involving the chemical redox process involving the oxidized form of the complex and the reduced form of dipyrone. The best voltammetric response was observed for a paste composition of 25% (m/m) [VO(Salen)], KCl solution pH from 5.5 to 8.0 as the electrolyte and potential scan rate of 10 mV s-1 in the presence of dipyrone. A linear voltammetric response for dipyrone was obtained in the concentration range from 9.9 × 106 to 2.8 × 10 -3 mol L-1, with a detection limit of 7.2 × 10 -6 mol L-1. Among of several compounds tested as potential interference, only ascorbic acid presented some interference. The proposed electrode is useful for the quality control and routine analysis of dipyrone in pharmaceutical formulations.

Substitution-newton-Raphson method for the solution of electric network equations

The power flow problem, in transmission networks, has been well solved, for most cases, using Newton-Raphson method (NR) and its decoupled versions. Generally speaking, the solution of a non-linear system of equations refers to two methods: NR and Successive Substitution. The proposal of this paper is to evaluate the potential of the Substitution-Newton-Raphson Method (SNR), which combines both methods, on the solution of the power flow problem. Simulations were performed using a two-bus test network in order to observe the characteristics of these methods. It was verified that the NR is faster than SNR, in terms of convergence, considering non-stressed scenarios. For those cases where the power flow in the network is closed to the limits (stressed system), the SNR converges faster. This paper presents the power flow formulation of the SNR and describes its potential for its application in special cases such as stressed scenarios. © 2006 IEEE.