Extension of a combined analytical/numerical initial value problem solver for unsteady periodic flow

Here we describe analytical and numerical modifications that extend the Differential Reduced Ejector/ mixer Analysis (DREA), a combined analytical/numerical, multiple species ejector/mixing code developed for preliminary design applications, to apply to periodic unsteady flow. An unsteady periodic flow modelling capability opens a range of pertinent simulation problems including pulse detonation engines (PDE), internal combustion engine ICE applications, mixing enhancement and more fundamental fluid dynamic unsteadiness, e.g. fan instability/vortex shedding problems. Although mapping between steady and periodic forms for a scalar equation is a classical problem in applied mathematics, we will show that extension to systems of equations and, moreover, problems with complex initial conditions are more challenging. Additionally, the inherent large gradient initial condition singularities that are characteristic of mixing flows and that have greatly influenced the DREA code formulation, place considerable limitations on the use of numerical solution methods. Fortunately, using the combined analytical–numerical form of the DREA formulation, a successful formulation is developed and described. Comparison of this method with experimental measurements for jet flows with excitation shows reasonable agreement with the simulation. Other flow fields are presented to demonstrate the capabilities of the model. As such, we demonstrate that unsteady periodic effects can be included within the simple, efficient, coarse grid DREA implementation that has been the original intent of the DREA development effort, namely, to provide a viable tool where more complex and expensive models are inappropriate. Copyright © 2002 John Wiley & Sons, Ltd.     

Flow in a waterfall with large Froude number

An approximate method is developed to solve the full nonlinear equations governing two-dimensional irrotational flow in a free waterfall, falling under the influence of gravity, at high Froude number based on conditions far upstream. Schwarz—Christoffel transformation is used to map the region, in the complex potential-plane, onto the upper half-plane. The Hilbert transformation as well as the perturbation technique, for large Froude number, are used as a basis for the approximate solution of the problem. A complete solution, up to second-order approximation, for the downstream free-surfaces profiles, for different Froude number, is discussed and illustrated. The obtained approximate solutions are compared with those of other authors. Favourable agreement with other results suggests that this method is effective in dealing with flow problems strongly influenced by gravity and high Froude number. The results obtained by this method are sufficiently accurate for practical purposes.     

Effect of particle charging on momentum and heat transfer from rarefied plasma flow

The features of interaction of a spherical metallic particle with a rarefied thermal plasma flow due to the presence o charges-electrons and ions in the gaseous phase-are considered. Analytical expressions describing charge, momentum, and energy exchange between the plasma and the particle für the cases of strong and weak Debye screening are obtained. It is illustrated that the efficiency of particle heating in the plasma considerably grows as compared with a hot molecular gas due to participation of electrons and ions in file transfer processes.     

A numerical study of an unsteady laminar flow in a doubly constricted 3D vessel

Unsteady flow dynamics in doubly constricted 3D vessels have been investigated under pulsatile flow conditions for a full cycle of period T. The coupled non‐linear partial differential equations governing the mass and momentum of a viscous incompressible fluid has been numerically analyzed by a time accurate Finite Volume Scheme in an implicit Euler time marching setting. Roe's flux difference splitting of non‐linear terms and the pseudo‐compressibility technique employed in the current numerical scheme makes it robust both in space and time. Computational experiments are carried out to assess the influence of Reynolds' number and the spacing between two mild constrictions on the pressure drop across the constrictions. The study reveals that the pressure drop across a series of mild constrictions can get physiologically critical and is also found to be sensitive both to the spacing between the constrictions and the oscillatory nature of the inflow profile. The flow separation zone on the downstream constriction is seen to detach from the diverging wall of the constriction leading to vortex shedding with 3D features earlier than that on the wall in the spacing between the two constrictions. Copyright © 2002 John Wiley & Sons, Ltd.     

Quantification of mass flow effects in capillary SFC

A model for restrictor flow produces accurate predictions of flow that can be used to optimize restrictor design. The relative amplitudes of restrictor and other flows and their effect on efficiency are discussed.     

CSR外靶实验终端中子流测量的模拟

随着兰州重离子加速器冷却储存环(HIRFL-CSR)外靶实验终端中子墙的建立, 为实验测量高能中子提供了机遇. 为确定CSR外靶实验终端对中子流测量的可行性, 基于BUU理论模型分别对对称系统(Ni+Ni, Pb+Pb)和非对称系统(Pb+Ni)进行了模拟计算, 发现当系统能量达到几百MeV/u时, 中子流信号相当明显, 并与碰撞参数有明显的依赖关系. 模拟结果表明, 在前角20°的覆盖范围内, 可以较好地实现中子流测量所需要的反应平面确定及碰撞参数选择. 对双击事件及其对中子流的影响进行了简单的讨论.     

流动注射化学发光抑制法测定尼群地平

基于尼群地平在碱性条件下对铁氰化钾鲁米诺化学发光体系有较强的抑制作用,联用流动注射技术建立了测定尼群地平的新方法.本法的线形范围为0.1～100.0 μg/mL,检出限0.02 μg/mL ,此法用于尼群地平片剂含量测量和药典法对照结果满意.     

Precise lattice parameter measurements in CsCl-Br solid solutions between room temperature and 90°K by powder x-ray diffraction

Solid solutions of CsCl-Br in five different concentrations were prepared in sealed quartz tubes by heating the mixture to 1123°K for 6–8 hr and quenching to room temperature. X-ray diffractograms were taken at eight different temperatures between room temperature and 90°K for these solid solutions using the YPC50NM powder diffractometer and a continuous flow cryostat. The observed lattice parameters for each sample at each temperature obtained from the powder diffractograms were then extrapolated to give the true lattice parameters using the least square method with Nelson-Riley extrapolation scheme. The values of the true lattice parameters at each concentration and at each temperature were tabulated and the results discussed. It is shown that the lattice parameters vs temperature for some concentrations exhibit an anomalous behaviour. Contribution No. 691     

Modeling plane turbulent Couette flow

Among the salient features of shear-driven plane Couette flow is the constancy of the total shear stress (viscous and turbulent) across the flow. This constancy gives rise to a quasi-homogenous core region, which makes the bulk of the flow substantially different from pressure-driven Poiseuille flow. The present second-moment closure study addresses the conflicting hypotheses relating to turbulent Couette flow. The inclusion of a new wall-proximity function in the wall-reflection part of the pressure-strain model seems mandatory, and the greement with recent experimental and direct numerical simulation (DNS) results is encouraging. Analysis of model computations in the range 750 ≤ Re ≤ 35,000 and comparisons with low-Re DNS data suggest that plane Couette flow exhibits a local-equilibrium core region, in which anisotropic, homogeneous turbulence prevails. However, the associated variation of the mean velocity in the core, as obtained by the model, conflicts with the intuitively appealing assumption of homogeneous mean shear. The constancy of the velocity gradient exhibited by the DNS therefore signals a deficiency in the modeled transport equation for the energy dissipation rate.     

Optimal control of unsteady compressible viscous flows

The control of complex, unsteady flows is a pacing technology for advances in fluid mechanics. Recently, optimal control theory has become popular as a means of predicting best case controls that can guide the design of practical flow control systems. However, most of the prior work in this area has focused on incompressible flow which precludes many of the important physical flow phenomena that must be controlled in practice including the coupling of fluid dynamics, acoustics, and heat transfer. This paper presents the formulation and numerical solution of a class of optimal boundary control problems governed by the unsteady two‐dimensional compressible Navier–Stokes equations. Fundamental issues including the choice of the control space and the associated regularization term in the objective function, as well as issues in the gradient computation via the adjoint equation method are discussed. Numerical results are presented for a model problem consisting of two counter‐rotating viscous vortices above an infinite wall which, due to the self‐induced velocity field, propagate downward and interact with the wall. The wall boundary control is the temporal and spatial distribution of wall‐normal velocity. Optimal controls for objective functions that target kinetic energy, heat transfer, and wall shear stress are presented along with the influence of control regularization for each case. Copyright © 2002 John Wiley & Sons, Ltd.