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

learn more… | top users | synonyms

2
votes
0answers
52 views

Increasing the accuracy of numerical discretization

In order to numerically solve the following differential equation: \begin{equation} \text{Fr}\{f\} := v(k)\dfrac{\partial f(z,k)}{\partial z} - F(z) \dfrac{\partial f(z,k)}{\partial k} = ...
0
votes
0answers
79 views
+50

Concave Cells and Finite Volume Method

I would like to know how to compute the "best" cell center for pressure calculation in case of concave polyhedral cells, if it makes any sense. My face centers are defined as faces barycenters and I'm ...
3
votes
1answer
46 views

Is it necessary to do characteristic decomposition to solve 2D Euler Equation?

Recently, I am trying to solve 2D Euler Equation using FV-WENO reconstruction. For the smooth initial problem, component-wise WENO reconstruction works well. My problem is, if I want to solve double ...
1
vote
1answer
94 views

Outflow boundary condition

I know that in outflow boundary we assume a zero normal gradient condition and use upwind scheme for approximation. However, I saw this sentence in a book which I do not understand; "Convective fluxes ...
1
vote
1answer
46 views

Is OPENFOAM a good idea for direction/operator splitting in parabolic equation

I would like to know if I should switch to OPENFOAM for my task. I work only with Cartesian grids, right now in 2d, rectangular domains only. If $\mathbf{w} = (w_1, w_2)$ and suppose I want to solve ...
4
votes
1answer
55 views

Boltzmann equation and equilibrium distribution function

The free-streaming term in stationary 1D Boltzmann equation satisfies the equilibrium distribution function, that is: \begin{equation} \text{Fr}\{f\} := v(k)\dfrac{\partial f(x,k)}{\partial x} - ...
3
votes
0answers
44 views

Reflecting boundary condition posed as a Riemann problem

I am trying to implement a solver for the Euler/Navier Stokes equations. I have a problem implementing boundary conditions for the wall. I am using an unstructured solver. A lot of literature says ...
3
votes
1answer
101 views

Is the Finite Volume Method conservative when the source term depends on the variable?

I trying to do a simulation where there are two quantities, $\Delta\theta(x,t) = \theta(x,t) - \theta_{o}$ and $\Delta\nu(x,t) = \nu(x,t) - \nu_{o}$. These quantities are chemical concentrations ...
3
votes
2answers
139 views

Are the Finite Difference and Finite Volume Methods different after the application of the Gauss Divergence theorem on the FVM?

I have a rudimentary question about the differences between finite difference (FDM) and finite volume methods (FVM). In FDM we concentrate on the nodes (points) in space while in FVM we concentrate ...
1
vote
1answer
107 views

1D uniform flow test case for compressible flow

Are there any one dimensional hyperbolic test cases for a 1D uniform flow? (I need the test case for the purpose of showing the satisfaction of geometric conservation law on a moving mesh).
3
votes
1answer
107 views

a few questions on understanding geometric conservation law

I was wondering if anyone could help with understanding the geometric conservation law for moving domains. I came across Link1, and have tried to understood the paper by Farhat et al Link2. So far ...
0
votes
1answer
67 views

Order of accuracy of DGFEM or FEM

I know that it is possible to determine the theoretical order of accuracy (B) of numerical solutions in FVM (for instance for a steady problem in which only the central differencing scheme (CDS) is ...
1
vote
1answer
136 views

what is the difference between non-conformal and conformal?

So far what I understand is that two neighbouring elements are conformal if their edges and faces match exactly, whereas with non-conformal elements this is not the case. For instance, h-refinement ...
0
votes
1answer
77 views

Is it possible to show global conservative properties FEM as it is done in FVM?

I know that in FVM, it is possible to show that a discretisation scheme is conservative by adding the discrete terms over a few control volumes and showing that all terms cancel apart from those ...
0
votes
0answers
51 views

Do collocated grid arrangements definitely result in the checkerboard effect?

I understand the checkerboard effect due to the use of collocated grid arrangements in FVM. However, I wanted to know whether this problem is definitely bound to effect the results? For instance, I ...
1
vote
1answer
101 views

Are FEM or DGFEM methods based on integrals or PDEs?

I know that FVM is based on the integral form of conservation laws, and FDS is based on PDEs. What I'm confused by, is whether FEM and DGFEM formulations are based on integral or pde form of ...
0
votes
0answers
53 views

problem with understanding the fluid boundary conditions of a 1D probelm

I am having problems understanding the boundary conditions of the problem described in this paper on researchgate Essentially the problem consists of a one dimensional fluid chamber in contact with a ...
0
votes
0answers
54 views

difference/similarities between DGFEM and FVM and FEM methods

I have been to trying to understand the discontinous galerkin method. My understanding is that in comparison to FVM, the essential difference is in the weighting function used and by setting the ...
2
votes
1answer
41 views

Use of structs in Axisymmetric Finite Volume method

This might be better somewhere else, but I'll give it a try here first. I'm implementing a finite volume scheme for an axisymmetric problem in C, and am looking for a more efficient way to handle all ...
1
vote
1answer
103 views

Diffusion with space dependent drift in Fipy

I need to solve a diffusion equation in periodic boundary conditions using fipy but I would like to have a drift term that depends on the position so like this: $$ \partial_t u(x,t) = \partial_x^2 ...
2
votes
1answer
137 views

MAC Projection in Projection method?

My question concerns the following paper: A Conservative Adaptive Projection Method for the Variable Density Incompressible Navier–Stokes Equations ...
0
votes
0answers
41 views

Restarting the Integrator

I'm solving a diffusion-reaction PDE with discontinuous variable coefficient and source term (like a step function) at a point. I use FVM with harmonic average on the coefficient, but no special ...
1
vote
2answers
193 views

Dirichlet BCs - alternative implementation methods

I am having problems with solving a hyperbolic wave problem with Dirichlet BCs. I have tried reducing the time step sizes, which does not affect the results, and notices increasing the number of nodes ...
1
vote
0answers
45 views

should boundary conditions be effecting moving mesh results?

I have a question on the use of moving mesh to solve the inviscid euler equations. I have solved the following equations: $$\frac{\partial}{\partial t}\left[\begin{array}{c} \rho\\ \rho u ...
0
votes
1answer
86 views

How to define fluxes for two dimensional convection-diffusion equation?

I want to solve the following differential equation using control volume approach on a Cartesian mesh: $$\frac{\partial T}{\partial t} + \frac{\partial T}{\partial x} + \frac{\partial T}{\partial y}= ...
1
vote
1answer
107 views

FVM - virtual node discretisation

I have come across the paper titled: A monolithic fluid structure interaction algorithm ... For a 1D grid, at the boundaries the paper uses virtual nodes $x_{0}$ and $x_{N+1}$ (page 372) and for ...
1
vote
0answers
76 views

Discontinuity at Interface

The equation at the left of the interface is \begin{equation} \displaystyle\frac{\partial C_i}{\partial t} = D_i \nabla^2 C_i - z_i \frac{D_i}{RT}F \nabla \cdot (C_i \nabla \phi_2) \end{equation} ...
3
votes
1answer
70 views

linear stability analysis using spectral radius

I am analysing the stability of a series of 1D linear equations of the form \begin{equation} \frac{d}{dt} x = A x \end{equation} discretised using upwind and central finite volume methods, etc, with ...
5
votes
1answer
139 views

Flux at coarse-fine mesh grid interface?

I am trying to solve one dimensional inviscid Burger's equation using adaptive mesh refinement. This is the PDE: $$\frac{\delta U}{\delta t} + \frac{\delta F}{\delta x} = 0$$ where the flux F of the ...
3
votes
0answers
82 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 ...
4
votes
1answer
93 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 ...
2
votes
2answers
424 views

The meaning of conservative discretization in Galerkin FEM and Discontinuous Galerkin

I do understand the meanning of "conservative discretization" within the FVM/FDM framework, indeed it is well explained in this post. Now, according to the table in this slide (pp.4), it concludes: ...
4
votes
1answer
54 views

Estimating the local compression/expansion ratio for a transformation on a point cloud

Let's say we have an unorganized point cloud P1 with N points, each with coordinates {x,y,z}. We apply non-rigid transformation to P1 (translation + rotation + warping), to obtain point cloud P2. ...
2
votes
1answer
222 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
178 views

What is the difference between SSPRK3 and RK3 time discretization methods?

Here, SSPRK3 refers to third order strong stability preserving Runge-Kutta and RK3 refres to regular third order Runge-Kutta method. The meaning of the method is obvious from the name. However there ...
0
votes
1answer
189 views

Temperature dependent 1-d conduction in Python?

I'm trying to write a Python code that is a numerical solver for 1-d heat conduction (using FVM) with a temperature dependent thermal conductivity. The solver has three functions I need to iterate ...
1
vote
1answer
110 views

What does “strongly conservative” mean in the context of numerical methods?

I have a homework problem that asks me to show that 1st order unwinding or central differencing can give a strongly conservative, consistent scheme for the 1-D Burger's Equation using a finite volume ...
4
votes
2answers
366 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 ...
2
votes
0answers
107 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 ...
0
votes
0answers
33 views

How to implement two-point viscous fluxes in FVM context

I am trying to implement the two-point viscous fluxes in an in-house unstructured code (FV, all hex cells). It will be really helpful if anyone can explain the algorithm for non-orthogonal meshes. ...
1
vote
1answer
66 views

Application of CLAWPACK to Richards' equation

I'm looking to solve the Richards' equation. This models water flow in porous media and is a nonlinear, possibly degenerative, parabolic differential equation that takes the form $\partial_t ...
2
votes
1answer
109 views

Unable to validate the Roe matrix for the Shallow Water Equations

In LeVeque's Finite Volume Methods for Hyperbolic Problems, p. 320-321, one may find the derivation of the Roe matrix to the 1D Shallow Water Equations (SWEs). It is $$ ...
2
votes
0answers
163 views

Finite Volume/difference scheme for 2d continuity equation

I am trying to solve the following 2D continuity equation on the rectangular domain [0,1]x[0,1] using the finite volume method. $\rho_t + \nabla \cdot (\rho v) = 0$ where the velocity $v =(v1,v2)$ is ...
2
votes
1answer
172 views

Finite volume method implementation issues

I am trying to write a simple finite volume method code but there are some concepts I'm still not really getting right (perhaps I'm overcomplicating things) Given a uniform grid, the idea is to ...
1
vote
1answer
82 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} ...
1
vote
1answer
191 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( ...
2
votes
0answers
131 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 ...
2
votes
1answer
232 views

Energy Conservation in Conservation Laws with Source Terms

I'm wondering if anyone can help me understand energy conservation when using conservation law methods (i.e. Riemann solver, High-Resolution Wave-Propagation Methods) with the addition of source ...
3
votes
1answer
444 views

Numerically solve a PDE in Python with a term calculated by coarse-graining

I'm trying to solve a PDE in Python of the form, $\dfrac{\partial c(\mathbf{x}, t)}{\partial t} = \mathrm{D} \nabla^2 c(\mathbf{x}, t) -\gamma \rho(\mathbf{x}, t) c(\mathbf{x}, t)$ where $c$ ...
0
votes
1answer
509 views

a circular plot from a vector which represents the temperature along the radius surface, which is the same for every radius

I have calculated the temperature of the section of a cylinder, which is subjected to a heat flow on its upper surface. Getting the temperature distribution in the 2D section. As shown in the ...