Questions tagged [discretization]

The process of representing a continuum space with a finite set of points/elements

43 questions with no upvoted or accepted answers
Filter by
Sorted by
Tagged with
5
votes
0answers
61 views

Why does the naive barycentric hodgestar fail?

The discrete exterior calculus is defined first using circumcentric dual cells, because the primal and dual edges are orthogonal and thus the dual cells are convex. This leads to a diagonal hodge star ...
5
votes
0answers
548 views

Order of accuracy of FVM discretization

I've recently got interested in CFD and started a small project by solving the radial Reynoldsequation. Why the Reynoldsequation? I recently encountered it through my studies and somehow got stuck :) ...
3
votes
0answers
109 views

How to account for the interface between two different phases in a discretized diffusion model?

I have tried to set up a model for the diffusion of a gas into a liquid. The two media are next to each other and the geometry is spherical because the system should simulate the diffusion out of a ...
3
votes
0answers
50 views

Lax equivalence theorem for integro-differential equation

Can the Lax equivalence theorem (http://en.wikipedia.org/wiki/Lax_equivalence_theorem) be applied to the discretization of integro-differential equations, or does a similar theorem exist for them?
3
votes
1answer
208 views

Which discretization scheme to use for elliptic PDE?

While simulating motion of nonlinear inelastic wire one meets the following equations \begin{align} &{\partial^2\varphi\over\partial t^2}=2{\partial F\over\partial s}{\partial\varphi\over\partial ...
2
votes
0answers
41 views

In sights into why higher order finite differencing method leads faster to instability

I was playing around with numerically solving the 1D wave equation with density and stiffness varying with position using central differencing methods and noticed that for certain discretization steps ...
2
votes
0answers
81 views

Literature on numerical solving based on multiple meshes?

Consider solving a differential equation system for 1D, 2D, or 3D. It involves various input and output "field" variables, which, more often than not, correspond to various physical quantities ...
2
votes
0answers
589 views

Time discretization of wave equation

I am trying to model the seismic wave equation and have therefore been reading about discretization schemes and their stability. I recently came across an insightful paper on 'Galerkin FEM methods for ...
2
votes
0answers
120 views

Finite difference aproximation - Darcy law

I am solving following problem: Filtration of water can be described in bi-dimensional case by $$- \partial_x(K(x,y) \partial_x u ) - \partial_y (K(x,y) \partial_y u ) = 0, $$ where $u$ - water ...
2
votes
0answers
352 views

Finite difference scheme for solving nonlinear least-squares problem

I am dealing with following problem: $$ \min_{u,\gamma}\Bigg\{ \frac{1}{1000} \iint_{S_2} {\gamma (x,y)^2 dxdy} + \iint_{S_2} {[u(x,y) - u_0 (x,y)]^2 dxdy} + \iint_{S_2} {[\Delta u(x,y) - \gamma (...
2
votes
0answers
65 views

2nd order differences for modelling a physically discrete phenomenon?

Consider the diffusion of some unknown $R$ on a physically discrete, 1-D ring: x---------x---------x---------x n = 0 1 2 3 We let ...
2
votes
0answers
185 views

boundary conditions with non-constant coefficients in cell centred finite volume method

Suppose am solving the heat conduction equation in 1d with Dirichlet boundary conditions. The thermal conductivity $k$ is a non constant function. So $-(k(x)u'(x))' = f(x)$ The value of $k$ enters the ...
2
votes
0answers
71 views

What are Implications of Commutative Diagrams?

This question may be too broad. But I really want to know some concrete explanations. I often find various commutations appear here and there, which concerns the application order of two operators. ...
1
vote
0answers
33 views

Procedure to convert continuous equations of motion to discrete version

Let's take a mobile robot and suppose we know its continuous equations of motion, for example this car-like simple model. Now if I am simulating this robot in a continuous 2D plane, then coding the ...
1
vote
0answers
39 views

discretization of advection diffusion with variable coefficients

I am looking for help to find a somewhat stable FD numerical scheme for the advection diffusion equation posed on a curve $(x(r),y(r))$. The equation becomes $$u_t=\alpha(r) u_r +\beta(r) u_{rr}+f(r,t)...
1
vote
0answers
60 views

Can the standard multigrid performance be used for time-dependent PDEs?

Consider a time dependent pde(i.e u(x,t)).I know when only space-coarsening is used the standard multigrid performance can be applied but what if instead we use only time-coarsening?Can we apply the ...
1
vote
0answers
108 views

FDM discretization of equation on the boundary

In order to simulate the following equation using FDM $$u_t(t,x)-u_{xx}(t,x)=0, \quad (t,x) \in (0,1)\times (0,1)$$ $$(u_t(t,x)-u_{x}(t,x))\rvert_{x=0}=0, \quad t \in (0,1)$$ $$(u_t(t,x)+u_{x}(t,x))\...
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
81 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
vote
0answers
114 views

PDE discretization on triangular domain

Given the 2D Poisson equation $$\Delta u = f\\ u(x,0) = g_1(x), 0<x<1\\u(0,y) = g_2(y), 0<y<1\\ \partial_n u (x, 1-x) =0, 0<x<1$$ defined on the domain $\Omega := \{(x,y) \in \...
1
vote
0answers
52 views

Decrease in slope during convergence analysis

I am using the method of manufactured solutions to perform the order of accuracy testing. I am using a cube for the testing. The cube is size 1m on all sides. I used 5 refinements: $dx = dy = dz = ...
1
vote
0answers
41 views

Limit to volume change in a discretized mathematical model?

I have set up a mathematical model describing the diffusion of ozone out of a gas bubble. The bubble is surrounded by a thin gas film. So actually, the model describes the diffusion of ozone through ...
1
vote
0answers
331 views

Global truncation error behavior at fixed time step

I am trying to solve the following diffusion equation problem: $\frac{\partial f}{\partial t}=\frac{\partial (D\frac{\partial f}{\partial x})}{\partial x}+S$ $D=1+x^{2}+\sin(x)$ $f(x,0)=1 , f(0,t)...
1
vote
0answers
385 views

Implement Robin boundary condition (finite volume)

I have a PDE equation with Robin Boundary condition in an annulus system and I should solve it by finite volume method: \begin{align} \frac{\partial T_f}{\partial t} - k \left(\frac{\...
1
vote
0answers
82 views

Discrete operator textbooks

This will be a vague question. When I was writing a finite element matrix assembly routine, a colleague noticed that I had a bug in my code because the sparsity pattern of the one of the blocks didn'...
1
vote
0answers
447 views

Rhie and Chow Pressure Velocity Coupling

In a collocated grid, which velocity is used in the convective term in the momentum equation? Is the Rhie and Chow constructed face velocity or an average of the adjacent cell center values(in a ...
1
vote
0answers
79 views

Three steps of pde numerical solution and nonlinear equation

I'm very new here. I'm trying to solve nonlinear elliptic equation $$ (n(u)u')' = f(u) $$ and face with crucial misunderstanding. As I suppose, the procedure of solving some nonlinear equation ...
1
vote
0answers
37 views

Inner Products vs Discretizations of Functions

Let $f:(0,1)\rightarrow [0,1]$ be a $\mathcal{C}^{\infty}$ function. We say that a function $D_r^f:\{1,...,r\}\rightarrow [0,1]$ is an $r$-discretization of $f$ if $$D_r^f(j) = \frac{1}{|a(j)-b(j)|}\...
1
vote
0answers
55 views

Second and Higher Order Order Corrector in Spectral Deferred Correction

I am trying to work out a second order or higher order correction for the method of Spectral deferred Correction (SDC). Specifically using as a corrector a second order or third order multi-step. In ...
1
vote
0answers
70 views

Does scaling factor affect discretization?

Suppose I want to solve the below equation numerically. $$ \frac{dy}{dx}=y $$ I'd like to normalize the space discretization by choosing $$ a\bar{x}=x $$ where I assume $\bar{x}$ is unity. Then the ...
1
vote
0answers
41 views

Stability in Discretization of 1D Stationary Boltzmann equation

I want to discretize and numerically solve the following PDE: \begin{equation} v(k)\dfrac{\partial f}{\partial x} + E(x)\dfrac{\partial f}{\partial k} = S\{f\} \end{equation} using finite volume (box ...
1
vote
0answers
189 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
0answers
92 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 ...
1
vote
0answers
54 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 \end{array}\...
1
vote
0answers
105 views

How to compute this double integral?

Let $$T=1, K=100, S_0=100, \sigma=0.05, r=0.15. $$ Define $\nu:=\frac{2r}{\sigma^2}-1$ and $$H(y,z)=\frac{z e^{\pi^2 /4y}}{\pi \sqrt{\pi y}}\int_0^{\infty} e^{-z \cosh(u) -u^2/(4y)} \sinh(u) \sin(...
1
vote
0answers
554 views

Equilateral triangle based mesh generation by intersection

In work I am currently working on I need to mesh some structure with equilateral triangles to study it using a kind of discrete element method known as spring networks or Lattice model. To mesh the ...
0
votes
0answers
47 views

Discretizing Multi-species Ion Exchange Equations by Finite Volume Method

I'm solving a system of multispecies ion exchange equations (diffusion+drift fluxes) in 1-d spherical domain using finite volume method to obtain the ion concentrations at the next time step. After ...
0
votes
0answers
51 views

What is the rationale of second-order finite volume discretization?

When it comes to a second-order accurate finite volume discretization of Navier-Stokes equations, which one of the two following rationales is adopted? 1- Second-order accuracy is a direct consequence ...
0
votes
0answers
44 views

How to solve for underlying function from discrete data set containing integral of that function

New to Computational Science, I hope I'm on the right exchange network for this question. I have a time series data set that contains the sum of a source data set representing an exponential decay ...
0
votes
0answers
65 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 ...
0
votes
0answers
75 views

Analytic vs discrete understanding of PDE

The PDE I am working with: $$\partial_tu = \nabla \cdot (a(x)\nabla u)-\beta(x)u\\ \partial_nu=0, x \in \Omega \subset \mathbb{R}^2\\ \beta(x)>0$$ Integrate the PDE: $$\int_\Omega \partial_t u=\...
0
votes
0answers
53 views

Existence and uniquness of solution of FVM for Poisson equation

I'm discretizing the following Poisson equation using FVM where the domain $\Omega$ of the solution is a regular hexagon of side $1$ centered about the origin. $$\Delta u =k,\text{ $k$ constant}\\ \...
0
votes
0answers
47 views

Order of accuracy of linearised vs non-linear system

Does the order of accuracy of a combination of schemes applied to solve a system of non-linear equations, match those of the same schemes applied to the linearised version of the system? In other ...