Energy stable and momentum conserving hybrid finite element method for the incompressible Navier-Stokes equations

A hybrid method for the incompressible Navier-Stokes equations is presented. The method inherits the attractive stabilizing mechanism of upwinded discontinuous Galerkin methods when momentum advection becomes significant, equal-order interpolations can be used for the velocity and pressure fields, and mass can be conserved locally. Using continuous Lagrange multiplier spaces to enforce flux continuity across cell facets, the number of global degrees of freedom is the same as for a continuous Galerkin method on the same mesh. Different from our earlier investigations on the approach for the Navier-Stokes equations, the pressure field in this work is discontinuous across cell boundaries. It is shown that this leads to very good local mass conservation and, for an appropriate choice of finite element spaces, momentum conservation. Also, a new form of the momentum transport terms for the method is constructed such that global energy stability is guaranteed, even in the absence of a pointwise solenoidal velocity field. Mass conservation, momentum conservation, and global energy stability are proved for the time-continuous case and for a fully discrete scheme. The presented analysis results are supported by a range of numerical simulations. © 2012 Society for Industrial and Applied Mathematics.

Thermal transport in high porosity cellular metal foams

The heat dissipation capability of highly porous cellular metal foams with open cells subject to forced air convection is studied using a combined experimental and analytical approach. The cellular morphologies of six FeCrAlY (an iron-based alloy) foams and six copper alloy foams with a range of pore sizes and porosities are quantified with the scanning electronic microscope and image analysis. Experimental measurements on pressure drop and heat transfer for copper foams are carried out. A numerical model for forced convection across open-celled metal foams is subsequently developed, and the predictions are compared with those measured. Reasonably good agreement with test data is obtained, given the complexity of the cellular foam morphology and the associated momentum/energy transport. The results show that cell size has a more significant effect on the overall heat transfer than porosity. An optimal porosity is obtained based on the balance between pressure drop and overall heat transfer, which decreases as the Reynolds number is increased.

ZnO film based surface acoustic wave micro-pump

In this study, a micro-pump unit based on surface acoustic wave (SAW) on piezoelectric ZnO film is designed and fabricated as a micro-fluidic device. It employs a mechanical wave, which is generated electrically using an aluminum interdigital transducer (IDT), and propagates on the surface of the ZnO film. The ZnO film was used in this study because it has a high electromechanical coefficient and an excellent bonding with various substrate materials, in particular silicon. The sputtering parameters for ZnO film deposition have been optimized, and the ZnO films with different thickness from 1 micron to 5.5 microns were prepared. The film properties have been characterized using different methods, such as scanning electron microscopy, X-ray diffraction and atomic force microscopy. Aluminum IDT with a finger width and spacing of 8 microns was patterned on the ZnO film using a lift-off process. The frequency generated was measured using a network analyzer, and it varies from 130 MHz to 180 MHz as a function of film thickness. A signal generator was used to generate the frequency with a power amplifier to amplify the signal, which was then applied to aluminum IDT to generate the surface acoustic wave. If a liquid droplet exists on the surface carrying the acoustic wave, the energy and the momentum of the SAW will be coupled into the fluid, causing the liquid to vibrate and move on film surface. The strength of this movement is determined by the applied voltage and frequency. The volume of the liquid drop loaded on the SAW device in this study is of several hundreds of nanoliters. The movement of the liquid inside the droplet and also on the ZnO film surface can be demonstrated. The performance of ZnO SAW device was characterized as a function of film thickness. © 2007 IOP Publishing Ltd.

Evolution of localized blobs of swirling or buoyant fluid with and without an ambient magnetic field.

We investigate the evolution of localized blobs of swirling or buoyant fluid in an infinite, inviscid, electrically conducting fluid. We consider the three cases of a strong imposed magnetic field, a weak imposed magnetic field, and no magnetic field. For a swirling blob in the absence of a magnetic field, we find, in line with others, that the blob bursts radially outward under the action of the centrifugal force, forming a thin annular vortex sheet. A simple model of this process predicts that the vortex sheet thins exponentially fast and that it moves radially outward with constant velocity. These predictions are verified by high-resolution numerical simulations. When an intense magnetic field is applied, this phenomenon is suppressed, with the energy and angular momentum of the blob now diffusing axially along the magnetic field lines, converting the blob into a columnar structure. For modest or weak magnetic fields, there are elements of both types of behavior, with the radial bursting dominating over axial diffusion for weak fields. However, even when the magnetic field is very weak, the flow structure is quite distinct to that of the nonmagnetic case. In particular, a small but finite magnetic field places a lower bound on the thickness of the annular vortex sheet and produces an annulus of counter-rotating fluid that surrounds the vortex core. The behavior of the buoyant blob is similar. In the absence of a magnetic field, it rapidly develops the mushroomlike shape of a thermal, with a thin vortex sheet at the top and sides of the mushroom. Again, a simple model of this process predicts that the vortex sheet at the top of the thermal thins exponentially fast and rises with constant velocity. These predictions are consistent with earlier numerical simulations. Curiously, however, it is shown that the net vertical momentum associated with the blob increases linearly in time, despite the fact that the vertical velocity at the front of the thermal is constant. As with the swirling blob, an imposed magnetic field inhibits the formation of a vortex sheet. A strong magnetic field completely suppresses the phenomenon, replacing it with an axial diffusion of momentum, while a weak magnetic field allows the sheet to form, but places a lower bound on its thickness. The magnetic field does not, however, change the net vertical momentum of the blob, which always increases linearly with time.