PARTITIONED ANALYSIS FOR DIMENSIONALLY-HETEROGENEOUS HYDRAULIC NETWORKS

In this work an iterative strategy is developed to tackle the problem of coupling dimensionally-heterogeneous models in the context of fluid mechanics. The procedure proposed here makes use of a reinterpretation of the original problem as a nonlinear interface problem for which classical nonlinear solvers can be applied. Strong coupling of the partitions is achieved while dealing with different codes for each partition, each code in black-box mode. The main application for which this procedure is envisaged arises when modeling hydraulic networks in which complex and simple subsystems are treated using detailed and simplified models, correspondingly. The potentialities and the performance of the strategy are assessed through several examples involving transient flows and complex network configurations.

Determination of the minimum value of the safety factor from geodesic Alfven eigenmodes in Joint European Torus

An analysis of the experimental conditions under which low-frequency (70-150 kHz) Alfven eigertmodes (AE) are excited during the monster sawtooth in Joint European Torus [F Romanelli et al, Proceedings of the 22nd IAEA Fusion Energy Conference, Geneva, Switzerland, 2008] is presented for the specific case of a discharge with ion cyclotron heating (5 MW) Using a simplified AE model for modes excited at the Alfven wave continuum maximum with geodesic corrections taken into account, the temporal evolution of the value of the safety factor q(0) at the magnetic axis is determined We describe a new scheme to determine the time variation of q(0) that works under conditions in which other standard diagnostics, such as the motional Stark effect do not give reliable results such as during a monster sawtooth [doi 10 1063/1 3494212]

Cloud condensation nuclei in pristine tropical rainforest air of Amazonia: size-resolved measurements and modeling of atmospheric aerosol composition and CCN activity

Atmospheric aerosol particles serving as cloud condensation nuclei (CCN) are key elements of the hydrological cycle and climate. We have measured and characterized CCN at water vapor supersaturations in the range of S=0.10-0.82% in pristine tropical rainforest air during the AMAZE-08 campaign in central Amazonia. The effective hygroscopicity parameters describing the influence of chemical composition on the CCN activity of aerosol particles varied in the range of kappa approximate to 0.1-0.4 (0.16+/-0.06 arithmetic mean and standard deviation). The overall median value of kappa approximate to 0.15 was by a factor of two lower than the values typically observed for continental aerosols in other regions of the world. Aitken mode particles were less hygroscopic than accumulation mode particles (kappa approximate to 0.1 at D approximate to 50 nm; kappa approximate to 0.2 at D approximate to 200 nm), which is in agreement with earlier hygroscopicity tandem differential mobility analyzer (H-TDMA) studies. The CCN measurement results are consistent with aerosol mass spectrometry (AMS) data, showing that the organic mass fraction (f(org)) was on average as high as similar to 90% in the Aitken mode (D <= 100 nm) and decreased with increasing particle diameter in the accumulation mode (similar to 80% at D approximate to 200 nm). The kappa values exhibited a negative linear correlation with f(org) (R(2)=0.81), and extrapolation yielded the following effective hygroscopicity parameters for organic and inorganic particle components: kappa(org)approximate to 0.1 which can be regarded as the effective hygroscopicity of biogenic secondary organic aerosol (SOA) and kappa(inorg)approximate to 0.6 which is characteristic for ammonium sulfate and related salts. Both the size dependence and the temporal variability of effective particle hygroscopicity could be parameterized as a function of AMS-based organic and inorganic mass fractions (kappa(p)=kappa(org) x f(org)+kappa(inorg) x f(inorg)). The CCN number concentrations predicted with kappa(p) were in fair agreement with the measurement results (similar to 20% average deviation). The median CCN number concentrations at S=0.1-0.82% ranged from N(CCN,0.10)approximate to 35 cm(-3) to N(CCN,0.82)approximate to 160 cm(-3), the median concentration of aerosol particles larger than 30 nm was N(CN,30)approximate to 200 cm(-3), and the corresponding integral CCN efficiencies were in the range of N(CCN,0.10/NCN,30)approximate to 0.1 to N(CCN,0.82/NCN,30)approximate to 0.8. Although the number concentrations and hygroscopicity parameters were much lower in pristine rainforest air, the integral CCN efficiencies observed were similar to those in highly polluted megacity air. Moreover, model calculations of N(CCN,S) assuming an approximate global average value of kappa approximate to 0.3 for continental aerosols led to systematic overpredictions, but the average deviations exceeded similar to 50% only at low water vapor supersaturation (0.1%) and low particle number concentrations (<= 100 cm(-3)). Model calculations assuming aconstant aerosol size distribution led to higher average deviations at all investigated levels of supersaturation: similar to 60% for the campaign average distribution and similar to 1600% for a generic remote continental size distribution. These findings confirm earlier studies suggesting that aerosol particle number and size are the major predictors for the variability of the CCN concentration in continental boundary layer air, followed by particle composition and hygroscopicity as relatively minor modulators. Depending on the required and applicable level of detail, the information and parameterizations presented in this paper should enable efficient description of the CCN properties of pristine tropical rainforest aerosols of Amazonia in detailed process models as well as in large-scale atmospheric and climate models.

Hidden vorticity in binary Bose-Einstein condensates

We consider a binary Bose-Einstein condensate (BEC) described by a system of two-dimensional (2D) Gross-Pitaevskii equations with the harmonic-oscillator trapping potential. The intraspecies interactions are attractive, while the interaction between the species may have either sign. The same model applies to the copropagation of bimodal beams in photonic-crystal fibers. We consider a family of trapped hidden-vorticity (HV) modes in the form of bound states of two components with opposite vorticities S(1,2) = +/- 1, the total angular momentum being zero. A challenging problem is the stability of the HV modes. By means of a linear-stability analysis and direct simulations, stability domains are identified in a relevant parameter plane. In direct simulations, stable HV modes feature robustness against large perturbations, while unstable ones split into fragments whose number is identical to the azimuthal index of the fastest growing perturbation eigenmode. Conditions allowing for the creation of the HV modes in the experiment are discussed too. For comparison, a similar but simpler problem is studied in an analytical form, viz., the modulational instability of an HV state in a one-dimensional (1D) system with periodic boundary conditions (this system models a counterflow in a binary BEC mixture loaded into a toroidal trap or a bimodal optical beam coupled into a cylindrical shell). We demonstrate that the stabilization of the 1D HV modes is impossible, which stresses the significance of the stabilization of the HV modes in the 2D setting.

Quantum statistical correlations in thermal field theories: Boundary effective theory

We show that the one-loop effective action at finite temperature for a scalar field with quartic interaction has the same renormalized expression as at zero temperature if written in terms of a certain classical field phi(c), and if we trade free propagators at zero temperature for their finite-temperature counterparts. The result follows if we write the partition function as an integral over field eigenstates (boundary fields) of the density matrix element in the functional Schrodinger field representation, and perform a semiclassical expansion in two steps: first, we integrate around the saddle point for fixed boundary fields, which is the classical field phi(c), a functional of the boundary fields; then, we perform a saddle-point integration over the boundary fields, whose correlations characterize the thermal properties of the system. This procedure provides a dimensionally reduced effective theory for the thermal system. We calculate the two-point correlation as an example.

Pressure of massless hot scalar theory in the boundary effective theory framework

We use the boundary effective theory approach to thermal field theory in order to calculate the pressure of a system of massless scalar fields with quartic interaction. The method naturally separates the infrared physics, and is essentially nonperturbative. To lowest order, the main ingredient is the solution of the free Euler-Lagrange equation with nontrivial (time) boundary conditions. We derive a resummed pressure, which is in good agreement with recent calculations found in the literature, following a very direct and compact procedure.

Magnetoconfined levels in a parabolic quantum dot: An analytical solution of a three-dimensional Fock-Darwin problem in a tilted magnetic field

The energy spectrum of an electron confined in a quantum dot (QD) with a three-dimensional anisotropic parabolic potential in a tilted magnetic field was found analytically. The theory describes exactly the mixing of in-plane and out-of-plane motions of an electron caused by a tilted magnetic field, which could be seen, for example, in the level anticrossing. For charged QDs in a tilted magnetic field we predict three strong resonant lines in the far-infrared-absorption spectra.

f t value of the 0(+)-> 0(+) beta(+) decay of (32)Ar: A measurement of isospin symmetry breaking in a superallowed decay

We determined the absolute branch of the T=2 superallowed decay of (32)Ar by detecting the beta(+)-delayed protons and gamma decays of the daughter state. We obtain b(SA)(beta)=(22.71 +/- 0.16)%, which represents the first determination of a proton branch to better than 1%. Using this branch along with the previously determined (32)Ar half-life and energy release, we determined ft=(1552 +/- 12) s for the superallowed decay. This ft value, together with the corrected Ft value extracted from previously known T=1 superallowed decays, yields a measurement of the isospin symmetry breaking correction in (32)Ar decay delta(exp)(C)=(2.1 +/- 0.8)%. This can be compared to a theoretical calculation delta(C)=(2.0 +/- 0.4)%. As by-products of this work, we determined the gamma and proton branches for the decay of the lowest T=2 state of (32)Cl, made a precise determination of the total proton branch and relative intensities of proton groups that leave (31)S in its first excited state and deduced an improved value for the (32)Cl mass.

gamma Cassiopeiae: an X-ray Be star with personality

An exciting unsolved problem in the study of high energy processes of early type stars concerns the physical mechanism for producing X-rays near the Be star gamma Cassiopeiae. By now we know that this source and several ""gamma Cas analogs"" exhibit an unusual hard thermal X-ray spectrum, compared both to normal massive stars and the non-thermal emission of known Be/X-ray binaries. Also, its light curve is variable on almost all conceivable timescales. In this study we reanalyze a high dispersion spectrum obtained by Chandra in 2001 and combine it with the analysis of a new (2004) spectrum and light curve obtained by XMM-Newton. We find that both spectra can be fit well with 3-4 optically thin, thermal components consisting of a hot component having a temperature kT(Q) similar to 12-14 keV, perhaps one with a value of similar to 2.4 keV, and two with well defined values near 0.6 keV and 0.11 keV. We argue that these components arise in discrete (almost monothermal) plasmas. Moreover, they cannot be produced within an integral gas structure or by the cooling of a dominant hot process. Consistent with earlier findings, we also find that the Fe abundance arising from K-shell ions is significantly subsolar and less than the Fe abundance from L-shell ions. We also find novel properties not present in the earlier Chandra spectrum, including a dramatic decrease in the local photoelectric absorption of soft X-rays, a decrease in the strength of the Fe and possibly of the Si K fluorescence features, underpredicted lines in two ions each of Ne and N (suggesting abundances that are similar to 1.5-3x and similar to 4x solar, respectively), and broadening of the strong NeXLy alpha and OVIII Ly alpha lines. In addition, we note certain traits in the gamma Cas spectrum that are different from those of the fairly well studied analog HD110432 - in this sense the stars have different ""personalities."" In particular, for gamma Cas the hot X-ray component remains nearly constant in temperature, and the photoelectric absorption of the X-ray plasmas can change dramatically. As found by previous investigators of gamma Cas, changes in flux, whether occurring slowly or in rapidly evolving flares, are only seldomly accompanied by variations in hardness. Moreover, the light curve can show a ""periodicity"" that is due to the presence of flux minima that recur semiregularly over a few hours, and which can appear again at different epochs.

TESTING STATISTICAL HYPOTHESIS ON RANDOM TREES AND APPLICATIONS TO THE PROTEIN CLASSIFICATION PROBLEM

Efficient automatic protein classification is of central importance in genomic annotation. As an independent way to check the reliability of the classification, we propose a statistical approach to test if two sets of protein domain sequences coming from two families of the Pfam database are significantly different. We model protein sequences as realizations of Variable Length Markov Chains (VLMC) and we use the context trees as a signature of each protein family. Our approach is based on a Kolmogorov-Smirnov-type goodness-of-fit test proposed by Balding et at. [Limit theorems for sequences of random trees (2008), DOI: 10.1007/s11749-008-0092-z]. The test statistic is a supremum over the space of trees of a function of the two samples; its computation grows, in principle, exponentially fast with the maximal number of nodes of the potential trees. We show how to transform this problem into a max-flow over a related graph which can be solved using a Ford-Fulkerson algorithm in polynomial time on that number. We apply the test to 10 randomly chosen protein domain families from the seed of Pfam-A database (high quality, manually curated families). The test shows that the distributions of context trees coming from different families are significantly different. We emphasize that this is a novel mathematical approach to validate the automatic clustering of sequences in any context. We also study the performance of the test via simulations on Galton-Watson related processes.

Eddy Formation in the Tropical Atlantic Induced by Abrupt Changes in the Meridional Overturning Circulation

The variability of the meridional overturning circulation (MOC) in the upper tropical Atlantic basin is investigated using a reduced-gravity model in a simplified domain. Four sets of idealized numerical experiments are performed: (i) switch-on of the MOC until a fixed value when a constant northward flow is applied along the western boundary; (ii) MOC with a variable flow; (iii) MOC in a quasi-steady flow; and (iv) shutdown of the MOC in the Northern Hemisphere. Results from experiments (i) show that eddies are generated at the equatorial region by shear instability and detached northward; eddies are responsible for an enhancement of the mean flow and the variability of the MOC. Results from experiments (ii) show a transitional behavior of the MOC related to the eddy generation in interannual-decadal time scales as the Reynolds number varies due to the variations in the MOC. In experiments (iii), a critical Reynolds number Re(c) around 30 is found, above which eddies are generated. Experiments (iv) demonstrate that even after the collapse of MOC in the Northern Hemisphere, eddies can still be generated and carry energy across the equator into the Northern Hemisphere; these eddies act to attenuate the impact of the MOC shutdown on short time scales. The results described here may be particularly pertinent to ocean general circulation models in which the Reynolds number lies close to the bifurcation point separating the laminar and turbulent regimes.

Homogenization in a thin domain with an oscillatory boundary

In this paper we analyze the behavior of the Laplace operator with Neumann boundary conditions in a thin domain of the type R(epsilon) = {(x(1), x(2)) is an element of R(2) vertical bar x(1) is an element of (0, 1), 0 < x(2) < epsilon G(x(1), x(1)/epsilon)} where the function G(x, y) is periodic in y of period L. Observe that the upper boundary of the thin domain presents a highly oscillatory behavior and, moreover, the height of the thin domain, the amplitude and period of the oscillations are all of the same order, given by the small parameter epsilon. (C) 2011 Elsevier Masson SAS. All rights reserved.

Analytical calculation of the transition to complete phase synchronization in coupled oscillators

Here we present a system of coupled phase oscillators with nearest neighbors coupling, which we study for different boundary conditions. We concentrate at the transition to the total synchronization. We are able to develop exact solutions for the value of the coupling parameter when the system becomes completely synchronized, for the case of periodic boundary conditions as well as for a chain with fixed ends. We compare the results with those calculated numerically.

An alternative to the Pythagorean rule? Reevaluating Problem 1 of cuneiform tablet BM 34 568

The first problem of the Seleucid mathematical cuneiform tablet BM 34 568 calculates the diagonal of a rectangle from its sides without resorting to the Pythagorean rule. For this reason, it has been a source of discussion among specialists ever since its first publication. but so far no consensus in relation to its mathematical meaning has been attained. This paper presents two new interpretations of the scribe`s procedure. based on the assumption that he was able to reduce the problem to a standard Mesopotamian question about reciprocal numbers. These new interpretations are then linked to interpretations of the Old Babylonian tablet Plimpton 322 and to the presence of Pythagorean triples in the contexts of Old Babylonian and Hellenistic mathematics. (C) 2007 Elsevier Inc. All rights reserved.

Stability region bifurcations of nonlinear autonomous dynamical systems: Type-zero saddle-node bifurcations

The behavior of stability regions of nonlinear autonomous dynamical systems subjected to parameter variation is studied in this paper. In particular, the behavior of stability regions and stability boundaries when the system undergoes a type-zero sadle-node bifurcation on the stability boundary is investigated in this paper. It is shown that the stability regions suffer drastic changes with parameter variation if type-zero saddle-node bifurcations occur on the stability boundary. A complete characterization of these changes in the neighborhood of a type-zero saddle-node bifurcation value is presented in this paper. Copyright (C) 2010 John Wiley & Sons, Ltd.