Skip to main content
Share Your Experience: Take the 2024 Developer Survey

Questions tagged [fluid-dynamics]

The study of the properties of fluids and gases in motion

Filter by
Sorted by
Tagged with
1 vote
1 answer
14k views

Python Finite Difference Schemes for 1D Heat Equation: How to express for loop using numpy expression

Hello all, I've recently been introduced to Python and Numpy, and am still a beginner in applying it for numerical methods. I've been performing simple 1D diffusion computations. I suppose my ...
Leonardo Parra García's user avatar
22 votes
2 answers
5k views

A good finite difference for the continuity equation

What would be a good finite difference discretization for the following equation: $\frac{\partial \rho}{\partial t} + \nabla \cdot \left(\rho u\right)=0$? We can take the 1D case: $\frac{\partial \...
tiam's user avatar
  • 864
4 votes
1 answer
3k views

Comparison of Lattice Boltzmann Method vs Traditional Navier-Stokes based Methods

I have a choice of two options, analysing and implementing Lattice Boltzmann methods or traditional Navier Stokes based methods. I'm a CFD newbie and I have a rough idea (though not rigorous enough to ...
Aditya Kumar Praharaj's user avatar
33 votes
4 answers
2k views

Why is local conservation important when solving PDEs?

Engineers often insist on using locally conservative methods such as finite volume, conservative finite difference, or discontinuous Galerkin methods for solving PDEs. What can go wrong when using a ...
Jed Brown's user avatar
  • 25.7k
7 votes
3 answers
702 views

C - OpenMP, MPI, Serial Program

I'm part of a Computational Science course and come from a non-programming background, so please forgive me my ignorance. I'm working on a set of code in C to numerically solve the Navier Stokes ...
kfkhalili's user avatar
  • 319
7 votes
1 answer
5k views

computing turbulent energy spectrum from isotropic turbulence flow field in a box

I have my 3 dimensional velocity flow-field u, v and w at a given instant of time from DNS using pseudo-spectral method. I need to calculate the energy spectrum ( in Fourier space ) as a function of ...
user4557934's user avatar
6 votes
3 answers
5k views

How to find QR decomposition of a rectangular matrix in overdetermined linear system solution?

While trying to find cell-centered gradients in finite volume method computation of incompressible fluid flow I get over-determined linear system. This is a well known "cell based least-square" ...
Johntra Volta's user avatar
3 votes
1 answer
452 views

Finite Elements Weak Formulation generalization

I am struggling with an equation that represents the Weak form of Galerkin method: $ \phi^{T}F(\textbf{u})\sim \int_{\Omega}^{ } \phi.f_{0}(\mathit{u},\nabla \mathit{u}) + \nabla\phi:f_{1}(\mathit{u},...
Pedro R.'s user avatar
  • 157
2 votes
2 answers
508 views

In practice, what are the most useful ways to visualize 2d fluid flow, to tell what is happening in the simulation? Esp for verification and debugging

My simulation creates a 2d grid of vectors and scalars (EDIT of velocity, depth etc), at 60 frames per second. Is it correct? It looks sort of right... but who knows? How can I tell what's happening - ...
hyperpallium's user avatar
2 votes
2 answers
541 views

2d Euler manufactured solutions

Where can I find manufactured solutions for the 2d Euler equations, with the complete analytical terms, including the Jacobian of the source term ?
Gonzague's user avatar
2 votes
1 answer
451 views

Physical interpretation of divergence theorem

In a diverging pipe section like the following, the pipe of radius $r$ splits into two pipes of radius $r/2$. Consider a solute transported by convection from node 1. $$\frac{\partial C}{\partial t}...
Natasha's user avatar
  • 421
1 vote
1 answer
552 views

FFT Poisson Solver for non-uniform grid

I have a 3D solver for the incompressible Navier-Stokes equations which uses a FFT library for the Poisson equation with a uniform grid on all directions. In 2D the Poisson equation is given by: $$ ...
user26633's user avatar
24 votes
2 answers
6k views

Discontinuous Galerkin: Nodal vs Modal advantages and disadvantages

There are two general approaches to representing solutions in the discontinuous galerkin method: nodal and modal. Modal: Solutions are represented by sums of modal coefficients multiplied by a set of ...
Aurelius's user avatar
  • 2,375
17 votes
4 answers
3k views

Is there a good, easy-to-use, high quality open source CFD solver out there?

My thesis is on developing numerical methods for model reduction in combustion. I run my methods purely on the chemistry part of combustion simulations, and I have plenty of case studies for 0-D ...
Geoff Oxberry's user avatar
16 votes
2 answers
3k views

(how to) write simulations that run faster?

I have started using python as the programming language for doing all my assignments in CFD. I have a very little experience in programming. I am from mechanical engineering background and am pursuing ...
Subodh's user avatar
  • 1,460
16 votes
2 answers
4k views

Disadvantages of common discretization schemes for CFD simulations

The other day, my computational fluid dynamics instructor was absent and he sent in his PhD candidate to substitute for him. In the lecture he gave, he seemed to indicate several disadvantages ...
Paul's user avatar
  • 12k
14 votes
1 answer
460 views

What spatial discretizations work for incompressible flow with anisotropic boundary meshes?

High Reynolds number flows produce very thin boundary layers. If wall resolution is used in Large Eddy Simulation, the aspect ratio may be on the order of $10^6$. Many methods become unstable in this ...
Jed Brown's user avatar
  • 25.7k
12 votes
3 answers
2k views

Can compressible flow solvers be used to solve incompressible flow?

I know that incompressible and compressible flow solvers are specifically designed to solve different types of problems with different fluid properties/flow conditions. Clearly, among the advantages ...
Paul's user avatar
  • 12k
11 votes
4 answers
3k views

Constructing random divergence-free velocity fields

I am trying to simulate decaying homogeneous isotropic turbulence. As initial condition I want a divergence-free vector-field, i.e, $\mathrm{div} U = 0$. How do I initialize random velocity field in ...
verito's user avatar
  • 129
8 votes
7 answers
8k views

What is the difference between MATLAB and FORTRAN?

In our university some Ph.D students for computational methods prefer FORTRAN over MATLAB. I can't understand why? What is the difference between them when are used in computational methods like ...
MohammadSh's user avatar
8 votes
2 answers
364 views

Stabilization of convection-dominated flow and turbulence modeling

Are stabilization techniques for convection-dominated flows like SUPG+PSPG, interior penalty methods, etc. able to handle turbulent flows without tubulence model being employed, at least up to some ...
Jan Blechta's user avatar
8 votes
2 answers
6k views

simple MHD simulation code for (self) education and play with

I would like some super simple computational code for solving magnetohydrodynamics problems. High accuracy nor performance is not my concern. I wan't it just to visually explore qualitative behavior ...
Prokop Hapala's user avatar
7 votes
1 answer
418 views

Potential flow around a non-symmetric obstacle using stream functions

I've seen that there is a way to use the finite differences method, on a cartesian orthogonal grid, to perform calculations on potential flow about an obstacle without using the Neumann conditions, ...
john's user avatar
  • 71
6 votes
1 answer
250 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$$ ...
2Napasa's user avatar
  • 362
6 votes
1 answer
181 views

A simple PDE solution question

I need to ask a question about partial derivatives. I want to solve this equation (steady state, one dimensional continuity equation): $$\frac{\partial (\rho u)}{\partial z}=0$$ which is equivalent to:...
ASK22's user avatar
  • 71
6 votes
1 answer
944 views

ENO/WENO component-wise vs characteristic-wise

Can someone give some references to understand what's the differences between a component-wise and a characteristic-wise ENO scheme? If I'm right, the characteristic variables come from the ...
BRabbit27's user avatar
  • 1,029
6 votes
3 answers
2k views

Is lattice Boltzmann suitable for simulation of incompressible Stokes flow?

We have a flow that is dominated by adhesion forces from the substrate and surface tension from the free surface. The material is nearly solid and at rest first, and gets a bit less solid by heating. ...
Thomas Klimpel's user avatar
5 votes
3 answers
3k views

Flow past a cylinder - Projection Method - Boundary Conditions

I plan to write a code that solves the flow past a cylinder and try to see the Von-Karman vortex street. I will solve the 2-D viscous, incompressible Navier-Stokes using the projection method. The ...
l3win's user avatar
  • 205
5 votes
2 answers
333 views

Preconditioner for the GMRES method in the Uzawa algorithm

I'm trying to solve \begin{equation}\left\{ \begin{split} \frac{\partial u}{\partial t}+(u\cdot\nabla)u-\nu\Delta u+\frac1\rho\nabla p&=f\;\;\;\text{in }\Lambda\\ u&=0\;\;\;\text{on }\partial\...
0xbadf00d's user avatar
  • 283
5 votes
1 answer
4k views

Stability of numerical method for 1D Burger's equation

I am trying to solve 1D viscous Burger's equation numerically and I cannot apply von Neumann analysis because the equation is non-linear. How do I predict the stability criteria for my system? I also ...
tumchaaditya's user avatar
4 votes
1 answer
4k views

Finite difference methods in cylindrical and spherical co-ordinate systems

I am quite familiar with finite difference schemes in cartesian coordinates. The key point here is that every point in the cartesian grid is treated equally as the spacing between consecutive points ...
Babaji's user avatar
  • 195
4 votes
1 answer
1k views

Effecient CFD programming techniques

I'm trying to make highly efficient CFD programming complex for solving combustion problems. I've finished writing core which realises mathematical model, and now I'm concerned about code performance. ...
Yakovenko Ivan's user avatar
4 votes
4 answers
12k views

Estimating the Courant number for the Navier-Stokes Equations under differing Reynolds number regimes

I am familiar with the Courant-Friedrich-Lewy Condition in as far as it applies to the stability of explicit finite difference schemes for standard parabolic and hyperbolic PDEs. However, when ...
Paul's user avatar
  • 12k
3 votes
2 answers
551 views

Energy Conservation in Conservation Laws with Source Terms

I'm wondering if anyone can help me understand energy conservation when using conservation law methods (i.e. Riemann solver, High-Resolution Wave-Propagation Methods) with the addition of source terms:...
Wes Lowrie's user avatar
3 votes
2 answers
2k views

GPU-enabled Lattice Boltzmann solvers?

Is anybody aware of any GPU-enabled Lattice Boltzmann solvers (preferably on C++/OpenCL and open-source) that would be recommended? I have found Advanced Simulation Library, but it seems to be very ...
user1876484's user avatar
3 votes
3 answers
1k views

Solving Stokes flow with walls using Oseen tensor

Introduction I've developed a code to solve for generalised, incompressible 2D Stokes flow $\eta \nabla^2 \mathbf{v} - \nabla p + \mathbf{S} = 0$ $\nabla . \mathbf{v} = 0$ where $\mathbf{S}$ can ...
Hemmer's user avatar
  • 265
3 votes
0 answers
109 views

numerical instabilities in Fluid Dynamics, Finite Element Method

I'm looking for references to understand where the numerical instabilities come from in hydrodynamics in general, and notably when the Péclet number: $Pe>1$. I'm using the finite element method. ...
J.A's user avatar
  • 171
3 votes
1 answer
985 views

Add User-defined/custom differential equations in OpenFoam (CFD)

I am new to OpenFoam. And I am trying to add a set (user defined) of differential equations to OpenFoam. I want to solve this user defined set of equations at each time point in addition to standard ...
Ivan M's user avatar
  • 31
3 votes
0 answers
401 views

Euler Equation Eigensystem with Gravity in the Energy Flux

I am modifying a conservative form of the Euler equations with gravity in the energy flux (see previous question: Energy Conservation in Conservation Laws with Source Terms) for use in a Riemann ...
Wes Lowrie's user avatar
3 votes
1 answer
190 views

Numerical computation of the velocity in the steady Navier-Stokes equation

I've asked this question on Math.SE too. Let $d\in\left\{1,\ldots,4\right\}$ $\Lambda\subseteq\mathbb R^d$ be bounded, nonempty and open and $\partial\Lambda$ be Lipschitz $V:=\left\{u\in H_0^1(\...
0xbadf00d's user avatar
  • 283
2 votes
1 answer
167 views

Confusion about Zabusky and Kruskal's stepper for the KdV equation

In Zabusky and Kruskal's paper about solitons, they derive the following update for the Korteweg de Vries equation (their footnote 6): \begin{align*} u_{i}^{j+1} = u_{i}^{j-1} - \frac{1}{3} \frac{k}{...
user14717's user avatar
  • 2,155
2 votes
5 answers
4k views

Why are upwind schemes stable in convection flow calculation?

It is well known that upwind schemes are stable when calculating convection flows with $|\text{Pe}|>2$, $\text{Pe}$ is the Peclet number. Why is that, and why is central difference unstable? Is ...
Kozuki's user avatar
  • 99
1 vote
0 answers
85 views

PDE discretization (via finite difference sheme) question

So after posting this question and reading all your comments I would like to make this new question (update). If you consider the three equations presented here: $$\frac{\partial \rho}{\partial t} +\...
ASK22's user avatar
  • 71
1 vote
0 answers
92 views

Vortex Lattice Method: better basis than horse-shoe vortex?

I read several introductory texts about potential flow and vortex lattice method. Basically, it is fitting of some velocity field described by conditions on velocity at some control points using basis ...
Prokop Hapala's user avatar
1 vote
1 answer
347 views

How to compute turbulent energy cascade

I need to compute the kinetic energy cascade using a finite volume solution in an equally spaced grid. I wonder if it is more correct to first compute the kinetic energy in the space (or time) domain, ...
Adr's user avatar
  • 173
1 vote
2 answers
2k views

New to CFD, Lattice Boltzmann or Navier-Stokes?

I apologize if some of my questions are naive; I am very new to computer simulations and fluid-dynamics. I am going to start a PhD in early 2017, and I would like to bone-up on some Computational ...
user109527's user avatar
0 votes
1 answer
759 views

boundary conditions of linear advection problem

I am solving the 1D advection problem given by: $$\frac{\partial u}{\partial t}=-c\frac{\partial u}{\partial x}$$ where c is the wave speed, and u is the unknown field variable, and x and t are time ...
Hooman's user avatar
  • 169