Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAn isogeometric boundary element method for electromagnetic scattering with compatible B-spline discretizations
76
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We outline the construction of compatible B-splines on 3D surfaces that satisfy the continuity requirements for electromagnetic scattering analysis with the boundary element method (method of moments). Our approach makes use of Non-Uniform Rational B-splines to represent model geometry and compatible B-splines to approximate the surface current, and adopts the isogeometric concept in which the basis for analysis is taken directly from CAD (geometry) data. The approach allows for high-order approximations and crucially provides a direct link with CAD data structures that allows for efficient design workflows. After outlining the construction of div- and curl-conforming B-splines defined over 3D surfaces we describe their use with the electric and magnetic field integral equations using a Galerkin formulation. We use Bezier extraction to accelerate the computation of NURBS and B-spline terms and employ H-matrices to provide accelerated computations and memory reduction for the dense matrices that result from the boundary integral discretization. The method is verified using the well known Mie scattering problem posed over a perfectly electrically conducting sphere and the classic NASA almond problem. Finally, we demonstrate the ability of the approach to handle models with complex geometry directly from CAD without mesh generation.
rate research
Read More
Understanding the flow of deformable particles such as liquid drops, synthetic capsules and vesicles, and biological cells confined in a small channel is essential to a wide range of potential chemical and biomedical engineering applications. Computer simulations of this kind of fluid-structure (membrane) interaction in low-Reynolds-number flows raise significant challenges faced by an intricate interplay between flow stresses, complex particles interfacial mechanical properties, and fluidic confinement. Here, we present an isogeometric computational framework by combining the finite-element method (FEM) and boundary-element method (BEM) for an accurate prediction of the deformation and motion of a single soft particle transported in microfluidic channels. The proposed numerical framework is constructed consistently with the isogeometric analysis paradigm; Loops subdivision elements are used not only for the representation of geometry but also for the membrane mechanics solver (FEM) and the interfacial fluid dynamics solver (BEM). We validate our approach by comparison of the simulation results with highly accurate benchmark solutions to two well-known examples available in the literature, namely a liquid drop with constant surface tension in a circular tube and a capsule with a very thin hyperelastic membrane in a square channel. We show that the numerical method exhibits second-order convergence in both time and space. To further demonstrate the accuracy and long-time numerically stable simulations of the algorithm, we perform hydrodynamic computations of a lipid vesicle with bending stiffness and a red blood cell with a composite membrane in capillaries. The present work offers some possibilities to study the deformation behavior of confining soft particles, especially the particles shape transition and dynamics and their rheological signature in channel flows.
In this article, we propose an exponential B-spline collocation method to approximate the solution of the fractional sub-diffusion equation of Caputo type. The present method is generated by use of the Gorenflo-Mainardi-Moretti-Paradisi (GMMP) scheme in time and an efficient exponential B-spline based method in space. The unique solvability is rigorously discussed. Its stability is well illustrated via a procedure closely resembling the classic von Neumann approach. The resulting algebraic system is tri-diagonal that can rapidly be solved by the known algebraic solver with low cost and storage. A series of numerical examples are finally carried out and by contrast to the other algorithms available in the literature, numerical results confirm the validity and superiority of our method.
In this work we present an adaptive boundary element method for computing the electromagnetic response of wave interactions in hyperbolic metamaterials. One unique feature of hyperbolic metamaterial is the strongly directional wave in its propagating cone, which induces sharp transition for the solution of the integral equation across the cone boundary when wave starts to decay or grow exponentially. In order to avoid a global refined mesh over the whole boundary, we employ a two-level a posteriori error estimator and an adaptive mesh refinement procedure to resolve the singularity locally for the solution of the integral equation. Such an adaptive procedure allows for the reduction of the degree of freedom significantly for the integral equation solver while achieving desired accuracy for the solution. In addition, to resolve the fast transition of the fundamental solution and its derivatives accurately across the propagation cone boundary, adaptive numerical quadrature rules are applied to evaluate the integrals for the stiff matrices. Finally, in order to formulate the integral equations over the boundary, we also derive the limits of layer potentials and their derivatives in the hyperbolic media when the target points approach the boundary.
We developed a fast numerical algorithm for solving the three dimensional vectorial Helmholtz equation that arises in electromagnetic scattering problems. The algorithm is based on electric field integral equations and is essentially a boundary element method. Nystroms quadrature rule with a triangular grid is employed to linearize the integral equations, which are then solved by using a right-preconditioned iterative method. We apply the fast multipole technique to accelerate the matrix-vector multiplications in the iterations. We demonstrate the broad applications and accuracy of this method with practical examples including dielectric, plasmonic and metallic objects. We then apply the method to investigate the plasmonic properties of a silver torus and a silver split-ring resonator under the incidence of an electromagnetic plane wave. We show the silver torus can be used as a trapping tool to bind small dielectric or metallic particles.
We develop a numerical method for solving the boundary value problem of The Linear Seventh Ordinary Boundary Value Problem by using seventh degree B-Spline function. Formulation is based on particular terms of order of seventh order boundary value problem. We obtain Septic B-Spline formulation and the Collocation B-spline Method is formulated as an approximation solution. We apply the presented method to solve an example of seventh-order boundary value problem which the results show that there is an agreement between approximate solutions and exact solutions. Resulting low absolute errors show that the presented numerical method is effective for solving high order boundary value problems. Finally, a general conclusion has been included.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
