Referring to the discretization of derivatives by Finite differences, and its applications to numerical solutions of partial differential equations.

learn more… | top users | synonyms

2
votes
0answers
22 views

Stability analysis for a hyperbolic PDE on staggered grid

I am trying to understand the stability of a finite difference equation on the staggered grid. I could understand the Von Neumann stability analysis for the collocated grid for a simple acoustic ...
6
votes
1answer
64 views

Best practice for dealing with Dirichlet boundary conditions in finite-difference schemes: add artificial unknowns?

I know of at least two ways of dealing with Dirichlet boundary conditions in finite-difference schemes (and, to a lesser extent, finite-element schemes). Here I'm thinking of solving Poisson's ...
1
vote
0answers
45 views

Can I model laminar incompressible fluid flow and heat transfer in MATLAB's PDE toolbox?

I have a system of PDEs in cylindrical coordinates that needs to be solved: 1. Continuity equation 2. Incompressible Navier stokes ( in r & z coordinates) 3. Heat transfer equation with both ...
-1
votes
0answers
23 views

nonlinear coupled pde by finite difference

I want to solve nonlinear coupled PDE by finite difference method. I have done the discretization. the number of variables and number of equations are not same. How to deal with this situation?
1
vote
0answers
42 views

Oscillations in Chorin's method due to the BC

I am pretty new to the CFD and I wanted to start with Chorin's projection. The starting problem is just a free jet flowing in the investigated area. I got terrible oscillations almost immediately and ...
0
votes
1answer
36 views

How to determine the truncation error with products and quotients

If I have an equation given by $$\displaystyle Y = \frac{a^2}{d^2}\frac{(1-c^2\frac{c}{a})}{(1-b^2)}$$ and I expand $a,b,c,d$ in a Taylor series, where $a$ is truncated at the $A^{th}$ order, $b$ is ...
0
votes
1answer
51 views

Heat Equation in 3D mass Matrix set-up

I am solving a 3D heat transfer equation with variable boundaries (insulated, convective, radiative or free) using a F.D.M. technique. My geometry of choice is a cube. The purpose of my work is to get ...
0
votes
2answers
87 views

Help implementing finite difference scheme for heat equation

I am trying to solve the following problem via a finite difference approximation: $u_t = k \, u_{xx}$, on $0 < x < L$ and $t > 0$; $u(0,t) = u(L,t) = 0$; $u(x,0) = f(x)$. I take $u(x,0) = ...
2
votes
1answer
94 views

finite difference frequency domain Eigenvalue matrix to get eigenmodes

I am using the FDFD method to calculate eigenmodes for an empty wave guide. It's a 1x1meter structure with a PEC boundary. (Here I have 6x6 points to make it simple). Sounds simple, but I can't get it ...
0
votes
0answers
20 views

monitor functions for mesh generation: error estimate by FD or by FEM?

I am using a local truncation error estimate as the monitor function for adaptive mesh refinement that comes from a finite difference(FD) scheme and its values are available only at nodal points. ...
0
votes
2answers
63 views

What is wrong with my code for solving Poisson equation with one side Neumann boundary condition?

I wrote a Matlab code for solving 2D Poisson equation $u_{xx} + u_{yy} + f(x,y) = 0$ on $[a,b]\times [c,d]$ with neumann boundary condition on $x = b$ and the other boundary conditions are dirichlet,...
0
votes
0answers
27 views

Radial wave-function using Finite Differences

I want to solve an equation of the form $$\left[\frac{d^2}{d r^2} + \frac{1}{r}\frac{d}{d r}\right]\psi = 0$$ where $\psi$ is a wave function using finite difference method. The equation is more ...
0
votes
1answer
57 views

Exact finite difference method for advection equation

I want to solve the advection equation: $u_t+au_x=0, a > 0$. Here is our method: $$U_j^{n+1}=U_{j-1}^n-\left(\frac{ak}{h}-1\right)(U_{j-1}^n - U_{j-2}^n)$$ I am trying to answer the following ...
6
votes
2answers
131 views

PetSc vs Sundials for serial numerical computations?

I am currently working on a physics problem that turns into a non-linear boundary value problem. I need an efficient numerical solver that I could run on my laptop with i5 dual core CPU. I am ...
3
votes
0answers
48 views

Corner Transport Upwind for Linear Advection in Arbitrary Velocity Field

I need to implement a 3D version of the Corner Transport Upwind (CTU) finite volume method (in python); and so I've been reading Leveque, "Finite Volume Methods for Hyperbolic Problems" which I think ...
0
votes
0answers
46 views

Multigrid for cell-centered finite difference

Let's say I'm trying to solve a problem of the form: $$ \frac{\partial^2\psi}{\partial x^2} + C \psi = f $$ The problem is discretized with cell-centered finite difference (not something I can ...
1
vote
2answers
32 views

Approximate Neumann BVP operator by a matrix

I am considering approximating the operator of the BVP \begin{cases} -u''+u'=g,&\quad x\in [a,b]\\ u'(a)=-1, u'(b)=1, \end{cases} by a matrix. I tried to use the idea of finite difference ...
0
votes
1answer
55 views

Implementing temperature depending viscosity in a finite-difference scheme

I have a little question that might be basic for some experts, but right now, its not clear for me. I want to implement temperature depending viscosity in a finite difference scheme (incompressible ...
0
votes
1answer
45 views

Discretized matrix from the integral kernel function

Recently, I read a paper [1] and then I want to handle the two-dimensional linear integro-differential equation \begin{equation*} -\triangle u + q\Big(\frac{\partial u}{\partial x} + \frac{\partial u}{...
1
vote
0answers
35 views

Finite differences for incompressible viscous fluid equations

I am working with the equations for incompressible viscous fluid: $$ \partial_t \vec{\omega} + (\vec{u}\cdot\nabla)\vec{\omega} = \nu\nabla^2\vec{\omega} $$ $$ \nabla^2 \vec{\psi} = -\vec{\omega} $$ $...
0
votes
0answers
50 views

Collocated Grid Navier Stokes Solver

I want to solve Navier Stokes equations on a collocated grid. Earlier, I was using a MacCormick scheme based solver where I discretized predictor step in forward differences and corrector step in ...
1
vote
1answer
54 views

How to set the temperature at the vertices points for a rectangular domain?

Suppose I have to solve the 2-D heat equation in a rectangular domain using the finite difference method, for the boundary conditions say: $T_1$ is the temperature of the right side of the rectangle, ...
0
votes
1answer
133 views

Three body problem in C++

I am in a begginers programming course and we got a little project. I chose to simulate the three body problem using the Euler method. Even though the system is chaotic there are some special cases ...
3
votes
1answer
101 views

Second order interpolation scheme

On a grid I am having the values of a physical quantity say for example Temperature, at the E,W,N,S and P node all of them being calculated using a second order discretization scheme. I want a second ...
0
votes
0answers
51 views

Gauss-Seidel iteration weighted by change

In general, I can use Gauss-Seidel iteration for finite difference solution of partial differential equations. In my case I am only solving an analog of steady-state heat transfer, so there is no ...
1
vote
1answer
120 views

How to define residual in multigrid approach?

I wish to solve the two-dimensional Navier Stokes equations using multigrid method on a collocated grid using the Predictor-Corrector method mentioned below. But first, let me elaborate on what I had ...
3
votes
1answer
116 views

Discretize Poisson equation with derivative of delta function as source

Consider the PDE \begin{equation} \frac{d^2}{dx^2} g(x) = \frac{d}{dx} \delta(x-x_0), \end{equation} with $x, x_0 \in [0,1]$ and $g(0)=g(1)=0$. What is the best method to discretise the derivative of ...
1
vote
0answers
117 views

Unwanted Oscillation in FDM simulation of elastic wave equation

I am using staggered grid FDTD for solving elastic wave equation. A description of which can be found at (geodynamics.usc.edu/~becker/teaching/557/reading/Virieux1987.pdf). I have generated a ...
3
votes
1answer
125 views

Is this finite difference approach correct?

I am solving incompressible 2D Navier-Stokes equations with zero y-component velocity. The geometry is a simple 2D pipe of a length $L$ and diameter $W$ and there is only two boundary conditions: ...
3
votes
1answer
100 views

Stability Criterion for this Explicit Scheme

I am solving an unsteady flow using the dual-time Navier-Stokes equation in which I write my momentum equation as: $$\frac{\partial u}{\partial \tau} + \frac{\partial u}{\partial t} + \frac{\partial u^...
1
vote
0answers
72 views

CFD implementation in software [closed]

I've been learning CFD theory for 4 months (which includes FEA, FDM AND FVM) and now I want to start running simulations. As a novice to CFD software, I would appreciate some advice on where to start, ...
0
votes
0answers
139 views

how to solve a 2D non-linear Poisson equation?

I am trying to solve the following equation for $P(x,y)$: $$I = P \nabla \cdot \frac{1}{P^2} \nabla{P}$$ where $P$ and $I$ are functions of x and y. $I(x,y)$ and $P(x,y)$ I want to understand ...
1
vote
1answer
24 views

Suitable method for simulation of in-fiber interferometer

I am trying to simulate an optic-fiber sensor (in-fiber interferometer) to study its respond to temperature. The method I am using is finite-difference time-domain (FDTD), and I come out with a large ...
0
votes
0answers
75 views

Solving Lid-driven Cavity on a Collocated grid

I want to solve the two dimensional lid driven cavity flow on a collocated grid. I already have used a staggered grid and validated my results with Ghia (1982). Recently, I came across this paper (...
0
votes
1answer
63 views

Boundary conditions generalized eigenvalue problem

Consider the following eigenvalue problem \begin{equation} \mathcal {L} x(s) = \lambda x(s), \end{equation} where \begin{equation} \mathcal {L} = \alpha \partial^4_s + (s^2-1)\partial^2_s + s \...
2
votes
1answer
123 views

How do I program periodic boundary conditions? [duplicate]

Hi I have a code below that solves non linear coupled PDE's given Dirichlet boundary conditions. However I need to implement periodic boundary conditions. The periodic boundary conditions are ...
1
vote
0answers
37 views

Representing a 3D system in 2D (Electromagnetic modelling)

Ok so I'm a complete beginner in computational modelling (I use analytical methods of physics typically) but I would like to model an anisotropic, aperiodic (but not random) finite array of metallic ...
-1
votes
1answer
89 views

Stability analysis for coupled nonlinear system of partial differential equations

I'm trying to solve a nonlinear partial differential equation \begin{equation} L(u_{xxtt},u_{xx}u_{tt},u_{xt}^2,u_{xt},u_{tt})=0 \end{equation} using finite difference methods. In order to remove the ...
6
votes
1answer
172 views

What does the Von Neumann's stability analysis tell us about non-linear finite difference equations?

I am reading a paper [1] where they solve the following non-linear equation \begin{equation} u_t + u_x + uu_x - u_{xxt} = 0 \end{equation} using finite difference methods. They also analyse the ...
4
votes
1answer
58 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} - \...
1
vote
1answer
75 views

Finite difference scheme for Webster equation

Webster equation is a popular generalization of 1D wave equation used for ducts of variable cross-section $S \equiv S(x)$. Assuming harmonicity in time, the spatial equation for propagation of ...
5
votes
0answers
75 views

The numerical solution of a (very ugly) set of integro-diferential equations

I've been having fun playing around with a simple model for reaction in a co-flow fluidized bed (I'm afraid I'm a chemical engineer, not a computer scientist). Unfortunately, simple physical ...
7
votes
1answer
132 views

Finite difference recursion and higher order

This may be a trivial question, but I've always wondered... For the classical, central finite difference schemes, if I'm interested in determining the second derivative, does applying the first ...
0
votes
0answers
41 views

Using Periodic BC for Taylor Green Vortex Flow Simulation

I want to solve the very simple case of two dimensional Taylor Green Vortex Flow. I am using incompressible Navier Stokes equations and Artificial Compressibility Method as my solver where I introduce ...
0
votes
0answers
54 views

Periodic boundary conditions for solving Navier Stokes Equations on a Staggered Grid

I want to solve two dimensional Navier Stokes equations on a staggered grid for the case of Taylor-Green Vortex. My initial conditions are standard sine and cosine functions. As I am aware, I should ...
0
votes
0answers
128 views

Heat equation from implicit scheme with Neumann B.C

To solve the heat equation from implicit scheme subject to Neumann boundary condition we can write: $$ T_i^{j+1}-T_i^{j}=\alpha (T_{i+1}^{j+1}-2T_{i}^{j+1}+T_{i-1}^{j+1}) $$ $$ \textbf{A} T^{n+1} = T^...
3
votes
1answer
86 views

Stability in discretization of a PDE

Suppose I want to numerically solve for $f(x,k)$ the one-dimensional Boltzmann equation for electrons in steady-state condition, given as: \begin{equation} \left( \dfrac{\hbar k}{m} \right) \dfrac{\...
0
votes
3answers
314 views

Numerical simulation of a reaction-diffusion system on MATLAB with finite difference discretization of spatial derivative

I have a model of a system which consists of diffusing reactants and intermediates. For each variable, $u_i$, the final representative equation looks like the standard reaction-diffusion form: $$\...
1
vote
0answers
52 views

PDE discretization (via finite difference sheme) question

So after posting this question and reading all your comments I would like to make this new question (update). If you consider the three equations presented here: $$\frac{\partial \rho}{\partial t} +\...
1
vote
1answer
59 views

Non-uniform finite difference Adaptive Mesh Refinement

Assuming that the crosses in the figure below are unknowns in a vertex-centered finite difference scheme in an Adaptive Mesh, how can I calculate the double derivate (Laplacian) at the Red x ? The Red ...