High-accuracy adaptive modeling of the energy distribution of a meniscus-shaped cell culture in a Petri dish

Cylindrical Petri dishes embedded in a rectangular waveguide and exposed to a polarized electromagnetic wave are often used to grow cell cultures. To guarantee the success of these cultures, it is necessary to enforce that the specific absorption rate distribution is sufficiently high and uniform over the Petri dish. Accurate numerical simulations are needed to design such systems. These simulations constitute a challenge due to the strong discontinuity of electromagnetic material properties involved, the relative low field value within the dish cultures compared with the rest of the domain, and the presence of the meniscus shape developed at the liquid boundary. The latter greatly increases the level of complexity of the model in terms of geometry and intensity of the gradients/singularities of the field solution. In here, we employ a three-dimensional (3D) hp-adaptive finite element method using isoparametric elements to obtain highly accurate simulations. We analyze the impact of the geometrical modeling of the meniscus shape cell culture in the hp-adaptivity. Numerical results showing the error convergence history indicate the numerical difficulties arisen due to the presence of a meniscus-shaped object. At the same time, the resulting energy distribution shows that to consider such meniscus shape is essential to guarantee the success of the cell culture from the biological point of view

Goal-oriented self-adaptive hp-strategies for finite element analysis of electromagnetic scattering and radiation problems

In this paper, a fully automatic goal-oriented hp-adaptive finite element strategy for open region electromagnetic problems (radiation and scattering) is presented. The methodology leads to exponential rates of convergence in terms of an upper bound of an user-prescribed quantity of interest. Thus, the adaptivity may be guided to provide an optimal error, not globally for the field in the whole finite element domain, but for specific parameters of engineering interest. For instance, the error on the numerical computation of the S-parameters of an antenna array, the field radiated by an antenna, or the Radar Cross Section on given directions, can be minimized. The efficiency of the approach is illustrated with several numerical simulations with two dimensional problem domains. Results include the comparison with the previously developed energy-norm based hp-adaptivity

An interface between an hp-adaptive finite element package and the pre- and post-processor GiD

An interface between GiD, the interactive graphical user interface used for numerical simulations, developed at the International Center for Numerical Methods in Engineering (CIMNE) of the Universidad Politécnica de Cataluña and the Geometrical Modeling Package (GMP) of the fully automatic hp-adaptive finite element (FE) software, developed at the Institute for Computational Engineering and Sciences (ICES) of the University of Texas at Austin, is presented. GiD is used to construct a tessellation of the problem domain into FE-like regions (blocks in GMP terminology), and the interface obtains and transfers all the topological and geometrical information to GMP. Then, GMP automatically constructs a parameterization for each FE-like region of the GMP mesh, which later can be used to generate the actual FE-mesh and support geometry updates during mesh refinements

Método de Elementos Finitos hp con Adaptabilidad Automática Orientada a un Objetivo para Problemas Abiertos en 2D

In this paper, we describe a fully automatic goaloriented hp-adaptive Finite Element strategy, which is applied to open problems (radiation and scattering). The methodology produces exponential convergence rates in terms of an upper bound of an user-prescribed quantity of interest (in our case, the S-parameter, the far radiated field or far scattering field) against the problem size (number of degrees of freedom). We illustrate the efficiency of the method with 2D numerical simulations of open problems (radiation and scattering). Applications include the far scattering (radiated) field by an object (antenna) and the computation of mutual coupling of the antennas (S-parameters). Results show that self-adaptive goal-oriented hp obtains more accuracy in the quantity of interest than self-adaptive energynorm hp with the same number of degrees of freedom

A comparison between PML, infinite elements and an iterative BEM as mesh truncation methods for HP self-adaptive procedures in electromagnetics

Finite element hp-adaptivity is a technology that allows for very accurate numerical solutions. When applied to open region problems such as radar cross section prediction or antenna analysis, a mesh truncation method needs to be used. This paper compares the following mesh truncation methods in the context of hp-adaptive methods: Infinite Elements, Perfectly Matched Layers and an iterative boundary element based methodology. These methods have been selected because they are exact at the continuous level (a desirable feature required by the extreme accuracy delivered by the hp-adaptive strategy) and they are easy to integrate with the logic of hp-adaptivity. The comparison is mainly based on the number of degrees of freedom needed for each method to achieve a given level of accuracy. Computational times are also included. Two-dimensional examples are used, but the conclusions directly extrapolated to the three dimensional case

High-accuracy adaptive simulations of a Petri dish exposed to electromagnetic radiation

This paper analyses numerically the electric field distribution of a liquid contained in a Petri dish when exposed to electromagnetic waves excited in a rectangular waveguide. Solutions exhibit high-gradients due to the presence of the dielectric liquid contained in the dish. Furthermore, electromagnetic fields within the dielectric have a dramatically lower value than on the remaining part of the domain, which difficults its simulation. Additionally, various singularities of different intensity appear along the boundary of the Petri dish. To properly reproduce and numerically study those effects, we employ a highly-accurate hp-adaptive finite element method. Results of this study demonstrate that the electric field generated within the circular Petri dish is non-homogeneous, and thus, a better shape, size, or location of the dish is needed to achieve an equally distributed radiation enabling the uniform growth of cell cultives

Cálculo de Dosimetría Mediante Elementos Finitos con Adaptabilidad Automática hp en Tres Dimensiones

In this communication the effect of the electromagnetic radiation on in vitro cell cultures is analyzed using a self-adaptive hp-Finite Element Method (hp-FEM) in three dimensions. Computer dosimetry is a challenging problem as it involves complex geometries with high contrast electromagnetic materials. hp-FEM produces exponential convergence rates in terms of the energy-norm error of the solution against the problem size (number of degrees of freedom), even in the presence of singularities. Thus, accurate electromagnetic solutions of complex problems, as computer dosimetry for in vitro experiments, can be obtained. Furthermore, the use of self-adaptive techniques provides solutions, with a user pre-specified degree of accuracy from a initial very coarse mesh without any a priori knowledge of the problem solutio

A three-dimensional self-adaptive hp finite element method for the characterization of waveguide discontinuities

We propose the use of a highly-accurate three-dimensional (3D) fully automatic hp-adaptive finite element method (FEM) for the characterization of rectangular waveguide discontinuities. These discontinuities are either the unavoidable result of mechanical/electrical transitions or deliberately introduced in order to perform certain electrical functions in modern communication systems. The proposed numerical method combines the geometrical flexibility of finite elements with an accuracy that is often superior to that provided by semi-analytical methods. It supports anisotropic refinements on irregular meshes with hanging nodes, and isoparametric elements. It makes use of hexahedral elements compatible with high-order H(curl)H(curl) discretizations. The 3D hp-adaptive FEM is applied for the first time to solve a wide range of 3D waveguide discontinuity problems of microwave communication systems in which exponential convergence of the error is observed

3D hp-adaptive finite element simulations of a magic-T electromagnetic waveguide structure

This paper employs a 3D hp self-adaptive grid-refinement finite element strategy for the solution of a particular electromagnetic waveguide structure known as Magic-T. This structure is utilized as a power divider/combiner in communication systems as well as in other applications. It often incorporates dielectrics, metallic screws, round corners, and so on, which may facilitate its construction or improve its design, but significantly difficult its modeling when employing semi-analytical techniques. The hp-adaptive finite element method enables accurate modeling of a Magic-T structure even in the presence of these undesired materials/geometries. Numerical results demonstrate the suitability of the hp-adaptive method for modeling a Magic-T rectangular waveguide structure, delivering errors below 0.5% with a limited number of unknowns. Solutions of waveguide problems delivered by the self-adaptive hp-FEM are comparable to those obtained with semi-analytical techniques such as the Mode Matching method, for problems where the latest methods can be applied. At the same time, the hp-adaptive FEM enables accurate modeling of more complex waveguide structures

Medidas experimentales de la complejidad computacional de un código autoadaptativo HP para problemas abiertos acelerado mediante ACA

The aim of the novel experimental measures presented in this paper is to show the improvement achieved in the computation time for a 2D self-adaptive hp finite element method (FEM) software accelerated through the Adaptive Cross Approximation (ACA) method. This algebraic method (ACA) was presented in an previous paper in the hp context for the analysis of open region problems, where the robust behaviour, good accuracy and high compression levels of ACA were demonstrated. The truncation of the infinite domain is settled through an iterative computation of the Integral Equation (IE) over a ficticious boundary, which, regardless its accuracy and efficiency, turns out to be the bottelneck of the code. It will be shown that in this context ACA reduces drastically the computational effort of the problem