Questions tagged [fluid-dynamics]

The study of the properties of fluids and gases in motion

Filter by
Sorted by
Tagged with
-3 votes
0 answers
32 views

Will we ever be able to control weather? [closed]

More specifically, will we be able to understand weather dynamics to such a computational detail that we will be able to influence the weather patterns? E.g. inducing rain during droughts over a ...
user avatar
0 votes
0 answers
29 views

How does Tannehill impose boundary conditions when coding the Parabolized Navier Stokes on an Implicit Finite Differences Scheme?

I'm trying to implement the scheme he describes on his book "Computational Fluid Mechanics and Heat Transfer" on Chap.9 and I'm having trouble imposing BC. I don’t get how he imposes them. I ...
user avatar
1 vote
0 answers
32 views

Adding wall friction to a 3D model - Stokes equation

Maybe that is not the best place for the following question, if such, let me know and sorry. I would like to model the flow of expanding bentonite in a fracture following the given reference. In order ...
user avatar
  • 89
0 votes
0 answers
75 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
195 views

CFD in Python: Mesh Generation

I am trying to set up a 3D CFD scheme for thermal and flow modelling in Python using the finite volume method. The first concern is to build the geometry and an accompanying mesh that is efficient for ...
user avatar
-1 votes
1 answer
58 views

The local and average Nusselt number in a square cavity

I am in the process of programming the local & average Nusselt number in a left vertical wall but my Matlab script gives me inappropriate values and it doesn't change with changing of Rayleigh ...
user avatar
1 vote
0 answers
36 views

Semi-lagrangian method for compressible fluid (non divergence free)

I am looking for a semi-Lagrangian method for advection with a non divergence free velocity field. The equation is \begin{align} \frac{\partial C(x,t)}{\partial t} &= - \nabla \cdot (\vec{v}(x,t) ...
user avatar
0 votes
0 answers
56 views

Passing boundary conditions to solver

Quite broad question, Currently building own Poisson solver subroutine for CFD solver. Works smooth, the goal is to generalise the input and make it flexible. Description of also: Memory allocation. ...
user avatar
  • 326
0 votes
1 answer
47 views

How can I correctly determine velocity of a point inside a grid after using mixed finite element method to solve Poisson equation?

I am using the mixed finite element method (MFEM) to solve the Poisson equation: $$\Delta h = 0,$$where $h$ denotes hydraulic pressure. The MFEM could determine the normal flux rate, $q_n$, through ...
user avatar
2 votes
1 answer
69 views

Shock Capturing Methods for Shallow Water Equations

I am looking for some help finding a numerical solution to the shallow water equations: $\partial_tu+\partial_x(u^2/2+g\eta)=0$ $\partial_t \eta+\partial_x(u\eta)=0$. where $u$ is the depth averaged ...
user avatar
1 vote
1 answer
162 views

How to implement Lax-Friedrich flux splitting with WENO scheme

I'm working on modeling a shock wave using the Euler equation with an advanced Equation of state and the fifth order WENO scheme. The equation are set up on the form: \begin{equation} \frac{\partial U}...
user avatar
  • 55
2 votes
0 answers
147 views

Recommended CPU for CFD running parallel, FDS (Fire Dynamics Simulator)

I am purchasing a new workstation to run FDS (Fire Dynamics Simulator) simulations, a CFD code for thermally driven fluid flow. Currently, I am using an Ubuntu Linux build with a Xeon E5-2630 v3 @ 2....
user avatar
  • 29
1 vote
0 answers
40 views

TVD slope / flux limiters formulation

Even if the formulation is the same the TVD slope limiter can be applied: to state reconstruction at the interface, in 1D FV formulation, we reconstruct the $Q^*_{j+1/2}$ and the $Q^*_{j-1/2}$ in the ...
user avatar
  • 382
0 votes
1 answer
143 views

Can I use periodic boundary conditions for `U` but not for `p`?

Cross-posted from Stack Overflow. (https://stackoverflow.com/questions/70686368/can-i-use-periodic-boundary-conditions-for-u-but-not-for-p) I am trying to numerically compute the drag force around a ...
user avatar
1 vote
1 answer
73 views

Direct integration of 2D Euler Equations with Runge Kutta shows oscillating Courant-Friedrichs-Lewy coefficient. Stiff or Bug?

By writing the direct integration of the 2D Euler Equations in a wide and short box where the fluid enters and exits through the horizontal faces using the Runge Kutta O(4) method I have found that ...
user avatar
  • 113
2 votes
1 answer
57 views

The effect of grid size on the total flux when solving Darcy flow with mixed finite element method

I am solving Darcy flow now with mixed finite element method. The Dary flow is $$\begin{equation}\begin{aligned}k^{-1}\mathbf{q} + \nabla h=0, \text{ in } \Omega\\ % \nabla\cdot \mathbf{q} = 0, \text{ ...
user avatar
0 votes
0 answers
25 views

Need help implementing constrained transport method for MHD riemann slver

I tried implementing constrained transport method to 1st order Godunov solver and when I run it on 2D rotor test case, the solver becomes unstable and crashed. The solver is able to maintain ...
user avatar
  • 101
1 vote
1 answer
66 views

What are simple ways to simulate gas expansion from a pulsed nozzle into a vacuum?

I am doing a physics experiment that involves releasing gas from a reservoir into a vacuum chamber via a pulsed nozzle and I am interested in knowing what would be a simple (and reasonably accurate) ...
user avatar
0 votes
1 answer
102 views

Implementation of mixed hybrid finite element method

The mixed hybrid finite element method (MHFEM) is based on the mixed finite element method (MFEM). So, I'd recall the implementation of MFEM. The mixed formulation of Poisson equation reads $$\begin{...
user avatar
2 votes
2 answers
70 views

Possible to use Iterative FD methods to solve a transformed non square domain [matlab]?

For the 2-D Poisson equation $$-(u_{xx}+u_{yy}) = f \ \ \text{where} f = 1$$ For boundary conditions $$\frac{\partial u}{\partial n} = 0 \ \text{on AB and AD}$$ $$ u = 0 \ \ \ \text{on BC and CD no-...
user avatar
  • 23
2 votes
1 answer
118 views

Mineral dissolution and solute transport around a solid

I am trying to simulate solute transport of acid (HCl) and consequent mineral dissolution around a grain (calcite). The governing equation for transport is the advection-diffusion equation, given as: ...
user avatar
  • 47
2 votes
1 answer
129 views

Elementary matrix of Raviart-Thomas elements

We can use the $RT0$ to solve the Darcy equation, i.e. $$k^{-1}\mathbf{u}+\nabla p = 0, \text{ in } \Omega,$$ $$-\nabla \cdot \mathbf{u} = 0, \text{ in } \Omega,$$ $$p = p_D \text{ on } \partial\Omega,...
user avatar
3 votes
2 answers
207 views

Test functions of Raviart-Thomas elements?

The test functions of general finite elements are like interpolation functions (if my understanding is correct). But how about test functions of Raviart-Thomas elements? Let's raise the $RT0$ element ...
user avatar
2 votes
2 answers
138 views

What is the most suitable numerical approach for modelling multiphase flow with particle interactions?

If I want to build a solver for this following problem: 1. There is stagnant water governed by the Navier-Stokes equation in the domain. ...
user avatar
3 votes
2 answers
201 views

How to enforce fluid and solid dynamic coupling in fluid-structure interactions using the finite element method?

I apologize in advance if the question has been posted before or if it sounds a bit naive. I am writing my own code in MATLAB for a staggered finite element solver for fluid-structure interaction ...
user avatar
2 votes
1 answer
66 views

Derivation of compressible volume-of-fluids formulation

I am trying to derivative the equations from [1] for a compressible Volume-of-fluids formulation but I am stuck in one of the last steps and would like to request some help to solve it. The governing ...
user avatar
2 votes
0 answers
72 views

Solute transport around a solid obstacle

I am a newbie in CFD and single/multiphase flow and transport in general. As part of my quest to learn, I am trying to model solute transport around a solid object in the center of a 2D domain. The ...
user avatar
  • 47
3 votes
1 answer
72 views

Instability at the boundary of a finite difference simulation of a hyperbolic PDE

I want to simulate the hyperbolic partial differential equation $$W_{tt} + V W_{tx} + k_E V W_x + k W_t = 0,$$ but I am having trouble finding a discrete analog of this equation which is numerically ...
user avatar
3 votes
1 answer
122 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
  • 590
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
57 views

Jos Stam's stable fluids — why is the timestep multiplied by the number of grid cells?

One more question about Jos Stam's GDC tutorial on stable fluids: in the advection step on page 8, the timestep for each dimension is implemented as dt * N, where <...
user avatar
  • 145
1 vote
1 answer
41 views

Jos Stam's stable fluids — why multiply by the grid size when subtracting out the curl-free part?

Considering page 10 of Jos Stam's GDC tutorial on stable fluids, the function project() first obtains the divergence of the velocity field: ¹ ...
user avatar
  • 145
0 votes
1 answer
89 views

Simulating the behavior of falling peas

I hope this is the right place for my question. I am a university student who designed a pea-pouring machine.Now I‘m really interested in simulating the behavior of those peas. Basically, they are ...
user avatar
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
1 vote
1 answer
150 views

Non-reflective boundary condition

I'm currently solving incompressible Navier-Stokes system of equations with periodic flow and high viscosity. Is there any outlet boundary types that avoids the reflection of flow from the outlet back ...
user avatar
  • 326
2 votes
0 answers
66 views

Multigrid method: linear solver and modified residual

I am trying to better understand the FAS multigrid algorithm for Euler equation in FV discretization. The usage of the modified residual (the residual with forcing) inside the different cases: ...
user avatar
  • 382
2 votes
2 answers
81 views

Filter size of large eddy simulation with Smagorinsky sub-grid stress model

I am trying to implement Large eddy simulation for solving air flow simulation with large Reynolds number using Smagorinsky sub-grid stress model. I read that the filter size should be calculated ...
user avatar
  • 23
1 vote
0 answers
100 views

Weird "oscillatory" modes appearing in FEM simulations

I am using COMSOL to solve a mathematical model involving thermoelectric hydrodynamic (TEMHD) flow. I am running a very large parameter sweep and using the solutions obtained to make some plots. ...
user avatar
5 votes
1 answer
187 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
3 votes
1 answer
209 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
1 vote
0 answers
53 views

How to numerically solve PDE that governs the free vortex wake model?

Crossposted at Math SE I am reading a paper on the free vortex wake model for a helicopter rotor blade, which is described by the following PDE: $$\frac{\partial \vec{r}}{\partial \psi} (\psi, \zeta) ...
user avatar
1 vote
2 answers
148 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
58 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
4 votes
1 answer
91 views

Why do we have to resort to Higher order schemes for solving the 1-D advection equation/ continuity equation?

\begin{equation} \begin{aligned} \frac{\partial N}{\partial t} &+ \frac{\partial J}{\partial r} = 0, \\ \frac{\partial N}{\partial t} &+ \frac{\partial }{\partial r}(N \upsilon ...
user avatar
0 votes
0 answers
36 views

Insert a boundary condition without removing periodicity assumption

I've to perform a multi-fluid internal flow simulation with a code which intrinsically assumes periodicity on a given direction (say z) on a cartesian grid. Nonetheless, the problem I'm trying to ...
user avatar
  • 1
-3 votes
1 answer
93 views

How to avoid negative concentration from numerical solution using FDM scheme?

$\frac{\partial C}{\partial t} + u \frac{\partial C}{\partial x} + w \frac{\partial C}{\partial x} = D \left(\frac{\partial^2C}{\partial x^2}+\frac{\partial^2C}{\partial y^2}\right)-C \cdot \left(\...
user avatar
0 votes
0 answers
335 views

Tool that calculates total drag for a submarine

Can you please suggest a tool (software or online tool) that can analyze a submarine hull shape and give the total drag at a specific velocity. A simplified tool would work too, where I select a ...
user avatar
3 votes
0 answers
110 views

Galerkin Least-Squares stabilization for FEM solution advection (hyperbolic) equations

I am playing with Galerkin Least-Squares stabilization to solve advection diffusion problem in the context of the finite element method. This works very well for steady-state advection-diffusion ...
user avatar
  • 1,127
7 votes
1 answer
202 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
  • 326
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
  • 590

1
2 3 4 5
12