A compressible near-wall turbulence model for boundary layer calculations

A compressible near-wall two-equation model is derived by relaxing the assumption of dynamical field similarity between compressible and incompressible flows. This requires justifications for extending the incompressible models to compressible flows and the formulation of the turbulent kinetic energy equation in a form similar to its incompressible counterpart. As a result, the compressible dissipation function has to be split into a solenoidal part, which is not sensitive to changes of compressibility indicators, and a dilational part, which is directly affected by these changes. This approach isolates terms with explicit dependence on compressibility so that they can be modeled accordingly. An equation that governs the transport of the solenoidal dissipation rate with additional terms that are explicitly dependent on the compressibility effects is derived similarly. A model with an explicit dependence on the turbulent Mach number is proposed for the dilational dissipation rate. Thus formulated, all near-wall incompressible flow models could be expressed in terms of the solenoidal dissipation rate and straight-forwardly extended to compressible flows. Therefore, the incompressible equations are recovered correctly in the limit of constant density. The two-equation model and the assumption of constant turbulent Prandtl number are used to calculate compressible boundary layers on a flat plate with different wall thermal boundary conditions and free-stream Mach numbers. The calculated results, including the near-wall distributions of turbulence statistics and their limiting behavior, are in good agreement with measurements. In particular, the near-wall asymptotic properties are found to be consistent with incompressible behavior; thus suggesting that turbulent flows in the viscous sublayer are not much affected by compressibility effects

A near-wall two-equation model for compressible turbulent flows

A near-wall two-equation turbulence model of the K - epsilon type is developed for the description of high-speed compressible flows. The Favre-averaged equations of motion are solved in conjunction with modeled transport equations for the turbulent kinetic energy and solenoidal dissipation wherein a variable density extension of the asymptotically consistent near-wall model of So and co-workers is supplemented with new dilatational models. The resulting compressible two-equation model is tested in the supersonic flat plate boundary layer - with an adiabatic wall and with wall cooling - for Mach numbers as large as 10. Direct comparisons of the predictions of the new model with raw experimental data and with results from the K - omega model indicate that it performs well for a wide range of Mach numbers. The surprising finding is that the Morkovin hypothesis, where turbulent dilatational terms are neglected, works well at high Mach numbers, provided that the near wall model is asymptotically consistent. Instances where the model predictions deviate from the experiments appear to be attributable to the assumption of constant turbulent Prandtl number - a deficiency that will be addressed in a future paper

Dissipative chaotic scattering

We show that weak dissipation, typical in realistic situations, can have a metamorphic consequence on nonhyperbolic chaotic scattering in the sense that the physically important particle-decay law is altered, no matter how small the amount of dissipation. As a result, the previous conclusion about the unity of the fractal dimension of the set of singularities in scattering functions, a major claim about nonhyperbolic chaotic scattering, may not be observable.Comment: 4 pages, 2 figures, revte

Atmospheres and Spectra of Strongly Magnetized Neutron Stars -- III. Partially Ionized Hydrogen Models

We construct partially ionized hydrogen atmosphere models for magnetized neutron stars in radiative equilibrium with surface fields B=10^12-5 \times 10^14 G and effective temperatures T_eff \sim a few \times 10^5-10^6 K. These models are based on the latest equation of state and opacity results for magnetized, partially ionized hydrogen plasmas that take into account various magnetic and dense medium effects. The atmospheres directly determine the characteristics of thermal emission from isolated neutron stars. For the models with B=10^12-10^13 G, the spectral features due to neutral atoms lie at extreme UV and very soft X-ray energy bands and therefore are difficult to observe. However, the continuum flux is also different from the fully ionized case, especially at lower energies. For the superstrong field models (B\ga 10^14 G), we show that the vacuum polarization effect not only suppresses the proton cyclotron line as shown previously, but also suppresses spectral features due to bound species; therefore spectral lines or features in thermal radiation are more difficult to observe when the neutron star magnetic field is \ga 10^14 G.Comment: 12 pages, 10 figures; ApJ, accepted (v599: Dec 20, 2003

General-relativistic coupling between orbital motion and internal degrees of freedom for inspiraling binary neutron stars

We analyze the coupling between the internal degrees of freedom of neutron stars in a close binary, and the stars' orbital motion. Our analysis is based on the method of matched asymptotic expansions and is valid to all orders in the strength of internal gravity in each star, but is perturbative in the ``tidal expansion parameter'' (stellar radius)/(orbital separation). At first order in the tidal expansion parameter, we show that the internal structure of each star is unaffected by its companion, in agreement with post-1-Newtonian results of Wiseman (gr-qc/9704018). We also show that relativistic interactions that scale as higher powers of the tidal expansion parameter produce qualitatively similar effects to their Newtonian counterparts: there are corrections to the Newtonian tidal distortion of each star, both of which occur at third order in the tidal expansion parameter, and there are corrections to the Newtonian decrease in central density of each star (Newtonian ``tidal stabilization''), both of which are sixth order in the tidal expansion parameter. There are additional interactions with no Newtonian analogs, but these do not change the central density of each star up to sixth order in the tidal expansion parameter. These results, in combination with previous analyses of Newtonian tidal interactions, indicate that (i) there are no large general-relativistic crushing forces that could cause the stars to collapse to black holes prior to the dynamical orbital instability, and (ii) the conventional wisdom with respect to coalescing binary neutron stars as sources of gravitational-wave bursts is correct: namely, the finite-stellar-size corrections to the gravitational waveform will be unimportant for the purpose of detecting the coalescences.Comment: 22 pages, 2 figures. Replaced 13 July: proof corrected, result unchange

Scaling Theory of Polyelectrolyte Adsorption on Repulsive Charged Surface

We studied polyelectrolyte adsorption on a repulsive charged surface by scaling analysis. At low ionic strength and low surface charge density in which a single polyelectrolyte is able to be adsorbed onto the surface, different regimes in the phase diagram are identified. The possibility of multi-layer structure formed by polyelectrolytes of like charge is also investigated.Comment: 4 pages, 2 figure