流固耦合力学的数值研究方法的发展及软件应用概述

流固耦合力学是多种工程领域研究的重要课题。本文对流固耦合力学数值研究方法发展过程进行了概述,阐述了各发展阶段使用的三类数值研究方法,并简单介绍了常见的计算流固耦合商业软件,最后对流固耦合力学的发展趋势进行了展望。

许昌学院“分数阶偏微分方程高效数值算法”研究团队科研进展

分数阶微分方程能更精确地描述记忆和遗传现象或过程,已被广泛应用于生物化学、系统控制、金融等诸多科学工程领域。由于分数阶微分方程的精确解不易求得,所以关于该方程的数值方法研究备受关注。许昌学院数理学院“分数阶偏微分方程高效数值算法”研究团队在2019年获批国家自然科学基金面上项目,在研究时间分数阶偏微分方程的高精度有限元逼近方面取得新进展。

地表及岩体移动研究进展

从国内外理论与试验研究的角度总结了地表及岩体移动的预测方法。地表及岩体移动的预测方法可分为唯象法、力学方法、数值方法和物理模拟法4类;讨论了地表及岩体移动研究领域的一些热点问题,并从工程应用、数值方法及本构理论等方面探讨了其发展方向。

MHC耦合作用下岩体力学性质研究现状

随着工程岩体开挖的不断加深,应力(M)-水流(H)-化学(C)耦合作用对岩体的影响越来越大,研究MHC耦合作用下岩体的力学性质具有重要的工程意义。目前,国内外学者对MHC耦合下岩体力学性质研究主要包括试验研究、模型研究和数值方法研究3个方面。其中,MHC耦合过程的试验研究方面主要集中在MH、MC两种因素的耦合研究上,而对于M、H、C 3种因素的耦合研究较少;模型研究方面还需要考虑化学反应过程的引入,以及如何用更精准的数学模型表达岩体MHC耦合中的水化学损伤过程;数值方法研究方面还需要在计算参数取值、本构关系选取和化学反应过程与应力水流的融合上进一步完善。

用复合插值提高离散导热热流场计算精度

提出了一种将中心差分与拉格朗日抛物线插值相结合的复合插值法来确定离散导热热流场，并与其它方法进行比较，结果表明该方法具有更高的精度．

AR5395活动区中的旋转储能(摘要)

本文根据太阳活动区AR5395中存在旋转运动特征,建立了一个扭转的共生磁流管模型,从一组完整的磁流体动力学方程组研究活动区底部的旋转运动储存能量。通过数值方法研究了三维流场和磁场的耦合效应。

Numerical simulation of steady flow past a liquid sphere immersed in simple shear flow at low and moderate Re

This work presents a numerical investigation on steady internal, external and surface flows of a liquid sphere immersed in a simple shear flow at low and intermediate Reynolds numbers. The control volume formulation is adopted to solve the governing equations of two-phase flow in a 3-D spherical coordinate system. Numerical results show that the streamlines for Re = 0 are closed Jeffery orbits on the surface of the liquid sphere, and also closed curves outside and inside the liquid sphere. However, the streamlines have intricate and non-closed structures for Re ≠ 0. The flow structure is dependent on the values of Reynolds number and interior-to-exterior viscosity ratio.

Numerical Investigations on Transient Behaviours of Two 3-D Freely Floating Structures by Using a Fully Nonlinear Method

Two floating structures in close proximity are very commonly seen in offshore engineering. They are often subjected to steep waves and, therefore, the transient effects on their hydrodynamic features are of great concem. This paper uses the quasi arbitrary Lagrangian Eulerian finite element method （QALE-FEM）, based on the fully nonlinear potential theory （FNPT）, to numerically investigate the interaction between two 3-D floating structures, which undergo motions with 6 degrees of freedom （DOFs）, and are subjected to waves with different incident angles. The transient behaviours of floating structures, the effect of the accompanied structures, and the nonlinearity on the motion of and the wave loads on the structures are the main focuses of the study. The investigation reveals an important transient effects causing considerably larger structure motion than that in steady state. The results also indicate that the accompanied structure in close proximity enhances the interaction between different motion modes and results in stronger nonlinearity causing 2hal-order component to be of similar significance to the fundamental one.

New Periodic Wave Solutions and Their Interaction for （2＋1）-dimensional KdV Equation

A class of new doubly periodic wave solutions for （2＋1）-dimensional KdV equation are obtained by introducing appropriate Jacobi elliptic functions and Weierstrass elliptic functions in the general solution（contains two arbitrary functions） got by means of multilinear variable separation approach for （2＋1）-dimensional KdV equation. Limiting cases are considered and some localized excitations are derived, such as dromion, multidromions, dromion-antidromion, multidromions-antidromions, and so on. Some solutions of the dromion-antidromion and multidromions-antidromions are periodic in one direction but localized in the other direction. The interaction properties of these solutions, which are numerically studied, reveal that some of them are nonelastic and some are completely elastic. Furthermore, these results are visualized.

宽色域视频标准ITU-R BT.1361与IEC 61966-2-4的分析和比较

ITU-RBT.1361和IEC61966-2-4(xvYCC)是关于宽色域视频系统的2个标准,它们扩展色域的方法不尽相同。以电视系统中广泛使用的彩条测试信号为例,分析了两者间的差异、计算了实际可用信号幅度范围并比较了它们的色域大小。结果表明,在与常规色域兼容传输情况下,2种标准虽均可扩展色域,但数值不够大,进一步研究扩展色域方法的空间仍较大。

Numerical research on mixing characteristics of different injection schemes for supersonic transverse jet

A numerical method using AUSMDV scheme and k-ω SST turbulence model with an explicit compressibility correction was developed,and a 3-D numerical simulation of a supersonic flow field with a vertical sonic jet of hydrogen was performed.Good agreement between numerical results and experimental data validated the reliability of the numerical method.Whereafter,two parameters,mass-weighted average total pressure and mixing efficiency,were defined to evaluate the mixing performance of different injection schemes.Based on the numerical method and evaluation criterion,the mixing characteristics of different injection schemes were studied in detail.It was found that for the mixing field of supersonic transverse jet,the near-field mixing is controlled by convection transport while the far-field mixing is controlled by mass diffusion;the circular-hole injection causes a loss of total pressure comparable to the slot injection,but can induce a much higher mixing efficiency because of its 3-D flow characteristic;the variation of injection angle under circular-hole injection mainly affects the near-field mixing degree,and among the five injection angles studied in the present paper,angle 120° is the optimal one;with the increase of the ratio between injector space and diameter,the induced mixing efficiency increases while the caused loss of total pressure can grow greatly;the two-stage injection method designed through reducing the injector area to keep the same hydrogen mass flowrate can induce a much higher mixing efficiency while only a bit larger loss of total pressure when compared to the single-stage injection,and hence the two-stage injection is superior to the single-stage injection.The research results can direct the design of the fuel injection method in the combustor of scramjet engine.

Numerical Investigation of Active Tip-clearance Control through Tip Cooling Injection in an Axial Turbine Cascade

This paper presents a numerical investigation of an active tip-clearance control method based on cooling injectionfrom the blade tip surface. It aims to study the influences of air injection on controlling tip clearance flow, withemphasis on the effects of the injection location on secondary flow and the potential thermal benefits from thecooling jet. The results show that injection location plays an important role in the redistribution of secondary flowwithin the cascade passage. Injection located much closer to the pressure-side comer performs better in reducingtip clearance massflow and its associated losses. However, it also intensifies tip passage vortex, due to less restraintderiving from the reduced tip clearance vortex. Lower plenum total pressure is required to inject equivalentamount of cooling air, but the heat transfer condition on the blade tip surface is a bit worse than that with injectionfrom the reattachment region. Thus the optimum location of air injection should be at the tip separation vortex region.

Numerical study on two-phase flow through fractured porous media

Based on the flux equivalent principle of a single fracture, the discrete fracture concept was developed, in which the macroscopic fractures are explicitly described as （n-l） dimensional geometry element. On the fundamental of this simplification, the discrete-fractured model was developed which is suitable for all types of fractured porous media. The principle of discrete-fractured model was introduced in detail, and the general mathematical model was expressed subsequently. The fully coupling discrete-fractured mathematical model was implemented using Galerkin finite element method. The validity and accuracy of the model were shown through the Buckley-Leverett problem in a single fracture. Then the discrete-fractured model was applied to the two different type fractured porous media. The effect of fractures on the water flooding in fractured reservoirs was investigated. The numerical results showed that the fractures made the porous media more heterogeneous and anisotropic, and that the orientation, size, type of fracture and the connectivity of fractures network have important impacts on the two-phase flow.

Numerical simulation of immiscible liquid-liquid flow in microchannels using lattice Boltzmann method

Immiscible kerosene-water two-phase flows in microchannels connected by a T-junction were numerically studied by a Lattice Boltzmann (LB) method based on field mediators.The two-phase flow lattice Boltzmann model was first validated and improved by several test cases of a still droplet.The five distinct flow regimes of the kerosene-water system,previously identified in the experiments from Zhao et al.,were reproduced.The quantitative and qualitative agreement between the simulations and the experimental data show the effectiveness of the numerical method.The roles of the interfacial tension and contact angle on the flow patterns and shapes of droplets were discussed and highlighted according to the numerical results based on the improved two-phase LB model.This work demonstrated that the developed LBM simulator is a viable tool to study immiscible two-phase flows in microchannels,and such a tool could provide tangible guidance for the design of various microfluidic devices that involve immiscible multi-phase flows.

Continuous Method Development and Numerical Study of HHV Water Gas Production by Pulverized Coal

In energy industries, it is always of difficulty to produce high heat value（HHV） gas continuously using pulverized coal. In this paper, a new type furnace for partitioned alternative gasification using pulverized coal is developed, in which the oxidation and reduction reaction occur alternatively with the dropping of pulverized coal and finally HHV gas could be continuously obtained at the reduction zone exit and low heat value（LHV） gas at the oxidation zone exit. Furthermore, the gasification characteristics and their factors in furnace are numerically simulated under two dimensional model with a self-coded program, based on heat, mass and energy transfer as well as reaction principles. It is found that the producing rate of HHV gas is 1.10Nm3/kg with heat value of ll.72MJ/Nm3, how- ever, that of LHV gas is 2.58Nm3/kg with heat value of 5.30MJ/Nm3, and the coal gas efficiency is 81.16% under optimized conditions.

Numerical study on flexural behaviors of steel reinforced engineered cementitious composite(ECC) and ECC/concrete composite beams

Engineered cementitious composite(ECC)is a class of high performance cementitious composites with pseudo strain-hardening behavior and excellent crack control capacity.Substitution of concrete with ECC can largely reduce the cracking and durability problems associated with brittleness of concrete.In this paper,a simplified constitutive model of the ECC material was applied to simulate the flexural behaviors of the steel reinforced ECC and ECC/concrete composite beams with finite element method.The simulation results are found to be in good agreement with test results,indicating that the finite element model is reasonably accurate in simulating the flexural behaviors of the steel reinforced ECC flexural members.The effects of the ECC modulus,ECC tensile ductility,ECC thickness and ECC position on flexural behaviors in terms of ultimate moment,deflection and the maximum crack width of the steel reinforced ECC or ECC/concrete composite beam are hence evaluated.

A mesoscale model based on Monte-Carlo method for concrete fracture behavior study

An aggregate generation and packing algorithm based on Monte-Carlo method is developed to express the aggregate random distribution in cement concrete. A mesoscale model is proposed on the basis of the algorithm. In this model, the concrete con- sists of three parts, namely coarse aggregate, cement matrix and the interracial transition zone （ITZ） between them. To verify the proposed model, a three-point bending beam test was performed and a series of two-dimensional mesoscale concrete mod- els were generated for crack behavior investigation. The results indicate that the numerical model proposed in this study is helpful in modeling crack behavior of concrete, and that treating concrete as heterogeneous material is very important in frac- ture modeling.

Numerical Study on Spectral Domain Optical Coherence Tomography Spectral Calibration and Re-Sampling Importance

A spectral calibration technique, a data processing method and the importance of calibration and re-sampling methods for the spectral domain optical coherence tomography system were numerically studied, targeted to optical coherence tomography （OCT） signal processing implementation under graphics processing unit （GPU） architecture. Accurately, assigning the wavelength to each pixel of the detector is of paramount importance to obtain high quality images and increase signal to noise ratio （SNR）. High quality imaging can be achieved by proper calibration methods, here performed by phase calibration and interpolation. SNR was assessed employing two approaches, single spectrum moving window averaging and consecutive spectra data averaging, to investigate the optimized method and factor for background noise reduction. It was demonstrated that the consecutive spectra averaging had better SNR performance.

Analytical and Numerical Studies of Quantum Plateau State in One Alternating Heisenberg Chain

By using the coupled duster method and the numerical density matrix renormalization group method, we investigate the properties of the quantum plateau state in an alternating Heisenberg spin chain. In the absence of a magnetic field, the results obtained from the coupled cluster method and density matrix renormalization group method both show that the ground state of the aiternating chain is a gapped dimerized state when the parameter a exceeds a critical point ac. The value of the critical points can be determined precisely by a detailed investigation of the behavior of the spin gap. The system therefore possesses an m = 0 plateau state in the presence of a magnetic field When a 〉 ac. In addition to the m = 0 plateau state, the results of density matrix renormaiization group indicate that there is an m = 1/4 plateau state that occurs between two critical fields in the alternating chain if a 〉 1. The mechanism for the m = 1/4 plateau state and the critical behavior of the magnetization as one approaches this plateau state are also discussed.

A numerical study on matching relationships of gravity waves in nonlinear interactions

Applying a fully nonlinear numerical scheme with second-order temporal and spatial precision,nonlinear interactions of gravity waves are simulated and the matching relationships of the wavelengths and frequencies of the interacting waves are discussed.In resonant interactions,the wavelengths of the excited wave are in good agreement with the values derived from sum or difference resonant conditions,and the frequencies of the three waves also satisfy the matching condition.Since the interacting waves obey the resonant conditions,resonant interactions have a reversible feature that for a resonant wave triad,any two waves are selected to be the initial perturbations,and the third wave can then be excited through sum or difference resonant interaction.The numerical results for nonresonant triads show that in nonresonant interactions,the wave vectors tend to approximately match in a single direction,generally in the horizontal direction.The frequency of the excited wave is close to the matching value,and the degree of mismatching of frequencies may depend on the combined effect of both the wavenumber and frequency mismatches that should benefit energy exchange to the greatest extent.The matching and mismatching relationships in nonresonant interactions differ from the results of weak interaction theory that the wave vectors are required to satisfy the resonant matching condition but the frequencies are permitted to mismatch and oscillate with amplitude of half the mismatching frequency.Nonresonant excitation has an irreversible characteristic,which is different from what is found for the resonant interaction.For specified initial primary and secondary waves,it is difficult to predict the values of the mismatching wavenumber and frequency for the excited wave owing to the complexity.