Numerical Simulation of R-M Instability

In order to capture shock waves and contact discontinuities in the field and easy to program with parallel computation a new algorithm is developed to solve the N-S equations for simulation of R-M instability problems. The method with group velocity control is used to suppress numerical oscillations, and an adaptive non-uniform mesh is used to get fine resolution. Numerical results for cylindrical shock-cylindrical interface interaction with a shock Mach number Ms=1.2 and Atwood number A=0.818, 0.961, 0.980 (the interior density of the interface/outer density p(1)/p(2) = 10, 50, 100, respectively), and for the planar shock-spherical interface interaction with Ms=1.2 and p(1)/p(2) = 14.28are presented. The effect of Atwood number and multi-mode initial perturbation on the R-M instability are studied. Multi-collisions of the reflected shock with the interface is a main reason of nonlinear development of the interface instability and formation of the spike-bubble structures In simulation with double mode perturbation vortex merging and second instability are found. After second instability the small vortex structures near the interface produced. It is important factor for turbulent mixing.

谈镐生文集

本文集收入了著名力学家、应用数学家、中国科学院院士谈镐生先生在流体力学、稀薄气体动力学和应用数学研究领域的论文和研究报告26篇，谈镐生先生倡导和支持力学基础研究的有关文章和论述21篇，谈先生的学术活动和生活图片多幅，以及谈镐生先生生平，最后附有谈先生生平年表。

 

 

Slender-body hypervelocity boundary-layer instability

With novel application of optical techniques, the slender-body hypervelocity boundary-layer instability is characterized in the previously unexplored regime where thermo-chemical effects are important. Narrowband disturbances (500-3000~kHz) are measured in boundary layers with edge velocities of up to 5~km/s at two points along the generator of a 5 degree half angle cone. Experimental amplification factor spectra are presented. Linear stability and PSE analysis is performed, with fair prediction of the frequency content of the disturbances; however, the analysis over-predicts the amplification of disturbances. The results of this work have two key implications: 1) the acoustic instability is present and may be studied in a large-scale hypervelocity reflected-shock tunnel, and 2) the new data set provides a new basis on which the instability can be studied.

Planar reflection of gaseous detonations

Pipes containing flammable gaseous mixtures may be subjected to internal detonation. When the detonation normally impinges on a closed end, a reflected shock wave is created to bring the flow back to rest. This study built on the work of Karnesky (2010) and examined deformation of thin-walled stainless steel tubes subjected to internal reflected gaseous detonations. A ripple pattern was observed in the tube wall for certain fill pressures, and a criterion was developed that predicted when the ripple pattern would form. A two-dimensional finite element analysis was performed using Johnson-Cook material properties; the pressure loading created by reflected gaseous detonations was accounted for with a previously developed pressure model. The residual plastic strain between experiments and computations was in good agreement.

During the examination of detonation-driven deformation, discrepancies were discovered in our understanding of reflected gaseous detonation behavior. Previous models did not accurately describe the nature of the reflected shock wave, which motivated further experiments in a detonation tube with optical access. Pressure sensors and schlieren images were used to examine reflected shock behavior, and it was determined that the discrepancies were related to the reaction zone thickness extant behind the detonation front. During these experiments reflected shock bifurcation did not appear to occur, but the unfocused visualization system made certainty impossible. This prompted construction of a focused schlieren system that investigated possible shock wave-boundary layer interaction, and heat-flux gauges analyzed the boundary layer behind the detonation front. Using these data with an analytical boundary layer solution, it was determined that the strong thermal boundary layer present behind the detonation front inhibits the development of reflected shock wave bifurcation.

Boundary-layer transition on a slender cone in hypervelocity flow with real gas effects

The laminar to turbulent transition process in boundary layer flows in thermochemical nonequilibrium at high enthalpy is measured and characterized. Experiments are performed in the T5 Hypervelocity Reflected Shock Tunnel at Caltech, using a 1 m length 5-degree half angle axisymmetric cone instrumented with 80 fast-response annular thermocouples, complemented by boundary layer stability computations using the STABL software suite. A new mixing tank is added to the shock tube fill apparatus for premixed freestream gas experiments, and a new cleaning procedure results in more consistent transition measurements. Transition location is nondimensionalized using a scaling with the boundary layer thickness, which is correlated with the acoustic properties of the boundary layer, and compared with parabolized stability equation (PSE) analysis. In these nondimensionalized terms, transition delay with increasing CO₂ concentration is observed: tests in 100% and 50% CO₂, by mass, transition up to 25% and 15% later, respectively, than air experiments. These results are consistent with previous work indicating that CO₂ molecules at elevated temperatures absorb acoustic instabilities in the MHz range, which is the expected frequency of the Mack second-mode instability at these conditions, and also consistent with predictions from PSE analysis. A strong unit Reynolds number effect is observed, which is believed to arise from tunnel noise. N_Tr for air from 5.4 to 13.2 is computed, substantially higher than previously reported for noisy facilities. Time- and spatially-resolved heat transfer traces are used to track the propagation of turbulent spots, and convection rates at 90%, 76%, and 63% of the boundary layer edge velocity, respectively, are observed for the leading edge, centroid, and trailing edge of the spots. A model constructed with these spot propagation parameters is used to infer spot generation rates from measured transition onset to completion distance. Finally, a novel method to control transition location with boundary layer gas injection is investigated. An appropriate porous-metal injector section for the cone is designed and fabricated, and the efficacy of injected CO₂ for delaying transition is gauged at various mass flow rates, and compared with both no injection and chemically inert argon injection cases. While CO₂ injection seems to delay transition, and argon injection seems to promote it, the experimental results are inconclusive and matching computations do not predict a reduction in N factor from any CO₂ injection condition computed.

Modern developments in hypersonic wind tunnels

The development of new methods of producing hypersonic wind-tunnel flows at increasing velocities during the last few decades is reviewed with attention to airbreathing propulsion, hypervelocity aerodynamics and superorbital aerodynamics. The role of chemical reactions in these flows leads to use of a binary scaling simulation parameter, which can be related to the Reynolds number, and which demands that smaller wind tunnels require higher reservoir pressure levels for simulation of flight phenomena. The use of combustion heated vitiated wind tunnels for propulsive research is discussed, as well as the use of reflected shock tunnels for the same purpose. A flight experiment validating shock-tunnel results is described, and relevant developments in shock tunnel instrumentation are outlined. The use of shock tunnels for hypervelocity testing is reviewed, noting the role of driver gas contamination in determining test time, and presenting examples of air dissociation effects on model flows. Extending the hypervelocity testing range into the superorbital regime with useful test times is seen to be possible by use of expansion tube/tunnels with a free piston driver.

Hemodynamics and outcome of septic shock

Septic shock is a common killer in intensive care units (ICU). The most crucial issue concerning the outcome is the early and aggressive start of treatment aimed at normalization of hemodynamics and the early start of antibiotics during the very first hours. The optimal targets of hemodynamic treatment, or impact of hemodynamic treatment on survival after first resuscitation period are less known. The objective of this study was to evaluate different aspects of the hemodynamic pattern in septic shock with special attention to prediction of outcome. In particular components of early treatment and monitoring in the ICU were assessed. A total of 401 patients, 218 with septic shock and 192 with severe sepsis or septic shock were included in the study. The patients were treated in 24 Finnish ICUs during 1999-2005. 295 of the patients were included in the Finnish national epidemiologic Finnsepsis study. We found that the most important hemodynamic variables concerning the outcome were the mean arterial pressures (MAP) and lactate during the first six hours in ICU and the MAP and mixed venous oxygen saturation (SvO2) under 70% during first 48 hours. The MAP levels under 65 mmHg and SvO2 below 70% were the best predictive thresholds. Also the high central venous pressure (CVP) correlated to adverse outcome. We assessed the correlation and agreement of SvO2 and mean central venous oxygen saturation (ScvO2) in septic shock during first day in ICU. The mean SvO2 was below ScvO2 during early sepsis. Bias of difference was 4.2% (95% limits of agreement 8.1% to 16.5%) by Bland-Altman analysis. The difference between saturation values correlated significantly to cardiac index and oxygen delivery. Thus, the ScvO2 can not be used as a substitute of SvO2 in hemodynamic monitoring in ICU. Several biomarkers have been investigated for their ability to help in diagnosis or outcome prediction in sepsis. We assessed the predictive value of N-terminal pro brain natriuretic peptide (NT-proBNP) on mortality in severe sepsis or septic shock. The NT-proBNP levels were significantly higher in hospital nonsurvivors. The NT-proBNP 72 hrs after inclusion was independent predictor of hospital mortality. The acute cardiac load contributed to NTproBNP values at admission, but renal failure was the main confounding factor later. The accuracy of NT-proBNP, however, was not sufficient for clinical decision-making concerning the outcome prediction. The delays in start of treatment are associated to poorer prognosis in sepsis. We assessed how the early treatment guidelines were adopted, and what was the impact of early treatment on mortality in septic shock in Finland. We found that the early treatment was not optimal in Finnish hospitals and this reflected to mortality. A delayed initiation of antimicrobial agents was especially associated with unfavorable outcome.

Detonation-driven-shock wave interactions with perforated plates

The study of detonations and their interactions is vital for the understanding of the high-speed flow physics involved and the ultimate goal of controlling their detrimental effects. However, producing safe and repeatable detonations within the laboratory can be quite challenging, leading to the use of computational studies which ultimately require experimental data for their validation. The objective of this study is to examine the induced flow field from the interaction of a shock front and accompanying products of combustion, produced from the detonation taking place within a non-electrical tube lined with explosive material, with porous plates with varying porosities, 0.7-9.7%. State of the art high-speed schlieren photography alongside high-resolution pressure measurements is used to visualise the induced flow field and examine the attenuation effects which occur at different porosities. The detonation tube is placed at different distances from the plates' surface, 0-30 mm, and the pressure at the rear of the plate is recorded and compared. The results indicate that depending on the level of porosity and the Mach number of the precursor shock front secondary reflected and transmitted shock waves are formed through the coalescence of compression waves. With reduced porosity, the plates act almost as a solid surface, therefore the shock propagates faster along its surface.

The Interaction Between Shock Waves and Foam in a Shock Tube

An experimental study and a numerical simulation were conducted to investigate the mechanical and thermodynamic processes involved in the interaction between shock waves and low density foam. The experiment was done in a stainless shock tube (80mm in inner diameter, 10mm in wall thickness and 5360mm in length). The velocities of the incident and reflected compression waves in the foam were measured by using piezo-ceramic pressure sensors. The end-wall peak pressure behind the reflected wave in the foam was measured by using a crystal piezoelectric sensor. It is suggested that the high end-wall pressure may be caused by a rapid contact between the foam and the end-wall surface. Both open-cell and closed-cell foams with different length and density were tested. Through comparing the numerical and experimental end-wall pressure, the permeability coefficients a and 0 are quantitatively determined.

Shock wave reflection from a wedge in a dusty gas

The present paper contains a detailed study of shock wave reflection from a wedge placed in various suspensions. In past works, the incident shock propagated initially in pure gas and the suspension started only at the leading edge of the deflecting wedge. However, in the present case the entire flow field is filled with a gas-dust suspension and the initial shock wave has steady-state structure relative to the shock front. In former studies the transmitted shock wave starts its propagation into the suspension and is reflected from the wedge at the same time. It is therefore obvious that the two unrelated processes of (2D) reflection and (1D) "transitional" relaxation occur simultaneously. In the present case the suspension behind the incident shock wave has reached steady state (i.e., it is a traveling wave) before the shock reaches the wedge leading edge. The reflection process from the deflecting wedge is studied for different dust mass loadings and different dust-particle diameter. It is shown that when the dust loading is low and the dust particle diameter is small the wave reflection pattern is similar to that observed in a similar pure gas case. In addition, an equilibrium state is reached, behind the evolved waves, very quickly. On the other hand, when the dust loading is relatively high and/or the dust particle diameter is relatively large, the observed reflection wave pattern is very different from that seen in a similar pure gas case. In such cases it takes much longer time to reach an equilibrium state behind the reflecting waves. It is also shown that the dust presence significantly affects the (gas) pressure on the wedge surface. The higher the dust loading is, the higher the pressure on the wedge surface. Suspensions composed of solid particle of different size, but having the same dust mass loading, will approach the same equilibrium pressure. However, it will take longer time to reach an equilibrium state for suspensions having large diameter particles. (C) 2004 Elsevier Ltd. All rights reserved.

Effect of plasma temperature on electrostatic shock generation and ion acceleration by laser

The effect of plasma temperature on electrostatic shock generated by a circularly polarized laser pulse in overdense plasma is studied by particle-in-cell simulation. Ion reflection and transmission in the collisionless electrostatic shock (CES) are investigated analytically. As the initial ion temperature is varied, a distinct transition from the laser-driven piston scenario with all ions being reflected to the CES scenario with partial ion reflection is found. The results show that at low but finite temperatures the ions are much more accelerated than if they were cold.

Control of normal shock wave / turbulent boundary-layer interaction using streamwise slots

The effect of streamwise slots on the interaction of a normal shock wave / turbulent boundary layer has been investigated experimentally at a Mach number of 1.3. The surface pressure distribution for the controlled interaction was found to be significantly smeared, featuring a distinct plateau. This was due to a change in shock structure from a typical unseparated normal shock wave boundary layer interaction to a large bifurcated Lambda type shock pattern. Boundary layer velocity measurements downstream of the slots revealed a strong spanwise variation of boundary layer properties whereas the modified shock structure was relatively twodimensional. Oil flow visualisation indicated that in the presence of slots the boundary layer surface flow was highly three dimensional and confirmed that the effect of slots was mainly due to suction and blowing similar to that for passive control with uniform surface ventilation. Three hole probe measurements confirmed that the boundary layer was three dimensional and that the slots introduced vortical motion into the flowfield. Results indicate that when applied to an aerofoil, the control device has the potential to reduce wave drag while incurring only small viscous penalties. The introduction of streamwise vorticity may also be beneficial to delay trailing edge separation and the device is thought to be capable of postponing buffet onset. © 2001 by A N Smith.

Effects of atrazine on female Wistar rats: Morphological alterations in ovarian follicles and immunocytochemical labeling of 90 kDa heat shock protein

Fertility in female mammals may be affected by a variety of endocrine disrupters present in the environment. Herbicide atrazine is an example of endocrine disrupter employed in agriculture, which disrupts estrous cyclicity in rats. Aiming to characterize morphologically the effect of low and sublethal doses of atrazine on the ovaries of Wistar rats, in an effort to determine the possible intrafollicular target site through which this herbicide acts adult females were submitted to both subacute and subchronic treatments. Additionally, immunocytochemical labeling of 90 kDa heat shock protein (HSP90) was performed in order to evaluate the role played by this protein in the ovary, under stressed conditions induced by herbicide exposure. The results indicated that atrazine induced impaired folliculogenesis, increased follicular atresia and HSP90 depletion in female rats submitted to subacute treatment, while the subchronic treatment with low dose of atrazine could compromise the reproductive capacity reflected by the presence of multioocytic follicle and stress-inducible HSP90. © 2007 Elsevier Ltd. All rights reserved.