A cost‐effective curvature calculation approach for interfacial flows on unstructured meshes

We present a simple and cost‐effective curvature calculation approach for simulations of interfacial flows on structured and unstructured grids. The interface is defined using volume fractions, and the interface curvature is obtained as a function of the gradients of volume fractions. The gradient computation is based on a recently proposed gradient recovery method that mimicks the least squares approach without the need to solve a system of equations and is quite easy to implement on arbitrary polygonal meshes. The resulting interface curvature is used in a continuum surface force formulation within the framework of a well‐balanced finite‐volume algorithm to simulate multiphase flows dominated by surface tension. We show that the proposed curvature calculation is at least as accurate as some of the existing approaches on unstructured meshes while being straightforward to implement on any mesh topology. Numerical investigations also show that spurious currents in stationary problems that are dependent on the curvature calculation methodology are also acceptably low using the proposed approach. Studies on capillary waves and rising bubbles in viscous flows lend credence to the ability of the proposed method as an inexpensive, robust, and reasonably accurate approach for curvature calculation and numerical simulation of multiphase flows.     

A mixed hp FEM for the approximation of fourth‐order singularly perturbed problems on smooth domains

We consider fourth‐order singularly perturbed problems posed on smooth domains and the approximation of their solution by a mixed Finite Element Method on the so‐called Spectral Boundary Layer Mesh. We show that the method converges uniformly, with respect to the singular perturbation parameter, at an exponential rate when the error is measured in the energy norm. Numerical examples illustrate our theoretical findings.     

The local ultraconvergence of high‐order finite element method for second‐order elliptic problems with constant coefficients over a rectangular partition

In this article, we will discuss the local ultraconvergence of high‐degree finite element method based on a rectangular partition for the second‐degree elliptic problem with constant coefficients in Ω ? ?² , u( y ) = 0 on ?Ω . Based on suitable regularity, ultraconvergence of the displacement of the extrapolated kth (k ≥ 3) degree finite element solution has been obtained by an extrapolation technique. Finally, numerical experiments are applied to demonstrate our theoretical findings.     

Evaluation of RBFs for volume data interpolation in CFD

A greedy method for choosing an optimum reduced set of control points is integrated with RBF interpolation and evaluated for the purpose of interpolating large‐volume data sets in CFD. Given a function defined at a set of points, the greedy method selects a small subset of these points that is sufficient to keep the interpolation error at all the remaining points below a chosen bound. This is equivalent to a type of data compression and would have useful storage, post‐processing, and computational applications in CFD. To test the method in terms of both the point selection scheme and the suitability of reduced control point volume interpolation, a trial application of the interpolation to velocity fields in CFD volume meshes is considered. To optimise the point selection process, and attempt to be able to capture multiple length scales, a variable support radius formulation has also been included. Structured and unstructured mesh cases are considered for aerofoils, a wing case and a wing‐body case. For smooth volume functions, the method is shown to work well, producing accurate velocity interpolations using a very small number of the cells in the mesh. For general complex fields including large gradients, the method is still shown to be effective, although large gradients require more interpolation points to be used.Copyright © 2013 John Wiley & Sons, Ltd.     

Exact pressure integrations on submerged bodies in waves using a quadtree adaptive mesh algorithm

The development of an adaptive free surface, mesh cutting, methodology, in order to analytically integrate pressures on varying wet parts of partially submerged surfaces in the presence of waves, is presented. Given a function of free‐surface elevation, the algorithm checks for the intersection of the body with the free surface and, based on user‐defined parameters, modifies the initial mesh, by subdividing the elements where necessary and eliminating others, via a quadtree approach. Redundant sub‐divisions, generated in the quad‐division process, are partially eliminated, but the quadrilateral nature of the elements is always kept. The free‐surface function must be single‐valued and its definition domain simply connected. Hydrostatic and Froude–Krylov forces are computed exactly on each panel by means of analytical formulations, which are derived and presented, based on the theory of linear gravity waves and from applying Green's theorem. Copyright © 2014 John Wiley & Sons, Ltd.     

Direct and adjoint global stability analysis of turbulent transonic flows over a NACA0012 profile

In this work, various turbulent solutions of the two‐dimensional (2D) and three‐dimensional compressible Reynolds averaged Navier–Stokes equations are analyzed using global stability theory. This analysis is motivated by the onset of flow unsteadiness (Hopf bifurcation) for transonic buffet conditions where moderately high Reynolds numbers and compressible effects must be considered. The buffet phenomenon involves a complex interaction between the separated flow and a shock wave. The efficient numerical methodology presented in this paper predicts the critical parameters, namely, the angle of attack and Mach and Reynolds numbers beyond which the onset of flow unsteadiness appears. The geometry, a NACA0012 profile, and flow parameters selected reproduce situations of practical interest for aeronautical applications. The numerical computation is performed in three steps. First, a steady baseflow solution is obtained; second, the Jacobian matrix for the RANS equations based on a finite volume discretization is computed; and finally, the generalized eigenvalue problem is derived when the baseflow is linearly perturbed. The methodology is validated predicting the 2D Hopf bifurcation for a circular cylinder under laminar flow condition. This benchmark shows good agreement with the previous published computations and experimental data. In the transonic buffet case, the baseflow is computed using the Spalart–Allmaras turbulence model and represents a mean flow where the high frequency content and length scales of the order of the shear‐layer thickness have been averaged. The lower frequency content is assumed to be decoupled from the high frequencies, thus allowing a stability analysis to be performed on the low frequency range. In addition, results of the corresponding adjoint problem and the sensitivity map are provided for the first time for the buffet problem. Finally, an extruded three‐dimensional geometry of the NACA0012 airfoil, where all velocity components are considered, was also analyzed as a Triglobal stability case, and the outcoming results were compared to the previous 2D limited model, confirming that the buffet onset is well detected. Copyright © 2014 John Wiley & Sons, Ltd.     

A high‐order element based adaptive mesh refinement strategy for three‐dimensional unstructured grid

Adaptive mesh refinement (AMR) shows attractive properties in automatically refining the flow region of interest, and with AMR, better prediction can be obtained with much less labor work and cost compared to manually remeshing or the global mesh refinement. Cartesian AMR is well established; however, AMR on hybrid unstructured mesh, which is heavily used in the high‐Reynolds number flow simulation, is less matured and existing methods may result in degraded mesh quality, which mostly happens in the boundary layer or near the sharp geometric features. User intervention or additional constraints, such as freezing all boundary layer elements or refining the whole boundary layer, are required to assist the refinement process. In this work, a novel AMR strategy is developed to handle existing difficulties. In the new method, high‐order unstructured elements are first generated based on the baseline mesh; then the refinement is conducted in the parametric space; at last, the mesh suitable for the solver is output. Generating refined elements in the parametric space with high‐order elements is the key of this method and this helps to guarantee both the accuracy and robustness. With the current method, 3‐dimensional hybrid unstructured mesh of huge size and complex geometry can be automatically refined, without user intervention nor additional constraints. With test cases including the 2‐dimensional airfoil and 3‐dimensional full aircraft, the current AMR method proves to be accurate, simple, and robust.     

Synthesis of Stimuli-responsive Nanoparticles by Solution Polymerization

In this study, stimuli-responsive nanoparticles were prepared by solution polymerization. Two synthesis routes are proposed to synthesize the particles, the monomer route and the polymer/monomer route. For the monomer route, pH and thermal sensitive nanoparticles were synthesized from acrylic acid and N-isopropylacrylamide. For the polymer/monomer route, the pH sensitive nanoparticles were synthesized from chitosan and acrylic acid. The effect of reaction time, initiator concentration and agitation rate on the particle size and the size distribution were investigated. The stimuli-responsive nanoparticles could be directly blended with other polymers to prepare stimuli-responsive functional membranes.     

Search for a Steric Penultimate Effect: Copolymerization of 2,3,4-Trimethyl-3-pentyl Methacrylate with Styrene

The hindered monomer, 2,3,4-trimethyl-3-pentyl methacrylate (I), was synthesized for penultimate effect studies. Since it readily homopoiymerized (km¹¹¹≠ 0) and readily copolymerized with styrene, copolymerizations of I with styrene were carried out at 60°C in benzene with AIBN as initiator. The conversion to copolymer and the copolymer composition were determined by using GLC techniques. Composition-conversion data was analyzed by performing a computerized nonlinear least-squares fitting to the integrated form of the penultimate model equation. The experimental design included the use of optimized M¹°/M²° ratios. The penultimate reactivity ratios calculated from these data were r¹′ = 0.23, r¹′= 0.59, r² = 0.59, r²′ = 1.34. Thus, when I is the penultimate unit, a terminal styryl radical prefers to add styrene, whereas when styrene is the penultimate unit, terminal styryl radicals prefer to add I. These results constitute the best evidence for a steric penultimate effect yet available in the literature from composition-conversion studies. However, the case is not yet proved. Further studies to strengthen this conclusion are proposed.     

The Copolymerization of Alkenes with Maleic Anhydride

The investigations presented deal with the experimental results of the copolymerization of maleic anhydride (MAn) with alkenes. The course of the reaction is explained by the overall rate of the copolymerization (v _Br), which correlates with the solution viscosity of the copolymer, and the dependence of the v _Br maximum on the mole ratio of the monomers at constant total monomer concentration. The use of solvents with increasing donor power leads to increased complexing of the free MAn molecules and of the MAn radical chain ends. The results demonstrate that, for low 1-alkenes, the addition of the MAn chain radical is the rate-determining step of the copolymerization. As the substituents on the olefinic double bond become larger or the double bond shifts to the 1,2-position, the addition of MAn to the hydrocarbon radical becomes more and more the rate-determining step. On the other hand, an increase of the CT complexation of the MAn polymer radical by use of donor solvents decreases the alkene addition rate.