# Questions tagged [fluid-dynamics]

The tag has no usage guidance.

104 questions
Filter by
Sorted by
Tagged with
357 views

### Shallow-water equations on a rotating sphere by using the NDSolve method with a purpose of earthquake-generated tsunami simulation

I'm trying to numerically solve a system of Shallow Water Equations in Wolfram Mathematica 14.0.0 on the rotating sphere by using the NDSolve method with a purpose of earthquake-generated Tsunami ...
• 71
235 views

### Speed up finite element method solver with small features

I post a question on making mesh for small features. @user21 gave a useful suggestion which resolves the small features. To improve the mesh, I have also tried to increase the resolution near the ...
• 641
1 vote
56 views

### How to constrain two PDE variables to be the same value on boundary?

I have a 4 pdes for 4 variables: v, v2, ua, and p. I would like to specify that v==v2 on the left and right boundaries at all times. Here are the equations: I discretise the system using PDEtoODE. ...
1 vote
69 views

### Why is NDSolve ignoring two PDEs out of three ones I am solving?

I am solving the following coupled system of 3 PDEs modelling a 1D membrane coupled to a 1D fluid flow field underneath. However, on putting them into NDSolveValue (and trying with FEM), it says the &...
174 views

• 47
311 views

### How to improve FDM solver for unsteady viscous flow?

To solve the problem that is discussed in the paper Finite Difference Analysis of Time-Dependent Viscous Nanofluid Flow Between Parallel Plates we developed FDM solver based on the code from the blog ...
• 47.7k
199 views

• 47.7k
148 views

### Infinite expresssion 1/0 encountered when use NDSolve for 3D axisymmetric Navier-Stokes (Euler) equations

The PDEs we are interested in solving using NDSolve is the vorticity-stream formulation of the 3D axisymmetric Navier-Stokes (Euler) equations (Ref.1 :T. Y. Hou, Potential singularity of the 3D Euler ...
• 313
171 views

### Invalid PeriodicBoundaryCondition for Navier Stokes problem using NDSolveValue

Background: I'm looking to have a 2D re-entrant channel of a geophysical flow that is forced by a wind stress at the surface ("taux"; below) and experiences planetary rotation ("cor&...
• 33
96 views

### How to extract gauss coordinates in AceGen/ AceFEM?

I am trying to implement a 2D turbulence problem with the model of Spalart-Allmaras. Therefore, I need to extract the global gauss coordinates to define parameter 'd' which stands for the distance to ...
1 vote
88 views

### NDSolve stops solving past a certain time

This is a continuation of the post I've made Unable to solve Delay PDEs Error in Boussinesq Approximation. I apologise if I shouldn't have posted a seperate question for this but I think that the ...
• 103
230 views

### Steady state fluid flow for downward flow past an obstacle

I was attempting to simulate fluid flow past a circular obstacle. The following is the code which I used ...
1 vote
85 views

### Unable to solve Delay PDEs Error in Boussinesq Approximation

I'm trying to solve the set of equations below describing the flow of a pot of water being heated slightly. The equations are 2D axisymmetric in nature. ...
• 103
149 views

### Spalart-Allmaras turbulence model

In this post Alex gives an implementation of the Spalart-Allmaras turbulence model [1, 2]. The example produces reasonable results, as far as I can tell. However, the implementation Alex uses deviates ...
• 40.4k
464 views

### Implementation of FEM to steady-state coupled fluid flow and heat transfer

This is a steady-state conjugate heat transfer problem (the time-independent version of this problem). The problem is conjugate as the energy equation is being solved in thermally connected solid and ...
• 297
200 views

### NDSolveValue gives up without warning or solution on PDEtoODE system

I have the following pair of PDEs that I want to solve in the half-space x>=0:  I followed this post to decompose the 4th-order height equation into two 2nd-order ones to respect the hyperbolic ...
372 views

### Refining mesh size leads to absurd results for a coupled heat transfer FEM model

I have been recently solving a conjugate heat transfer problem, which involves fully-reversing or reciprocating flow of fluid over a heated block of solid. The problem is 2D and the temperature field ...
• 297
271 views

• 21
146 views

### How to find non-trivial solution to homogeneous linear equations by Mathematica using code (Input form)?

I have many homogeneous equations to which I am interested to find the non-trivial solution using Mathematica, But I am totally new to Mathematica, so I do no know, how to find the non-trivial ...
• 223
872 views

### Reciprocating flow in a channel over a heated surface

The following is a coupled heat transfer and fluid flow problem. A thick plane channel is being heated with a constant flux from the bottom (at $y=-e$) with a constant heat flux $q$ as shown in the ...
• 297
549 views

### 3D stable fluids algorithm to simulate transition from laminar to turbulent flow

This algorithm is 3D extension of our 2D algorithm published on this page and here. We suppose that with this code we can simulate transition from laminar to turbulent flow. In this example we compute ...
• 47.7k
173 views

### How to display warning message inside manipulate?

Using the code shown, Reynolds number is calculated for given parameters. Is it possible to display a warning message to the user if a parameter is outside applicable range? e.g., if the fluid ...
• 141
93 views

### Can mathematica numerically solve this system of coupled integro differential equations?

Here's my system: R(x,y,t) is the concentration of a chemical that's diffusing across the box. The disk is a mobile sink for the chemical, which absorbs it in proportion to its concentration. The ...
1 vote
176 views

### Interaction between fluid flow and elastic obstacle (fluid-structure interaction) [closed]

I would like to model the velocity and pressure distribution of an elastic obstacle in a channel where the fluid flows. There are Navier-Stokes examples in Mathematica webpage with rigid obstacles but ...
• 19
1 vote
96 views

### How to make a calculator to display the Non-Dimensional number

I am a beginner here and want to learn this tool. At first, I want to make one calculator that will display the value of the Non-Dimensional number of Fluid dynamics based on user inputs. In this code,...
• 141
335 views

### Modified Heat Transfer in Fluid Flow

I am trying to simulate Modified Heat Transfer in Fluid Flow (based on Buoyancy-Driven Flow in a Square Cavity ). The modified heat transfer takes the form: with the solid volume fraction: The ...
• 1,854
669 views

### Stable fluids code for electromagnetic mixture application

This code has been translated from the original Jos Stam code and improved with some Mathematica functions. It solves problem of viscous incompressible flow with electromagnetic force in a rectangle ...
• 47.7k
1k views

### Stokes equations in 2D with traction boundary conditions

This is a follow-up to a previous question (see here). We would like to solve the two-dimensional Stokes equations using the FEM package in Mathematica, when we prescribe traction boundary conditions. ...
• 163
586 views

### Stokes equations in 2D with Neumann conditions

I would like to solve the 2D Stokes equations within a unit disk, say $\Omega$, by using the finite element method (FEM) as it is implemented in NDSolve (by loading ...
• 41
1 vote
71 views

### How to find invariant sets of nonlinear operators below

How to find invariant sets of nonlinear operators below. $$x[n+1]=0.9∗y[n]^2/(y[n]+2)+x[n]∗(1−0.8/(1+x[n]))$$ $$y[n+1]=0.8∗x[n]/(1+x[n])+0.6∗y[n]$$
• 61
341 views

### 2D incompressible flow solver with known non-zero initial velocity distribution

I know that there are lots of questions on this forum regarding mathematica FEM, but I haven't quite found what I am looking for. I will be working with data obtained through PIV (particle image ...
I want to find a numerical solution of $\omega(t,q)$, where $t\in[0,T]$ and $q\in[0,N]$. $q_t$ is a determinstic w.r.t $t$, other unmentioned parameters can be fullfilled or given by the user. I met ...