A numerical method for simulating the dynamics of 3D axisymmetric vesicles suspended in viscous flows

We extend [Shravan K. Veerapaneni, Denis Gueyffier, Denis Zorin, George Biros, A boundary integral method for simulating the dynamics of inextensible vesicles suspended in a viscous fluid in 2D, Journal of Computational Physics 228(7) (2009) 2334–2353] to the case of three-dimensional axisymmetric vesicles of spherical or toroidal topology immersed in viscous flows. Although the main components of the algorithm are similar in spirit to the 2D case—spectral approximation in space, semi-implicit time-stepping scheme—the main differences are that the bending and viscous force require new analysis, the linearization for the semi-implicit schemes must be rederived, a fully implicit scheme must be used for the toroidal topology to eliminate a CFL-type restriction and a novel numerical scheme for the evaluation of the 3D Stokes single layer potential on an axisymmetric surface is necessary to speed up the calculations. By introducing these novel components, we obtain a time-scheme that experimentally is unconditionally stable, has low cost per time step, and is third-order accurate in time. We present numerical results to analyze the cost and convergence rates of the scheme. To verify the solver, we compare it to a constrained variational approach to compute equilibrium shapes that does not involve interactions with a viscous fluid. To illustrate the applicability of method, we consider a few vesicle-flow interaction problems: the sedimentation of a vesicle, interactions of one and three vesicles with a background Poiseuille flow.     

A numerical method for acoustic normal modes for shear flows

The normal modes and their propagation numbers for acoustic propagation in wave guides with flow are the eigenvectors and eigenvalues of a boundary value problem for a non-standard Sturm-Liouville problem. It is non-standard because it depends non-linearly on the eigenvalue parameter. (In the classical problem for ducts with no flow, the problem depends linearly on the eigenvalue parameter.) In this paper a method is presented for the fast numerical solution of this problem. It is a generalization of a method that was developed for the classical problem. A finite difference method is employed that combines well known numerical techniques and a generalization of the Sturm sequence method to solve the resulting algebraic eigenvalue problem. Then a modified Richardson extrapolation method is used that dramatically increases the accuracy of the computed eigenvalues. The method is then applied to two problems. They correspond to acoustic propagation in the ocean in the presence of a current, and to acoustic propagation in shear layers over flat plates.     

A Fictitious Domain,parallel numerical method for rigid particulate flows

In this paper, we develop a Fictitious Domain, parallel numerical method for the Direct Numerical Simulation of the flow of rigid particles in an incompressible viscous Newtonian fluid. A Simultaneous Directions Implicit algorithm is employed which gives the model a high level of parallelization. The projection of the fluid velocity onto rigid motion on the particles is based on a fast computational technique which relies on the conservation of linear and angular momenta. Numerical results are presented which confirm the ability of the proposed method to simulate the sedimentation of one and many particles; the parallel efficiency of the algorithm is also assessed.     

A multislope MUSCL method on unstructured meshes applied to compressible Euler equations for axisymmetric swirling flows

A finite volume method for the numerical solution of axisymmetric inviscid swirling flows is presented. The governing equations of the flow are the axisymmetric compressible Euler equations including swirl (or tangential) velocity. A first-order scheme is introduced where the convective fluxes at cell interfaces are evaluated by the Rusanov or the HLLC numerical flux while the geometric source terms are discretizated to provide a well-balanced scheme i.e. the steady-state solutions with null velocity are preserved. Extension to the second-order space approximation using a multislope MUSCL method is then derived. To test the numerical scheme, a stationary solution of the fluid flow following the radial direction has been established with a zero and nonzero tangential velocity. Numerical and exact solutions are compared for classical Riemann problems where we employ different limiters and effectiveness of the multislope MUSCL scheme is demonstrated for strongly shocked axially symmetric flows like in spherical bubble compression problem. Two other tests with axisymmetric geometries are performed: the supersonic flow in a tube with a cone and the axisymmetric blunt body with a free stream.     

Some considerations on the nonlinear stability of stationary planar Euler flows

We give sufficient conditions for the nonlinear stability of possibly nonsmooth stationary solutions of the two-dimensional Euler equation in symmetric bounded domains. We use, as Lyapunov functions, first integrals due to the symmetry of the problem. Moreover, we investigate the stability of smooth solutions under perturbations of the boundary. The last result is based on a generalization of the well known Arnold approach.Research partially supported by Italian CNR and Ministero della Pubblica Istruzione     

A new method based on the harmonic balance method for nonlinear oscillators

The harmonic balance (HB) method as an analytical approach is widely used for nonlinear oscillators, in which the initial conditions are generally simplified by setting velocity or displacement to be zero. Based on HB, we establish a new theory to address nonlinear conservative systems with arbitrary initial conditions, and deduce a set of over-determined algebraic equations. Since these deduced algebraic equations are not solved directly, a minimization problem is constructed instead and an iterative algorithm is employed to seek the minimization point. Taking Duffing and Duffing-harmonic equations as numerical examples, we find that these attained solutions are not only with high degree of accuracy, but also uniformly valid in the whole solution domain.     

A numerical method for one-dimensional nonlinear sine-Gordon equation using multiquadric quasi-interpolation

In this paper, we use a univariate multiquadric quasi-interpolation scheme to solve the one-dimensional nonlinear sine-Gordon equation that is related to many physical phenomena. We obtain a numerical scheme by using the derivative of the quasi-interpolation to approximate the spatial derivative and a difference scheme to approximate the temporal derivative. The advantage of the obtained scheme is that the algorithm is very simple so that it is very easy to implement. The results of numerical experiments are presented and compared with analytical solutions to confirm the good accuracy of the presented scheme.     

A nonlinear PSE method for two-fluid shear flows with complex interfacial topology

A nonlinear stability method is developed for laminar two-fluid shear flows which undergo changes in the interface topology. The method is based on the nonlinear parabolized stability equations (PSE) and incorporates a scalar-based interface capturing (IC) scheme in order to track complex deformations of the fluid interface. In doing so, the formulation retains the flexibility and physical insight of instability-wave based methods, while providing hydrodynamic modeling capabilities similar to direct numerical calculations: the new formulation, referred to as the IC-PSE, can capture the nonlinear physical mechanisms responsible for generating large-scale, two-fluid structures, without incurring heavy computational costs. This approach is valid for spatially developing, laminar two-fluid shear flows which are convectively unstable, and can naturally account for the growth of finite amplitude interfacial waves, along with changes to the interfacial topology. We demonstrate the accuracy of the IC-PSE against direct Navier–Stokes calculations for two-fluid mixing layers with density and viscosity stratification. The comparisons show that the IC-PSE can predict the dynamics of the instability waves and capture the formation of Kelvin–Helmholtz vortex rolls and large scale liquid structures, at an order of magnitude less computational cost than direct calculations. The role of surface tension in the IC-PSE formulation is shown to be valid for flows in which Re/We ? 1, and the method accurately predicts the formation and non-linear evolution of flow structures in this limit. This is demonstrated for spatially developing mixing layers which lead to vortex roll-up and ligaments, prior to droplet formation. The pinch-off process itself is a high surface tension phenomenon and in not considered herein. The method also accurately captures the effect of interfacial waves on the mean flow, and the topology changes during the non-linear evolution of the two-fluid structures.     

A high-order numerical method for the nonlinear Helmholtz equation in multidimensional layered media

We present a novel computational methodology for solving the scalar nonlinear Helmholtz equation (NLH) that governs the propagation of laser light in Kerr dielectrics.     

A two-dimensional analytical method for the transformation crystallography based on vector analysis

A simple 2-D analytical method has been developed for calculating the crystallographic features, including the orientation relationship (OR), habit plane (HP) and the direction and spacing of dislocations in the HP, under the condition that a pair of Burgers vectors from both phases are parallel to each other lying in the HP. In contrast to previous 2-D models for the OR and HP, which are on the plane normal to the parallel Burgers vectors, the present 2-D model is based on vector analysis in the HP. Closed-form expressions for the crystallographic features have been derived as functions of lattice parameter ratios. While the results are in agreement with those from a 2-D method based on the O-line model, the present method is simpler and improves our understanding of crystallographic features.     

非线性波动方程的新数值迭代方法

提出了一种新的求解非线性波动方程的数值迭代法,它是一种半解析的方法.与完全的数值计算方法扰法相比,它能够考虑各阶谐波的相互作用,且能够满足能量守恒定律.用它研究了非线性声波在液体中的传播性质,结果表明,在微扰法适用的声强范围内迭代法也适用,在微扰法不适用的一个较宽的声强范围内迭代法依然适用.     

基于矢量水听器的一种时间方位历程目标跟踪方法

文章简要介绍了时间方位历程显示方式的原理，针对时间方位历程中最大值自动跟踪方法存在的弊端，提出了应用矢量水听器解决多目标跟踪问题的一种新方法，并应用于海试数据处理中，取得了良好的效果。     

A numerical method for the study of nucleation of ordered phases

A numerical approach based on the string method is developed to study nucleation of ordered phases in first-order phase transitions. Among other things, this method allows an efficient computation of the minimum energy path (MEP) during the nucleation process. The MEP provides information about the size, shape and free energy barrier of the critical nucleus. To improve the efficiency of the string method, a special initialization process is proposed. Constraints from physical models are treated using two methods, a generalized coordinates method and a projection method. Strategies for choosing the computational domain and defining the nucleus boundary are also introduced. The validity of our approach is illustrated by two nontrivial examples from soft condensed matter physics, namely the nematic–isotropic transition of liquid crystals and the ordered-to-ordered phase transition of diblock copolymers.     

A spatial correlation method for reconstructing blood flow velocity vector in the nonlinear tomography

A method for reconstructing the total blood flow velocity vector distribution is proposed. Experimental data are obtained with the help of the second-order nonlinearity parameter tomography scheme that employs wideband coded signals and a small number of transducers. The data are processed using the selective matched filtering technique in combination with spatial correlation of fragments of the speckle structure that appears in the process of reconstruction of the spatial nonlinearity parameter distribution. Results of numerical simulations are presented.     

基于Burton-Miller方程的轴对称结构声学边界元方法

提出了一种求取轴对称结构任意边界条件下声辐射特性的边界元方法。采用Burton和Miller改进型公式将高阶奇异项转化为弱奇异项之和,保证声辐射参数的唯一性,且计算简单精确。将结构表面声压与振速按照旋转轴角度进行Fourier级数展开,利用级数的正交性建立各项待定系数的求解公式;然后转化格林函数的法向偏导为切向偏导,方便直接计算各项积分,并将面积分公式表示为沿结构边界的线积分和沿旋转角度的积分;进一步采用二次等参单元离散结构边界线,建立声压与振速的关系矩阵,从而确定结构声辐射参数。以脉动球源和横向振动球源为例计算,与解析解和传统边界元法结果作对比,说明该方法的有效精确性。     

A numerical algorithm for viscous incompressible interfacial flows

We present a new algorithm to numerically simulate two-dimensional viscous incompressible flows with moving interfaces. The motion is updated in time by using the backward difference formula through an iterative procedure. At each iteration, the pseudo-spectral technique is applied in the horizontal direction. The resulting semi-discretized equations constitute a boundary value problem in the vertical coordinate which is solved by decoupling growing and decaying solutions. Numerical tests justify that this method achieves fully second-order accuracy in both the temporal variable and vertical coordinate. As an application of this algorithm, we study the motion of Stokes waves in the presence of viscosity. Our numerical results are consistent with the recently published asymptotic solution for Stokes waves in slightly viscous fluids.     

A method based on a nonlinear generalized Heisenberg algebra to study the molecular vibrational spectrum

We propose a method, based on a generalized Heisenberg algebra (GHA), to reproduce the anharmonic spectrum of diatomic molecules. The theoretical spectrum generated by GHA allows us to fit the experimental data and to obtain the dissociation energy for the carbon monoxide molecule. Our outcomes are more accurate than the standard models used to study molecular vibrations, namely the Morse and the q-oscillator models and comparable to the perturbed Morse model proposed by Huffaker [CITE], for the first experimental levels. The dissociation energy obtained here is more accurate than all previous models.     

用于光谱分类的光谱峰谷沿匹配算法

针对超短脉冲激光与气体相互作用产生的非线性荧光光谱的特点,提出了一种新的光谱去噪、特征提取、特征匹配的有效方法.通过对三种气体不同浓度的荧光光谱进行小波压缩去噪,得到5层664个小波系数,重构后获得除噪效果满意而特征信息保持不变的光谱.对重构后的光谱进行光谱峰/谷沿相似度的分析,构造了三种气体的特征谱.在此基础上,经对三种气体光谱进行特征峰/谷沿的匹配分析,提出一种利用同种物质不同浓度下的光谱总体相似度来构造标准特征谱新方法,通过与标准特征谱的匹配分析得到分类结果.实验结果表明,该算法的分类结果受浓度变化影响小,对低至0.05%浓度条件下的全部27条光谱均获得100%的分类识别结果.     

A numerical study of the applicability of Melnikov's method

Melnikov's method can be used to predict the appearance of the homoclinic and also the heteroclinic tangency. The applicability of the method is studied by using numerically calculated invariant manifolds of the Poincaré map of a soft spring oscillator.     

A numerical study of the Fermi hypernetted chain method

We compare several elaborations of the Fermi hypernetted chain (FHNC) method for the calculation of ground-state properties of dense Fermi systems. A variety of simple nuclear potentials is employed. Consistency tests give estimates for the convergence of the resulting FHNC expansions with respect to elementary diagrams, and to antisymmetry corrections.