Questions tagged [finite-volume]

Referring to the discretization of partial differential equations using Finite Volume Method.

Filter by
Sorted by
Tagged with
3
votes
1answer
55 views

How does the diffusion of a finite volume method with a WENO scheme compare with that of spectral methods?

I know that, in general, finite volume (FV) methods are more (numerically) diffusive than spectral methods. However, I can't find any information on how the advection scheme changes that. For example, ...
1
vote
2answers
103 views

How to apply central difference to viscous fluxes in 2D Navier-Stokes equations?

I'm trying to solve 2D unsteady compressible Navier-Stokes equations with finite-difference or finite-volume method. Here is the system, it's pretty standard: $$ \frac{\partial U}{\partial t} + \frac{...
1
vote
1answer
37 views

How to evaluate the average value of a polynomial inside the triangle area in finite volume sense?

Consider we have a linear bivariate polynomial: $$p(x,y)=ax+by+c.$$ To construct the linear polynomial using least square method, we need to evaluate the value of the average polynomial $p$ in at ...
1
vote
1answer
59 views

Where could error terms that blow up in SWE come from?

I have been working on a solver for shallow water equations with reflective boundary conditions. I have found that it diverges very fast. As a workaround I noticed yesterday that if I smooth the ...
1
vote
0answers
52 views

Finite volume reconstruction techniques

For a cell-centered finite-volume calculation, with a nontrivial stencil (either using an unstructured grid or a strongly non-uniform structured grid), what are the main techniques for reconstruction ...
3
votes
1answer
114 views

Is upwinding needed for slope limiter / flux limiter and numerical flux?

I have a cell centered cartesian grid and am trying to implement the flux inside the divergence term using numerical flux with a flux limiter. I found different formulas for MUSCL flux limiter, where ...
2
votes
2answers
111 views

Simplest way to “upgrade” from Euler equations to Navier-Stokes equations in FV or FD framework

I have quite a lot of experience solving unsteady Euler equations, including multi-component ones, with in house-coded finite-difference and finite-volume methods, including MacCormack and MUSCLE ...
1
vote
0answers
30 views

Solution of non-linear Poisson equation does not match reference

I'm trying to solve the non-linear Poisson equation as a first step to solve the drift-diffusion equations for semiconductors. For reference, I'm using a preprint from the Weierstrass Institut (which ...
0
votes
0answers
51 views

Simple Finite Volume method for Stokes equations

I'm trying to understand how to implement fluid problems using Finite Volume Elements, for example a simple Finite Volume Elements for the Stokes problem: $-\nu\Delta u+\nabla p =f$ in a bounded ...
2
votes
1answer
106 views

How to compute gradient of a cell having a boundary face?

In many situations in unstructured mesh solvers, one needs to compute gradient of arbitrary variable $\phi$ such as temperature or velocity at face centers (one of such situations is correction for ...
2
votes
1answer
181 views

How to calculate skewness for a mesh?

I am writing a code to calculate mesh quality stats such as: cell volume, face areas and non-orthogonality between faces (basically something like OpenFOAM's ...
0
votes
1answer
102 views

Is it possible to predict solution oscillation before solving the system by looking at coefficient matrix?

Question When it is about solving a system of equations, is it possible to predict that whether high-frequency noise (e.g. checker-boarding) is likely to appear in the converged solution by looking at ...
2
votes
3answers
158 views

Flux sign and face normal confusion in finite volume method

I implemented a solver for the 2D steady-state heat equation (without heat generation and homogeneous material) $\nabla. (k\nabla T) = 0$, using finite volume method, however, I am having some ...
1
vote
2answers
156 views

$P0$ elements for $H1$

Are there $P0$ (zero degree/constant element) nonconforming methods for approximating solutions in $H1$? More specifically, I have the equation: $$u-f - T\Delta u = 0$$ Which can be interpreted as ...
1
vote
0answers
37 views

Cell-centered DG extension to the two-point flux approximation scheme

A current problem that I am working on requires me to compute the solution from the heat diffusion evolution on a discontinuous function. More precisely - I have a Delaunay triangulation and within ...
2
votes
0answers
48 views

Solving Stokes Equations in 3D - Do I need to treat pressure-velocity coupling iteratively?

I need to develop a code to solve Stokes Equations in 3D in cubic geometries (structured grid, uniform mesh spacing). My code needs to take a pressure gradient in one direction as a BC (pinlet=p1, ...
1
vote
0answers
35 views

Effect of reducing flux consistency order at boundary on convergence order

Consider the 1D nonstationary convection-diffusion PDE $$ \begin{alignat}{2} \partial_t u &= -a \partial_x u + D \partial_{xx}u, &\qquad x \in (0,1), t \in (0,T), \\ f(t) &= \left.\left( a ...
4
votes
1answer
154 views

Without positive definiteness, does an iterative solver work?

Question Does lacking positive definiteness of the matrix of coefficients in a system of equations, make using iterative solvers impractical? Description Using the finite volume method, I have ...
1
vote
0answers
150 views

Two-dimensional modelling of a flate-plate reactor

I am trying to simulate the unsteady pollutant concentration along the reactor by solving the 2nd order PDE below with the stated BC and IC. which method is appropriate to solve convection-diffusion-...
0
votes
0answers
41 views

Discretisation by FVM with triangulare mesh

I work on electromagnetic modeling, and I use the finite volume method as a tool, so I want to develop code by the modified finite volume method (with triangular mesh), my problem studied is in a ...
1
vote
0answers
24 views

Convection equation expanded in Legendre polynomials

In some areas of physics, for example radiation transport in plasmas, it is beneficial to re-express the convection equation in terms of a Legendre polynomial expansion. This results in a set of ...
0
votes
0answers
74 views

Import stl file into Gmsh and create a 3d mesh for CFD

I am trying to develop a case study for teaching CFD. To this end, I am using OpenVSP, Gmsh and SU2. I have created a wing in OpenVSP, exported the CFD mesh as .stl file. Opened the .stl file in Gmsh ...
11
votes
8answers
6k views

What programming language should I choose and why?

I am a mechanical engineer, intermediated/advanced level in MATLAB and MATHEMATICA, and beginner in Python. I intend to get a PhD in aeroelasticity (FEM + CFD) and coding my own program. I intend to ...
2
votes
0answers
45 views

How does the “Stable Fluids” algorithm by Jos Stam relate to the SIMPLE and PISO algorithms?

The "Stable Fluids" paper (*) by Jos Stam starts by acknowledging that "Our method cannot be found in the computational fluids literature, since it is custom made for computer graphics applications. ...
0
votes
0answers
62 views

The relation between PDE order and discretization order

In Jasak's Ph.D. thesis (2000), a notion is given about discretization of a transport equation: For good accuracy, it is necessary for the order of the discretization to be equal to or higher than the ...
2
votes
3answers
185 views

How is central difference scheme second-order accurate?

In an arbitrarily unstructured mesh, shown in the figure below, in the context of finite volume method, I want to obtain an approximation of $\phi_f$, where $N$ and $P$ are cell centers of adjacent ...
2
votes
1answer
78 views

How do I find the portion of a cell/voxel lying within a defined surface?

We have a 3-dimensional grid of voxels (or cells), with individual voxels being of volume $dx\,dy\,dz$ where $dx=dy=dz=1$. A cone-like surface is defined by some function, $z = f(x, y)$, which in ...
7
votes
2answers
232 views

Numerical flux and source term in FVM (Burger's like equation)

I'm trying to solve the following equation with FVM $$u_t + f(u)_x = g(u)$$ where $g$ is some smooth function of $u$ and $f(u) = \frac{u^2}{2}$. This is really similar to Burger's equation, except ...
0
votes
0answers
52 views

2D diffusion equation using Finite Volume Method

i am working on an assignment problem: Consider a two-dimensional rectangular plate of dimension L = 1 m in the x direction and H = 2 m in the y direction. The plate material has constant thermal ...
3
votes
2answers
145 views

Mass conservation for hyperbolic relaxation problem

I have solved numerically the following system: \begin{cases} \partial_t{u} + \partial_x{v} = 0 \\ \partial_t{v} + \frac{1}{\varepsilon^2}\partial_x{u} = -\frac{1}{\varepsilon^2}(v-f(u)) \end{cases} ...
3
votes
2answers
209 views

Why FVM can handle unstructured meshes while FDM cannot?

How come Finite Volume Method(FVM) handle the unstructured meshes and Finite difference Method cannot, whereas in FVM to approximate the fluxes at the boundary we use the central differencing? My ...
0
votes
0answers
46 views

Pressure boundary conditions in Stokes Equation in 2D (Finite Volumes)

I am solving the steady-state incompressible Stokes equations in 2D: \begin{equation} \frac{\partial u_x}{\partial x} + \frac{\partial u_y}{\partial y} = 0, \end{equation} \begin{equation} \mu\left[\...
2
votes
1answer
182 views

Testing a block tridiagonal system of equations

In 1D problems, tridiagonal systems of equations are obtained when we use finite-difference or finite-volumes in a structured mesh. A wide solver is the TDMA algorithm here. In two-dimensional ...
1
vote
0answers
37 views

How to avoid density getting “deleted” in two way rigid body coupling with LBM CFD?

I've been reading this paper recently, which talks about using Lattice Boltzmann methods and two way coupling. Specifically, it outlines fluid solid coupling, and solid fluid coupling, and how simply ...
0
votes
1answer
39 views

Produce vertex displacements from volumetric shrinkage data on unstructured meshes

I was wondering what would be an efficient way to produce compatible displacements for mesh nodes/vertices if the computed data is volume shrinkage of each element/cell in the unstructured mesh? ...
1
vote
2answers
127 views

FVM vs FDM vs Conservative form vs Non conservative form

My question is regarding solving the conservative form and the non-conservative form of the governing-equations (GE), like continuity or the navier stokes equation, using finite difference method (FDM)...
2
votes
0answers
127 views

CFD and finite volume method: Dirichlet boundary conditions for the Euler equations

Please point me to an answer if one already exists, but after some searching, I still can't find the answer to what seems like a very simple question. There are plenty of references out there for ...
1
vote
1answer
56 views

Obtain velocity from imposed energy spectrum using the inverse FFT

I am trying to obtain the spatial representation of $u(x)$ (e.g. velocity) from its energy spectrum $E(k)=k^4\exp(-(k/k_0)^2)$, which is given in the frequency domain, provided $|u(k)|=\sqrt{2E(k)}$. ...
1
vote
1answer
131 views

How to compute turbulent energy cascade

I need to compute the kinetic energy cascade using a finite volume solution in an equally spaced grid. I wonder if it is more correct to first compute the kinetic energy in the space (or time) domain, ...
0
votes
0answers
37 views

Grids for atmosphere simulation with finite volumes on the globe

I am currently in the early construction process of building a simple CFD model of a rotating planetary atmosphere. The planet should be allowed to tilt significantly, so that a time-dependent source ...
5
votes
1answer
108 views

Approximate $\|\Delta f\|^2_{L^2(\Omega)}$ in finite element context

I have minimization problem of the form $$ G(f) + \|\Delta f\|^2_{L^2(\Omega)} \to \min $$ over all $f\in C^2(\Omega)$, $\Omega$ being closed and bounded. Let us forgot about $G$; I'm interested in ...
1
vote
1answer
105 views

Abaqus, ANSYS, and FVM solver for thermal expansion problem converges to different values

Is it reasonable for a FEM and FVM code to converge to slightly different solutions for the same physical problem (identical BCs, geometry, properties, etc...), provided stability constraints are ...
2
votes
1answer
78 views

Computing face fluxes in FVM

In FVM, we have to compute fluxes at some face of a cell. There are many ways to compute this face flux value, but the most common and easiest way involves some simple averaging at the face. However, ...
1
vote
1answer
235 views

Is there any fundamental difference between meshing for FEM, FVM and FDM?

I am a novice to the field of computational science and have just started studying the FDM and FEM (haven't started on FVM yet). While trying the understand the subject I got this question and trying ...
2
votes
2answers
179 views

Automatic timestep adjustment in a CFD solver

I have developed my own 3D Finite Volume Navier-Stokes solver based on projection method for nonuniform grid. I am looking to incorporate automatic timestep adjustment at each time step based on ...
0
votes
1answer
82 views

Is “Gradient Computation” in Finite Volume Discretization Really 2nd order accurate?

Based on this, pp 245, we go through these steps to discretize a gradient statement, namely $\nabla\phi$: 1- Gauss theorem reads, $$ \int_V\nabla \phi dV = \oint_{\partial V}\phi dS $$ 2- Integral ...
2
votes
1answer
42 views

Discretization with non-constant matrix containg entries form unknown vector

Consider a system of PDEs $$ \begin{cases} u_t = \nabla \cdot (D(u)\nabla u) + \frac{c}{K_U+c}u-ku\\ c_t = d_c\Delta c -\frac{\nu_U c}{K_U + c}u \end{cases} $$ with some boundary conditions. Here, $D(...
0
votes
0answers
124 views

Double mach reflection at a inclined wedge

I am running into a strange problem when solving the 2D compressible Euler equations on a inclined wedge. To elaborate, my top boundary condition seems to emitting some type of instability. I have ...
1
vote
0answers
56 views

Discretizing a parabolic PDE with finite volume method

I want to discretize the following parabolic PDE: $$u_t = \nabla\cdot(\alpha(x)\nabla u)- \beta u\\ x\in\Omega \subset \mathbb{R}^2\\ \partial_n u = 0\\ u(t,0) = u_0(x)\ge 0, \alpha(x)>0$$ Given ...
1
vote
0answers
74 views

Determine truncation error of PDE discretization

The equation is $$\frac{\partial}{\partial x}\left(u\frac{\partial u}{\partial x}\right)=f(x)\\ 0<x<1, u(0)=u(1)=0$$ I'm discretizing this PDE using FVM as follows: $0=x_0=x_{1/2}<x_1<x_{...

1
2 3 4 5