Questions tagged [fluid-dynamics]
The study of the properties of fluids and gases in motion
631
questions
1
vote
1
answer
82
views
Symmetry axis boundary condition
I was wondering about the symmetry axis boundary condition in commercial CFD solvers such as ANSYS Fluent.
If the problem is the flow through a round pipe or out of a round nozzle, it is natural to ...
2
votes
3
answers
199
views
Approximation of derivatives in the finite volume method
I am looking into the finite volume method and I have come to a problem with discritisation. Suppose I am looking at a particular cell(I'm dealing) with cartesian grid. at the points $(X_{i},Y_{j}),(...
1
vote
1
answer
77
views
Is the weak form on this book about a level set FSI problem wrong?
I'm trying to repeat a level set FSI problem on the book <Level Set Methods for Fluid-Structure Interaction>, on the page 89, the provided freefem code define a weak form of the discretized ...
0
votes
0
answers
18
views
Is there any papers about simulating material damage under supersonic combustion?
I'm quite interested in this field, Is there any papers about simulating the damage of material(such as alloy) under supersonic combustion, after searching for a while, I didn't find anything
0
votes
0
answers
25
views
Time step not converging in Transient simulation with RPI Wall boiling
I am trying to simulate subcooled flow boiling in horizontal channel with Non-equilibrium RPI wall boiling model. It is transient simulation with Implicit scheme. Steady state simulation are not ...
0
votes
0
answers
109
views
Where has my implementation of the SHAS2D (Book and Boris, 1975) model gone wrong?
I am new to computational fluid dynamics, and am having a tough time implementing the Book and Boris flux limiter SHAS2D model.
I am attempting to re-create the SHAS2D fluid code to solve the Euler ...
2
votes
1
answer
100
views
FENICSx or deal.ii for modeling fluid structure interaction of cylindrical energy harvesters in various arrangements
I'm trying to model fluid-structure interaction of cylindrical energy harvesters in various arrangements. I have chosen OpenFOAM for the fluid part but I have no clue which framework to use for the ...
2
votes
1
answer
101
views
AMR framework for efficient simulation of PDEs, potentially with boundary layers
I am aiming at developing a simulation software mostly aimed at reactive fluid flows applications. Here are a simplified list of points to drive the initial choices (code structure, libraries...):
I ...
0
votes
1
answer
75
views
In level set fluid structure interaction method why can we rewrite the elastic force in this form while there is no shear force?
when we consider a Immersed Membrane Without Shear, we can define a regularized elastic energy as $$\mathcal{E}(\varphi)=\int_{\Omega} E(|\nabla \varphi|) \frac{1}{\varepsilon} \zeta\left(\frac{\...
0
votes
0
answers
28
views
OpenFOAM- write interval changed after restarting simulation
I am running a pimpleFoam simulation on a supercomputer. I have to submit it as a job which has an allocated time. The simulation I was running ran out of this allocated time. So then I restarted it ...
0
votes
1
answer
111
views
Lattice Boltzmann method parallelization
I'm learning about the lattice Boltzmann method (LBM) and how it can be parallelized. On Wikipedia, they say, "The LBM was designed from scratch to run efficiently on massively parallel ...
0
votes
0
answers
57
views
What happens at the interface between a solid and a fluid?
I am doing an MPM simulation of water colliding against a solid, I am currently encoding the model for collision forces and I am realising something isn't clicking.
Let us assume the fluid experiences ...
0
votes
0
answers
26
views
Error - Getting negative temperature while solving Quasi 1D flow in a Nozzle (JD Anderson Book)
I am currently trying to solve Nozzle Flow problem from JD Anderson Book.
Here's my entire python code -
...
1
vote
2
answers
189
views
How should I solve generalized eigenvalue problems in Python? (Orr-Sommerfeld equation)
I am trying to solve the Orr-Sommerfeld equation numerically, using the techniques given in this article. This leads to solving a generalized eigenvalue problem, that is, given two matrices $\mathbf A,...
0
votes
0
answers
28
views
How can i model upward natural convection in various angles (0 - 90)?
I am looking for a way to numerically solve the Naiver-Stokes equations for steady incompressible flow using FDM over a surface with various angles?
0
votes
1
answer
43
views
How to address the element face adjacent to boundaries when the finite difference method and marker-and-cell scheme are used to solve the Stokes flow?
The Stokes equations are
$$-\Delta \mathbf u + \nabla p = f \text{, in }\Omega,$$ and
$$ -\nabla \cdot \mathbf u = g, \text{ in } \Omega$$
where $\mathbf u = \left( u, v \right)$ is the flow ...
1
vote
1
answer
120
views
How can we use the two boundary conditions for the Taylor-Maccoll Equations for Cone Shock Waves?
The Taylor-Maccoll equations below govern the gas dynamics around a supersonic, axially oriented cone.
Using the notation from Anderson's Fundamentals of Aerodynamics, which uses a dimensionless ...
0
votes
0
answers
19
views
MPM course, implicit updat clarifification
In the 2016 SIGGRAPH notes on MPM the authors explain how to do an implicit update in section 11.4 throguh by this formula:
It is not clear to me what the $j$ is over in this case (for eq 202), is it ...
0
votes
0
answers
14
views
Resultant Lift Direction Confusion
Quite a small scale question but I have gotten myself confused. Inside Ansys Fluent I am simulating a wing at different angle of attacks, instead of rotating the airfoil I have rotated the inlet ...
1
vote
0
answers
116
views
Why do we use modified pressure in incompressible multiphase solvers with gravity?
The context of my question is two-phase incompressible solvers such as interFoam in OpenFOAM, but I have seen this trick used ...
0
votes
1
answer
60
views
Fluent Induced Drag Aerofoil Convergence Error
I am currently modelling what happens, to lift induced drag, when you apply washout to a wing (Twist). However I have a few issues
the drag values I calculate are into the 1000s
the drag values ...
0
votes
0
answers
39
views
Energy potential involvign pressure and viscoscity?
Making a really long explanation short. I am looking for an energy density function for an incompressible isotropic fluid that involves both viscoscity and pressure.
The reason is, there is a way to ...
0
votes
0
answers
54
views
Solving a minimization problem without flattening inputs?
I am reading this paper on improving time steps for solving simulation problems:
https://www.math.ucla.edu/~jteran/papers/GSSJT15.pdf
The authors developed this energy function:
$$E(x) = \frac{1}{2\...
0
votes
0
answers
30
views
MPM implicit integration
Setup
I am having a hard time translating the formulas in the 2016 MPM course from siggraph into algorithmic computations.
In particular we have equations 200 and 201
$$h(v^{n+1}) = Mv^{n+1} - \Delta ...
3
votes
1
answer
163
views
Projection (or fractional-step) methods Vs coupled method for incompressible Navier-Stokes
My question is in the context of the finite element method.
Incompressible Navier-Stokes equations can be solved using the coupled method or projection/fractional step methods. Each method has its ...
0
votes
0
answers
39
views
Refluxing step on Finite difference AMR
Hi I am a computer scientist working on MHD code for astrophysics simulation. We use a finite difference scheme where we first solve the spatial derivatives and with them solve the right hand side and ...
2
votes
1
answer
120
views
How do the navier stoke equations model materials who "forget" their original form?
Sorry for the screenshot but I don't want to try to format this on latex:
We have this annotation of the Navier-Stokes equations:
I am particularly puzzled by the viscosity/stress term. For an ...
0
votes
0
answers
91
views
Adding stability to MPM simulation?
I am writing a 2D implementation of MLS-MPM, I have fluids working perfetly fine, solids technically work as well, at low time steps.
This is the fluid simulation at a large time step:
https://i.stack....
0
votes
0
answers
36
views
Matrix for Marker and Cell grid?
I have an assignment question that reads:
Show that the combined matrix for the Marker and Cell (MAC) grid
for velocity and pressure for the steady Stokes equations is symmetric.
You can consider the ...
0
votes
0
answers
53
views
Help with code for FDM Couette flow of two fluids on python
I'm trying to solve numerically with a python code a very simple BVP that represent the velocity field of two fluids between two plates where the driving force of the momentum transport is only for ...
1
vote
0
answers
60
views
SPH example codes for solving differential equations
I'm quite new to the smoothed particle hydrodynamics method, so I'm trying to find Matlab or Python code to help me understand it. Specifically, I want to know how to solve the one-dimensional ...
0
votes
0
answers
55
views
Solving Laplace for Velocity Potential in Constricted Channel
I am trying to solve the 2D Laplace equation numerically to give the velocity potential of a fluid flowing in a channel with a constriction:
$u_{xx}+u_{yy}=0$
There is a constriction in the channel at ...
1
vote
0
answers
52
views
How to calculate the reynolds number for airborne wind energy system?
I am asking for some help on how to calculate the reynolds number for an AWES. I know the formula of Re=rhovl/viscosity. I have selected the air foil as Clark y with chord length 3.72m (which will be ...
0
votes
0
answers
119
views
discrete definition of curl $ \nabla \times \vec{A}$ on a 3D grid?
I have the data for 3D vector field $\vec{A}$ (with components $\vec{A_1}$, $\vec{A_2}$ and $\vec{A_3}$) sampled on a 3D grid with integer indices i, j and k.
Assuming that only the third component $\...
1
vote
0
answers
69
views
Fortran - Lid-Driven Cavity Boundary Conditions Error when using SIMPLE method
I am studying Numerical Methods for incompressible flows. part of the tasks is to model the lid driven cavity problem in 2D using the SIMPLE method.
I have been provided with Fortran code that is ...
0
votes
0
answers
39
views
books/paper recommendation on computational thermal-turbulence by using FEM
I have just learned basic FEM for 2D N-S euqation, now my teacher let me to do the following problem, the document of this problem is in large fluid problem, the system of equations is listed in that ...
0
votes
0
answers
39
views
What is the boundary condition of $\varepsilon $ for $k-\varepsilon $ turbulence model?
I'm working on the numerical solution of systems of equations, you can use this link to access it, the system of equations is
$\begin{aligned} \partial_t \theta+u \nabla \theta-\nabla \cdot\left(\...
0
votes
0
answers
76
views
Book recommendation on multiphysics
I want to learn multiphysics such as fluid-structure interaction where the simulation is performed for heat transfer, fluid dynamics, and electrodynamics.
Could you please recommend some books about ...
1
vote
1
answer
54
views
How to implement pressure stabilization in matlab while solving steady stokes equation
the equation is $$
\left\{\begin{array}{l}
-\nabla \cdot \mathbb{T}(\mathbf{u}, p)=\mathbf{f} \text { in } \Omega, \\
\nabla \cdot \mathbf{u}=0 \text { in } \Omega, \\
\mathbf{u}=\mathbf{g} \text { on ...
0
votes
0
answers
49
views
How to specify the value of pressure in a point in solving steady stokes equation numerically?
I try to solve the steady-Stokes equation numerically on , that is
\begin{aligned}
-\nabla \cdot \mathbb{T}(\mathbf{u}, p) & =\mathbf{f} \quad \text { on } \Omega=[0,1]*[-0.25,0], \\
\nabla \cdot \...
0
votes
0
answers
63
views
Navier Stokes Equation Discretization by upwind scheme
I've been reading this paper regarding the SIMPLER algorithm in CFD and I don't understand exactly how the discretized equation mentioned by the author is arrived at. It says that the convective terms ...
0
votes
0
answers
49
views
Parallel Plate Flow simulation using SIMPLER Alogortithm
I am new to the community and CFD itself and am working on my first CFD program in Scilab based on the SIMPLER algorithm in Scilab i.e. modeling a flow between parallel plates with entrance region. I ...
0
votes
0
answers
110
views
Once a method is implemented in a computer program, then does learning the theory become less useful use of time?
Once a method is implemented in a computer program, then does learning the theory become less useful use of time?
This has confused me. I used to view that all things should be studied from first ...
0
votes
1
answer
60
views
How to combine multigrid preconditioner with jacobi preconditioner?
I have not found any relevant information in the literature on the following rather simple problem:
How to combine (geometric) multigrid preconditioned conjugate gradient (MGPCG) with an additional ...
2
votes
0
answers
111
views
Starting point to use Trilinos linear algebra packages?
I have experience with classical linear algebra packages in C++ like Eigen, Blaze, etc.
I have never wrote my own PetSC back-end solver but I used/modified several of them.
In a new project I would ...
0
votes
1
answer
119
views
Problem with my Octave code (unsteady heat equation with FEM)
I want help with my Octave code regarding the unsteady heat equation.
My geometry and mesh are generated with FreeFEM++, so there is no problem with that (I tried it with the steady problem with no ...
3
votes
0
answers
934
views
What EXACTLY are the "outflow" and "pressure outlet" BCs in ANSYS Fluent?
I am just starting my MSc and we have been given an assignment in which we must model two-dimensional, laminar, incompressible Poiseuille flow in a channel. We are supposed to use a VELOCITY INLET ...
0
votes
1
answer
106
views
Aerofoil study using CFD, struggling to find aerofoil coordinates
I’ve been messing around with Ansys and I’m struggling to find the aerofoil coordinates for a NACA 66-012?
I looked on Airfoil tools, but it doesn’t allow you to generate a 6 series aerofoil, only 4 ...
1
vote
1
answer
192
views
Differential equation for radioactive cooling in fortran
Today I'm trying to evaluate this differential equation for internal energy in a gas in Fortran:
$$
\frac{du}{dt} = - \frac{n_H^2}{\rho}\frac{\Lambda}{n_H^2}
$$
Where nH is the density of hydrogen in ...
1
vote
0
answers
58
views
Anisotropic lines identification algorithm
I am implementing an algorithm to identify the anisotropic lines in an unstructured mesh. This is done in the framework of an Agglomeration algorithm for Non-Structured multigrid applied to FV ...