Coupling model uncertainty for coupled rainfall/runoff and surface water quality models in river problems

Coupled hydrology and water quality models are an important tool today, used in the understanding and management of surface water and watershed areas. Such problems are generally subject to substantial uncertainty in parameters, process understanding, and data. Component models, drawing on different data, concepts, and structures, are affected differently by each of these uncertain elements. This paper proposes a framework wherein the response of component models to their respective uncertain elements can be quantified and assessed, using a hydrological model and water quality model as two exemplars. The resulting assessments can be used to identify model coupling strategies that permit more appropriate use and calibration of individual models, and a better overall coupled model response. One key finding was that an approximate balance of water quality and hydrological model responses can be obtained using both the QUAL2E and Mike11 water quality models. The balance point, however, does not support a particularly narrow surface response (or stringent calibration criteria) with respect to the water quality calibration data, at least in the case examined here. Additionally, it is clear from the results presented that the structural source of uncertainty is at least as significant as parameter-based uncertainties in areal models. © 2012 John Wiley & Sons, Ltd.

Search for maximal flavor violating scalars in same-charge lepton pairs in pp̄ collisions at s=1.96TeV

Models of maximal flavor violation (MxFV) in elementary particle physics may contain at least one new scalar SU(2) doublet field ΦFV=(η0,η+) that couples the first and third generation quarks (q1, q3) via a Lagrangian term LFV=ξ13ΦFVq1q3. These models have a distinctive signature of same-charge top-quark pairs and evade flavor-changing limits from meson mixing measurements. Data corresponding to 2fb-1 collected by the Collider Dectector at Fermilab II detector in pp̄ collisions at s=1.96TeV are analyzed for evidence of the MxFV signature. For a neutral scalar η0 with mη0=200GeV/c2 and coupling ξ13=1, ∼11 signal events are expected over a background of 2.1±1.8 events. Three events are observed in the data, consistent with background expectations, and limits are set on the coupling ξ13 for mη0=180-300GeV/c2. © 2009 The American Physical Society.

An explicit scheme for coupling temperature and concentration fields in solidification models

A numerical scheme for coupling temperature and concentration fields in a general solidification model is presented. A key feature of this scheme is an explicit time stepping used in solving the governing thermal and solute conservation equations. This explicit approach results in a local point-by-point coupling scheme for the temperature and concentration and avoids the multi-level iteration required by implicit time stepping schemes. The proposed scheme is validated by predicting the concentration field in a benchmark solidification problem. Results compare well with an available similarity solution. The simplicity of the proposed explicit scheme allows for the incorporation of complex microscale models into a general solidification model. This is demonstrated by investigating the role of dendrite coarsening on the concentration field in the solidification benchmark problem.

Mixed-dimensional coupling in finite element models

In many finite element analysis models it would be desirable to combine reduced- or lower-dimensional element types with higher-dimensional element types in a single model. In order to achieve compatibility of displacements and stress equilibrium at the junction or interface between the differing element types, it is important in such cases to integrate into the analysis some scheme for coupling the element types. A novel and effective scheme for establishing compatibility and equilibrium at the dimensional interface is described and its merits and capabilities are demonstrated. Copyright (C) 2000 John Wiley & Sons, Ltd.

Meson properties in an extended nonlocal NJL model

A nonlocal version of the NJL model is investigated. It is based on a separable quark-quark interaction, as suggested by the instanton liquid picture of the QCD vacuum. The interaction is extended to include terms that bind vector and axial-vector mesons. The nonlocality means that no further regulator is required. Moreover the model is able to confine the quarks by generating a quark propagator without poles at real energies. Features of the continuation of amplitudes from Euclidean space to Minkowski energies are discussed. These features lead to restrictions on the model parameters as well as on the range of applicability of the model. Conserved currents are constructed, and their consistency with various Ward identities is demonstrated. In particular, the Gell-Mann-Oakes-Renner relation is derived both in the ladder approximation and at meson loop level. The importance of maintaining chiral symmetry in the calculations is stressed throughout. Calculations with the model are performed to all orders in momentum. Meson masses are determined, along with their strong and electromagnetic decay amplitudes. Also calculated are the electromagnetic form factor of the pion and form factors associated with the processes gamma gamma* --> pi0 and omega --> pi0 gamma*. The results are found to lead to a satisfactory phenomenology and demonstrate a possible dynamical origin for vector-meson dominance. In addition, the results produced at meson loop level validate the use of 1/Nc as an expansion parameter and indicate that a light and broad scalar state is inherent in models of the NJL type.

Stratosphere‐troposphere coupling and annular mode variability in chemistry‐climate models

The internal variability and coupling between the stratosphere and troposphere in CCMVal‐2 chemistry‐climate models are evaluated through analysis of the annular mode patterns of variability. Computation of the annular modes in long data sets with secular trends requires refinement of the standard definition of the annular mode, and a more robust procedure that allows for slowly varying trends is established and verified. The spatial and temporal structure of the models’ annular modes is then compared with that of reanalyses. As a whole, the models capture the key features of observed intraseasonal variability, including the sharp vertical gradients in structure between stratosphere and troposphere, the asymmetries in the seasonal cycle between the Northern and Southern hemispheres, and the coupling between the polar stratospheric vortices and tropospheric midlatitude jets. It is also found that the annular mode variability changes little in time throughout simulations of the 21st century. There are, however, both common biases and significant differences in performance in the models. In the troposphere, the annular mode in models is generally too persistent, particularly in the Southern Hemisphere summer, a bias similar to that found in CMIP3 coupled climate models. In the stratosphere, the periods of peak variance and coupling with the troposphere are delayed by about a month in both hemispheres. The relationship between increased variability of the stratosphere and increased persistence in the troposphere suggests that some tropospheric biases may be related to stratospheric biases and that a well‐simulated stratosphere can improve simulation of tropospheric intraseasonal variability.

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.

(D)over-barN interaction from meson-exchange and quark-gluon dynamics

We investigate the (D) over barN interaction at low energies using a meson exchange model supplemented with a short-distance contribution from one-gluon exchange. The model is developed in close analogy to the meson-exchange KN interaction of the Julich group utilizing SU(4) symmetry constraints. The main ingredients of the interaction are provided by vector meson (rho, omega) exchange and higher-order box diagrams involving (D) over bar *N , (D) over bar Delta, and (D) over bar*Delta intermediate states. The short-range part is assumed to receive additional contributions from genuine quark-gluon processes. The predicted cross-sections for (D) over barN for excess energies up to 150MeV are of the same order of magnitude as those for KN but with average values of around 20mb, roughly a factor two larger than for the latter system. It is found that the omega-exchange plays a very important role. Its interference pattern with the rho-exchange, which is basically fixed by the assumed SU(4) symmetry, clearly determines the qualitative features of the (D) over barN interaction - very similiar to what happens also for the KN system.

Strength of the trilinear Higgs boson coupling in technicolor models

We discuss the strength of the trilinear Higgs boson coupling in technicolor (or composite) models in a model independent way. The coupling is determined as a function of a very general ansatz for the technicolor self-energy, and turns out to be equal or smaller than the one of the Standard Model Higgs boson depending on the dynamics of the theory. (c) 2006 Published by Elsevier B.V.