The Influence of Quadrature Errors on Isogeometric Mortar Methods

Mortar methods have recently been shown to be well suited for isogeometric analysis. We review the recent mathematical analysis and then investigate the variational crime introduced by quadrature formulas for the coupling integrals. Motivated by finite element observations, we consider a quadrature rule purely based on the slave mesh as well as a method using quadrature rules based on the slave mesh and on the master mesh, resulting in a non-symmetric saddle point problem. While in the first case reduced convergence rates can be observed, in the second case the influence of the variational crime is less significant

Orbital HP-Clouds for Solving Schrödinger Equation in Quantum Mechanics

Solving Schroedinger equation in quantum mechanics presents a challenging task in numerical methods due to the high order behavior and high dimension characteristics in the wave functions, in addition to the highly coupled nature between wave functions. This work introduces orbital and polynomial enrichment functions to the partition of unity for solution of Schroedinger equation under the framework of HP-Clouds. An intrinsic enrichment of orbital function and extrinsic enrichment of monomial functions are proposed. Due to the employment of higher order basis functions, a higher order stabilized conforming nodal integration is developed. The proposed methods are implemented using the density functional theory for solution of Schroedinger equation. Analysis of several single and multi-electron/nucleus structures demonstrates the effectiveness of the proposed method

A parameter-free total Lagrangian smooth particle hydrodynamics algorithm applied to problems with free surfaces

This paper presents a new Smooth Particle Hydrodynamics computational framework for the solution of inviscid free surface flow problems. The formulation is based on the Total Lagrangian description of a system of first-order conservation laws written in terms of the linear momentum and the Jacobian of the deformation. One of the aims of this paper is to explore the use of Total Lagrangian description in the case of large deformations but without topological changes. In this case, the evaluation of spatial integrals is carried out with respect to the initial undeformed configuration, yielding an extremely efficient formulation where the need for continuous particle neighbouring search is completely circumvented. To guarantee stability from the SPH discretisation point of view, consistently derived Riemann-based numerical dissipation is suitably introduced where global numerical entropy production is demonstrated via a novel technique in terms of the time rate of the Hamiltonian of the system. Since the kernel derivatives presented in this work are fixed in the reference configuration, the non-physical clumping mechanism is completely removed. To fulfil conservation of the global angular momentum, a posteriori (least-squares) projection procedure is introduced. Finally, a wide spectrum of dedicated prototype problems is thoroughly examined. Through these tests, the SPH methodology overcomes by construction a number of persistent numerical drawbacks (e.g. hour-glassing, pressure instability, global conservation and/or completeness issues) commonly found in SPH literature, without resorting to the use of any ad-hoc user-defined artificial stabilisation parameters. Crucially, the overall SPH algorithm yields equal second order of convergence for both velocities and pressure

The Viscoelastic Properties of Passive Eye Muscle in Primates. I: Static Forces and Step Responses

The viscoelastic properties of passive eye muscles are prime determinants of the deficits observed following eye muscle paralysis, the root cause of several types of strabismus. Our limited knowledge about such properties is hindering the ability of eye plant models to assist in formulating a patient's diagnosis and prognosis. To investigate these properties we conducted an extensive in vivo study of the mechanics of passive eye muscles in deeply anesthetized and paralyzed monkeys. We describe here the static length-tension relationship and the transient forces elicited by small step-like elongations. We found that the static force increases nonlinearly with length, as previously shown. As expected, an elongation step induces a fast rise in force, followed by a prolonged decay. The time course of the decay is however considerably more complex than previously thought, indicating the presence of several relaxation processes, with time constants ranging from 1 ms to at least 40 s. The mechanical properties of passive eye muscles are thus similar to those of many other biological passive tissues. Eye plant models, which for lack of data had to rely on (erroneous) assumptions, will have to be updated to incorporate these properties

Maximum-Entropy Meshfree Method for Compressible and Near-Incompressible Elasticity

Numerical integration errors and volumetric locking in the near-incompressible limit are two outstanding issues in Galerkin-based meshfree computations. In this paper, we present a modified Gaussian integration scheme on background cells for meshfree methods that alleviates errors in numerical integration and ensures patch test satisfaction to machine precision. Secondly, a lockingfree small-strain elasticity formulation for meshfree methods is proposed, which draws on developments in assumed strain methods and nodal integration techniques. In this study, maximum-entropy basis functions are used; however, the generality of our approach permits the use of any meshfree approximation. Various benchmark problems in two-dimensional compressible and near-incompressible small strain elasticity are presented to demonstrate the accuracy and optimal convergence in the energy norm of the maximumentropy meshfree formulation

Improved robustness for nearly-incompressible large deformation meshfree simulations on Delaunay tessellations

Artículo de publicación ISIA displacement-based Galerkin meshfree method for large deformation analysis of nearly-incompressible elastic solids is presented. Nodal discretization of the domain is defined by a Delaunay tessellation (three-node triangles and four-node tetrahedra), which is used to form the meshfree basis functions and to numerically integrate the weak form integrals. In the proposed approach for nearly-incompressible solids, a volume-averaged nodal projection operator is constructed to average the dilatational constraint at a node from the displacement field of surrounding nodes. The nodal dilatational constraint is then projected onto the linear approximation space. The displacement field is constructed on the linear space and enriched with bubble-like meshfree basis functions for stability. The new procedure leads to a displacement-based formulation that is similar to F-bar methodologies in finite elements and isogeometric analysis. We adopt maximum-entropy meshfree basis functions, and the performance of the meshfree method is demonstrated on benchmark problems using structured and unstructured background meshes in two and three dimensions. The nonlinear simulations reveal that the proposed methodology provides improved robustness for nearlyincompressible large deformation analysis on Delaunay meshes.Chilean National Fund for Scientific and Technological Development (Fondecyt) 11110389 U.S. National Science Foundation CMMI-133478