A Particle Numerical Model for Wall Film Dynamics in Port-Injected Engines

To help predict hydrocarbon emissions during cold-start conditions the authors are developing a numerical model for the dynamics and vaporization of the liquid wall films formed in port-injected spark-ignition engines and incorporating this model in the KIVA-3 code for complex geometries. This paper summarizes the current status of the project and presents illustrative example calculations. The dynamics of the wall film is influenced by interactions with the impinging spray, the wall, and the gas flow near the wall. The spray influences the film through mass, tangential momentum, and energy addition. The wall affects the film through the no-slip boundary condition and heat transfer. The gas alters film dynamics through tangential stresses and heat and mass transfer in the gas boundary layers above the films. New wall functions are given to predict transport in the boundary layers above the vaporizing films. It is assumed the films are sufficiently thin that film flow is laminar and that liquid inertial forces are negligible. Because liquid Prandtl numbers are typically about then, unsteady heating of the film should be important and is accounted for by the model. The thin film approximation breaks down near sharp corners, where an inertial separation criterion is used. A particle numerical method is used for the wall film. This has the advantages of compatibility with the KIVA-3 spray model and of very accurate calculation of convective transport of the film. The authors have incorporated the wall film model into KIVA-3, and the resulting combined model can be used to simulate the coupled port and cylinder flows in modern spark-ignition engines. They give examples by comparing computed fuel distributions with closed- and open-valve injection during the intake and compression strokes of a generic two-valve engine

An integrated approach to modelling the fluid-structure interaction of a collapsible tube

The well known collapsible tube experiment was conducted to obtain flow, pressure and materials property data for steady state conditions. These were then used as the boundary conditions for a fully coupled fluid-structure interaction (FSI) model using a propriety computer code, LS-DYNA. The shape profiles for the tube were also recorded. In order to obtain similar collapse modes to the experiment, it was necessary to model the tube flat, and then inflate it into a circular profile, leaving residual stresses in the walls. The profile shape then agreed well with the experimental ones. Two departures from the physical properties were required to reduce computer time to an acceptable level. One of these was the lowering of the speed of sound by two orders of magnitude which, due to the low velocities involved, still left the mach number below 0.2. The other was to increase the thickness of the tube to prevent the numerical collapse of elements. A compensation for this was made by lowering the Young's modulus for the tube material. Overall the results are qualitatively good. They give an indication of the power of the current FSI algorithms and the need to combine experiment and computer models in order to maximise the information that can be extracted both in terms of quantity and quality

Hydrodynamical analysis of symmetric nucleus-nucleus collisions at CERN/SPS energies

We present a coherent theoretical study of ultrarelativistic heavy-ion data obtained at the CERN/SPS by the NA35/NA49 Collaborations using 3+1-dimensional relativistic hydrodynamics. We find excellent agreement with the rapidity spectra of negative hadrons and protons and with the correlation measurements in two experiments: $S+S$ at 200 $AGeV$ and $Pb+Pb$ at 160 $AGeV$ (preliminary results). Within our model this implies that for $Pb+Pb$ ($S+S$) a quark-gluon-plasma of initial volume 174 $fm^3$ (24 $fm^3$) with a lifetime 3.4 $fm/c$ (1.5 $fm/c$) was formed. It is found that the Bose-Einstein correlation measurements do not determine the maximal effective radii of the hadron sources because of the large contributions from resonance decay at small momenta. Also within this study we present an NA49 acceptance corrected two-pion Bose-Einstein correlation function in the invariant variable, $Q_{inv}$.Comment: 21 pages, 11 Postscript figures (1 File, 775654 Bytes, has to be requested for submission via e.mail from [email protected]

Thermal photons as a measure for the rapidity dependence of the temperature

The rapidity distribution of thermal photons produced in Pb+Pb collisions at CERN-SPS energies is calculated within scaling and three-fluid hydrodynamics. It is shown that these scenarios lead to very different rapidity spectra. A measurement of the rapidity dependence of photon radiation can give cleaner insight into the reaction dynamics than pion spectra, especially into the rapidity dependence of the temperature.Comment: 3 Figure

K-TIF: a two-fluid computer program for downcomer flow dynamics. [PWR]

The K-TIF computer program has been developed for numerical solution of the time-varying dynamics of steam and water in a pressurized water reactor downcomer. The current status of physical and mathematical modeling is presented in detail. The report also contains a complete description of the numerical solution technique, a full description and listing of the computer program, instructions for its use, with a sample printout for a specific test problem. A series of calculations, performed with no change in the modeling parameters, shows consistent agreement with the experimental trends over a wide range of conditions, which gives confidence to the calculations as a basis for investigating the complicated physics of steam-water flows in the downcomer

SALE: a simplified ALE computer program for fluid flow at all speeds

A simplified numerical fluid-dynamics computing technique is presented for calculating two-dimensional fluid flows at all speeds. It combines an implicit treatment of the pressure equation similar to that in the Implicit Continuous-fluid Eulerian (ICE) technique with the grid rezoning philosophy of the Arbitrary Lagrangian-Eulerian (ALE) method. As a result, it can handle flow speeds from supersonic to the incompressible limit in a grid that may be moved with the fluid in typical Lagrangian fashion, or held fixed in an Eulerian manner, or moved in some arbitrary way to give a continuous rezoning capability. The report describes the combined (ICEd-ALE) technique in the framework of the SALE (Simplified ALE) computer program, for which a general flow diagram and complete FORTRAN listing are included. A set of sample problems show how to use or modify the basic code for a variety of applications. Numerical listings are provided for a sample problem run with the SALE program

Freeze-out in hydrodynamical models

We study the effects of strict conservation laws and the problem of negative contributions to final momentum distribution during the freeze-out through 3-dimensional hypersurfaces with spacelike normal. We study some suggested solutions for this problem, and demonstrate it in one example