Relaxation system based sub-grid scale modelling for large eddy simulation of Burgers' equation

An implicit sub-grid scale model for large eddy simulation is presented by utilising the concept of a relaxation system for one dimensional Burgers' equation in a novel way. The Burgers' equation is solved for three different unsteady flow situations by varying the ratio of relaxation parameter (epsilon) to time step. The coarse mesh results obtained with a relaxation scheme are compared with the filtered DNS solution of the same problem on a fine mesh using a fourth-order CWENO discretisation in space and third-order TVD Runge-Kutta discretisation in time. The numerical solutions obtained through the relaxation system have the same order of accuracy in space and time and they closely match with the filtered DNS solutions.

A Numerical Scheme for Three-Dimensional Front Propagation and Control of Jordan Mode

As an example of a front propagation, we study the propagation of a three-dimensional nonlinear wavefront into a polytropic gas in a uniform state and at rest. The successive positions and geometry of the wavefront are obtained by solving the conservation form of equations of a weakly nonlinear ray theory. The proposed set of equations forms a weakly hyperbolic system of seven conservation laws with an additional vector constraint, each of whose components is a divergence-free condition. This constraint is an involution for the system of conservation laws, and it is termed a geometric solenoidal constraint. The analysis of a Cauchy problem for the linearized system shows that when this constraint is satisfied initially, the solution does not exhibit any Jordan mode. For the numerical simulation of the conservation laws we employ a high resolution central scheme. The second order accuracy of the scheme is achieved by using MUSCL-type reconstructions and Runge-Kutta time discretizations. A constrained transport-type technique is used to enforce the geometric solenoidal constraint. The results of several numerical experiments are presented, which confirm the efficiency and robustness of the proposed numerical method and the control of the Jordan mode.

SuSeFLAV 1.2: Program for supersymmetric mass spectra with seesaw mechanism and rare lepton flavor violating decays

Accurate supersymmetric spectra are required to confront data from direct and indirect searches of supersymmetry. SuSeFLAV is a numerical tool capable of computing supersymmetric spectra precisely for various supersymmetric breaking scenarios applicable even in the presence of flavor violation. The program solves MSSM RGEs with complete 3 x 3 flavor mixing at 2-loop level and one loop finite threshold corrections to all MSSM parameters by incorporating radiative electroweak symmetry breaking conditions. The program also incorporates the Type-I seesaw mechanism with three massive right handed neutrinos at user defined mass scales and mixing. It also computes branching ratios of flavor violating processes such as l(j) -> l(i)gamma, l(j) -> 3 l(i), b -> s gamma and supersymmetric contributions to flavor conserving quantities such as (g(mu) - 2). A large choice of executables suitable for various operations of the program are provided. Program summary Program title: SuSeFLAV Catalogue identifier: AEOD_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEOD_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GNU General Public License No. of lines in distributed program, including test data, etc.: 76552 No. of bytes in distributed program, including test data, etc.: 582787 Distribution format: tar.gz Programming language: Fortran 95. Computer: Personal Computer, Work-Station. Operating system: Linux, Unix. Classification: 11.6. Nature of problem: Determination of masses and mixing of supersymmetric particles within the context of MSSM with conserved R-parity with and without the presence of Type-I seesaw. Inter-generational mixing is considered while calculating the mass spectrum. Supersymmetry breaking parameters are taken as inputs at a high scale specified by the mechanism of supersymmetry breaking. RG equations including full inter-generational mixing are then used to evolve these parameters up to the electroweak breaking scale. The low energy supersymmetric spectrum is calculated at the scale where successful radiative electroweak symmetry breaking occurs. At weak scale standard model fermion masses, gauge couplings are determined including the supersymmetric radiative corrections. Once the spectrum is computed, the program proceeds to various lepton flavor violating observables (e.g., BR(mu -> e gamma), BR(tau -> mu gamma) etc.) at the weak scale. Solution method: Two loop RGEs with full 3 x 3 flavor mixing for all supersymmetry breaking parameters are used to compute the low energy supersymmetric mass spectrum. An adaptive step size Runge-Kutta method is used to solve the RGEs numerically between the high scale and the electroweak breaking scale. Iterative procedure is employed to get the consistent radiative electroweak symmetry breaking condition. The masses of the supersymmetric particles are computed at 1-loop order. The third generation SM particles and the gauge couplings are evaluated at the 1-loop order including supersymmetric corrections. A further iteration of the full program is employed such that the SM masses and couplings are consistent with the supersymmetric particle spectrum. Additional comments: Several executables are presented for the user. Running time: 0.2 s on a Intel(R) Core(TM) i5 CPU 650 with 3.20 GHz. (c) 2012 Elsevier B.V. All rights reserved.

Nonrotating Beams Isospectral to Tapered Rotating Beams

In this paper, we seek to find nonrotating beams that are isospectral to a given tapered rotating beam. Isospectral structures have identical natural frequencies. We assume the mass and stiffness distributions of the tapered rotating beam to be polynomial functions of span. Such polynomial variations of mass and stiffness are typical of helicopter and wind turbine blades. We use the Barcilon-Gottlieb transformation to convert the fourth-order governing equations of the rotating and the nonrotating beams, from the (x, Y) frame of reference to a hypothetical (z, U) frame of reference. If the coefficients of both the equations in the (z, U) frame match with each other, then the nonrotating beam is isospectral to the given rotating beam. The conditions on matching the coefficients lead to a pair of coupled differential equations. Wesolve these coupled differential equations numerically using the fourth-order Runge-Kutta scheme. We also verify that the frequencies (given in the literature) of standard tapered rotating beams are the frequencies (obtained using the finite-element analysis) of the isospectral nonrotating beams. Finally, we present an example of beams having a rectangular cross-section to show the application of our analysis. Since experimental determination of rotating beam frequencies is a difficult task, experiments can be easily conducted on these isospectral nonrotating beams to calculate the frequencies of the rotating beam.

LEVEL SET输运方程的求解方法及其对气-液两相流运动界面数值模拟的影响

用数值模拟方法来研究气-液两相流动与传热现象是当今多相流领域的一个热门课题.由于两相流固有的复杂性,气-液两相流界面迁移现象的数值模拟一直是两相流研究中的一大难点.本文介绍了捕捉气-液两相流相界面运动的水平集方法(Level Set)及其研究进展,介绍了求解Level Set输运方程的3种方法,即一般差分格式、Superbee-TVD格式和Runge-Kutta法-5阶WENO组合格式.结合主流场的求解,分别用这3种方法对4种典型相界面在5种流场中的迁移特性进行了模拟计算,并对计算结果进行了比较和分析.结果表明,Runge-Kutta法-5阶WENO组合格式求解Level Set输运方程的效果最好,在以后的计算中将主要采用这种组合格式来进行气-液相界面输运方程的求解.

Chaotic dynamic analysis of viscoelastic plates

The dynamic buckling of viscoelastic plates with large deflection is investigated in this paper by using chaotic and fractal theory. The material behavior is given in terms of the Boltzmann superposition principle. in order to obtain accurate computation results, the nonlinear integro-differential dynamic equation is changed into an autonomic four-dimensional dynamical system. The numerical time integrations of equations are performed by using the fourth-order Runge-Kutta method. And the Lyapunov exponent spectrum, the fractal dimension of strange attractors and the time evolution of deflection are obtained. The influence of geometry nonlinearity and viscoelastic parameter on the dynamic buckling of viscoelastic plates is discussed.

Numerical solution of the incompressible Navier-Stokes equations with an upwind compact difference scheme

A new finite difference method for the discretization of the incompressible Navier-Stokes equations is presented. The scheme is constructed on a staggered-mesh grid system. The convection terms are discretized with a fifth-order-accurate upwind compact difference approximation, the viscous terms are discretized with a sixth-order symmetrical compact difference approximation, the continuity equation and the pressure gradient in the momentum equations are discretized with a fourth-order difference approximation on a cell-centered mesh. Time advancement uses a three-stage Runge-Kutta method. The Poisson equation for computing the pressure is solved with preconditioning. Accuracy analysis shows that the new method has high resolving efficiency. Validation of the method by computation of Taylor's vortex array is presented.

Analysis of super compact finite difference method and application to simulation of vortex-shock interaction

Turbulence and aeroacoustic noise high-order accurate schemes are required, and preferred, for solving complex flow fields with multi-scale structures. In this paper a super compact finite difference method (SCFDM) is presented, the accuracy is analysed and the method is compared with a sixth-order traditional and compact finite difference approximation. The comparison shows that the sixth-order accurate super compact method has higher resolving efficiency. The sixth-order super compact method, with a three-stage Runge-Kutta method for approximation of the compressible Navier-Stokes equations, is used to solve the complex flow structures induced by vortex-shock interactions. The basic nature of the near-field sound generated by interaction is studied.

Surface waves in a vertically excited circular cylindrical container

The nonlinear free surface amplitude equation, which has been derived from the inviscid fluid by solving the potential equation of water waves with a singular perturbation theory in a vertically oscillating rigid circular cylinder, is investigated successively in the fourth-order Runge-Kutta approach with an equivalent time-step. Computational results include the evolution of the amplitude with time, the characteristics of phase plane determined by the real and imaginary parts of the amplitude, the single-mode selection rules of the surface waves in different forced frequencies, contours of free surface displacement and corresponding three-dimensional evolution of surface waves, etc. In addition, the comparison of the surface wave modes is made between theoretical calculations and experimental measurements, and the results are reasonable although there are some differences in the forced frequency.

变截面流化床理论模型及数值模拟

基于双流体模型,在一定假设的条件下推导得到了用于描述变截面床流化过程的局部平衡模型。对于局部平衡模型,空间导数项的离散采用五阶精度的WENO有限差分格式,时间导数项的离散采用TVD Runge-Kutta型的离散格式,对流量突变后的瞬态过程进行了模拟,得到固相体积分数在整个变化过程中沿床高的分布以及床高变化规律和床层表面颗粒速度变化曲线。对于流量突增过程,在床内不同位置形成了一系列的连续波,当所有的连续波到达床面整个变化过程结束。而对于流量突然减小过程,将会有固相体积分数间断在分布板处形成,当所有间断到达床面时,塌落过程结束。

可压缩尖锥边界层湍流的直接数值模拟

采用高精度差分方法对来流马赫数0.7,来流Reynolds数250000/Inch,锥角为20°的尖锥边界层的整个空间转捩过程进行了直接数值模拟.对流项采用了7阶迎风格式离散,黏性项采用6阶中心格式离散,时间推进为3阶Runge-Kutta方法.对转捩形成的充分发展湍流进行了统计分析,包括平均速度分布,近壁湍流强度和雷诺应力等统计数据与平板边界层理论及实验吻合很好,验证了结果的正确性.显示了近壁湍流的典型拟序结构——高、低速条带结构并根据流向速度的周向相关量确定了条带的间距,以当地壁面尺度度量的条带间距沿流向并没有显著变化.给出了柱坐标下的可压湍动能发展方程,并据此对近壁湍动能的生成、耗散和输运机制进行了分析.

微薄梁三点弯曲的尺度效应研究

应用高灵敏度的力传感器以及时间序列电子散斑干涉法,同时测出了不同厚度纯镍薄片三点弯曲试件的抗力与变形,得到薄梁中心点处的载荷与挠度曲线.应用Fleck和Hutchinson的偶应力理论,结合平面应变弯曲模型,建立了薄梁处于弹性状态和弹塑性状态的控制方程,应用Runge-Kutta法进行数值求解,并将计算得到的载荷-挠度曲线以及无量纲化弯矩-表面应变曲线和实验结果进行了比较.在理论计算过程中,没有拟合任何材料参数,所有的材料参数均来自实验测量的结果,材料特征尺度也是根据Stolken和Evans的工作给出的.结果表明:应用偶应力理论预测的结果和实验结果符合良好,而经典理论的预测结果与实验不相符合.

蒸汽-冷流体接触冷凝流动的数值模拟

介绍了关于蒸汽-冷流体直接接触冷凝流动与传热的数值计算模型与部分研究结果。用Level Set方法确定蒸汽-冷流体接触界面的位置和形状,建立了对蒸汽和冷流体普遍适用的动量、能量和质量守恒方程,在能量和质量寺恒方程中增加了部分项用于计算蒸汽冷凝所产生的影响。用有限差分法在交错网格上离散控制方程,用Runge-Kutta法-五阶WENO组合格式求解Level Set输运方程,用压力修正的迭代Projection方法求解动量方程,而用SIMPLE方法求解温度控制方程。对算例的计算结果表明,本文所建立的数值计算模型能反映物理现象的宏观特性。根据计算结果,分析了本文模型的优缺点,并指出了今后改进的方向。

不同重力条件下不相溶流体层间界面波演化规律的数值模拟研究

报道了关于不相溶流体层间界面波演化规律的数值模拟研究及结果，重点考察了重力条件对界面波演化特性的影响。考虑在深度方向无限扩展的互不相容的两个流体层，上层流体比下层的轻，但比下层的运动速度快；两层流体间的界面上存在正弦波形的初始扰动，并随流体流动而不断变化。本文采用Level Set方法来实现对运动的相界面的追踪，用有限差分法来离散控制方程组。为了提高数值算法的稳定性，采用三阶的Runge-Kutta法来离散时间导数，而采用五阶的WENO（Weighted Essentially Non-oscillatory）格式来离散一阶对流输运项，并用压力修正投影法（Pressure Correction Projection Method）来实现离散控制方程组的求解。为了提高对复杂非稳态过程的解的准确度，采用了嵌套的三层迭代循环。本文对一系列工况条件下的界面波演化过程进行了计算；除了研究重力的作用之外，还考察了流体密度、粘性、表面张力、初始界面波频率、振幅及波数对界面波演化特性的影响。其中，上下流体层的最大密度比和粘性比可达3000/1，而重力加速度在0~5g0（g0=9.8m/s^2）之间变化，上下流体层间的最大速度差为8m/s。研究结果表明，随着重力、流体密度比、流体粘性比及表面张力的增加，界面波的演化受到不同程度的抑制，而界面波的传播速度也与重力及流体的密度、粘性和表面张力等因素相关。