An all-speed relaxation scheme for interface flows with surface tension

We consider interface flows where compressibility and capillary forces (surface tension) are significant. These flows are described by a non-conservative, unconditionally hyperbolic multiphase model. The numerical approximation is based on finite-volume method for unstructured grids. At the discrete level, the surface tension is approximated by a volume force (CSF formulation). The interface physical properties are recovered by designing an appropriate linearized Riemann solver (Relaxation scheme) that prevents spurious oscillations near material interfaces. For low-speed flows, a preconditioning linearization is proposed and the low Mach asymptotic is formally recovered. Numerical computations, for a bubble equilibrium, converge to the required Laplace law and the dynamic of a drop, falling under gravity, is in agreement with experimental observations.     

An adaptive SPH method for strong shocks

We propose an alternative SPH scheme to usual SPH Godunov-type methods for simulating supersonic compressible flows with sharp discontinuities. The method relies on an adaptive density kernel estimation (ADKE) algorithm, which allows the width of the kernel interpolant to vary locally in space and time so that the minimum necessary smoothing is applied in regions of low density. We have performed a von Neumann stability analysis of the SPH equations for an ideal gas and derived the corresponding dispersion relation in terms of the local width of the kernel. Solution of the dispersion relation in the short wavelength limit shows that stability is achieved for a wide range of the ADKE parameters. Application of the method to high Mach number shocks confirms the predictions of the linear analysis. Examples of the resolving power of the method are given for a set of difficult problems, involving the collision of two strong shocks, the strong shock-tube test, and the interaction of two blast waves.     

An algorithm and its application for obtaining some kind of infinite-dimensional Hamiltonian canonical formulation

Using factorization viewpoint of differential operator, this paper discusses how to transform a nonlinear evolution equation to infinite-dimensional Hamiltonian linear canonical formulation. It proves a sufficient condition of canonical factorization of operator, and provides a kind of mechanical algebraic method to achieve canonical evolution equation, infinite-dimensional Hamiltonian canonical system, factorization of differential operator, commutator,evolution equation, infinite-dimensional Hamiltonian canonical system, factorization of differential operator, commutatorProject supported by the National Natural Science Foundation of China (Grant No~10562002) and the Natural Science Foundation of Nei Mongol, China (Grant No~200508010103).2007-05-09{\partial}/{\partial x}Corresponding author．E-mail：alatanca@imu．edu.cn/qk/85823A/200711/25754042.html0200, 03405/9/2007 12:00:00 AMUsing factorization viewpoint of differential operator, this paper discusses how to transform a nonlinear evolution equation to infinite-dimensional Hamiltonian linear canonical formulation. It proves a sufficient condition of canonical factorization of operator, and provides a kind of mechanical algebraic method to achieve canonical $`{\partial}/{\partial x}'$-type expression, correspondingly. Then three examples are given, which show the application of the obtained algorithm. Thus a novel idea for inverse problem can be derived feasibly.http://cpb.iphy.ac.cn/CN/10.1088/1009-1963/16/11/002https://cpb.iphy.ac.cn/CN/article/downloadArticleFile.do?attachType=PDF&id=1088382007-11-20'-type expression, correspondingly. Then three examples are given, which show the application of the obtained algorithm. Thus a novel idea for inverse problem can be derived feasibly.     

Һ̬�������ɱ���������������ѧЧӦ����ֵģ��

分析了磁流体力学效应对液态金属自由表面射流稳定性的影响。从射流的感应电势、电流、速度等方面,解释射流在磁场中稳定的原因。数值计算结果验证了理论分析。     

液态金属自由表面膜流MHD效应的数值模拟

对液态金属自由表面膜流在强磁场下的磁流体力学效应进行了数值模拟研究,获得了液态金属自由表面的形状、截面流速分布及截面上的电动势分布,从而能对膜流的一些磁流体动力学行为作出解释。数值计算结果与理论分析和实验结果符合较好。由实验和数值模拟结果可以得出,液态金属膜流通过强磁场时,磁场会阻碍膜流的运动。     

A semi-coherent interface with nanoledges intersecting the free surface of an elastic half-space

Abstract

The elastic field generated by a semi-coherent interface formed by a succession of parallel terraces and nanoledges is investigated when the interface intersects the free surface of an isotropic half-space medium. The method uses the combination of semi-infinite Somigliana dislocations placed on each terrace and semi-infinite translation (or Volterra) dislocations placed at one or two corners of the nanoledges. Numerical examples depict the deformation fields nearby emerging regions of two semi-coherent interfaces; one is flat, while the other includes nanoledges.     

An interface capturing method for the simulation of multi-phase compressible flows

A novel finite-volume interface (contact) capturing method is presented for simulation of multi-component compressible flows with high density ratios and strong shocks. In addition, the materials on the two sides of interfaces can have significantly different equations of state. Material boundaries are identified through an interface function, which is solved in concert with the governing equations on the same mesh. For long simulations, the method relies on an interface compression technique that constrains the thickness of the diffused interface to a few grid cells throughout the simulation. This is done in the spirit of shock-capturing schemes, for which numerical dissipation effectively preserves a sharp but mesh-representable shock profile. For contact capturing, the formulation is modified so that interface representations remain sharp like captured shocks, countering their tendency to diffuse via the same numerical diffusion needed for shock-capturing. Special techniques for accurate and robust computation of interface normals and derivatives of the interface function are developed. The interface compression method is coupled to a shock-capturing compressible flow solver in a way that avoids the spurious oscillations that typically develop at material boundaries. Convergence to weak solutions of the governing equations is proved for the new contact capturing approach. Comparisons with exact Riemann problems for model one-dimensional multi-material flows show that the interface compression technique is accurate. The method employs Cartesian product stencils and, therefore, there is no inherent obstacles in multiple dimensions. Examples of two- and three-dimensional flows are also presented, including a demonstration with significantly disparate equations of state: a shock induced collapse of three-dimensional van der Waal’s bubbles (air) in a stiffened equation of state liquid (water) adjacent to a Mie-Grüneisen equation of state wall (copper).     

SPH二阶粒子近似光滑函数的充分条件及数值验证

为实现SPH二阶粒子近似,讨论光滑函数应满足的充分条件,给出对应的光滑函数形式,分析二阶粒子近似解的唯一性,并进行数值验证.结果表明：和标准SPH方法相比,在周期边界条件下,随着粒子数加密,SPH方法二阶粒子近似的L₁误差减小且L₁阶趋近2.0.     

The potential energy surface and chaos in 2D Hamiltonian systems

We provide a new insight into the relationship between the geometric property of the potential energy surface and chaotic behavior of 2D Hamiltonian dynamical systems, and give an indicator of chaos based on the geometric property of the potential energy surface by defining Mean Convex Index (MCI). We also discuss a model of unstable Hamiltonian in detail, and show our results in good agreement with HBLSL's (Horwitz, Ben Zion, Lewkowicz, Schiffer and Levitan) new Riemannian geometric criterion.     

An integrable Hamiltonian hierarchy and associated integrable couplings system

This paper establishes a new isospectral problem. By making use of the Tu scheme, a new integrable system is obtained. It gives integrable couplings of the system obtained. Finally, the Hamiltonian form of a binary symmetric constrained flow of the system obtained is presented.     

An immersed interface method for Stokes flows with fixed/moving interfaces and rigid boundaries

We present an immersed interface method for solving the incompressible steady Stokes equations involving fixed/moving interfaces and rigid boundaries (irregular domains). The fixed/moving interfaces and rigid boundaries are represented by a number of Lagrangian control points. In order to enforce the prescribed velocity at the rigid boundaries, singular forces are applied on the fluid at these boundaries. The strength of singular forces at the rigid boundary is determined by solving a small system of equations. For the deformable interfaces, the forces that the interface exerts on the fluid are calculated from the configuration (position) of the deformed interface. The jumps in the pressure and the jumps in the derivatives of both pressure and velocity are related to the forces at the fixed/moving interfaces and rigid boundaries. These forces are interpolated using cubic splines and applied to the fluid through the jump conditions. The positions of the deformable interfaces are updated implicitly using a quasi-Newton method (BFGS) within each time step. In the proposed method, the Stokes equations are discretized via the finite difference method on a staggered Cartesian grid with the incorporation of jump contributions and solved by the conjugate gradient Uzawa-type method. Numerical results demonstrate the accuracy and ability of the proposed method to simulate incompressible Stokes flows with fixed/moving interfaces on irregular domains.     

Edge misfit dislocation formation at the interface of a nanopore and infinite substrate with surface/interface effects

The model of an edge misfit dislocation at the interface of the hollow nanopore and the infinite substrate with surface/interface stress is investigated. Using the complex variable method, analytical solutions for complex potentials of a film due to an edge misfit dislocation located in the film with surface/interface effect are derived, and the stress fields of the film and the edge misfit dislocation formation energy can be obtained. The critical conditions for edge misfit dislocation formation are given at which the generation of an edge misfit dislocation is energetically favourable. The influence of the ratio of the shear modulus between the film and the infinite substrate, the misfit strain, the radius of the nanopore and the surface/interface stress on the critical thickness of the film is discussed.     

An energy-consistency-preserving large eddy simulation-scalar filtered mass density function (LES-SFMDF) method for high-speed flows

An energy-consistency-preserving large eddy simulation-scalar filtered mass density function (LES-SFMDF) method is developed to improve the existing LES-SFMDF method in high-speed flows, especially supersonic flows. The high-speed source term in the SFMDF transport equation is analysed theoretically from a new point of view, and then several primary principles are proposed for the LES-SFMDF to achieve a good consistency even along discontinuities. Based on these principles and further theoretical analysis, the high-speed source term of the enthalpy in the SFMDF is modelled and computed from both MC and LES variables rather than the usually used solely LES variables. This new LES-SFMDF method is used for simulating the flows in a shock tube and in a subsonic temporally developing mixing layer. This method shows a better particle energy consistency than the existing method when applied across discontinuities in supersonic laminar flows. Unlike the existing method, with this energy-consistency-preserving LES-SFMDF method, particle energy consistency is consistent with particle mass consistency so that particle energy consistency can benefit from particle velocity correction. This method also demonstrates robustness for various numbers of particles.     

An integrable Hamiltonian hierarchy, a high-dimensional loop algebra and associated integrable coupling system

A subalgebra of loop algebra ?_2 is established. Therefore, a new isospectral problem is designed. By making use of Tu's scheme, a new integrable system is obtained, which possesses bi-Hamiltonian structure. As its reductions, a formalism similar to the well-known Ablowitz－Kaup－Newell－Segur (AKNS) hierarchy and a generalized standard form of the Schr?dinger equation are presented. In addition, in order for a kind of expanding integrable system to be obtained, a proper algebraic transformation is supplied to change loop algebra ?_2 into loop algebra ?_1. Furthermore, a high-dimensional loop algebra is constructed, which is different from any previous one. An integrable coupling of the system obtained is given. Finally, the Hamiltonian form of a binary symmetric constrained flow of the system obtained is presented.     

锂离子电池中的尺寸效应与表界面问题研究

文章介绍了锂离子电池中的尺寸效应与表界面问题的研究.从热力学与动力学方面讨论了小尺寸材料与体材料性质的区别以及作为储能材料的优缺点.从实际考虑,建议发展动力学稳定的纳米结构电极材料.     

Canonical transformations and exact invariants for time‐dependent Hamiltonian systems

An exact invariant is derived for n‐degree‐of‐freedom non‐relativistic Hamiltonian systems with general time‐dependent potentials. To work out the invariant, an infinitesimalcanonical transformation is performed in the framework of the extended phase‐space. We apply this approach to derive the invariant for a specific class of Hamiltonian systems. For the considered class of Hamiltonian systems, the invariant is obtained equivalently performing in the extended phase‐space a finitecanonical transformation of the initially time‐dependent Hamiltonian to a time‐independent one. It is furthermore shown that the invariant can be expressed as an integral of an energy balance equation. The invariant itself contains a time‐dependent auxiliary function ξ (t) that represents a solution of a linear third‐order differential equation, referred to as the auxiliary equation. The coefficients of the auxiliary equation depend in general on the explicitly known configuration space trajectory defined by the system's time evolution. This complexity of the auxiliary equation reflects the generally involved phase‐space symmetry associated with the conserved quantity of a time‐dependent non‐linear Hamiltonian system. Our results are applied to three examples of time‐dependent damped and undamped oscillators. The known invariants for time‐dependent and time‐independent harmonic oscillators are shown to follow directly from our generalized formulation.     

Lie symmetries and non-Noether conserved quantities for Hamiltonian canonical equations

This paper focuses on studying Lie symmetries and non-Noether conserved quantities of Hamiltonian dynamical systems in phase space. Based on the infinitesimal transformations with respect to the generalized coordinates and generalized momenta, we obtain the determining equations and structure equation of the Lie symmetry for Hamiltonian dynamical systems. This work extends the research of non-Noether conserved quantity for Hamilton canonical equations, and leads directly to a new type of non-Noether conserved quantities of the systems. Finally, an example is given to illustrate these results.     

Stability analysis on the free surface phenomena of a magnetic fluid for general use

This paper presents an analysis for elucidating a variety of physical processes on the interface (free surface) of magnetic fluid. The present analysis is composed of the magnetic and the fluid analysis, both of which have no limitations concerning the interface elevation or its profile. The magnetic analysis provides rigorous interface magnetic field under arbitrary distributions of applied magnetic field. For the fluid analysis, the equation for interface motion includes all nonlinear effects. Physical quantities such as the interface magnetic field or the interface stresses, obtained first as the wavenumber components, facilitate confirming the relations with those by the conventional theoretical analyses. The nonlinear effect is formulated as the nonlinear mode coupling between the interface profile and the applied magnetic field. The stability of the horizontal interface profile is investigated by the dispersion relation, and summarized as the branch line. Furthermore, the balance among the spectral components of the interface stresses are shown, within the sufficient range of the wavenumber space.     

Dynamic behavior of a single bubble between the free surface and rigid wall

The objective of this paper is to apply high-speed photography and schlieren method to investigate the bubble dynamics between the free surface and a rigid wall. The temporal evolution of the bubble shape and the free surface motion are recorded by two synchronous high-speed cameras. Experiments are carried out for a single bubble generated at various normalized stand-off distances from bubble center to the free surface and to the rigid wall. The results show that (1) three distinctive patterns are identified with the morphology of the bubble and free surface, namely single toroidal bubble without spike (STB), single toroidal bubble with a spike (STBS) and double toroidal bubbles with a spike (DTBS). (2) The dynamic characteristics of the bubble at collapse and rebound stage vary evidently at different patterns, including the bubble shape variations and free surface motion. In detail, the schlieren images show the formation and propagation of shock waves, which explains the radiative process of bubble collapse energy. (3) Qualitative comparisons are carried out for the bubble and free surface at the same pattern. And quantitative analyses are conducted for the jet velocity, bubble collapse position, bubble collapse time and spike height, etc. for different values of bubble-rigid wall distance.