A non-conforming monolithic finite element method for problems of coupled mechanics

In this study, a Lagrange multiplier technique is developed to solve problems of coupled mechanics and is applied to the case of a Newtonian fluid coupled to a quasi-static hyperelastic solid. Based on theoretical developments in [57], an additional Lagrange multiplier is used to weakly impose displacement/velocity continuity as well as equal, but opposite, force. Through this approach, both mesh conformity and kinematic variable interpolation may be selected independently within each mechanical body, allowing for the selection of grid size and interpolation most appropriate for the underlying physics. In addition, the transfer of mechanical energy in the coupled system is proven to be conserved. The fidelity of the technique for coupled fluid–solid mechanics is demonstrated through a series of numerical experiments which examine the construction of the Lagrange multiplier space, stability of the scheme, and show optimal convergence rates. The benefits of non-conformity in multi-physics problems is also highlighted. Finally, the method is applied to a simplified elliptical model of the cardiac left ventricle.     

The backward phase flow and FBI-transform-based Eulerian Gaussian beams for the Schrödinger equation

We propose the backward phase flow method to implement the Fourier–Bros–Iagolnitzer (FBI)-transform-based Eulerian Gaussian beam method for solving the Schrödinger equation in the semi-classical regime. The idea of Eulerian Gaussian beams has been first proposed in [12]. In this paper we aim at two crucial computational issues of the Eulerian Gaussian beam method: how to carry out long-time beam propagation and how to compute beam ingredients rapidly in phase space. By virtue of the FBI transform, we address the first issue by introducing the reinitialization strategy into the Eulerian Gaussian beam framework. Essentially we reinitialize beam propagation by applying the FBI transform to wavefields at intermediate time steps when the beams become too wide. To address the second issue, inspired by the original phase flow method, we propose the backward phase flow method which allows us to compute beam ingredients rapidly. Numerical examples demonstrate the efficiency and accuracy of the proposed algorithms.     

欧拉游历间K1—变换下的算法复杂性

本文主要讨论了欧拉图方面的国际权威ＨｅｒｂｅｒｔＦｌｅｉｃｈｎｅｒ教授所著书［１］中的一个问题；对欧拉图Ｇ的任意两个欧拉游历要经过多少Ｋ－变换或Ｋ＾＊－变换才能从一个游历得到另一个游历？我们得到的结论是：对欧拉图Ｇ中的任两个欧拉游历Ｔ、Ｔ＇最多经过‖Ｅ（Ｇ）‖－‖Ｖ（Ｇ）‖－变换可以使Ｔ变换成Ｔ＇。且此结果不能再改进。进一步我们分别对Ｋ－变换和Ｋ＆－变换的算法复杂性进行了讨论。     

Three‐dimensional computation for flow‐induced vibrations in in‐line and cross‐flow directions of a circular cylinder

We present numerical results for in‐line and cross‐flow vibrations of a circular cylinder, which is immersed in a uniform flow and is elastically supported by damper‐spring systems to compute vibrations of a rigid cylinder. In the case of a circular cylinder with a low Scruton number, it is well‐known that two types of self‐excited vibrations appear in the in‐line direction in the range of low reduced velocities. On the other hand, a cross‐flow vibration of the circular cylinder can be excited in the range of high reduced velocities. Therefore, we compute the flow‐induced vibrations of the circular cylinder in the wide range of the reduced velocities at low and high Scruton numbers and discuss about excitation mechanisms in the in‐line and cross‐flow directions. Copyright © 2011 John Wiley & Sons, Ltd.     

Eulerian numbers: A spline perspective

In this paper, the spline interpretations of Eulerian numbers and refined Eulerian numbers are presented. Many classical results about Eulerian numbers can follow from the properties of B-splines directly, and some new results about the refined Eulerian numbers and descent polynomials are also derived. Specifically, the explicit and recurrence formulas for the refined Eulerian numbers and descent polynomials are obtained. This paper also provides a new approach to study Eulerian numbers.     

The number of rooted Eulerian planar maps

In this paper we provide a solution of the functional equation unsolved in the paper, by the second author, "On functional equations arising from map enumerations" that appeared in Discrete Math, 123: 93-109 (1993). It is also the number of combinatorial distinct rooted general eulerian planar maps with the valency of root-vertex, the number of non-root vertices and non-root faces of the maps as three parameters. In particular, a result in the paper, by the same author, "On the number of eulerian planar map...     

Joint Distribution of Rises and Falls

The marginal distributions of the number of rises and the number of falls have been used successfully in various areas of statistics, especially in non-parametric statistical inference. Carlitz (1972, Duke Math. J. 39, 268–269) showed that the generating function of the joint distribution for the numbers of rises and falls satisfies certain complex combinatorial equations, and pointed out that he had been unable to derive the explicit formula for the joint distribution from these equations. After more than two decades, this latter problem remains unsolved. In this article, the joint distribution is obtained via the probabilistic method of finite Markov chain imbedding for random permutations. A numerical example is provided to illustrate the theoretical results and the corresponding computational procedures.     

A characteristic nonoverlapping domain decomposition method for multidimensional convection‐diffusion equations

We develop a quasi‐two‐level, coarse‐mesh‐free characteristic nonoverlapping domain decomposition method for unsteady‐state convection‐diffusion partial differential equations in multidimensional spaces. The development of the domain decomposition method is carried out by utilizing an additive Schwarz domain decomposition preconditioner, by using an Eulerian‐Lagrangian method for convection‐diffusion equations and by delicately choosing appropriate interface conditions that fully respect and utilize the hyperbolic nature of the governing equations. Numerical experiments are presented to illustrate the method. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Assessment of an Eulerian CFD model for prediction of dilute droplet dispersion in a turbulent jet

The turbulent dispersion of non-evaporating droplets in an axisymmetric round jet issuing from a nozzle is investigated both experimentally and theoretically. The experimental data set has a well-defined inlet boundary with low turbulence intensity at the nozzle exit, so that droplet dispersion is not affected by the transport of nozzle-generated fluctuating motion into the jet, and is influenced solely by turbulence in the gas phase produced in the shear layer of the jet. This data set is thus ideal for testing algebraic models of droplet fluctuating motion that assume local equilibrium with the turbulence in the gas phase. Moreover, the droplet flux measurements are sufficiently accurate that conservation of the total volume flow of the droplet phase has been demonstrated. A two-fluid turbulence modelling approach is adopted, which uses the k–ε turbulence model and a simple algebraic model that assumes local equilibrium to predict the fluid and droplet turbulent correlations, respectively. We have shown that the k–ε turbulence model lacks generality for predicting the spread of momentum in jets with and without a potential core. However, in general, the model predicts the radial dispersion of droplets in the considered turbulent jet with reasonable accuracy over a broad range of droplet sizes, once deficiencies in the k–ε turbulence model are taken into account.