skip to main content


Title: Anisotropic Locally-Conformal Perfectly Matched Layer for Higher Order Curvilinear Finite-Element Modeling
A perfectly matched layer (PML) method is proposed for electrically large curvilinear meshes based on a higher order finite-element modeling paradigm and the concept of transformation electromagnetics. The method maps the non-Maxwellian formulation of the locally conformal PML to a purely Maxwellian implementation using continuously varying anisotropic and inhomogeneous material parameters. An approach to the implementation of a conformal PML for higher order meshes is also presented, based on a method of normal projection for PML mesh generation around an already existing convex volume mesh of a dielectric scatterer, with automatically generated constitutive material parameters. Once the initial mesh is generated, a PML optimization method based on gradient descent is implemented to most accurately match the PML material parameters to the geometrical interface. The numerical results show that the implementation of a conformal PML in the higher order finite-element modeling paradigm dramatically reduces the reflection error when compared to traditional PMLs with piecewise constant material parameters. The ability of the new PML to accurately and efficiently model scatterers with a large variation in geometrical shape and those with complex material compositions is demonstrated in examples of a dielectric almond and a continuously inhomogeneous and anisotropic transformation-optics cloaking structure, respectively.  more » « less
Award ID(s):
1646562
NSF-PAR ID:
10076263
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
IEEE transactions on antennas and propagation
Volume:
65
Issue:
12
ISSN:
0018-926X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Additive manufacturing, no longer reserved exclusively for prototyping components, can create parts with complex geometries and locally tailored properties. For example, multiple homogenous material sources can be used in different regions of a print or be mixed during printing to define properties locally. Additionally, heterogeneous composites provide an opportunity for another level of tuning properties through processing. For example, within particulate-filled polymer matrix composites before curing, the presence of an applied electric and/or magnetic fields can reorient filler particles and form hierarchical structures depending on the fields applied. Control of particle organization is important because effective material properties are highly dependent on the distribution of filler material within composites once cured. While previous work in homogenization and effective medium theories have determined properties based upon ideal analytic distributions of particle orientations and spatial location, this work expands upon these methods generating discrete distributions from quasi-Monte Carlo simulations of the electromagnetic processing event. Results of simulations provide predicted microarchitectures from which effective properties are determined via computational homogenization.

    These particle dynamics simulations account for dielectric and magnetic forces and torques in addition to hydrodynamic forces and hard particle separation. As such, the distributions generated are processing field dependent. The effective properties for a composite represented by this distribution are determined via computational homogenization using finite element analysis (FEA). This provides a path from constituents, through processing parameters to effective material properties. In this work, we use these simulations in conjunction with a multi-objective optimization scheme to resolve the relationships between processing conditions and effective properties, to inform field-assisted additive manufacturing processes.

    The constituent set providing the largest range of properties can be found using optimization techniques applied to the aforementioned simulation framework. This key information provides a recipe for tailoring properties for additive manufacturing design and production. For example, our simulation results show that stiffness for a 10% filler volume fraction can increase by 34% when aligned by an electric field as compared to a randomly distributed composite. The stiffness of this aligned sample is also 29% higher in the direction of the alignment than perpendicular to it, which only differs by 5% from the random case [1]. Understanding this behavior and accurately predicting composite properties is key to producing field processed composites and prints. Material property predictions compare favorably to effective medium theory and experimentation with trends in elastic and magnetic effective properties demonstrating the same anisotropic behavior as a result of applied field processing. This work will address the high computational expense of physics simulation based objective functions by using efficient algorithms and data structures. We will present an optimization framework using nested gradient searches for micro barium hexaferrite particles in a PDMS matrix, optimizing on composite magnetization to determine the volume fraction of filler that will provide the largest range of properties by varying the applied electric and magnetic fields.

     
    more » « less
  2. SUMMARY

    Combining finite element methods for the incompressible Stokes equations with particle-in-cell methods is an important technique in computational geodynamics that has been widely applied in mantle convection, lithosphere dynamics and crustal-scale modelling. In these applications, particles are used to transport along properties of the medium such as the temperature, chemical compositions or other material properties; the particle methods are therefore used to reduce the advection equation to an ordinary differential equation for each particle, resulting in a problem that is simpler to solve than the original equation for which stabilization techniques are necessary to avoid oscillations.

    On the other hand, replacing field-based descriptions by quantities only defined at the locations of particles introduces numerical errors. These errors have previously been investigated, but a complete understanding from both the theoretical and practical sides was so far lacking. In addition, we are not aware of systematic guidance regarding the question of how many particles one needs to choose per mesh cell to achieve a certain accuracy.

    In this paper we modify two existing instantaneous benchmarks and present two new analytic benchmarks for time-dependent incompressible Stokes flow in order to compare the convergence rate and accuracy of various combinations of finite elements, particle advection and particle interpolation methods. Using these benchmarks, we find that in order to retain the optimal accuracy of the finite element formulation, one needs to use a sufficiently accurate particle interpolation algorithm. Additionally, we observe and explain that for our higher-order finite-element methods it is necessary to increase the number of particles per cell as the mesh resolution increases (i.e. as the grid cell size decreases) to avoid a reduction in convergence order.

    Our methods and results allow designing new particle-in-cell methods with specific convergence rates, and also provide guidance for the choice of common building blocks and parameters such as the number of particles per cell. In addition, our new time-dependent benchmark provides a simple test that can be used to compare different implementations, algorithms and for the assessment of new numerical methods for particle interpolation and advection. We provide a reference implementation of this benchmark in aspect (the ‘Advanced Solver for Problems in Earth’s ConvecTion’), an open source code for geodynamic modelling.

     
    more » « less
  3. Abstract The embedded finite element technique provides a unique approach for modeling of fiber-reinforced composites. Meshing fibers as distinct bundles represented by truss elements embedded in a matrix material mesh allows for the assignment of more specific material properties for each component rather than homogenization of all of the properties. However, the implementations of the embedded element technique available in commercial software do not replace the material of the matrix elements with the material of the embedded elements. This causes a redundancy in the volume calculation of the overlapping meshes leading to artificially increased stiffness and mass. This paper investigates the consequences in the energy calculations of an explicit dynamic model due to this redundancy. A method for the correction of the edundancy within a finite element code is suggested which removes extra energy and is shown to be effective at correcting the energy calculations for large amounts of redundant volume. 
    more » « less
  4. Summary

    This work describes a domain embedding technique between two nonmatching meshes used for generating realizations of spatially correlated random fields with applications to large‐scale sampling‐based uncertainty quantification. The goal is to apply the multilevel Monte Carlo (MLMC) method for the quantification of output uncertainties of PDEs with random input coefficients on general and unstructured computational domains. We propose a highly scalable, hierarchical sampling method to generate realizations of a Gaussian random field on a given unstructured mesh by solving a reaction–diffusion PDE with a stochastic right‐hand side. The stochastic PDE is discretized using the mixed finite element method on an embedded domain with a structured mesh, and then, the solution is projected onto the unstructured mesh. This work describes implementation details on how to efficiently transfer data from the structured and unstructured meshes at coarse levels, assuming that this can be done efficiently on the finest level. We investigate the efficiency and parallel scalability of the technique for the scalable generation of Gaussian random fields in three dimensions. An application of the MLMC method is presented for quantifying uncertainties of subsurface flow problems. We demonstrate the scalability of the sampling method with nonmatching mesh embedding, coupled with a parallel forward model problem solver, for large‐scale 3D MLMC simulations with up to 1.9·109unknowns.

     
    more » « less
  5. Abstract

    A fundamental requirement in standard finite element method (FEM) over four‐node quadrilateral meshes is that every element must be convex, else the results can be erroneous. A mesh containing concave element is said to betangled, and tangling can occur, for example, during: mesh generation, mesh morphing, shape optimization, and/or large deformation simulation. The objective of this article is to introduce a tangled finite element method (TFEM) for handling concave elements in four‐node quadrilateral meshes. TFEM extends standard FEM through two concepts. First, the ambiguity of the field in the tangled region is resolved through a careful definition, and this naturally leads to certain correction terms in the FEM stiffness matrix. Second, an equality condition is imposed on the field at re‐entrant nodes of the concave elements. When the correction terms and equality conditions are included, we demonstrate that one can achieve accurate results, and optimal convergence, even over severely tangled meshes. The theoretical properties of the proposed TFEM are established, and the implementation, that requires minimal changes to standard FEM, is discussed in detail. Several numerical experiments are carried out to illustrate the robustness of the proposed method.

     
    more » « less