Modeling and simulation of martensitic phase transitions with a triple point

A framework for modeling complex global energy landscapes in a piecewise manner is presented. Specifically, a class of strain-dependent energy functions is derived for the triple point of Zirconia (ZrO₂), where tetragonal, orthorhombic (orthoI) and monoclinic phases are stable. A simple two-dimensional framework is presented to deal with this symmetry breaking. An explicit energy is then fitted to the available elastic moduli of Zirconia in this two-dimensional setting. First, we use the orbit space method to deal with symmetry constraints in an easy way. Second, we introduce a modular (piecewise) approach to reproduce or model elastic moduli, energy barriers and other characteristics independently of each other in a sequence of local steps. This allows for more general results than the classical Landau theory (understood in the sense that the energy is a polynomial of invariant polynomials). The class of functions considered here is strictly larger. Finite-element simulations for the energy constructed here demonstrate the pattern formation in Zirconia at the triple point.     

Wave propagation in an elastic tube: a numerical study

In this paper, a fluid–wall interaction model, called the elastic tube model, is introduced to investigate wave propagation in an elastic tube and the effects of different parameters. The unsteady flow was assumed to be laminar, Newtonian and incompressible, and the vessel wall to be linear-elastic, isotropic and incompressible. A fluid–wall interaction scheme is constructed using a finite element method. The results demonstrate that the elastic tube plays an important role in wave propagation. It is shown that there is a time delay between the velocity waveforms at two different locations and that the peak velocity increases while the low velocity decreases in the elastic tube model, contrary to the rigid tube model where velocity waveforms overlap each other. Compared with the elastic tube model, the increase of the wall thickness makes wave propagation faster and the time delay cannot be observed clearly, however, the velocity amplitude is reduced slightly due to the decrease of the internal radius. The fluid–wall interaction model simulates wave propagation successfully and can be extended to study other mechanical properties considering complicated geometrical and material factors.     

Wave induced thermal boundary layers in a compressible fluid: analysis and numerical simulations

The response of a semi-infinite compressible fluid to a step-wise change in temperature of its boundary is investigated analytically and numerically. Numerical results of the boundary layer structure are compared with Clarke’s analytical solution for a gas with thermal conductivity proportional to temperature. To avoid unwanted numerical dissipation in the numerical analysis, the space-time conservation element and solution element (CESE) method has been adopted to solve the unsteady 1-D Navier-Stokes equations. Good agreement between analytical and numerical results has been found for the development of the thermal boundary layer on a long time scale. Weak shock waves and expansion waves induced by the thermal boundary layer due to its compressibility, are observed in the numerical simulation. Finally, the numerical method has been applied to the reflection of a non-linear expansion wave and to a shock wave from an isothermal wall, thereby illustrating the effect of the boundary layer on the external flow field.     

The elastic fields of misfit cylindrical particles: a dislocation-based numerical approach

Analytical solutions for the displacement and stress fields due to a misfit particle in a host matrix can be difficult to obtain, especially for complex particle geometries. In this work, we present a numerical method for finding such fields in the case of infinitely-long particles. The method is based on discretizing the continuous misfit region between the particle and matrix into local misfit regions consisting of interstitial dislocation loops. The results presented here indicate very good agreement with analytical solutions and better convergence with increasing loop density.     

The elastic response of a cohesive aggregate—a discrete element model with coupled particle interaction

A model is presented for the deformation of a cohesive aggregate of elastic particles that incorporates two important effects of large-sized inter-particle junctions. A finite element model is used to derive a particle response rule, for both normal and tangential relative deformations between pairs of particles. This model agrees with the Hertzian contact theory for small junctions, and is valid for junctions as large as half the nominal particle size. Further, the aggregate model uses elastic superposition to account for the coupled force–displacement response due to the simultaneous displacement of all of the neighbors of each particle in the aggregate. A particle stiffness matrix is developed, relating the forces at each junction to the three displacement degrees of freedom at all of the neighboring-particle junctions. The particle response satisfies force and moment equilibrium, so that the model is properly posed to allow for rigid rotation of the particle without introducing rotational degrees of freedom. A computer-simulated sintering algorithm is used to generate a random particle packing, and the stiffness matrix is derived for each particle. The effective elastic response is then estimated using a mean field or affine displacement calculation, and is also found exactly by a discrete element model, solving for the equilibrium response of the aggregate to uniform-strain boundary conditions. Both the estimate and the exact solution compare favorably with experimental data for the bulk modulus of sintered alumina, whereas Hertzian contact-based models underestimate the modulus significantly. Poisson's ratio is, however, accurately determined only by the full equilibrium discrete element solution, and shown to depend significantly on whether or not rigid particle rotation is permitted in the model. Moreover, this discrete element model is sufficiently robust, so it can be applied to problems involving non-homogeneous deformations in such cohesive aggregates.     

Buckling of nano-fibre reinforced composites: a re-examination of elastic buckling

Elastic buckling of layered/fibre reinforced composites is investigated. Assuming the existence of both shear and transverse modes of failure, the fibre is analysed as a layer embedded in a matrix. Interacting stresses, acting at the interfaces are determined from an exact derived stress field in the matrix. It is shown that buckling can occur only in the shear buckling mode and that the transverse buckling mode is spurious. As opposed to the well known Rosen shear buckling mode solution (predicated on an infinite buckling wavelength), shear buckling is shown to exist under two régimes: buckling of dilute composites with finite wavelengths and buckling of non-dilute composites with infinite wavelengths. Based on the analysis, a model is constructed which defines the fibre concentration at which the transition between the two régimes occurs. The buckling strains are shown to be (approximately) constant for dilute composites and, in the case of very stiff fibres, to have realistic values compatible with elastic behaviour. For the case of non-dilute composites, the strains are found to be in agreement with those given by the Rosen shear buckling solution. Numerical results for the buckling strains and stresses are presented and compared with the Rosen solution. These reveal that the Rosen solution is valid only for the case of non-dilute composites. The investigation demonstrates that elastic buckling may be a dominant failure mechanism of composites consisting of very stiff fibres fabricated in the framework of nano-technology.     

Extraction of cohesive-zone laws from elastic far-fields of a cohesive crack tip: a field projection method

A solution method of an inverse problem is developed to extract cohesive-zone laws from elastic far-fields surrounding a crack-tip cohesive zone. The solution method is named the “field projection method (FPM).” In the process of developing the method a general form of cohesive-crack-tip fields is obtained and used for eigenfunction expansions of the plane elastic field in a complex variable representation. The closing tractions and the separation-gradients at the cohesive zone are expressed in terms of orthogonal polynomial series expansions of the general-form complex functions. The series expansion forms a set of cohesive-crack-tip eigenfunctions, which is complete and orthogonal in the sense of the interaction J-integral in the far field as well as at the cohesive-zone faces. The coefficients of the eigenfunctions in the J-orthogonal representation are extracted directly, using interaction J-integrals in the far field between the physical field of interest and auxiliary probing fields. The path-independence of the interaction J-integral enables us to identify the cohesive-zone variables, i.e. tractions and separations, and thus the cohesive-zone constitutive laws uniquely from the far-field data. A set of numerical algorithms is developed for the inversion method and the results from numerical experiments suggest that the proposed algorithms are well suited for extracting cohesive-zone laws from the far-field data. The set includes methods to find the position and size of a cohesive zone. Further included are discussions on error analysis and stability of the inversion scheme.     

The merging of two co-rotating vortices: a numerical study

The merging of two-dimensional co-rotating vortices is analysed through direct numerical simulations at large Reynolds numbers. It is shown how the Reynolds number affects each of the three phases that characterise this phenomenon. In the first phase, we examine the merging onset and focus on its definition. During the second rapid phase, the contributions of various flow regions upon the dynamics of a vortex are quantitatively studied. These regions are respectively the companion vortex, the filaments and an intermediate zone between vortices and filaments. The third phase is interpreted in terms of an advection diffusion process. Finally the final profile and circulation of the merged vortex is determined: the two thirds of the total circulation of the two initial vortices is contained in the newly formed vortex.     

Translation of a sphere in a rotating viscous fluid: A numerical study

The translation of a sphere moving along the axis of a rotating viscous fluid is studied by the finite difference method at moderate Reynolds (up to R = 500) and Taylor (up to T = 100) numbers. Suppression of the separation is observed with increasing rotation parameter T. The drag coefficient is also presented. It is observed that the drag coefficient is less than that with no rotation in the range 0<N<0·7, where N = 2T/R is the inverse Rossby number. The same phenomenon was observed experimentally by Maxworthy in the range 0<N<0·75±0·03.     

Height of burst explosions: a comparative study of numerical and experimental results

In the current work, we use the Constant Volume model and the numerical method, Regularized Smoothed Particle Hydrodynamics (RSPH) to study propagation and reflection of blast waves from detonations of the high explosives C-4 and TNT. The results from simulations of free-field TNT explosions are compared to previously published data, and good agreement is found. Measurements from height of burst tests performed by the Norwegian Defence Estates Agency are used to compare against numerical simulations. The results for shock time of arrival and the pressure levels are well represented by the numerical results. The results are also found to be in good agreement with results from a commercially available code. The effect of allowing different ratios of specific heat capacities in the explosive products are studied. We also evaluate the effect of changing the charge shape and height of burst on the triple point trajectory.      

Magnetohydrodynamic flow of a Sisko fluid in annular pipe: A numerical study

This paper presents a numerical study for the unsteady flow of a magnetohydrodynamic (MHD) Sisko fluid in annular pipe. The fluid is assumed to be electrically conducting in the presence of a uniform magnetic field. Based on the constitutive relationship of a Sisko fluid, the non‐linear equation governing the flow is first modelled and then numerically solved. The effects of the various parameters especially the power index n, the material parameter of the non‐Newtonian fluid b and the magnetic parameter B on the flow characteristics are explored numerically and presented through several graphs. Moreover, the shear‐thinning and shear‐thickening characteristics of the non‐Newtonian Sisko fluid are investigated and a comparison is also made with the Newtonian fluid. Copyright © 2009 John Wiley & Sons, Ltd.     

A numerical study of a rectangular jet along a shear free boundary: Surface jet

A numerical simulation of a rectangular surface jet is performed at a Reynolds number of ${\mathit{Re}}_j=4400$. The global parameters of the jet e.g. maximum velocity decay, jet surface normal and lateral spread rates, entrainment, jet momentum flux and turbulent momentum flux are in agreement with several other studies reported in the literature. It is shown that the mean velocity and Reynolds stress profiles scale with the maximum local streamwise velocity and jet half width in the surface normal and lateral directions. The current simulation provides balance, explicitly calculated budgets for the turbulence kinetic energy, Reynolds normal and shear stresses. The surface jet develops a thin layer of fast moving fluid in the lateral direction near the surface. This layer is called the ‘surface current’. It has been suggested that the surface current arises due to the Reynolds stress anisotropy in the near surface region. The current study shows that this explanation is incomplete. The turbulence production for the Reynolds stress in the lateral direction is negative, which can drive the mean flow in the lateral direction. The higher level of negative production in the near surface region is responsible for the development of the surface current.     

Wave scattering by a floating porous elastic plate of arbitrary shape: A semi-analytical study

In this paper, a semi-analytical model based on linear potential flow theory and an eigenfunction expansion method is developed to study wave scattering by a porous elastic plate with arbitrary shape floating in water of finite depth. The water domain is divided into the interior and exterior regions, corresponding to the domain beneath the plate and the rest extending towards infinite distance horizontally, respectively. The unknown coefficients in the potential expressions are determined by satisfying the continuity conditions for pressure and velocity at the interface of the two regions, together with the conditions for the motion/force at the edge of the plate, where the Fourier series expansion method is employed to deal with the terms associated with the radius function. A plate with three shapes – circular, cosine and elliptical – and three edge conditions are considered. We find that the largest deflection of the plate with a simply supported edge and a clamped edge is more sensitive to the change in porosity when the porosity is small. The power dissipated by an elliptical plate with its major axis perpendicular to the incident wave direction is the largest among the case studies for the majority of the porosity values tested.     

How a flexible tail improves the power extraction efficiency of a semi-activated flapping foil system: A numerical study

The improvement of power extraction of a semi-activated flapping foil system via the use of a flexible tail is numerically investigated in this work. A NACA0015 airfoil arranged in a two-dimensional laminar flow synchronously executes a forced pitching motion and an induced plunging motion. A flat plate attached to the trailing edge of the foil is utilized to model a tail, and thereby they are considered as a unit for the purpose of power extraction. The tail is either rigid or deformable. At a Reynolds number of 1100 and the position of the pitching axis at third chord, the effects of the mass and flexibility of the tail as well as the frequency of pitching motion on the net power extraction are systematically examined. It is found that compared to the foil with a rigid tail, the efficiency of net power extraction for the foil with a deformable tail can be improved. Based on the numerical analysis, it is indicated that the enhanced power extraction, which is caused by the increased lift force, directly contributes to the net efficiency improvement. Moreover, owing to high enhancement of power extraction, a flexible tail with high flexibility is recommended in the semi-activated flapping foil based power extraction system.     

Effect of incidence angle with wake passing on a film cooled leading edge: A numerical study

This work presents the numerical study of a film‐cooled blade under the influence of wake passing at different incidence angles. The film cooling technology has been proven to be effective to increase the blade life of first turbine stages. However, the leading edge is affected by an high heat transfer rate and cooling this region is difficult. Moreover, separated regions downstream the coolant injection increases the local heat transfer coefficient and can have a detrimental effect in terms of airfoil life. This work analyses how the flow field is affected by the wake passing at different incidence angles (?5, 0, 5) and the impact on heat transfer coefficient. The test case is a linear cascade with two rows of cylindrical holes at the leading edge. Two different holes arrangements are compared in terms of film cooling structures, namely AGTB‐B1 and AGTB‐B2 with 0 and 45^° spanwise inclination. The numerical results show a good agreement with the experiments. A deeper investigation is carried out on AGTB‐B1. The results obtained show that the wake passing and the incidence angle have a strong effect on coolant jets. In particular, there is a significative impact on coolant redistribution near the leading edge. The wake passing has a stronger effect on pressure side, mainly at negative incidence. The predictive approach is based on an U‐RANS in‐house CFD solver using a conventional two‐equations closure. In order to avoid extra turbulence production, critical in the leading edge region, the turbulence model incorporates an extra algebraic equation that enforces a realizability constraint. The unsteady formulation is based on a dual time stepping approach with a sliding plane between the moving bars and the cascade. Copyright © 2009 John Wiley & Sons, Ltd.     

An efficient solution methodology to study the response of a beam on viscoelastic and nonlinear unilateral foundation: Static response

Many engineering materials and foundations such as soils demonstrate nonlinear and viscoelastic behaviour. Yet, it is challenging to develop static and dynamic models of systems that include these materials and are able to predict the behaviour over a wide range of loading conditions. This research is focused on a specific example: a pinned–pinned beam interacting with polyurethane foam foundation. Two cases, when the foundation can react in tension and compression as well as only in compression, are considered. The model developed here is capable of predicting the response to static as well as dynamic forces, whether concentrated or distributed. Galerkin’s method is used to derive modal amplitude equations. In the tensionless foundation case, the contact region changes with beam motion and the estimation of the co-ordinates of the lift-off points is embedded into the solution procedure. An efficient solution technique is proposed that is capable of handling cases where there are multiple contact and non-contact regions. Depending on the loading profiles a high number of modes may need to be included in the solution and to speed up computation time, a convolution method is used to evaluate the integral terms in the model. The adaptability of the solution scheme to complicated loading patterns is demonstrated via examples. The solution approach proposed is applicable to dynamic loadings as well and in these cases the automated treatment of complicated response patterns makes the convolution approach particularly attractive. The influence of various parameters on the static response is discussed.