Packed beds find extensive application in a wide variety of industries. The objective of the present work is to analyze and evaluate the effects of the wall on structural characteristics, hydrodynamics and heat transfer in packed beds of spheres. As a first attempt, spheres of uniform size are considered. The cylindrical wall of the bed confines the location of the particles thus leading to significant radial variations in void fraction and specific lateral surface area. The two characteristics at any given radial position r* are estimated by defining a concentric cylindrical channel (CCC) of an arbitrary thickness such that its boundaries are equidistant from the cylindrical surface passing through r* and accounting for the solid volumes or lateral surface areas of the segments of spheres (cap, slice, rod and annular ring) contained in the CCC and with centers lying within a distance of a particle radius from r*.The curved boundaries of the sphere segments are rigorously accounted for. The low aspect ratio beds (aspect ratio less than or equal to 2) show three distinct types of behavior. In beds of aspect ratio 2, the void fraction starts from a value of unity at the wall and decreases to a minimum and then increases to unity at the center of the bed. In beds with aspect ratio between l\/¯3/2, there is a continuous decrease in void fraction from unity at the wall to a fairly low value towards the axis and then a slight increase followed by another decrease. The profiles for aspect ratio less than l\/¯3/2 show a continuous decrease from a value of unity at the wall to zero towards the axis. In contrast, beds of high aspect ratio show heavily damped oscillations in the void fraction up to about five particle diameters from the wall and then a constant value. The lateral surface area variations in low aspect ratio beds show a steep fall from a very high value near the wall, and in high aspect ratio beds an oscillatory nature, though not as strong as in the corresponding void fraction profiles. The distribution of flow in packed beds for steady flow of an incompressible Newtonian fluid under isothermal conditions is modeled by using Ergun equation with Brinkman-type correction to account for the viscous effects in the region close to the wall. The confining effect of the wall is incorporated through the radial variations in void fraction and specific lateral surface area. The hydraulic radius in the region next to the wall is modified to take into account the resistance of the wall surface to flow. The resulting model equations with appropriate boundary conditions are solved numerically by collocation technique. The influence of aspect ratio in the range 1.25 to 20.3 and Reynolds number from 0.1 to 1000, the two most important factors affecting the flow behavior, is evaluated. The velocity profiles show a peak in the region close to the wall thus indicating severe channeling effect in this region. The magnitude and location of the peak depend on aspect ratio and Reynolds number. The model predictions agree remarkably with reported experimental data on velocity profiles in a bed of aspect ratio 10.7, and on the effect of Reynolds number on friction factors in beds of low aspect ratio. The radial variations in void fraction, velocity and effective thermal conductivity are incorporated in the two-dimensional pseudo-homogeneous steady-state model to analyze the wall effects on heat transfer in packed beds. Both constant wall temperature and constant wall flux boundary conditions are adopted. The equations are solved numerically using finite difference technique. The radial temperature profiles are seen to be fairly uniform in beds of low aspect ratio thus showing that the often made assumption of complete radial thermal mixing in low aspect ratio beds is valid. Beds of high aspect ratio show strong radial gradients. For constant heat flux condition the slope of the temperature profile remains constant after a small distance from the Inlet thus leading to thermally fully-developed flow. For this condition the heat transfer equations are solved analytically to obtain expressions for Nusselt number and the radial temperature profiles. There is a significant difference in the temperature profiles evaluated in the presence and absence of wall effects. Good agreement is found between the Nusselt numbers obtained from the model and reported experimental data.

<p>A numerical high order difference method is developed for solution of the incompressible Navier-Stokes equations. The solution is determined on a staggered curvilinear grid in two dimensions and by a Fourier expansion in the third dimension. The description in curvilinear body-fitted coordinates is obtained by an orthogonal mapping of the equations to a rectangular grid where space derivatives are determined by compact fourth order approximations. The time derivative is discretized with a second order backward difference method in a semi-implicit scheme, where the nonlinear terms are linearly extrapolated with second order accuracy.</p><p>An approximate block factorization technique is used in an iterative scheme to solve the large linear system resulting from the discretization in each time step. The solver algorithm consists of a combination of outer and inner iterations. An outer iteration step involves the solution of two sub-systems, one for prediction of the velocities and one for solution of the pressure. No boundary conditions for the intermediate variables in the splitting are needed and second order time accurate pressure solutions can be obtained.</p><p>The method has experimentally been validated in earlier studies. Here it is validated for flow past a circular cylinder as an example of a physical test case and the fourth order method is shown to be efficient in terms of grid resolution. The method is applied to external flow past a parabolic body and internal flow in an asymmetric diffuser in order to investigate the performance in two different curvilinear geometries and to give directions for future development of the method. It is concluded that the novel formulation of boundary conditions need further investigation.</p><p>A new iterative solution method for prediction of velocities allows for larger time steps due to less restrictive convergence constraints.</p>

Im Rahmen des Modellierungsseminars wurde die Kopplung einer Strömung mit der Ausbreitung einer reaktiven Substanz im Strömungsgebiet untersucht. Die Strömung wurde dabei durch die inkompressiblen Navier-Stokes-Gleichungen beschrieben. Zusätzlich wurde ein mathematisches Modell für die Ausbreitung der Substanz durch eine Diffusions-Konvektions-Gleichung bestimmt. Beide wurden durch die FEM- Sofware NAVIER berechnet und simuliert.

Thesis (M.S.) -- Mississippi State University. Department of Computational Engineering. / Title from title screen. Includes bibliographical references.

Thesis (Ph. D.)--Mississippi State University. Department of Engineering. / Title from title screen. Includes bibliographical references.

Not available / text

Not available / text

Disertacijos santraukoje apžvelgiami Navjė-Stokso lygčių periodiniai laiko atžvilgiu uždaviniai srityse su cilindriniais išėjimais į begalybę. / In this PhD dissertation summary is considered time periodic Navier-Stokes equations in domains with cylindrical outlets to infinity.

An unstructured explicit, Reynolds averaged Navier-Stokes solver is developed to operate on inviscid flows, laminar flows and turbulent flows and one equation Spalart-Allmaras turbulence modeling is implemented to the solver. A finite volume formulation, which is cell-center based, is used for numerical discretization of Navier-Stokes equations in conservative form. This formulation is combined with one-step, explicit time marching upwind numerical scheme that is the first order accurate in space. Turbulent viscosity is calculated by using one equation Spalart-Allmaras turbulence transport equation. In order to increase the convergence of the solver local time stepping technique is applied. Eight test cases are used to validate the developed solver,for inviscid flows, laminar flows and turbulent flows. All flow regimes are tested on NACA-0012 airfoil. The results of NACA-0012 are compared with the numerical and experimental data.

ABSTRACT SIMULATION OF SURFACE WAVES GENERATED BY A RAPID RISE OF A BLOCK AT THE SEA BOTTOM SENOL, Nalan M.Sc., Department of Civil Engineering, Supervisor: Assoc. Prof. Dr. ismail AYDIN July 2005, 74 Pages A mathematical model is developed for investigating time dependent surface deformations of a hydrostatic water volume, when it is subjected to a sudden partial rise of the sea bottom. In the model, 2-dimensional, compressible, and viscous Navier-Stokes equations are solved by Marker and Cell (MAC) method. Variable mesh size in both horizontal and vertical directions with a staggered grid arrangement is used. Limited compressibility model is utilized for pressure. Various computational tests are done for the selection of computational parameters of the model. It is found that the amplitude of surface waves generated by vertical displacements of the sea bottom depends on size and speed of bottom displacements.