On reflection of a planar solitary wave at a vertical wall

The reflection of a planar solitary wave at a vertical wall is investigated by solving the Boussinesq equations analytically as well as numerically. The analytical solution is obtained by means of the inner-outer expansions technique, while the numerical solution is based on a finite-difference scheme. The maximum wave amplitude at the wall and the time at which this maximum amplitude is reached are presented. It is also found that the incident wave does not reflect immediately at the wall as predicted by the linear wave theory. Rather, the wave suffers a time delay, called the phase lag, during the reflection process. This phase lag is found to be inversely proportional to the square root of the initial wave amplitude. As the reflected wave eventually propagates away from the wall, it has a phase shift in comparison with that obtained by the linear wave theory. The analytical results obtained in this paper are in good agreement with the numerical results, and they also agree fairly well with the existing experimental data.     

Quasiperiodic and exponential transient phase waves and their bifurcations in a ring of unidirectionally coupled parametric oscillators

Phase waves rotating in a ring of unidirectionally coupled parametric oscillators are studied. The system has a pair of spatially uniform stable periodic solutions with a phase difference and an unstable quasiperiodic traveling phase wave solution. They are generated from the origin through a period doubling bifurcation and the Neimark?CSacker bifurcation, respectively. In transient states, phase waves rotating in a ring are generated, the duration of which increases exponentially with the number of oscillators (exponential transients). A power law distribution of the duration of randomly generated phase waves and the noise-sustained propagation of phase waves are also shown. These properties of transient phase waves are well described with a kinematical equation for the propagation of wave fronts. Further, the traveling phase wave is stabilized through a pitchfork bifurcation and changes into a standing wave through pinning. These bifurcations and exponential transient rotating waves are also shown in an autonomous system with averaging and a coupled map model, and they agree with each other.     

Stability of a stationary steam-water front in a porous medium

The instability of a plane front between two phases of the same fluid (steam and water) in a porous medium is considered. The configuration is taken to be initially stationary with the more dense phase overlying the less dense phase. The frontal region is assumed sharp, so that macroscopic boundary conditions can be utilized. This assumption precludes the existence of dispersion instabilities. The stabilizing influence of phrase transition as well as the implication of different macroscopic pressure boundary conditions on the stability of the front are discussed and illustrated.     

The maximal Lyapunov exponent of a co-dimension two-bifurcation system excited by a bounded noise

In the present paper,the maximal Lyapunov exponent is investigated for a co-dimension two bifurcation system that is on a three-dimensional central manifold and subjected to parametric excitation by a bounded noise.By using a perturbation method,the expressions of the invariant measure of a one-dimensional phase diffusion process are obtained for three cases,in which different forms of the matrix B,that is included in the noise excitation term,are assumed and then,as a result,all the three kinds of singular boundaries for one-dimensional phase diffusion process are analyzed.Via Monte-Carlo simulation,we find that the analytical expressions of the invariant measures meet well the numerical ones.And furthermore,the P-bifurcation behaviors are investigated for the one-dimensional phase diffusion process.Finally,for the three cases of singular boundaries for one-dimensional phase diffusion process,analytical expressions of the maximal Lyapunov exponent are presented for the stochastic bifurcation system.     

Darcy–Brinkman Flow Through a Corrugated Channel

A perturbation analysis is carried out to the second order to give effective equations for Darcy–Brinkman flow through a porous channel with slightly corrugated walls. The flow is either parallel or normal to the corrugations, and the corrugations of the two walls are either in phase or half-period out of phase. The present study is based on the assumptions that the corrugations are periodic sinusoidal waves of small amplitude, and the channel is filled with a sparse porous medium so that the flow can be described by the Darcy–Brinkman model, which approaches the Darcian or Stokes flow limits for small or large permeability of the medium. The Reynolds number is also assumed to be so low that the nonlinear inertia can be ignored. The effects of the corrugations on the flow are examined, quantitatively and qualitatively, as functions of the flow direction, the phase difference, and the wavelength of the corrugations, as well as the permeability of the channel. It is found that the corrugations will have greater effects when it is nearer the Stokes’ flow limit than the Darcian flow limit, and when the wavelength is shorter. For the same wavelength and phase difference, cross flow is more affected than longitudinal flow by the corrugations. Opposite effects can result from 180° out-of-phase corrugations, depending on the flow direction, the wavelength, as well as the permeability.     

Large Eddy Simulation of a Novel Gas-Assisted Coal Combustion Chamber

In this work a recently presented combustion chamber that is specifically designed for the investigation of gas-assisted coal combustion and the validation of models is simulated under reactive conditions for the first time. In the configuration coal combustion is assisted and stabilized by a methane flame. In the course of the investigation, the configuration’s complexity is increased successively. Results of the isothermal single-phase flow are discussed first. Subsequently, reproducibility of the single-phase methane flame by means of the applied modeling approach is evaluated. In a further step, coal particles having the same thermal power as the methane flame are injected into the configuration. Particle histories, the conversion of the coal particles as well as its retroactive effect on the gas phase are investigated. Experimental results based on laser diagnostics are provided for all operating points and used for comparison with numerical results. Gas phase velocity fields for all operating points are available. In order to identify the reaction in the reactive single-phase case planar laser induced fluorescence of the OH-radical (OH-PLIF) was applied. Overall good agreement between numerical and experimental results could be obtained. In the Large Eddy Simulation (LES) a Flamelet Generated Manifold (FGM) based model is utilized. The four-dimensional manifold is spanned by two mixture fractions, a reaction progress variable and the enthalpy on which the gas phase chemistry gets mapped onto. Thereby, the model accounts for both, volatiles reaction and char conversion. Furthermore, finite rate chemistry effects as well as non-adiabatic physics are considered.     

Hamiltonian structures of dynamics of a gyrostat in a gravitational field

The dynamics of a gyrostat in a gravitational field is a fundamental problem in celestial mechanics and space engineering. This paper investigates this problem in the framework of geometric mechanics. Based on the natural symplectic structure, non-canonical Hamiltonian structures of this problem are derived in different sets of coordinates of the phase space. These different coordinates are suitable for different applications. Corresponding Poisson tensors and Casimir functions, which govern the phase flow and phase space structures of the system, are obtained in a differential geometric method. Equations of motion, as well as expressions of the force and torque, are derived in terms of potential derivatives. We uncover the underlying Lie group framework of the problem, and we also provide a systemic approach for equations of motion. By assuming that the gravitational field is axis-symmetrical and central, SO(2) and SO(3) symmetries are introduced into the general problem respectively. Using these symmetries, we carry out two reduction processes and work out the Poisson tensors of the reduced systems. Our results in the central gravitational filed are in consistent with previous results. By these reductions, we show how the symmetry of the problem affects the phase space structures. The tools of geometric mechanics used here provide an access to several powerful techniques, such as the determination of relative equilibria on the reduced system, the energy-Casimir method for determining the stability of equilibria, the variational integrators for greater accuracy in the numerical simulation and the geometric control theory for control problems.     

Axi-symmetric wave propagation in a cylinder coated with a piezoelectric layer

This paper presents the wave propagation in a cylinder coated with a thin piezoelectric layer. The piezoelectric coupling effects are fully modeled in the mechanics model for this piezoelectric coupled cylindrical shell with bending resistance. The decoupled torsional wave velocity and the dispersion curves for the two- mode shell model are obtained theoretically. The cut-off frequency and phase velocities at limit wave number are also derived. The numerical simulations are conducted to present the results of wave propagation in this cylindrical shell and as well as to compare the results by the current bending theory and the membrane shell theory. From the comparisons, the results display obvious deference of wave propagations in terms of dispersion characteristics by different shell theories when thicker piezoelectric layer are used and when higher wave number is considered. The results of this paper can serve as a reference for future study on wave propagation in coupled structures as well as in the design of smart structures incorporating piezoelectric materials.     

Internal waves around a disturbance in a fluid with arbitrary stratification and background shear flow

A three-dimensional ray theory is presented for calculating the phase configuration of internal waves around a moving disturbance in a flow with arbitrary stratification and background shear. The theory is applied to two-dimensional stratified shear flows which have been produced in the laboratory and good agreement is shown between the theoretical and experimental phase configurations. Good agreement is also shown when caustics and critical levels are present. This paper includes the wave systems arising from a combined translation and oscillatory motion of a source and it shows how distinct systems of waves arise from each, type of motion under different background conditions. This paper shows that for two-dimensional steady wave systems the critical level is at a well defined height which is independent of wavenumber but in three dimensions the critical height can in general vary with wavenumber.     

Nonlinear resonance of a subsatellite on a short constant tether

The nonlinear resonant behavior of a subsatellite on a short constant tether during station-keeping phase is investigated in this paper. The nonlinear dynamic equations of in-plane motion of the system are derived based on Kane’s method first. Then an approach of multiple scales expressed in matrix form is employed in solving the simplified nonlinear system of cubic nonlinearity near its local equilibrium position. Analysis shows that there exists a three-to-one resonance in such a nonlinear system with two degrees of freedom. Afterward, the approximate solution up to third order determined analytically by the Weierstrass elliptic function is obtained and the comparison between the approximate and numerical solutions presented as well. The results show that the approximate solution is coincide well with the numerical solution of original system. The nonlinear resonance of the subsatellite on short tether exhibits coexistent quasiperiodic motions or a quasiperiodic oscillation near local equilibrium position.     

Obtaining phase averaged turbulence properties in the near wake of a circular cylinder at high Reynolds number using POD

The flow past a circular cylinder at high Reynolds number is studied by means of PIV, 3C-PIV and Time-Resolved PIV techniques. One of the goals of this study was to allow comparisons with numerical simulations on a domain identical to that of the experiment. For this reason, the cylinder was placed in a confined environment, with a high blockage and a low aspect ratio, thereby allowing computations on a mesh of reasonable size, and avoiding “infinite conditions”. This paper deals with the decomposition of the flow in a coherent and random parts. To this aim, phase averaged quantities were first obtained using the wall pressure signal on the cylinder as a trigger signal. This was achieved using both conditional sampling and LSE with similar results. This decomposition is then analysed using the Time Resolved PIV measurements, as well as by comparison of the contributions of the organised and turbulent fluctuations to the time-independent Reynolds stress tensor with those estimated from velocity spectra by interpolation and integration of the continuous part. In agreement with other studies, it is found that the contribution of the turbulent motion is overestimated as a result of the occurence of phase jitter between the trigger and velocity signal. A POD analysis was then performed to extract the coherent motion and to compare this decomposition with that obtained by phase averaging. Similarly to the phase averaging, the POD allows the decomposition of the time-independent stress tensor as the sum of two contributions corresponding to the first N modes, and the rest of the modes. This decomposition is then analysed by comparing these contributions to those obtained from the velocity spectra, according to the value N chosen. It is found that these contributions are strongly dependent on N, and the contribution of the first modes greatly overestimate the coherent motion if N is too large. In order to obtain a good decomposition of the flow in coherent and random parts, the difficulty in this case lies in the choice of the modes. Finally, the POD coefficients of the first two modes are used instead of the pressure signal to determine the phase of the vortex shedding, and the phase averaging is reconsidered. It is found that the phase averaged vortices are less smeared by the averaging process, the turbulent stresses better follow the evolution of the vortices, and the contributions of both coherent and turbulent fluctuations are found to agree well with those evaluated from the velocity spectra. This enhancement is obtained because the phase angle is determined directly from the velocity fields to be averaged, thereby reducing the phase-jitter effect. A comparison with a detached eddy simulation is also briefly shown and demonstrates the high level of agreement obtainable between simulation and experiment, as well as confirming the enhancement of the phase averaging using this procedure.     

Buckling of a circular plate made of a shape memory alloy due to a reverse thermoelastic martensite transformation

We solve the axisymmetric buckling problem for a circular plate made of a shape memory alloy undergoing reverse martensite transformation under the action of a compressing load, which occurs after the direct martensite transformation under the action of a generally different (extending or compressing) load. The problem was solved without any simplifying assumptions concerning the transverse dimension of the supplementary phase transition region related to buckling. The mathematical problem was reduced to a nonlinear eigenvalue problem. An algorithm for solving this problem was proposed. It was shown that the critical buckling load under the reverse transition, which is obtained by taking into account the evolution of the phase strains, can be many times lower than the same quantity obtained under the assumption that the material behavior is elastic even for the least (martensite) values of the elastic moduli. The critical buckling force decreases with increasing modulus of the load applied at the preliminary stage of direct transition and weakly depends on whether this load was extending or compressing. In shape memory alloys (SMA), mutually related processes of strain and direct (from the austenitic into the martensite phase) or reverse thermoelastic phase transitions may occur. The direct transition occurs under cooling and (or) an increase in stresses and is accompanied by a significant decrease (nearly by a factor of three in titan nickelide) of the Young modulus. If the direct transition occurs under the action of stresses with nonzero deviator, then it is accompanied by accumulation of macroscopic phase strains, whose intensity may reach 8%. Under the reverse transition, which occurs under heating and (or) unloading, the moduli increase and the accumulated strain is removed. For plates compressed in their plane, in the case of uniform temperature distribution over the thickness, one can separate trivial processes under which the strained plate remains plane and the phase ratio has a uniform distribution over the thickness. For sufficiently high compressing loads, the trivial process of uniform compression may become unstable in the sense that, for small perturbations of the plate deflection, temperature, the phase ratio, or the load, the difference between the corresponding perturbed process and the unperturbed process may be significant. The results of several experiments concerning the buckling of SMA elements are given in [1, 2], and the statement and solution of the corresponding boundary value problems can be found in [3–11]. The experimental studies [2] and several analytic solutions obtained for the Shanley column [3, 4], rods [5–7], rectangular plates under direct [8] and reverse [9] transitions showed that the processes of thermoelastic phase transitions can significantly (by several times) decrease the critical buckling loads compared with their elastic values calculated for the less rigid martensite state of the material. Moreover, buckling does not occur in the one-phase martensite state in which the elastic moduli are minimal but in the two-phase state in which the values of the volume fractions of the austenitic and martensite phase are approximately equal to each other. This fact is most astonishing for buckling, studied in the present paper, under the reverse transition in which the Young modulus increases approximately half as much from the beginning of the phase transition to the moment of buckling. In [3–9] and in the present paper, the static buckling criterion is used. Following this criterion, the critical load is defined to be the load such that a nontrivial solution of the corresponding quasistatic problem is possible under the action of this load. If, in the problems of stability of rods and SMA plates, small perturbations of the external load are added to small perturbations of the deflection (the critical force is independent of the amplitude of the latter), then the critical forces vary depending on the value of perturbations of the external load [5, 8, 9]. Thus, in the case of small perturbations of the load, the problem of stability of SMA elements becomes indeterminate. The solution of the stability problem for SMA elements also depends on whether the small perturbations of the phase ratio and the phase strain tensor are taken into account. According to this, the problem of stability of SMA elements can be solved in the framework of several statements (concepts, hypotheses) which differ in the set of quantities whose perturbations are admissible (taken into account) in the process of solving the problem. The variety of these statements applied to the problem of buckling of SMA elements under direct martensite transformation is briefly described in [4, 5]. But, in the problem of buckling under the reverse transformation, some of these statements must be changed. The main question which we should answer when solving the problem of stability of SMA elements is whether small perturbations of the phase ratio (the volume fraction of the martensite phase q) are taken into account, because an appropriate choice significantly varies the results of solving the stability problem. If, under the transition to the adjacent form of equilibrium, the phase ratio of all points of the body is assumed to remain the same, then we deal with the “fixed phase atio” concept. The opposite approach can be classified as the “supplementary phase transition” concept (which occurs under the transition to the adjacent form of equilibrium). It should be noted that, since SMA have temperature hysteresis, the phase ratio in SMA can endure only one-sided small variations. But if we deal with buckling under the inverse transformation, then the variation in the volume fraction of the martensite phase cannot be positive. The phase ratio is not an independent variable, like loads or temperature, but, due to the constitutive relations, its variations occur together with the temperature variations and, in the framework of connected models for a majority of SMA, together with variations in the actual stresses. Therefore, the presence or absence of variations in q is determined by the presence or absence of variations in the temperature, deflection, and load, as well as by the system of constitutive relations used in this particular problem. In the framework of unconnected models which do not take the influence of actual stresses on the phase ratio into account, the “fixed phase ratio” concept corresponds to the case of absence of temperature variations. The variations in the phase ratio may also be absent in connected models in the case of specially chosen values of variations in the temperature and (or) in the external load, as well as in the case of SMA of CuMn type, for which the influence of the actual stresses on the phase compound is absent or negligible. In the framework of the “fixed phase ratio” hypothesis, the stability problem for SMA elements has a solution coinciding in form with the solution of the corresponding elastic problem, with the elastic moduli replaced by the corresponding functions of the phase ratio. In the framework of the supplementary phase transition” concept, the result of solving the stability problem essentially depends on whether the small perturbations of the external loads are taken into account in the process of solving the problem. The point is that, when solving the problem in the connected setting, the supplementary phase transition region occupies, in general, not the entire cross-section of the plate but only part of it, and the location of the boundary of this region depends on the existence and the value of these small perturbations. More precisely, the existence of arbitrarily small perturbations of the actual load can result in finite changes of the configuration of the supplementary phase transition region and hence in finite change of the critical values of the load. Here we must distinguish the “fixed load” hypothesis where no perturbations of the external loads are admitted and the “variable load” hypothesis in the opposite case. The conditions that there no variations in the external loads imply additional equations for determining the boundary of the supplementary phase transition region. If the “supplementary phase transition” concept and the “fixed load” concept are used together, then the solution of the stability problem of SMA is uniquely determined in the same sense as the solution of the elastic stability problem under the static approach. In the framework of the “variable load” concept, the result of solving the stability problem for SMA ceases to be unique. But one can find the upper and lower bounds for the critical forces which correspond to the cases of total absence of the supplementary phase transition: the upper bound corresponds to the critical load coinciding with that determined in the framework of the “fixed phase ratio” concept, and the lower bound corresponds to the case where the entire cross-section of the plate experiences the supplementary phase transition. The first version does not need any additional name, and the second version can be called as the "all-round supplementary phase transition" hypothesis. In the present paper, the above concepts are illustrated by examples of solving problems about axisymmetric buckling of a circular freely supported or rigidly fixed plate experiencing reverse martensite transformation under the action of an external force uniformly distributed over the contour. We find analytic solutions in the framework of all the above-listed statements except for the case of free support in the “fixed load” concept, for which we obtain a numerical solution.     

Observation of Stress Field Around an Oscillating Crack Tip in a Quenched Thin Glass Plate

The variation of stress field around an oscillating crack tip in a quenched thin glass plate is observed using instantaneous phase-stepping photoelasticity. The successive images around the propagating crack are recorded by a CCD camera that is equipped with a pixelated micro-retarder array. Then, the phase maps of the principal stress difference and the principal direction are easily obtained even though the photoelastic fringes cannot be visualized. The path of the crack growth as well as the stress intensity factors and the crack tip constraint are obtained from these phase distributions. Results show that the mode I stress intensity factor and the crack tip constraint vary remarkably with the crack growth. In addition, the results show that the mode-II stress intensity factor exists even though the crack propagates smoothly.     

Measurements of concentration per size class in a dense polydispersed jet using planar laser-induced fluorescence and phase Doppler techniques

In dense two-phase flows, it is well known that phase Doppler anemometry is not well suited for the measurement of concentration and mass flux. Laser diagnostics based on fluorescence can provide the dispersed phase concentration but without discrimination between size classes. We present a new method of coupling the two techniques, in order to extract the local value of concentration and flux per size class. The method is applied to an axisymmetric turbulent jet, laden with polydispersed droplets 1–90 μm. Droplet concentration profiles are obtained in the development zone (x/d ₀ < 20) of the dense jet and are used to study droplet dispersion. The results are then introduced into the momentum transport equations to analyze the influence of droplets on the carrier phase. We show that the local decrease of the rate of variation of mean momentum with mass loading is due both to an increase in interfacial transfer rate and to a decrease in turbulent diffusion effects. Received: 20 November 2000 / Accepted: 3 April 2001     

Modeling the acidizing of a horizontal well opening a carbon-bearing oil reservoir

The model of two-phase seepage flow occurring when oil is displaced by an acid solution in the neighborhood of a horizontal well is considered with account for flow in the well bore. The effect of acid activation and acid fracturing is studied. The influence of the phase permeability variations accompanying the rock break-up on the treatment result is estimated. The nonuniformity of the inflow distribution along the horizontal bore is shown to increase as a result of the acidizing. Numerical calculations are carried out for actual situations.     

Pressure Dependent Properties of a Compressible Polymer

This paper deals with the investigations of a porous carbon black-filled rubber, tested with regard to its pressure and tension behaviour. In the tension range only uniaxial tests are performed while in the pressure range uniaxial as well as hydrostatic tests are performed. The uniaxial experiments are carried out in a custom-made uniaxial device and the hydrostatic tests in a pressure chamber which is specially developed for this application. The construction and use of the pressure chamber is clearly described in this paper. All experiments are related to the basic elasticity of the material. The viscoelastic behaviour is completely disregarded at this point. Not only the experiments are discussed, also the modelling of the material is looked at. The tested cellular rubber is composed of an incompressible solid phase and a compressible gas phase. For that reason a so-called structural compressibility is observed. The compressible behaviour of cellular rubber is an important property. So the main focus of the paper is on the pressure tests and the simulation of these. The existing material models for rubber like materials only deal with incompressible rubber structures. To represent the compressible behaviour, the Theory of Porous Media is used. The constitutive model is based on a polynomial approach for an incompressible material. This is complemented by a volumetric expansion term with a point of compaction to model the structural compressibility.     

X-ray based measurements of the local solid phase content in a three-phase flow of a bubble column: statistical significance

Information regarding the flow properties in bubble columns with a three-phase flow is of great interest for research into its performance, as well as for the validation of computer models. In an earlier paper, we proposed an X-ray based particle tracking velocimetry (PTV), called XPTV. This method allows not only the measurement of the velocities of both solid and fluid phases in a three-phase flow, but also the assessment of the time-averaged local solid phase content. In this paper, we are concerned with the statistical significance of the obtained measurements.