Flow physics of normal and abnormal bioprosthetic aortic valves

Flow physics of transvalvular flows in the aorta with bioprosthetic valves are investigated using computational modelling. For the efficient simulations of flow-structure-interaction in transvalvular flows, a simplified, reduced degree of freedom valve model is employed with a sharp interface immersed boundary based incompressible flow solver. Simulations are performed for normal as well as abnormal valves with reduced leaflet motion that models the effect of early leaflet thrombosis. The structure of the aortic jet and the hemodynamic stresses on the aortic wall are analysed to understand the hemodynamic impacts and possible long-term clinical implications of sub-clinical, reduced leaflet motion. The simulation results have shown that the reduced leaflet motion tilts the direction of aortic jet and generates stronger flow separation and re-attachment on the aortic wall downstream from the reduced motion leaflets. The modified flow pattern increases the wall pressure fluctuation and average wall shear stress on the downstream aortic wall, and results in the asymmetric oscillatory shear index distributions, which may have long-term clinical implications such as aortic wall damage and remodelling.     

Spectral estimates for saddle point matrices arising in weak constraint four‐dimensional variational data assimilation

We consider the large sparse symmetric linear systems of equations that arise in the solution of weak constraint four‐dimensional variational data assimilation, a method of high interest for numerical weather prediction. These systems can be written as saddle point systems with a 3 × 3 block structure but block eliminations can be performed to reduce them to saddle point systems with a 2 × 2 block structure, or further to symmetric positive definite systems. In this article, we analyse how sensitive the spectra of these matrices are to the number of observations of the underlying dynamical system. We also obtain bounds on the eigenvalues of the matrices. Numerical experiments are used to confirm the theoretical analysis and bounds.     

砂岩方梁断裂过程区范围的实验确定方法

裂纹前端的断裂过程区是引起岩石非线性断裂及尺寸效应的主要原因。利用数字图像相关技术对砂岩开展了三点弯曲梁实验,获得观测区域高精度的全场位移和应变数据,根据断裂韧带区域水平位移和水平应变的分布特征,结合裂尖岩石颗粒变化的微观分析,提出采用裂纹尖端水平位移波动性和水平应变突变性所得到的波动系数和水平应变突变值,确定断裂过程区形状和临界尺寸的方法。结果表明:砂岩断裂过程区的形状为不规则的狭长带状区域,断裂过程区的临界长度为11~13mm,临界宽度为1.58~2.36mm。断裂过程区区域内形变在趋向裂尖时呈指数增加,但其单位区域内的形变增量呈波动状态。该方法能够更加准确判断岩石断裂过程区的范围,有助于分析岩石的非线性断裂特性。     

QEPAS系统中石英音叉模态仿真计算与研究

对石英音叉增强型光声光谱(QEPAS)系统中常用的石英音叉进行了有限元模态计算，获得石英音叉前6阶振型与模态频率，认知了第4阶对称摆动振型为有效振动，利用单因素法分析了石英音叉的音臂长度l₁、音臂宽度w₁、音臂厚度t、音臂切角θ、音臂圆孔直径d及音臂圆孔高度h对低阶有效共振频率(Fre)的影响，敏感度依次为： l₁> w₁>d>θ>t>h，考虑实际设计情形，筛选出了l₁，w₁，d与h四个石英音叉设计变量，采用Box-Behnken实验设计方案与RSM(response surface methodology)方法，以Fre为函数目标，建立l₁，w₁，d与h的二次回归响应面模型，得到了参数之间的交互作用，利用Design-Expert软件对响应面模型进行设计参数反求，结果表明，在15 000 Hz≤Fre≤25 000 Hz计算区域内误差较小，基本满足QEPAS系统的计算需求，所提出的研究与设计方法具有一定通用性，可为QEPAS系统中石英音叉结构参数设计提供参考。     

The Determination of Protonation Constants of Some Amino Acids and Their Esters by Potentiometry in Different Media

In this study, stoichiometric protonation constants of L-tyrosine, L-cysteine, L-tryptophane, L-lysine, and L-histidine, and their methyl and ethyl esters in water and ethanol–water mixtures of 30, 50, and 70% ethanol (v/v), were determined potentiometrically using a combined pH electrode system calibrated as the concentration of hydrogen ion. Titrations were performed at 25^∘C and the ionic strength of the medium was maintained at 0.10 mol⋅L⁻¹ using sodium chloride. Protonation constants were calculated by using the BEST computer program. The effect of solvent composition on the protonation constants is discussed. The log₁₀ K₂ values of esters generally decreased with increasing ethanol content. However, the log₁₀ K₁ values of the esters of L-tyrosine, L-cysteine, and L-tryptophane were found to increase with increasing ethanol content in contrast those of L-lysine and L-histidine esters.     

非线性互补问题的一种全局收敛的显式光滑Newton方法

本针对Po函数非线性互补问题，给出了一种显式光滑Newton方法，该方法将光滑参数μ进行显式迭代而不依赖于Newton方向的搜索过程，并在适当的假设条件下，证明了算法的全局收敛性。     

电路过渡过程的计算机辅助分析

电路过渡过程所列方程是微分方程，本文中采用的是方框图模型分析法，即将微分方程的复杂示解分解成最基本的加（减）、乘（除）、积分（微分）、增益等运算，采用VB设计用户界面产进行计算，并给出了一算例。     

The number of powers of 2 in a representation of large even integers II

It is proved that for any integerk≥ 54 000, there isN _k >0 depending onk only such that every even integer ≥N _k is a sum of two odd prime numbers andk powers of 2. Project partially supported by RGC Research Grant (No.HKU 7122/97P) and Post-Doctoral Fellowship of the University of Hong Kong.     

Surfaces of Constant Mean Curvature Bounded by Two Planar Curves

In this paper we study constant mean curvature compact surfaces with two Jordan curves in parallel planes as boundary and we investigate the point at which the surface inherits the symmetries of its boundary.     

A semi-analytical numerical method for time-dependent radiative transfer problems in slab geometry with coherent isotropic scattering

We describe a semi-analytical numerical method for coherent isotropic scattering time-dependent radiative transfer problems in slab geometry. This numerical method is based on a combination of two classes of numerical methods: the spectral methods and the Laplace transform (LTS_N) methods applied to the radiative transfer equation in the discrete ordinates (S_N) formulation. The basic idea is to use the essence of the spectral methods and expand the intensity of radiation in a truncated series of Laguerre polynomials in the time variable and then solve recursively the resulting set of “time-independent” S_N problems by using the LTS_N method. We show some numerical experiments for a typical model problem.