A Smooth Finite Element Method Based on Reproducing Kernel DMS-Splines

The element-based piecewise smooth functional approximation in the conventional finite element method (FEM) results in discontinuous first and higher order derivatives across element boundaries Despite the significant advantages of the FEM in modelling complicated geometries, a motivation in developing mesh-free methods has been the ease with which higher order globally smooth shape functions can be derived via the reproduction of polynomials There is thus a case for combining these advantages in a so-called hybrid scheme or a `smooth FEM' that, whilst retaining the popular mesh-based discretization, obtains shape functions with uniform C-p (p >= 1) continuity One such recent attempt, a NURBS based parametric bridging method (Shaw et al 2008b), uses polynomial reproducing, tensor-product non-uniform rational B-splines (NURBS) over a typical FE mesh and relies upon a (possibly piecewise) bijective geometric map between the physical domain and a rectangular (cuboidal) parametric domain The present work aims at a significant extension and improvement of this concept by replacing NURBS with DMS-splines (say, of degree n > 0) that are defined over triangles and provide Cn-1 continuity across the triangle edges This relieves the need for a geometric map that could precipitate ill-conditioning of the discretized equations Delaunay triangulation is used to discretize the physical domain and shape functions are constructed via the polynomial reproduction condition, which quite remarkably relieves the solution of its sensitive dependence on the selected knotsets Derivatives of shape functions are also constructed based on the principle of reproduction of derivatives of polynomials (Shaw and Roy 2008a) Within the present scheme, the triangles also serve as background integration cells in weak formulations thereby overcoming non-conformability issues Numerical examples involving the evaluation of derivatives of targeted functions up to the fourth order and applications of the method to a few boundary value problems of general interest in solid mechanics over (non-simply connected) bounded domains in 2D are presented towards the end of the paper

The Reproducing Kernel DMS-FEM: 3D Shape Functions and Applications to Linear Solid Mechanics

We propose a family of 3D versions of a smooth finite element method (Sunilkumar and Roy 2010), wherein the globally smooth shape functions are derivable through the condition of polynomial reproduction with the tetrahedral B-splines (DMS-splines) or tensor-product forms of triangular B-splines and ID NURBS bases acting as the kernel functions. While the domain decomposition is accomplished through tetrahedral or triangular prism elements, an additional requirement here is an appropriate generation of knotclouds around the element vertices or corners. The possibility of sensitive dependence of numerical solutions to the placements of knotclouds is largely arrested by enforcing the condition of polynomial reproduction whilst deriving the shape functions. Nevertheless, given the higher complexity in forming the knotclouds for tetrahedral elements especially when higher demand is placed on the order of continuity of the shape functions across inter-element boundaries, we presently emphasize an exploration of the triangular prism based formulation in the context of several benchmark problems of interest in linear solid mechanics. In the absence of a more rigorous study on the convergence analyses, the numerical exercise, reported herein, helps establish the method as one of remarkable accuracy and robust performance against numerical ill-conditioning (such as locking of different kinds) vis-a-vis the conventional FEM.

Smooth DMS-FEM: A new approach to solving nearly incompressible nonlinear elasto-static problems

The DMS-FEM, which enables functional approximations with C(1) or still higher inter-element continuity within an FEM-based meshing of the domain, has recently been proposed by Sunilkumar and Roy [39,40]. Through numerical explorations on linear elasto-static problems, the method was found to have conspicuously superior convergence characteristics as well as higher numerical stability against locking. These observations motivate the present study, which aims at extending and exploring the DMS-FEM to (geometrically) nonlinear elasto-static problems of interest in solid mechanics and assessing its numerical performance vis-a-vis the FEM. In particular, the DMS-FEM is shown to vastly outperform the FEM (presently implemented through the commercial software ANSYS (R)) as the former requires fewer linearization and load steps to achieve convergence. In addition, in the context of nearly incompressible nonlinear systems prone to volumetric locking and with no special numerical artefacts (e.g. stabilized or mixed weak forms) employed to arrest locking, the DMS-FEM is shown to approach the incompressibility limit much more closely and with significantly fewer iterations than the FEM. The numerical findings are suggestive of the important role that higher order (uniform) continuity of the approximated field variables play in overcoming volumetric locking and the great promise that the method holds for a range of other numerically ill-conditioned problems of interest in computational structural mechanics. (C) 2011 Elsevier Ltd. All rights reserved.

Wrinkled and slack membranes: nonlinear 3D elasticity solutions via smooth DMS-FEM and experiment

The smooth DMS-FEM, recently proposed by the authors, is extended and applied to the geometrically nonlinear and ill-posed problem of a deformed and wrinkled/slack membrane. A key feature of this work is that three-dimensional nonlinear elasticity equations corresponding to linear momentum balance, without any dimensional reduction and the associated approximations, directly serve as the membrane governing equations. Domain discretization is performed with triangular prism elements and the higher order (C1 or more) interelement continuity of the shape functions ensures that the errors arising from possible jumps in the first derivatives of the conventional C0 shape functions do not propagate because the ill-conditioned tangent stiffness matrices are iteratively inverted. The present scheme employs no regularization and exhibits little sensitivity to h-refinement. Although the numerically computed deformed membrane profiles do show some sensitivity to initial imperfections (nonplanarity) in the membrane profile needed to initiate transverse deformations, the overall patterns of the wrinkles and the deformed shapes appear to be less so. Finally, the deformed profiles, computed through the DMS FEM-based weak formulation, are compared with those obtained through an experiment on an ultrathin Kapton membrane, wherein wrinkles form because of the applied boundary displacement conditions. Comparisons with a reported experiment on a rectangular membrane are also provided. These exercises lend credence to the feasibility of the DMS FEM-based numerical route to computing post-wrinkled membrane shapes. Copyright (c) 2012 John Wiley & Sons, Ltd.

Molecular weight dependence of phase Behavior of PEO/P(EO-b-DMS) blends: Application of Sanchez-Lacombe lattice fluid theory

Molecular weight dependence of phase separation behavior of the Poly (ethylene oxide) (PEO)/Poly(ethylene oxide-block-dimethylsiloxane) (P(EO-b-DMS)) blends was investigated by both experimental and theoretical methods. The cloud point curves of PEO/P(EO-b-DMS) blends were obtained by turbidity method. Based on Sanchez-Lacombe lattice fluid theory (SLLFT), the adjustable parameter, epsilon*(12)/k (quantifying the interaction energy between different components), was evaluated by fitting the experimental data in phase diagrams. To calculate the spinodals, binodals, and the volume changes of mixing for these blends, three modified combining rules of the scaling parameters for the block copolymer were introduced.

Abnormal pressure dependence of the phase boundaries in TL/PEO/P(EO-b-DMS) ternary mixtures

The cloud-point temperatures (T-cl's) of both binary poly(ethylene oxide) (PEO)-poly(ethylene oxide-b-dimethylsiloxane) [P(EO-b-DMS)] and ternary[toluene/PEO/P(EO-b-DMS)] systems were determined by light scattering measurements at atmospheric pressure. The phase separation behavior upon cooling in the ternary system has been investigated at atmospheric pressure and under high pressure and compared to the phase behavior in the binary system. The phase transition temperatures have been obtained for all of the samples. As a result, the pressure induces compatibility in the binary mixtures, but for the ternary system, pressure not only can induce mixing but also can induce phase separation.

Thermodynamics of phase behavior in PEO/P(EO-b-DMS) homopolymer and block co-oligomer mixtures under pressure

The cloud-point temperatures (T-cl's) of poly(ethylene oxide) (PEO) and poly(ethylene oxide)-block-polydimethylsiloxane (P(EO-b-DMS)) homopolymer and block-oligomer mixtures were determined by turbidity measurements over a range of temperatures (105 to 130degrees), pressures (1 to 800 bar), and compositions (10-40 wt.-% PEO). The system phase separates upon cooling and T-cl was found to decrease with an increase in pressure for a constant composition. In the absence of special effects, this finding indicates negative excess volumes. Special attention was paid to the demixing temperatures as a function of the pressure for the different polymer mixtures and the plots in the T-phi plane (where phi signifies volume fractions). The cloud-point curves of the polymer mixture under pressures were observed for different compositions. The Sanchez-Lacombe (SL) lattice fluid theory was used to calculate the spinodals, the binodals, the Flory-Huggins (FH) interaction parameter, the enthalphy of mixing, and the volume changes of mixing. The calculated results show that modified P(EO-b-DMS) scaling parameters with the new combining rules can describe the thermodynamics of the PEO/P(EO-b-DMS) system well with the SL theory.

Effect of Grazing-Mediated Dimethyl Sulfide (DMS) Production on the Swimming Behavior of the Copepod Calanus helgolandicus

Chemical interactions play a fundamental role in the ecology of marine foodwebs. Dimethyl sulfide (DMS) is a ubiquitous marine trace gas that acts as a bioactive compound by eliciting foraging behavior in a range of marine taxa including the copepod Temora longicornis. Production of DMS can rapidly increase following microzooplankton grazing on phytoplankton. Here, we investigated whether grazing-induced DMS elicits an increase in foraging behavior in the copepod Calanus helgolandicus. We developed a semi-automated method to quantify the effect of grazing-mediated DMS on the proportion of the time budget tethered females allocate towards slow swimming, typically associated with feeding. The pooled data showed no differences in the proportion of the 25 min time budget allocated towards slow swimming between high (23.6 +/- 9.74%) and low (29.1 +/- 18.33%) DMS treatments. However, there was a high degree of variability between behavioral responses of individual copepods. We discuss the need for more detailed species-specific studies of individual level responses of copepods to chemical signals at different spatial scales to improve our understanding of chemical interactions between copepods and their prey.