Bounding the influence of domain parameterization and knot spacing on numerical stability in isogeometric analysis. (English) Zbl 1295.65116

Pilgerstorfer, Elisabeth; Jüttler, Bert

Comput. Methods Appl. Mech. Eng. 268, 589-613 (2014).

Summary: Isogeometric Analysis (IGA) was introduced by T. J. R. Hughes et al. [Comput. Methods Appl. Mech. Eng. 194, No. 39–41, 4135–4195 (2005; Zbl 1151.74419)] as a new method to bridge the gap between the geometry description and numerical analysis. Similar to the finite element approach, the IGA concept to solve a partial differential equation leads to a (linear) system of equations. The condition number of the coefficient matrix is a crucial factor for the stability of the system. It depends strongly on the domain parameterization, which provides the isogeometric discretization. In this paper we derive a bound for the condition number of the stiffness matrix of the Poisson equation. In particular, we investigate the influence of the domain parameterization and the knot spacing on the stability of the numerical system. The factors appearing in our bound reflect the stability properties of a given domain parameterization.

MSC:

65N30	Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs
65D17	Computer-aided design (modeling of curves and surfaces)

Keywords:

numerical stability; isogeometric analysis; domain parameterization; knot spacing; condition number; quality measure

Citations:

Zbl 1151.74419

Software:

FEAPpv; ISOGAT

References:

[1]	Hughes, T. J.R.; Cottrell, J. A.; Bazilevs, Y., Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement, Computer Methods in Applied Mechanics and Engineering, 194, 4135-4195 (2005) · Zbl 1151.74419
[2]	Bazilevs, Y.; Beirão da Veiga, L.; Cottrell, J. A.; Hughes, T. J.R.; Sangalli, G., Isogeometric analysis: approximation, stability and error estimates for h-refined meshes, Mathematical Models and Methods in Applied Sciences, 16, 1031-1090 (2006) · Zbl 1103.65113
[3]	Beirão da Veiga, L.; Buffa, A.; Rivas, J.; Sangalli, G., Some estimates for h-p-k-refinement in isogeometric analysis, Numerische Mathematik, 118, 271-305 (2011) · Zbl 1222.41010
[4]	Beirão da Veiga, L.; Cho, D.; Sangalli, G., Anisotropic NURBS approximation in isogeometric analysis, Computer Methods in Applied Mechanics and Engineering, 209-212, 1-11 (2012) · Zbl 1243.65027
[5]	Hughes, T. J.R.; Reali, A.; Sangalli, G., Efficient quadrature for NURBS-based isogeometric analysis, Computer Methods in Applied Mechanics and Engineering, 199, 301-313 (2010) · Zbl 1227.65029
[6]	Strang, W. G.; Fix, G., An Analysis of the Finite Element Method (1973), Prentice-Hall: Prentice-Hall New Jersey · Zbl 0278.65116
[7]	Zienkiewicz, O. C.; Taylor, R. L.; Zhu, J. Z., The Finite Element Method: Its Basis and Fundamentals (2005), Elsevier Science: Elsevier Science Oxford · Zbl 1307.74005
[8]	Brenner, S. C.; Scott, R. L., The Mathematical Theory of Finite Element Methods (2005), Springer
[9]	Bank, R. E.; Scott, L. R., On the conditioning of finite element equations with highly refined meshes, SIAM Journal on Numerical Analysis, 26, 1383-1394 (1989) · Zbl 0688.65062
[10]	Hu, N.; Guo, X.-Z.; Katz, I. N., Bounds for eigenvalues and condition numbers in the p-version of the finite element method, Mathematics of Computation, 67, 1423-1450 (1998) · Zbl 0907.65112
[11]	Du, Q.; Wang, D.; Zhu, L., On mesh geometry and stiffness matrix conditioning for general finite element spaces, SIAM Journal on Numerical Analysis, 47, 1421-1444 (2009) · Zbl 1191.65152
[12]	Aigner, M.; Heinrich, C.; Jüttler, B.; Pilgerstorfer, E.; Simeon, B.; Vuong, A.-V., Swept volume parameterization for isogeometric analysis, (Hancock, E.; Martin, R., The Mathematics of Surfaces XIII. The Mathematics of Surfaces XIII, Lecture Notes in Computer Science, vol. 5654 (2009), Springer), 19-44 · Zbl 1253.65182
[13]	Cohen, E.; Martin, T.; Kirby, R. M.; Lyche, T.; Riesenfeld, R. F., Analysis-aware modeling: understanding quality considerations in modeling for isogeometric analysis, Computer Methods in Applied Mechanics and Engineering, 199, 334-356 (2010) · Zbl 1227.74109
[14]	Martin, T.; Cohen, E.; Kirby, R. M., Volumetric parameterization and trivariate B-spline fitting using harmonic functions, Computer Aided Geometric Design, 26, 648-664 (2009) · Zbl 1205.65094
[15]	Xu, G.; Mourrain, B.; Duvigneau, R.; Galligo, A., Optimal analysis-aware parameterization of computational domain in isogeometric analysis, Proceedings of Geometric Modeling and Processing, 236-254 (2010)
[16]	Xu, G.; Mourrain, B.; Duvigneau, R.; Galligo, A., Parameterization of computational domain in isogeometric analysis: methods and comparison, Computer Methods in Applied Mechanics and Engineering, 200, 2021-2031 (2011) · Zbl 1228.65232
[17]	Xu, G.; Mourrain, B.; Duvigneau, R.; Galligo, A., Analysis-suitable volume parameterization of multi-block computational domain in isogeometric applications, Computer-Aided Design, 45, 395-404 (2013)
[18]	Takacs, T.; Jüttler, B., Existence of stiffness matrix integrals for singularly parameterized domains in isogeometric analysis, Computer Methods in Applied Mechanics and Engineering, 200, 3568-3582 (2011) · Zbl 1239.65014
[19]	Takacs, T.; Jüttler, B., \(H^2\) regularity properties of singular parameterizations in isogeometric analysis, Graphical Models, 74, 361-372 (2012)
[21]	Nguyen, D. M.; Evgrafov, A.; Gersborg, A. R.; Gravesen, J., Isogeometric shape optimization of vibrating membranes, Computer Methods in Applied Mechanics and Engineering, 200, 1343-1353 (2011) · Zbl 1228.74062
[22]	Nguyen, D. M.; Evgrafov, A.; Gravesen, J., Isogeometric shape optimization for electromagnetic scattering problems, Progress in Electromagnetic Research (PIER) B, 45, 117-146 (2012)
[23]	Qian, X., Full analytical sensitivities in NURBS based isogeometric shape optimization, Computer Methods in Applied Mechanics and Engineering, 199, 2059-2071 (2010) · Zbl 1231.74352
[24]	Wall, W. A.; Frenzel, M. A.; Cyron, C., Isogeometric structural shape optimization, Computer Methods in Applied Mechanics and Engineering, 197, 2976-2988 (2008) · Zbl 1194.74263
[25]	Seo, Y.-D.; Kim, H.-J.; Youn, S.-K., Shape optimization and its extension to topological design based on isogeometric analysis, International Journal of Solids and Structures, 47, 1618-1640 (2010) · Zbl 1194.74275
[26]	Pospíšilová, A.; Lepš, M.; Rypl, D.; Patzák, B., Shape optimization using isogeometric analysis and particle swarm optimization, (Topping, B. H.V., Proceedings of the 11th International Conference on Computational Structures Technology (2012), Civil-Comp Press: Civil-Comp Press Stirlingshire, UK), (Paper 220)
[27]	Lipton, S.; Evans, J. A.; Bazilevs, Y.; Elguedj, T.; Hughes, T. J.R., Robustness of isogeometric structural discretizations under severe mesh distortion, Computer Methods in Applied Mechanics and Engineering, 199, 357-373 (2010) · Zbl 1227.74112
[29]	de Boor, C., A Practical Guide to Splines (1978), Springer Verlag: Springer Verlag New York · Zbl 0406.41003
[30]	Prautzsch, H.; Boehm, W.; Paluszny, M., Bézier and B-Spline Techniques (2002), Springer: Springer New York · Zbl 1033.65008
[31]	Hoschek, J.; Lasser, D., Fundamentals of Computer Aided Geometric Design (1993), A.K. Peters Ltd.: A.K. Peters Ltd. Natick · Zbl 0788.68002
[32]	Piegl, L.; Tiller, W., The NURBS Book (1995), Springer: Springer London · Zbl 0828.68118
[33]	Cottrell, J. A.; Hughes, T. J.R.; Bazilevs, Y., Isogeometric Analysis. Towards Integration of CAD and FEA (2009), Wiley & Sons Ltd.: Wiley & Sons Ltd. Chichester · Zbl 1378.65009
[34]	Ikebe, Y.; Inagaki, T.; Miyamoto, S., The monotonicity theorem. Cauchy’s interlace theorem, and the Courant-Fischer theorem, The American Mathematical Monthly, 94, 352-354 (1987) · Zbl 0623.15010
[35]	Nedelec, J., Acoustic and Electromagnetic Equations, Integral Representations for Harmonic Problems (2001), Springer: Springer New York · Zbl 0981.35002
[36]	Gahalaut, K.; Kraus, J.; Tomar, S., Multigrid methods for isogeometric discretization, Computer Methods in Applied Mechanics and Engineering, 253, 413-425 (2012) · Zbl 1297.65153