Scribd
Upload
47 views1 page

Problema

This 3-sentence summary provides the key details about the document: The document outlines Homework Problem #3 which involves solving a one-dimensional wave equation using various numerical methods, including higher order time-stepping schemes, compact spatial differencing, and pseudo-spectral methods with FFT, and examining the effect of adding weak spatial smoothing consistent with the spatial discretization. Students are asked to implement the problem using MATLAB, solve it using different grid sizes and numerical methods, and plot and compare the solutions to the exact solution.

Uploaded by

César Augusto Niño Castro
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
47 views1 page

Problema

This 3-sentence summary provides the key details about the document: The document outlines Homework Problem #3 which involves solving a one-dimensional wave equation using various numerical methods, including higher order time-stepping schemes, compact spatial differencing, and pseudo-spectral methods with FFT, and examining the effect of adding weak spatial smoothing consistent with the spatial discretization. Students are asked to implement the problem using MATLAB, solve it using different grid sizes and numerical methods, and plot and compare the solutions to the exact solution.

Uploaded by

César Augusto Niño Castro
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 1

Numerical Methods for Coastal and Oceanographic Engineering Slinn Spring 2004

Home Work Problem # 3. Due March 4, 2004.

This problem set will give you experience with, 1) periodic boundary conditions, 2) the
compact scheme, 3) higher order explict time step methods, 4) fast Fourier transform
(FFT) methods, and 5) numerical diffusion and dissipation.

u u
Again we examine the one way wave equation, +C = 0.
t x

1. Solve using the third order Adams-Bashforth time step method. Use a Courant
2
number of 0.2 and initial conditions u ( x,0) = 3e 5( x 2) , with C = 1 m s-1, on the
domain, 0 < x < 10 m. Use periodic boundary conditions, and plot solutions with
t = 10 s and t = 20 s. Use three values of x, with 64, 128, and 256 grid points.
Use 3 types of spatial differences; 1) second-order central differences, 2) 4th order
compact scheme, 3) pseudo-spectral methods using FFTs to calculate the spatial
derivatives. Plot compared to the exact solution, which is the same as the initial
conditions.

2. Add weak spatial smoothing to each solution, consistent with the type of spatial
derivatives being used, plot compared to the exact solution.

A) second order spatial smoothing with the coefficient chosen to damp the
2x wave. (use 64 and 256 grid points)
B) 4th-order compact smoothing with = 0, = 0.4, a = 0.925, b = 0.9,
c = -0.025 (see page 57 of thesis notes). Use 64 and 128 grid points.
C) Spectral cutoff filter, for de-aliasing kcutoff = 2/3 kmax. Use 32 and 64 grid
points only.

A periodic tri-diagonal matrix solver and a FFT package is provided for the assignment.
The FFT package provided uses only powers of 2 so the number of grid points needs to
be 128, 256, 512, a power of 2 (or 129, 257, 513,) depending on method of
periodicity.

Plot the results with Tecplot.

You might also like