In this paper we discuss the numerical solution on a simple 2D domain of the Helmoltz equation with mixed boundary conditions. The so called radiation problem depends on the wavenumber constant parameter k and it is inspired here by medical applications, where a transducer emits a pulse at a given frequency. This problem has been successfully solved in the past with the classical Finite Element Method (FEM) for relative small values of k. But in modern applications the values of k can be of order of thousands and FEM faces up several numerical difficulties. To overcome these difficulties we solve the radiation problem using the Isogeometric Analysis (IgA), a kind of generalization of FEM. Starting with the variational formulation of the radiation problem, we show with details how to apply the isogeometric approach in order to compute the coefficients of the approximated solution of radiation problem in terms of the B-spline basis functions. Our implementation of IgA using GeoPDEs software shows that isogeometric approach is superior than FEM, since it is able to reduce substantially the pollution error, especially for high values of k, producing additionally smoother solutions which depend on less degrees of freedom.
翻译:在本文中,我们讨论Helmoltz 方程式简单 2D 域的数值解决方案,并使用混合边界条件。所谓的辐射问题取决于波数常数参数 k, 并在此受到医疗应用的启发, 即一个传感器以给定频率发出脉冲。 这个问题过去曾成功地用传统微量元素法( FEM) 来计算相对小的k。 但在现代应用中, k 的数值可以按千位顺序排列, FEM 面临数方面的困难。 要克服这些困难,我们用Isogeologic 分析(IgA) 解决辐射问题, 这是一种FEM 的概括化。 从辐射问题的变式配方开始, 我们展示了如何应用异度测量方法, 以便用B- spline 基函数来计算辐射问题近似解决办法的系数。 我们使用 GeoPDES软件对IGA 的运用情况表明, 等量法方法优于 FEM, 因为它能够大幅减少污染错误, 特别是高值的K, 产生更多自由度, 以较低程度为依次。