Finite difference Definition and 144 Threads

I thought I solved the problem in answering my own post a few days ago, but the tunneling probability vs. energy trend is clearly wrong. I've remade the post because I have totally changed my approach and need a better understanding of the boundary setup. Overall description: a plane wave...

I recently made a Python library for modelling (very basic) finite difference problems. The Github readme goes into details of what it does and how it works, and I put together a Google Colab with some examples (diffusion, advection, water wave refraction) with interactive visuals. I'd love to...

I am trying to apply finite difference scheme for Beam propagation method by following this paper. I was wondering if anyone can share their code if they have implemented this method. I can share my code which is not working as expected and can get some insights if possible.

It seems no finite difference scheme is stable for a>0, dt>0, correct?

u'(t)=u(t)^2, u(0)=1, 0<t<2 What finite difference scheme can overcome the difficulty when t->1+ and help the solution jump to the negative branch?

Let's discuss whether the energy under a finite difference (FD) scheme is conserved. Take the simplest vibration eq mx''+kx=0, which one will use a FD scheme to solve. The energy is mx'^2/2+kx^2/2. Whether the energy is conserved doesn't depend on the FD scheme for the ODE but upon the FD scheme...

Homework Statement:: Discuss the limitation of the Explicit Finite Difference Model. Relevant Equations:: no formula Hello there, I have to discuss the limitations of using the Explicit Finite Difference model to calculate a 2D Heat Diffusion through an aluminium place, however, I really...

Here is the paper again: https://www.mdpi.com/2218-2004/6/2/22?type=check_update&version=2#related_content For a class project I need to calculate the energy levels of atoms using the Hartree Fock method as presented in this paper which essentially brute forces the calculation using finite...

So for my scheme I obtained ##\frac{\mu}{h^2} U_{p}+(\frac{v_{1}}{2 h}-\frac{\mu}{h^2})U_{E}+(\frac{v_{2}}{2 h} - \frac{\mu}{h^2})U_{N} - (\frac{v_{1}}{2 h}+\frac{\mu}{h^2})U_{W} - (\frac{v_{2}}{2 h} + \frac{\mu}{h^2})U_{N} + \tau = f## however I am not sure this is correct. I am quite new to...

As part of my project I was asked to use the finite difference method to solve Schrodinger equation. I see how you can turn it into a matrix equation, but I don't know how to solve it if the energy eigenvalues are unknown. Are there any recommended methods I can use to determine those...

I have lately been working with Numerical Analysis and I am using Finite Difference Methods for Ordinary and Partial Differential Equations by Randall J. LeVeque. It was recommended to me by a friend of mine (physicist) https://epubs.siam.org/doi/book/10.1137/1.9780898717839?mobileUi=0&...

Homework Statement Homework Equations Finite difference method The Attempt at a Solution I have tried two different approaches, but still i am wrong in the question. Can anyone guide me how to attempt this question? Thank you

I am working on modelling of a heat storage tank. More specifically, I need to find out the transient temperature variation through the tank height. My question is about the plume entrainment considering the case that hot water is injected from the bottom part of a heat storage tank. This...

Hi. I have been trying to solve this problem that has been keeping me up at night for a coupe weeks at least. If anyone can help me, I would be greatly appreciated. Hot air enters a cylindrical duct. The duct has some R-value and radiation and convection is being accounted for on the outside...

Hi PF, initially I would like you to focus on that link https://books.google.com.tr/books?id=Dkp6CwAAQBAJ&pg=PA389&lpg=PA389&dq=runge+kutta+method++is+tvd+proof&source=bl&ots=47ULQDVwcC&sig=e2zjdnXENJ7WxBbrf6hXkSouvLI&hl=tr&sa=X&ved=0ahUKEwjU5Z2XsbXZAhUMCMAKHWpnATQ4ChDoAQhKMAQ#v=onepage&q=runge...

Homework Statement Determine the Finite Difference Method stencil for approximating a second derivative u''(x) at a discrete set of nodes with maximum accuracy for stencil of sizes (0,4) (off-centered). My questions: I think I am able to answer the question I am just not sure about what is...

For example, when we solve simple 1D Poisson equation by finite difference method, why three point central difference scheme on uniform grid (attached image) is second order method for solution convergence? I understand why approximation of first derivative is second order (and that second...

Hi PF! I'm trying to finite difference a FTCS PDE $$h_t+\partial_x\left[h^3(h_{xxx}+(Kf'(h)-G)h_x\right] = 0$$ where ##K## and ##G## are constants. ##f'(h) = -(n(h^*/h)^n-m(h^*/h)^m)/h##. Boundary conditions are ##h_x=h_{xxx}=0## at both ends of the ##x## domain (however long you want to make...

Given a function ##\psi## having ##N## components, how would you (as fast as possible) construct ##\psi'## also having ##N## components? I thought about taking a forward finite difference approach on the first ##N-1## components of ##\psi## to generate ##\psi'## and then a backward finite...

What is a finite difference discretization for the fourth-order partial differential terms \frac{\partial u}{\partial x}k\frac{\partial u}{\partial x}\frac{\partial u}{\partial x}k(x,y)\frac{\partial u}{\partial x} and \frac{\partial u}{\partial x}k(x,y) \frac{\partial u}{\partial y}...

Hi. I am trying to simulate this paper since apparently I have a lot of time. Scrolling down to the last page, he simulated a transient 2D heat conduction plate with composite slabs on it. Darkest one is copper, lighter one is steel, lightest one is glass. If you look closely, the authors said...

So I have this PDE: d2T/dr2 + 1/r dT/dr + d2T/dθ2 = 0. How do I implement dT/dr || [r = 0] = 0? Also, what should I do about 1/r? This is actually the first time I am going to attack FDF in polar/cylindrical coordinates. I can finite-difference the base equation fairly decently; I am just...

I have completed a 2D finite difference code in MATLAB that has a domain of (0,1)x(0,1) and has Dirichlet Boundary Conditions of value zero along the boundary. I get convergence rates of 2 for second order and 4 for fourth order. My issue now is that I'm now wanting to change the domain to a...

Hi, Physics forum! Just a little push of my doubts I hope somebody could help me with my confusion of one of our home works. I know that all boundary conditions are zero. My doubt is how do I interpret (x,y,0)=0.01 source in the figure? Where is it located in the grid. I am hoping someone...

Homework Statement Homework Equations I could really use a push on how to approach this problem. My primary problem is it asks for the heat flux into the page, which makes no sense to me as that is the z direction and this is in the x/y plane. If anyone could explain this problem and maybe...

I am using the Finite Difference Method to solve Poisson's equation \frac{\partial \phi}{\partial z^2} = \frac{\rho}{\epsilon} To do it is discretized according to the Finite Difference Approximation of the second order derivative yielding the following set of equations for each grid point...

Homework Statement Which algebraic expressions must be solved when you use finite difference approximation to solve the following Possion equation inside of the square : $$U_{xx} + U_{yy}=F(x,y)$$[/B] $$0<x<1$$ $$0<y<1$$ Boundary condition $$U(x,y)=G(x,y)$$ Homework Equations Central...

Homework Statement Solve the time dependent 1D heat equation using the Crank-Nicolson method. The conditions are a interval of length L=1, initial distribution of temperature is u(x,0) = 2-1.5x+sin(pi*x) and the temperature in the ends of the interval are u(0,t) = 2; u(1,t) = 0.5. Homework...

Hi, I'm attempting to solve the 3D poisson equation ∇ ⋅ [ ε(r) ∇u ] = -ρ(r) Using a finite difference scheme. The scheme is simple to implement in 3D when ε(r) is constant, and I have found an algorithm that solves for a non-constant ε(r) in 2D. But I am having trouble finding an algorithm...

Suppose I want to solve the time-independent Schrödinger equation (ħ2/2m ∂2/∂x2 + V)ψ = Eψ using a numerical approach. I then discretize the equation on a lattice of N points such that x=(x1,x2,...,xN) etc. Finally I approximate the second order derivative with the well known central difference...

This question is inspired by Gilbert Strang's Course on Computational Science and Engineering, MIT 18.085. Consider the three matrices Fixed-Fixed $$K=\begin{bmatrix} 2 &-1 & 0 &0 \\ -1&2 & -1 &0 \\ 0 & -1 &2 & -1 \\ 0 & 0 & -1 & 2 \\ \end{bmatrix} $$ Free-Fixed $$T=\begin{bmatrix} 1 &-1 & 0 &0...

I am solving the simple 2nd-order wave equation: $$ \frac {\partial ^2 E}{\partial t^2} = c^2 \frac {\partial ^2 E}{\partial z^2} $$ Over a domain of (in SI units): ## z = [0,L=10]##m, ##t = [0,t_{max} = 10]##s and boundary/initial conditions: $$ E(z=0) = E(z=L) = 0 $$ $$ E(t=0) =...

Suppose I am given some 1D Hamiltonian: H = ħ2/2m d2/dx2 + V(x) (1) Which I want to solve on the interval [0,L]. I think most of you are familiar with the standard approach of discretizing the interval [0,L] in N pieces and using the finite difference formulas for V and the...

Homework Statement d2T/dx2 = 5*(dT/dx) - 0.1*x = 0 T(0) = 50 T(10) = 400 (Δx) = 2I've figured out how to do these problems when Δx = 1, but when it equals any other number it goes wrong. I know you start by plugging in the algebraic approximations for the differential elements, I think maybe my...

I'm trying to numerically solve the time dependent Schrödinger equation and I've been told that the best approach is to numerically integrate using a finite difference method, however I don't understand why I couldn't just use ode45 to solve it?! Is the finite difference (interpolation) method...

I found this part of forum the most relevant to this theme so excuse me if I missed. This year I'm doing a high school summer project related to quantum mechanics. Anyway I'm using finite difference method to solve Schrodinger equation. Before starting to work on a project I decided to get some...

Hi I've been trying to get a simple solution to the 2D Navier-Lame equations using finite difference on a rectangular grid. I want to see the displacements, u and v, when a simple deformation is imposed - e.g. top boundary is displaced by 10%. The equations are as follows: \begin{eqnarray*}...

Hello! (Wave) We consider the finite difference method for the approximation $\left\{\begin{matrix} -u''(x)+q(x)u(x)=f(x)\\ u'(a)=u'(b)=0 \end{matrix}\right.$ and let $K$ be the $(N+2) \times (N+2)$ matrix of the method. Let $v \in \mathbb{R}^{N+2}, v=\begin{pmatrix} v_0\\ v_1\\ \dots\\...

I'm trying to replicate the model presented in this [paper](http://www.sciencedirect.com/science/article/pii/S1359431103000474), which is basically to model heat and mass transfer along a one-dimensional duct. There are four characteristic equations for this problem : Momentum conservation...

Hello! (Wave) I want to solve numerically the following boundary value problem: $\left\{\begin{matrix} -u''+qu=f & , x \in [a,b]\\ -u'(a)+d_1 u(a)=0 & \\ u'(b)+d_2 u(b)=0 & \end{matrix}\right.$ where $q(x) \geq 0 \forall x \in [a,b], d_1, d_2 \geq 0$. We consider the uniform partition of...

I am simulating electrons inside a cylindrical well like the one shown on the first figure. My current work has been on solving the Schrodinger equation numerically for the above potential and then finding corrections to the solution such that it is consistent with Poissons equation. To do so...

Hi, I have written some codes for the finite difference solution of diffusion equation (\frac{\partial c}{\partial t}= D {\nabla^2 c}, where c is the species concentration and D is the diffusion coefficient) as follows: DO k= 1, tsteps+1 DO i = 2, zsteps DO j = 2, rsteps...

I am trying to solve the following eigenvalue differential equation numerically: ∇2ψ = Eψ , where the coordinate system is polar coordinates and the boundary condition is ψ(R,Φ)=0, where R is the radius of the disk i am working on. To solve it I am using a finite difference scheme, but there...

Hi PF! I was wondering if anyone could help me with a finite difference question? The problem I am doing is a 1-D space and time problem, so ##z## (space variable, from left to right) and ##t## (time) are my independent variables and my dependent variable is ##h##, the height, governed by a PDE...

Hi, Let's consider the heat equation as \frac{\partial T}{\partial t}=\alpha \frac{{{\partial }^{2}}T}{\partial {{x}^{2}}} In order to have a second accuracy system, one can use the Crank-Nicolson method as \frac{{{\partial }^{2}}T}{\partial {{x}^{2}}}\approx \frac{1}{2}\left(...

I had these code in this forum but comes out error as below, any suggestion? Error 1 error C4430: missing type specifier - int assumed. Note: C++ does not support default-int c:\users\username\documents\visual studio 2010\projects\fdm 001\fdm 001\explicit 001.cpp 27 Error 2...

hi dear i have a question. i have equation (1/α ) dT/dt =d2T/dr2 +1/r dT/dr +d2T/dz T=T(r,z) T(Ri,z)=Ti T(Ro,z)=To T(r,0)=To dT(r,L)/dz =0 by finite difference method O(h^3) and this question's MATLAB program. is there anyone who can do it ? it is very important for me tnx

I would like to reproduce results from a much older code to test a new one. I only have the old code's results in the form of plots, not data, but I need data. The older code solves the unsteady vorticity transport equation in 2D with a constant kinematic viscosity term. I'm interested in 2-D...

I want to solve a Laplace PDE in a polar coordinate system with finite difference method. and the boundary conditions: Here that I found in the internet: and the analytical result is: The question is how its works? Can I give an example or itd?Thanks

Homework Statement I need to (computationally) solve the following linear elliptic problem for the function u(x,y): \Delta u(x,y) = u_{x,x} + u_{y,y} = k u(x,y) on the domain \Omega = [0,1]\times[0,1] with u(x,y) = 1 at all points on the boundary.Homework Equations [/B] I know that I...