Stack Exchange Network

Stack Exchange network consists of 174 Q&A communities including Stack Overflow, the largest, most trusted online community for developers to learn, share their knowledge, and build their careers.

Visit Stack Exchange

Questions tagged [fluid-dynamics]

The study of the properties of fluids and gases in motion

4
votes
2answers
720 views

How to discretize Burger's equation?

I am trying to solve the very simple one dimensional burgers equation which is: $$\frac{\delta U}{\delta t} + \frac{\delta F}{\delta x} = 0$$ where the flux F of some variable U is defined as$$ F= \...
0
votes
2answers
197 views

3D Divergence-free Stokes equations

I would like to have a divergence-free formulation of the Navier-Stokes equations for creeping flow (a.k.a. Stokes equation, N.-S. without the inertial terms) in 3D, for purposes of tracking particle ...
4
votes
0answers
144 views

Is it generally unstable to use in multidimensional simulations finite difference schemes with higher orders than 2?

I'm part of a team trying to generalize a 1D Advection-Diffussion-Reaction code we inherited by extending it to 2D by using dimensional splitting, i.e. solving advection and diffusion for x and y ...
6
votes
2answers
945 views

9-point stencil finite difference Laplacian with variable diffusion coefficients

So I'm trying to implement a 9-point stencil discretization to the 2D difussion equation. The stencil is here. However, most of the literature deals with a Laplacian that has a constant diffusion ...
4
votes
1answer
122 views

transverse component for multidimensional advection in method of lines

So I inherited from some people a code that solves the advection-diffusion-reaction equation for a particular system. The original code was first implemented in 1D which worked fine in cartesian ...
3
votes
1answer
365 views

Darcy flow finite elements

The Darcy equations for porous media flow are given by: $\frac{\mu}{\kappa}\mathbf{u} - \nabla p = \mathbf{0}$ $\nabla\cdot\mathbf{u} = 0$ where $\kappa$ is the permeability and can in general be ...
1
vote
0answers
341 views

ANSYS Fluent: Defining External Force on a Cell Zone [closed]

I am using ANSYS Fluent 15 to simulate a case of fluid flow. In my case, I have three cell zones in my simulation, let's call them top, middle and bottom. My requirement is to put a vertical constant ...
2
votes
1answer
107 views

CFD turbulence modelling mean pressures vs peak pressures

I have been experimenting with using Autodesk CFD to investigate facade/ cladding pressures on a (rectangular) building, comparing results with cladding/ facades pressures pressures from design codes ...
2
votes
1answer
429 views

What are acceptable boundary conditions for porous media flow?

I am attempting to simulate fluid flow through a porous foam. I would like to have no-slip boundary conditions on part of the boundary and free flow conditions on the inlet and outlet. Right now I am ...
3
votes
3answers
631 views

Where do I begin learning CFD?

I am working on some CFD codes - generating adjoint of a CFD model, for example. My background however is in computer science. I have absolutely no clue what the model represents or what it is doing. ...
0
votes
1answer
83 views

Visualise data under Tecplot

I have a set of data that I need to visualise under Tecplot, the file format is for example: ...
1
vote
0answers
99 views

Conservation at grid interface in adaptive mesh refinement

I am using adaptive mesh refinement to solve one dimensional inviscid Burgers equation. However I am facing some difficulty to handle grid interfaces which are not uniform (coarse-fine grid interface)....
2
votes
1answer
472 views

openfoam - Programming customized PDEs

I am looking for a method to automatically solve custom PDEs on a custom control volume. Specifically I would like to solve equations similar (but not exactly alike) to: $$\frac{\partial y}{\partial ...
3
votes
2answers
1k 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 ...
4
votes
3answers
2k views

How to deal with nonlinear term in Navier Stokes equations (finite element code)

I am trying to solve the Navier Stokes equations using the finite element method. I plan on using the pressure correction method to deal with the pressure and an implicit time stepping scheme for ...
0
votes
1answer
243 views

Adaptive mesh refinement basic conceptual problem

I am a beginner in adaptive mesh refinement (AMR). After I am done with the first two papers by Dr. Marsha Berger, I was trying to write my own code for a problem which has a parabolic partial ...
0
votes
1answer
79 views

Representation of polynomial order in CFD codes

I am currently working on a CFD code over a cubic grid. Now, the number of elements used in the simulation is decomposed among the number of processors. Each of those processors (a section of the cube)...
2
votes
4answers
1k 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 ...
5
votes
2answers
643 views

A Question About the Rhie-Chow Interpolation Used for Solving the Incompressible Navier-Stokes Equations on Unstructured Grids

When using the SIMPLE method on a mesh with a collocated variable arrangement, the following interpolation is used for the advecting velocities: \begin{equation} u_f = \overline{u}_f - \overline{D}_f\...
3
votes
2answers
285 views

Direct Numerical Simulation of reacting flows

I need to write a DNS code for simulating reacting flows in turbulent conditions. The code has to be highly scalable, because the computational cost of simulations is expected very large. My idea is ...
3
votes
1answer
557 views

2D Poisson Solver for Taylor Green Vortex Problem

I am trying to write a 2D Navier Stokes solver using an RK3 for time advancement in python. For debugging, I have converted the RK3 to an Euler step for simplicity. Checking my divergence for my ...
6
votes
2answers
405 views

applications of computational geometry in fields such as CFD?

Out of curiosity, I was recently trying to search what skills are required to be successful as developer in scientific computing field (e.g. CFD or similar). And to do so, I was going to through ...
1
vote
1answer
97 views

alternatives to moving mesh technique

If one has to simulate oscillating plate/solid wall ( sinusoidal function of time) in a domain (a simple piston movement in 'y' direction), the obvious way would be scenario 1. For brevity, assume ...
1
vote
3answers
1k views

2D Stokes equation Code

Does anyone know where could I find a code (in Matlab or Mathematica, for example) for he Stokes equation in 2D? It has been solved numerically by so many people and referenced in so many paper that I ...
3
votes
1answer
623 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 ...
1
vote
0answers
61 views

Help about Fluid-Fluid coupling techniques

I need a little help and advice with a project I want to do: the idea is to "couple" (I don't know whether I can call it like this) a conservative Navier-Stokes Solver (Fractional-step, 2nd order FDM) ...
0
votes
2answers
158 views

The “standard map” and Jupiter's red spot [closed]

Has there been an attempt to "explain" Jupiter's Red Spot via what is known as the standard map in chaos theory? I cannot be the first to notice the visual similarities...       &...
2
votes
2answers
439 views

How to interpret/convert density in 2D simulations?

Density has units of $kg/m^3$. In simulations often scaling is done to deal with non-dimensional values, so the physical density is converted to the non-dimensional density so that it can be used in ...
2
votes
2answers
131 views

Value of density when there are no or very few neighbours in SPH simulation

Sorry for the noob question. I am trying to implement SPH using the directions shown in this paper. The density needs to be updated using the formula The smoothing kernel is If there are no ...
4
votes
2answers
522 views

Manufactured solution for pressure based 3d incompressible Navier-Stokes solver with wall boundaries

I already successfully verified my solver (SIMPLE-type FVM-method) with the following manufactured solution (3d Taylor-Green vortex) on the solution domain $[-1,1]^3$ with Dirichlet boundary ...
0
votes
1answer
580 views

How to solve Energy Balance equation by numerical method

Good Day I am new to heat transfer technique please give me some suggestion on solving energy balance equation $$a \frac{\partial T_p}{\partial t}=\frac{\partial}{\partial x}\left(b\frac{\partial ...
2
votes
2answers
341 views

How to prove the strong formulation for Discontinuous Galerkin is conservative?

I've been trying to figure out how to prove that the following equation is conservative $$\int_{D^k} \partial_t(u^k) v_j^k + \partial_x(f) v_j^k dx = [(f -f^*) v^k_j]_{x^k}^{x^{k+1}}$$ where $v_j^k$ ...
2
votes
0answers
128 views

Arbitrary Choosing of the Solution Domain - Navier Stokes and Manufactured Solutions

I want to verify a finite-volume solver (SIMPLE-Algorithm) for the incompressible Navier-Stokes equations by using a manufactured solution. I use Dirichlet boundary conditions for the velocity at all ...
5
votes
3answers
376 views

Variable viscosity Stokes equation

One very efficient way to solve Stokes equation with periodic boundary conditions \begin{equation} -\eta \nabla^{2} \bf{v} + \nabla p = f \\ \nabla \cdot \bf{v} = 0 \end{equation} is using the ...
6
votes
2answers
1k views

Solution blows up when using Runge-Kutta to solve simultaneous ODEs for liquid film equations

This is a continuation of a problem I asked over at physics exchange and math exchange. Basically I have two ODEs that I am solving in order to calculate the radial and tangential velocity of liquid ...
2
votes
1answer
3k views

What's the difference between grid-based and mesh-based methods for PDEs?

I am into computational fluid dynamics and so far I've found that the most common approaches to solve for the governing equations are Eulerian and Lagrangian. The former samples the domain at fixed ...
4
votes
1answer
91 views

Zero-k mode in Pseudo-spectral solution of Stokes Flow

I'm trying to solve a Stokes flow problem with a pseudo-spectral method in periodic boundary conditions. The equations of interest are $-\nabla^2 \bf{v} + \nabla p = \bf{f} \\ \nabla \cdot \bf{v} = ...
6
votes
1answer
609 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 ...
2
votes
2answers
122 views

Expected computational time for DNS computation of fluid flow

Using an established criterion involving capturing eddies down to the Kolmogorov length scale it can be reasoned that the order of grid points in the computational mesh needs to be $N^3 \ge Re^{9/4}$ ...
3
votes
1answer
356 views

Finite Element integration with tensor notation

While I was studying discontinuous finite element methods I found an integration of a Navier Stokes equation using tensorial notation. The equation is the following: $\mathbf{\bar {u}}_{t} + (\...
0
votes
3answers
688 views

What Is The Intersection Between CFD and Computer Science

What Is The Intersection Between CFD and Computer Science... And.. How hard would it be to undertake a PHD in CFDs after an undergrad and masters in CompSci Many Thanks....
2
votes
2answers
250 views

Can the conservative form of the advection equation be re-written by replacing the velocity term with an integral over all other points in space?

Suppose I have a 1D advection equation in conservation (divergence) form $\partial_t u(x,t) = -\partial_x [v(x)u(x,t)],$ where $u$ is a conserved quantity in space, and $v$ gives the velocity of the ...
1
vote
1answer
114 views

Finite Volume Method flux integration

I was reading through this document about FVM. I understood all up to the point where we have on page 15 the following $$(\bar u^{n+1}_i - \bar u^n_i) \Delta x + \int^{t^{n+1}}_{t^n} f(u(x_{i+1/2},t))...
1
vote
1answer
247 views

What kind of test cases are convenient to use for testing the code for Euler equations of gas dynamics in polar coordinates?

Consider Euler equation of gas dynamics in polar coordinates as $$ \left( \begin{array}{ccc} \rho \\ \rho u_{r} \\ \rho u_{\theta} \\ E \end{array} \right)_{t} + \frac{1}{r} \left( \begin{array}{ccc}...
0
votes
1answer
394 views

CFD: Mesh problem [closed]

I wish to study the 2D heat transfer of the soda can using Fluent. For my geometry I took a circle. My problem while meshing the circle; Because i need to mesh it from inside. I used the mapped faced ...
1
vote
0answers
280 views

Two Dimensional Nozzle Flow

I've been looking at the book "Computational Fluid Dynamics - The Basics With Applications", by John D. Anderson, and his solution for the one-dimensional subsonic-supersonic nozzle flow using the ...
2
votes
1answer
109 views

Trignometric/Fourier spectral collocation with zero Dirichlet BC in 2D

I am concerned with numerical solution to the following problem on $[0,1]\times[0,1]$. $\dfrac{\partial\theta}{\partial t}+u(x,t).\nabla \theta(x,t)=\kappa \nabla^2\theta(t,x)$ with Dirichlet ...
1
vote
2answers
5k views

Fluent time step size (transient implicit solver)

I am studying fluid-structure interaction across integrated circuits packaging by using Ansys Fluent (I'm a beginner). Currently I had an issue about the time step size. When I set time step size of ...
3
votes
0answers
287 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 ...
3
votes
2answers
670 views

How to set up a shock tube problem such that the solution includes a shock with a specified Mach number

One of the famous and convenient test cases for shock wave modeling is the 1D Sod's shock tube. This is a Riemann problem for the compressible Euler equations of gas dynamics. The initial set up has ...