Airfoils &deg / apping in pitch and plunge are studied, and the &deg / apping motion parameters are op- timized to maximize thrust generation and the e&plusmn / ciency of the thrust generation. Unsteady viscous &deg / ow&macr / elds over &deg / apping airfoils are computed on overset grids using a Navier-Stokes solver. Computations are performed in parallel using Parallel Virtual Machine library routines in a computer cluster. A single &deg / apping airfoil and dual airfoils &deg / apping in a biplane con- &macr / guration are considered. A gradient based optimization algorithm is employed. The thrust production and the e&plusmn / ciency of the thrust production are optimized with respect to &deg / apping parameters / the plunging and pitching amplitudes, the &deg / apping frequency, and the phase shift between the pitch and plunge motions. It is observed that thrust generation of &deg / apping airfoils strongly depends on the phase shift and high thrust values may be obtained at the expense of reduced e&plusmn / ciency. For a high e&plusmn / ciency in thrust generation, the e&reg / ective angle of attack of the airfoil is reduced and large scale vortex formations at the leading edge are prevented. At a &macr / xed reduced &deg / apping frequency of 1, a single &deg / apping airfoil in pitch and plunge motion produces the maximum average thrust coe&plusmn / cient of 1:41 at the plunge amplitude of 1:60, the pitch amplitude of 23:5o, and the phase shift of 103:4o whereas the maximum e&plusmn / ciency of 67:5% is obtained at the plunge amplitude of 0:83, the pitch amplitude of 35:5o and the phase shift of 86:5o.

High explosive charges when detonated ensue in a flow field characterized by several physical phenomena that include blast wave propagation, hydrodynamic instabilities, real gas effects, fluid mixing and afterburn effects. Solid metal particles are often added to explosives to augment the total impulsive loading, either through direct bombardment if inert, or through afterburn energy release if reactive. These multiphase explosive charges, termed as heterogeneous explosives, are of interest from a scientific perspective as they involve the confluence and interplay of various additional physical phenomena such as shock-particle interaction, particle dispersion, ignition, and inter-phase mass, momentum and energy transfer. In the current research effort, chemical explosions in multiphase environments are investigated using a robust, state-of-the-art Eulerian-gas, Lagrangian-solid methodology that can handle both the dense and dilute particle regimes. Explosions into ambient air as well as into aluminum particle clouds are investigated, and hydrodynamic instabilities such as Rayleigh- Taylor and Richtmyer-Meshkov result in a mixing layer where the detonation products mix with the air and afterburn. The particles in the ambient cloud, when present, are observed to pick up significant amounts of momentum and heat from the gas, and thereafter disperse, ignite and burn. The amount of mixing and afterburn are observed to be independent of particle size, but dependent on the particle mass loading and cloud dimensions. Due to fast response times, small particles are observed to cluster as they interact with the vortex rings in the mixing layer, which leads to their preferential ignition/ combustion. The total deliverable impulsive loading from heterogeneous explosive charges containing inert steel particles is estimated for a suite of operating parameters and compared, and it is demonstrated that heterogeneous explosive charges deliver a higher near-field impulse than homogeneous explosive charges containing the same mass of the high explosive. Furthermore, particles are observed to introduce significant amounts of hydrodynamic instabilities in the mixing layer, resulting in augmented fluctuation intensities and fireball size, and different growth rates for heterogeneous explosions compared to homogeneous explosions. For aluminized explosions, the particles are observed to burn in two regimes, and the average particle velocities at late times are observed to be independent of the initial solid volume fraction in the explosive charge. Overall, this thesis provides useful insights on the role played by solid particles in chemical explosions.

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.

Incompressible flows -- flows in which variations in the density of a fluid are negligible -- arise in a wide variety of applications, from hydraulics to aerodynamics. The incompressible Navier-Stokes equations which govern such flows are also of fundamental physical and mathematical interest. They are believed to hold the key to understanding turbulent phenomena; precise conditions for the existence and uniqueness of solutions remain unknown -- and establishing such conditions is the subject of one of the Clay Mathematics Institute's Millennium Prize Problems. Typical solutions of incompressible flow problems involve both fine- and large-scale phenomena, so that a uniform finite element mesh of sufficient granularity will at best be wasteful of computational resources, and at worst be infeasible because of resource limitations. Thus adaptive mesh refinements are required. In industry, the adaptivity schemes used are ad hoc, requiring a domain expert to predict features of the solution. A badly chosen mesh may cause the code to take considerably longer to converge, or fail to converge altogether. Typically, the Navier-Stokes solve will be just one component in an optimization loop, which means that any failure requiring human intervention is costly. Therefore, I pursue technological foundations for a solver of the incompressible Navier-Stokes equations that provides robust adaptivity starting with a coarse mesh. By robust, I mean both that the solver always converges to a solution in predictable time, and that the adaptive scheme is independent of the problem -- no special expertise is required for adaptivity. The cornerstone of my approach is the discontinuous Petrov-Galerkin (DPG) finite element methodology developed by Leszek Demkowicz and Jay Gopalakrishnan. For a large class of problems, DPG can be shown to converge at optimal rates. DPG also provides an accurate mechanism for measuring the error, and this can be used to drive adaptive mesh refinements. Several approximations to Navier-Stokes are of interest, and I study each of these in turn, culminating in the study of the steady 2D incompressible Navier-Stokes equations. The Stokes equations can be obtained by neglecting the convective term; these are accurate for "creeping" viscous flows. The Oseen equations replace the convective term, which is nonlinear, with a linear approximation. The steady-state incompressible Navier-Stokes equations approximate the transient equations by neglecting time variations. Crucial to this work is Camellia, a toolbox I developed for solving DPG problems which uses the Trilinos numerical libraries. Camellia supports 2D meshes of triangles and quads of variable polynomial order, allows simple specification of variational forms, supports h- and p-refinements, and distributes the computation of the stiffness matrix, among other features. The central contribution of this dissertation is design and development of mathematical techniques and software, based on the DPG method, for solving the 2D incompressible Navier-Stokes equations in the laminar regime (Reynolds numbers up to about 1000). Along the way, I investigate approximations to these equations -- the Stokes equations and the Oseen equations -- followed by the steady-state Navier-Stokes equations. / 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.