Iterative strong coupling of dimensionally heterogeneous models

In this article we address decomposition strategies especially tailored to perform strong coupling of dimensionally heterogeneous models, under the hypothesis that one wants to solve each submodel separately and implement the interaction between subdomains by boundary conditions alone. The novel methodology takes full advantage of the small number of interface unknowns in this kind of problems. Existing algorithms can be viewed as variants of the `natural` staggered algorithm in which each domain transfers function values to the other, and receives fluxes (or forces), and vice versa. This natural algorithm is known as Dirichlet-to-Neumann in the Domain Decomposition literature. Essentially, we propose a framework in which this algorithm is equivalent to applying Gauss-Seidel iterations to a suitably defined (linear or nonlinear) system of equations. It is then immediate to switch to other iterative solvers such as GMRES or other Krylov-based method. which we assess through numerical experiments showing the significant gain that can be achieved. indeed. the benefit is that an extremely flexible, automatic coupling strategy can be developed, which in addition leads to iterative procedures that are parameter-free and rapidly converging. Further, in linear problems they have the finite termination property. Copyright (C) 2009 John Wiley & Sons, Ltd.

D (sJ) (2317) meson production at RHIC

We study the production of D (sJ) (2317) mesons in relativistic heavy ion collisions using the quark coalescence model. The predicted D (sJ) (2317) abundance depends sensitively on the quark structure of the D (sJ) (2317) meson. We have also evaluated the absorption cross sections of the D (sJ) (2317) meson by pi, rho, kaon and K* in a phenomenological hadronic model. We find that the final yield of D (sJ) (2317) mesons remains sensitive to its initial number produced from the quark-gluon plasma, providing thus the possibility of studying the quark structure of the D (sJ) (2317) meson and its production mechanism in relativistic heavy ion collisions.

EXACT AND QUASI-EXACT MODELS OF STRANGE STARS

We construct and compare in this work a variety of simple models for strange stars, namely, hypothetical self-bound objects made of a cold stable version of the quark-gluon plasma. Exact, quasi-exact and numerical models are examined to find the most economical description for these objects. A simple and successful parametrization of them is given in terms of the central density, and the differences among the models are explicitly shown and discussed. In particular, we present a model starting with a Gaussian ansatz for the density profile that provides a very accurate and almost complete analytical integration of the problem, modulo a small difference for one of the metric potentials.

NU-PROCESS IN EXOTIC MODELS

The exact physical conditions generating the abundances of r-elements in environments such as supernovae explosions are still under debate. We evaluated the characteristics expected for the neutrino wind in the proposed model of type-II supernova driven by conversion of nuclear matter to strange matter. Neutrinos will change the final abundance of elements after freeze out of r-process nucleosynthesis, specially those close to mass peaks.

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.

A coupling cutting stock-lot sizing problem in the paper industry

An important production programming problem arises in paper industries coupling multiple machine scheduling with cutting stocks. Concerning machine scheduling: how can the production of the quantity of large rolls of paper of different types be determined. These rolls are cut to meet demand of items. Scheduling that minimizes setups and production costs may produce rolls which may increase waste in the cutting process. On the other hand, the best number of rolls in the point of view of minimizing waste may lead to high setup costs. In this paper, coupled modeling and heuristic methods are proposed. Computational experiments are presented.

Exclusive processes in electron-ion collisions

The exclusive processes in electron-ion (eA) interactions are an important tool to investigate the QCD dynamics at high energies as they are in general driven by the gluon content of the target which is strongly subject to parton saturation effects. In this Letter we compute the cross sections for the exclusive vector meson production as well as the deeply virtual Compton scattering (DVCS) relying on the color dipole approach and considering the numerical solution of the Balitsky-Kovchegov equation including running coupling corrections (rcBK). The production cross sections obtained with the rcBK solution and bCGC parametrization are very similar, the former being slightly larger. (C) 2011 Elsevier B.V. All rights reserved.

QCD sum rules study of the meson Z(+)(4430)

We use QCD sum rules to study the recently observed meson Z(+)(4430), considered as a D*D-1 molecule with J(P) = 0(-). We consider the contributions of condensates up to dimension eight and work at leading order in alpha(s). We get m(Z) = (4.40 +/- 0.10) GeV in a very good agreement with the experimental value. We also make predictions for the analogous mesons Z(s) and Z(bb) considered as D-s*D-1 and B*B-1 molecules, respectively. For Z(s) we predict mZ(s) = (4.70 +/- 0.06) GeV, which is above the D-s* D-1 threshold, indicating that it is probably a very broad state and, therefore, difficult to observe experimentally. For Z(bb) we predict m(Zbb) = (10.74 +/- 0.12) GeV, in agreement with quark model predictions. (c) 2008 Elsevier B.V. All rights reserved.

Stability of the curvature perturbation in dark sectors` mutual interacting models

We consider perturbations in a cosmological model with a small coupling between dark energy and dark matter. We prove that the stability of the curvature perturbation depends on the type of coupling between dark sectors. When the dark energy is of quintessence type, if the coupling is proportional to the dark matter energy density, it will drive the instability in the curvature perturbations: however if the coupling is proportional to the energy density of dark energy, there is room for the stability in the curvature perturbations. When the dark energy is of phantom type, the perturbations are always stable, no matter whether the coupling is proportional to the one or the other energy density. (C) 2008 Elsevier B.V. All rights reserved.

Observational constraints on the dark energy and dark matter mutual coupling

We examine different phenomenological interaction models for Dark Energy and Dark Matter by performing statistical joint analysis with observational data arising from the 182 Gold type la supernova samples, the shift parameter of the Cosmic Microwave Background given by the three-year Wilkinson Microwave Anisotropy Probe observations, the baryon acoustic oscillation measurement from the Sloan Digital Sky Survey and age estimates of 35 galaxies. Including the time-dependent observable, we add sensitivity of measurement and give complementary results for the fitting. The compatibility among three different data sets seem to imply that the coupling between dark energy and dark matter is a small positive value, which satisfies the requirement to solve the coincidence problem and the second law of thermodynamics, being compatible with previous estimates. (c) 2008 Elsevier B.V. All rights reserved.

Energy and system size dependence of phi meson production in Cu plus Cu and Au plus Au collisions

We study the beam-energy and system-size dependence of phi meson production (using the hadronic decay mode phi -> K(+) K(-)) by comparing the new results from Cu + Cu collisions and previously reported Au + Au collisions at root s(NN) = 62.4 and 200 GeV measured in the STAR experiment at RHIC. Data presented in this Letter are from mid-rapidity (vertical bar y vertical bar < 0.5) for 0.4 < p(T) < 5 GeV/c. At a given beam energy, the transverse momentum distributions for phi mesons are observed to be similar in yield and shape for Cu + Cu and Au + Au colliding systems with similar average numbers of participating nucleons. The phi meson yields in nucleus-nucleus collisions, normalized by the average number of participating nucleons, are found to be enhanced relative to those from p + p collisions. The enhancement for phi mesons lies between strange hadrons having net strangeness = 1 (K(-) and <(A)over bar>) and net strangeness = 2 (Xi). The enhancement for phi mesons is observed to be higher at root s(NN) = 200 GeV compared to 62.4 GeV. These observations for the produced phi(s (s) over bar) mesons clearly suggest that, at these collision energies, the source of enhancement of strange hadrons is related to the formation of a dense partonic medium in high energy nucleus-nucleus collisions and cannot be alone due to canonical suppression of their production in smaller systems. (C) 2009 Elsevier B.V. All rights reserved.