Implementation and efficiency analysis of parallel computation using OpenACC: a case study using flow field simulations

The Open Accelerator (OpenACC) application programming interface is a relatively new parallel computing standard. In this paper, particle-based flow field simulations are examined as a case study of OpenACC parallel computation. The parallel conversion process of the OpenACC standard is explained, and further, the performance of the flow field parallel model is analysed using different directive configurations and grid schemes. With careful implementation and optimisation of the data transportation in the parallel algorithm, a speedup factor of 18.26× is possible. In contrast, a speedup factor of just 11.77× was achieved with the conventional Open Multi-Processing (OpenMP) parallel mode on a 20-kernel computer. These results demonstrate that optimised feature settings greatly influence the degree of speedup, and models involving larger numbers of calculations exhibit greater efficiency and higher speedup factors. In addition, the OpenACC parallel mode is found to have good portability, making it easy to implement parallel computation from the original serial model.     

Flexural-wave-generation using a phononic crystal with a piezoelectric defect

This paper proposes a method to amplify the performance of a flexural-wave-generation system by utilizing the energy-localization characteristics of a phononic crystal(PnC) with a piezoelectric defect and an analytical approach that accelerates the predictions of such wave-generation performance. The proposed analytical model is based on the Euler-Bernoulli beam theory. The proposed analytical approach, inspired by the transfer matrix and S-parameter methods, is used to perform band-structure an...     

Fractional Fringe Order Determination using Digital Moire Fringe Shifting

Ｆｒａｃｔｉｏｎａｌ　Ｆｒｉｎｇｅ　Ｏｒｄｅｒ　Ｄｅｔｅｒｍｉｎａｔｉｏｎ　ｕｓｉｎｇ　Ｄｉｇｉｔａｌ　Ｍｏｉｒｅ　Ｆｒｉｎｇｅ　ＳｈｉｆｔｉｎｇＦｒａｃｔｉｏｎａｌＦｒｉｎｇｅＯｒｄｅｒＤｅｔｅｒｍｉｎａｔｉｏｎｕｓｉｎｇＤｉｇｉｔａｌＭｏｉｒｅＦ...     

Hole-drilling method using grating rosette and Moiré interferometry

The hole-drilling method is one of the most wellknown methods for measuring residual stresses. To identify unknown plane stresses in a specimen, a circular hole is first drilled in the infinite plate under plane stress, then the strains resulting from the hole drilling is measured. The strains may be acquired from interpreting the Moire signature around the hole. In crossed grating Moire interferometry, the horizontal and vertical displacement fields （u and v） can be obtained to determinate two strain fields and one shearing strain field. In this paper, by means of Moire interferometry and three directions grating （grating rosette） developed by the authors, three displacement fields （u, v and s） are obtained to acquire three strain fields. As a practical application, the hole-drilling method is adopted to measure the relief strains for aluminum and fiber reinforced composite. It is a step by step method; in each step a single laminate or equivalent depth is drilled to find some relationships between the drilling depth and the residual strains relieved in the fiber reinforced composite materials.     

Comparative study of in-cylinder tumble flows in an internal combustion engine using different piston shapes—an insight using particle image velocimetry

This paper deals with the experimental investigations of the in-cylinder tumble flows in a single-cylinder engine with five different piston crown shapes at an engine speed of 1,000 rev/min., during suction and compression strokes under motoring conditions using particle image velocimetry. Two-dimensional in-cylinder tumble flow measurements and analysis are carried out in combustion space on a vertical plane passing through cylinder axis. Ensemble average velocity vectors are used to analyze the tumble flow structure. Tumble ratio and average turbulent kinetic energy are evaluated and used to characterize the tumble flows. From results, it is found that at end of compression, pentroof-offset-bowl piston shows about 41 and 103% improvement in tumble ratio and average turbulent kinetic energy respectively, compared to that of flat piston. The present study will be useful in understanding effect of piston crown shapes on nature of the in-cylinder fluid tumble flows under real engine conditions.     

Wrinkling in graded core/shell systems using symplectic formulation

Wrinkles in flat graded elastic layers have been recently described as a timevarying Hamiltonian system by the energy method. Cylindrical core/shell structures can also undergo surface instabilities under the external pressure. In this study, we show that by treating the radial direction as a pseudo-time variable, the graded core/shell system with radially decaying elastic properties can also be described within the symplectic framework. In combination with the shell buckling equation, the prese...     

Investigation of liquid–liquid drop coalescence using tomographic PIV

High-speed tomographic PIV was used to investigate the coalescence of drops placed on a liquid/liquid interface; the coalescence of a single drop and of a drop in the presence of an adjacent drop (side-by-side drops) was investigated. The viscosity ratio between the drop and surrounding fluids was 0.14, the Ohnesorge number (Oh = μ_d/(ρ_dσD)^1/2) was 0.011, and Bond numbers (Bo = (ρ _d  − ρ _s )gD ²/σ) were 3.1–7.5. Evolving volumetric velocity fields of the full coalescence process allowed for quantification of the velocity scales occurring over different time scales. For both single and side-by-side drops, the coalescence initiates with an off-axis film rupture and film retraction speeds an order of magnitude larger than the collapse speed of the drop fluid. This is followed by the formation and propagation of an outward surface wave along the coalescing interface with wavelength of approximately 2D. For side-by-side drops, the collapse of the first drop is asymmetric due to the presence of the second drop and associated interface deformation. Overall, tomographic PIV provides insight into the flow physics and inherent three-dimensionalities in the coalescence process that would not be achievable with flow visualization or planar PIV only.     

Modeling of oil–water flow using energy minimization concept

A new model is developed to predict flow behaviors including flow pattern, pressure gradient and holdup for oil–water flow in horizontal and slightly inclined pipes. The model is based on the universal principle that a system stabilizes to its minimum total energy. The structural configurations observed in two-phase flow systems can be interpreted in terms of total energy minimization. Performance of the developed model is tested against several experimental data, and comparisons with existing models are presented. It is evident from the results and comparisons that the model estimates the pressure gradient and flow pattern very well. The model provides extensive information about oil–water flow characteristics.     

Vertical vibration control using nonlinear energy sink with inertial amplifier

To reduce additional mass, this work proposes a nonlinear energy sink(NES)with an inertial amplifier(NES-IA) to control the vertical vibration of the objects under harmonic and shock excitations. Moreover, this paper constructs pure nonlinear stiffness without neglecting the gravity effect of the oscillator. Both analytical and numerical methods are used to evaluate the performance of the NES-IA. The research findings indicate that even if the actual mass is 1% of the main oscillator, the NES-IA...     

Determining the Mie potential parameters using Poisson’s ratio

The known value of Poisson’s ratio specifying the relation between the strains along the principal directions in the case of uniaxial strain is used to propose an approach to derive an equation relating this ratio to the exponents of the Mie pair potential. An example of determining one of these exponents is discussed when the other exponent is given.     

A multilevel treatment of moiré fringe data using finite elements

Finite element procedures are applied to the modeling, analysis and visualization of experimental moiré data. Smoothing elements are introduced and evaluated with respect to data sparseness and error. A one-dimensional smoothing element is uniquely coupled with the method of principal curves to extract moiré fringe centers. A two-dimensional smoothing element is then used to produce a full-field representation given the fringe locations. The moiré technique is applied to the four-point bend experiment, and the surface-modeling technique is used to obtain displacement and gradient (strain) information.     

Block-structured compressible Navier–Stokes solution using the OPS high-level abstraction

abstract

In this paper, we report the development and validation of a compressible solver with shock capturing, using a domain-specific high-level abstraction framework, OPS, that is being developed at the University of Oxford. OPS uses an active library approach for block-structured meshes, capable of generating codes for a variety of parallel implementations with different parallelisation strategies. Performance results on various architectures are reported for the 1D Shu–Osher test case.     

Visualization of convection patterns near an evaporating meniscus using μPIV

Experimental visualizations of the three dimensional (3D) convection patterns generated near an evaporating meniscus in horizontally oriented capillary tubes are presented. These patterns are caused due to the differential evaporation along the meniscus. In this study, transparent capillary tubes with refractive index close to that of the evaporating liquid were used to minimize refraction effects and obtain velocity vectors near the walls. Polystyrene fluorescent particles of 0.5 μm diameter suspended in methanol were used to make the measurements in tubes of 75, 200 and 400 μm diameter. For the 75 μm tube, gravity was observed to have no effect on the flow patterns and an axisymmetric counter-rotating vortex pair was present along the horizontal and vertical center planes, suggesting the presence of a toroidal vortex near the meniscus. With an increase in tube size, buoyancy effects became apparent as the axisymmetric pattern broke down. A counter-rotating symmetric vortex pair was observed in the horizontal center plane, whereas in the vertical center plane, a single vortex dominated the flow and pushed the secondary vortex to a corner. Particle streak and μPIV images were obtained in multiple horizontal planes and a vertical center plane to understand this 3D flow behavior.     

Three-dimensional crack tip deformation measurement using combined Moiré-Sagnac interferometry

A combined Moiré-Sagnac interferometry method is developed for in-plane (u andv) and out-of-plane (w) surface deformation measurement. The combined optical setup is used to measure three-dimensional crack tip deformations of AI 2024-0 and AI 2024-T4 specimens at room temperature and an inconel 909 specimen at 570°C. Measured displacements near the crack tip region of the AI 2024-T4 specimen are used as input nodal displacements to determine stress intensity factors based on two-dimensional and three-dimensional Jacobian derivative method. The values compare favorably with theoretical calculations. The extent of the three-dimensional crack tip deformation zone is also discussed.     

Cellular-level near-wall unsteadiness of high-hematocrit erythrocyte flow using confocal μPIV

In hemodynamics, the inherent intermittency of two-phase cellular-level flow has received little attention. Unsteadiness is reported and quantified for the first time in the literature using a combination of fluorescent dye labeling, time-resolved scanning confocal microscopy, and micro-particle image velocimetry (μPIV). The near-wall red blood cell (RBC) motion of physiologic high-hematocrit blood in a rectangular microchannel was investigated under pressure-driven flow. Intermittent flow was associated with (1) the stretching of RBCs as they passed through RBC clusters with twisting motions; (2) external flow through local obstacles; and (3) transitionary rouleaux formations. Velocity profiles are presented for these cases. Unsteady flow clustered in local regions. Extra-cellular fluid flow generated by individual RBCs was examined using submicron fluorescent microspheres. The capabilities of confocal μPIV post-processing were verified using synthetic raw PIV data for validation. Cellular interactions and oscillating velocity profiles are presented, and 3D data are made available for computational model validation.     

Viscoelastic flow past confined objects using a micro–macro approach

This paper is concerned with the numerical prediction of viscoelastic flow past a cylinder in a channel and a sphere in a cylinder using molecular-based models. The basis of the numerical method employed is a micro–macro model in which the polymer dynamics is described by the evolution of an ensemble of Brownian configuration fields. The spectral element method is used to discretize the equations in space. Comparisons are made between the macroscopic simulations based on the Oldroyd B constitutive model and microscopic simulations based on Hookean dumbbells, and excellent agreement is found. The micro–macro approach can be used to simulate models, such as the finitely extensible nonlinear elastic (FENE) dumbbell model, which do not possess a closed-form constitutive equation. Numerical simulations are performed for the FENE model. The influence of the model parameters on the flow is described and, in particular, the dependence of the drag as a function of the Weissenberg number.     

Shockwave�Cturbulent boundary layer interaction control using magnetically driven surface discharges

This study demonstrates the potential for shockwave?Cturbulent boundary layer interaction control in air using low current DC constricted surface discharges forced by moderate strength magnetic fields. An analytical model describing the physics of magnetic field forced discharge interaction with boundary layer flow is developed and compared to experiments. Experiments are conducted in a Mach 2.6 indraft air tunnel with discharge currents up to 300?mA and magnetic field strengths up to 5?Tesla. Separation- and non-separation-inducing shocks are generated with diamond-shaped shockwave generators located on the wall opposite to the surface electrodes, and flow properties are measured with schlieren imaging, static wall pressure probes and acetone flow visualization. The effect of plasma control on boundary layer separation depends on the direction of the Lorentz force (j × B). It is observed that by using a Lorentz force that pushes the discharge upstream, separation can be induced or further strengthened even with discharge currents as low as 30?mA in a 3-Tesla magnetic field. If shock-induced separation is present, it is observed that by using Lorentz force that pushes the discharge downstream, separation can be suppressed, but this required higher currents, greater than 80?mA. Acetone planar laser scattering is used to image the flow structure in the test section and the reduction in the size of recirculation bubble and its elimination are observed experimentally as a function of actuation current and magnetic field strength.     

Modeling finite thickness slab perforation using a coupled Eulerian–Lagrangian approach

A coupled Eulerian–Lagrangian (CEL) method can be used to model many types of dynamic events. Projectile penetration through solids is particularly well-suited to a CEL method. In this study the CEL method in the commercially-available code Abaqus was used to model a near rigid projectile perforating finite thickness concrete slabs. A near rigid projectile can be modeled as a Lagrangian material with distinct material interfaces, while the solid target can be modeled as an Eulerian material capable of large deformations. An improved concrete constitutive model is also described that was implemented into Abaqus as a user material model. A simplified stochastic model was also implemented to capture some of the heterogeneous nature of concrete. The CEL simulations are compared to experimental data to demonstrate the utility of this method for this type of perforation event.     

Modelling cell wall growth using a fibre-reinforced hyperelastic–viscoplastic constitutive law

A fibre-reinforced hyperelastic–viscoplastic model using a finite strain Finite Element (FE) analysis is presented to study the expansive growth of cell walls. Based on the connections between biological concepts and plasticity theory, e.g. wall-loosening and plastic yield, wall-stiffening and plastic hardening, the modelling of cell wall growth is established within a framework of anisotropic viscoplasticity aiming to represent the corresponding biology-controlled behaviour of a cell wall. In order to model in vivo growth, special attention is paid to the differences between a living cell and an isolated wall. The proposed hyperelastic–viscoplastic theory provides a unique framework to clarify the interplay between cellulose microfibrils and cell wall matrix and how this interplay regulates sustainable growth in a particular direction while maintaining the mechanical strength of the cell walls by new material deposition. Moreover, the effect of temperature is taken into account. A numerical scheme is suggested and FE case studies are presented and compared with experimental data.