Elastic wave scattering around cavities in inhomogeneous continua by the BEM

This work examines elastic wave scattering around cavities embedded in a continuum with depth-dependent shear modulus and under conditions of plane strain. A restricted case of inhomogeneity is considered, where the Poisson ratio is fixed at 0.25 and where the density profile also varies, but proportionally to the shear modulus. For this specific case, the wave speeds remain macroscopically constant and it becomes possible to recover the exact Green functions by using an algebraic transformation method. These functions are subsequently used as kernels in a standard 2D boundary element formulation defined in the Laplace transform domain. The final step involves an inverse Laplace transformation, whereby the transient behavior of cavities in the aforementioned inhomogeneous continuum can be recovered. Two basic examples are solved, namely the circular cylindrical cavity under sudden internal explosion and under a pressure wave sweep. In the latter case, it is possible to investigate the effect that the angle of wave incidence has on the displacement and stress that develop along the cavity's perimeter, given the fact that the shear modulus is changing along the vertical direction. These examples serve to illustrate the present approach and to reveal some interesting differences that are observed in transient wave scattering phenomena between homogeneous and continuously inhomogeneous models, where the latter models yield a more realistic representation of geological formations.     

Determination of flutter derivatives via a neural network approach

A neural-network-based method is offered to determine the flutter derivatives of section models under smooth and turbulent flows. The approach uses the observed dynamic responses to train an appropriate neural network. Subsequently, the modal parameters of the model for different mean velocities of wind flow are directly estimated using weight matrices in the neural network. The flutter derivatives can then be determined accurately. The validity of the present method is verified through numerical studies. Finally, the procedure is employed to process experimental data from an inverted-U-type section model, obtained from wind tunnel tests.     

Direct measurement of fluid velocity gradients at a wall by PIV image processing with stereo reconstruction

Velocity gradient is typically estimated in Particle Image Velocimetry (PIV) by differentiating a measured velocity field, which amplifies noise in the measured velocities. If gradients near a boundary are sought, such noise is usually greater than in bulk fluid, because of small tracer displacement, uncertainty in the effective positions of velocity vectors, intense deformation of tracer patterns, and laser reflection. We consider here a modified form of the Particle Image Distortion (PID) method todirectly calculate velocity gradients at a fixed wall, and refer it as “PIV/IG” (“Interface Gradiometry”). Results from synthetic 2D PIV images suggest our method achieves higher SNR and accuracy than velocity differentiation. Also, we have developed a procedure to reconstruct three-dimensional velocity gradient at a fixed wall the two non-zero components from PIV/IG data obtained in stereo views; these equations simplify considerably thanks to the no-slip condition. Experimental data from the bottom wall of turbulent open channel flow appear to suffer from a form of pixel locking. As with standard PIV, this underlines the importance of adequate tracer diameter in the images, sufficient seeding density, and of dynamic range of the camera sensor.     

Die-swell,splashing drop and a numerical technique for solving the Oldroyd B model for axisymmetric free surface flows

This work deals with the development of a numerical method for simulating viscoelastic axisymmetric free surface flow of an Oldroyd B fluid. A novel formulation is developed for the computation of the non-Newtonian extra-stress components on rigid boundaries and on the symmetry axis. The full free surface stress conditions are employed. The resulting governing equations are solved by finite differences on a Marker-and-cell (MAC) type grid. Validation is provided by simulating a pipe flow problem. The classical die-swell problem is solved and swelling ratios are provided. The height of the splash caused by a falling liquid drop for various Reynolds and Weissenberg numbers is then studied, and the height of the splash is shown to diminish with increasing viscoelasticity.     

Error estimation,adaptivity and data transfer in enriched plasticity continua to analysis of shear band localization

In this paper, an adaptive FE analysis is presented based on error estimation, adaptive mesh refinement and data transfer for enriched plasticity continua in the modelling of strain localization. As the classical continuum models suffer from pathological mesh-dependence in the strain softening models, the governing equations are regularized by adding rotational degrees-of-freedom to the conventional degrees-of-freedom. Adaptive strategy using element elongation is applied to compute the distribution of required element size using the estimated error distribution. Once a new mesh is generated, state variables and history-dependent variables are mapped from the old finite element mesh to the new one. In order to transfer the history-dependent variables from the old to new mesh, the values of internal variables available at Gauss point are first projected at nodes of old mesh, then the values of the old nodes are transferred to the nodes of new mesh and finally, the values at Gauss points of new elements are determined with respect to nodal values of the new mesh. Finally, the efficiency of the proposed model and computational algorithms is demonstrated by several numerical examples.     

Boundary value problems for first order impulsive delay differential equations with a parameter

In this paper, by means of the method of upper and lower solutions and monotone iterative technique, the existence of maximal and minimal solutions of the boundary value problems for first order impulsive delay differential equations is established.     

Exact Solutions of Extended Boussinesq Equations

The problem addressed in this paper is the verification of numerical solutions of nonlinear dispersive wave equations such as Boussinesq-like system of equations. A practical verification tool for numerical results is to compare the numerical solution to an exact solution if available. In this work, we derive some exact solitary wave solutions and several invariants of motion for a wide range of Boussinesq-like equations using Maple software. The exact solitary wave solutions can be used to specify initial data for the incident waves in the Boussinesq numerical model and for the verification of the associated computed solution. The invariants of motions can be used as verification tools for the conservation properties of the numerical model.