Numerical simulation of DNA sample preconcentration in microdevice electrophoresis

A numerical model is presented for the accurate and efficient prediction of preconcentration and transport of DNA during sample introduction and injection in microcapillary electrophoresis. The model incorporates conservation laws for the different buffer ions, salt ions, and DNA sample, coupled through a Gaussian electric field to account for the field modifications that cause electromigration. The accuracy and efficiency required to capture the physics associated with such a complex transient problem are realized by the use of the finite element-flux corrected transport (FE-FCT) algorithm in two dimensions. The model has been employed for the prediction of DNA sample preconcentration and transport during electrophoresis in a double-T injector microdevice. To test its validity, the numerical results have been compared with the corresponding experimental data under similar conditions, and excellent agreement has been found. Finally, detailed results from a simulation of DNA sample preconcentration in electrophoretic microdevices are presented using as parameters the electric field strength and the other species concentrations. The effect of the Tris concentration on sample stacking is also investigated. These results demonstrate the great potential offered by the model for future optimization of such microchip devices with respect to significantly enhanced speed and resolution of sample separation.     

神经网络的函数逼近能力分析

本文综述了多层前传网络（ＭＬＰ）及径向基函数网络（ＲＢＦ）对函数任意精度逼近的能力，比较了两种网络的最佳逼近特性。对激活函数类的扩充作了介绍，并说明有限数值精度对函数逼近能力实现的影响。     

高空核电磁脉冲模拟波形的双指数函数拟合法

 许多标准和公开出版物中都用双指数函数描述高空核电磁脉冲典型波形。通过数值方法，研究了双指数函数一项重要的性质。根据该性质，详细讨论了高空核电磁脉冲模拟波形数值拟合中，双指数函数特征参数与脉冲峰值、前沿、后沿以及半宽等物理参数的关系，从而提出一种简单有效的脉冲参数计算方法。通过IEC标准中规定的高空核电磁脉冲参数的估计与一个实测高空核电磁脉冲模拟波形的数值拟合，验证了该方法的有效性和可靠性。该方法能够适应多种峰值、前沿和半宽的高宽核电磁脉冲模拟波形的数值拟合。     

Numerical solutions for jump-diffusions with regime switching

This paper is devoted to numerical solutions for a class of jump-diffusions with regime switching. After briefly reviewing the notion of jump-diffusions with regime switching, finite-difference procedures are constructed. Under simple conditions, it is proved that the algorithm converges to the desired limit by means of a martingale problem formulation. Numerical experiments are carried out to demonstrate the performance of the algorithm.     

Vibration of strings with nonlinear supports

The dynamic string motion, which displacement is unilaterally constrained by the rigid termination condition of an arbitrary geometry has been simulated and analyzed. The treble strings of a grand piano usually terminate at a capo bar, which is situated above the strings. The apex of a V-shaped section of the capo bar defines the end of the speaking length of the strings. A numerical calculation based on the traveling wave solution is proposed for modeling the nonlinearity inducing interactions between the vibrating string and the contact condition at the point of string termination. It was shown that the lossless string vibrates in two distinct vibration regimes. In the beginning the string starts to interact in a nonlinear fashion with the rigid terminator, and the resulting string motion is aperiodic. Consequently, the spectrum of the string motion depends on the amplitude of string vibrations, and its spectral structure changes continuously with the passage of time. The duration of that vibration regime depends on the geometry of the terminator. After some time of aperiodic vibration, the string vibrations settle in a periodic regime where the resulting spectrum remains constant.     

数值分析课程教学的实践与认识

本文针对数值分析课程的特点，结合自身的教学实践和体会，给出了实际教学过程中的三个方面的具体做法．     

Thermomechanical formulation of ductile damage coupled to nonlinear isotropic hardening and multiplicative viscoplasticity

In this paper, we present a thermomechanical framework which makes use of the internal variable theory of thermodynamics for damage-coupled finite viscoplasticity with nonlinear isotropic hardening. Damage evolution, being an irreversible process, generates heat. In addition to its direct effect on material's strength and stiffness, it causes deterioration of the heat conduction. The formulation, following the footsteps of Simó and Miehe (1992), introduces inelastic entropy as an additional state variable. Given a temperature dependent damage dissipation potential, we show that the evolution of inelastic entropy assumes a split form relating to plastic and damage parts, respectively. The solution of the thermomechanical problem is based on the so-called isothermal split. This allows the use of the model in 2D and 3D example problems involving geometrical imperfection triggered necking in an axisymmetric bar and thermally triggered necking of a 3D rectangular bar.     

Numerical analysis of the initial conditions in fractional systems

Fractional dynamics is a growing topic in theoretical and experimental scientific research. A classical problem is the initialization required by fractional operators. While the problem is clear from the mathematical point of view, it constitutes a challenge in applied sciences. This paper addresses the problem of initialization and its effect upon dynamical system simulation when adopting numerical approximations. The results are compatible with system dynamics and clarify the formulation of adequate values for the initial conditions in numerical simulations.     

Numerical study on combustion of H2—O2 in a supersonic reactive expansion in a nozzle

A study on expansion flow inside a nozzle considering full mechanism chemistry of hydrogen and oxygen was carried out. In this study, a full implicit scheme for turbulent reactive flow was obtained by combining the second order TVD scheme of Yee and Harten (1987, Implicit TVD schemes for hyperbolic conservation laws in curvilinear coordinates. American Institute of Aeronautics and Astronautics Journal, 25(2), 266–274) with the efficient implicit lower-upper scheme of Shuen and Yoon (1989, Numerical study of chemically reacting flows using a lower-upper symmetric successive overrelaxation scheme. American Institute of Aeronautics and Astronautics Journal, 27(12), 1752–1756). The species equations, Navier–Stokes equations and turbulence model were implemented in the numerical scheme and solved in conjunction with full detailed finite rate chemistry. The numerical scheme is verified by comparison with experimental results of a converging–diverging nozzle. Effects of inlet pressure, temperature and fuel-oxidant mass ratio on nozzle flow field were studied. Variation of chemical species under different conditions was investigated by considering a chemical mechanism. Results show that increasing inlet pressure increases the rate of reactions due to increasing the concentration of reactants. For lower inlet pressure the radical H increases slightly in the diverging part of the nozzle, while for higher pressures it decreases along the nozzle. Inlet fuel–oxidant mass ratio affects the variation of all species with a greater effect for a near stoichiometric ratio. It was also shown that a higher inlet temperature provides a more enhanced reaction zone in the diverging part of the nozzle.