Modified kinetic flux vector splitting schemes for compressible flows

We investigate the traditional kinetic flux vector splitting (KFVS) and BGK schemes for the compressible Euler equations. First, based on a careful study of the behavior of the discrete physical variables across the contact discontinuity, we analyze quantitatively the mechanism of inducing spurious oscillations of the velocity and pressure in the vicinity of the contact discontinuity for the first-order KFVS and BGK schemes. Then, with the help of this analysis, we propose a first-order modified KFVS (MKFVS) scheme which is oscillation-free in the vicinity of the contact discontinuity, provided certain consistent conditions are satisfied. Moreover, by using piecewise linear reconstruction and van Leer’s limiter, the first-order MKFVS scheme is extended to a second-order one, consequently, a nonoscillatory second-order MKFVS scheme is constructed. Finally, by combing the MKFVS schemes with the γ-model, we successfully extend the MKFVS schemes to multi-flows, and propose therefore a first- and second-order MKFVS schemes for multi-fluid computations, which are nonoscillatory across fluid interfaces. A number of numerical examples presented in this paper validate the theoretic analysis and demonstrate the good performance of the MKFVS schemes in simulation of contact discontinuities for both single- and multi-fluids.     

Population dynamics in compressible flows

Organisms often grow, migrate and compete in liquid environments, as well as on solid surfaces. However, relatively little is known about what happens when competing species are mixed and compressed by fluid turbulence. In these lectures we review our recent work on population dynamics and population genetics in compressible velocity fields of one and two dimensions. We discuss why compressible turbulence is relevant for population dynamics in the ocean and we consider cases both where the velocity field is turbulent and when it is static. Furthermore, we investigate populations in terms of a continuos density field and when the populations are treated via discrete particles. In the last case we focus on the competition and fixation of one species compared to another.     

Multifractal clustering in compressible flows

A quantitative relationship is found between the multifractal properties of the asymptotic mass distribution in a random dissipative system and the long-time fluctuations of the local stretching rates of the dynamics. It captures analytically the fine aspects of the strongly intermittent clustering of dynamical trajectories. Applied to a simple compressible hydrodynamical model with known stretching-rate statistics, the relation produces a nontrivial spectrum of multifractal dimensions that is confirmed numerically.     

Population genetics in compressible flows

We study competition between two biological species advected by a compressible velocity field. Individuals are treated as discrete Lagrangian particles that reproduce or die in a density-dependent fashion. In the absence of a velocity field and fitness advantage, number fluctuations lead to a coarsening dynamics typical of the stochastic Fisher equation. We investigate three examples of compressible advecting fields: a shell model of turbulence, a sinusoidal velocity field and a linear velocity sink. In all cases, advection leads to a striking drop in the fixation time, as well as a large reduction in the global carrying capacity. We find localization on convergence zones, and very rapid extinction compared to well-mixed populations. For a linear velocity sink, one finds a bimodal distribution of fixation times. The long-lived states in this case are demixed configurations with a single interface, whose location depends on the fitness advantage.     

Three-dimensional boundary conditions for numerical simulations of reactive compressible flows with complex thermochemistry

The Navier–Stokes characteristic boundary conditions (NSCBC) is a very efficient numerical strategy to treat boundary conditions in fully compressible solvers. The present work is an extension of the 3D-NSCBC method proposed by Yoo et al. and Lodato et al. in order to account for multi-component reactive flows with detailed chemistry and complex transport. A new approach is proposed for the outflow boundary conditions which enables clean exit of non-normal flows, and the specific treatment of all kinds of edges and corners is carefully addressed. The proposed methodology is successfully validated on various challenging multi-component reactive flow configurations.     

On the relation between finite element and finite volume schemes for compressible flows with cylindrical and spherical symmetry

A complete set of equivalence conditions, connecting the node-centered finite volume and the mass-lumped finite element schemes, is derived for the first time for the one-dimensional compressible Euler equations with cylindrical and spherical symmetry. Analytical expressions for the evaluation of the equivalent cell volumes and interface normals in terms of the finite element integrals are presented. Numerical experiments for compressible unsteady flows, including expanding and converging shock problems, are carried out using the new approach and the differences with the results from a finite volume scheme violating the equivalence conditions are discussed.     

可压流体Rayleigh-Taylor不稳定性的离散Boltzmann模拟

使用离散Boltzmann模型模拟了可压流体系统中多模初始情况下的Rayleigh-Taylor不稳定性.该离散Boltzmann模型等效于一个Navier-Stokes模型外加一个关于热动非平衡行为的粗粒化模型.通过模拟Riemann问题:Sod激波管、冲击波碰撞和热Couette流问题验证模型的有效性,所得数值结果与解析解一致.利用该模型对界面间断随机多模初始扰动的可压Rayleigh-Taylor不稳定性进行数值模拟研究,得到不稳定性界面演化过程的基本图像.由于黏性和热传导共同作用,一开始扰动界面被"抹平",演化较慢;随着模式互相耦合而减少,演化开始加速,并经历非线性小扰动阶段和不规则非线性阶段,而后发展成典型的"蘑菇状",后期进入湍流混合阶段.由于扰动模式的耦合与发展,轻重流体的重力势能、压缩能与动能相互转化,系统先是趋于热动平衡态,而后偏离热动平衡态以线性形式增长,接着再次趋于热动平衡态,最后慢慢远离热动平衡态.     

Low-frequency components and modulation processes in compressible cavity flows

The modulated character of cavity-flow oscillations is investigated through a so-called modulation analysis of spectral distributions. The approach is based on a recently proposed viewpoint on the Rossiter formula and concerns mid- to high-subsonic flows in shallow cavities. For this type of configuration, the spectra are mainly characterized by the presence of several dominant peaks (Rossiter modes) that are not in a harmonic relation but uniformly spaced at a distance equal to the fundamental frequency of the oscillation mechanism (aero-acoustic feedback loop). This feature is interpreted as the result of an amplitude modulation process and related to variations in the vortex-corner interaction in the downstream part of the cavity (γ modulation). A lower frequency modulation is identified through the secondary peak distribution. A detailed analysis of the spectral structure confirms the presence of a component at the corresponding frequency value (Δ_f mode). The assumption of a specific coupling between the two modulation processes is investigated. It leads to a new approximate form for the γ-modulation ratio that allows an explicit expression of the Rossiter constant γ related to aspect ratio of the cavity (L/D).     

A discrete model for compressible flows in heterogeneous media

This work deals with the building of a discrete model able to describe and to predict the evolution of complex gas flows in heterogeneous media. In many physical applications, large scales numerical simulation is no longer possible because of a lack of computing resources. Indeed the medium topology may be complex due to the presence of many obstacles (walls, pipes, equipments, geometric singularities etc.). Aircraft powerplant compartments are examples where topology is complex due to the presence of pipes, ducts, coolers and other equipment. Other important examples are gas explosions and large scale dispersion of hazardous materials in urban places, cities or underground involving obstacles such as buildings and various infrastructures. In all cases efficient safety responses are required. Then a new discrete model is built and solved in reasonable execution times for large cells volumes including such obstacles. Quantitative comparisons between experimental and numerical results are shown for different significant test cases, showing excellent agreement.     

Sticky particles in compressible flows: aggregation and Richardson's law

The evolution of the distance r(t) between a pair of passive tracer particles in rough compressible velocity fields is studied. The scaling behavior depends on the stickiness of the particles. Sticky particles start aggregating in moderately compressible flows, which can be realized on the free-slip surface of a turbulent fluid; nonsticky particles can aggregate only in less common strongly compressible flows (even then, the aggregation rate remains lower). Aggregation gives rise to an anomalous scaling law for the mean-square-distance growth rate, slower than Richardson's law. These findings help understand the results of recent experiments.     

A Brinkman penalization method for compressible flows in complex geometries

To simulate flows around solid obstacles of complex geometries, various immersed boundary methods had been developed. Their main advantage is the efficient implementation for stationary or moving solid boundaries of arbitrary complexity on fixed non-body conformal Cartesian grids. The Brinkman penalization method was proposed for incompressible viscous flows by penalizing the momentum equations. Its main idea is to model solid obstacles as porous media with porosity, ϕ, and viscous permeability approaching zero. It has the pronounced advantages of mathematical proof of error bound, strong convergence, and ease of numerical implementation with the volume penalization technique. In this paper, it is extended to compressible flows. The straightforward extension of penalizing momentum and energy equations using Brinkman penalization with respective normalized viscous, η, and thermal, η_T, permeabilities produces unsatisfactory results, mostly due to nonphysical wave transmissions into obstacles, resulting in considerable energy and mass losses in reflected waves. The objective of this paper is to extend the Brinkman penalization technique to compressible flows based on a physically sound mathematical model for compressible flows through porous media. In addition to penalizing momentum and energy equations, the continuity equation for porous media is considered inside obstacles. In this model, the penalized porous region acts as a high impedance medium, resulting in negligible wave transmissions. The asymptotic analysis reveals that the proposed Brinkman penalization technique results in the amplitude and phase errors of order O((ηϕ)^1/2) and O((η/η_T)^1/4ϕ^3/4), when the boundary layer within the porous media is respectively resolved or unresolved. The proposed method is tested using 1- and 2-D benchmark problems. The results of direct numerical simulation are in excellent agreement with the analytical solutions. The numerical simulations verify the accuracy and convergence rates.     

A Brinkman penalization method for compressible flows in complex geometries

To simulate flows around solid obstacles of complex geometries, various immersed boundary methods had been developed. Their main advantage is the efficient implementation for stationary or moving solid boundaries of arbitrary complexity on fixed non-body conformal Cartesian grids. The Brinkman penalization method was proposed for incompressible viscous flows by penalizing the momentum equations. Its main idea is to model solid obstacles as porous media with porosity, , and viscous permeability approaching zero. It has the pronounced advantages of mathematical proof of error bound, strong convergence, and ease of numerical implementation with the volume penalization technique. In this paper, it is extended to compressible flows. The straightforward extension of penalizing momentum and energy equations using Brinkman penalization with respective normalized viscous, η, and thermal, η_T, permeabilities produces unsatisfactory results, mostly due to nonphysical wave transmissions into obstacles, resulting in considerable energy and mass losses in reflected waves. The objective of this paper is to extend the Brinkman penalization technique to compressible flows based on a physically sound mathematical model for compressible flows through porous media. In addition to penalizing momentum and energy equations, the continuity equation for porous media is considered inside obstacles. In this model, the penalized porous region acts as a high impedance medium, resulting in negligible wave transmissions. The asymptotic analysis reveals that the proposed Brinkman penalization technique results in the amplitude and phase errors of order O((η)^1/2) and O((η/η_T)^1/43/4), when the boundary layer within the porous media is respectively resolved or unresolved. The proposed method is tested using 1- and 2-D benchmark problems. The results of direct numerical simulation are in excellent agreement with the analytical solutions. The numerical simulations verify the accuracy and convergence rates.     

Flow visualization and density measurement of compressible flows by a laser speckle method

Laser speckle method is a well known technique that is useful for both visualization and quantitative measurement. This technique was applied to the density measurement of Mach reflection of shock waves in the present experiment. The Object of the measurement is the density field of simple Mach reflection in relatively low shock Mach number. The non-uniform flow field is divided into three regions by incident, Mach and reflected shock waves. A shock tube was employed in the present experiment. Wedges of 20 degrees and 45 degrees were placed in the test section. YAG laser was employed as a light source. Speckle photograph was taken by a digital still camera. Simple subtraction between the reference and flow images shows a shock pattern and a degree of the correlation of speckle pattern in the flow field. Thus, we can obtain a visualized flow image showing a configuration of Mach reflection from speckle photograph. Speckle photographs which was obtained in the experiments were processed with cross-correlation method. A reconstructed density gradient vector map of Mach reflection was obtained. Comparing the experimental result with numerical one, the measured density gradient shows a good agreement with theoretical prediction.     

Assessment of measuring errors in DIC using deformation fields generated by plastic FEA

In this article, systematic errors that arise from different implementations of digital image correlation (DIC) techniques are analyzed. In particular, we investigate the influence of the adopted correlation function, the interpolation order, the shape function and the subset size on the derived displacements. These errors are estimated using numerically deformed images that were obtained by imposing finite element (FE) displacement fields on an undeformed image yielding plastic deformation of the specimen. This FE procedure simulates realistic experimental heterogeneous deformations at various load steps. It is shown that DIC is able to reproduce these displacements up to a satisfactory level if conscious choices in the above-mentioned implementations are made.     

A front-tracking/ghost-fluid method for fluid interfaces in compressible flows

A front-tracking/ghost-fluid method is introduced for simulations of fluid interfaces in compressible flows. The new method captures fluid interfaces using explicit front-tracking and defines interface conditions with the ghost-fluid method. Several examples of multiphase flow simulations, including a shock–bubble interaction, the Richtmyer–Meshkov instability, the Rayleigh–Taylor instability, the collapse of an air bubble in water and the breakup of a water drop in air, using the Euler or the Navier–Stokes equations, are performed in order to demonstrate the accuracy and capability of the new method. The computational results are compared with experiments and earlier computational studies. The results show that the new method can simulate interface dynamics accurately, including the effect of surface tension. Results for compressible gas–water systems show that the new method can be used for simulations of fluid interface with large density differences.     

Direct numerical simulation of passive scalar in decaying compressible turbulence

1IntroductionDirectnumericalsimulation(DNS)becomesanimportanttoolinrecentresearchofturbulence[1].DNSofcompressibleturbulenceismoredifficultthanthatoftheincompressibleturbulence.WhentheturbulentMachnumberisgreaterthan0.3theshockletsmayappearinthecompressibleturbulentflowfields.Thereasonandmechanismofshockletsexistencearenotclearyet.TheturbulentMachnumberinDNScannotbeveryhighwiththepresentexistingnumericalmethodsandcomputerresource.Fortheproblemofcompressibleisotropicturbulencewiththeinitia…     

Numerical investigation of cavitation generated by an ultrasonic dental scaler tip vibrating in a compressible liquid

Bacterial biofilm accumulation around dental implants is a significant problem leading to peri-implant diseases and implant failure. Cavitation occurring in the cooling water around ultrasonic scaler tips can be used as a novel solution to remove debris without any surface damage. However, current clinically available instruments provide insufficient cavitation around the activated tip surface. To solve this problem a critical understanding of the vibro-acoustic behaviour of the scaler tip and the associated cavitation dynamics is necessary. In this research, we carried out a numerical study for an ultrasound dental scaler with a curved shape tip vibrating in water, using ABAQUS based on the finite element method. We simulated the three-dimensional, nonlinear and transient interaction between the vibration and deformation of the scaler tip, the water flow around the scaler and the cavitation formation and dynamics. The numerical model was well validated with the experiments and there was excellent agreement for displacement at the free end of the scaler. A systematic parametric study has been carried out for the cavitation volume around the scaler tip in terms of the frequency, amplitude and power of the tip vibration. The numerical results indicate that the amount of cavitation around the scaler tip increases with the frequency and amplitude of the vibration. However, if the frequency is far from the natural frequency, the cavitation volume around the free end decreases due to reduced free end vibration amplitude.     

静电放电火花产生的电磁场树枝模型

 静电放电火花产生的电磁脉冲会对电子系统的正常工作造成严重的干扰，甚至造成系统的损伤。利用时域有限差分法建立了静电放电火花产生的电磁场的数值模型，模型中充分考虑了放电电极上的静电荷对电场的影响。把由此模型放电计算的电磁场值与由此解析方法得到的场值进行了比较，结果吻合良好。由此可以用此模型来研究静电放电火花产生的电磁场与电子系统的能量耦合问题。     

静电放电火花产生的电磁场树枝模型

静电放电火花产生的电磁脉冲会对电子系统的正常工作造成严重的干扰,甚至造成系统的损伤。利用时域有限差分法建立了静电放电火花产生的电磁场的数值模型,模型中充分考虑了放电电极上的静电荷对电场的影响。把由此模型放电计算的电磁场值与由此解析方法得到的场值进行了比较,结果吻合良好。由此可以用此模型来研究静电放电火花产生的电磁场与电子系统的能量耦合问题。     

Out-of-plane bipolar and quadrupolar magnetic fields generated by shear flows in two-dimensional resistive reconnection

Shear flows perpendicular to the anti-parallel reconnecting magnetic field are often observed in magnetosphere and interplanetary plasmas, and in laboratory plasmas toroidal differential rotations can also be generated in magnetic confinement devices. Our study finds that such shear flows can generate bipolar or quadrupolar out-of-plane magnetic field perturbations in a two-dimensional resistive MHD reconnection without the Hall effects. The quadrupolar structure has otherwise been thought a typical Hall MHD reconnection feature caused by the in-plane electron convection. The results will challenge the conventional understanding and satellite observations of the signature of reconnection evidences in space plasmas.