Questions tagged [fluid-dynamics]

The tag has no usage guidance.

Filter by
Sorted by
Tagged with
2 votes
0 answers
46 views

Can Mathematica's FEM solve coupled Dirichlet Boundary Conditions?

I am solving a coupled system of PDEs using Mathematica's FEM capabilities. Specifically, the Navier-Stokes equations with a no-flux stress-free boundary. To do this, I need to specify a Dirichlet ...
Eric Hester's user avatar
7 votes
0 answers
251 views

Code for quasi 1D nozzle flows

The quasi-one-dimensional model describing the flow of compressible gas in rocket nozzles is very common. The corresponding equations have a divergent non dimensional form $ \frac{\partial \mathbf{U}}...
Alex Trounev's user avatar
  • 41.1k
2 votes
1 answer
126 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 ...
mike's user avatar
  • 303
3 votes
1 answer
145 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&...
Chris's user avatar
  • 33
2 votes
1 answer
77 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 ...
user1234321's user avatar
1 vote
0 answers
74 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 ...
Lucas's user avatar
  • 103
3 votes
1 answer
129 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 ...
XU KANGYOU's user avatar
1 vote
0 answers
77 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. ...
Lucas's user avatar
  • 103
5 votes
1 answer
103 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 ...
user21's user avatar
  • 39.1k
0 votes
0 answers
56 views

How to solve thin-film PDE with pressure singularity?

I would like to solve the following PDE representing the height of a thin fluid film: However, due to the D[p,x] term in the denominator, in regions where ...
Ariana Fenris's user avatar
0 votes
0 answers
54 views

PDE System solves in Cartesian but not in Cylindrical Axisymmetry

As a follow-up to my previous question, I am now solving the same system in axisymmetry. This introduces a singularity at the origin. Following various posts on solving in cylindrical coordinates, I ...
Ariana Fenris's user avatar
4 votes
1 answer
368 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 ...
Avrana's user avatar
  • 297
4 votes
1 answer
174 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 ...
Ariana Fenris's user avatar
7 votes
1 answer
333 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 ...
Avrana's user avatar
  • 297
2 votes
1 answer
207 views

Speeding up NDSolve to reasonable speeds to solve a coupled PDE system

Problem Statement I am planning to solve a PDE system which consists of a fluid droplet spreading on a non-Newtonian substrate. The system consists of the following equations: $$\frac{\partial p_1}{\...
FLP's user avatar
  • 475
1 vote
1 answer
166 views

How to use VectorPlot3D to visualize a solution of Navier-Stokes equations in cylindrical coordinates?

I try to solve axially symmetric Navier-Stokes equations in cylindrical coordinates. I have successfuly solved the equations, but I have troubles when I try to plot the solution, as VectorPlot3D does ...
lemurman's user avatar
  • 402
5 votes
0 answers
169 views

Resolving singularity in convection-diffusion equation using pdetoode

Building on the system of equations in this post, I attempted to solve an additional convection-diffusion equation describing the concentration of solute in the lens, which affects its spreading. ...
FLP's user avatar
  • 475
2 votes
1 answer
140 views

Directly calculating the cyclic steady state of a time-periodic conjugate heat transfer problem

Context The following transient problem is the reciprocating (i.e., fully reversing) flow of a fluid $0<x<L, 0<y<d$ over a thick heated block $0<x<L, -e<y<0$ until the system ...
Avrana's user avatar
  • 297
3 votes
1 answer
136 views

Conjugate heat transfer modelling of reciprocating flow crashes for long flow times

The following transient problem is essentially the reciprocating (i.e., fully reversing) flow of a fluid over a thick heated block until the system reaches a cyclic steady-state (i.e., the system ...
Avrana's user avatar
  • 297
20 votes
4 answers
1k views

What's wrong with this FFT-based Von Kármán vortex street simulation?

About 9 years ago, I came across this interesting website, and found the following paragraph with a broken Mathematica code sample: When fluid passes an object, it can leave a trail of vortices ...
xzczd's user avatar
  • 62.8k
2 votes
1 answer
113 views

Preventing Mathematica from applying product rule

I am programming viscous terms of compressible Navier-Stokes equations ...
dylewiczk's user avatar
4 votes
1 answer
137 views

NDSolve for Fluid Flow - Monitor Residuals

I'm using NDSolve for fluid flow and would like to monitor the convergence of the solution more closely. I have a simple stationary case set up. NDSolve should return the solutions for u,v and p. ...
kienappel's user avatar
  • 411
9 votes
1 answer
542 views

Fourth order, nonlinear PDE depicting dynamics of liquid lenses

PROBLEM STATEMENT Recently, I was trying to verify the solution of a set of thin film spreading equations given by this paper. Where $$p1=-(h_{1xx}+\sigma h_{2xx})$$ $$p2=-\sigma h_{2xx}-\Pi[h2-h1]$$ ...
FLP's user avatar
  • 475
1 vote
1 answer
101 views

Transient darcy flow equation

I can solve the following stationary Darcy problem: $div(\nabla u)=1$ with $ u = 0 \in \Gamma_d$ ...
Stratus's user avatar
  • 2,882
0 votes
0 answers
73 views

Plotting contours of a two-variable function containing a sum

I'm trying to use Mathematica to plot contours of a rather intricate two-variable function. The equation describes the velocity profile for laminar flow in a tube of rectangular section, namely: $$ {w^...
Lucas's user avatar
  • 21
0 votes
0 answers
120 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 ...
Noor Aslam's user avatar
10 votes
2 answers
818 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 ...
Avrana's user avatar
  • 297
20 votes
1 answer
487 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 ...
Alex Trounev's user avatar
  • 41.1k
4 votes
1 answer
168 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 ...
ASHISH's user avatar
  • 141
0 votes
0 answers
90 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 ...
dr_strangeloop's user avatar
1 vote
0 answers
149 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 ...
metudust's user avatar
1 vote
1 answer
91 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,...
ASHISH's user avatar
  • 141
5 votes
2 answers
301 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 ...
ABCDEMMM's user avatar
  • 1,816
17 votes
1 answer
599 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 ...
Alex Trounev's user avatar
  • 41.1k
16 votes
2 answers
990 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. ...
Alex R's user avatar
  • 163
4 votes
1 answer
520 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 ...
Alex R's user avatar
  • 41
1 vote
1 answer
70 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]$$
Zafar's user avatar
  • 61
4 votes
2 answers
297 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 ...
Adrian del Ser's user avatar
0 votes
0 answers
95 views

How to solve this linear complementarity problem with NDSolve or other functions? The arguments should be ordered consistently

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 ...
Hot Pizza's user avatar
0 votes
0 answers
80 views

Nonlinear Differential Equation with Paramaters in Sqrt Function

For a university project, I am trying to see if my system will have choked flow and also plot the resulting pressure spike. I set up the system below to try to model the transient response. I am able ...
OK1991's user avatar
  • 37
0 votes
1 answer
102 views

Using Solution to Differential Equation in New Function

I am trying to feed the solutions from my differential equations into another function. I am using NDSolve to solve two differential equations simultaneously. It ...
OK1991's user avatar
  • 37
6 votes
0 answers
216 views

System of coupled nonlinear pdes with FEM (magneto-hydrodynamic mass transfer)

EDIT 01 See end of post for an update. Original Question I'm trying to reproduce the magneto-hydrodynamic flow studied here using the nonlinear FEM functionality, and have been having trouble ...
George Varnavides's user avatar
0 votes
1 answer
204 views

FEM analyis of flow in a hydrocyclone [closed]

I am interested in solving the flow in a hydrocyclone (steady Navier Stokes equation). I have an STL file of the 3D geometry. My question is how to extract the internal geometry of the hydrocyclone ...
user2088899's user avatar
1 vote
0 answers
222 views

Fluid dynamics collapse 1D Euler equations with NDSolve

I am getting a solution with NDSolve for studying the collapse of an interstellar isothermal cloud. But I have two problems and questions: I am getting some warnings even if I am obtaining a good ...
Hergé's user avatar
  • 73
0 votes
2 answers
123 views

Better method for computing streamlines

Can someone suggest a better way to solve for streamlines? P.S. StreamPlot is not providing accurate results. Here, velocity is...
Anubhav Kamal's user avatar
0 votes
0 answers
91 views

please help me .. can any one solve this problem usingg NDSolve [duplicate]

solve for dependent variable f and theta. where m and n are parameters please help me. ...
muhammad dawood's user avatar
0 votes
1 answer
146 views

How can i solve and plot this non linear ODEs usin NDSOLVE?

I am trying to solve the following non-linear ODEs ...
muhammad dawood's user avatar
7 votes
2 answers
1k views

Cylindrical coordinates in FEM

I am trying to solve the Stokes equation for fluid flow in a 3d cylinder. All boundaries are no-slip, apart from the top boundary, which enforces flow in the x-direction. My problem is that I can't ...
Oscillon's user avatar
  • 1,231
12 votes
4 answers
841 views

Area/Volume of a 2D/3D object as it is filled up with water

How could I get the volume of water on a region as it is filled up with water from below? Assume gravity points in some appropriate downward axis like -y (or -z in 3D) so the water fills upward from ...
flinty's user avatar
  • 23.1k
6 votes
1 answer
405 views

FEM: periodic solution of 2D Navier-Stokes equations

Let’s consider a horizontal channel with a round obstacle in the middle. ...
Rodion Stepanov's user avatar