Downloads & Free Reading Options - Results

http://uf.catalog.fcla.edu/uf.jsp?st=UF003163679&ix=nu&I=0&V=D&pm=1

11

We present a cut finite element method for shape optimization in the case of linear elasticity. The elastic domain is defined by a level-set function, and the evolution of the domain is obtained by moving the level-set along a velocity field using a transport equation. The velocity field is the largest decreasing direction of the shape derivative that satisfies a certain regularity requirement and the computation of the shape derivative is based on a volume formulation. Using the cut finite element method no re--meshing is required when updating the domain and we may also use higher order finite element approximations. To obtain a stable method, stabilization terms are added in the vicinity of the cut elements at the boundary, which provides control of the variation of the solution in the vicinity of the boundary. We implement and illustrate the performance of the method in the two--dimensional case, considering both triangular and quadrilateral meshes as well as finite element spaces of different order.

We consider the model Poisson problem $-\Delta u = f \in \Omega$, $u = g$ on $\pa \Omega$, where $ \Omega $ is a bounded polyhedral domain in $\RR^n$. The objective of the paper is twofold. The first objective is to review the well posedness and the regularity of our model problem using appropriate weighted spaces for the data and the solution. We use these results to derive the domain of the Laplace operator with zero boundary conditions on a concave domain, which seems not to have been fully investigated before. We also mention some extensions of our results to interface problems for the Elasticity equation. The second objective is to illustrate how anisotropic weighted regularity results for the Laplace operator in 3D are used in designing efficient finite element discretizations of elliptic boundary value problems, with the focus on the efficient discretization of the Poisson problem on polyhedral domains in $\RR^3$, following {\em Numer. Funct. Anal. Optim.}, 28(7-8):775--824, 2007. The anisotropic weighted regularity results described and used in the second part of the paper are a consequence of the well-posedness results in (isotropically) weighted Sobolev spaces described in the first part of the paper. The paper is based on the talk by the last named author at the Congress of Romanian Mathematicians, Brasov 2011, and is largely a survey paper.

Most discretizations of the Helmholtz equation result in a system of linear equations that has an indefinite coefficient matrix. Much effort has been put into solving such systems of equations efficiently. In a previous work, the current author and D.C. Dobson proposed a numerical method for solving the complex Helmholtz equation based on the minimization variational principles developed by Milton, Seppecher, and Bouchitte. This method results in a system of equations with a symmetric positive definite coefficient matrix, but at the same time requires solving simultaneously for the solution and its gradient, resulting in a much larger system of equations. Herein is presented a method based on the saddle point variational principles of Milton, Seppecher, and Bouchitte, which produces symmetric positive definite systems of equations, but eliminates the necessity of solving for the gradient of the solution.

In the present work, we propose to extend to the Stokes problem a fictitious domain approach inspired by eXtended Finite Element Method and studied for Poisson problem in [Renard]. The method allows computations in domains whose boundaries do not match. A mixed finite element method is used for fluid flow. The interface between the fluid and the structure is localized by a level-set function. Dirichlet boundary conditions are taken into account using Lagrange multiplier. A stabilization term is introduced to improve the approximation of the normal trace of the Cauchy stress tensor at the interface and avoid the inf-sup condition between the spaces for velocity and the Lagrange multiplier. Convergence analysis is given and several numerical tests are performed to illustrate the capabilities of the method.

Bibliography: l. 128-129

This paper studies the nonconforming Morley finite element approximation of the eigenvalues of the biharmonic operator. A new $C^1$ conforming companion operator leads to an $L^2$ error estimate for the Morley finite element method which directly compares the $L^2$ error with the error in the energy norm and, hence, can dispense with any additional regularity assumptions. Furthermore, the paper presents new eigenvalue error estimates for nonconforming finite elements that bound the error of (possibly multiple or clustered) eigenvalues by the approximation error of the computed invariant subspace. An application is the proof of optimal convergence rates for the adaptive Morley finite element method for eigenvalue clusters.

Rigorous computer simulations of propagating electromagnetic fields have become an important tool for optical metrology and optics design of nanostructured components. As has been shown in previous benchmarks some of the presently used methods suffer from low convergence rates and/or low accuracy of the results and exhibit very long computation times which makes application to extended 2D layout patterns impractical. We address 3D simulation tasks by using a finite-element solver which has been shown to be superior to competing methods by several orders of magnitude in accuracy and computational time for typical microlithography simulations. We report on the current status of the solver, incorporating higher order edge elements, adaptive refinement methods, and fast solution algorithms. Further, we investigate the performance of the solver in the 3D simulation project of light diffraction off an alternating phase-shift contact-hole mask.

In this paper we consider the numerical approximation of a general second order semi-linear parabolic partial differential equation. Equations of this type arise in many contexts, such as transport in porous media. Using finite element method for space discretization and the exponential Rosenbrock-Euler method for time discretization, we provide a rigorous convergence proof in space and time under only the standard Lipschitz condition of the nonlinear part for both smooth and nonsmooth initial solution. This is in contrast to very restrictive assumptions made in the literature, where the authors have considered only approximation in time so far in their convergence proofs. The optimal orders of convergence in space and in time are achieved for smooth and nonsmooth initial solution.

In contrast to the commonly used lattice Boltzmann method, off-lattice Boltzmann methods decouple the velocity discretization from the underlying spatial grid, thus allowing for more efficient geometric representations of complex boundaries. The current work combines characteristic-based integration of the streaming step with the free-energy based multiphase model by Lee et. al. [Journal of Computational Physics, 206 (1), 2005 ]. This allows for simulation time steps more than an order of magnitude larger than the relaxation time. Unlike previous work by Wardle et. al. [Computers and Mathematics with Applications, 65 (2), 2013 ] that integrated intermolecular forcing terms in the advection term, the current scheme applies collision and forcing terms locally for a simpler finite element formulation. A series of thorough benchmark studies reveal that this does not compromise stability and that the scheme is able to accurately simulate flows at large density and viscosity contrasts.

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 point-wise 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.

For the stationary advection-diffusion problem the standard continuous Galerkin method is unstable without some additional control on the mesh or method. The interior penalty discontinuous Galerkin method is stable but at the expense of an increased number of degrees of freedom. The hybrid method proposed in [5] combines the computational complexity of the continuous method with the stability of the discontinuous method without a significant increase in degrees of freedom. We discuss the implementation of this method using the finite element library deal.ii and present some numerical experiments.

Final report for period October 1982 - September 1983.--pt.1

The search for more efficient gas turbine engines has led to the proposal for the replacement of metal high temperature components with ceramic components . Essential to this effort is the numerical analysis of proposed designs . This thesis report describes the model discretization of a proposed blade design, the development of pre- and post-processors for the ADINA finite element code and the initial stress analysis of the modeled blade.

The research accomplishments in this project include: parallel implementation of the p-version, development of multi-p methods and multi-p conditioning, study of the conditioning of global stiffness matrices, investigation of the problem of model selection, numerical analysis of material interface singularities in two dimensions, and extraction of generalized stress intensity factors using the principle of complementary energy. (AN)

An axisymmetric boundary element flow algorithm is coupled with a finite element structural analyzer to perform interactive calculations of the growth and subsequent collapse of an explosion bubble near a submerged compliant structure. The validity of the program is established through direct comparisons with published experimental data. A parametric study of the interaction between a growing and collapsing bubble and a spherical shell is presented. The results show that if the stiffness of the shell is sufficiently low, then the mass of the shell is a critical parameter in the collapse problem. If the mass of the structure is high, a reentrant jet forms and is then directed towards the shell. As the mass of the spherical shell is decreased, the collapse becomes spherical with no jet formation. At the lowest structural mass for which calculations are performed, a jet directed away from the structure begins to form. The ratio of the depth of submergence to bubble maximum radius was also found to be a critical parameter in the collapse problem. When this ratio is large (greater than 100), the collapse is driven by interaction forces. However, for shallow submergence, buoyancy effects become more important than interaction forces. Fluid-Structure Interaction, Boundary Element, Fluid Loading, Bubble Collapse, Numerical Prediction Method, Whipping, Finite Element, Cavitation.

A fatigue damage computational algorithm utilizing a multiaxial, isothermal, continuum based fatigue damage model for unidirectional metal matrix composites has been implemented into the commercial finite element code MARC using MARC user subroutines. Damage is introduced into the finite element solution through the concept of effective stress which fully couples the fatigue damage calculations with the finite element deformation solution. An axisymmetric stress analysis was performed on a circumferentially reinforced ring, wherein both the matrix cladding and the composite core were assumed to behave elastic-perfectly plastic. The composite core behavior was represented using Hill's anisotropic continuum based plasticity model, and similarly, the matrix cladding was represented by an isotropic plasticity model. Results are presented in the form of S-N curves and damage distribution plots.

Aim: The aim of this study was to compare plate positions in Lapidus arthrodesis in terms of displacement and plate stress.

An investigation was made of the buckling of rectangular plates with free, simply supported, or clamped edges and loaded in uniform compression or pure shear. The finite element method with a rectangular, sixteen-degree-of -freedom element was used. Results show that, for the same order of error, this element allows much coarser plate divisions than a widely used twelve-degree-of-freedom element. An order of magnitude reduction in the required order of the assembled stiffness and stability coefficient matrices results.

The theory of the computable a-posteriori error estimate for a finite element method is developed. Among other things, it is shown that the error estimate is very reliable and the ratio (called effectivity index) between the estimator and the true error approaches one. Numerical examples computed by program FEARS (Finite Element Adaptive Research Solver) of the University of Maryland, illustrate the effectivity and reliability of the estimators. (Author)

The theory of the computable a-posteriori error estimate for a finite element method is developed. Among other things, it is shown that the error estimate is very reliable and the ratio (called effectivity index) between the estimator and the true error approaches one. Numerical examples computed by program FEARS (Finite Element Adaptive Research Solver) of the University of Maryland, illustrate the effectivity and reliability of the estimators. (Author)

There are three versions of the finite element method. The classical h- version achieves accuracy by refining the mesh while using low degrees p of the elements, usually p = 1, 2. The p-version keeps the mesh fixed and the accuracy is achieved by increasing the degree p. The h-p version combines both approaches. Although the p and h-p versions of the finite element method are relatively new developments, many basic results are available. The aim of this paper is to give a survey of the basic available results and directions for further development. The paper tries to survey basic theoretical implementational and computational aspects of the method as of today.

The nonlinear analysis of an orthotropic, laminated shell of revolution with transverse shear deformations is presented. The finite element method employing a curved shell element was the method used to perform the analysis. Each element has two nodes with four degrees of freedom at each node, viz., two translation, one bending rotation and one transverse shear rotation. The classical Kirchhoff-Love assumption for normals to the midsurface was relaxed in favor of the shear deformation mode. A curved shell element was developed that matches slopes and curvatures as well as displacements at the nodal circles. This significantly reduces the meridional bending moments present at element junctures as compared to the case in which straight line elements (conical frusta) are used to represent a shell of revolution having meridional curvatures. A computer program was written in Fortran IV to implement the theory. Several example problems, both linear and nonlinear, were solved and the results compared with solutions from the literature.

The primary aim of this study is to develop a finite element model for elastic scattering by axisymmetric bodies submerged in a heterogeneous ocean acoustic waveguide.

The rate of convergence of the finite element method is a function of the strategy by which the number of degrees of freedom are increased. Alternative strategies are examined in the light of recent theoretical results and computational experience. (Author)

In the past two decades computational structural analysis and design have developed into very extensive subject areas. A major contributing factor certainly is the evolution of the finite element method into a powerful tool for the solution of a broad range of problems. In this article we address two such questions which may be broadly phrased as follows: (a) How can we calculate, at reasonable cost, reliable estimates of the error of computed finite element solutions? (b) How should we set up the overall solution process so as to achieve, with reasonable certainty and efficiency, a solution for which the accuracy either falls into a prescribed range or is best possible for the allotted computational cost?

AUTOMESH-2D is a computer program specifically designed as a preprocessor for the scattering analysis of two dimensional bodies by the finite element method. This program was developed due to a need for reproducing the effort required to define and check the geometry data, element topology, and material properties. There are six modules in the program: (1) Parameter Specification; (2) Data Input; (3) Node Generation; (4) Element Generation; (5) Mesh Smoothing; and (5) Data File Generation.

ADA038759

ADA038759

A study was undertaken to investigate the measurement of wing deformation and internal loads using measured strain data. Future aerospace vehicle research depends on the ability to accurately measure the deformation and internal loads during ground testing and in flight. The approach uses the inverse Finite Element Method (iFEM). The iFEM is a robust, computationally efficient method that is well suited for real-time measurement of real-time structural deformation and loads. The method has been validated in previous work, but has yet to be applied to a large-scale test article. This work is in preparation for an upcoming loads test of a half-span test wing in the Flight Loads Laboratory at the National Aeronautics and Space Administration Armstrong Flight Research Center (Edwards, California). The method has been implemented into an efficient MATLAB® (The MathWorks, Inc., Natick, Massachusetts) code for testing different sensor configurations. This report discusses formulation and implementation along with the preliminary results from a representative aerospace structure. The end goal is to investigate the modeling and sensor placement approach so that the best practices can be applied to future aerospace projects.

The least squares (L sub 2) finite element method is introduced for 2-D steady state pure convection problems with smooth solutions. It is proven that the L sub 2 method has the same stability estimate as the original equation, i.e., the L sub 2 method has better control of the streamline derivative. Numerical convergence rates are given to show that the L sub 2 method is almost optimal. This L sub 2 method was then used as a framework to develop an iteratively reweighted L sub 2 finite element method to obtain a least absolute residual (L sub 1) solution for problems with discontinuous solutions. This L sub 1 finite element method produces a nonoscillatory, nondiffusive and highly accurate numerical solution that has a sharp discontinuity in one element on both coarse and fine meshes. A robust reweighting strategy was also devised to obtain the L sub 1 solution in a few iterations. A number of examples solved by using triangle and bilinear elements are presented.

The least squares (L sub 2) finite element method is introduced for 2-D steady state pure convection problems with smooth solutions. It is proven that the L sub 2 method has the same stability estimate as the original equation, i.e., the L sub 2 method has better control of the streamline derivative. Numerical convergence rates are given to show that the L sub 2 method is almost optimal. This L sub 2 method was then used as a framework to develop an iteratively reweighted L sub 2 finite element method to obtain a least absolute residual (L sub 1) solution for problems with discontinuous solutions. This L sub 1 finite element method produces a nonoscillatory, nondiffusive and highly accurate numerical solution that has a sharp discontinuity in one element on both coarse and fine meshes. A robust reweighting strategy was also devised to obtain the L sub 1 solution in a few iterations. A number of examples solved by using triangle and bilinear elements are presented.

The least squares (L sub 2) finite element method is introduced for 2-D steady state pure convection problems with smooth solutions. It is proven that the L sub 2 method has the same stability estimate as the original equation, i.e., the L sub 2 method has better control of the streamline derivative. Numerical convergence rates are given to show that the L sub 2 method is almost optimal. This L sub 2 method was then used as a framework to develop an iteratively reweighted L sub 2 finite element method to obtain a least absolute residual (L sub 1) solution for problems with discontinuous solutions. This L sub 1 finite element method produces a nonoscillatory, nondiffusive and highly accurate numerical solution that has a sharp discontinuity in one element on both coarse and fine meshes. A robust reweighting strategy was also devised to obtain the L sub 1 solution in a few iterations. A number of examples solved by using triangle and bilinear elements are presented.

The least squares (L sub 2) finite element method is introduced for 2-D steady state pure convection problems with smooth solutions. It is proven that the L sub 2 method has the same stability estimate as the original equation, i.e., the L sub 2 method has better control of the streamline derivative. Numerical convergence rates are given to show that the L sub 2 method is almost optimal. This L sub 2 method was then used as a framework to develop an iteratively reweighted L sub 2 finite element method to obtain a least absolute residual (L sub 1) solution for problems with discontinuous solutions. This L sub 1 finite element method produces a nonoscillatory, nondiffusive and highly accurate numerical solution that has a sharp discontinuity in one element on both coarse and fine meshes. A robust reweighting strategy was also devised to obtain the L sub 1 solution in a few iterations. A number of examples solved by using triangle and bilinear elements are presented.

Pustaka Online Punyo Den Free

We present a locally adapted parametric finite element method for interface problems. For this adapted finite element method we show optimal convergence for elliptic interface problems with a discontinuous diffusion parameter. The method is based on the adaption of macro elements where a local basis represents the interface. The macro elements are independent of the interface and can be cut by the interface. A macro element which is a triangle in the triangulation is divided into four subtriangles. On these subtriangles, the basis functions of the macro element are interpreted as linear functions. The position of the vertices of these subtriangles is determined by the location of the interface in the case a macro element is cut by the interface. Quadrature is performed on the subtriangles via transformations to a reference element. Due to the locality of the method, its use is well suited on distributed architectures.

The goal of this paper is to introduce a simple finite element method to solve the Stokes and the Navier-Stokes equations. This method is in primal velocity-pressure formulation and is so simple such that both velocity and pressure are approximated by piecewise constant functions. Implementation issues as well as error analysis are investigated. A basis for a divergence free subspace of the velocity field is constructed so that the original saddle point problem can be reduced to a symmetric and positive definite system with much fewer unknowns. The numerical experiments indicate that the method is accurate and robust.

We consider the time-dependent Gross-Pitaevskii equation describing the dynamics of rotating Bose-Einstein condensates and its discretization with the finite element method. We analyze a mass conserving Crank-Nicolson-type discretization and prove corresponding a priori error estimates with respect to the maximum norm in time and the $L^2$- and energy-norm in space. The estimates show that we obtain optimal convergence rates under the assumption of additional regularity for the solution to the Gross-Pitaevskii equation. We demonstrate the performance of the method in numerical experiments.

We present a robust computational framework for advective-diffusive-reactive systems that satisfies maximum principles, the non-negative constraint, and element-wise species balance property. The proposed methodology is valid on general computational grids, can handle heterogeneous anisotropic media, and provides accurate numerical solutions even for very high P\'eclet numbers. The significant contribution of this paper is to incorporate advection (which makes the spatial part of the differential operator non-self-adjoint) into the non-negative computational framework, and overcome numerical challenges associated with advection. We employ low-order mixed finite element formulations based on least-squares formalism, and enforce explicit constraints on the discrete problem to meet the desired properties. The resulting constrained discrete problem belongs to convex quadratic programming for which a unique solution exists. Maximum principles and the non-negative constraint give rise to bound constraints while element-wise species balance gives rise to equality constraints. The resulting convex quadratic programming problems are solved using an interior-point algorithm. Several numerical results pertaining to advection-dominated problems are presented to illustrate the robustness, convergence, and the overall performance of the proposed computational framework.

We consider the P1/P1 or P1b/P1 finite element approximations to the Stokes equations in a bounded smooth domain subject to the slip boundary condition. A penalty method is applied to address the essential boundary condition $u\cdot n = g$ on $\partial\Omega$, which avoids a variational crime and simultaneously facilitates the numerical implementation. We give $O(h^{1/2} + \epsilon^{1/2} + h/\epsilon^{1/2})$-error estimate for velocity and pressure in the energy norm, where $h$ and $\epsilon$ denote the discretization parameter and the penalty parameter, respectively. In the two-dimensional case, it is improved to $O(h + \epsilon^{1/2} + h^2/\epsilon^{1/2})$ by applying reduced-order numerical integration to the penalty term. The theoretical results are confirmed by numerical experiments.

In this paper we consider the initial-boundary value problem for the time-dependent Maxwell-Schr\"{o}dinger equations, which arises in the interaction between the matter and the electromagnetic field for the semiconductor quantum devices. A Crank-Nicolson finite element method for solving the problem is presented. The optimal energy-norm error estimates for the numerical algorithm without any time-step restrictions are derived. Numerical tests are then carried out to confirm the theoretical results.

Fluid flow applications can involve a number of coupled problems. One is the simulation of free-surface flows, which require the solution of a free-boundary problem. Within this problem, the governing equations of fluid flow are coupled with a domain deformation approach. This work reviews five of those approaches: interface tracking using a boundary-conforming mesh and, in the interface capturing context, the level-set method, the volume-of-fluid method, particle methods, as well as the phase-field method. The history of each method is presented in combination with the most recent developments in the field. Particularly, the topics of extended finite elements (XFEM) and NURBS-based methods, such as Isogeometric Analysis (IGA), are addressed. For illustration purposes, two applications have been chosen: two-phase flow involving drops or bubbles and sloshing tanks. The challenges of these applications, such as the geometrically correct representation of the free surface or the incorporation of surface tension forces, are discussed.

In this paper, we replace the asymptotic enrichments around the crack tip in the extended finite element method (XFEM) with the semi-analytical solution obtained by the scaled boundary finite element method (SBFEM). The proposed method does not require special numerical integration technique to compute the stiffness matrix and it improves the capability of the XFEM to model cracks in homogeneous and/or heterogeneous materials without a priori knowledge of the asymptotic solutions. A heaviside enrichment is used to represent the jump across the discontinuity surface. We call the method as the extended scaled boundary finite element method (xSBFEM). Numerical results presented for a few benchmark problems in the context of linear elastic fracture mechanics show that the proposed method yields accurate results with improved condition number. A simple MATLAB code is annexed to compute the terms in the stiffness matrix, which can easily be integrated in any existing FEM/XFEM code.

An adaptive finite element method is presented for the elastic scattering of a time-harmonic plane wave by a periodic surface. First, the unbounded physical domain is truncated into a bounded computational domain by introducing the perfectly matched layer (PML) technique. The well-posedness and exponential convergence of the solution are established for the truncated PML problem by developing an equivalent transparent boundary condition. Second, an a posteriori error estimate is deduced for the discrete problem and is used to determine the finite elements for refinements and to determine the PML parameters. Numerical experiments are included to demonstrate the competitive behavior of the proposed adaptive method.

Report for period April 1984 - June 1984.--pt.2

Final report for period October 1982 - September 1983.--pt.1

Bibliography: l. 109

This article is from Journal of the Royal Society Interface , volume 11 . Abstract Significant advances in the functional analysis of musculoskeletal systems require the development of modelling techniques with improved focus, accuracy and validity. This need is particularly visible in the fields, such as palaeontology, where unobservable parameters may lie at the heart of the most interesting research questions, and where models and simulations may provide some of the most innovative solutions. Here, we report on the development of a computational modelling method to generate estimates of the mechanical properties of vertebral bone across two living species, using elderly human and juvenile porcine specimens as cases with very different levels of bone volume fraction and mineralization. This study is presented in two parts; part I presents the computational model development and validation, and part II the virtual loading regime and results. This work paves the way for the future estimation of mechanical properties in fossil mammalian bone.