A discontinuous Galerkin method/HLLC solver for the Euler equations

This paper proposes a fully three‐dimensional non‐linear Euler methodology for solving aerodynamic and acoustic problems in the presence of strong shocks and rarefactions. It uses a discontinuous Galerkin method (DGM) within the element, and a Riemann solver (HLLC) at the boundaries to propagate rarefactions while preserving the entropy condition and capturing shocks with no spurious oscillations. This approach is thought to marry the best aspects of finite element and finite volume methods, achieving conservation while not requiring the solution of a large matrix. Examples in which shock and rarefaction waves are well captured are presented and the propagation of acoustic pulses is well demonstrated. Copyright © 2003 John Wiley & Sons, Ltd.     

Stability analysis of abnormal multiplication of plankton using parameter identification technique

This paper presents the mathematical approach for the abnormal multiplication of plankton. An abnormal multiplication can be expressed as an unstable problem and the stability of the system is investigated by introducing eigenvalues of a mathematical equation. The stability of the system can be judged by an eigenvalue based on the Lyapunov's stability theory. In this paper, the Arnoldi‐QR method is used to obtain eigenvalues and eigenvectors of the system. The mode superposition method is employed to create spatial distribution needed to analyse the stability. To obtain the objective eigenvalue, the parameter identification technique is employed. The finite element method is used for the discretization in space. Lake Kasumigaura, which is located in Ibaraki Prefecture in Japan, is selected and actual data in 1975, 1976, 1991 and 2000 are used in order to investigate the stability of the specified lake in Japan. Copyright © 2004 John Wiley & Sons, Ltd.     

Indirect boundary element method for unsteady linearized flow over prolate and oblate spheroids and hemispheroidal protuberances

The indirect boundary element method was used to study the hydrodynamics of oscillatory viscous flow over prolate and oblate spheroids, and over hemispheroidal bodies hinged to a plate. Analytic techniques, such as spheroidal coordinates, method of images, and series representations, were used to make the numerical methods more efficient. A novel method for computing the hydrodynamic torque was used, since for oscillatory flow the torque cannot be computed directly from the weightings. Instead, a Green's function for torque was derived to compute the torque indirectly from the weightings. For full spheroids, the method was checked by comparing the results to exact solutions at low and high frequencies, and to results computed using the singularity method. For hemispheroids hinged to a plate, the method for low frequencies was checked by comparing the results to previous results, and to exact solutions at high frequencies. Copyright © 2004 John Wiley & Sons, Ltd.     

Numerical analysis of the rotating viscous flow approaching a solid sphere

A numerical simulation is performed to investigate the flow induced by a sphere moving along the axis of a rotating cylindrical container filled with the viscous fluid. Three‐dimensional incompressible Navier–Stokes equations are solved using a finite element method. The objective of this study is to examine the feature of waves generated by the Coriolis force at moderate Rossby numbers and that to what extent the Taylor–Proudman theorem is valid for the viscous rotating flow at small Rossby number and large Reynolds number. Calculations have been undertaken at the Rossby numbers (R_o) of 1 and 0.02 and the Reynolds numbers (R_e) of 200 and 500. When R_o=O(1), inertia waves are exhibited in the rotating flow past a sphere. The effects of the Reynolds number and the ratio of the radius of the sphere and that of the rotating cylinder on the flow structure are examined. When R_o ? 1, as predicted by the Taylor–Proudman theorem for inviscid flow, the so‐called ‘Taylor column’ is also generated in the viscous fluid flow after an evolutionary course of vortical flow structures. The initial evolution and final formation of the ‘Taylor column’ are exhibited. According to the present calculation, it has been verified that major theoretical statement about the rotating flow of the inviscid fluid may still approximately predict the rotating flow structure of the viscous fluid in a certain regime of the Reynolds number. Copyright © 2004 John Wiley & Sons, Ltd.     

有限元节点编号优化

从分析最简单最优的矩形单元节点编号着手，分析其特点，得出各个节点的相邻单元节点编号总和除以该节点相邻单元数所得的节点商以及相邻单元最小最大节点编号和都为升序排列。总结出一种新型的有限元节点重编号方法，并举例说明该算法能减少有限元的带宽。     

Least‐squares finite element method on adaptive grid for PDEs with shocks

We apply the least‐squares finite element method with adaptive grid to nonlinear time‐dependent PDEs with shocks. The least‐squares finite element method is also used in applying the deformation method to generate the adaptive moving grids. The effectiveness of this method is demonstrated by solving a Burgers' equation with shocks. Computational results on uniform grids and adaptive grids are compared for the purpose of evaluation. The results show that the adaptive grids can capture the shock more sharply with significantly less computational time. For moving shock, the adaptive grid moves correctly with the shock. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Conformation of Tax-response elements in the human T-cell leukemia virus type I promoter

Background: HTLV I Tax is believed to activate viral gene expression by binding bZIP proteins (such as CRIB) and increasing their affinities for proviral THE target sites. Each 21 by THE target site contains an imperfect copy of the intrinsically bent CRE target site (the TRE core) surrounded by highly conserved flanking sequences. These flanking sequences are essential for maximal increases in DNA affinity and transactivation, but they are not, apparently, contacted by protein. Here we employ non-denaturing gel electrophoresis to evaluate TRE conformation in the presence and absence of bZIP proteins, and to explore the role of DNA conformation in viral transactivation.Results: Our results show that the TRE-1 flanking sequences modulate the structure and modestly increase the affinity of a CREB bZIP peptide for the TRE-1 core recognition sequence. These flanking sequences are also essential for a maximal increase in stability of the CREB-DNA complex in the presence of Tax.Conclusions: The CRE-like TRE core and the TRE flanking sequences are both essential for formation of stable CRIB-TRE-1 and Tax-CREB-TRE-1 complexes. These two DNA segments may have co-evolved into a unique structure capable of recognizing Tax and a bZIP protein.     

三重四极杆电感耦合等离子体质谱(ICP-MS)法测定氧化铕中13种稀土杂质元素

稀土杂质元素直接影响高纯单金属稀土材料的整体性能,是高科技领域许多材料的重要组成部分。通过考察最佳的消解酸量、温度、时间、氧气反应气流量、稀释气流量,建立了基于三重四极杆电感耦合等离子体质谱仪(ICP-MS/MS)直接测定氧化铕中13种稀土杂质元素分析方法。该方法采用0.1%基体直接进样,可以很大程度提高前处理分析效率。利用碰撞模式测定氧化铕稀土中的Y、La、Pr、Nd、Sm、Gd、Tb、Dy、Ho、Er、Yb、Lu元素,氧气质量转移模式测定氧化铕中的Tm,两种模式结合可以有效去除多原子干扰,实现氧化铕的稳定测试分析。通过对氧化铕标准物质（GBW02902）直接测定分析,结果表明,在碰撞和氧气质量转移模式下,各元素线性相关系数（r）均大于 0.9999,方法检出限为0.001~0.023 mg/kg,测试精密度优于1.99%,13种元素的测试值都在认定值的不确定度范围之内。该分析方法操作简单,测试稳定,效率高,为实验室进行氧化铕材料中REE杂质的准确测试分析提供思路和借鉴。     

激光诱导击穿光谱(LIBS)技术在低碳钢铁冶金行业的应用

激光诱导击穿光谱技术(Laser induced breakdown spectroscopy,简称LIBS)是一种基于激光等离子体的发射光谱技术,具有样品制备简单,分析迅速,多元素同时分析,现场在线和远程分析等独特优势,已在钢铁冶金领域取得了越来越广泛的应用与发展。综述了近年来LIBS技术在钢铁冶金领域的应用研究进展,从冶金生产的全流程进行介绍,包括矿石、熔融钢水/铁水、钢铁产品及炉渣和废气分析。同时,还总结了LIBS技术在钢铁冶金领域应用的优点和缺点,并对该技术的应用前景及未来发展方向进行展望。     

蝙蝠喉部发声组织建模及发声机理*

建立蝙蝠发声组织模型对超声机理研究及在智能设备的应用具有重要意义。根据蝙蝠喉部发声组织结构特点,通过有限元方法构建了蝙蝠的3种不同发声组织模型,分析了尺寸、材料力学参数、组织结构和张力4个因素对发声组织特征频率的影响。结果表明,如果用人类声带,按比例缩小构建蝙蝠喉部组织模型,蝙蝠无法发出超声波。构建组织结构含甲状软骨和声带的半鼓状模型和只含声带的条状模型,两种模型的特征频率相近且在合理的参数域内均无法达到超声范围。而含膜条状模型的特征频率可以通过张力进行超声频率的调节,这与文献的实验结果一致。因此,可基于含膜条状模型对蝙蝠喉管发声组织进行建模及其发声机理研究。