# Linsolve vs. PETSc for transient laminar Navier - Stokes

**URL:** https://discourse.julialang.org/t/linsolve-vs-petsc-for-transient-laminar-navier-stokes/116876
**Category:** Specific Domains
**Tags:** question
**Created:** [July 10, 2024, 3:48pm UTC](https://discourse.julialang.org/t/linsolve-vs-petsc-for-transient-laminar-navier-stokes/116876 "2024-07-10T15:48:22Z")
**Posts on this page:** 2
**Page:** 1

<div class="post-metadata">

### Author: ![Domenico\_Lahaye](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/domenico_lahaye/32/203728_2.png) [@Domenico\_Lahaye](https://discourse.julialang.org/u/Domenico_Lahaye)
#### Post date: [July 10, 2024, 3:48pm UTC](https://discourse.julialang.org/t/linsolve-vs-petsc-for-transient-laminar-navier-stokes/116876/1 "2024-07-10T15:48:22Z")

</div>

The question that follows was phrased at JuliaCon:

given: suppose solving a 2D / 3D transient laminar Navier-Stokes problems (lid-driven cavity, flow part cylinder or airfoil) using Gridap, Trixi or similar

how do the following 2 options compare: a/ using linsolve within DiffEQ or b/ using TS/SNES/KPS/PC within PETSc ?

does is make sense to employ the DAE formulation to treat the pressure equation?

---

<div class="post-metadata">

### Author: ![DanDoe](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/dandoe/32/52717_2.png) [@DanDoe](https://discourse.julialang.org/u/DanDoe)
#### Post date: [July 12, 2024, 9:47am UTC](https://discourse.julialang.org/t/linsolve-vs-petsc-for-transient-laminar-navier-stokes/116876/2 "2024-07-12T09:47:18Z")

</div>

So the thing is that Trixi.jl is fundamentally designed as a hyperbolic PDE solver, so the main scope are always phenomena with finite speed of information propagation.

That having said, stuff like Stokes Flow is really out of scope. For such scenarios I would forward you to [Waterlily.jl](https://github.com/WaterLily-jl/WaterLily.jl), which tackle explicitly incompressible flows, thus probably more suited.

Back to the question:

In Trixi.jl, there is no need to solve a linear system due to the nature of the DG solver we employ.
