Stable loosely-coupled-type algorithm for fluid–structure interaction in blood flow

We introduce a novel loosely coupled-type algorithm for fluid–structure interaction between blood flow and thin vascular walls. This algorithm successfully deals with the difficulties associated with the “added mass effect”, which is known to be the cause of numerical instabilities in fluid–structure interaction problems involving fluid and structure of comparable densities. Our algorithm is based on a time-discretization via operator splitting which is applied, in a novel way, to separate the fluid sub-problem from the structure elastodynamics sub-problem. In contrast with traditional loosely-coupled schemes, no iterations are necessary between the fluid and structure sub-problems; this is due to the fact that our novel splitting strategy uses the “added mass effect” to stabilize rather than to destabilize the numerical algorithm. This stabilizing effect is obtained by employing the kinematic lateral boundary condition to establish a tight link between the velocities of the fluid and of the structure in each sub-problem. The stability of the scheme is discussed on a simplified benchmark problem and we use energy arguments to show that the proposed scheme is unconditionally stable. Due to the crucial role played by the kinematic lateral boundary condition, the proposed algorithm is named the “kinematically coupled scheme”.     

Scattering in thick multifractal clouds, Part I: Overview and single scattering

Over the last twenty years, many studies have been made of radiative transfer in scaling cloud fields, the vast majority of which have been limited to numerical studies in clouds with relatively small optical thickness. An exception to this was the development of a formalism for treating single scattering in optically thick but conservative multifractal clouds without significant holes. In part I of this paper we show how these results can be extended to general “universal” multifractal clouds dominated by low density “Lévy holes”. In part II, we demonstrate how the analytic single scattering results can be generalized to multiple scattering including the case of very thick clouds as well as to realistic nonconservative clouds.     

Two mechanisms of spiral-pair-source creation in excitable media

Two new modes of generating spiral pairs in an excitable medium have been found. They depend on a geometrical structure (GS) inside the medium. This may be formed e.g. as a result of scars or fibrosis in the heart tissue, or artificially built in a chemical reaction substrate. Both sources involve a GS composed of a circular “convergent lens” bounded by two opaque “walls”. One mode can be induced by a single wave and behaves as a “flip-flop” type of a limit cycle. The other mode is generated by a train of plane waves impinging on the GS, and is created at the focus of the converging wave-fragments.     

Kerr black holes are not fragile

Certain AdS black holes are “fragile”, in the sense that, if they are deformed excessively, they become unstable to a fundamental non-perturbative stringy effect analogous to Schwinger pair-production [of branes]. Near-extremal topologically spherical AdS-Kerr black holes, which are natural candidates for string-theoretic models of the very rapidly rotating black holes that have actually been observed to exist, do represent a very drastic deformation of the AdS-Schwarzschild geometry. One therefore has strong reason to fear that these objects might be “fragile”, which in turn could mean that asymptotically flat rapidly rotating black holes might be fragile in string theory. Here we show that this does not happen: despite the severe deformation implied by near-extremal angular momenta, brane pair-production around topologically spherical AdS-Kerr-Newman black holes is always suppressed.     

Emergent gravity in two dimensions

We explore models with emergent gravity and metric by means of numerical simulations. A particular type of two-dimensional non-linear sigma-model is regularized and discretized on a quadratic lattice. It is characterized by lattice diffeomorphism invariance which ensures in the continuum limit the symmetry of general coordinate transformations. We observe a collective order parameter with properties of a metric, showing Minkowski or Euclidean signature. The correlation functions of the metric reveal an interesting long-distance behavior with power-like decay. This universal critical behavior occurs without tuning of parameters and thus constitutes an example of “self-tuned criticality” for this type of sigma-models. We also find a non-vanishing expectation value of a “zweibein” related to the “internal” degrees of freedom of the scalar field, again with long-range correlations. The metric is well described as a composite of the zweibein. A scalar condensate breaks Euclidean rotation symmetry.     

Frequency modulation indicator, Arnold’s web and diffusion in the Stark-Quadratic-Zeeman problem

We notice that the fundamental frequencies of a slightly perturbed integrable Hamiltonian system are not time-constant inside a resonance but frequency modulated, as is evident from pendulum models and wavelet analysis. Exploiting an intrinsic imprecision inherent to the numerical frequency analysis algorithm itself, hence transforming a drawback into an opportunity, we define the Frequency Modulation Indicator, a very sensitive tool in detecting where fundamental frequencies are modulated, localizing so the resonances without having to resort, as in other methods, to the integration of variational equations. For the Kepler problem, the space of the orbits with a fixed energy has the topology of the product of two 2-spheres. The perturbation Hamiltonian, averaged over the mean anomaly, has surely a maximum and a minimum, to which correspond two periodic orbits in physical space. Studying the neighbourhood of these two elliptic stable points, we are able to define adapted action-angle variables, for example, the usual but “SO(4)-rotated” Delaunay variables. The procedure, implemented in the program KEPLER, is performed transparently for the user, providing a general scheme suited for generic perturbation. The method is then applied to the Stark-Quadratic-Zeeman problem, displaying very clearly the Arnold web of the resonances. Sectioning transversely one of the resonance strips so highlighted and performing a numerical frequency analysis, one is able to locate with great precision the thin stochastic layer surrounding a separatrix. Another very long (10⁸ revolutions) frequency analysis on an orbit starting here reveals, as expected, a well defined pattern, which ensures that the integration errors do not eject the point out of the layer, and moreover a very slow drift in the frequency values, clearly due to Arnold diffusion.     

Understanding the time dependence of atomic level populations in evolving plasmas

The time dependence of atomic level populations in evolving plasmas is studied using an eigenfunction expansion of the non-LTE rate equations. The work aims to develop understanding without the need for, and as an aid to, numerical solutions. The discussion is mostly limited to linear systems, especially those for optically thin plasmas, but the implicitly non-linear case of non-LTE radiative transfer is briefly discussed. Eigenvalue spectra for typical atomic systems are examined using results compiled by Hearon. Diagonal dominance and sign symmetry of rate matrices show that just one eigenvalue is zero (corresponding to the equilibrium state), that the remaining eigenvalues have negative real parts, and that oscillations, if any, are necessarily damped. Gershgorin's theorems are used to show that many eigenvalues are determined by the radiative lifetimes of certain levels, because of diagonal dominance. With other properties, this demonstrates the existence of both “slow” and “fast” time-scales, where the “slow” evolution is controlled by properties of meta-stable levels. It is shown that, when collisions are present, Rydberg states contribute only “fast” eigenvalues. This justifies use of the quasi-static approximation, in which atoms containing just meta-stable levels can suffice to determine the atomic evolution on time-scales long compared with typical radiative lifetimes. Analytic solutions for two- and three-level atoms are used to examine the basis of earlier intuitive ideas, such as the “ionizing plasma” approximation. The power and limitations of Gershgorin's theorems are examined through examples taken from the solar atmosphere. The methods should help in the planning and interpretation of both experimental and numerical experiments in which atomic evolution is important. While the examples are astrophysical, the methods and results are applicable to plasmas in general.     

Study of mutual information multimodality medical image registration based on modified simplex optimization method

Because of a different imaging mechanism and highly complexity of body tissues and structures. Different modality medical images provide non-overlay complementary information. This has very important significance for multimodal medical image registration. Image registration is the first and key part of problem to be solved in the integrations. When the spatial position of two medical images is same, the registration could be achieved. For two CT and PET images, the principal axis method is adopted to achieve the rough registration. The modified simplex algorithm is employed to implement global search using the mutual information as similarity measure. The initial registration parameters are achieved through principal axis Based on the results of test, improved simplex method can adjust reflecting distance. Stepped-up optimization algorithm on the new experimental points through the methods of “reflection”, “enlargement”, “shrinkage” or “global systolic”. A mutual information registration based on modified simplex optimization method is presented in this paper to improve the speed of medical image registration.Results indicate that the proposed registration method prevents the optimizing process from falling into local extremum and improves the convergence speed while keeping the precision. The accurate registration of multimodal image with different resolutions is achieved.     

An optical study of the correlation between growth kinetics and microstructure of μc-Si grown by SiH4-H2 PECVD

Fully microcrystalline silicon, μc-Si, thin films have been deposited on corning glass by plasma enhanced chemical vapor deposition (PECVD) using SiH₄-H₂. The effects of the surface treatment and of the deposition temperature on microstructure of μc-Si films are investigated by “in situ” laser reflectance interferometry (LRI), “ex situ” spectroscopic ellipsometry (SE) and Raman spectroscopy. LRI indicated the existence of a “crystalline seeding time”, which is indicative of the crystallite nucleation, and depends on substrate treatments. Longer “crystalline seeding time” results in a lower density of crystalline nuclei, which grow laterally, yielding to complete suppression of the amorphous incubation layer and to growth of very dense, fully crystalline layer at a growth temperature as low as 120 °C.     

Exchange rate regimes, saving glut and the Feldstein-Horioka puzzle: The East Asian experience

This paper investigates whether the recent experience of the emerging East Asian countries with current account surpluses is consistent with the “saving glut” hypothesis and the Feldstein and Horioka puzzle. The evidence suggests that the saving retention coefficients declined substantially in most of the countries after an endogenous break date coinciding with a major exchange rate regime change with the 1997-1998 crisis. Exchange rate flexibility appears to be enhancing financial integration. The results are consistent with an “investment slump” explanation rather than the “saving glut” postulation.     

Image segmentation using adaptive loopy belief propagation

Loopy belief propagation (LBP) algorithm over pairwise-connected Markov random fields (MRFs) has become widely used for low-level vision problems. However, Pairwise MRF is often insufficient to capture the statistics of natural images well, and LBP is still extremely slow for application on an MRF with large discrete label space. To solve these problems, the present study proposes a new segmentation algorithm based on adaptive LBP. The proposed algorithm utilizes local region information to construct a local region model, as well as a local interaction region MRF model for image segmentation. The adaptive LBP algorithm maximizes the global probability of the proposed MRF model, which employs two very important strategies, namely, “message self-convergence” and “adaptive label pruning”. Message self-convergence can improve the reliability of a pixel in choosing a label in local region, and label pruning can dismiss impossible labels for every pixel. Thus, the most reliable information messages transfer through the LBP algorithm. The experimental results show that the proposed algorithm not only obtains more accurate segmentation results but also greater speed.     

Sequential estimation of velocity fields using episodic proper orthogonal decomposition

In this paper, the problem of approximating velocity fields at future and past times based on information available at the current time is addressed. A novel method called “episodic POD” is described and developed that enables us to achieve our objective. Application of episodic POD to an ensemble of flow data results in a set of spatio-temporal eigenfunctions and a set of coefficients associated with the eigenfunctions. From these eigenfunctions, we develop two models called the “forward model” and “inverse model” that enable us to approximate the velocity fields at future and past times, based on information provided at the current time. A second set of models, the forward and inverse sequential models are also developed that enable the dynamic update of approximated velocity fields when new information is made available, making these models more adept at on-line estimation. The various properties associated with these models are described in detail, and four different examples are used to validate the models and show the different properties associated with episodic POD. It is shown through numerical validation, that the episodic POD model has a form that is dynamically consistent with the original system. It is also shown that episodic POD outperforms linear Kalman filters in the presence of noise.     

Elastic modes frequencies in the mirrors with “ears” for parametric oscillatory instability in advanced LIGO interferometer

We discuss the importance of accurate numerical calculation of elastic modes in the mirrors with suspension “ears” in advanced LIGO interferometer to enable precise predictions of parametric oscillatory instability problem. We show that “ears” of test masses produce additional shift of elastic modes frequencies which is larger than relaxation rate of optical modes. These shifts may increase the possibility of parametric oscillatory instability.     

A heuristic approach to automated molecular line assignment

A first investigation on the feasibility of automated molecular line assignment is presented. Dense rovibrational molecular spectra are normally assigned by strongly interactive computer methods, ranging from commercial spreadsheets to dedicated programs, like Loomis-Wood or Ritz. While a general-purpose, fully automated assignment procedure seems to be out of reach for the near future, we show that a thorough investigation of the problem can lead to new, more efficient and less interactive methods, at least in reasonably favorable conditions. Interesting suggestions are provided by some modern “heuristic” problem-solving algorithms, which mimic natural processes. As a first step, we have developed a “transgenic-evolutionary” algorithm, which has successfully assigned artificial spectra of up to almost 3500 lines. We discuss also its performance on an experimental, but “filtered,” methanol spectrum. Possible future improvements and developments of this method, as well as its limits, are discussed.     

Hydrodynamics and phases of flocks

We review the past decade’s theoretical and experimental studies of flocking: the collective, coherent motion of large numbers of self-propelled “particles” (usually, but not always, living organisms). Like equilibrium condensed matter systems, flocks exhibit distinct “phases” which can be classified by their symmetries. Indeed, the phases that have been theoretically studied to date each have exactly the same symmetry as some equilibrium phase (e.g., ferromagnets, liquid crystals). This analogy with equilibrium phases of matter continues in that all flocks in the same phase, regardless of their constituents, have the same “hydrodynamic”—that is, long-length scale and long-time behavior, just as, e.g., all equilibrium fluids are described by the Navier-Stokes equations. Flocks are nonetheless very different from equilibrium systems, due to the intrinsically nonequilibrium self-propulsion of the constituent “organisms.” This difference between flocks and equilibrium systems is most dramatically manifested in the ability of the simplest phase of a flock, in which all the organisms are, on average moving in the same direction (we call this a “ferromagnetic” flock; we also use the terms “vector-ordered” and “polar-ordered” for this situation) to exist even in two dimensions (i.e., creatures moving on a plane), in defiance of the well-known Mermin-Wagner theorem of equilibrium statistical mechanics, which states that a continuous symmetry (in this case, rotation invariance, or the ability of the flock to fly in any direction) can not be spontaneously broken in a two-dimensional system with only short-ranged interactions. The “nematic” phase of flocks, in which all the creatures move preferentially, or are simply oriented preferentially, along the same axis, but with equal probability of moving in either direction, also differs dramatically from its equilibrium counterpart (in this case, nematic liquid crystals). Specifically, it shows enormous number fluctuations, which actually grow with the number of organisms faster than the “law of large numbers” obeyed by virtually all other known systems. As for equilibrium systems, the hydrodynamic behavior of any phase of flocks is radically modified by additional conservation laws. One such law is conservation of momentum of the background fluid through which many flocks move, which gives rise to the “hydrodynamic backflow” induced by the motion of a large flock through a fluid. We review the theoretical work on the effect of such background hydrodynamics on three phases of flocks—the ferromagnetic and nematic phases described above, and the disordered phase in which there is no order in the motion of the organisms. The most surprising prediction in this case is that “ferromagnetic” motion is always unstable for low Reynolds-number suspensions. Experiments appear to have seen this instability, but a quantitative comparison is awaited. We conclude by suggesting further theoretical and experimental work to be done.     

Spooky action-at-a-distance in general probabilistic theories

We call a probabilistic theory “complete” if it cannot be further refined by no-signaling hidden-variable models, and name a theory “spooky” if every equivalent hidden-variable model violates Shimony?s Outcome Independence. We prove that a complete theory is spooky if and only if it admits a pure steering state in the sense of Schrödinger. Finally we show that steering of complementary states leads to a Schrödinger?s-cat-like paradox.     

A leap-frog discontinuous Galerkin time-domain method of analyzing electromagnetic scattering problems

Several major challenges need to be faced for efficient transient multiscale electromagnetic simulations, such as flexible and robust geometric modeling schemes, efficient and stable time-stepping algorithms, etc. Fortunately, because of the versatile choices of spatial discretization and temporal integration, a discontinuous Galerkin time-domain(DGTD) method can be a very promising method of solving transient multiscale electromagnetic problems. In this paper, we present the application of a leap-frog DGTD method to the analyzing of the multiscale electromagnetic scattering problems. The uniaxial perfect matching layer(UPML) truncation of the computational domain is discussed and formulated in the leap-frog DGTD context. Numerical validations are performed in the challenging test cases demonstrating the accuracy and effectiveness of the method in solving transient multiscale electromagnetic problems compared with those of other numerical methods.     

A numerical scheme for the solution of axisymmetric radiative transfer problems in irregular domains filled by media with opaque and transparent diffuse and specular boundaries

This paper presents a new numerical scheme of the discrete ordinates method for the solution of axisymmetric radiative transfer problems in irregular domains filled by media with opaque and transparent diffuse and specular (Fresnel) boundaries and interfaces. New test problems of radiative transfer, which describe radiative transfer in domains with Fresnel interfaces, are proposed in this paper. These problems admit analytic solutions and can be used as benchmark ones. The proposed scheme is applied to the solution of the problems. Numerical results show that the presence of Fresnel interfaces leads to an appreciably larger error in numerical solution. This is connected with the “discontinuity” of the Fresnel reflectivity, which, through numerical diffusion, leads to the distortion of numerical solution. Modification of the scheme allows to reduce the numerical error.     

A novel fast numerical algorithm for cascaded Raman fiber laser using the analytic approximate solution

In this paper, we proposed a novel numerical algorithm for cascaded Raman fiber laser (CRFL) using the approximate analytic results as the initial values for shooting method, which effectively reduces the calculating time from several hours to a few minutes. With the algorithm, we obtained a numerical solution of the bilateral-pumping Ge-doped fifth-order CRFL which can avoid the so-called “end face damage” phenomenon efficiently. At the same time, we also simulated the unilateral-pumping one as a comparison, which showed both the characteristics of them are similar with each other.     

Statistical physics of media processes: Mediaphysics

The processes of mass communications in complicated social or sociobiological systems such as marketing, economics, politics, animal populations, etc. as a subject for the special scientific subbranch—“mediaphysics”—are considered in its relation with sociophysics. A new statistical physics approach to analyze these phenomena is proposed. A keystone of the approach is an analysis of population distribution between two or many alternatives: brands, political affiliations, or opinions. Relative distances between a state of a “person's mind” and the alternatives are measures of propensity to buy (to affiliate, or to have a certain opinion). The distribution of population by those relative distances is time dependent and affected by external (economic, social, marketing, natural) and internal (influential propagation of opinions, “word of mouth”, etc.) factors, considered as fields. Specifically, the interaction and opinion-influence field can be generalized to incorporate important elements of Ising-spin-based sociophysical models and kinetic-equation ones. The distributions were described by a Schrödinger-type equation in terms of Green's functions. The developed approach has been applied to a real mass-media efficiency problem for a large company and generally demonstrated very good results despite low initial correlations of factors and the target variable.