A self-similar solution of exponential shock waves in non-ideal magnetogasdynamics

A self similar method is used to analyze numerically the one-dimensional, unsteady flow of a strong cylindrical shock wave driven by a piston moving with time according to an exponential law in a plasma of constant density. The plasma is assumed to be a non-ideal gas with infinite electrical conductivity permeated by an axial magnetic field. Numerical solutions in the region between the shock and the piston are presented for the cases of adiabatic and isothermal flow. The general behaviour of density, velocity, and pressure profiles remains unaffected due to presence of magnetic field in non-ideal gas. However, there is a decrease in values of density, velocity and pressure in case of magnetogasdynamics as compared to non-magnetic case. It may be noted that the effect of magnetic field on the flow pattern is more significant in case of isothermal flow as compared to adiabatic flow. The effect of non-idealness, specific heat exponent and magnetic field strength on the variation of shock strength across the shock front is also investigated.     

考虑活塞回复的侧向后喷武器两相流数值模拟

为研究活塞回复运动对火药燃气流动的影响,基于两相流理论对活塞控制侧向后喷武器的发射过程进行了数值模拟研究。考虑控制侧向后喷通道开闭的活塞-弹簧系统的往复运动,建立了结合膛内气固两相流、活塞腔内流固耦合和侧向排气管内气体瞬态流动的武器发射过程数学模型,并将数值模拟结果与相关文献进行了比较验证。得到了该武器发射过程中膛内流场分布与稀疏波传播特性,并与普通武器的膛内流场进行了对比分析。进一步研究了活塞回复运动对火药燃气流动和减后坐效率的影响。结果表明:相对于不考虑活塞的回复运动,在弹丸初速都降低1.52%的情况下,因为活塞回复关闭后喷通道,其减后坐效率由38.86%下降到32.88%,说明在此类武器研究中,不可忽视活塞回复运动。     

Panel flutter prediction in two dimensional flow with enhanced piston theory

Piston theory may be used in the high Mach number supersonic flow region and/or in very high frequency subsonic or supersonic flow. In this flow model, the pressure at a point on the fluid-solid interface only depends on the downwash at the same point. However the classical piston theory may not be sufficient for some phenomena in aeroelasticity and aeroacoustics (far field prediction). Dowell and Bliss have created an extension of piston theory that allows for higher order effects that take into account the effect the distribution of downwash on pressure at any point. For simple harmonic motion, expansions in reduced frequency, inverse reduced frequency and/or inverse (square of) Mach number have all been created; The effects of higher order terms in these several expansion in creating an enhanced piston theory was illustrated for plunge and pitch motion of an airfoil (discrete system) by Ganji and Dowell. In the present paper, flutter prediction for a flexible panel in two –dimensional flow is investigated using enhanced piston theory. The goal of the present paper is to demonstrate that an enhance version of piston theory can analyze single degree of freedom flutter of a panel as compared to the classical piston theory and quasi-steady aerodynamic models which can only treat coupled mode flutter.     

Unsteady unidirectional flow of an Oldroyd-B fluid in a circular duct with different given volume flow rate conditions

In this paper, the velocity profile and pressure gradient of the unsteady state unidirectional flow of an Oldroyd-B fluid in a circular duct are considered. The flow motion in the duct is induced by a given but arbitrary inlet volume flow rate which varies with time. Based on the flow conditions described, two basic flow situations are solved, which are a suddenly started, and a constant acceleration, flow respectively. These two results are applied to a practical case that is a trapezoidal piston motion which contains three phases of piston motion, the constant acceleration from the rest to a fixed velocity, then keeping at this velocity, following with the constant deceleration to a stop. In addition, oscillatory flow is also considered.     

The extrusion of a model yield stress fluid imaged by MRI velocimetry

Extrusion tests were performed by forcing a well-characterized model yield stress fluid from a cylindrical cartridge through various cylindrical extrusion dies using a variety of different piston velocities. In this study the Bingham number within the die ranged from 0.1 to 10. MRI techniques allowed for the non-invasive determination of the local velocity within the extruded material in the range [0.015; 20 mm s^?1]. The velocity profile within a very long die was determined by MRI and agreed very well with the analytical results for the flow of a Herschel–Bulkley fluid within a conduit using parameters determined from independent rheometrical tests, validating both the rheological approach and the accuracy of the MRI techniques. Although the velocity was determined by MRI in the upper and lower zones separately, the intersection of these zones showed great agreement, providing an entire view of the extrusion process. In the range of Bingham number studied, the velocity field for a given contraction ratio appeared similar when scaled by the piston velocity, with a dimpled acceleration zone above the die and lateral dead zones varying negligibly with the piston velocity. For a further analysis the experimental results were compared with the results of numerical simulations. Finite element simulations using an elastic solids model were performed to provide this comparison. It was found that this model did well in representing the characteristics of extrusion flow seen in the experiments; an aspect that was not present in the biviscous simulations. The MRI results show that for the range of values studied, both the piston velocity and the contraction ratio have little effect on the characteristics of the flow, including the size and location of the apparent dead zones. It was found that with an appropriate scaling the central, longitudinal velocity follows a master curve. A decreasing contraction ratio, on the other hand, appears to increase the size of the weak velocity region, in contrast with the simulation results.     

Numerical study of a piston-driven unsteady flow in a pipe with sudden expansion

This paper presents results of the numerical study of a piston-driven unsteady flow in a pipe with sudden expansion. The piston closes the larger-diameter pipe and moves between two limiting positions with strong acceleration or deceleration at the beginning and end of each stroke and constant velocity in between. The piston velocity in the exhaust stroke is about four times higher than in the intake stroke. Periodic piston movement in this fashion creates a complex unsteady flow between the piston head and the plane of sudden expansion. The numerical method is implicit and of finite volume type, using a moving grid and a collocated arrangement of variables. Second-order spatial discretization, fine grids and a multigrid solution method were used to ensure high accuracy and good efficiency. Spatial and temporal discretization errors were of the order of 1% and 0.1% respectively. The features of the flow are discussed and the velocity profiles are compared with experimental data, showing good qualitative and quantitative agreement.     

Comparative study of in-cylinder tumble flows in an internal combustion engine using different piston shapes—an insight using particle image velocimetry

This paper deals with the experimental investigations of the in-cylinder tumble flows in a single-cylinder engine with five different piston crown shapes at an engine speed of 1,000 rev/min., during suction and compression strokes under motoring conditions using particle image velocimetry. Two-dimensional in-cylinder tumble flow measurements and analysis are carried out in combustion space on a vertical plane passing through cylinder axis. Ensemble average velocity vectors are used to analyze the tumble flow structure. Tumble ratio and average turbulent kinetic energy are evaluated and used to characterize the tumble flows. From results, it is found that at end of compression, pentroof-offset-bowl piston shows about 41 and 103% improvement in tumble ratio and average turbulent kinetic energy respectively, compared to that of flat piston. The present study will be useful in understanding effect of piston crown shapes on nature of the in-cylinder fluid tumble flows under real engine conditions.     

NUMERICAL MODELING FOR FLUID-STRUCTURE INTERACTION UNDER BLAST AND IMPACT LOADING RESPONSE

建立了基于位移相等条件的流固耦合数值模拟程序框架,通过串联弱耦合方式对爆炸冲击波作用下流固耦合效应进行数值模拟. 其中非线性固体位移场采用基于Lagrange方法描述的时域间断伽辽金有限元方法进行处理. 基于修正弹簧近似的非结构动网格新技术,非定常流场采用格心格式的有限体积方法进行求解. 数值模拟结果表明: 该文所发展的弱耦合分析程序在流固耦合数值模拟过程中具有可靠的计算精度. 同时,程序对于爆炸强间断气动冲击载荷冲击作用下固体变形引起流场的反射和叠加效应, 以及流场变化引起复杂的固体响应具有良好的耦合求解能力.     

基于当地流活塞理论的气动弹性计算方法研究

发展了一种高效、高精度的超音速、高超音速非定常气动力计算 方法------基于定常CFD技术的当地流活塞理论. 运用当地流活塞理论计算非定常 气动力，耦合结构运动方程，实现超音速、高超音速气动弹性的时域模拟. 运用这 种方法计算了一系列非定常气动力算例和颤振算例，并和原始活塞理论、非定 常Euler方程结果作了比较. 由于局部地使用活塞理论假设，这种方法大大地克服 了原始活塞理论对飞行马赫数、翼型厚度和飞行迎角的 限制. 与非定常Euler方程方法相比，当地流活塞理论的效率很高.     

A semi‐implicit scheme for large Eddy simulation of piston engine flow and combustion

A semi‐implicit scheme is presented for large eddy simulation of turbulent reactive flow and combustion in reciprocating piston engines. First, the governing equations in a deforming coordinate system are formulated to accommodate the moving piston. The numerical scheme is made up of a fourth‐order central difference for the diffusion terms in the transport equations and a fifth‐order weighted essentially nonoscillatory (WENO) scheme for the convective terms. A second‐ order Adams–Bashforth scheme is used for time integration. For higher density ratios, it is combined with a predictor–corrector scheme. The numerical scheme is explicit for time integration of the transport equations, except for the continuity equation which is used together with the momentum equation to determine the pressure field and velocity field by using a Poisson equation for the pressure correction field. The scheme is aimed at the simulation of low Mach number flows typically found in piston engines. An efficient multigrid method that can handle high grid aspect ratio is presented for solving the pressure correction equation. The numerical scheme is evaluated on two test engines, a laboratory four‐stroke engine with rectangular‐shaped engine geometry where detailed velocity measurements are available, and a modified truck engine with practical cylinder geometry where lean ethanol/air mixture is combusted under a homogeneous charge compression ignition (HCCI) condition. Copyright © 2012 John Wiley & Sons, Ltd.     

Effect of pressure-dependent slip on flow curve multiplicity

Various microstructural pictures for slip at polymer/solid interfaces lead to relations which have a region where multiple values of slip velocity are predicted for the same shear stress. This leads to the expectation of multivalued flow curves, which has been verified in specific cases by numerous researchers. We study the effect of pressure dependence on flow curve multiplicity using a simple multivalued slip relation to model the phenomena of hysteresis and spurt flow in polymer extrusion. A continuation technique is used to trace out the boundaries of the region of flow curve multiplicity as pressure drop and die length to diameter (L/D) ratio are changed. Results for Newtonian, shear thinning and viscoelastic constitutive equations show that, despite the multivalued nature of the slip model, multiplicity (and thus hysteresis) is absent at high L/D.  For the sake of completeness, we also carry out time-dependent simulations at constant piston speed taking fluid compressibility into account. These simulations show that oscillations in the pressure drop and exit volumetric flow rate result only if the system is operated in the multiplicity region of the steady state flow curve, in agreement with the results of similar simulations by researchers using various multivalued slip models without pressure dependence. The results demonstrate that a multivalued slip model does not guarantee multiplicity in the flow curve for the constant pressure drop operation, nor oscillations for constant piston speed operation. Received: 18 August 1997 Accepted: 30 March 1998     

A Combined Numerical and Experimental Study of the 3D Tumble Structure and Piston Boundary Layer Development During the Intake Stroke of a Gasoline Engine

Due to its positive effect on flame propagation in the case of a well-defined breakdown, the formation of a large-scale tumble motion is an important goal in engine development. Cycle-to-cycle variations (CCV) in the tumble position and strength however lead to a fluctuating tumble breakdown in space and time and therefore to combustion variations, indicated by CCV of the peak pressure. This work aims at a detailed investigation of the large-scale tumble motion and its interaction with the piston boundary layer during the intake stroke in a state-of-the-art gasoline engine. To allow the validation of the flow near the piston surface obtained by simulation, a new measurement technique called “Flying PIV” is applied. A detailed comparison between experimental and simulation results is carried out as well as an analysis of the obtained flow field. The large-scale tumble motion is investigated based on numerical data of multiple highly resolved intake strokes obtained using scale-resolving simulations. A method to detect the tumble center position within a 3D flow field, as an extension of previously developed 2D and 3D algorithms, is presented and applied. It is then used to investigate the phase-averaged tumble structure, its characteristics in terms of angular velocity and the CCV between the individual intake strokes. Finally, an analysis is presented of the piston boundary layer and how it is influenced by the tumble motion during the final phase of the intake stroke.     

LDV measurements of the flow inside the combustion chamber of a 4-valve D.I. diesel engine with axisymmetric piston-bowls

 LDV measurements of the tangential and radial velocity components have been performed to investigate the flow patterns inside a 4-valve single cylinder diesel engine for various engine speeds, swirl numbers and piston-bowl geometries, under conditions similar to those of a production engine. The work focused on the near-wall region of the axisymmetric combustion chambers where the strong swirl-squish and spray-wall interactions take place. The different axial flow stratification observed in both velocity components has confirmed the influence of piston shape and intake configuration on the flow patterns and turbulence levels near compression TDC. Measured values are in reasonable agreement with those provided by a two-zone phenomenological model. Received: 6 November 1995 / Accepted: 12 April 1996     

The unsteady drag on a spherical bubble at large Reynolds numbers

The flow past a spherical bubble undergoing a rectilinear motion in the unsteady flow of an unbounded liquid medium is investigated. The liquid velocity field at infinity is assumed to be uniform and the Reynolds number to be large. The Strouhal number is taken to be of order unity. The velocity distribution is sought by superposition of a perturbation field on the potential flow past the bubble so that the flow field is divided into four regions, i.e. the external flow field where the potential flow holds, the boundary layer, the rear stagnation point region and the wake. The flow in the rear stagnation point region and the wake is assumed to be essentially inertial. The unsteady drag experienced by the bubble is calculated from the mechanical energy balance of the liquid.     

Application of a Synthetic Jet to Control Boundary Layer Separation under Ultra-High-Lift Turbine Pressure Distribution

The transition and separation processes of the boundary layer developing on a flat plate under a prescribed adverse pressure gradient typical of Ultra-High-Lift low-pressure turbine profiles have been investigated, with and without the application of a synthetic jet (zero net mass flow rate jet). A mechanical piston has been adopted to produce an intermittent flow with zero net mass flow rate. The capability of the device to suppress or reduce the large laminar separation bubble occurring under steady inflow condition at low Reynolds numbers has been experimentally investigated by means of hot-wire measurements. Wall static pressure measurements complement the hot-wire time-resolved velocity results. The paper reports the investigations performed for both steady and controlled conditions. The active device is able to control the laminar separation bubble induced at low Reynolds number conditions by the strong adverse pressure gradient. An overall view of the time-dependent evolution of the controlled boundary layer is provided by the phase-locked ensemble averaging technique, triggered at the synthetic jet frequency. The separated flow transition process, which is detected for the uncontrolled condition, is modified by the synthetic jet in different ways during the blowing and suction phases. Overall, the phase-locked velocity distributions show a reduced separated flow region for the whole jet cycle as compared to the uncontrolled condition. The phase-locked distributions of the random unsteadiness allow the identification of vortical structures growing along the shear layer mainly during the blowing phase.     

On the formulation of the problem of strip oscillations and stability in supersonic gas flow

The formulation of the strip flutter problem based on a refined expression for the excess pressure is proposed. Two cases are considered: almost transverse and almost longitudinal flow past the strip. In the first case, an exact expression for the excess pressure, fundamentally different from the formulas of piston theory, which leads to novel, little studied eigenvalue problems is obtained. In the second case, an exact expression is obtained for the flow potential and, for purely longitudinal flow, for the excess pressure also. It is shown that for purely longitudinal flow at high supersonic velocities the critical flutter velocity is equal to the phase velocity of perturbation propagation along the strip, which coincides with the results of piston theory.     

Development of a hybrid finite volume/element solver for incompressible flows

In this paper, we report our development of an implicit hybrid flow solver for the incompressible Navier–Stokes equations. The methodology is based on the pressure correction or projection method. A fractional step approach is used to obtain an intermediate velocity field by solving the original momentum equations with the matrix‐free implicit cell‐centred finite volume method. The Poisson equation derived from the fractional step approach is solved by the node‐based Galerkin finite element method for an auxiliary variable. The auxiliary variable is closely related to the real pressure and is used to update the velocity field and the pressure field. We store the velocity components at cell centres and the auxiliary variable at cell vertices, making the current solver a staggered‐mesh scheme. Numerical examples demonstrate the performance of the resulting hybrid scheme, such as the correct temporal convergence rates for both velocity and pressure, absence of unphysical pressure boundary layer, good convergence in steady‐state simulations and capability in predicting accurate drag, lift and Strouhal number in the flow around a circular cylinder. Copyright © 2007 John Wiley & Sons, Ltd.     

In-cylinder flows of a motored four-stroke engine with flat-crown and slightly concave-crown pistons

The temporal and spatial evolution processes of the in-cylinder flow structures and turbulence intensities in the symmetry and offset planes of a motored four-valve, four-stroke engine during the intake and compression strokes are diagnosed by using a particle image velocimeter. Two pistons of different crown shapes (flat-crown and slightly concave-crown pistons) are studied. The inception, establishment, and evolution of the tumbling vortical flow structures during the intake and compression strokes are clearly depicted. Quantitative strengths of the rotating vortical flow motions are presented by a dimensionless parameter, the tumble ratio, which can represent the mean angular velocity of the vortices in the target plane. The turbulence intensity of the in-cylinder flow is also calculated by using the measured time-varying velocity data. The results show that the flat-crown piston induces higher bulk-averaged tumble ratio and turbulence intensity than the slightly concave-crown piston does because the tumble ratio and turbulence generated by the flat-crown piston in the offset planes during the compression stroke are particularly large. The engine with the flat-crown piston also presents larger torque and power outputs and lower hydrocarbon emission than that with the slightly concave-crown piston. This might be caused by the enhanced combustion in the engine cylinder due to the stronger tumble ratio and turbulence intensity.     

Development of pulsation regimes in one-dimensional unsteady MHD flows with switching off of the electrical conductivity

This paper investigates the gas flow in an electromagnetic field when the conductivity, being a function of the thermodynamic gas parameters, vanishes during the flow (switching off of the conductivity). In the case of steady supersonic flows in an expanding nozzle it was first shown analytically [1] and then confirmed by numerical experiment [2] that stable steady flow is not possible for all the problem parameters (for example, the values of the magnetic field at the exit). Instead of a steady flow a periodic regime is realized when narrow regions of conducting gas with currents flowing through them detach from the conducting region and propagate down the channel. In these papers the conductivity was assumed to be a function of only the temperature, such that for T T* (T) = 0. In [3, 4] the flows of conducting gas in the channels were calculated both with the given dependence of the gas conductivity on the temperature and on the basis of a three-component model by means of the Saha equation. At the same time, the development of periodic regimes in the flow in the nozzle was observed in both cases, but the mechanism of the origin of the current layers was not explained. The self-similar problem of the withdrawal of a nonconducting piston from a half-space occupied by a conducting gas with a magnetic field was investigated in [5] in a linear formulation. At the same time, regions of the problem parameters (the velocity of the piston and the magnetic field on it) were found when, in spite of the self-similar formulation of the problem, there is no self-similar solution. At the same time, regions exist where several solutions are possible. The possibility of the formation of isothermal rarefaction zones with low electrical conductivity when the Joule heating is balanced by the cooling of the gas on expansion (Butler waves) [6] was not taken into account in this paper, since they are unstable with respect to superheating. However, in the case of flow in a nozzle it was shown [2] that precisely the development of instabilities in these zones leads to the formation of the periodic regime. In the present paper the solution of the self-similar problem is constructed in a nonlinear formulation. The reason for the occurrence of regions in which the solution is multiply valued, which is associated with the process of arrival at self-similar boundary conditions, is explained. It is shown that a quasiperiodic regime can arise in the solution, occurring, in particular, in the regions of the problem parameters where there is no self-similar solution.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 115–122, July–August, 1986.     

Detached eddy simulation of cyclic large scale fluctuations in a simplified engine setup

Computational Fluid Dynamics using RANS-based modelling approaches have become an important tool in the internal combustion engine development and optimization process. However, these models cannot resolve cycle to cycle variations, which are an important aspect in the design of new combustion systems. In this study the feasibility of using a Detached Eddy Simulation (DES) SST model, which is a hybrid RANS/LES model, to predict cycle to cycle variations is investigated. In the near wall region or in regions where the grid resolution is not sufficiently fine to resolve smaller structures, the two-equation RANS SST model is used. In the other regions with higher grid resolution an LES model is applied. The case considered is a geometrically simplified engine, for which detailed experimental data for the ensemble averaged and single cycle velocity field are available from Boreé et al. [Boreé, J., Maurel, S., Bazile, R., 2002. Disruption of a compressed vortex, Physics of Fluids 14 (7), 2543–2556]. The fluid flow shows a strong tumbling motion, which is a major characteristic for modern turbo-charged, direct-injection gasoline engines. The general flow structure is analyzed first and the extent of the LES region and the amount of resolved fluctuations are discussed. Multiple consecutive cycles are computed and turbulent statistics of DES SST, URANS and the measured velocity field are compared for different piston positions. Cycle to cycle variations of the velocity field are analyzed for both computation and experiment with a special emphasis on the useability of the DES SST model to predict cyclic variations.