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.

Global perturbation of the marine Ca isotopic composition in the aftermath of the Marinoan global glaciation

Ca isotopic compositions of Marinoan post-glacial carbonate successions in Brazil and NW Canada were measured Both basal dolostones display delta(44/40)Ca values between 1 and 0 7 parts per thousand overlying limestones show a negative Ca isotope excursion to values around 0 1 parts per thousand and delta(44/40)Ca values rapidly increase up-section to near 2 0 parts per thousand In the Brazilian successions those high delta(44/40)Ca values rapidly decrease and stabilize to values between 0 6 and 0 9 parts per thousand These Ca isotope secular variation trends are unlike those of Sturtian post-glacial carbonate successions but similar to those of Marinoan post-glacial carbonate successions in Namibia suggesting that the perturbation of the marine Ca cycle was global This recommends Ca isotope stratigraphy as a tool to correlate Neoproterozoic post-glacial carbonate successions worldwide While the lowermost and uppermost strata have delta(44/40)Ca values typical of Phanerozoic carbonates the extremes 0 1 and 2 0 parts per thousand have not been thus far reported for other marine carbonates These extreme values suggest a short-lived non-actualistic perturbation in the marine Ca cycle Simple box modelling of the Marinoan post-glacial marine Ca cycle can reproduce the extreme values only by postulating a two-step process with Ca input initially exceeding Ca removal trough carbonate precipitation followed by precipitation overtaking a decreased Ca Input (C) 2010 Elsevier B V All rights reserved