An efficient edge based data structure for the compressible Reynolds-averaged Navier–Stokes equations on hybrid unstructured meshes

An efficient edge based data structure has been developed in order to implement an unstructured vertex based finite volume algorithm for the Reynolds-averaged Navier–Stokes equations on hybrid meshes. In the present approach, the data structure is tailored to meet the requirements of the vertex based algorithm by considering data access patterns and cache efficiency. The required data are packed and allocated in a way that they are close to each other in the physical memory. Therefore, the proposed data structure increases cache performance and improves computation time. As a result, the explicit flow solver indicates a significant speed up compared to other open-source solvers in terms of CPU time. A fully implicit version has also been implemented based on the PETSc library in order to improve the robustness of the algorithm. The resulting algebraic equations due to the compressible Navier–Stokes and the one equation Spalart–Allmaras turbulence equations are solved in a monolithic manner using the restricted additive Schwarz preconditioner combined with the FGMRES Krylov subspace algorithm. In order to further improve the computational accuracy, the multiscale metric based anisotropic mesh refinement library PyAMG is used for mesh adaptation. The numerical algorithm is validated for the classical benchmark problems such as the transonic turbulent flow around a supercritical RAE2822 airfoil and DLR-F6 wing-body-nacelle-pylon configuration. The efficiency of the data structure is demonstrated by achieving up to an order of magnitude speed up in CPU times.     

Bifurcation analysis of visual angle model with anticipated time and stabilizing driving behavior

In the light of the visual angle model (VAM), an improved car-following model considering driver's visual angle, anticipated time and stabilizing driving behavior is proposed so as to investigate how the driver's behavior factors affect the stability of the traffic flow. Based on the model, linear stability analysis is performed together with bifurcation analysis, whose corresponding stability condition is highly fit to the results of the linear analysis. Furthermore, the time-dependent Ginzburg-Landau (TDGL) equation and the modified Korteweg-de Vries (mKdV) equation are derived by nonlinear analysis, and we obtain the relationship of the two equations through the comparison. Finally, parameter calibration and numerical simulation are conducted to verify the validity of the theoretical analysis, whose results are highly consistent with the theoretical analysis.     

桉树杂交种与其亲本的近红外光谱判别

研究桉树控制授粉后目标性状的基因作用方式是探索其基因重组规律的重要内容。常规的数量统计分析精度往往不高,而DNA分析的专业要求高,且费时费力。该研究利用近红外光谱(NIRs)研究不同基因型桉树杂交种、亲本及杂交种与亲本间近红外光谱信息的关系,探索NIRs用于桉树杂交种与其亲本判别的可行性和准确性。以控制授粉的桉树亲本及其杂交F1代材料为对象,每种基因型从各自田间试验分别选取10个单株,采集树冠中上部新鲜健康叶片。用手持式近红外仪Phazir Rx(1624)采集桉树杂交种与其亲本叶片的NIRs信息。每单株选10片完全生理成熟的健康叶片,避开叶脉扫描其正面光谱5次,以50条NIRs信息的均值代表单个叶片的NIRs信息,最终每个基因型获得10条NIRs信息。对原始NIRs采用二阶多项式S.G一阶导数预处理。预处理后的NIRs用于多元统计分析,首先对桉树杂交亲本和子代样本进行主成分分析(PCA),直观展示不同基因型的分类情况。然后运用簇类独立软模式(SIMCA)和偏最小二乘判别分析(PLS-DA)两种有监督的判别模式验证NIRs用于桉树杂交种与其亲本树种的分类判别效果。PCA结果显示,不同的亲本间、杂交种间及杂交种与亲本间样本的主因子得分可以清晰地将各基因型分开。SIMCA模式判别分析中,桉树杂交种样本到亲本PCA模型的样本距离显示,待判别样本能够形成单独的聚类,且能直观反映两者的遗传相似。PLS-DA判别结果显示,桉树杂交亲本的PLS模型能通过预测其杂交子代的响应变量将其与亲本准确分开。结果表明,桉树叶片的NIRs信息可以准确地反映桉树杂交子代遗传信息的传递规律,NIRs判别模型可以准确地将各种基因型予以区分。因此,NIRs信息不仅可用于桉树杂交种和纯种的定性判别,还可以分析桉树基因重组过程中加性遗传效应的大小,从而为桉树遗传基础分析及其育种改良研究提供理论支撑。     

A probabilistic approach to spectral analysis of growth-fragmentation equations

The growth-fragmentation equation describes a system of growing and dividing particles, and arises in models of cell division, protein polymerisation and even telecommunications protocols. Several important questions about the equation concern the asymptotic behaviour of solutions at large times: at what rate do they converge to zero or infinity, and what does the asymptotic profile of the solutions look like? Does the rescaled solution converge to its asymptotic profile at an exponential speed? These questions have traditionally been studied using analytic techniques such as entropy methods or splitting of operators. In this work, we present a probabilistic approach: we use a Feynman–Kac formula to relate the solution of the growth-fragmentation equation to the semigroup of a Markov process, and characterise the rate of decay or growth in terms of this process. We then identify the Malthus exponent and the asymptotic profile in terms of a related Markov process, and give a spectral interpretation in terms of the growth-fragmentation operator and its dual.     

Mass distribution and skewness for passive scalar transport in pipes with polygonal and smooth cross sections

We extend our previous results characterizing the loading properties of a diffusing passive scalar advected by a laminar shear flow in ducts and channels to more general cross‐sectional shapes, including regular polygons and smoothed corner ducts originating from deformations of ellipses. For the case of the triangle and localized, cross‐wise uniform initial distributions, short‐time skewness is calculated exactly to be positive, while long‐time asymptotics shows it to be negative. Monte Carlo simulations confirm these predictions, and document the timescale for sign change. The equilateral triangle appears to be the only regular polygon with this property—all others possess positive skewness at all times. Alternatively, closed‐form flow solutions can be constructed for smooth deformations of ellipses, and illustrate how both nonzero short‐time skewness and the possibility of multiple sign switching in time is unrelated to domain corners. Exact conditions relating the median and the skewness to the mean are developed which guarantee when the sign for the skewness implies front (more mass to the right of the mean) or back (more mass to the left of the mean) “loading” properties of the evolving tracer distribution along the pipe. Short‐ and long‐time asymptotics confirm this condition, and Monte Carlo simulations verify this at all times. The simulations are also used to examine the role of corners and boundaries on the distribution for short‐time evolution of point source , as opposed to cross‐wise uniform, initial data.     

Triesterase and Promiscuous Diesterase Activities of a Di‐CoII‐Containing Organophosphate Degrading Enzyme Reaction Mechanisms

The reaction mechanism for the hydrolysis of trimethyl phosphate and of the obtained phosphodiester by the di‐Co^II derivative of organophosphate degrading enzyme from Agrobacterium radiobacter P230(OpdA), have been investigated at density functional level of theory in the framework of the cluster model approach. Both mechanisms proceed by a multistep sequence and each catalytic cycle begins with the nucleophilic attack by a metal‐bound hydroxide on the phosphorus atom of the substrate, leading to the cleavage of the phosphate‐ester bond. Four exchange‐correlation functionals were used to derive the potential energy profiles in protein environments. Although the enzyme is confirmed to work better as triesterase, as revealed by the barrier heights in the rate‐limiting steps of the catalytic processes, its promiscuous ability to hydrolyze also the product of the reaction has been confirmed. The important role played by water molecules and some residues in the outer coordination sphere has been elucidated, while the binuclear Co^II center accomplishes both structural and catalytic functions. To correctly describe the electronic configuration of the d shell of the metal ions, high‐ and low‐spin arrangement jointly with the occurrence of antiferromagnetic coupling, have been herein considered.     

Influence of Polarity and Activation Energy in Microwave–Assisted Organic Synthesis (MAOS)

The aim of this work was to determine the parameters that have decisive roles in microwave-assisted reactions and to develop a model, using computational chemistry, to predict a priori the type of reactions that can be improved under microwaves. For this purpose, a computational study was carried out on a variety of reactions, which have been reported to be improved under microwave irradiation. This comprises six types of reactions. The outcomes obtained in this study indicate that the most influential parameters are activation energy, enthalpy, and the polarity of all the species that participate. In addition to this, in most cases, slower reacting systems observe a much greater improvement under microwave irradiation. Furthermore, for these reactions, the presence of a polar component in the reaction (solvent, reagent, susceptor, etc.) is necessary for strong coupling with the electromagnetic radiation. We also quantified that an activation energy of 20–30 kcal mol⁻¹ and a polarity (μ) between 7–20 D of the species involved in the process is required to obtain significant improvements under microwave irradiation.     

Feasibility of application of conductometric biosensor based on acetylcholinesterase for the inhibitory analysis of toxic compounds of different nature

This study was aimed at the development of a conductometric biosensor based on acetylcholinesterase considering the feasibility of its application for the inhibitory analysis of various toxicants. In this paper, the optimum conditions for enzyme immobilization on the transducer surface are selected as well as the optimum concentration of substrate for inhibitory analysis. Sensitivity of the developed biosensor to different classes of toxic compounds (organophosphorus pesticides, heavy metal ions, surfactants, aflatoxin, glycoalkaloids) was tested. It is shown that the developed biosensor can be successfully used for the analysis of pesticides and mycotoxins, as well as for determination of total toxicity of the samples. A new method of biosensor analysis of toxic substances of different classes in complex multicomponent aqueous samples is proposed.     

稳定分层湍流中逆梯度输运的关联分析研究

本文采用关联分析方法研究了稳定温度分层湍流中的结构特性、输运特性，以及热量、动量逆梯度输运现象的尺度效应及其参数演化．首先采用大涡模拟方法对稳定分层湍流中的结构特性和输运特性进行了分析，将逆梯度输运发生的时间尺度作为已知条件，结合关联量分析方法在波数空间中的解析解，对逆梯度输运现象的尺度效应进行了分析研究．结果发现，稳定分层强度较大的流动中发生垂向热量及动量逆梯度输运现象，发生的结构尺度与关联分析所发现垂向热量、动量逆梯度输运的波数形成了呼应．随着分层强度增加，热量、动量的输运强度均受抑制，与逆梯度输运关联的流场结构尺度减小，同样的效应也发生在流场结构向下游演化的过程中．