Efficient discretization of pressure‐correction equations on non‐orthogonal grids

An alternative discretization of pressure‐correction equations within pressure‐correction schemes for the solution of the incompressible Navier–Stokes equations is introduced, which improves the convergence and robustness properties of such schemes for non‐orthogonal grids. As against standard approaches, where the non‐orthogonal terms usually are just neglected, the approach allows for a simplification of the pressure‐correction equation to correspond to 5‐point or 7‐point computational molecules in two or three dimensions, respectively, but still incorporates the effects of non‐orthogonality. As a result a wide range (including rather high values) of underrelaxation factors can be used, resulting in an increased overall performance of the underlying pressure‐correction schemes. Within this context, a second issue of the paper is the investigation of the accuracy to which the pressure‐correction equation should be solved in each pressure‐correction iteration. The scheme is investigated for standard test cases and, in order to show its applicability to practical flow problems, for a more complex configuration of a micro heat exchanger. Copyright © 2003 John Wiley & Sons, Ltd.     

Axisymmetric frictionless contact of functionally graded materials

The main interest of this study is a new method to solve the axisymmetric frictionless contact problem of functionally graded materials (FGMs). Based on the fact that an arbitrary curve can be approached by a series of continuous but piecewise linear curves, the FGM is divided into a series of sub-layers with shear modulus varying linearly in each sub-layer and continuous at the sub-interfaces. With this model, the axisymmetric frictionless contact problem of a functionally graded coated half-space is investigated. By using the transfer matrix method and Hankel integral transform technique, the problem is reduced to a Cauchy singular integral equation. The contact pressure, contact region and indentation are calculated for various indenters by solving the equations numerically. An erratum to this article can be found at     

Numerical method and models for anti-plane strain of a system with a thin elastic wedge

The paper presents a general method to find asymptotics for a (multi-)wedge system containing a thin wedge. It employs separation of the symmetric and anti-symmetric parts of the boundary displacements and tractions of the wedge. The method is applicable when the angle of the thin wedge turns to zero. A physical interpretation of the derived equations is obtained by using power expansions of non-polynomial functions, which appear after the Mellin transform. We establish that the first term in the expansion of the symmetric part corresponds to shear, while the first term of the anti-symmetric part describes deflection of the wedge axis. Numerical experiments, performed by using a code developed on the basis of the theory, show that using only the first terms of the expansions insignificantly influence accuracy: the approximate results coincide with the exact values of roots to the third significant digit even for the wedge angle of 30°.     

A new VOF‐based numerical scheme for the simulation of fluid flow with free surface. Part II: application to the cavity filling and sloshing problems

Finite element analysis of fluid flow with moving free surface has been performed in 2‐D and 3‐D. The new VOF‐based numerical algorithm that has been proposed by the present authors (Int. J. Numer. Meth. Fluids, submitted) was applied to several 2‐D and 3‐D free surface flow problems. The proposed free surface tracking scheme is based on two numerical tools; the orientation vector to represent the free surface orientation in each cell and the baby‐cell to determine the fluid volume flux at each cell boundary. The proposed numerical algorithm has been applied to 2‐D and 3‐D cavity filling and sloshing problems in order to demonstrate the versatility and effectiveness of the scheme. The proposed numerical algorithm resolved successfully the free surfaces interacting with each other. The simulated results demonstrated applicability of the proposed numerical algorithm to the practical problems of large free surface motion. It has been also demonstrated that the proposed free surface tracking scheme can be easily implemented in any irregular non‐uniform grid systems and can be extended to 3‐D free surface flow problems without additional efforts. Copyright © 2003 John Wiley & Sons, Ltd.     

Structure–electrical property study of anisotropic polyaniline films

The relationships of the structure and electrical properties of anisotropic HCl‐doped polyaniline (PANI) films cast from N,N′‐dimethylpropylene urea (DMPU) solutions and stretched to different draw ratios were studied. The anisotropic structure of the stretched PANI films was examined by X‐ray diffraction, near‐infrared wave‐guide coupling, and polarized infrared measurements. The PANI emeraldine base (EB) films cast from DMPU solutions had a single‐phase noncrystalline structure, and stretching of the films did not cause crystallization to occur. The transition moment angles of two weakly absorbing infrared bands were determined, and the Hermans' orientation functions for the PANI EB films were calculated. The PANI films were then doped with HCl, and the electrical properties were determined by impedance spectroscopy. With a specially designed test fixture, the in‐plane and through‐plane impedance was obtained. The conductivity along the stretch direction increased with orientation. The in‐plane conductivity was significantly higher than the through‐plane conductivity. © 2003 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 41: 823–841, 2003     

Structure induced by irradiation of femtosecond laser pulse in dyed polymeric materials

We investigated the structures induced by an irradiation of a near‐infrared (NIR) femtosecond laser pulse in dye‐doped polymeric materials {poly(methyl methacrylate) (PMMA), thermoplastic epoxy resin (Epoxy), and a block copolymer of methyl methacrylate and ethyl acrylate‐butyl acrylate [p(MMA/EA‐BA) block copolymer]}. Dyes used were classified into two types—type 1 with absorption at 400 nm and type 2 with no absorption at 400 nm. The 400‐nm wavelength corresponds to the two‐photon absorption region by the irradiated NIR laser pulse at 800 nm. Type 1 dye‐doped PMMA and p(MMA/EA‐BA) block copolymer showed a peculiar dye additive effect for the structures induced by the line irradiation of a NIR femtosecond laser pulse. On the contrary, dye‐doped Epoxy did not exhibit a dye additive effect. The different results among PMMA, p(MMA/EA‐BA) block copolymer, and Epoxy matrix polymers are supposed to be related to the difference of electron‐acceptor properties. The mechanism of this type 1 dye‐additive‐effect phenomenon for PMMA and p(MMA/EA‐BA) block copolymer is discussed on the basis of two‐photon absorption of type 1 dye at 400 nm by the irradiation of a femtosecond laser pulse with 800 nm wavelength and the dissipation of the absorbed energy to the polymer matrix among various transition processes. Dyes with a low‐fluorescence quantum yield favored the formation of thicker grating structures. © 2002 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 40: 2800–2806, 2002     

Diffusion via splitting and remeshing via merging in vortex methods

The technique of splitting a fat vortex element (with a core width larger than some threshold) into some thin ones in order to fix the convergence problem of the core‐spreading vortex methods is convenient and efficient. In particular, it keeps the method purely Lagrangian. In the present investigation, the splitting process is further viewed as part of the physical diffusion process. A new splitting method in which several weaker child vortices surround a thinned but still strong parent vortex is proposed. It is found that because of the survival of the parent vortex, the error arising from the splitting events can be largely reduced. The computational amount on the other hand is kept reasonably large by merging similar and close‐by vortices. The merging scheme designed herein not only involves fewer restrictions but also allows merging vortices of opposite rotations through the viewpoint of remeshing. The validity and accuracy of these techniques, proposed particularly for simulations undergoing lots of splitting and merging events, are verified by successfully simulating the interactions between two Burgers vortices under an external straining field. Copyright © 2005 John Wiley & Sons, Ltd.     

On the use of the Galerkin method for 3D numerical modelling of the general circulation: the South Atlantic eperiment

A linear three‐dimensional hydrodynamical numerical model, with the application of the Galerkin Method for the vertical dependence, is here presented. The spherical coordinate system is used, in order to allow large‐scale simulations. The equations and mathematical development of the model are shown in detail, together with the boundary and initial conditions, and the sequence of equations' solution. The model is applied to the South Atlantic Ocean, for estimating typical seasonal circulations, and the results are summarized in maps of currents at surface and 1000 m depth, and in transport values of the Brazil Current between 30°S and 40°S. Copyright © 2002 John Wiley & Sons, Ltd.     

Least‐squares finite element method on adaptive grid for PDEs with shocks

We apply the least‐squares finite element method with adaptive grid to nonlinear time‐dependent PDEs with shocks. The least‐squares finite element method is also used in applying the deformation method to generate the adaptive moving grids. The effectiveness of this method is demonstrated by solving a Burgers' equation with shocks. Computational results on uniform grids and adaptive grids are compared for the purpose of evaluation. The results show that the adaptive grids can capture the shock more sharply with significantly less computational time. For moving shock, the adaptive grid moves correctly with the shock. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Synthesis and characterization of novel perfluorinated polyimides with 3D‐controlled structures for photonic applications

We have designed and synthesized novel perfluorinated polyimides with 3D controlled structure. At first, we successfully developed the new synthetic routes to diamines with pendant bulky perfluorinated aromatic units in a multi‐step synthetic procedure. Novel perfluorinated polyimides were prepared in a two‐step reaction in N‐methyl‐2‐pyrrolidinone (NMP) solution: The first step was for the synthesis of polyamic acids (PAAs) and the second reaction was for the imidization of PAA. The polymer yield was over 89% and the inherent viscosity of PAAs was in the range of 0.24–0.36 dL/g. The thin films were prepared by spin‐coating the PAA solution in NMP onto various substrates such as a Si wafer or a KBr pellet, dried at 80 °C and further cured at 230 °C. The resultant polyimides are thermally stable over 400 °C. The refractive index and birefringence of the resultant polyimides are 1.5858–1.6452 and 0.01–0.005, respectively. The refractive index of polyimide decreases with increasing the fluorine content. The copolymerization and the ether linkages into the backbone reduce the birefringence of polyimides. Surprisingly, the pendant ether linkage is not a crucial factor in reducing the polyimide birefringence. Their birefringence is comparatively very low, compared with that of previous polyimides. © 2005 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 44: 1326–1342, 2006