Determination of the effective elastic–plastic response of metal–ceramic composites

The mechanical response of metal–ceramic composites is analysed through a homogenization model accounting for the mechanical behaviour of the constituent materials. In order to achieve this purpose a nonlinear homogenization method based on the phase field approach has been suitably implemented into a numerical code. A prescribed homogenized strain state is applied to a unit volume element of a metal–ceramic composite with proportional loading in which all components of the strain tensor are proportional to one scalar parameter. The mechanical response of the material has been modeled by considering a von Mises plasticity model for the metal phase and a Drucker–Prager associative elastic–plastic material model for the ceramic phase. A two stages plasticity has been obtained in which inelastic strain develops in the metal phase followed by a fully plastic response. A comparison with a finite element model of the stress–strain response of an axisymmetric unit cell has been carried out with the purpose to validate the homogenization based modeling presented in the paper. Plastic parameters of a Drucker–Prager yield surface for the homogenized composite have been calculated at different materials compositions. Associative Drucker–Prager plasticity has been found to be accurate for high ceramic content.     

Experimental and numerical analysis of nonlinear dynamics of rotor–bearing–seal system

Nonlinear dynamics and stability of the rotor–bearing–seal system are investigated both theoretically and experimentally. An experimental rotor–bearing–seal device is designed and corresponding tests are carried out. The experimental rotor system is simplified as the Jeffcott rotor. The nonlinear oil–film forces are obtained under the short bearing theory and Muszynska nonlinear seal force model is used. Numerical method is utilized to solve the nonlinear governing equations. Bifurcation diagrams, waterfall plots, Poincaré maps, spectrum plots and rotor orbits are drawn to analyze various nonlinear phenomena and system unstable processes. Theoretical results from numerical analysis are in good agreement with results from experiments. Conclusions are drawn and prove that this study will contribute to the further understanding of nonlinear dynamics and stability of the rotor system with the fluid-induced forces from oil–film bearings and the seals.     

An example of stick–slip and stick–slip–separation waves

The dynamical problem of a brake-like mechanical system composed of an elastic cylindrical tube with Coulomb's friction in contact with a rigid and rotating cylinder is considered. This model problem enables us to give an example of non-trivial periodic solutions in the form of stick–slip or stick–slip–separation waves propagating on the contact surface. A semi-analytical analysis of stick–slip waves is obtained when the system of governing equations is reduced by condensation to a simpler system involving only the contact displacements. This reduced system, of only one space variable in addition to time, can be solved almost analytically and gives some interesting informations on the existence and the characteristics of stick–slip waves such as the wave numbers on the circumference, stick and slip proportions, wave celerities, tangential and normal forces. It is shown in particular that the stick–slip–separation solutions would occur for small normal pressures or high rotational speeds. Since the analytical discussion becomes cumbersome in this case, a second approach based on numerical analysis by the finite element method is performed. The existence and the characteristics of stick–slip and stick–slip–separation waves are discussed numerically.     

Dynamics of cavitation–structure interaction

Cavitation–structure interaction has become one of the major issues for most engineering applications. The present work reviews recent progress made toward developing experimental and numerical investigation for unsteady turbulent cavitating flow and cavitation–structure interaction. The goal of our overall efforts is to(1) summarize the progress made in the experimental and numerical modeling and approaches for unsteady cavitating flow and cavitation–structure interaction,(2) discuss the global multiphase structures for different cavitation regimes, with special emphasis on the unsteady development of cloud cavitation and corresponding cavitating flow-induced vibrations,with a high-speed visualization system and a structural vibration measurement system, as well as a simultaneous sampling system,(3) improve the understanding of the hydroelastic response in cavitating flows via combined physical and numerical analysis, with particular emphasis on the interaction between unsteady cavitation development and structural deformations. Issues including unsteady cavitating flow structures and cavitation–structure interaction mechanism are discussed.     

Group Properties of 2–Submodels for the Evolutionary Class of Gas–Dynamic Equations

Invariant 2–submodels (submodels with two independent variables) of the evolutionary class are considered for the equations of gas dynamics with an equation of state of general form. Group analysis of these submodels is performed. Allowable operators and transformations of equivalence are indicated, and group classification is performed.     

Detonability of THDCPD-exo–air mixtures

In the frame of industrial risk and propulsive application, the detonability study of JP10–air mixtures was performed. The simulation and measurements of detonation parameters were performed for THDCPD-exo/air mixtures at various initial pressure (1 bar < P ₀ < 3 bar) and equivalence ratio (0.8 < Φ < 1.6) in a heated tube (T ₀ ~ 375 K). Numerical simulations of the detonation were performed with the STANJAN code and a detailed kinetic scheme of the combustion of THDCPD. The experimental study deals with the measurements of detonation velocity and cell size λ. The measured velocity is in a good agreement with the calculated theoretical values. The cell size measurements show a minimum value for Φ ~ 1.2 at every level of initial pressure studied and the calculated induction length L _i corresponds to cell size value with a coefficient k = λ/L _i = 24 at P ₀ = 1 bar. Based on the comparison between the results obtained during this study and those available in the literature on the critical initiation energy E _c, critical tube diameter d _c and deflagration to detonation transition length L _DDT, we can conclude that the detonability of THDCPD–air mixtures corresponds to that of hydrocarbon–air mixtures.

---

This paper is based on the work presented at the 33rd International Pyrotechnics Seminar, IPS 2006, Fort Collins, July 16–21, 2006.     

Yielding and intrinsic plasticity of Ti–Zr–Ni–Cu–Be bulk metallic glass

Bulk metallic glass with composition Ti₄₀Zr₂₅Ni₈Cu₉Be₁₈ exhibits considerably high compressive yield stress, significant plasticity (with a concomitant vein-like fracture morphology) and relatively low density. Yielding and intrinsic plasticity of this alloy are discussed in terms of its thermal and elastic properties. An influence of normal stresses acting on the shear plane is evidenced by: (i) the fracture angle (<45°) and (ii) finite-element simulations of nanoindentation curves, which require the use of a specific yield criterion, sensitive to local normal stresses acting on the shear plane, to properly match the experimental data. The ratio between hardness and compressive yield strength (constraint factor) is analyzed in terms of several models and is best adjusted using a modified expanding cavity model incorporating a pressure-sensitivity index defined by the Drucker–Prager yield criterion. Furthermore, comparative results from compression tests and nanoindentation reveal that deformation also causes strain softening, a phenomenon which is accompanied with the occurrence of serrated plastic flow and results in a so-called indentation size effect (ISE). A new approach to model the ISE of this metallic glass using the free volume concept is presented.     

A combined experimental–numerical investigation of ductile fracture in bending of a class of ferritic–martensitic steel

We present a combined experimental–numerical study on fracture initiation at the convex surface and its propagation during bending of a class of ferritic–martensitic steel. On the experimental side, so-called free bending experiments are conducted on DP1000 steel sheets until fracture, realizing optical and scanning electron microscopy analyses on the post mortem specimens for fracture characterization. A blended Mode I – Mode II fracture pattern, which is driven by cavitation at non-metallic inclusions as well as martensitic islands and resultant softening-based intense strain localization, is observed. Phenomena like crack zig-zagging and crack alternation at the bend apex along the bending axis are introduced and discussed. On the numerical side, based on this physical motivation, the process is simulated in 2D plane strain and 3D, using Gurson’s dilatant plasticity model with a recent shear modification, strain-based void nucleation, and coalescence effects. The effect of certain material parameters (initial porosity, damage at coalescence and failure, shear modification term, etc.), plane strain constraint and mesh size on the localization and the fracture behavior are investigated in detail.     

Existence of Quasipattern Solutions of the Swift–Hohenberg Equation

We consider the steady Swift–Hohenberg partial differential equation, a one-parameter family of PDEs on the plane that models, for example, Rayleigh–Bénard convection. For values of the parameter near its critical value, we look for small solutions, quasiperiodic in all directions of the plane, and which are invariant under rotations of angle ${\pi/q, q \geqq 4}$ . We solve an unusual small divisor problem and prove the existence of solutions for small parameter values, then address their stability with respect to quasi-periodic perturbations.     

A nonlocal phase-field model of Ginzburg–Landau–Korteweg fluids

A thermodynamic model of Korteweg fluids undergoing phase transition and/or phase separation is developed within the framework of weakly nonlocal thermodynamics. Compatibility with second law of thermodynamics is investigated by applying a generalized Liu procedure recently introduced in the literature. Possible forms of the free energy and of the stress tensor, which generalize some earlier ones proposed by several authors in the last decades, are carried out. Owing to the new procedure applied for exploiting the entropy principle, the thermodynamic potentials are allowed to depend on the whole set of variables spanning the state space, including the gradients of the unknown fields, without postulating neither the presence of an energy or entropy extra-flux, nor an additional balance law for microforce.     

Invariant submodels and exact solutions of Khokhlov–Zabolotskaya–Kuznetsov model of nonlinear hydroacoustics with dissipation

We study three-dimensional Khokhlov–Zabolotskaya–Kuznetsov (KZK) model of the nonlinear hydroacoustics with dissipation. This model is described by third order quasilinear partial differential equation of the (KZK). We obtained that the (KZK) equation admits an infinite Lie group of the transformations, depending on the three arbitrary functions. This is due to the fact that in the (KZK) model the main direction of the wave’s propagation is singled out. The submodels of the (KZK) model.are described by the invariant solutions of the (KZK) equation. We studied essentially distinct, not linked by means of the point transformations, invariant solutions of rank 0 and 1 of this equation. Also we considered the invariant solutions of rank 2 and 3. The invariant solutions of rank 0 and 1 are found either explicitly, or their search is reduced to the solution of the nonlinear integro-differential equations. For example, we obtained the invariant solutions that we called by “Ultrasonic knife” and “Ultrasonic destroyer”. The submodel “Ultrasonic knife” have the following property: at each fixed moment of the time in the field of the existence of the solution near a some plane the pressure increases indefinitely and becomes infinite on this plane. The submodel “Ultrasonic destroyer” contains a countable number of “Ultrasonic knives”. The presence of the arbitrary constants in the integro-differential equations, that determine invariant solutions of rank 1 provides a new opportunities for analytical and numerical study of the boundary value problems for the received submodels, and, thus, for the original (KZK) model. With a help of these invariant solutions we researched a propagation of the intensive acoustic waves (one-dimensional, axisymmetric and planar) for which the acoustic pressure, speed and acceleration of its change, or the acoustic pressure , speed and acceleration of its change in the radial direction, or the acoustic pressure, speed and acceleration of its change in the direction of one of the axes are specified at the initial moment of the time at a fixed point. Under the certain additional conditions, we established the existence and the uniqueness of the solutions of boundary value problems, describing these wave processes. Mechanical relevance of the obtained solutions is as follows: (1) these solutions describe nonlinear and diffraction effects in ultrasonic fields of a special kind, (2) these solutions can be used as a test solutions in the numerical calculations performed in studies of ultrasonic fields generated by powerful emitters. Application of the obtained formula generating the new solutions for the found solutions gives families of the solutions containing three arbitrary functions.     

Influence of nanoprecipitates on the creep strength and ductility of a Fe–Ni–Al alloy

Creep studies of a duplex Fe–Ni–Al intermetallic alloy, in two microstructural states, have been carried out at temperatures between 725 and 800 °C (about 0.6 T_m). In the as-cast state, the alloy contains a large volume fraction of nanoprecipitates (50–100 nm) which confer a very high creep strength with a stress exponent of 3 and an activation energy of 280 kJ/mol. The different microstructure obtained in the second state of the alloy, obtained after annealing at 1000 °C for 24 h, leads to a much lower creep strength with a higher stress exponent as well as a large value of the apparent activation energy. While volume diffusion appears to control creep in the as-cast state, both thermal and athermal processes seem to contribute to the different creep rate of material in the annealed state. The latter also exhibits a much larger ductility (12%) relative to that observed in the as-cast material (3%), due to the presence of large numbers of interfaces between the two phases present where strain incompatibilities can be accommodated.     

A computational study of the wing–wing and wing–body interactions of a model insect

The aerodynamic interaction between the contralateral wings and between the body and wings of a model insect are studied, by using the method of numerically solving the Navier-Stokes equations over moving overset grids, under typical hovering and forward flight conditions. Both the interaction between the contralateral wings and the interaction between the body and wings are very weak, e.g. at hovering, changes in aerodynamic forces of a wing due to the present of the other wing are less than 3% and changes in aerodynamic forces of the wings due to presence of the body are less than 2%. The reason for this is as following. During each down- or up-stroke, a wing produces a vortex ring, which induces a relatively large jet-like flow inside the ring but very small flow outside the ring. The vortex rings of the left and right wings are on the two sides of the body. Thus one wing is outside vortex ring of the other wing and the body is outside the vortex rings of the left and right wings, resulting in the weak interactions.     

Dynamics of stochastic Lorenz–Stenflo system

This paper discusses the Lorenz–Stenflo system under the influence of L \(\acute{\hbox {e}}\) vy noise. We find conditions under which the solution to stochastic Lorenz–Stenflo system is exponentially stable. We then investigate the estimation of the global attractive set and stochastic bifurcation behavior of the stochastic Lorenz–Stenflo system. Results show that the jump noise can make the solution stable, the bounds and bifurcation to undergo change under some conditions. Numerical results show the effectiveness and advantage of our methods.     

An example of stick–slip waves

An analytical solution representing a family of stick–slip waves is obtained in a simple example modelling the dynamic behaviour of an elastic cylindrical tube in contact with Coulomb's friction with a rigid and rotating cylinder. This family of waves, representing the non-trivial periodic responses of a continuous system of one space variable, is not classical in the literature.     

Stability of the Front under a Vlasov–Fokker–Planck Dynamics

We consider a kinetic model for a system of two species of particles interacting through a long range repulsive potential and a reservoir at given temperature. The model is described by a set of two coupled Vlasov–Fokker–Plank equations. The important front solution, which represents the phase boundary, is a stationary solution on the real line with given asymptotic values at infinity. We prove the asymptotic stability of the front for small symmetric perturbations.     

Behaviour of the Fokker–Planck–Boltzmann equation near a Maxwellian

We consider the initial value problem for the Fokker–Planck–Boltzmann equation namely, viewed as the Boltzmann equation with an additional diffusion term in velocity space to describe, for instance, the transport in thermal baths of binary elastic collisional particles. The strong solution for initial data near an absolute Maxwellian is proved to exist globally in time and tends asymptotically in the -norm to another time dependent self-similar Maxwellian in large time. The effect of the diffusion in phase space is investigated. It produces a diffusion process in velocity space and results in a heating process on the macroscopic fluid-dynamic observable, accelerating the convergence of solutions to the equilibrium of a self-similar Maxwellian at a faster time-decay rate than the Boltzmann equation. This phenomena is also observed for homogeneous Fokker–Planck–Boltzmann equations, where the time-decay rate in the -norm to the self-similar Maxwellian is proved to be faster than exponential. Moreover, the Fokker–Planck–Boltzmann equation is shown to converge (under an appropriate scaling) strongly to the Boltzmann equation in the process of the zero diffusion limit.     

Lumps and rogue waves of generalized Nizhnik–Novikov–Veselov equation

In this paper, via generalized bilinear forms, we consider the (\(2+1\))-dimensional bilinear p-Sawada–Kotera (SK) equation. We derive analytical rational solutions in terms of positive quadratic functions. Through applying the dependent transformation, we present a class of lump solutions of the (\(2+1\))-dimensional SK equation. Those rationally decaying solutions in all space directions exhibit two kinds of characters, i.e., bright lump wave (one peak and two valleys) and bright–dark lump wave (one peak and one valley). In addition, we also obtain three families of bright–dark lump wave solutions to the nonlinear p-SK equation for \(p=3\).     

Analysis of unsteady gas–liquid flows in a rectangular tank: Comparison of Euler–Eulerian and Euler–Lagrangian simulations

We study the dynamics of gas–liquid flows experimentally and computationally in a rectangular bubble column where the gas source is introduced at the corner. The flow in this reactor is complex and inherently unsteady in nature. The two-dimensional liquid phase velocity field is calculated by an Eulerian approach solving the unsteady Reynolds Averaged Navier Stokes equations. The conservation equations are closed using a two parameter turbulence model. The two-way coupling was accounted for by adding source terms in the conservation equations of the continuous phase to take into account the interaction with the dispersed phase. Bubble tracking is achieved through a Lagrangian approach. Here the equations of motion are solved taking into account the drag, pressure, buoyancy and gravity forces. The time-averaged flows along with the variables which characterize turbulence are analyzed for a wide range of gas flow-rates using Euler–Lagrangian simulations. These simulation predictions are validated with Euler–Eulerian simulations where the gas-phase distribution is captured as a void fraction and PIV experiments. The motion of bubbles induces turbulence in the flow. The applicability of two parameter models for turbulence like the standard k–ε model on time-averaged flow properties is addressed. From the results of the time averaged velocity field, turbulence intensity, turbulent viscosity and gas hold-up profiles, it is concluded that the Euler–Lagrangian model is applicable at lower gas flow-rates. The Euler–Eulerian approach was found to be valid at lower as well as higher gas flow-rates.