Scribd
Upload
14 views2 pages

Wave Motion Simulation and Analysis

The document discusses wave motion and its mathematical modeling, particularly focusing on simulating one-dimensional waves on a string using the wave equation. It outlines the initial and boundary conditions necessary for solving the wave equation and describes a finite difference method for numerical approximation. Additionally, it presents an example with numerical errors and convergence accuracy related to the method used.

Uploaded by

amareshbej
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
14 views2 pages

Wave Motion Simulation and Analysis

The document discusses wave motion and its mathematical modeling, particularly focusing on simulating one-dimensional waves on a string using the wave equation. It outlines the initial and boundary conditions necessary for solving the wave equation and describes a finite difference method for numerical approximation. Additionally, it presents an example with numerical errors and convergence accuracy related to the method used.

Uploaded by

amareshbej
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

Project 1

Rakhi Bala Roy


April 2025

1 Introduction
A very wide range of physical processes lead to wave motion, where signals are propagated through a
medium in space and time,normally with the tittle or no permanent movement of the medium itself. The
shape of the signals may undergo changes as they travel through matter,but usually not so much that
the signals cannot be recognized at some later point in space and [Link] types of wave motion can
be described by the equation utt = ∇.(c2 ∇u) + f , which we will solve in the forthcoming text by finite
difference methods [1].

1.1 Simulation of waves on a string


We begin our study of wave equations by simulating one-dimensional waves on a string, say on a guitar
or violin. Let the string in the deformed state coincide with the interval [0,L] on the x axis, and let
u(x,t) be the displacement at time t in the y direction of a point initially at x. The displacement function
u is governed by the mathematical model

∂2u ∂2u
2
= c2 2 x ∈ (0, L), t ∈ (0, T ] (1)
∂t ∂x

u(x, 0) = I(x), x ∈ [0, L] (2)

∂
u(x, 0) = 0, x ∈ [0, L] (3)
∂t

u(0, t) = 0, t ∈ (0, T ] (4)

u(L, t) = 0, t ∈ (0, T ] (5)


The constant c and the function I(x) must be prescribed.
Equation (1) is known as the one-dimensional wave equation. Since this PDE contains a second-order
derivative in time, we need two initial conditions. The condition (2) specifies the initial shape of the
string, I(x), and (3) expresses that the initial velocity of the string is zero. In addition, PDEs need
boundary conditions, given here as (4) and (5). These two conditions specify that the string is fixed at
the ends, i.e.,that the displacement u is zero. The solution u(x,t) varies in space and time and describes
waves that move with velocity c to the left and right. Sometimes we will use a more compact notation
for the partial derivatives to save space;
∂2u
ut = ∂u
∂t , utt = ∂t2 ,
and similar expressions for derivatives with respect to other variables. Then the wave equation can
be written compactly as utt = c2 uxx .
The PDE problem (1)-(5) will now be discretized in space and time by a finite difference method.
24
Example 1. We take λ1 = λ2 = 1, T = 1, f (t) = 12t2 + r(5−α) t4 − α
The exact solution of Example 1 is u(t) = t4 . In our computation, we use the difference scheme to solve
this problem. The numerical error was computed by
E∞ (τ ) = max0≤n≤N |u(tn ) − un  
E∞ (τ )
and the corresponding convergence order was calculated by Order=log2 E∞ ( τ2 ) .
We chose α = 1.3, 1.5, 1.7, respectively, to rest the convergence accuracy.

1
Table 1 shows the numerical errors and the corresponding convergence accuracy for each fixed α. As
expected, our scheme generates the numerical approximation accuracy with Q(τ 3−α ).
Table 1.

Students Name Roll Numbers Reg. No.


Rakhi Bala Roy 16 1825

Table 1: Numerical approximation error and convergence accuracy in Example 1.

References
[1] Manruo Cui, Cui-Cui Ji, and Weizhong Dai. A finite difference method for solving the wave equation
with fractional damping. Mathematical and Computational Applications, 29(1), 2024.

Common questions

Powered by AI

Finite difference methods are used to discretize the one-dimensional wave equation in both space and time, enabling its numerical solution. Key aspects include defining the initial conditions, such as the initial shape and velocity of the string, and boundary conditions specifying that the string is fixed at both ends, resulting in zero displacement at these points . The numerical error of these methods is evaluated using specific schemes, where the convergence order is determined by the formula Order=log2(E∞(τ)/E∞(τ/2)), and the scheme's accuracy varies with different α values .

The convergence accuracy of a numerical scheme can be tested by calculating numerical errors using a specific formula, E∞(τ) = max0≤n≤N|u(tn) −un, and then determining the convergence order using Order=log2(E∞(τ)/E∞(τ/2)). In the discussed example, different values of α (1.3, 1.5, 1.7) were tested, and it was found that the scheme's accuracy aligns with Q(τ^(3−α)), an expected result which confirms the validity and efficiency of the finite difference approach for solving wave equations .

Calculating numerical error is crucial because it provides a quantitative measure of how closely the finite difference scheme approximates the true solution of the wave equation. It helps in assessing the accuracy and reliability of the method. If the numerical error is significant, it indicates potential issues in the discretization approach or parameter selections. By evaluating the numerical error and its convergence, one can ensure that the scheme is adequately capturing the dynamics of wave motion and adjust the model as needed to improve precision .

Setting boundary conditions like u(0, t) = 0 and u(L, t) = 0 ensures that the ends of the string remain fixed, meaning the displacement is zero at these points throughout the simulation. These conditions effectively confine the wave to travel within the defined length of the string, impacting how the wave reflects and interacts with these boundaries, thereby shaping the wave patterns observed in the simulation. This setup simulates a realistic scenario for wave motion on instruments such as guitars or violins, influencing the wave dynamics and harmonics .

Establishing initial and boundary conditions is crucial for solving the wave equation as they define the specific constraints and starting points of the problem, allowing for an accurate representation of the physical system. The initial conditions set the initial shape and velocity of the string, while the boundary conditions ensure that the string's displacement is zero at its fixed endpoints. These conditions are necessary to solve the partial differential equation uniquely and accurately capture the vibrations and resultant wave patterns over time .

The condition ∂tu(x, 0) = 0 specifies that the initial velocity of the string is zero across its length. This is significant as it defines a resting state of the system before any wave motion occurs, ensuring that any subsequent movement is purely a result of the wave propagation rather than initial motion—providing a clear starting point for analyzing the wave dynamics generated by subsequent external forces or initial displacements .

The term f in the equation represents external forces or influences acting upon the wave system. Its presence in the wave equation adds a component that accounts for these forces, modifying the wave dynamics beyond what would occur due solely to inherent properties of the medium and initial conditions. By including f, the model can capture a broader range of real-world scenarios where external disturbances or forces alter wave behavior .

The velocity parameter 'c' determines the speed at which waves propagate through the medium along the string. It influences the shape and timing of wave patterns as they travel left and right from their initial positions. In the one-dimensional wave equation, c appears in the term c^2uxx, affecting the rate at which changes in position lead to changes in displacement over time .

The equation utt = ∇.(c^2∇u) + f describes the propagation of waves through a medium while considering external forces denoted by f. The equation reflects changes over time and space, capturing the essence of wave motion where signals travel through a medium without permanent displacement of the medium itself. It is applied in simulations by discretizing the problem using finite difference methods to solve for wave behavior over time and space, as demonstrated in the study of waves on a string .

Discretizing partial derivatives in space and time is essential for applying finite difference methods to solve wave equations because it transforms continuous problems into a set of algebraic equations that can be solved computationally. This discretization is important for approximating the values of derivatives at discrete grid points, allowing the simulation of waves to account for changes across both space and time in a manageable manner. It enables the practical computation of solutions on a finite grid and is a fundamental process that underpins the numerical analysis of the wave equation's solutions .

You might also like