First- and second-order sensitivities of beams with respect to cross-sectional cracks

This paper deals with a variational formulation for the sensitivity problem of beam systems in the context of deformable solids with cracks. Natural frequencies are defined as state variables involved in the energy functional of the system, while the cracks length and position are treated as system parameters. The hierarchical equation system is formed and solved for the first and second derivatives of the natural frequency functions with respect to the cracks length and position. An analytical procedure for calculations of the second-order sensitivities of natural frequencies is proposed for the non-symmetrical equation system operator. Numerical algorithms are worked out and implemented computationally. Analytical and numerical aspects of the problem are discussed in detail through a number of illustrative results.The support of this work by the State Committee for Scientific Research (KBN) under Grant No. 4-050-0148/17-98-00 is gratefully acknowledged.     

Plane Problem of Vibrations of an Elastic Floating Plate under Periodic External Loading

The Wiener–Hopf technique is used to construct an analytical solution of the problem of vibrations of a semiinfinite elastic floating plate under periodic external loading. The solution is obtained in explicit form ignoring draft. The dependences of the amplitudes of surface waves and iceplate deflection on the loading distribution and frequency, ice thickness, and liquid depth are studied numerically. It is established that for some types of acting load, no waves propagate in the plate and liquid and the plate vibrations are standing waves localized near the loading region. An example of such vibrations is given and a condition for the occurrence of localized vibrations is found.     

On the Displacement of NAPL Lenses and Plumes in a Phreatic Aquifer

The movement of an LNAPL lens above a sloping or horizontal water table, and a DNAPL lens above an impermeable surface is discussed. The governing equations are derived, using the vertical equilibrium approach and assuming the water mobility to be much greater than that of the NAPL. Analytical solutions are obtained for onedimensional movement of a lens along a sloping water table. They describe the lens movement with the formation of a jump at the leading front (large-scale approximation), and the distribution of NAPL in the transition zone near the jump (smallscale approximation). A model,describing the movement of a lens, taking into account NAPL retention,is proposed. Approximate onedimensional solutions for the movement of a NAPL lens along sloping or horizontal surfaces under such conditions are presented for this model. Some approximate analytical solutions for twodimensional lens (plume) formation and movement are obtained for the case of a point source at a sloping surface.     

Critical Review of Foamy Oil Flow

Foamy oil flow is a term coined to describe a form of twophase oilgas flow that appears to occur during solution gas drive in some heavy oil reservoirs and does not fit the classical models of twophase flow. Most of the evidence supporting the presence of this unusual flow mechanism is circumstantial and comes from attempts to explain much higher than expected well productivity and primary recovery factors in several heavy oil reservoirs. This paper is a review of the available literature on foamy oil flow in primary production of heavy oils under solution gas drive.The mechanisms operating in solution gas drive in heavy oil reservoirs are briefly discussed. The issues related to supersaturation in oil phase, bubble nucleation, critical gas saturation, and relative permeability are discussed. The possible role of rate processes related to the release of solution gas and the formation of a segregated gas phase is reviewed. The porescale mechanisms involved in creation and propagation of dispersed gas flow are discussed. Several published mathematical models of foamy solution gas drive are reviewed with focus on their limitations.The review shows that the theoretical and experimental investigations of foamy oil flow are still in early stages. Although the occurrence of foamy oil flow has been verified in laboratory experiments, its existence at the reservoir scale has not been confirmed. The theoretical understanding of the mechanisms underlying foamy oil flow remains poor.     

Disruption of Shaped‐Charge Jets by a Pulsed Current

This paper gives results from experimental studies of the effect of electricpulse parameters on the development of current instability and disruption of shapedcharge jets. A simple physical model for the development of current instability and the decrease in the cavern depth in the target is proposed. Notations are introduced and analytical expressions are obtained for the critical current (critical linear current density) and the ideal shape of the current pulse required for the disruption of shapedcharge jets. It is shown that the estimate of the final cavern depths in steel target for the proposed model is in fairly good agreement with both experimental values of cavern depths in targets and with the results of numerical calculations using quasitwodimensional unsteady models of MHD instability and volume disruption.     

Analytical solutions of one-dimensional discrete dynamical systems with chaotic behaviour

This paper studies exact analytical solutions in closed form of the difference equation

A K-ε TWO-EQUATION TURBULENCE MODEL FOR THE SOLID-LIQUID TWO-PHASE FLOWS

A two-equation turbulence model has been dereloped for predicting two-phase flow the two equations describe the conserration of turbulence kinetic energy and dissipation rate of that energy for the incompressible carrier fluid in a two-phase flow The continuity, the momentum, K and ε equations are modeled. In this model,the solid-liquid slip veloeites, the particle-particte interactions and the interactions between two phases are considered,The sandy water pipe turbulent flows are sueeessfuly predicted by this turbulince model.     

Evolution of the Diffusion Mixing Layer of Two Gases upon Interaction with Shock Waves

A mathematical model of mechanics of a twovelocity twotemperature mixture of gases is developed. Based on this model, evolution of the mixing layer of two gases with different densities under the action of shock and compression waves is considered by methods of mathematical simulation in the onedimensional unsteady approximation. In the asymptotic approximation of the full model, a solution of an initialboundary problem is obtained, which describes the formation of a diffusion layer between two gases. Problems of interaction of shock and compression waves with the diffusion layer are solved numerically in the full formulation. It is shown that the layer is compressed as the shock wave traverses it; the magnitude of compression depends on shockwave intensity. As the shock wave passes from the heavy gas to the light gas, the mixing layer becomes overcompressed and expands after shockwave transition. The wave pattern of the flow is described in detail. The calculated evolution of the mixinglayer width is in good agreement with experimental data.     

A correlation for yield stress fluid flow through packed beds

Relatively few correlations are available for non-Newtonian fluid flows through packed beds, even though such fluids are frequently used in industry. In this paper, a correlation is presented for yield stress fluid flow through packed beds. The correlation is developed by introducing the yield stress model in place of the Newtonian model used in deriving Erguns equation. The resulting model has three parameters that are functions of the geometry and roughness of the particle surfaces. Two of the parameters can be deduced in the limit as the yield stress becomes negligible and the model reduces to Erguns equation for Newtonian fluids. The third model parameter is determined from experimental data. The correlation relates a defined friction factor to the dimensionless Reynolds and Hedstrom numbers and can be used to predict pressure drop for flow of a yield stress fluid through a packed bed of spherical particles. Conditions for flow or no-flow are also determined in the correlation. Comparison of model calculations, between a Newtonian and a yield stress fluid for flow penetration into a packed bed of spheres, shows the yield stress fluid initially performs similar to the Newtonian fluid, at large Reynolds numbers. At lower Reynolds numbers the yield stress effect becomes important and the flow rate significantly decreases when compared to the Newtonian fluid.     

The ‘Missing’ Self-Similar Free Convection Boundary-Layer Flow Over a Vertical Permeable Surface in a Porous Medium

The free convection boundary-layer flow of a Darcy–Boussinesq fluid from a vertical permeable plate with an inverse-linear temperature distribution is considered. The outstanding characteristics of this self-similar flow which, according to the usual reduction procedure of pseudo-similarity to full similarity, should not exist at all, are analyzed in detail. Thus it is shown that this flow only exists if a lateral suction with a sufficiently large suction parameter _min = 1.079131 is applied. For the threshold value _min the solution is unique but above it multiple solutions are encountered for every given value of .     

Stability of Steady‐State Shear Jet Flows of an Ideal Fluid with a Free Boundary in an Azimuthal Magnetic Field Against Small Long‐Wave Perturbations

The problem of the linear stability of steadystate axisymmetric shear jet flows of a perfectly conducting inviscid incompressible fluid with a free surface in an azimuthal magnetic field is studied. The necessary and sufficient condition for the stability of these flows against small axisymmetric longwave perturbations of special form is obtained by the direct Lyapunov method. It is shown that if this stability condition is not satisfied, the steadystate flows considered are unstable to arbitrary small axisymmetric longwave perturbations. A priori exponential estimates are obtained for the growth of small perturbations. Examples are given of the steadystate flows and small perturbations imposed on them which evolve in time according to the estimates obtained.     

On the Propagation of Long‐Wave Perturbations in a Two‐Layer Free‐Boundary Rotational Fluid

A mathematical model for the propagation of longwave perturbations in a freeboundary shear flow of an ideal stratified twolayer fluid is considered. The characteristic equation defining the velocity of perturbation propagation in the fluid is obtained and studied. The necessary hyperbolicity conditions for the equations of motion are formulated for flows with a monotonic velocity profile over depth, and the characteristic form of the system is calculated. It is shown that the problem of deriving the sufficient hyperbolicity conditions is equivalent to solving a system of singular integral equations. The limiting cases of weak and strong stratification are studied. For these models, the necessary and sufficient hyperbolicity conditions are formulated, and the equations of motion are reduced to the Riemann integral invariants conserved along the characteristics.     

Influence of leading-edge lateral injection angles on the film cooling effectiveness of a gas turbine blade

Typical film-cooling configuration of a symmetrical turbine blade leading edge is investigated using a three-dimensional finite volume method and a multi-block technique. The computational domain includes the curved blade surface as well as the coolant regions and the plenum. The turbulence is approximated by a two layer k– model. The computations have been performed using the TLV two-layer and the TLVA models. However, the utilization of the TLV and TLVA models has not improved the prediction of the lateral averaged film cooling effectiveness of gas turbine blades when compared with those obtained using wall function strategy.The general features of film cooling such as jet blow-off, high turbulence intensity in the shear layer, and secondary rotating vortices are captured in the present study. Comparison between predicted and experimental results indicates that the trends of the thermal field are well predicted in most cases. In the second part of this study, the influence of lateral injection angle on lateral averaged adiabatic film cooling effectiveness is investigated by varying the lateral injection angle around the experimental value ( = 25°, 30°, 35° and 60° spanwise to the blade surface). It was found that the best coverage and consequently, the maximum film cooling effectiveness are provided by the most extremely inclined injection angle, which is 25° in this investigation.     

Nonaqueous Phase Liquid Dissolution in Porous Media: Current State of Knowledge and Research Needs

Our understanding of nonaqueous phase liquid (NAPL) dissolution in the subsurface environment has been increasing rapidly over the past decade. This knowledge has provided the basis for recent developments in the area of NAPL recovery, including cosolvent and surfactant flushing. Despite these advances toward feasible remediation technologies, there remain a number of unresolved issues to motivate environmental researchers in this area. For example, the lack of an effective NAPLlocation methodology precludes effective deployment of NAPL recovery technologies. The objectives of this paper are to critically review the state of knowledge in the area of stationary NAPL dissolution in porous media and to identify specific research needs. The review first compares NAPL dissolutionbased mass transfer correlations reported for environmental systems with more fundamental results from the literature involving model systems. This comparison suggests that our current understanding of NAPL dissolution in smallscale (on the order of cm) systems is reasonably consistent with fundamental mass transfer theory. The discussion then expands to encompass several issues currently under investigation in NAPL dissolution research, including: characterizing NAPL morphology (i.e. effective size and surface area); multicomponent mixtures; scale-related issues (dispersion, flow by-passing); locating NAPL in the subsurface and enhanced NAPL recovery. Research needs and potential approaches are discussed throughout the paper. This review supports the following conclusions: (1) Our knowledge related to local dissolution and remediation issues is maturing, but should be brought to closure with respect to the link between NAPL emplacement theory (as it impacts NAPL morphology) and NAPL dissolution; (2) The role of nonideal NAPL mixtures, and intra-NAPL mass transfer processes must be clarified; (3) Valid models for quantifying and designing NAPL recovery schemes with chemical additives need to be refined with respect to chemical equilibria, mass transfer and chemical delivery issues; (4) Computational and large-scale experimental studies should begin to address parameter up-scaling issues in support of model application at the field scale; and (5) Inverse modeling efforts aimed at exploiting the previous developments should be expanded to support field-scale characterization of NAPL location and strength as a dissolving source.     

The Experimental Approach to Effective Stress Law of Coal Mass by Effect of Methane

Effective stress law of all kinds of coal samples, including steam coal, fat coal, corking coal, thin coal and anthracite, under pore pressure of gas, is experimentally studied using a newly developed test machine. These samples are taken from Coal Mines in Wuda, Hebi, Yanzhou, Yangquan, Qingshui, and Gujiao in China. The experiment results show that, under pore pressure of gas, the tested coal samples comply with Biots effective stress law, where the Biots coefficient is not a constant, and is bilinear function of volumetric stress () and pore pressure (p), that is, We define four areas according to the numerical feature of , that is, functionless area of pore pressure, normal function area, fracturing function area, and quasi-soil function area. The effective stress law of coal mass introduced by this paper is a constitutive equation in the study of coupled solid and fluid. This has significance in the drainage and outburst of methane in coal seam.     

Methods of Small and Large δ in the Nonlinear Dynamics – A Comparative Analysis

New asymptotic approaches for dynamical systems containing a power nonlinear term x ⁿ are proposed and analyzed. Two natural limiting cases are studied: n 1 + , 1 and n . In the firstcase, the 'small method' (SM)is used and its applicability for dynamical problems with the nonlinearterm sin as well as the usefulness of the SMfor the problem with small denominators are outlined. For n , a new asymptotic approach is proposed(conditionally we call it the 'large method' –LM). Error estimations lead to the followingconclusion: the LM may be used, even for smalln, whereas the SM has a narrow application area. Both of the discussed approaches overlap all values ofthe parameter n.     

On the Quasi‐Steadiness Hypothesis as Applied to Gas Exhaustion from a Receiver

A semiempirical method of determining the stabilization time for a quasisteady mode of gas exhaustion from a receiver after sudden opening of the nozzle and the time evolution of the real flow rate at the stage of the transitional process are considered. The numerical solution of the equations of exhaustion gas dynamics in a twodimensional formulation and the results of model experiments demonstrated that the method can be used to estimate the conditions of applicability of the quasisteadiness hypothesis and to determine the discharge coefficient of the nozzle with controlled accuracy.     

Ši’lnikov Chaos in the Generalized Lorenz Canonical Form of Dynamical Systems

This paper studies the generalized Lorenz canonical form of dynamical systems introduced by elikovský and Chen [International Journal of Bifurcation and Chaos 12(8), 2002, 1789]. It proves the existence of a heteroclinic orbit of the canonical form and the convergence of the corresponding series expansion. The ilnikov criterion along with some technical conditions guarantee that the canonical form has Smale horseshoes and horseshoe chaos. As a consequence, it also proves that both the classical Lorenz system and the Chen system have ilnikov chaos. When the system is changed into another ordinary differential equation through a nonsingular one-parameter linear transformation, the exact range of existence of ilnikov chaos with respect to the parameter can be specified. Numerical simulation verifies the theoretical results and analysis.     

Design of a Beam with Controllable Force Elements

A method for solving the problem of design of an intellectual structure formulated for the pair optimal position of actuators, optimal control of actuators is developed. In the method proposed, physical and logical objects are treated as equivalent.     

Upscaling in Subsurface Transport Using Cluster Statistics of Percolation

Transport/flow problems in soils have been treated in random resistor network representations (RRNs). Two lines of argument can be used to justify such a representation. Solute transport at the pore-space level may probably be treated using a system of linear, first-order differential equations describing inter-pore probability fluxes. This equation is equivalent to a random impedance network representation. Alternatively, Darcys law with spatially variable hydraulic conductivity is equivalent to an RRN. Darcys law for the hydraulic conductivity is applicable at sufficiently low pressure head in saturated soils, but only for steady-state flow in unsaturated soils. The result given here will have two contributions, one of which is universal to any linear conductance problem, i.e., requires only the applicability of Darcys (or Ohms) law. The second contribution depends on the actual distribution of linear conductances appropriate. Although nonlinear effects in RRNs (including changes in resistance values resulting from current, analogous to changes in matric potential resulting from flow) have been treated within the framework of percolation theory, the theoretical development lags the corresponding development of the linear theory, which is, in principle, on a solid foundation. In practice, calculations of the nonlinear conductivity in relatively (compared with soils) well characterized solid-state systems such as amorphous or impure semiconductors, do not agree with each other or with experiment. In semiconductors, however, experiments do at least appear consistent with each other.In the limit of infinite system size the transport properties of a sufficiently inhomogeneous medium are best calculated through application of critical rate analysis with the system resistivity related to the critical (percolating) resistance value, R_c. Here well-known cluster statistics of percolation theory are used to derive the variability, W (R,x) in the smallest maximal resistance, R of a path spanning a volume x³ as well as to find the dependence of the mean value of the conductivity, (x). The functional form of the cluster statistics is a product of a power of cluster size, and a scaling function, either exponential or Gaussian, but which, in either case, cuts off cluster sizes at a finite value for any maximal resistance other than R_c. Either form leads to a maximum in W (R,x) at R=R_c. When the exponential form of the cluster statistics is used, and when individual resistors are exponential functions of random variables (as in stochastic treatments of the unsaturated zone by the McLaughlin group [see Graham and MacLaughlin (1991), or the series of papers by Yeh et al. (1985, 1995), etc.], or as is known for hopping conduction in condensed matter physics), then W (R,x) has a power law decay in R/R_c (or R_c/R, the power being an increasing function of x. If the statistics of the individual resistors are given by power law functions of random variables (as in Poiseiulles Law), then an exponential decay in R for W (R,x) is obtained with decay constant an increasing function of x. Use, instead, of the Gaussian cluster statistics alters the case of power law decay in R to an approximate power, with the value of the power a function of both R and x.