Implement the fdtd method in a matlab program

Reference no: EM131307706

1. The 2D diffusion equation u_t = D∇²u is discretized using the Finite Difference method as

(u_i,jⁿ⁺¹- u_i,jⁿ)/DΔt = γ[(u_i+1ⁿ⁺¹_,j - 2u_i,jⁿ⁺¹ + u_i-1ⁿ⁺¹_,j)/Δx² + u_i,j+1ⁿ⁺¹_,j - 2u_i,jⁿ⁺¹ + u_i-1ⁿ⁺¹/Δy²] + (1-γ)[(u_i+1ⁿ_,j - 2u_i,jⁿ + u_i-1ⁿ_,j)/Δx² + u_i,j+1ⁿ - 2u_i,jⁿ + u_i,j-1ⁿ/Δy²]

where γ is a parameter between 0 and 1 which determines the method of time discretization. Specifically, we have

γ = 0 : explicit method
γ = 1 : implicit method
γ = 1/2 : Crank-Nicholson method

Derive the stability condition for each of the above time discretization methods.

2. The propagation of electromagnetic waves in 2D is governed by the wave equation

n²/c^2.∂²u/∂²t = ∂²u/∂x² + ∂u/∂y²,

where n(x, y) is the refractive index of the medium and c = 3 × 10⁸ m/s is the speed of light in vacuum. We would like to solve the above equation using the Finite Difference Time Domain (FDTD) method. The computation domain is restricted to a rectangular region of size a × b and Radiating Boundary Conditions based on the one-way wave equation are applied to all four boundaries.

(a) Give a FD discretization of the wave equation for an interior node (i, j).

(b) Derive the FD equations for nodes on the left, right, top and bottom boundaries.

(c) Implement the FDTD method in a MATLAB program to solve the above wave equation. An outline of the program is given at the end of the assignment for your reference.

(d) Use your program to run the following simulations

Reference no: EM131307706

Questions Cloud

Create diagram of dsl and cable modem connection to your isp : Create a diagram of the DSL and Cable Modem connections to your ISP, cable organization, and telecom to your home router using Visio or its open source alternative software.

Define the security strategies of defense in depth security : In this paper, you will define the security strategies of Defense in Depth and Layered Security along with comparing and contrasting the strategies by explaining.

Write a program where you create a class called bankaccount : Write a program where you create a class called BankAccount. The BankAccount class should contain a String to store the customer name and a double to store the account balance.

What is conflict management : What is conflict management.watch the video lecture below to learn why conflict at work isn't always a bad thing.

Implement the fdtd method in a matlab program : ECE 452 - Assignment Give a FD discretization of the wave equation for an interior node (i, j) - Derive the FD equations for nodes on the left, right, top and bottom boundaries - Derive the stability condition for each of the above time discretizatio..

Find the price elasticity of demand for newtons donuts : Calculate the price elasticity of demand for Newton's Donuts and describe what it means. Describe your answer and show your calculations.

Explain effective approaches to the broad spectrum : Explain effective approaches to the broad spectrum of employee relations, including career development, fostering ethical behavior, discipline, labor relations, and dismissals.

Analyze and derive an expression for the marginal cost curve : Describe and derive an expression for the marginal cost (MC) curve. Describe and estimate the incremental costs of the extra 200 pairs per week (from 1,000 pairs to 1,200 pairs of shoes).

By how much was the number of references reduced : Using the appropriate psychology database at your library, enter the subject words short-term memory and see how many references you obtain.

User Account

All Pages