G-jitter effects on half floating-zone convection in intermediate-frequency range

The g-jitter influence on thermocapillary convection and critical Marangoni number in a liquid bridge of half-floating rone was discussed in the low frequency range of 0.4 to 1.5 Hz in a previous paper. This paper extended the experiments to the intermediate frequency range of 2 to 18 Hz, which htrs often been recorded as vibration environment of spacecrafts. The experiment was completed on the deck of a vibration machine, which gave a periodical applied acceleration to simulate the effects of g-jitter. The experimental results in the intermediate frequency range are different from that in the low frequency range. The velocity field and the shape of the free surface have periodical fluctuations in response to g-jitter. The amplitude of the periodical varying part of the temperature response decreases obviously with increasing frequency of g-jitter and vanishes almost when the frequency of g-jitter is high enough. The critical Marangoni number is defined to describe the transition from a periodical convection in response to g-jitter to an oscillatory convection due to internal instability, and will increase with increasing g-jitter frequency. According to the spectral analysis, it can be found that the oscillatory part of temperature is a superposition of two harmonic waves if the Marangoni number is larger than a critical value.

Fission of solitary wave in stratified fluid

In this paper, we study the fission of a solitary wave in the stratified fluid with a free surface. It has been discovered that there is no difference between the fissions of the internal solitary waves in odd or even modes, and the effect of the stratification on the fission of a surface solitary wave can almost be neglected

The interaction of two pairs of solitary waves in a two-fluid system

This paper deals with the interaction of solitary waves in a two-fluid system which consistsof two superimposed incompressible inviscid fluids with a free surface and a horizontal rigidbottom. Under the assumption of shallow water wave, we first derive the basic equationssuitable for the model considered, a generalized form of the Boussinesq equations, then usingthe PLK method and the reductive perturbation method, obtain the second-order approximatesolution for the head-on collision between two pairs of interface and surface solitary waves,and give their maximum amplitudes during the collision and the nonuniform phase shiftsafter the collision which lead to the distortion of the wave profiles.

自由曲面的二维自适应测量及测球半径的三维补偿

基于贝塞尔方法提出了针对自由曲面的自适应测量算法，并与传统测量方法进行了对比。在对自适应测量数据的处理中，提出了一种新的基于离散点三角剖分的测球半径三维补偿算法。通过仿真和试验验证，上述算法已成功应用于集成化柔性激光加工系统中。

Nonlinear Dynamic Response of Tension Leg Platform

Tension Leg Platform (TLP) is a typical compliant offshore structure for oil exploitation in deep water. Most of the existing mathematical models for analyzing the dynamic response of TLP are based on explicit or implicit assumptions that displacements (translations and rotations) are small magnitude. Herein a theoretical method for analyzing the nonlinear dynamic behavior of TLP with finite displacement is developed, in which multifold nonlinearities are taken into account, i.e. finite displacement, coupling of the six degrees of freedom, instantaneous position, instantaneous wet surface, free surface effects and viscous drag force. Using this theoretical model, we perform the numerical analysis of dynamic response of a representative TLP. The comparison between the degenerative linear solution of the proposed nonlinear model and the published one shows good agreements. Furthermore, numerical results are presented which illustrate that nonlinearities exert a distinct influence on the dynamic responses of the TLP.

Short-Term Nonlinear Response of Tension Leg Platform in Random Sea Waves

Most of the existing mathematical models for analyzing the dynamic response of TLP are based on explicit or implicit assumptions that motions (translations and rotations) are small magnitude. However, when TLP works in severe adverse conditions, the a priori assumption on small displacements may be inadequate. In such situation, the motions should be regarded as finite magnitude. This paper will study stochastic nonlinear dynamic responses of TLP with finite displacements in random waves. The nonlinearities considered are: large amplitude motions, coupling the six degrees-of-freedom, instantaneous position, instantaneous wet surface, free surface effects and viscous drag force. The nonlinear dynamic responses are calculated by using numerical integration procedure in the time domain. After the time histories of the dynamic responses are obtained, we carry out cycle counting of the stress histories of the tethers with rain-flow counting method to get the stress range distribution.

Floating zone convection

The microgravity research, as a branch of the advanced sciences and a spe- cialized field of high technology, has been made in China since the late 1980's. The research group investigating microgravity fluid physics consisted of our col- leagues and the authors in the Institute of Mechanics of the Chinese Academy of Sciences (CAS), and we pay special attention to the floating zone convection as our first research priority. Now, the research group has expanded and is a part of the National Microgravity Laboratory of the CAS, and the research fields have been extended to include more subjects related to microgravity science. Howev- er, the floating zone convection is still an important topic that greatly holds our research interests.

谈镐生文集

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

 

 

Effect of humidity on microstructure and properties of YBCO film prepared by TFA-MOD method

Epitaxial YBCO superconducting films were deposited on the single crystal LaAlO3 (001) substrate by metal organic deposition method. All YBCO films were fired at 820 degrees C in humidity range of 2.6%-19.7% atmosphere. Microstructure of YBCO thin films was analyzed by means of X-ray diffraction (XRD) and scanning electron microscopy (SEM). Superconducting properties of YBCO films were measured by four-probe method. XRD results showed that the second phase (such as BaF2)and a-axis-oriented grains existed in the films prepared at 2.6% humidity condition; a-axis-oriented grains increased in the film prepared at higher than 4.2% humidity condition; almost pure c-axias-oriented grains existed in the films fired at 4.2% humidity condition. Morphologies of the YBCO films showed that all films had a smooth and crack-free surface. YBCO film prepared at 4.2% humidity condition showed J(c) value of 3.3 MA/cm(2) at 77 K in self-field.

Influence of Rayleigh-Taylor Instability on Liquid Propellant Reorientation in a Low-Gravity Environment

A computational simulation is conducted to investigate the influence of Rayleigh-Taylor instability on liquid propellant reorientation flow dynamics for the tank of CZ-3A launch vehicle series fuel tanks in a low-gravity environment. The volume-of-fluid (VOF) method is used to simulate the free surface flow of gas-liquid. The process of the liquid propellant reorientation started from initially flat and curved interfaces are numerically studied. These two different initial conditions of the gas-liquid interface result in two modes of liquid flow. It is found that the Rayleigh-Taylor instability can be reduced evidently at the initial gas-liquid interface with a high curve during the process of liquid reorientation in a low-gravity environment.

Effect of interfacial heat transfer on the onset of oscillatory convection in liquid bridge

In present study, effect of interfacial heat transfer with ambient gas on the onset of oscillatory convection in a liquid bridge of large Prandtl number on the ground is systematically investigated by the method of linear stability analyses. With both the constant and linear ambient air temperature distributions, the numerical results show that the interfacial heat transfer modifies the free-surface temperature distribution directly and then induces a steeper temperature gradient on the middle part of the free surface, which may destabilize the convection. On the other hand, the interfacial heat transfer restrains the temperature disturbances on the free surface, which may stabilize the convection. The two coupling effects result in a complex dependence of the stability property on the Biot number. Effects of melt free-surface deformation on the critical conditions of the oscillatory convection were also investigated. Moreover, to better understand the mechanism of the instabilities, rates of kinetic energy change and "thermal" energy change of the critical disturbances were investigated (C) 2009 Elsevier Ltd. All rights reserved.

EVALUATION OF FILM-SUBSTRATE ADHESION VIA IMPACT USING COATED BULLETS

A method was devised to evaluate the adhesion between a film and a substrate. A front-end coated bullet is accelerated by a gas gun and hits the substrate of the specimen under test. The impact generates a compressive stress pulse that propagates toward the film. After transmission through the interface, part of the pulse is reflected on the free surface of the film, and tensile stress arises at the film-substrate interface, possibly inducing debonding of the film. This dynamic process was demonstrated analytically and simulated numerically by the finite element method. The results validate the initial concept and lay the foundation for further optimization of this method.

I. Rigid body penetration into brittle materials. II. Phase change effect on shock wave propagation

Part I.

We have developed a technique for measuring the depth time history of rigid body penetration into brittle materials (hard rocks and concretes) under a deceleration of ~ 10⁵ g. The technique includes bar-coded projectile, sabot-projectile separation, detection and recording systems. Because the technique can give very dense data on penetration depth time history, penetration velocity can be deduced. Error analysis shows that the technique has a small intrinsic error of ~ 3-4 % in time during penetration, and 0.3 to 0.7 mm in penetration depth. A series of 4140 steel projectile penetration into G-mixture mortar targets have been conducted using the Caltech 40 mm gas/ powder gun in the velocity range of 100 to 500 m/s.

We report, for the first time, the whole depth-time history of rigid body penetration into brittle materials (the G-mixture mortar) under 10⁵ g deceleration. Based on the experimental results, including penetration depth time history, damage of recovered target and projectile materials and theoretical analysis, we find:

1. Target materials are damaged via compacting in the region in front of a projectile and via brittle radial and lateral crack propagation in the region surrounding the penetration path. The results suggest that expected cracks in front of penetrators may be stopped by a comminuted region that is induced by wave propagation. Aggregate erosion on the projectile lateral surface is < 20% of the final penetration depth. This result suggests that the effect of lateral friction on the penetration process can be ignored.

2. Final penetration depth, P_max, is linearly scaled with initial projectile energy per unit cross-section area, e_s , when targets are intact after impact. Based on the experimental data on the mortar targets, the relation is P_max(mm) 1.15e_s (J/mm² ) + 16.39.

3. Estimation of the energy needed to create an unit penetration volume suggests that the average pressure acting on the target material during penetration is ~ 10 to 20 times higher than the unconfined strength of target materials under quasi-static loading, and 3 to 4 times higher than the possible highest pressure due to friction and material strength and its rate dependence. In addition, the experimental data show that the interaction between cracks and the target free surface significantly affects the penetration process.

4. Based on the fact that the penetration duration, t_max, increases slowly with e_s and does not depend on projectile radius approximately, the dependence of t_max on projectile length is suggested to be described by t_max(μs) = 2.08e_s (J/mm² + 349.0 x m/(πR²), in which m is the projectile mass in grams and R is the projectile radius in mm. The prediction from this relation is in reasonable agreement with the experimental data for different projectile lengths.

5. Deduced penetration velocity time histories suggest that whole penetration history is divided into three stages: (1) An initial stage in which the projectile velocity change is small due to very small contact area between the projectile and target materials; (2) A steady penetration stage in which projectile velocity continues to decrease smoothly; (3) A penetration stop stage in which projectile deceleration jumps up when velocities are close to a critical value of ~ 35 m/s.

6. Deduced averaged deceleration, a, in the steady penetration stage for projectiles with same dimensions is found to be a(g) = 192.4v + 1.89 x 10⁴, where v is initial projectile velocity in m/s. The average pressure acting on target materials during penetration is estimated to be very comparable to shock wave pressure.

7. A similarity of penetration process is found to be described by a relation between normalized penetration depth, P/P_max, and normalized penetration time, t/t_max, as P/P_max = f(t/t_max, where f is a function of t/t_max. After f(t/t_max is determined using experimental data for projectiles with 150 mm length, the penetration depth time history for projectiles with 100 mm length predicted by this relation is in good agreement with experimental data. This similarity also predicts that average deceleration increases with decreasing projectile length, that is verified by the experimental data.

8. Based on the penetration process analysis and the present data, a first principle model for rigid body penetration is suggested. The model incorporates the models for contact area between projectile and target materials, friction coefficient, penetration stop criterion, and normal stress on the projectile surface. The most important assumptions used in the model are: (1) The penetration process can be treated as a series of impact events, therefore, pressure normal to projectile surface is estimated using the Hugoniot relation of target material; (2) The necessary condition for penetration is that the pressure acting on target materials is not lower than the Hugoniot elastic limit; (3) The friction force on projectile lateral surface can be ignored due to cavitation during penetration. All the parameters involved in the model are determined based on independent experimental data. The penetration depth time histories predicted from the model are in good agreement with the experimental data.

9. Based on planar impact and previous quasi-static experimental data, the strain rate dependence of the mortar compressive strength is described by σ_f/σ⁰_f = exp(0.0905(log(έ/έ_0) ^1.14, in the strain rate range of 10^-7/s to 10³/s (σ⁰_f and έ are reference compressive strength and strain rate, respectively). The non-dispersive Hugoniot elastic wave in the G-mixture has an amplitude of ~ 0.14 GPa and a velocity of ~ 4.3 km/s.

Part II.

Stress wave profiles in vitreous GeO₂ were measured using piezoresistance gauges in the pressure range of 5 to 18 GPa under planar plate and spherical projectile impact. Experimental data show that the response of vitreous GeO₂ to planar shock loading can be divided into three stages: (1) A ramp elastic precursor has peak amplitude of 4 GPa and peak particle velocity of 333 m/s. Wave velocity decreases from initial longitudinal elastic wave velocity of 3.5 km/s to 2.9 km/s at 4 GPa; (2) A ramp wave with amplitude of 2.11 GPa follows the precursor when peak loading pressure is 8.4 GPa. Wave velocity drops to the value below bulk wave velocity in this stage; (3) A shock wave achieving final shock state forms when peak pressure is > 6 GPa. The Hugoniot relation is D = 0.917 + 1.711u (km/s) using present data and the data of Jackson and Ahrens [1979] when shock wave pressure is between 6 and 40 GPa for ρ₀ = 3.655 gj cm³ . Based on the present data, the phase change from 4-fold to 6-fold coordination of Ge⁺⁴ with O^-2 in vitreous GeO₂ occurs in the pressure range of 4 to 15 ± 1 GPa under planar shock loading. Comparison of the shock loading data for fused SiO₂ to that on vitreous GeO₂ demonstrates that transformation to the rutile structure in both media are similar. The Hugoniots of vitreous GeO₂ and fused SiO₂ are found to coincide approximately if pressure in fused SiO₂ is scaled by the ratio of fused SiO₂to vitreous GeO₂ density. This result, as well as the same structure, provides the basis for considering vitreous Ge0₂ as an analogous material to fused SiO₂ under shock loading. Experimental results from the spherical projectile impact demonstrate: (1) The supported elastic shock in fused SiO₂ decays less rapidly than a linear elastic wave when elastic wave stress amplitude is higher than 4 GPa. The supported elastic shock in vitreous GeO₂ decays faster than a linear elastic wave; (2) In vitreous GeO₂ , unsupported shock waves decays with peak pressure in the phase transition range (4-15 GPa) with propagation distance, x, as α 1/x^-3.35 , close to the prediction of Chen et al. [1998]. Based on a simple analysis on spherical wave propagation, we find that the different decay rates of a spherical elastic wave in fused SiO₂ and vitreous GeO₂ is predictable on the base of the compressibility variation with stress under one-dimensional strain condition in the two materials.

I. Impact spallation experiments: fracture patterns and spall velocities. II. Craters in carbonate rocks: an electron paramagnetic resonance analysis of shock damage

This work is divided into two independent papers.

PAPER 1.

Spall velocities were measured for nine experimental impacts into San Marcos gabbro targets. Impact velocities ranged from 1 to 6.5 km/sec. Projectiles were iron, aluminum, lead, and basalt of varying sizes. The projectile masses ranged from a 4 g lead bullet to a 0.04 g aluminum sphere. The velocities of fragments were measured from high-speed films taken of the events. The maximum spall velocity observed was 30 m/sec, or 0.56 percent of the 5.4 km/sec impact velocity. The measured velocities were compared to the spall velocities predicted by the spallation model of Melosh (1984). The compatibility between the spallation model for large planetary impacts and the results of these small scale experiments are considered in detail.

The targets were also bisected to observe the pattern of internal fractures. A series of fractures were observed, whose location coincided with the boundary between rock subjected to the peak shock compression and a theoretical "near surface zone" predicted by the spallation model. Thus, between this boundary and the free surface, the target material should receive reduced levels of compressive stress as compared to the more highly shocked region below.

PAPER 2.

Carbonate samples from the nuclear explosion crater, OAK, and a terrestrial impact crater, Meteor Crater, were analyzed for shock damage using electron para- magnetic resonance, EPR. The first series of samples for OAK Crater were obtained from six boreholes within the crater, and the second series were ejecta samples recovered from the crater floor. The degree of shock damage in the carbonate material was assessed by comparing the sample spectra to spectra of Solenhofen limestone, which had been shocked to known pressures.

The results of the OAK borehole analysis have identified a thin zone of highly shocked carbonate material underneath the crater floor. This zone has a maximum depth of approximately 200 ft below sea floor at the ground zero borehole and decreases in depth towards the crater rim. A layer of highly shocked material is also found on the surface in the vicinity of the reference bolehole, located outside the crater. This material could represent a fallout layer. The ejecta samples have experienced a range of shock pressures.

It was also demonstrated that the EPR technique is feasible for the study of terrestrial impact craters formed in carbonate bedrock. The results for the Meteor Crater analysis suggest a slight degree of shock damage present in the β member of the Kaibab Formation exposed in the crater walls.

Part I - The radiation of elastic waves from a spherical cavity in a half space. Part II - Precision determination of focal depths and epicenters of earthquakes.

Part I: The dynamic response of an elastic half space to an explosion in a buried spherical cavity is investigated by two methods. The first is implicit, and the final expressions for the displacements at the free surface are given as a series of spherical wave functions whose coefficients are solutions of an infinite set of linear equations. The second method is based on Schwarz's technique to solve boundary value problems, and leads to an iterative solution, starting with the known expression for the point source in a half space as first term. The iterative series is transformed into a system of two integral equations, and into an equivalent set of linear equations. In this way, a dual interpretation of the physical phenomena is achieved. The systems are treated numerically and the Rayleigh wave part of the displacements is given in the frequency domain. Several comparisons with simpler cases are analyzed to show the effect of the cavity radius-depth ratio on the spectra of the displacements.

Part II: A high speed, large capacity, hypocenter location program has been written for an IBM 7094 computer. Important modifications to the standard method of least squares have been incorporated in it. Among them are a new way to obtain the depth of shocks from the normal equations, and the computation of variable travel times for the local shocks in order to account automatically for crustal variations. The multiregional travel times, largely based upon the investigations of the United States Geological Survey, are confronted with actual traverses to test their validity.

It is shown that several crustal phases provide control enough to obtain good solutions in depth for nuclear explosions, though not all the recording stations are in the region where crustal corrections are considered. The use of the European travel times, to locate the French nuclear explosion of May 1962 in the Sahara, proved to be more adequate than previous work.

A simpler program, with manual crustal corrections, is used to process the Kern County series of aftershocks, and a clearer picture of tectonic mechanism of the White Wolf fault is obtained.

Shocks in the California region are processed automatically and statistical frequency-depth and energy depth curves are discussed in relation to the tectonics of the area.