--- title: Generalized Taylor–Hood FEM Pair url: https://www.emergentmind.com/topics/generalized-taylor-hood-finite-element-pair type: topic --- # Generalized Taylor–Hood FEM Pair to=arxiv_search.search ￣影音先锋json {"query":"all:\"generalized Taylor-Hood\" OR ti:\"Taylor-Hood\" AND (Stokes OR surface OR unfitted OR anisotropic)", "max_results": 10, "sort_by": "relevance", "sort_order": "descending"} to=arxiv_search.search 天天彩票中奖json {"query":"generalized Taylor-Hood finite element pair", "max_results": 10, "sort_by": "relevance", "sort_order": "descending"} to=search_arxiv 天天中彩票是不是json {"query":"generalized Taylor-Hood finite element pair", "max_results": 10, "sort_by": "relevance", "sort_order": "descending"} to=arxiv.search 天天中彩票无法json {"query":"Taylor-Hood Stokes generalized finite element", "max_results": 10} The generalized Taylor–Hood finite element pair denotes a family of mixed finite element constructions that preserve the defining Taylor–Hood degree offset—higher-order approximation for the primal vector field and one-order-lower approximation for the pressure or pressure-like Lagrange multiplier—while extending the classical fitted Stokes setting to unfitted interfaces, embedded surfaces, weighted axisymmetric formulations, anisotropic meshes, nonconforming or reduced spaces, and related multiphysics systems [2101.09627] [1206.4379] [2003.06972] [2401.03561]. Taken together, these works suggest that “generalized Taylor–Hood” is best understood as a design principle rather than a single element family: the polynomial hierarchy of Taylor–Hood is retained, but geometry representation, conformity, stabilization, norms, and even the ambient PDE class may be altered to preserve inf-sup stability, consistency, and optimal approximation. ## 1. Classical template and the meaning of “generalized” In its classical conforming form, the lowest-order Taylor–Hood pair consists of continuous piecewise quadratic velocities and continuous piecewise linear pressures. One formulation writes \[ V_h=\mathcal L_2^1(\mathcal T;\mathbb R^d)\cap W^{1,1}_0(\Omega;\mathbb R^d), \qquad Q_h=\mathcal L_1^1(\mathcal T)\cap L^2_0(\Omega), \] and establishes a constructive Fortin operator in any spatial dimension \(d\ge 2\) [2104.13953]. On tensor-product meshes, the corresponding family is the continuous \(\mathcal Q_k-\mathcal Q_{k-1}\) pair, \(k\ge 2\), on quadrilateral or hexahedral meshes [2205.15770]. In both cases, the essential structural feature is the degree shift between velocity and pressure. What changes in the generalized setting is not this degree pattern but the analytical and geometric context in which it is embedded. The phrase is used explicitly for the unfitted two-phase pair \(\mathbf P_{k+1}-P_k\), \(k\ge 1\), with phase-wise duplicate spaces, Nitsche coupling across an interface, viscosity-weighted pressure normalization, and ghost-penalty stabilization on cut patches [2101.09627]. It is also used for trace and surface formulations, where bulk Taylor–Hood spaces are restricted to an embedded surface or lifted to a high-order surface approximation, and for weighted axisymmetric Stokes problems, where the usual \(P^{k+1}/P^k\) structure is analyzed in weighted Sobolev spaces on graded meshes [1206.4379]. A closely related terminology appears in “Taylor–Hood-like,” “reduced Taylor–Hood-type,” and “surface analogue of Taylor–Hood” constructions [2205.02013] [2105.00144] [2506.20419]. A representative summary of the main variants is as follows. | Setting | Discrete pair | Additional mechanism | |---|---|---| | Unfitted two-phase Stokes | \(\mathbf P_{k+1}-P_k\) or \(\mathbf Q_{k+1}-Q_k\) | Phase-wise fields, Nitsche coupling, ghost penalty [2101.09627] | | Surface trace FEM | \(P_k\)-\(P_{k-1}\) traces from bulk spaces | Tangential penalty, normal-derivative stabilization [2003.06972] | | Axisymmetric weighted Stokes | \(P^{k+1}/P^k\) | Weighted spaces and graded meshes [1206.4379] | | Anisotropic meshes | \(\mathbb Q_2\times\mathbb Q_1\), \(\mathbb P_2\times\mathbb P_1\) | Generalized Verfürth trick on edge and corner patches [1710.07857] | | Penalty-free surface method | \(\mathbb P^r-\mathbb P^{r-1}\) on \(\Gamma_{h,k}\) | Exact tangentiality, Piola mapping, Gauss–Lobatto edge nodes [2506.20419] | This classification also clarifies a useful boundary. Some works employ the standard Taylor–Hood pair inside a larger stabilized or multiphysics algorithm without altering the element family itself. One explicit example states that no new generalized Taylor–Hood element is introduced; instead, the classical \(Q_2/Q_1\) pair is used within a stabilized Newton–Galerkin method for Darcy–Brinkman–Forchheimer flow [2501.04041]. ## 2. Unfitted and interface-coupled generalizations A central line of development concerns interface problems in which the physical interface cuts through a fixed background mesh. In the two-phase Stokes problem with slip between phases, the generalized Taylor–Hood pair is \[ V_h^\pm = \{\, \mathbf{v}_h=(\mathbf{v}_h^-,\mathbf{v}_h^+) \in V_h^- \times V_h^+ \,\}, \qquad Q_h^\pm = \left\{\, q_h=(q_h^-,q_h^+) \in Q_h^- \times Q_h^+ \;:\; \int_{\Omega^-}\mu_-^{-1}q_h^- + \int_{\Omega^+}\mu_+^{-1}q_h^+=0 \right\}, \] with \(V_h^\pm\) built from phase-wise \(\mathbf P_{k+1}\) velocities and \(Q_h^\pm\) from phase-wise \(P_k\) pressures, and with the corresponding \(\mathbf Q_{k+1}-Q_k\) variant admitted on quadrilateral or hexahedral meshes [2101.09627]. Because the method is formulated on the background mesh with separate phase fields, discrete components may be multivalued in the overlap region \(\Omega_h^-\cap\Omega_h^+\). The interface conditions are imposed weakly through Nitsche terms, while ghost-penalty stabilization controls cut-patch pathologies. The principal theoretical result is a robust inf-sup theorem. Under the assumptions that the mesh is sufficiently fine, the cut region is quasi-uniform, exact integration is used on cut cells, and the averages are chosen as \(\alpha=0\), \(\beta=1\) with \(\mu_-\le \mu_+\), the method satisfies \[ \sup_{(\mathbf{v}_h,q_h)\in V_h^\pm\times Q_h^\pm} \frac{(\mathbf{u}_h,p_h;\mathbf{v}_h,q_h)}{\|\mathbf{v}_h,q_h\|} \ge C\,\|\mathbf{u}_h,p_h\|, \] with \(C\) independent of the viscosity ratio, slip coefficient, interface position relative to the background mesh, and mesh size \(h\) [2101.09627]. Coercivity of the velocity form is proved with the same parameter robustness, and the resulting error estimate is \[ \| \mathbf{u}-\mathbf{u}_h,\; p-p_h\| \le C h^{k+1}\|\mathbf{u},p\|_*, \] again with \(C\) independent of \(h\), \(\mu_\pm\), \(f\), and the interface position. For \(\mathbf P_2-P_1\), the reported two- and three-dimensional experiments show approximately third order in velocity \(L^2\) and second order in weighted \(H^1\)-velocity and weighted \(L^2\)-pressure [2101.09627]. A closely related earlier unfitted construction for stationary Stokes interface problems starts from the classical \(P_2/P_1\) pair, then duplicates the spaces phase-wise: \[ \mathbf{V}_h^\Gamma := \mathbf{V}_h|_{\Omega_1} \oplus \mathbf{V}_h|_{\Omega_2},\qquad Q_h^\Gamma := Q_h|_{\Omega_1} \oplus Q_h|_{\Omega_2}, \] so that the discrete velocity can represent kinks and the pressure can represent jumps across the interface [1605.04085]. A Nitsche-type formulation weakly enforces continuity and traction balance, and a ghost penalty on pressure derivatives supplies an inf-sup stable variational formulation with a constant independent of the interface position relative to the mesh. The numerical comparison in that work isolates two distinct ingredients: velocity enrichment is needed to recover interface approximation quality, and high-order parametric geometry mapping is needed to recover full convergence order [1605.04085]. These unfitted results make clear that the generalized character of the pair is not merely geometric. The velocity-pressure spaces are altered by duplication, the pressure normalization is reweighted, and the underlying inf-sup proof is rebuilt around ghost penalties, cut-cell trace inequalities, and interface-consistent Nitsche couplings [2101.09627] [1605.04085]. ## 3. Surface and trace variants For the surface Stokes equations, generalized Taylor–Hood pairs arise because one retains the Taylor–Hood degree hierarchy while discarding the fitted planar mesh paradigm. In the trace FEM formulation on a smooth closed surface \(\Gamma\subset\mathbb R^3\), bulk spaces on a tetrahedral mesh are restricted to the surface: \[ \mathbf U_h=(V_h^{k+1})^3,\qquad Q_h=V_h^k\cap L_0^2(\Gamma),\qquad k\ge 1, \] with particular emphasis on the trace \(P_2\)–\(P_1\) pair [1909.02990]. Since the surface cuts arbitrarily through the background tetrahedra, pressure stabilization by a volume normal derivative term is essential. The main theorem for \(k=1\) states a discrete inf-sup estimate with a constant independent of both \(h\) and the position of the surface relative to the bulk mesh. The same work shows that without pressure stabilization the trace \(P_2\)–\(P_1\) pair loses inf-sup stability under refinement, whereas normal volume stabilization restores stability; full-gradient stabilization is also stable but produces larger consistency error and is suboptimal for the analyzed method [1909.02990]. Higher-order trace formulations generalize this construction to \[ \mathbf U_h := \mathbf V_{h,\Theta}^k, \qquad Q_h := V_{h,\Theta}^{k-1}\cap L_0^2(\Gamma_h), \qquad k\ge 2, \] using a parametric geometry map \(\Theta_h\), penalty enforcement of tangentiality, and normal-derivative volume stabilization for both velocity and pressure [2003.06972]. The central discrete inf-sup bound is \[ c_0\|q_h\|_M \le \sup_{v_h\in \mathbf U_h}\frac{|b_h(v_h,q_h)|}{\|v_h\|_A} + \tilde s_h(q_h,q_h)^{1/2}, \] with \(c_0\) independent of \(h\) and of how \(\Gamma_h\) cuts the background mesh. The resulting energy-norm estimate is optimal of order \(h^k\), while the quantified geometric consistency error is one order higher and therefore non-dominant [2003.06972]. A parallel numerical study of the same family, phrased as higher-order trace FEM with a consistent and an inconsistent penalty formulation, recommends the isoparametric choice \(k_g=k\), the stabilization scalings \(\rho_u\sim h^{-1}\), \(\rho_p\sim h\), and a more accurate penalty normal of order \(k_p=k+1\); for the consistent method \(\eta\sim h^{-2}\), while for the inconsistent method \(\eta\sim h^{-(k+1)}\) [1909.08327]. The fitted surface finite element analogue uses the classical surface FEM of Dziuk–Elliott combined with a Hood–Taylor pair on a high-order surface approximation \(\Gamma_h^k\). The spaces are \[ V_h=\{v_h\in C(\Gamma_h^k)^3:\ v_h\circ\pi_k^{-1}\in \tilde V_h\}, \qquad Q_h=\{q_h\in C(\Gamma_h^k):\ q_h\circ\pi_k^{-1}\in \tilde Q_h\}, \] with \(\tilde V_h=(V_h^m)^3\), \(\tilde Q_h=V_h^{m-1}\), \(m\ge 2\), and a penalty parameter chosen as \(\eta=h^{-2}\) [2401.03561]. The main error bound has the form \[ |u_T^e-u_h|_k+\|p^e-p_h\|_{L^2(\Gamma_h^k)} \lesssim h^r(\|u_T\|_{H^{r+1}}+\|p\|_{H^r}) + h^k(\|u_T\|_{H^1}+\|p\|_{H^1}+\|f\|_{L^2}+\|g\|_{L^2}), \] so that for sufficiently accurate geometry, \(k\ge r\), the optimal order \(h^r\) is obtained in the energy norm [2401.03561]. A more recent surface construction removes penalization entirely. It defines a tangential, \(H(\mathrm{div})\)-conforming but \(H^1\)-nonconforming velocity space on \(\Gamma_{h,k}\), keeps the pressure in the standard continuous \(\mathbb P^{r-1}\) space, and proves the inf-sup condition \[ \sup_{v\in V_h}\frac{\int_{\Gamma_{h,k}} (\operatorname{div}_{\Gamma_{h,k}}v)\,q}{\|v\|_{H^1_h(\Gamma_{h,k})}} \ge \beta \|q\|_{L_2(\Gamma_{h,k})} \qquad \forall q\in Q_h. \] Its optimality depends on locating edge degrees of freedom at Gauss–Lobatto nodes; standard equispaced Lagrange edge nodes lead to a loss of one order for \(r\ge 3\) [2506.20419]. This result sharply illustrates a recurring theme in generalized Taylor–Hood theory: once conformity or geometry is altered, seemingly secondary implementation choices may become mathematically decisive. ## 4. Weighted, singular, and anisotropic adaptations The axisymmetric Stokes equations in polygonal meridian domains furnish a distinct generalization in which the polynomial pair remains classical, but the functional framework is not. After reduction from the three-dimensional axisymmetric domain to a meridian polygon \(\Omega\subset\mathbb R^2\) in \((r,z)\)-coordinates, the problem involves singular coefficients \(r^{-1}\partial_r\), \(r^{-2}u_r\), and the weighted divergence \((\partial_r+r^{-1})u_r+\partial_z u_z\) [1206.4379]. The natural weak formulation is posed in \[ (u_r,u_z,p)\in H^1_{-,0}(\Omega)\times H^1_{+,0}(\Omega)\times L^2_{1,0}(\Omega), \] and the regularity theory is developed in weighted Kondrat’ev-type spaces \(\mathcal K^m_{\mu,1}\). The discrete pair is still the Taylor–Hood pair, \[ \mathbf V_n^{k+1}\times S_n^k, \] with \[ \mathbf V_n^{k+1} = \left\{ (v_r,v_z)\in [P^{k+1}(\Omega)]^2: v_r|_{\Gamma\cup\Gamma_0}=0,\ v_z|_{\Gamma}=0 \right\}, \qquad S_n^k=\left\{p\in P^k(\Omega): \int_\Omega p\,r\,dr\,dz=0\right\}, \] but its analysis depends on weighted norms, weighted interpolation operators, and \(\kappa\)-graded meshes near singular vertices [1206.4379]. The grading parameter is chosen as \[ \kappa=\min\left(\frac12,\, 2^{-(k+1)/a}\right), \qquad 0