Questions tagged [finite-volume]

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

Filter by
Sorted by
Tagged with
5
votes
1answer
95 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 ...
0
votes
0answers
8 views

Universal formulation of adiabatic equations of state in compresible finite-volume simultions

I code some finite element solver which should work for broad variety of materials (i.e. gas, liquid, solid, plasma) and large span of compressions resp. densities. I want to simulate things like ...
1
vote
1answer
60 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
52 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
64 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
143 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
64 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
36 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
75 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
52 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
60 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_{...
0
votes
1answer
132 views

How to make a less diffusive code to solve 2D advection equation?

I would like to solve the following differential equation numerically in 2D, $$\frac{\partial z^-}{\partial t}+(\vec{B}\cdot\vec{\nabla})z^-=0,$$ see Wikipedia if you are curious about what the ...
0
votes
1answer
88 views

Elliptic PDE finite volume method with Dirichlet boundary condition

I want to discretize the following equation using a Finite Volume Method $$\nabla \cdot (a(x)\nabla u)=f(x)\\x\in \Omega \subset \mathbb{R}^2 \\u_{|\partial\Omega}=g$$ I'm using Voronoi cells here: ...
6
votes
3answers
245 views

Why is the FVM traditionally used in CFD, and FEM in computational structures?

Most CFD codes use FVM. Most computational structures codes use FEM. Why is the FEM not frequently used in CFD, and why is FVM not frequently used in FEM?
3
votes
2answers
112 views

Finite volume discretization of non-conservative linear hyperbolic equation

Problem. Consider the one-dimensional adjoint Euler equations for $(x,t) \in \Omega \times [0,T]$ with $\Omega \subset \mathbb{R}$ and $T > 0$ $$ \varphi_t + \Big(\frac{\mathrm{d}F}{\mathrm{d} U}(x)...
1
vote
1answer
102 views

How to generate a face list from vertices?

I have a little background in writing toy finite volume CFD codes. In 2D Cartesian scenarios, I typically take $x_{\min}$, $x_{\max}$, $y_{\min}$, $y_{\max}$, and the number of points in $x$ and $y$ ...
1
vote
2answers
66 views

TVD for temporal dicretisation

I have come across schemes where TVD (with flux limiters) is used for spatial discretisation along with Runge-kutta for Temporal discretisation. Can TVD be used for Temporal discretisation? If so ...
1
vote
1answer
53 views

Strain from FEM simulations to strain gauge measurements

I am looking for some intuition in making comparisons of FEM simulations to experimental measurements. In particular, I am interested in comparisons to strain gauge readings, and perhaps even LVDTs. ...
1
vote
1answer
82 views

Grid dependence of a numerical model

Statement of the problem Suppose, we consider the following model $$ \begin{array}{l} (1)~\mathbf{u}_t + \mathbf{F(u)}_x = \mathbf{S}(\mathbf{u},\mathbf{w}), \\ (2)~\mathbf{w}_x = \mathbf{P}(\mathbf{...
1
vote
0answers
47 views

Thermal stress to displacement through finite volume method

I am attempting to solve for a displacement field where I know the thermal stresses on my discretized domain, which consists of hex cells. My first question is: is easier way to solve for the ...
1
vote
0answers
50 views

Roe Riemann solver for perfect gas mixture

I have working program for solving one-component 1D Euler equations with Roe's approximate Riemann solver constructed according to this pdf. My implementation of the algorithm is as follows ($\rho$ is ...
0
votes
1answer
74 views

Exact Riemann solver for perfect gas mixture: problem with Newton's method convergence

I'm trying to solve multicomponent Euler equations for perfect gas mixture with Godunov-like scheme using exact Riemann solver. Of course, some approximate solver would probably be more cost-effective,...
2
votes
2answers
83 views

Elliptic equation with finite volume and unstructured high order geometry

I have found that in unstructured mesh, discretizing the laplacian operator with finite volumes requires special care, as given in An Introduction to Computational Fluid Dynamics: The Finite Volume ...
2
votes
0answers
55 views

WENO methods: why the characteristic wise method resulting big errors?

I was doing my research/project using WENO as the limiter in finite volume methods to solve hyperbolic conservation law. I have no idea why the result in the characteristic wise method has a big error ...
1
vote
1answer
78 views

TDMA with 3rd order upwind scheme

I'm trying to implement a model I found in a paper, but there is something I do not understand. The authors say they use TDMA to solve their equations; however, they use a 3rd order upwind biased ...
2
votes
0answers
62 views

Harmonic average of Diffusion Tensors in Finite Volume Method

I want to implement a Bilinear Finite Volume discretisation of the anisotropic diffusion problem: $$\frac{du}{dt} = \nabla \cdot (\textbf{D} \nabla u)$$ Both my degrees of freedom as well as the ...
1
vote
0answers
82 views

Basic approach for numerical solution of PDE

I'm looking for some guidance on how to write a program to numerically solve a PDE. As an example for comparison in 1D: $$\frac{d^2u}{dx^2} = f\;\;\;\;u(0) = 0\;\;\;\;u(1) = 0$$ We could try 6 ...
0
votes
1answer
24 views

Calculating volume of a discretised diffuse interface object

Suppose I have a spherical object projected onto a discrete square mesh. The dicretised circle can be represented by filling a logical matrix such that voxels in the interior of the sphere are filled ...
0
votes
1answer
39 views

Solving the diffusion/heat equation for a randomly distributed set of points in 3D

In this problem I am trying to solve, I have a messy set of points distributed in 3D space, each with a defined temperature. If I would want to calculate the heat transfer scenario in this system, how ...
2
votes
1answer
86 views

Does the PDE hold at every cell in a FVM mesh?

If you solve a given PDE (Navier stoke's, Euler, heat eqn, advection eqn, etc...) using FVM, is this PDE supposed to be valid at every cell in the discretized domain, or only in the global domain as a ...
1
vote
0answers
48 views

Jacobian Elements for Coupled Drift-Diffusion System using Vertex-Centered Finite Volume

I'm trying to solve the fully coupled drift-diffusion system using Newton's Method. Although I eventually plan to potentially use a Jacobian-Free Newton-Krylov approach, this is still something that I ...
3
votes
0answers
140 views

Finite Element Stabilization for Drift-Diffusion/Advection-Diffusion Equations

I've tried my best to look through the relevant suggested similar questions when posting this, and hopefully this contains enough new material to not be considered a duplicate. I'm currently trying ...
3
votes
1answer
140 views

Numerical Lax-Wendroff scheme order of convergence on Burgers equation

I was suggested to move that question here. The question to be as follows. Statement of the problem Is it possible to achieve the second order of convergence (OOC) of Lax-Wendroff (LxW) scheme ...
1
vote
1answer
50 views

Problems with deriving an equation for a finite-difference scheme given in the journal paper

I'm reading this paper and trying to follow everything that the author has done. A Comparative Study of Finite Volume Method and Finite Difference Method for Convection-Diffusion Problem But there ...
2
votes
0answers
88 views

Monotone, monotonicity preserving, LED, TVD, NVD, bounded, stable and stability preserving discretization schemes [closed]

When it comes to discretization schemes for finite volume method, the following terms can be found in literature: monotone schemes monotonicity preserving schemes local extremum diminishing schemes ...
1
vote
0answers
42 views

Convergence rate assesment in space/time

I am solving a hyperbolic PDE (e.g. the shallow water equations) which depends upon $x$ and $t$. Typically, the overall convergence rate is calculated by comparing the numerical error in different ...
1
vote
5answers
388 views

Don't we care about the numerical diffusion in the diffusion term?

In the context of the solution of advection-diffusion equations by finite volume method, many numerical schemes, papers and book chapters are dedicated to address the numerical diffusion and/or ...
1
vote
0answers
56 views

Discretization used for steady linear elasticity - followup to previous question

This is a follow-up of my previous scicomp question (https://scicomp.stackexchange.com/posts/28863/edit). I figured I'd start a new thread on this as the question is a bit different from my previous ...
4
votes
1answer
92 views

Finite-volume method applied to a particular advection equation

I'm trying to apply the finite-volume method (FVM), with which I'm not so familiar, so a simple 1D PDE equation. I already asked this question in the math stackexchange, but was told that it could be ...
0
votes
1answer
338 views

Finite volume software packages

There are many software packages for the finite element method, of which the most popular are listed e.g. on Wikipedia. When it comes to the finite volume method, I'm not aware of any similarly ...
1
vote
0answers
158 views

Odd-even decoupling at faces of cells

I am currently solving PDEs using the finite volume method. The surface integrals of the equations that I am solving involves computing face gradients. The current algorithm that we use to compute ...
1
vote
0answers
109 views

FiPy with derivative source terms

I have a coupled nonlinear PDE system in 1 spatial dimension, in which I want to solve using FiPy. The dependent variables are $n$ and $T$: \begin{align} \frac{\partial n}{\partial t} \,&=\, D\,\...
0
votes
1answer
285 views

Crank-Nicolson method and mixed derivatives

I am curious if anyone had literature references or knowledge on how to apply the Crank-Nicolson (with approximate factorization) to the $$ \nabla \cdot (\nu (\nabla \mathbf{u} + \left(\nabla \mathbf{...
2
votes
2answers
301 views

Finite Difference and Finite Volume as special cases of Finite Element

I have heard this many times that "FD is a special case of FE" and separately, that "FV is a special case of FE". However I have never come across a brief document that substantiates that claim, ...
1
vote
0answers
75 views

Choice of velocity grid - staggered or not?

I'm trying to understand when and why one would use a staggered vs. a colocated grid in problems that have velocities and scalars that they transport (e.e. density). If scalars are defined cell-...
0
votes
0answers
82 views

Reconstruction of cells information when given mesh in polyMesh format [duplicate]

The polyMesh format used in OpenFOAM is an intelligent format for which 'cells' file, i.e. the file giving list of nodes that make a cell, is not needed. The way to reconstruct geometry data such as ...
6
votes
1answer
216 views

Why is the continuous Galerkin Finite Element Method a poor choice for the inverse problem for the Navier-Lame equation?

I work in a field (elasticity reconstruction) that frequently uses standard Finite Element methods to solve the first order PDE, where we are given u and solve for mu and lambda: $$ (\mu(u_{i,j} + u_{...
2
votes
2answers
198 views

Computing Roe's average density for General Equation of State

I am solving a 1d Shock tube problem of compressible fluid obeying Euler equations(Hyperbolic pde). I am trying to simulate it using Finite Volume Method using Roe's Scheme. Half of the tube contains ...
7
votes
0answers
211 views

What can be done with Finite Element Method and not with the Finite Volume Method, and vice versa?

What are some applications where you would absolutely go for either FEM, but not FVM, or vice versa? What are some applications where both methods are equally suited? I worked with the FEM so far and ...
0
votes
2answers
125 views

If a numerical solution remains constant for different grid sizes, what does it mean?

I'm testing a finite volume scheme, Godunov type solver, using a problem with analytical solution. I'm not able to reproduce the solution, which includes source terms. I tried with different mesh ...