Hi. I am trying to solve two-dimensional Stokes-Darcy equations using Gridap.jl. One of the terms in a weak formulation has the form {\int_{I}\lambda\left(u\cdot\tau\right)\left(v\cdot\tau\right)ds} where \tau is a unit tangent vector to the interface I between two subdomains.

Is there a way to specify a unit tangent vector in Gridap.jl? Specifically, I am looking for a function similar to `get_normal_vector`

.

Canâ€™t you use \int \lambda (u \times n) \cdot (v \times n) \, ds instead, where n is a unit normal?

@stevengj Thank you for your reply. However I am solving a two-dimensional problem so I cannot use the cross product.

I believe that Gridap defines the cross product of 2d vectors as just a scalar corresponding to the â€śzâ€ť component, which is exactly what you want here (u \times n is then equivalent to (u \times n) \cdot e_3 = u \cdot (n \times e_3) = u \cdot \tau).

The cross product also gives you a way to get a tangent vector if you want it â€” given a 2d normal vector n = (n_1, n_2), a tangent vector is \tau = n \times e_3 = (n_2, -n_1). But it doesnâ€™t seem like you actually need to construct this explicitly.

4 Likes