Effect of imperfections on the yielding of two-dimensional foams

The influence of each of the six different types of morphological imperfection - waviness, non-uniform cell wall thickness, cell-size variations, fractured cell walls, cell-wall misalignments, and missing cells - on the yielding of 2D cellular solids has been studied systematically for biaxial loading. Emphasis is placed on quantifying the knock-down effect of these defects on the hydrostatic yield strength and upon understanding the associated deformation mechanisms. The simulations in the present study indicate that the high hydrostatic strength, characteristic of ideal honeycombs, is reduced to a level comparable with the deviatoric strength by several types of defect. The common source of this large knock-down is a switch in deformation mode from cell wall stretching to cell wall bending under hydrostatic loading. Fractured cell edges produce the largest knock-down effect on the yield strength of 2D foams, followed in order by missing cells, wavy cell edges, cell edge misalignments, Γ Voronoi cells, δ Voronoi cells, and non-uniform wall thickness. A simple elliptical yield function with two adjustable material parameters successfully fits the numerically predicted yield surfaces for the imperfect 2D foams, and shows potential as a phenomenological constitutive law to guide the design of structural components made from metallic foams.

Computer generated hologram with geometric occlusion using GPU-accelerated depth buffer rasterization for three-dimensional display

We present a method of rapidly producing computer-generated holograms that exhibit geometric occlusion in the reconstructed image. Conceptually, a bundle of rays is shot from every hologram sample into the object volume.We use z buffering to find the nearest intersecting object point for every ray and add its complex field contribution to the corresponding hologram sample. Each hologram sample belongs to an independent operation, allowing us to exploit the parallel computing capability of modern programmable graphics processing units (GPUs). Unlike algorithms that use points or planar segments as the basis for constructing the hologram, our algorithm's complexity is dependent on fixed system parameters, such as the number of ray-casting operations, and can therefore handle complicated models more efficiently. The finite number of hologram pixels is, in effect, a windowing function, and from analyzing the Wigner distribution function of windowed free-space transfer function we find an upper limit on the cone angle of the ray bundle. Experimentally, we found that an angular sampling distance of 0:01' for a 2:66' cone angle produces acceptable reconstruction quality. © 2009 Optical Society of America.

The direct field boundary impedance of two-dimensional periodic structures with application to high frequency vibration prediction.

Large sections of many types of engineering construction can be considered to constitute a two-dimensional periodic structure, with examples ranging from an orthogonally stiffened shell to a honeycomb sandwich panel. In this paper, a method is presented for computing the boundary (or edge) impedance of a semi-infinite two-dimensional periodic structure, a quantity which is referred to as the direct field boundary impedance matrix. This terminology arises from the fact that none of the waves generated at the boundary (the direct field) are reflected back to the boundary in a semi-infinite system. The direct field impedance matrix can be used to calculate elastic wave transmission coefficients, and also to calculate the coupling loss factors (CLFs), which are required by the statistical energy analysis (SEA) approach to predicting high frequency vibration levels in built-up systems. The calculation of the relevant CLFs enables a two-dimensional periodic region of a structure to be modeled very efficiently as a single subsystem within SEA, and also within related methods, such as a recently developed hybrid approach, which couples the finite element method with SEA. The analysis is illustrated by various numerical examples involving stiffened plate structures.

Simulation of three-dimensional viscous flow in turbomachinery geometries using a solution-adaptive unstructured mesh methodology

This paper presents a numerical method for the simulation of flow in turbomachinery blade rows using a solution-adaptive mesh methodology. The fully three-dimensional, compressible, Reynolds-averaged Navier-Stokes equations with k-ε turbulence modeling (and low Reynolds number damping terms) are solved on an unstructured mesh formed from tetrahedral finite volumes. At stages in the solution, mesh refinement is carried out based on flagging cell faces with either a fractional variation of a chosen variable (like Mach number) greater than a given threshold or with a mean value of the chosen variable within a given range. Several solutions are presented, including that for the highly three-dimensional flow associated with the corner stall and secondary flow in a transonic compressor cascade, to demonstrate the potential of the new method.

Finite volume distance field solution applied to medial axis transform

Accurate and efficient computation of the nearest wall distance d (or level set) is important for many areas of computational science/engineering. Differential equation-based distance/ level set algorithms, such as the hyperbolic-natured Eikonal equation, have demonstrated valuable computational efficiency. Here, in the context, as an 'auxiliary' equation to the main flow equations, the Eikonal equation is solved efficiently with two different finite volume approaches (the cell vertex and cell-centered). Application of the distance solution is studied for various geometries. Moreover, a procedure using the differential field to obtain the medial axis transform (MAT) for different geometries is presented. The latter provides a skeleton representation of geometric models that has many useful analysis properties. As an alternative approach to the pure geometric methods (e.g. the Voronoi approach), the current d-MAT procedure bypasses many difficulties that are usually encountered by pure geometric methods, especially in three dimensional space. It is also shown that the d-MAT approach provides the potential to sculpt/control the MAT form for specialized solution purposes. Copyright © 2010 by the American Institute of Aeronautics and Astronautics, Inc.

Three-dimensional effects on sliding and waving wings

Like large insects, micro air vehicles operate at low Reynolds numbers O(1; 000 - 10; 000) in a regime characterized by separated flow and strong vortices. The leading-edge vortex has been identified as a significant source of high lift on insect wings, but the conditions required for the formation of a stably attached leading-edge vortex are not yet known. The waving wing is designed to model the translational phase of an insect wing stroke by preserving the unsteady starting and stopping motion as well as three-dimensionality in both wing geometry (via a finite-span wing) and kinematics (via wing rotation). The current study examines the effect of the spanwise velocity gradient on the development of the leading-edge vortex along the wing as well as the effects of increasing threedimensionalityby decreasing wing aspect ratio from four to two. Dye flow visualization and particle image velocimetry reveal that the leading-edge vortices that form on a sliding or waving wing have a very high aspect ratio. The structure of the flow is largely two-dimensional on both sliding and waving wings and there is minimal interaction between the leading-edge vortices and the tip vortex. Significant spanwise flow was observed on the waving wing but not on the sliding wing. Despite the increased three-dimensionality on the aspect ratio 2 waving wing, there is no evidence of an attached leading-edge vortex and the structure of the flow is very similar to that on the higher-aspect-ratio wing and sliding wing. © Copyright 2010.

Compressive failure of finite size unidirectional composite laminates with a region of fibre waviness

The compressive behaviour of finite unidirectional composites with a region of misaligned reinforcement is investigated via finite element analyses. Models with and without fibre bending stiffness are compared, confirming that compressive strength is accurately predicted without modelling fibre bending stiffness for real composite components which typically have waviness defects of several millimetres wavelength. Various defect parameters are investigated. Results confirm the well-known sensitivity of compressive strength to misalignment angle, and also show that compressive strength falls rapidly with the proportion of laminate width covered by the wavy region. A simple empirical equation is proposed to model the effect of a single patch of waviness in finite specimens. Other parameters such as length and position of the wavy region are found to have a smaller effect on compressive strength. The modelling approach is finally adapted to model distributed waviness and thus determine the compressive strength of composites with realistic waviness defects. © 2011 Elsevier Ltd. All rights reserved.

Polarisation diversity of conformal arrays based on geometric algebra via convex optimisation

A simple and general design procedure is presented for the polarisation diversity of arbitrary conformal arrays; this procedure is based on the mathematical framework of geometric algebra and can be solved optimally using convex optimisation. Aside from being simpler and more direct than other derivations in the literature, this derivation is also entirely general in that it expresses the transformations in terms of rotors in geometric algebra which can easily be formulated for any arbitrary conformal array geometry. Convex optimisation has a number of advantages; solvers are widespread and freely available, the process generally requires a small number of iterations and a wide variety of constraints can be readily incorporated. The study outlines a two-step approach for addressing polarisation diversity in arbitrary conformal arrays: first, the authors obtain the array polarisation patterns using geometric algebra and secondly use a convex optimisation approach to find the optimal weights for the polarisation diversity problem. The versatility of this approach is illustrated via simulations of a 7×10 cylindrical conformal array. © 2012 The Institution of Engineering and Technology.

Physical consequences of a nonparametric uncertainty model in structural dynamics

One of the main claims of the nonparametric model of random uncertainty introduced by Soize (2000) [3] is its ability to account for model uncertainty. The present paper investigates this claim by examining the statistics of natural frequencies, total energy and underlying dispersion equation yielded by the nonparametric approach for two simple systems: a thin plate in bending and a one-dimensional finite periodic massspring chain. Results for the plate show that the average modal density and the underlying dispersion equation of the structure are gradually and systematically altered with increasing uncertainty. The findings for the massspring chain corroborate the findings for the plate and show that the remote coupling of nonadjacent degrees of freedom induced by the approach suppresses the phenomenon of mode localization. This remote coupling also leads to an instantaneous response of all points in the chain when one mass is excited. In the light of these results, it is argued that the nonparametric approach can deal with a certain type of model uncertainty, in this case the presence of unknown terms of higher or lower order in the governing differential equation, but that certain expectations about the system such as the average modal density may conflict with these results. © 2012 Elsevier Ltd.

Direction and polarisation estimation using polarised cylindrical conformal arrays

In this study, the authors describe two-dimensional direction finding and signal polarisation estimation from a cylindrical conformal array consisting of directional and polarised antenna elements. Firstly, a simple and general transformation procedure, based on the mathematical framework of geometric algebra, is presented for arbitrary conformal arrays with polarised and directional antennas. Subsequently, the authors utilise the symmetry of cylindrical arrays to estimate signal parameters via rotational invariance techniques. The authors show how to iteratively estimate the azimuth and elevation angles of the incident signal, as well as its polarisation. To illustrate the versatility of this method, the results of simulations on a 3×4 cylindrical conformal array are shown and discussed. © 2012 The Institution of Engineering and Technology.

Fluid-structure interaction of three-dimensional magnetic artificial cilia

A numerical model is developed to analyse the interaction of artificial cilia with the surrounding fluid in a three-dimensional setting in the limit of vanishing fluid inertia forces. The cilia are modelled using finite shell elements and the fluid is modelled using a boundary element approach. The coupling between both models is performed by imposing no-slip boundary conditions on the surface of the cilia. The performance of the model is verified using various reference problems available in the literature. The model is used to simulate the fluid flow due to magnetically actuated artificial cilia. The results show that narrow and closely spaced cilia create the largest flow, that metachronal waves along the width of the cilia create a significant flow in the direction of the cilia width and that the recovery stroke in the case of the out-of-plane actuation of the cilia strongly depends on the cilia width. © 2012 Cambridge University Press.

3D modeling of high-T c superconductors by finite element software

A three-dimensional (3D) numerical model is proposed to solve the electromagnetic problems involving transport current and background field of a high-T c superconducting (HTS) system. The model is characterized by the E-J power law and H-formulation, and is successfully implemented using finite element software. We first discuss the model in detail, including the mesh methods, boundary conditions and computing time. To validate the 3D model, we calculate the ac loss and trapped field solution for a bulk material and compare the results with the previously verified 2D solutions and an analytical solution. We then apply our model to test some typical problems such as superconducting bulk array and twisted conductors, which cannot be tackled by the 2D models. The new 3D model could be a powerful tool for researchers and engineers to investigate problems with a greater level of complicity.