Pedestrian flows and non‐classical shocks

We present a model for the flow of pedestrians that describes features typical of this flow, such as the fall due to panic in the outflow of people through a door. The mathematical techniques essentially depend on the use of non‐classical shocks in scalar conservation laws. Copyright © 2005 John Wiley & Sons, Ltd.     

可反映射的性质及其应用

对可反映射的性质进行了研究 ,利用其性质对混沌系统 Hénon映射的同宿轨进行了讨论     

A convective weakly viscoelastic rotating flow with pressure Neumann condition

The objective of this work is to investigate through the numeric simulation, the effects of the weakly viscoelastic flow within a rotating rectangular duct subject to a buoyancy force due to the heating of one of the walls of the duct. A direct velocity–pressure algorithm in primitive variables with a Neumann condition for the pressure is employed. The spatial discretization is made with finite central differences on a staggered grid. The pressure field is directly updated without any iteration. Numerical simulations were done for several Weissemberg numbers (We) and Grashof numbers (Gr) . The numerical results show that for high Weissemberg numbers (We>7.4 × 10^?5) and for ducts with aspect ratio 2:1 and 8:1, the secondary flow is restabilized with a stretched double vortex configuration. It is also observed that when the Grashof number is increased (Gr>17 × 10^?4) , the buoyancy force neutralizes the effects of the Coriolis force for ducts with aspect ratio 8:1. Copyright © 2008 John Wiley & Sons, Ltd.     

A compact high‐order HLL solver for nonlinear hyperbolic systems

We present a new finite‐volume method for calculating complex flows on non‐uniform meshes. This method is designed to be highly compact and to accurately capture all discontinuities that may arise within the solution of a nonlinear hyperbolic system. In the first step, we devise a fourth‐degree Hermite polynomial to interpolate the solution. The coefficients defining this polynomial are calculated by using a least‐square method. To introduce monotonicity conditions within the procedure, two constraints are added into the least‐square system. Those constraints are derived by locally matching the high‐order Hermite polynomial with a low‐order TVD polynomial. To emulate these constraints only in regions of discontinuities, data‐depending weights are defined; these weights are based upon normalized indicators of smoothness of the solution and are parameterized by an O(1) quantity. The reconstruction so generated is highly compact and is fifth‐order accurate when the solution is smooth; this reconstruction becomes first order in regions of discontinuities. In the second step, this reconstruction is inserted in an HLL approximate Riemann solver. This solver is designed to correctly capture all discontinuities that may arise into the solution. To this aim, we introduce the contribution of a possible contact discontinuity into the HLL Riemann solver. Thus, a spatially fifth‐order non‐oscillatory method is generated. This method evolves in time the solution and its first derivative. In a one‐dimensional context, a linear spectral analysis and extensive numerical experiments make it possible to assess the robustness and the advantages of the method in computing multi‐scale problems with embedded discontinuities. Copyright © 2008 John Wiley & Sons, Ltd.     

About the compatibility between ansatzes and constraints for a local formulation of orbital‐free density functional theory

Functional properties that are exact for the Hohenberg–Kohn functional may turn into mutually exclusive constraints at a given level of ansatz. This is exemplarily shown for the local density approximation. Nevertheless, it is possible to reach exactly the Kohn–Sham data from an orbital‐free density functional framework based on simple one‐point functionals by starting from the Levy–Perdew–Sahni formulation. The energy value is obtained from the density‐potential pair, and therefore does not refer to the functional dependence of the potential expression. Consequently, the potential expression can be obtained from any suitable model and is not required to follow proper scaling behavior.     

Novel ‘Nano in Nano’ Composites for Sustained Drug Delivery: Biodegradable Nanoparticles Encapsulated into Nanofiber Non‐Wovens

Novel ‘nano in nano’ composites consisting of biodegradable polymer nanoparticles incorporated into polymer nanofibers may efficiently modulate drug delivery. This is shown here using a combination of model compound‐loaded biodegradable nanoparticles encapsulated in electrospun fibers. The dye coumarin 6 is used as model compound for a drug in order to simulate drug release from loaded poly(lactide‐co‐glycolide) nanoparticles. Dye release from the nanoparticles occurs immediately in aqueous solution. Dye‐loaded nanoparticles which are encapsulated by electrospun polymer nanofibers display a significantly retarded release.

     

Dust Self‐Organized Structures II. Solutions of Master Equations for Small Diffusion

Master equations for spherical dust structures are solved numerically using the asymptotic solutions at the center of the structures for the case of absence of external ionization and small diffusions. The structures are determined by a single parameter, the external plasma flux at the surface of the structure. The equilibrium states that are possible in a limited range of this parameter are investigated numerically. It is demonstrated that in the range of existence of equilibria the structures are changing their shapes and type of distributions inside the structures. For large external fluxes the ion and dust distributions can have peaks inside the structures while for low external fluxes the dust distribution has a single maximum at the structure center. The lower is the external flux supporting the structure the larger is its size. An increase of the external flux decreases the accumulation of dust and ions at the center. The total number of dust confined by the structure is larger for larger size structures. Estimates of dust crystallization inside structures are given. The role of diffusion is calculated by perturbations and is shown to be small in all structure regions except the structure edges. In the perturbation theory we use the exact expressions of the diffusion coefficients calculated previously numerically. The regions with dust density peaks inside the structures have been calculated with two order of magnitude larger precision that allows to resolve the structure parameter dependencies inside the peaks. It is shown that although in peaks the gradients of all parameters are increased the diffusion flux is still small and that the continuity and hydrodynamic approach are applicable within an accuracy about several %‐s (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

有限元仿真兰姆波评价胶接接头界面弹性模量的研究

金属薄板检测中,声-超声检测方法可激发多模式的兰姆波,兰姆波能够评价金属胶接后结构质量,正确认识兰姆波传播特性与胶接接头弹性模量的关系是有效应用该技术的基础。本文采用有限元仿真法分析这种关系,得出随着胶层弹性模量的降低,各模式兰姆波通过胶接接头后能量有不同程度的衰减,频率为0.5 MHz和1 MHz时,兰姆波主要模式A₀、S₀能量值与接头弹性模量具有一定的线性关系。实验结果表明,低频时随着胶接接头的人工老化,接头中主要模式兰姆波能量均发生一定程度的衰减,其能量变化曲线与有限元仿真结果趋势一致。本方法能为实验研究、实际检测提供有效参考依据。