Simulation of viscous flows with undulatory boundaries. Part I: Basic solver

A numerical method for the simulation of viscous flows with undulatory walls and free surfaces is presented. The simulation domain is discretized by a boundary-fitted and time-dependent grid. The Navier–Stokes equations, subject to fully nonlinear kinematic and dynamic boundary conditions at the free surface and no-slip boundary condition at the wall, are simulated by a hybrid pseudo-spectral and finite difference method in space and a semi-implicit fractional-step method in time. The performance of the method is demonstrated by a series of test cases including flows over wavy boundaries, various surface waves, and interaction between vortices and free surfaces. Validation by convergence test and extensive comparisons with previous theoretical, experimental, and numerical studies indicate the accuracy and efficiency of the method. Finally, a simulation example of turbulence and free surface interaction is presented. Results show that the rich features of the free surface such as surface waves, splats, anti-splats, dimples, and scars are captured accurately. Characteristic vortical structures and variation of turbulence statistics in the near-surface region are also elucidated.     

Kinetic theory based lattice Boltzmann equation with viscous dissipation and pressure work for axisymmetric thermal flows

A lattice Boltzmann equation (LBE) for axisymmetric thermal flows is proposed. The model is derived from the kinetic theory which exhibits several features that distinguish it from other previous LBE models. First, the present thermal LBE model is derived from the continuous Boltzmann equation, which has a solid foundation and clear physical significance; Second, the model can recover the energy equation with the viscous dissipation term and work of pressure which are usually ignored by traditional methods and the existing thermal LBE models; Finally, unlike the existing thermal LBE models, no velocity and temperature gradients appear in the force terms which are easy to realize in the present model. The model is validated by thermal flow in a pipe, thermal buoyancy-driven flow, and swirling flow in vertical cylinder by rotating the top and bottom walls. It is found that the numerical results agreed excellently with analytical solution or other numerical results.     

A special class of stationary flows for two-dimensional euler equations: A statistical mechanics description. Part II

We continue and conclude our analysis started in Part I (see [CLMP]) by discussing the microcanonical Gibbs measure associated to a N-vortex system in a bounded domain. We investigate the Mean-Field limit for such a system and study the corresponding Microcanonnical Variational Principle for the Mean-Field equation. We discuss and achieve the equivalence of the ensembles for domains in which we have the concentration at (–8)⁺ in the canonical framework. In this case we have the uniqueness of the solutions of the Mean-Field equation. For the other kind of domains, for large values of the energy, there is no equivalence, the entropy is not a concave function of the energy, and the Mean-field equation has more than one solution. In both situations, we have concentration when the energy diverges. The Microcanonical Mean Field Limit for the N-vortex system is proven in the case of equivalence of ensembles.Communicated by J.L. Lebowitz     

Kinetic theory of hydrodynamic flows. II. The drag on a sphere and on a cylinder

In this paper we continue the study of solutions of the extended Boltzmann equation started previously. In particular, we study an iterated solution of the equation that can be used to describe the flow of a rarefied gas around a macroscopic object. We discuss the rarefied flow and then show how the iterated solution can be extended into the hydrodynamic regime. The results for the drag force and for the distribution function of the gas molecules are shown to be identical to the results obtained in a previous paper by a generalization of the normal solution method. We also discuss the special properties of both rarefied and continuum flows around a cylinder and show that in both regions one must take into account Oseen-like terms which naturally appear in the extended Boltzmann equation. In the hydrodynamic regime we obtain Lamb's formula for the force on the cylinder. By relating the terms in the iterated expression to dynamical events taking place in the fluid, we are able to discuss the dynamical origin of the results obtained here.A preliminary report on the work described here and in Part I was given in Ref. 2.     

浅调制下双二波混频：II.耦合理论

通过耦合波方程的方法，考虑光折变晶体中存在基频光栅和最低阶的高频和组合光栅，本文给出了各种情形下的耦合波方程及其解，结果表明一对相干光束的耦合可以被另一束相干的总光强，相位和调制度等控制，理论上，利用这种效应，可以实现相位－振幅转换及振幅－振幅调制。     

Artificial compressibility, characteristics-based schemes for variable-density, incompressible, multispecies flows: Part II. Multigrid implementation and numerical tests

The paper presents an investigation of the accuracy and efficiency of artificial compressibility, characteristics-based (CB) schemes for variable-density incompressible flows. The CB schemes have been implemented in conjunction with a multigrid method for accelerating numerical convergence and a fourth-order, explicit Runge–Kutta method for the integration of the governing equations in time. The implementation of the CB schemes is obtained in conjunction with first-, second- and third-order interpolation formulas for calculating the variables at the cell faces of the computational volume. The accuracy and efficiency of the schemes are examined against analytical and experimental results for diffusion broadening in two- and three-dimensional microfluidic channels, a problem that has motivated the development of the present methods. Moreover, unsteady, inviscid simulations have been performed for variable-density mixing layer. The computations revealed that accuracy and efficiency depend on the CB scheme design. The best multigrid convergence rates were exhibited by the conservative CB scheme, which is obtained by the fully conservative formulation of the variable-density, incompressible equations.     

Dual nature of localization in guiding systems with randomly corrugated boundaries: Anderson-type versus entropic

A unified theory for the conductance of an infinitely long multimode quantum wire whose finite segment has randomly rough lateral boundaries is developed. It enables one to rigorously take account of all feasible mechanisms of wave scattering, both related to boundary roughness and to contacts between the wire rough section and the perfect leads within the same technical frameworks. The rough part of the conducting wire is shown to act as a mode-specific randomly modulated effective potential barrier whose height is governed essentially by the asperity slope. The mean height of the barrier, which is proportional to the average slope squared, specifies the number of conducting channels. Under relatively small asperity amplitude this number can take on arbitrary small, up to zero, values if the asperities are sufficiently sharp. The consecutive channel cut-off that arises when the asperity sharpness increases can be regarded as a kind of localization, which is not related to the disorder per se but rather is of entropic or (equivalently) geometric origin. The fluctuating part of the effective barrier results in two fundamentally different types of guided wave scattering, viz., inter- and intramode scattering. The intermode scattering is shown to be for the most part very strong except in the cases of (a) extremely smooth asperities, (b) excessively small length of the corrugated segment, and (c) the asperities sharp enough for only one conducting channel to remain in the wire. Under strong intermode scattering, a new set of conducting channels develops in the corrugated waveguide, which have the form of asymptotically decoupled extended modes subject to individual solely intramode random potentials. In view of this fact, two transport regimes only are realizable in randomly corrugated multimode waveguides, specifically, the ballistic and the localized regime, the latter characteristic of one-dimensional random systems. Two kinds of localization are thus shown to coexist in waveguide-like systems with randomly corrugated boundaries, specifically, the entropic localization and the one-dimensional Anderson (disorder-driven) localization. If the particular mode propagates across the rough segment ballistically, the Fabry–Pérot-type oscillations should be observed in the conductance, which are suppressed for the mode transferred in the Anderson-localized regime.     

浅调制下双二波混频: I.耦合理论

通过耦合波方程的方法，考虑光折变晶体中存在基频光栅和最低阶的高频和组合光栅，本文给出了各种情形下的耦合波方程及其解。结果表明一对相干光束的耦合可以被另一束相干的总光强，相位和调制度等控制。理论上，利用这种效应，可以实现相位－振幅转换及振幅－振幅调制     

Absorbing boundary conditions for scalar waves in anisotropic media. Part 1: Time harmonic modeling

With the ultimate goal of devising effective absorbing boundary conditions (ABCs) for general anisotropic media, we investigate the accuracy aspects of local ABCs designed for the scalar anisotropic wave equation in the frequency domain (time harmonic case). The ABC analyzed in this paper is the perfectly matched discrete layers (PMDL). PMDL is a simple variant of perfectly matched layers (PML) and is equivalent to rational approximation-based local ABCs. Specifically, we derive a sufficient condition for PMDL to accurately absorb wave modes with outgoing group velocities and this condition turns out to be a simple bound on the PMDL parameters. The reflection coefficient derived in this paper clearly reveals that the PMDL absorption is based on group velocities, and not phase velocities, and hence a PMDL can be designed to correctly identify and accurately absorb all outgoing wave modes (even those with opposing signs of phase and group velocities). The validity of the sufficient condition is demonstrated through a series of frequency domain simulations. In part 2 of this paper [S. Savadatti, M.N. Guddati, Absorbing boundary conditions for scalar waves in anisotropic media. Part 2: Time-dependent modeling, J. Comput. Phys. (2010), http://dx.doi.org/10.1016/j.jcp.2010.05.017], the accuracy condition presented here is shown to govern both the well-posedness and accuracy aspects of PMDL designed for transient (time-dependent) modeling of scalar waves in anisotropic media.     

Absorbing boundary conditions for scalar waves in anisotropic media. Part 2: Time-dependent modeling

With the ultimate goal of devising effective absorbing boundary conditions (ABCs) for general anisotropic media, we investigate the well-posedness and accuracy aspects of local ABCs designed for the transient modeling of the scalar anisotropic wave equation. The ABC analyzed in this paper is the perfectly matched discrete layers (PMDL), a simple variant of perfectly matched layers (PML) that is also equivalent to rational approximation based ABCs. Specifically, we derive the necessary and sufficient condition for the well-posedness of the initial boundary value problem (IBVP) obtained by coupling an interior and a PMDL ABC. The derivation of the reflection coefficient presented in a companion paper (S. Savadatti, M.N. Guddati, J. Comput. Phys., 2010, doi:10.1016/j.jcp.2010.05.018) has shown that PMDL can correctly identify and accurately absorb outgoing waves with opposing signs of group and phase velocities provided the PMDL layer lengths satisfy a certain bound. Utilizing the well-posedness theory developed by Kreiss for general hyperbolic IBVPs, and the well-posedness conditions for ABCs derived by Trefethen and Halpern for isotropic acoustics, we show that this bound on layer lengths also ensures well-posedness. The time discretized form of PMDL is also shown to be theoretically stable and some instability related to finite precision arithmetic is discussed.     

Surface flows of granular mixtures: II. Segregation with grains of different size

We present the generalization of a theoretical model for segregation of granular mixtures due to surface flows, published in J. Phys. I France 6, 1295 (1996). Our generalized model is valid for grains differing by their size and/or their surface properties; in the present paper, we describe the case of two species with the same surface properties but two different sizes. The rolling stream is assumed to be homogeneous. Exchanges between the grains at rest and the rolling stream are modelized via binary collisions. The model predicts that during the filling of a two-dimensional silo, continuous segregation appears inside the static phase: small (respectively large) grains tend to stop uphill (respectively downhill), although both species remain present everywhere. This fits the observations when the size difference between the species is small. When the size difference is large, a different regime is observed. We argue that in this case, segregation occurs directly inside the rolling stream. Received: 25 February 1998 / Received in final form and Accepted: 6 July 1998     

基于大密度差多相流SPH方法的二维液滴碰撞数值模拟

该文结合了Ott提出的修正连续性方程和Adami改进的动量方程, 对空气中的液滴碰撞问题进行了二维数值模拟. 为有效提高计算精度, 推导了适用于大密度差多相流的人工黏性和人工应力方程. 通过表面张力作用下方形液滴自然变化和空气中两液滴互溶的算例, 验证了算法的有效性; 对不同韦伯数 (8.8, 19.8)、不同碰撞参数 (0, 0.5)下的液滴碰撞过程进行了数值模拟, 并与VOF方法对比,取得了较为一致的结果; 进一步计算多个韦伯数、多个碰撞参数下的液滴碰撞, 得到了空气中二维液滴碰撞结果分布图,与实验结果相符合. 结果表明, 该算法对于求解涉及大密度差多相流的液滴碰撞破碎问题十分有效,而且该方法容易拓展到三维, 从而为进一步模拟火箭发动机的二次雾化过程奠定了基础. 关键词： 光滑粒子流体动力学 大密度差 多相流 液滴碰撞     

Characteristic boundary conditions for simulations of compressible reacting flows with multi-dimensional,viscous and reaction effects

A generalized formulation of the characteristic boundary conditions for compressible reacting flows is proposed. The new and improved approach resolves a number of lingering issues of spurious solution behaviour encountered in turbulent reacting flow simulations in the past. This is accomplished (a) by accounting for all the relevant terms in the determination of the characteristic wave amplitudes and (b) by accommodating a relaxation treatment for the transverse gradient terms with the relaxation coefficient properly determined by the low Mach number asymptotic expansion. The new boundary conditions are applied to a comprehensive set of test problems including: vortex-convection; turbulent inflow; ignition front propagation; non-reacting and reacting Poiseuille flows; and counterflow cases. It is demonstrated that the improved boundary conditions perform consistently superior to existing approaches, and result in robust and accurate solutions with minimal acoustic wave interactions at the boundary in hostile turbulent combustion simulation conditions.     

The Simulation of Electrostatic Spray Painting Process with High‐Speed Rotary Bell Atomizers. Part II: External Charging

The present contribution summarizes investigations aiming to completely model the electrostatically supported spray painting process with high‐speed rotary bells by means of CFD. In this part II, so‐called external charging atomizers, where high voltage is applied to emitting electrode needles, are considered. Here, charging of the droplets takes place due to free ions produced from corona discharge at the electrodes. Part I [1] dealt with direct charging atomizers, where potential is applied directly to the rotating bell. The commercial CFD‐code Fluent has been extended to account for the electrostatic field and the space charge effect due to the ions. Here, a model for the time‐dependent and inhomogeneous field charging of the droplets was applied. Furthermore, the direct interaction between the ion current and the flow field, i.e., the so‐called ion wind, could be calculated. As input conditions, the airflow from the shaping air orifices and measured droplet sizes close to the bell edge using Fraunhofer diffraction were taken. In general, numerical and experimental results are in good agreement, confirming the applicability of the chosen physical approach. This is especially true for the final film thickness on the target and the transfer efficiency, i.e., the amount of paint solids that finally deposited on the target. In the near future the calculations must be extended to true unsteady simulations including the dynamic meshing procedure.     

A current density conservative scheme for incompressible MHD flows at a low magnetic Reynolds number. Part II: On an arbitrary collocated mesh

A conservative formulation of the Lorentz force is given here for magnetohydrodynamic (MHD) flows at a low magnetic Reynolds number with the current density calculated based on Ohm’s law and the electrical potential formula. This conservative formula shows that the total momentum contributed from the Lorentz force is conservative when the applied magnetic field is constant. For the case with a non-constant applied magnetic field, the Lorentz force has been divided into two parts: a strong globally conservative part and a weak locally conservative part.The conservative formula has been employed to develop a conservative scheme for the calculation of the Lorentz force on an unstructured collocated mesh. Only the current density fluxes on the cell faces, which are calculated using a consistent scheme with good conservation, are needed for the calculation of the Lorentz force. Meanwhile, a conservative interpolation technique is designed to get the current density at the cell center from the current density fluxes on the cell faces. This conservative interpolation can keep the current density at the cell center conservative, which can be used to calculate the Lorentz force at the cell center with good accuracy. The Lorentz force calculated from the conservative current at the cell center is equivalent to the Lorentz force from the conservative formula when the applied magnetic field is constant, which can conserve the total momentum. We will further prove that the simple interpolation scheme used in the Part I [M.-J. Ni, R. Munipalli, N.B. Morley, P.Y. Huang, M. Abdou, A current density conservative scheme for MHD flows at a low magnetic Reynolds number. Part I. On a rectangular collocated grid system, Journal of Computational Physics, in press, doi:10.1016/j.jcp.2007.07.025] of this series of papers is conservative on a rectangular grid and can keep the total momentum conservative in a rectangular grid.     

Scalar conservation laws with discontinuous flux function: II. On the stability of the viscous profiles

The equation , whereH is Heaviside's step function, appears for example in continuous sedimentation of solid particles in a liquid, in two-phase flow, in traffic-flow analysis and in ion etching. The discontinuity of the flux function atx=0 causes a discontinuity of a solution, which is not uniquely determined by the initial data. By a viscous profile of this discontinuity we mean a stationary solution ofu _t +(F ) _x =u _xx , whereF is a smooth approximation of the discontinuous flux, i.e.,H is smoothed. We present some results on the stability of the viscous profiles, which means that a small disturbance tends to zero uniformly ast. This is done by weighted energy methods, where the weight (depending onf andg) plays a crucial role.     

An adaptive version of ghost-cell immersed boundary method for incompressible flows with complex stationary and moving boundaries

An adaptive version of immersed boundary method for simulating flows with complex stationary and moving boundaries is presented.The method employs a ghost-cell methodology which allows for a sharp representation of the immersed boundary.To simplify the implementation of the methodology,a volume-of-fluid method is introduced to identify the immersed boundary.In addition,the domain is spatially discretized using a tree-based discretization which is relatively simple to implement a fully flexible adaptive refi...     

Low temperature spectroscopy of proteins. Part II: Experiments with single protein complexes

In this part of the review we describe aspects of the physics of proteins at low temperature as they are reflected in the spectra of individual pigment–protein complexes. The focus of this review is on the spectral diffusion of chromophores that are naturally embedded in light-harvesting complexes from purple bacteria. From the spectral diffusion behaviour we can deduce details about the organisation of the energy landscape of the protein and discuss the implications for the motions of the protein in conformational phase space.     

Finite-aperture effects for Fourier transform systems with convergent illumination. Part II: 3-D system analysis

One of the most important optical signal processing operations is the optical Fourier transform (OFT). Of the arrangements for implementation of the OFT, perhaps the most flexible is that for the scaled optical Fourier transform (SOFT), as it allows control over the scale of the output Fourier transform distribution. By means of an analysis in cylindrical coordinates, we examine some of the practical limits introduced by the use of a thin lens of finite aperture in the implementation of the SOFT. Using simple rules of thumb that are based on an examination of the phase and magnitude deviations from the ideal (infinite-lens) diameter case, we define a volume inside the geometric shadow, which we refer to as a sub-geometric shadow. We then show that inside this sub-geometric shadow errors introduced by diffraction can be quantified.