Questions tagged [fluid-dynamics]

The study of the properties of fluids and gases in motion

Filter by
Sorted by
Tagged with
1
vote
0answers
30 views

Pseudospectral method for Rayleigh-Benard with constant temperature gradient

$$ \nabla\cdot \mathbf{u} = 0 \\ \frac{\partial \mathbf{u}}{\partial t}+\left(\mathbf{u}\cdot \nabla\right)\mathbf{u} = -\nabla p+\nu\nabla^2\mathbf{u}+\alpha g\theta\mathbf{e}_z\\ \frac{\partial\...
2
votes
0answers
37 views

Verification on pressure predictor method for CFD code

I have developed a python code for a lid-drive cavity model. However, my results are not converging. The algorithm of my code looks like this: Euler Momentum Equation looks like this: $$\frac{u^{n+1}...
1
vote
0answers
16 views

What is the exponential trick to include laplacian term in Rayleigh-Bernard simulation

I have come across a Rayleigh-Bernard simulation code which doesn't have the laplacian term but an integrating factor (in the exponential form) containing viscosity and diffusivity. I found out that ...
2
votes
0answers
52 views

Solving a simple Shallow Water model

Hi. I have a question at Mathematics and they suggested post here, once it's not common. I transcript as following. Many thanks I need to solve with basic methods this simple Shallow Water Model: $...
0
votes
0answers
30 views

Is semi-Lagrangian 1D advection identical to upwind Euler, when $|u|<\Delta x/\Delta t$?

This looked true to me, I worked through the algebra for 1D and confirmed it. But no one seems to mention it, so maybe I'm missing something... I'm using Bridson's SIGGRAPH 2007 course notes. [5MB ...
6
votes
0answers
116 views

What are some good debugging habits for numerical simulation?

I'm currently writing a lid drive cavity CFD code on python. Currently, my code has some issues (values jumping bear b.c). I was wondering what are some good habits in debugging numerical codes. ...
0
votes
0answers
8 views

Universal formulation of adiabatic equations of state in compresible finite-volume simultions

I code some finite element solver which should work for broad variety of materials (i.e. gas, liquid, solid, plasma) and large span of compressions resp. densities. I want to simulate things like ...
1
vote
1answer
59 views

What will be the impact of quantum computing on existing numerical techniques (e.g. CFD)?

Quantum computing seems to be a very active and promising development area in computer science. However, I am curious as to what impact (if any) quantum computing will have on existing classical ...
1
vote
0answers
44 views

WENO scheme on curvilinear coordinates

I've been developing a curvilinear FVM code. So far I've implemented the PPM scheme and am looking into adding WENO schemes. So far I've been discretizing the grid metrics using a second-order central....
0
votes
0answers
37 views

OpenFOAM and CFDEM on GPU

I have a simulation project written for OpenFOAM and CFDEM and would like to find an alternative to run it on GPU since raising the number of cores already provided a promising speed up and ...
1
vote
0answers
27 views

turbulant flow modeling [closed]

I have been trying to model the flow in a very thin channel using COMSOL Multyphysics. The laminar flow module does not converge no matter how I change the meshes, or use ramping loads and such. The k-...
1
vote
1answer
56 views

How to get started with numerically solving a Stochastic Navier Stokes equation

I originally posted the question on the math stackexchange, and was told I should try here. I’m researching Stochastic PDE, in particular the Navier Stokes Equation, and would like to estimate the ...
1
vote
0answers
31 views

Dealing with spurious oscillations in particle tracking methods

I work on modelling high intensity discharge xenon-filled lamps. The model governing the discharge is quite complex and sadly includes fluid dynamics. After some time, I managed to implement a finite-...
0
votes
0answers
10 views

Compute cell values and isosurface constant from density values of particles

I am trying to reconstruct the surface for a fluid simulation based on a list of particles using the Marching Cubes algorithm. From different resources, such as http://paulbourke.net/geometry/...
2
votes
2answers
143 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 ...
0
votes
1answer
63 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 ...
0
votes
0answers
74 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 ...
-1
votes
1answer
55 views

Simulating magnetic particles in a field free point generated by two opposing magnets

This is probably a long shot with such a short time, but I've been trying to get theoretical data for a project I'm working on. The project involves using a very simplified version of magnetic ...
1
vote
0answers
25 views

Instability in Lattice Boltzmann Solver

I wrote a Lattice Boltzmann Solver in Rust a little while back using both the BGK approximation and the TRT (two relaxation time) method on the D2Q9 lattice. In both cases I run into major stability ...
0
votes
1answer
72 views

Imposition of Dirichlet BC for Fourier pseudospectral in this paper

I was trying to implement the algorithm from the paper "Adapting a Fourier pseudospectral method to Dirichlet boundary conditions for Rayleigh–Benard convection". I am having a hard time to ...
0
votes
0answers
29 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 ...
0
votes
1answer
116 views

Finite Element - Flux Calculation

I am solving an advection-diffusion equation using the FEM and am having trouble calculating my fluxes. I start with the equation, $$\frac{\partial n}{\partial t} = \frac{\partial j_{n}}{\partial x}\...
6
votes
3answers
245 views

Why is the FVM traditionally used in CFD, and FEM in computational structures?

Most CFD codes use FVM. Most computational structures codes use FEM. Why is the FEM not frequently used in CFD, and why is FVM not frequently used in FEM?
3
votes
2answers
111 views

Finite volume discretization of non-conservative linear hyperbolic equation

Problem. Consider the one-dimensional adjoint Euler equations for $(x,t) \in \Omega \times [0,T]$ with $\Omega \subset \mathbb{R}$ and $T > 0$ $$ \varphi_t + \Big(\frac{\mathrm{d}F}{\mathrm{d} U}(x)...
1
vote
1answer
144 views

Normalization of polynomials for discontinuous Galerkin methods (DGM)

I was curious if someone could share their opinion on this matter. I have noticed that some people in literature normalize their Legendre polynomials, i.e. divide or multiply the polynomial by $$\...
1
vote
1answer
115 views

Can OpenFoam be used for interactive simulations (flying a plane)?

I'm a software engineer and I've just started to fly rc-planes. I'm currently building a self learning auto-pilot for my plane and I'm really interested in CFD. I installed OpenFOAM and ran through ...
0
votes
1answer
160 views

Heat diffusion - Is this the correct approach to include Newmann boundary conditions?

Thank you for looking at this problem. Is this the correct approach to include neumann boundary conditions? With this solution temperature is not correct, and there´s no diffusion. The model seems ...
1
vote
0answers
92 views

How to simulate water, falling under gravity, and impinging on a curved surface, which is kept/present in a domain, containing air?

TL;DR: How do I simulate a hole, at the bottom of a (full) water tank? I am attempting to simulate water, flowing out of a hole/slit, at the bottom of a tank (Water Domain) (under the influence of ...
0
votes
1answer
1k 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 ...
1
vote
2answers
66 views

TVD for temporal dicretisation

I have come across schemes where TVD (with flux limiters) is used for spatial discretisation along with Runge-kutta for Temporal discretisation. Can TVD be used for Temporal discretisation? If so ...
0
votes
1answer
66 views

What's a time centered Riemann problem?

I am trying to understand the meshless methods as described in https://arxiv.org/pdf/1409.7395.pdf. I'm having trouble understanding the following step: (Page 7, just after equation 17) Now, rather ...
2
votes
0answers
56 views

Why do stabilized formulations for the Navier-Stokes equation maintain the convergence rate for high order polynomial interpolation?

I have a quick questions which has been troubling me lately. When reading the FENICS Finite Element Book they assess various approaches to solver the Stokes equation. Obviously, they discuss the ...
3
votes
2answers
149 views

Asymptotic error of forward Euler

I'm trying to understand the asymptotic error behaviour of forward Euler (finite difference method), as timesteps are decreased (refined), so I feel trust in the method of manufactured solutions (MMS)....
0
votes
1answer
40 views

Modeling Diodes in Autodesk CFD

I'm extremely new to Autodesk CFD, and I'm working on a project that deals with diodes heating up and cooling down based on a fixed temperature regulated by a temperature switch. Basically, I have a ...
2
votes
1answer
89 views

Numerically solving a partial differential equation

I am trying to numerically solve the following PDE, $$\frac{\partial u^A}{\partial t} = c_1\frac{\partial^2 u^A}{\partial^2x} \,,$$ where $c_1$ is a constant. The above can be discretized using the ...
1
vote
2answers
73 views

Guaranteed equality between binary results with increasing MPI processes

Testing on an MPI scientific code for compressible flow dynamics I noticed that the results may depend on the number of processors used for the calculation. In fact, comparing the binary files they ...
1
vote
0answers
51 views

Has there been a comparison bewteen SIMPLE/SIMPLER and JFNK for steady CFD?

I'm looking for a comparison between the Jacobian-Free Newton-Krylov (JFNK) method performance compared to the conventional CFD nonlinear solution methodologies like SIMPLE. Does anyone know if such ...
2
votes
0answers
55 views

WENO methods: why the characteristic wise method resulting big errors?

I was doing my research/project using WENO as the limiter in finite volume methods to solve hyperbolic conservation law. I have no idea why the result in the characteristic wise method has a big error ...
3
votes
1answer
57 views

Plotting or Visualizing a Higher dimensional vector field

I am trying to visualize a higher dimensional vector field. Is there a way to do this. I asked this question here, I was told to post it here. As an example, one can use $$\begin{eqnarray} \dot{x}...
4
votes
1answer
98 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 ...
1
vote
0answers
38 views

How to remove sigularities from 3D Vortex Lattice Method

I'm trying to solve aerodynamics of whole aircraft by vortex lattice method (or deeper here). The problem is that sometimes trailing vortex filaments from Horse Shoe Vortexes of main wing hits panels ...
3
votes
0answers
56 views

Should I expect computational gains using a second-order splitting method here?

I am trying to solve a three-dimensional baroclinic transport problem. The hydrodynamic (three-dimensional shallow water) equations are: \begin{align} \nabla \cdot \vec{v} = 0, \tag{1} \\ \frac{\...
2
votes
1answer
84 views

When writing a research article, how to deal with deviations from theoretical expectations due to numerical errors?

We research the aerodynamics of the vehicle. Usually, there exists turbulence, which will cause the time-varying measurement of the force. The unsteady flow field determines we could not get a ...
1
vote
2answers
434 views

Demo example for OpenFOAM with CUDA

I am looking for a simple usage example/demo of OpenFOAM + CUDA and would like to understand how exactly OpenFOAM benefits from CUDA. The thing is I do not have any background in fluid dynamics and so ...
1
vote
1answer
44 views

Loss of energy when using Roe Solver for solving onedimensional Shallow Water Equations

I have written a Roe solver with Harten entropy fix code in Matlab to numerically solve the one-dimensional Shallow Water Equations. : \begin{eqnarray} \dfrac{\partial h(x,t)}{\partial t} + \dfrac{\...
0
votes
1answer
166 views

Poincare map for Arnold-Beltrami-Childress Magnetic Field in Python

I want to plot the Poincare map for Arnold-Beltrami-Childress magnetic field for parameters $A=1, B=0.816, C=0.5773$ in Python for the Poincare section $z=0$. Also, I am not able to understand what ...
3
votes
1answer
114 views

Comparison between FEM and FDM methods for flow simulations

What are the main differences between finite element and finite difference approach for incompressible flow simulations? I have a vague idea about how FE methods rely on minimizing the residual over ...
4
votes
0answers
167 views

A good 2D finite difference for the continuity equation

How could I go about solving the continuity equation below in 2D? $$\frac{\partial \rho}{\partial t} + \nabla \cdot \left(\rho u\right)=0$$ I saw that a similar question was posted here: A good ...
2
votes
1answer
110 views

build a simple incompressible solver with spectral method

I am trying to build a simple solver for the incompressible fluid in a periodic box to learn the spectral method. I am following the textbook (Peyret R. Spectral methods for incompressible viscous ...
2
votes
2answers
140 views

Algorithm to generate water flow map, given a terrain

I've posted the same question at GameDev Stack Exchange, but unfortunately I am not getting any response. So I am going to post ( and reword) it here. Hopefully I can get an answer! I have a terrain (...