A spectral projection method for the analysis of autocorrelation functions and projection errors in discrete particle simulation

Discrete particle simulation is a well‐established tool for the simulation of particles and droplets suspended in turbulent flows of academic and industrial applications. The study of some properties such as the preferential concentration of inertial particles in regions of high shear and low vorticity requires the computation of autocorrelation functions. This can be a tedious task as the discrete point particles need to be projected in some manner to obtain the continuous autocorrelation functions. Projection of particle properties on to a computational grid, for instance, the grid of the carrier phase, is furthermore an issue when quantities such as particle concentrations are to be computed or source terms between the carrier phase and the particles are exchanged. The errors committed by commonly used projection methods are often unknown and are difficult to analyse. Grid and sampling size limit the possibilities in terms of precision per computational cost. Here, we present a spectral projection method that is not affected by sampling issues and addresses all of the above issues. The technique is only limited by computational resources and is easy to parallelize. The only visible drawback is the limitation to simple geometries and therefore limited to academic applications. The spectral projection method consists of a discrete Fourier‐transform of the particle locations. The Fourier‐transformed particle number density and momentum fields can then be used to compute the autocorrelation functions and the continuous physical space fields for the evaluation of the projection methods error. The number of Fourier components used to discretize the projector kernel can be chosen such that the corresponding characteristic length scale is as small as needed. This allows to study the phenomena of particle motion, for example, in a region of preferential concentration that may be smaller than the cell size of the carrier phase grid. The precision of the spectral projection method depends, therefore, only on the number of Fourier modes considered. Copyright © 2008 John Wiley & Sons, Ltd.     

A stabilized incremental projection scheme for the incompressible Navier–Stokes equations

It is well known that any spatial discretization of the saddle‐point Stokes problem should satisfy the Ladyzhenskaya–Brezzi–Babuska (LBB) stability condition in order to prevent the appearance of spurious pressure modes. Particularly, if an equal‐order approximation is applied, the Schur complement (or, as called some times, the Uzawa matrix) of the pressure system has a non‐trivial null space that gives rise to such modes. An idea in the past was that all the schemes that solve a Poisson equation for the pressure rather than the Uzawa pressure equation (splitting/projection methods) should overcome this difficulty; this idea was wrong. There is numerical evidence that at least the so‐called incremental projection scheme still suffers from spurious pressure oscillations if an equal‐order approximation is applied. The present paper tries to distinguish which projection requires LBB‐compliant approximation and which does not. Moreover, a stabilized version of the incremental projection scheme is derived. Proper bounds for the stabilization parameter are also given. The numerical results show that the stabilized scheme does indeed achieve second‐order accuracy and does not produce spurious (node to node) pressure oscillations. Copyright © 2001 John Wiley & Sons, Ltd.     

A projection scheme for incompressible multiphase flow using adaptive Eulerian grid

This paper presents a finite element method for incompressible multiphase flows with capillary interfaces based on a (formally) second‐order projection scheme. The discretization is on a fixed Eulerian grid. The fluid phases are identified and advected using a level set function. The grid is temporarily adapted around the interfaces in order to maintain optimal interpolations accounting for the pressure jump and the discontinuity of the normal velocity derivatives. The least‐squares method for computing the curvature is used, combined with piecewise linear approximation to the interface. The time integration is based on a formally second order splitting scheme. The convection substep is integrated over an Eulerian grid using an explicit scheme. The remaining generalized Stokes problem is solved by means of a formally second order pressure‐stabilized projection scheme. The pressure boundary condition on the free interface is imposed in a strong form (pointwise) at the pressure‐computation substep. This allows capturing significant pressure jumps across the interface without creating spurious instabilities. This method is simple and efficient, as demonstrated by the numerical experiments on a wide range of free‐surface problems. Copyright © 2004 John Wiley & Sons, Ltd.     

A class of fully third-order accurate projection methods for solving the incompressible Navier-Stokes equations

In this paper, a fully third-order accurate projection method for solving the incompressible Navier-Stokes equations is proposed. To construct the scheme, a continuous projection procedure is firstly presented. We then derive a sufficient condition for the continuous projection equations to be temporally third-order accurate approximations of the original Navier-Stokes equations by means of the local- truncation-error-analysis technique. The continuous projection equations are discretized temporally and spatially to third-order accuracy on the staggered grids, resulting in a fully third-order discrete projection scheme. The possibility to design higher-order projection methods is thus demonstrated in the present paper. A heuristic stability analysis is performed on this projection method showing the probability of its being stable. The stability of the present scheme is further verified through numerical tests. The third-order accuracy of the present projection method is validated by several numerical test cases.     

A class of fully third-order accurate projection methods for solving the incompressible Navier-Stokes equations

In this paper, a fully third-order accurate projection method for solving the incompressible Navier-Stokes equations is proposed. To construct the scheme, a continuous projection procedure is firstly presented. We then derive a sufficient condition for the continuous projection equations to be temporally third-order accurate approximations of the original Navier-Stokes equations by means of the local- truncation-error-analysis technique. The continuous projection equations are discretized temporally and spatially to third-order accuracy on the staggered grids, resulting in a fully third-order discrete projection scheme. The possibility to design higher-order projection methods is thus demonstrated in the present paper. A heuristic stability analysis is performed on this projection method showing the probability of its being stable. The stability of the present scheme is further verified through numerical tests. The third-order accuracy of the present projection method is validated by several numerical test cases. The project supported by the China NKBRSF (2001CB409604) The English text was polished by Yunming Chen     

A projection scheme for Navier-Stokes with variable viscosity and natural boundary condition

Combiningthe Navier-Stokes systems with Neumann (or natural) boundary condition to characterize a fluid flow is frequent. The popular projection (or pressure correction) methods inspired by Chorin and Temam are not well adapted to such boundary condition, which translate in loss of accuracy. If some alternative projection methods have been proposed to reduce the accuracy loss due to the Neumann condition in case of Newtonian fluids, little has been proposed for generalized Newtonian fluids. In this work, we propose two methods derived from the incremental pressure correction projection that can be used for fluids with inhomogeneous or variable viscosity with natural boundary condition. Both time and space accuracy of the methods will be illustrated using a manufactured solution.     

A Finite Element approach with patch projection for strain gradient plasticity formulations

Several strain gradient plasticity formulations have been suggested in the literature to account for inherent size effects on length scales of microns and submicrons. The necessity of strain gradient related terms render the simulation with strain gradient plasticity formulation computationally very expensive because quadratic shape functions or mixed approaches in displacements and strains are usually applied. Approaches using linear shape functions have also been suggested which are, however, limited to regular meshes with equidistanced Finite Element nodes. As a result the majority of the simulations in the literature deal with plane problems at small strains. For the solution of general three dimensional problems at large strains an approach has to be found which has to be computationally affordable and robust.     

A convenient look-up-table based method for the compensation of non-linear error in digital fringe projection

Although the structured light system that uses digital fringe projection has been widely implemented in three-dimensional surface profile measurement, the measurement system is susceptible to non-linear error. In this work, we propose a convenient look-up-table-based(LUT-based) method to compensate for the non-linear error in captured fringe patterns. Without extra calibration, this LUT-based method completely utilizes the captured fringe pattern by recording the full-field differences. Then, a phase compensation map is established to revise the measured phase. Experimental results demonstrate that this method works effectively.     

Remark on the pressure boundary condition for the projection method

Much has been said about the pressure boundary condition for the projection method, which is different from the actual boundary condition satisfied by the pressure in the Navier-Stokes equations. In this short note we present a different point of view which resolves the difficulty and we show how this point of view agrees with previous results.     

A novel tensile device for in situ scanning electron microscope mechanical testing

Experimental Techniques - To facilitate the study of deformation mechanisms and mechanical properties of bulk materials with feature size of centimeter level, a novel tensile device compatible with...     

X-ray and electron microscope examination of hardened steel in bending: impact, rising and constant load

Impact bending fracture of hardened steel is examined by X-ray measurements and electron microscopic fractographic data. They are compared to those for bending with rising and constant amplitude load.Obtained are residual stresses and diffraction lines the width of which are integrated with reference to the depth under the failure surface of the hardened steel. The characteristics of the plastic zone formed under the fracture surfaces are determined for the three different types of loading. In particular, the zones of initial failure and final rupture are analyzed by electron microscopy. The results are quantified by the texture of microstructure surface appearance with reference to local regions failured in a brittle and ductile fashion.     

A critical comparison of projection microradiography and the 2X method

Summary It is shown that Cosslett's criticism on the 2X method is based partly on a misconception of the resolving power of the photographic film and partly on a fundamentally wrong comparison between projection and 2X method. A generally valid comparison is given for the cases where the resolving power is not limited by Fresnel diffraction. Finally it is shown that Fresnel fringes are smaller in the 2X method than in a comparable projection system.The authors wish to express their gratitude to Dr. J.Bouman for the discussions which contributed especially to § 2.     

A new hybrid target concept for multi-keV X-ray sources

A novel concept for using hybrid targets to create multi-keV X-ray sources was tested on the GEKKO XII facility of the OSAKA University and on the OMEGA facility of the University of Rochester. The sources were made via laser irradiation of a titanium foil placed at the end of a plastic cylinder, filled with a very low-density (2 and 5 mg/cm³) silicon-dioxide aerogel that was designed to control the longitudinal expansion of the titanium plasma. Preliminary calculations were used to determine optimal conditions for the aerogel density, cylinder diameter and length that maximize multi-keV X-ray emission. The X-ray emission power was measured on OMEGA using absolutely calibrated broad-band, diode-based CEA diagnostics, in addition to high resolution crystal spectrometers. On GEKKO XII, the heat wave propagation velocity in the aerogel was also measured with an X-ray framing camera. The advantage of using the thermal wave generated in the aerogel to heat a solid material to increase the conversion efficiency has not been fully demonstrated in these experiments. However, it was shown that a 5 mg/cm³ aerogel placed in front of a titanium foil can improve the x-ray conversion efficiency with respect to the case of 2 mg/cm³ for some target diameter and length.     

A projection method for solving incompressible viscous flows on domains with moving boundaries

In this paper, a projection method is presented for solving the flow problems in domains with moving boundaries. In order to track the movement of the domain boundaries, arbitrary‐Lagrangian–Eulerian (ALE) co‐ordinates are used. The unsteady incompressible Navier–Stokes equations on the ALE co‐ordinates are solved by using a projection method developed in this paper. This projection method is based on the Bell's Godunov‐projection method. However, substantial changes are made so that this algorithm is capable of solving the ALE form of incompressible Navier–Stokes equations. Multi‐block structured grids are used to discretize the flow domains. The grid velocity is not explicitly computed; instead the volume change is used to account for the effect of grid movement. A new method is also proposed to compute the freestream capturing metrics so that the geometric conservation law (GCL) can be satisfied exactly in this algorithm. This projection method is also parallelized so that the state of the art high performance computers can be used to match the computation cost associated with the moving grid calculations. Several test cases are solved to verify the performance of this moving‐grid projection method. Copyright © 2004 John Wiley Sons, Ltd.     

Use of compression wave focusing for the production of a superjet

The principle of focusing compression waves by their reflection from nonflat surfaces is well-known. The invention in [1] postulates for this purpose the use in a pulsed hydrojet of a piston with a specially shaped surface. In this paper, the use of focusing is investigated numerically for the production of superjets.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 153–156, September–October, 1976.     

A finite element variational multiscale method for incompressible flows based on the construction of the projection basis functions

In this article, we present a finite element variational multiscale (VMS) method for incompressible flows based on the construction of projection basis functions and compare it with common VMS method, which is defined by a low‐order finite element space L_h on the same grid as X_h for the velocity deformation tensor and a stabilization parameter α. The best algorithmic feature of our method is to construct the projection basis functions at the element level with minimal additional cost to replace the global projection operator. Finally, we give some numerical simulations of the nonlinear flow problems to show good stability and accuracy properties of the method. Copyright © 2011 John Wiley & Sons, Ltd.     

A new microtensile tester for the study of MEMS materials with the aid of atomic force microscopy

An apparatus has been designed and implemented to measure the elastic tensile properties (Young's modulus and tensile strength) of surface micromachined polysilicon specimens. The tensile specimens are “dog-bone” shaped ending in a large “paddle” for convenient electrostatic or, in the improved apparatus, ultraviolet (UV) light curable adhesive gripping deposited with electrostatically controlled manipulation. The typical test section of the specimens is 400 μm long with 2 μm×50 μm cross section. The new device supports a nanomechanics method developed in our laboratory to acquire surface topologies of deforming specimens by means of Atomic Force Microscopy (AFM) to determine (fields of) strains via Digital Image Correlation (DIC). With this tool, high strength or non-linearly behaving materials can be tested under different environmental conditions by measuring the strains directly on the surface of the film with nanometer resolution.     

Subspace projection filters for imaging in random media

We consider the problem of selective imaging extended reflectors in cluttered media. We propose a random travel time model for simulating the array response matrix in clutter and we compare it with the full wave solution. Our simplified model captures very well the full wave random medium behavior as this is illustrated by our numerical results. The algorithm for selective array imaging uses coherent interferometry on a filtered version of the data. The filter, which is based on the singular value decomposition of the response matrix, enhances the signal reflected by the edges of the reflector. We illustrate the performance of the imaging algorithm with numerical simulations in the regime of ultrasonic non-destructive testing in concrete.