An efficient ghost-cell immersed boundary method (GCIBM) for simulating turbulent flows in complex geometries is presented. A boundary condition is enforced through a ghost cell method. The reconstruction procedure allows systematic development of numerical schemes for treating the immersed boundary while preserving the overall second-order accuracy of the base solver. Both Dirichlet and Neumann boundary conditions can be treated. The current ghost cell treatment is both suitable for staggered and non-staggered Cartesian grids. The accuracy of the current method is validated using flow past a circular cylinder and large eddy simulation of turbulent flow over a wavy surface. Numerical results are compared with experimental data and boundary-fitted grid results. The method is further extended to an existing ocean model (MITGCM) to simulate geophysical flow over a three-dimensional bump. The method is easily implemented as evidenced by our use of several existing codes.

A Ghost-Cell Immersed Boundary Method for Flow in Complex Geometry

We discuss recent advances on robust unfitted finite element methods on cut meshes. These methods are designed to facilitate computations on complex geometries obtained, for example, from computer- ...

CutFEM: Discretizing geometry and partial differential equations

Thank you very much for downloading mathematical foundations of elasticity. As you may know, people have search numerous times for their favorite books like this mathematical foundations of elasticity, but end up in infectious downloads. Rather than reading a good book with a cup of tea in the afternoon, instead they are facing with some infectious bugs inside their laptop. mathematical foundations of elasticity is available in our book collection an online access to it is set as public so you can download it instantly. Our digital library saves in multiple countries, allowing you to get the most less latency time to download any of our books like this one. Kindly say, the mathematical foundations of elasticity is universally compatible with any devices to read.

Mathematical Foundations Of Elasticity

The finite cell method is an embedded domain method, which combines the fictitious domain approach with higher-order finite elements, adaptive integration, and weak enforcement of unfitted essential boundary conditions. Its core idea is to use a simple unfitted structured mesh of higher-order basis functions for the approximation of the solution fields, while the geometry is captured by means of adaptive quadrature points. This eliminates the need for boundary conforming meshes that require time-consuming and error-prone mesh generation procedures, and opens the door for a seamless integration of very complex geometric models into finite element analysis. At the same time, the finite cell method achieves full accuracy, i.e. optimal rates of convergence, when the mesh is refined, and exponential rates of convergence, when the polynomial degree is increased. Due to the flexibility of the quadrature based geometry approximation, the finite cell method can operate with almost any geometric model, ranging from boundary representations in computer aided geometric design to voxel representations obtained from medical imaging technologies. In this review article, we first provide a concise introduction to the basics of the finite cell method. We then summarize recent developments of the technology, with particular emphasis on the research topics in which we have been actively involved. These include the finite cell method with B-spline and NURBS basis functions, the treatment of geometric nonlinearities for large deformation analysis, the weak enforcement of boundary and coupling conditions, and local refinement schemes. We illustrate the capabilities and advantages of the finite cell method with several challenging examples, e.g. the image-based analysis of foam-like structures, the patient-specific analysis of a human femur bone, the analysis of volumetric structures based on CAD boundary representations, and the isogeometric treatment of trimmed NURBS surfaces. We conclude our review by briefly discussing some key aspects for the efficient implementation of the finite cell method.

The Finite Cell Method: A Review in the Context of Higher-Order Structural Analysis of CAD and Image-Based Geometric Models

The wrinkling and delamination of stiff thin films adhered to a polymer substrate have important applications in “flexible electronics.” The resulting periodic structures, when used for circuitry, have remarkable mechanical properties because stretching or twisting of the substrate is mostly accommodated through bending of the film, which minimizes fatigue or fracture. To date, applications in this context have used substrate patterning to create an anisotropic substrate-film adhesion energy, thereby producing a controlled array of delamination “blisters.” However, even in the absence of such patterning, blisters appear spontaneously, with a characteristic size. Here, we perform well-controlled experiments at macroscopic scales to study what sets the dimensions of these blisters in terms of the material properties and explain our results by using a combination of scaling and analytical methods. Besides pointing to a method for determining the interfacial toughness, our analysis suggests a number of design guidelines for the thin films used in flexible electronic applications. Crucially, we show that, to avoid the possibility that delamination may cause fatigue damage, the thin film thickness must be greater than a critical value, which we determine.

http://absimage.aps.org/image/MAR09/MWS_MAR09-2008-000237.pdf

The macroscopic delamination of thin films from elastic substrates

SUMMARY

We introduce a new method to triangulate planar, curved domains that transforms a specific collection of triangles in a background mesh to conform to the boundary. In the process, no new vertices are introduced, and connectivities of triangles are left unaltered. The method relies on a novel way of parameterizing an immersed boundary over a collection of nearby edges with its closest point projection. To guarantee its robustness, we require that the domain be C2-regular, the background mesh be sufficiently refined near the boundary, and that specific angles in triangles near the boundary be strictly acute. The method can render both straight-edged and curvilinear triangulations for the immersed domain. The latter includes curved triangles that conform exactly to the immersed boundary, and ones constructed with isoparametric mappings to interpolate the boundary at select points. High-order finite elements constructed over these curved triangles achieve optimal accuracy, which has customarily proven difficult in numerical schemes that adopt nonconforming meshes.



Aside from serving as a quick and simple tool for meshing planar curved domains with complex shapes, the method provides significant advantages for simulating problems with moving boundaries and in numerical schemes that require iterating over the geometry of domains. With no conformity requirements, the same background mesh can be adopted to triangulate a large family of domains immersed in it, including ones realized over several updates during the coarse of simulating problems with moving boundaries. We term such a background mesh as a universal mesh for the family of domains it can be used to triangulate. Universal meshes hence facilitate a framework for finite element calculations over evolving domains while using only fixed background meshes. Furthermore, because the evolving geometry can be approximated with any desired order, numerical solutions can be computed with high-order accuracy. We present demonstrative examples using universal meshes to simulate the interaction of rigid bodies with Stokesian fluids. Copyright © 2014 John Wiley & Sons, Ltd.

Universal meshes: A method for triangulating planar curved domains immersed in nonconforming meshes

Summary


We formulate a class of delicately controlled problems to model the kink-free evolution of quasistatic cracks in brittle, isotropic, linearly elastic materials in two dimensions. The evolving crack satisfies familiar principles—Griffith's criterion, local symmetry, and irreversibility. A novel feature of the formulation is that in addition to the crack path, the loading is also treated as an unknown. Specifically, a scaling factor for prescribed Dirichlet and Neumann boundary conditions is computed as part of the solution to yield an always-propagating and ostensibly kink-free crack and a continuous loading history beyond the initial step. A dimensionless statement of the problem depends only on the Poisson's ratio of a homogeneous material, and is in particular, independent of its Young's modulus and fracture toughness.



Numerical resolution of the formulated problem relies on two new ideas. The first is an algorithm to compute triangulations conforming to cracked domains by locally deforming a given background mesh in the vicinity of evolving cracks. The algorithm is robust under mild assumptions on the sizes and angles of triangles near the crack and its smoothness. Hence, a subset, if not the entire family, of cracked domains realized during the course of a simulation can be discretized with the same background mesh; we term the latter a universal mesh for such a family of domains. Universal meshes facilitate adopting a discrete representation for the crack (as splines in our examples), preclude the need for local splitting/retriangulation operations, liberate the crack from following directions prescribed by the mesh, and enable the adoption of standard FEMs to compute the elastic fields in the cracked solid. Second, we employ a method specifically designed to approximate the stress intensity factors for curvilinear cracks. We examine the performance of the resulting numerical method with detailed examples, including comparisons with an exact solution of a crack propagating along a circular arc (which we construct) and comparisons with experimental fracture paths. In all cases, we observe convergence of computed paths, their derivatives, and loading histories with refinement of the universal mesh. Copyright © 2014 John Wiley & Sons, Ltd.

/pdf/simulating-curvilinear-crack-propagation-in-two-dimensions-1qsmnjhmqw.pdf

Simulating curvilinear crack propagation in two dimensions with universal meshes

A discontinuous-Galerkin-based immersed boundary method with non-homogeneous boundary conditions and its application to elasticity

A finite element method to compute three-dimensional equilibrium configurations of fluid membranes

We describe a method for discretizing planar C2-regular domains immersed in non-conforming triangulations. The method consists in constructing mappings from triangles in a background mesh to curvilinear ones that conform exactly to the immersed domain. Constructing such a map relies on a novel way of parameterizing the immersed boundary over a collection of nearby edges with its closest point projection. By interpolating the mappings to curvilinear triangles at select points, we recover isoparametric mappings for the immersed domain defined over the background mesh. Indeed, interpolating the constructed mappings just at the vertices of the background mesh yields a fast meshing algorithm that involves only perturbing a few vertices near the boundary. 
For the discretization of a curved domain to be robust, we have to impose restrictions on the background mesh. Conversely, these restrictions define a family of domains that can be discretized with a given background mesh. We then say that the background mesh is a universal mesh for such a family of domains. The notion of universal meshes is particularly useful in free/moving boundary problems because the same background mesh can serve as the universal mesh for the evolving domain for time intervals that are independent of the time step. Hence it facilitates a framework for finite element calculations over evolving domains while using a fixed background mesh. Furthermore, since the evolving geometry can be approximated with any desired order, numerical solutions can be computed with high-order accuracy. We demonstrate these ideas with various numerical examples.

Ramsharan Rangarajan

Papers

Universal meshes: A method for triangulating planar curved domains immersed in nonconforming meshes

Simulating curvilinear crack propagation in two dimensions with universal meshes

A discontinuous-Galerkin-based immersed boundary method with non-homogeneous boundary conditions and its application to elasticity

A finite element method to compute three-dimensional equilibrium configurations of fluid membranes

Universal Meshes: A new paradigm for computing with nonconforming triangulations