Computation of three-dimensional shock wave and boundary-layer interactions

Computations of the impingement of an oblique shock wave on a cylinder and a supersonic flow past a blunt fin mounted on a plate are used to study three dimensional shock wave and boundary layer interaction. In the impingement case, the problem of imposing a planar impinging shock as an outer boundary condition is discussed and the details of particle traces in windward and leeward symmetry planes and near the body surface are presented. In the blunt fin case, differences between two dimensional and three dimensional separation are discussed, and the existence of an unique high speed, low pressure region under the separated spiral vortex core is demonstrated. The accessibility of three dimensional separation is discussed

Development of relaxation turbulence models

Relaxation turbulence models have been intensively studied. The complete time dependent mass averaged Navier-Stokes equations have been solved for flow into a two dimensional compression corner. A new numerical scheme has been incorporated into the developed computed code with an attendant order of magnitude reduction in computation time. Computed solutions are compared with experimental measurements of Law for supersonic flow. Details of the relaxation process have been studied; several different relaxation models, including different relaxation processes and varying relaxation length, are tested and compared. Then a parametric study has been conducted in which both Reynolds number and wedge angle are varied. To assess effects of Reynolds number and wedge angle, the parametric study includes the comparison of computed separation location and upstream extent of pressure rise; numerical results are also compared with the measurements of surface pressure, skin friction and mean velocity field

Dilepton/Photon Production in Heavy Ion Collisions and the QCD Phase Transition

We calculate electromagnetic production from highly excited hadronic matter created in heavy ion collisions. The rates include the usual lowest order processes in quark-gluon plasma plus the usual reactions in the hadronic phase, related with rho,a_1 mesons. The space-time integration is done using a hydrodynamical model. We have found that dilepton mass spectrum agrees with results of other previous works, but *disagrees* with the CERES dilepton data. In order to explain these data, some "unconventional" production mechanism need to be incorporated: we discuss especially the notion of modified rho,a_1 masses, which indeed may explain the data. Other suggestions (e.g. the longer-lived fireball) to increase the production of low mass dileptons seem to be insufficient for the task. The results for direct photon production are below the current WA80 experimental bounds, for all variants considered.Comment: 34 pages LaTeX, with encapsulated ps figures. WWW home at http://looney.physics.sunysb.edu/~daffy/phd.html One comment and one reference added to appendi

A heater made from graphite composite material for potential deicing application

A surface heater was developed using a graphite fiber-epoxy composite as the heating element. This heater can be thin, highly electrically and thermally conductive, and can conform to an irregular surface. Therefore it may be used in an aircraft's thermal deicing system to quickly and uniformly heat the aircraft surface. One-ply of unidirectional graphite fiber-epoxy composite was laminated between two plies of fiber glass-epoxy composite, with nickel foil contacting the end portions of the composite and partly exposed beyond the composites for electrical contact. The model heater used brominated P-100 fibers from Amoco. The fiber's electrical resistivity, thermal conductivity and density were 50 micro ohms per centimeter, 270 W/m-K and 2.30 gm/cubic cm, respectively. The electricity was found to penetrate through the composite in the transverse direction to make an acceptably low foil-composite contact resistance. When conducting current, the heater temperature increase reached 50 percent of the steady state value within 20 sec. There was no overheating at the ends of the heater provided there was no water corrosion. If the foil-composite bonding failed during storage, liquid water exposure was found to oxidize the foil. Such bonding failure may be avoided if perforated nickel foil is used, so that the composite plies can bond to each other through the perforated holes and therefore lock the foil in place

Calculation of a separated turbulent boundary layer

The properties of a Navier-Stokes solution of a shock-separated turbulent flow over a flat wall are investigated. Refinements of an algebraic relaxation turbulence model previously shown to be of value for the simulation of separated flows are presented. A simplified analysis applicable near an adiabatic wall is developed and used to help verify the accuracy of the numerical solution. Features of the time-dependent response of a turbulent boundary layer to shock impingement are presented

Geometrical interpretation and graphical solution to minimum energy discrete-data control Scientific report no. 9

Linear PAM regulator minimum energy design by geometric interpretation and graphical solutio

Mining and Incentive Concession Contracts

This paper studies the design of a mining concession contract as a multi-period autoselection problem where production is the depletion of a non renewable resource. As compared to symmetric information, we show that overproduction (resp. underproduction) is optimal in the initial phase (resp. terminal phase ) of the resource extraction program. Also, asymmetric information lengthens the contract duration but reduces the scarcity rent. Finally, when there are several agents competing for contract bid, we show that optimal auctioning could be used to award the concession, assigning the lowest cost agent to carry out the extraction.ADVERSE SELECTION; EXHAUSTIBILITY; OVERPRODUCTION

Numerical study of large-eddy breakup and its effect on the drag characteristics of boundary layers

The break-up of a field of eddies by a flat-plate obstacle embedded in a boundary layer is studied using numerical solutions to the two-dimensional Navier-Stokes equations. The flow is taken to be incompressible and unsteady. The flow field is initiated from rest. A train of eddies of predetermined size and strength are swept into the computational domain upstream of the plate. The undisturbed velocity profile is given by the Blasius solution. The disturbance vorticity generated at the plate and wall, plus that introduced with the eddies, mix with the background vorticity and is transported throughout the entire flow. All quantities are scaled by the plate length, the unidsturbed free-stream velocity, and the fluid kinematic viscosity. The Reynolds number is 1000, the Blasius boundary layer thickness is 2.0, and the plate is positioned a distance of 1.0 above the wall. The computational domain is four units high and sixteen units long