Questions tagged [navier-stokes]

Questions about solution methods of the Navier-Stokes equations, related physical constants and non-dimensional number. Also special methods to solve the equations including the assumptions and their implementation in order to simplify them. Also, questions regarding modelling of the non-linear term, coefficients of these model can be subjective of this title.

Filter by
Sorted by
Tagged with
0 votes
0 answers
72 views

The mathematical meaning of a zero gradient pressure boundary condition in the Navier-Stokes equations

I would like to solve the Navier-Stokes equations for the unsteady problem of the flow around a circular cylinder. I would like to understand how to write mathematically the boundary condition for the ...
user avatar
1 vote
0 answers
41 views

Manufactured solution for incompressible navier-stokes

I am using the following manufactured solution for testing a incompressible navier-stokes finite volume solver solver $$ u = 0.25 \sin{\left(\pi x \right)} \cos{\left(\pi y \right)} + 0.25 \\ v = - 0....
user avatar
  • 21
1 vote
0 answers
56 views

How to assign initial velocity field and handle pressure-velocity coupling in FVM?

I am trying to solve the 2D incompressible Navier-Stokes equations for laminar flow over a backward facing step using the finite volume method. This is the plot that I generated of a generic mesh ...
user avatar
0 votes
0 answers
66 views

Slow convergence of Stokes solver used with the Immersed Boundary method

I am using Immersed Boundary Method to simulate elastic particles in 3D Stokes flow. Specifically, one has $\nabla ^2 \mathbf{u}-\nabla p + \mathbf{f}(t) = 0$, $\nabla \cdot \mathbf{u} \; $, where $\...
user avatar
3 votes
1 answer
116 views

Discontinuous pressure elements for incompressible Navier-Stokes

I am looking for some LBB-stable velocity-pressure combinations for incompressible Navier-Stokes where the pressure space is element-wise discontinuous, preferably with a linear variation elementwise. ...
user avatar
  • 565
1 vote
0 answers
60 views

Good non oscilliatory derivatives for an exsisting grid

I'm calculating the entropy production of a shockwave by utilizing the equations: \begin{equation} \sigma = J'_q\frac{\partial}{\partial x}\left(\frac{1}{T}\right) +\frac{1}{T}\frac{4\eta}{3}\left(\...
user avatar
  • 55
3 votes
1 answer
97 views

Discretizing the viscous component in 1 - D Navier stokes compressive flow

I've been working on modelling the NS equations in order to simulate shock waves. The equations are set up on the form: \begin{equation} \frac{\partial U}{\partial t} + \frac{\partial F(U)}{\partial x}...
user avatar
  • 55
5 votes
1 answer
184 views

Write incompressible Navier Stokes as ODE in $(\mathbf{u},p)$

Consider the Navier stokes equation after the discretization with conforming finite elements with source term $f=0$. We have the algebraic structure of a saddle point problem: $$M \dot{u} = f- Au -B^...
user avatar
1 vote
1 answer
110 views

Locking phenomena for $P1 - P0$ elements

Consider the Stokes problem and the usual divergence operator $B:V \rightarrow Q'$, $\langle Bv, q\rangle = b(v,q)=(\operatorname{div} v,q)$ and its discrete versione $B_h : V_h \rightarrow Q_h'$. In ...
user avatar
5 votes
0 answers
147 views

About the condition $\ker(B_h) \subset \ker(B)$ in mixed finite elements formulation

I'm studying mixed finite elements. The problem is a classical saddle-point one: we seek for $(u,p)$ in $V \times Q$: $$A u + B^t p = f$$ $$Bu = g$$ where $A: V \rightarrow V', B:V \rightarrow Q'$ ...
user avatar
2 votes
0 answers
94 views

Understanding inf-sup conditions for classical saddle point problems

I'm studying the inf-sup conditions for saddle point problems. I'm referring to the usual one $$\begin{cases}Au + B^t p = f \\Bu=g \end{cases}$$ In the book I'm using (Ern - Guermond: Theory and ...
user avatar
4 votes
1 answer
194 views

Discrete divergence free functions

I'm studying the weak formulation of NS equations. During the analysis, the book I'm using (Quarteroni-Valli, page 301-302), defined $$Z_h=\{v_h \in V_h: (\operatorname{div}(v_h),q_h)=0 \quad \forall ...
user avatar
3 votes
1 answer
200 views

Time discretization Navier Stokes equation

This question is a follow-up of this one. The weak form of Navier Stokes equation is (assuming $v,q$ test functions for the velocity and the pressure, respectively) $$(\frac{du}{dt},v)_{\Omega} + (\...
user avatar
  • 289
2 votes
0 answers
101 views

Confusion about preconditioner for incompressible Navier-Stokes equation with implicit-explicit method

Consider the time-dependent Navier-Stokes equation $$u_t + (u \cdot \nabla) u - \Delta u + \nabla p = f$$ $$\operatorname{div}(u)=0$$ Looking at deal.ii tutorials, I've notice that there are ...
user avatar
  • 289
1 vote
1 answer
108 views

What is the rationale of second-order finite volume discretization?

When it comes to a second-order accurate finite volume discretization of Navier-Stokes equations, which one of the two following rationales is adopted? 1- Second-order accuracy is a direct consequence ...
user avatar
  • 235
1 vote
0 answers
56 views

Why does including the pressure in this FVM for Stokes 2nd Problem lead to wrong solutions?

I'm trying to learn how to use finite volume methods and I want to solve a more general case of Stokes' second problem i.e. an infinite half-plane oscillating harmonically with no-slip boundary ...
user avatar
  • 136
7 votes
1 answer
201 views

Projection method FVM poisson part, adding source term

The idea of the method is to decompose the Navier-Stokes equation into the solenoidal and irrotational parts. $$\frac{\partial u}{\partial t}+u(\nabla \cdot u)=-\frac{1}{\rho}\nabla p+\nabla ^2 u$$ ...
user avatar
  • 316
3 votes
2 answers
127 views

Efficient schemes for solving the extended Saddle point problem

I am interested in knowing some efficient techniques for solving the following extended Saddle point problem. \begin{align} \begin{bmatrix} A & B^T & C^T \\ B & 0 & 0 \\ C & ...
user avatar
  • 565
1 vote
0 answers
123 views

Lumped mass matrices for higher-order finite elements for CFD

Given that some of the mass lumping techniques, for example, row-sum lumping does not produce practically viable lumped mass matrices for all the element shapes, what are the techniques used for mass ...
user avatar
  • 565
2 votes
1 answer
79 views

Fix for FD WENO method for multi-component compressible flows

I'm solving two-dimensional four-component compressible Navier-Stokes equations with finite-difference WENO approach. The equations are pretty standard: $$ \frac{\partial U}{\partial t} + \frac{\...
user avatar
  • 285
1 vote
2 answers
332 views

How to apply central difference to viscous fluxes in 2D Navier-Stokes equations?

I'm trying to solve 2D unsteady compressible Navier-Stokes equations with finite-difference or finite-volume method. Here is the system, it's pretty standard: $$ \frac{\partial U}{\partial t} + \frac{...
user avatar
  • 285
2 votes
2 answers
114 views

Different form of the Navier--Stokes equations

Normally I write the incompressible Newtonian isothermal flow Navier--Stokes equations as follows: $$\displaystyle \frac{\partial v}{\partial t} -\nu\Delta v +\color{red}{(\nabla v)v} +...
user avatar
  • 493
3 votes
0 answers
86 views

Explicit DG time step restriction for compressible Navier-Stokes equations

Hesthaven's book 1 mentions the following time step restriction for Navier-Stokes equations (see (7.32) in 2008 edition) $$ \Delta t \approx \frac{h}{N^2} \frac{C}{|u| + |c| + \frac {N^2 \mu}{h}} $$ ($...
user avatar
2 votes
0 answers
65 views

Haw to apply central difference to viscous flux in energy equation?

In many modern papers Navier-Stokes equations are solved with finite-difference or finite-volume methods using WENO reconstruction for non-viscous fluxes and central differences for viscous ones. It ...
user avatar
  • 285
2 votes
2 answers
144 views

Simplest way to "upgrade" from Euler equations to Navier-Stokes equations in FV or FD framework

I have quite a lot of experience solving unsteady Euler equations, including multi-component ones, with in house-coded finite-difference and finite-volume methods, including MacCormack and MUSCLE ...
user avatar
  • 285
2 votes
0 answers
54 views

Solving Stokes Equations in 3D - Do I need to treat pressure-velocity coupling iteratively?

I need to develop a code to solve Stokes Equations in 3D in cubic geometries (structured grid, uniform mesh spacing). My code needs to take a pressure gradient in one direction as a BC (pinlet=p1, ...
user avatar
3 votes
0 answers
99 views

How to construct a Fortin Operator for Crouzeix-Raviart Element?

I want to prove the LBB condition for the Stokes Equations discretised by the Crouzeix-Raviart element. The continuous Stokes Equation in the weak formulation is Find $u \in H_0^1(\Omega, \mathbb{R}^...
user avatar
  • 329
1 vote
0 answers
97 views

Stokes problem with imposed acceleration on boundaries (projection scheme)

I am trying to solve FSI problems with finite elements and using a projection scheme (I am taking as reference the review of Guermond: Guermond, J. L.; Minev, P.; Shen, Jie, An overview of ...
user avatar
  • 11
0 votes
1 answer
137 views

FDM on nonlinear PDEs

I'm working with a 2D Navier Stokes PDE in the unstabilized version - the equation is a linear equation of the type $\frac{∂u}{∂t} = F(u,t)$. In order to perform time discretization with FDM (finite ...
user avatar
5 votes
1 answer
283 views

What kind of a researcher am I?

So far, I've worked a bit in modeling, simulations and simple lab experiments, and I've really enjoyed all three research methods to approach a single research question. I can write tricky (in terms ...
user avatar
0 votes
0 answers
74 views

Numerical flow visualization in 2D for a moving boundary,

I have a rigid body that moves according to a set of governing ODEs, and I'd like to numerically visualize the vortices that are shed by this object. How could I proceed? I've been reading up on the ...
user avatar
0 votes
0 answers
68 views

The relation between PDE order and discretization order

In Jasak's Ph.D. thesis (2000), a notion is given about discretization of a transport equation: For good accuracy, it is necessary for the order of the discretization to be equal to or higher than the ...
user avatar
  • 235
0 votes
1 answer
140 views

How to visualize the vorticity / flow for a rigid body moving through a fluid?

How can I write down two-dimensional Navier-Stokes equations for a simple rigid object immersed in a flow and freely falling due to gravity? I'm trying to view the vorticity that's induced by the ...
user avatar
0 votes
0 answers
82 views

Pressure boundary conditions in Stokes Equation in 2D

I am solving the steady-state incompressible Stokes equations in 2D: \begin{equation} \frac{\partial u_x}{\partial x} + \frac{\partial u_y}{\partial y} = 0, \end{equation} \begin{equation} \mu\left[\...
user avatar
2 votes
1 answer
184 views

What is the correct way to calculate deviatoric stress tensor in lattice Boltzmann method?

Due to my previous question, where I asked about flux calculation in lattice Boltzmann (LB) method here, I have more or less same question for deviatoric stress tensor calculation due to pseudo-...
user avatar
3 votes
3 answers
2k views

Lattice Boltzmann methods vs Navier stokes/ other eulerian methods for *water* simulation

Note, there is already a question here, however the answers don't answer the original question, let alone specific considerations when dealing with nearly in-compressible fluids (water). Another ...
user avatar
  • 151
1 vote
0 answers
24 views

Equi-order in pressure correction schme of Navier-Stokes equation

I am wondering if there is an stabilized equi-order scheme in pressure correction scheme in solving Navier-Stokes equation? Usually P2-P1 element combination is used to solve NS equation, and a PSPG ...
user avatar
2 votes
0 answers
103 views

Governing equations vs Transport equations

This is a basic question. But I did not find any explanations for this. How are governing equations, like mass, momentum, energy conservations equations, different from 'Transport equation'?. Is a ...
user avatar
  • 205
0 votes
0 answers
72 views

Stokes Equation fails to converge for an ellipse

This might be because of the mesh, but the following code blows up for all values of b not 1. Does anybody have any experience working with the ellipse mesh in Fenics? ...
user avatar
2 votes
2 answers
219 views

Automatic timestep adjustment in a CFD solver

I have developed my own 3D Finite Volume Navier-Stokes solver based on projection method for nonuniform grid. I am looking to incorporate automatic timestep adjustment at each time step based on ...
user avatar
1 vote
1 answer
95 views

Is "Gradient Computation" in Finite Volume Discretization Really 2nd order accurate?

Based on this, pp 245, we go through these steps to discretize a gradient statement, namely $\nabla\phi$: 1- Gauss theorem reads, $$ \int_V\nabla \phi dV = \oint_{\partial V}\phi dS $$ 2- Integral ...
user avatar
  • 235
0 votes
0 answers
150 views

Double mach reflection at a inclined wedge

I am running into a strange problem when solving the 2D compressible Euler equations on a inclined wedge. To elaborate, my top boundary condition seems to emitting some type of instability. I have ...
user avatar
  • 75
0 votes
0 answers
60 views

3D Tollmien-Schlichting Waves Imposed in a Channel Flow (Are Physics correct?, etc)

So I am trying to do some further tests on a 2nd-order code Incompressible Navier Stokes equations, by studying transition to turbulence in a Poiseuille flow. Specifically, I'm interested to see ...
user avatar
  • 261
0 votes
2 answers
155 views

Taylor-Hood finite hexahedral elements, pressure diverging

I am developing a FEM fluid solver using the Taylor-Hood algorithm, i.e. quadratic interpolation for velocity, and linear for pressure. I have developed the code for 2-D quadrilaterals and triangles, ...
user avatar
1 vote
0 answers
62 views

Combining fluid flow solver based on lattice Boltzmann method with a mechanical deformation solver based on finite element method

I'm thinking to couple my fluid flow solver based on lattice Boltzmann method with a mechanical deformation solver based on finite element method to take account for solid deformation in my models. In ...
user avatar
  • 173
1 vote
0 answers
427 views

Inflow and outflow boundary conditions for advection-diffusion equation

I'm trying to solve this advection-diffusion equation (ADE): $$\frac{\partial \phi}{\partial t} + \nabla \cdot (-D \nabla \phi + \mathbf{u} \phi) = 0$$ In fact, this ADE framework is coupled to a ...
user avatar
  • 173
4 votes
1 answer
406 views

How to go from turbulent RANS to laminar Navier-Stokes and Euler

SU2 is an open-source CFD suite that is built around a RANS-solver. The main PDE that is solved, is the following: $$ \frac{\partial}{\partial t} \mathbf{U} + \nabla \cdot \mathbf{F^c} - \nabla \cdot ...
user avatar
  • 141
5 votes
0 answers
140 views

How can Navier--Stokes equations have asymmetric solutions such as Karman vortex streets

The Navier--Stokes equations are axially symmetric, so with symmetric boundary conditions, how can features such as Karman vortex streets develop? I understand that in reality symmetry does never ...
user avatar
  • 717
2 votes
0 answers
122 views

Derive total energy balance equation from Chapman-Enskog analysis of lattice Boltzmann equation

I'm interested to derive the total energy balance from Chapman-Enskog analysis of lattice Boltzmann equation (LBE). I know, I should go to the second moment of LBE (zeroth moment gives mass ...
user avatar
  • 173
3 votes
0 answers
79 views

Solving numerically a linearized system of elliptic (?) Navier-Stokes equation (Shallow Water Derived)

For my PhD Thesis, my advisor asked me to build a solver inspired from the article "Optimal Control Theory Applied to an Objective Analysis of a Tidal Current Mapping by HF Radar, J-L Devenon, 1989". ...
user avatar