Lattice Boltzmann Simulation of Convection in a Porous Medium with Temperature Jump and Velocity Slip Boundary Conditions

To investigate the convection in a porous medium, a horizontal quiescent layer of one fluid saturating a porous medium heated from bottom is numerically studied using single lattice-Boltzmann method (LBM) and the generalized Navier-Stokes equation proposed by Nithiarasu et al. [P. Nithiarasu, K.M. Seetharamu, and T Sundararajan, Int. J. Heat Mass Trans. 40 (1997) 3955]. Due to the rarefaction, the boundary conditions are considered as both temperature jump and velocity slip. The computational results are validated against the analytical results, and excellent agreement has been obtained. The results have shown that the Rayleigh number is increased with increasing temperature jump, the stabilization effect of temperature is much more significant than that of velocity slip, and the computation stability of present model is better than that of Darcy and Brinkman models.     

缓解企业创新融资约束的途径：市场还是政府?

本文基于金融市场和政府公共政策影响企业创新活动资金来源的理论机制,对欧拉方程加以拓展构建研发投资的融资约束模型,以2010—2019年制造业上市公司为样本,就企业技术创新的融资约束问题及其缓解途径进行了实证检验。结果表明：制造业上市公司的创新活动普遍存在着融资约束,市场途径中的股权融资方式能够缓解资金不足问题,银行贷款和债券融资没有发挥相应的作用,政府途径中的直接财政补贴并没有解决资金短缺困境,采用事后返还方式的税收优惠则效果十分显著;从不同所有制来看,国有企业有着较为充足的融资途径,而民营企业仍然面临着较为严重的金融歧视,只有税收优惠能够起到激励效果;地区分组表明,市场化先进地区的企业融资渠道更为丰富,滞后地区的企业难以通过市场途径获得研发投资所需资金,公共政策也没有起到“雪中送炭”的功能;不同技术密集度行业的实证结果显示,高技术企业得到了金融市场和政府层面的共同支持,低技术企业的研发投资面临着严重的融资困境。     

巨灾互换对巨灾风险可保性边界的影响——以超强台风“天鸽”为例

巨灾衍生品对提升保险公司巨灾风险承保能力意义深远。本文探讨了巨灾互换这一新型巨灾衍生产品对巨灾风险可保性边界的影响。构建了利差连续给付型巨灾互换的定价模型,推导了保险公司参与巨灾互换情形下的破产概率模型和可保性边界函数。并以超强台风“天鸽”对我国造成的损失数据为样本,给出了巨灾互换价差的计算实例,讨论了巨灾互换对可保性边界的影响。进一步,分析了损失赔付比例与初始资本、安全负荷系数和破产概率的关系。结果表明,保险公司参与巨灾互换能够有效拓展巨灾风险的可保性边界。研究结论可为巨灾互换在中国的研发与应用提供数据支持和决策依据。     

Preservation of Linear Constraints in Approximation of Tensors

For an arbitrary tensor (multi-index array) with linear constraints at each direction, it is proved that the factors of any minimal canonical tensor approximation to this tensor satisfy the same linear constraints for the corresponding directions.     

Generalized variational principles for boundary value problem of electromagnetic field in electrodynamics

An expression of the generalized principle of virtual work for the boundary value problem of the linear and anisotropic electromagnetic field is given. Using Chien＇s method, a pair of generalized variational principles （GVPs） are established, which directly leads to all four Maxwell＇s equations, two intensity-potential equations, two constitutive equations, and eight boundary conditions. A family of constrained variational principles is derived sequentially. As additional verifications, two degenerated forms are obtained, equivalent to two known variational principles. Two modified GVPs are given to provide the hybrid finite element models for the present problem.     

快速小波边界元的矩阵后压缩方法

介绍了一种基于传统边界元单元划分的小波 Galerkin 边界元法,该方法具有几乎线性(即 ON,N 为自由度)的求解复杂度。在准消失矩小波的框架下介绍了非标准型系数矩阵的压缩问题,提出了一种后压缩算法以降低小波边界元法的内存消耗。求解 Stokes 方程的算例表明,后压缩算法在保证结果收敛特性的情况下可以将系数矩阵的内存占用量降低 5 倍以上。     

一种基于多矩的有限体积浸入边界法

基于多矩VSIAM3格式及浸入边界法,提出一套在复杂计算区域内求解不可压缩流动的数值格式.不可压N-S方程使用VSIAM3格式进行离散,引入浸入边界法处理复杂、移动边界,使用虚拟网格方法计算动量方程修正项,同时还考虑了对连续方程的修正.使用标准算例对数值模式进行验证.     

消声器声学性能预测的子结构快速多极子边界元法

鉴于快速多极子边界元法的应用主要局限于单区域声学问题计算,发展基于子结构技术的快速多极子边界元法以计算多区域声场问题,介绍基本原理、具体实施过程以及优缺点.以带有插进口管的膨胀腔消声器为例,应用子结构快速多极子边界元法和传统边界元法计算其传递损失,通过与实验测量结果的比较,验证方法的有效性和计算精度.研究表明,快速多极子边界元法与传统边界元法相比,节点数越多,其在节省计算时间,减少计算量等方面的优势越明显.     

TRANSFER MATRIX METHOD FOR ANALYZING VIBRATION AND DAMPING CHARACTERISTICS OF ROTATIONAL SHELL WITH PASSIVE CONSTRAINED LAYER DAMPING TREATMENT

The first order differential matrix equations of the host shell and constrained layer for a sandwich rotational shell are derived based on the thin shell theory.Employing the layer wise principle and first order shear deformation theory, only considering the shearing deformation of the viscoelastic layer, the integrated first order differential matrix equation of a passive constrained layer damping rotational shell is established by combining with the normal equilibrium equation of the viscoelastic layer.A highly precise transfer matrix method is developed by extended homogeneous capacity precision integration technology.The numerical results show that present method is accurate and effective.     

APPROXIMATE SOLUTIONS FOR TRANSIENT RESPONSE OF CONSTRAINED DAMPING LAMINATED CANTILEVER PLATE

The series composed by beam mode function is used to approximate the displacement function of constrained damping of laminated cantilever plates, and the transverse deformation of the plate on which a concentrated force is acted is calculated using the principle of virtual work.By solving Lagrange＇s equation, the frequencies and model loss factors of free vibration of the plate are obtained, then the transient response of constrained damping of laminated cantilever plate is obtained, when the concentrated force is withdrawn suddenly.The theoretical calculations are compared with the experimental data, the results show：both the frequencies and the response time of theoretical calculation and its variational law with the parameters of the damping layer are identical with experimental results.Also, the response time of steel cantilever plate, unconstrained damping cantilever plate and constrained damping cantilever plate are brought into comparison, which shows that the constrained damping structure can effectively suppress the vibration.