Problem 17
Question
The Helmholz equation, \(p_{1 x x}+p_{y y}+p, z_{r z}+(\omega / c)^{2} p=0\), governs acoustic modes of vibration in a cavity with rigid walls. Here \(p=p(x, y, z)\) represents the amplitude of sinusoidally varying pressure, \(\omega\) is the circular frequency, and \(c\) is the speed of sound in the medium. The boundary condition is \(p_{1 n}=0\), where \(n\) is a direction normal to the wall. Derive formulas for finite element matrices, using the assumed pressure amplitude field \(p=\) \(\left\lfloor\mathbf{N} \mid\left\\{\mathbf{p}_{e}\right\\} .\right.\)
Step-by-Step Solution
Verified Answer
The finite element matrices for the given Helmholz equation can be derived by first reformulating the equation to its weak form by integrating over the volume of the cavity. Applying Galerkin's method, the weak form of equation is then discretized using a finite set of basis functions. This results in a system of linear equations with a stiffness matrix and a mass matrix. The final system equation can be expressed as \( (K + \frac{\omega^{2}}{c^{2}}M)p_{e} = 0 \). Lastly, the boundary condition is applied by modifying the elements of \(K\) and \(M\) corresponding to the degrees of freedom at the boundary.
1Step 1: Formulating the Weak Form of the Equation
Infinite element analysis, instead of directly solving the differential equation, it is re-formulated into its weak form. The weak form is obtained in this case by multiplying the Helmholtz equation by a test function \( v \), and integrating over the volume \( V \) of the cavity. The weak form of the Helmholz equation can then be represented as: \( \int_{V} \left( \nabla v \cdot \nabla p + \frac{\omega^{2}}{c^{2}}v p \right) dV = 0 \).
2Step 2: Applying Galerkin's method
The next step would be to discretize the weak form of the equation using a finite set of basis functions. In this exercise, we are provided with a basis function, \( \mathbf{N} \). We substitute \( p \) and \( v \) with approximations, \( \mathbf{N}p_{e} \) and \( \mathbf{N} \), into the weak form equation. This application of the Galerkin method reduces the problem to a system of linear equations.
3Step 3: Assembling Element Matrices
By integrating and applying the divergence theorem, the weak form of the Helmholz equation is turn into matrix form which yields two finite element matrices: a stiffness matrix and a mass matrix. Denoting \( K_{ij} = \int_{V} \nabla N_{i} \cdot \nabla N_{j} dV \), and \( M_{ij} = \int_{V} N_{i}N_{j} dV \), the finite element matrices can be represented as \( K \) and \( M \). The final form of the system is then \( (K + \frac{\omega^{2}}{c^{2}}M)p_{e} = 0 \), where \( p_{e} \) contains the coefficients of the pressure field in the finite element basis.
4Step 4: Applying the Boundary Condition
To implement the boundary condition \(p_{1n} = 0\), the elements of \(K\) and \(M\) corresponding to the degrees of freedom at the boundary are modified. The nature of the modification depends on how the boundary is discretized, whether with elements that touch the boundary or with boundary elements.
Key Concepts
Helmholtz equationWeak formGalerkin methodStiffness matrixMass matrix
Helmholtz equation
The Helmholtz equation is a partial differential equation used to describe various physical phenomena, such as acoustic vibrations in a cavity with rigid walls.
This equation is represented as:\[ p_{1xx} + p_{yy} + p_{zz} + \left( \frac{\omega}{c} \right)^2 p = 0 \]where:
This equation is represented as:\[ p_{1xx} + p_{yy} + p_{zz} + \left( \frac{\omega}{c} \right)^2 p = 0 \]where:
- \( p = p(x, y, z) \) is the amplitude of sinusoidally varying pressure.
- \( \omega \) is the circular frequency.
- \( c \) is the speed of sound in the medium.
Weak form
In finite element analysis, the weak form of an equation is an alternative formulation that's more amenable to numerical methods.
Instead of solving the original Helmholtz equation directly, the weak form is obtained by multiplying the equation by a test function \( v \) and integrating over the entire volume \( V \).
Mathematically, this can be written as:\[ \int_{V} \left( abla v \cdot abla p + \frac{\omega^{2}}{c^{2}}v p \right) dV = 0 \]The weak form transforms the differential equation into an integral one, which can be easier to handle when constructing finite element matrices.
It also accommodates irregular geometries and boundary conditions more naturally.
Instead of solving the original Helmholtz equation directly, the weak form is obtained by multiplying the equation by a test function \( v \) and integrating over the entire volume \( V \).
Mathematically, this can be written as:\[ \int_{V} \left( abla v \cdot abla p + \frac{\omega^{2}}{c^{2}}v p \right) dV = 0 \]The weak form transforms the differential equation into an integral one, which can be easier to handle when constructing finite element matrices.
It also accommodates irregular geometries and boundary conditions more naturally.
Galerkin method
The Galerkin method is a technique used to solve the weak form of a differential equation.
Its essence lies in approximating the solution by a finite set of basis functions.
For the Helmholtz equation, these approximations for the solution \( p \) and the test function \( v \) involve basis functions given by \( \mathbf{N} \).
This approach leads to a system of linear equations, facilitating computation by breaking down the continuous problem into a discrete one.
The Galerkin method is crucial because it ensures that the errors for the governing equations' residuals are orthogonal to the chosen finite-dimensional subspace, which effectively minimizes the global error within that space.
Its essence lies in approximating the solution by a finite set of basis functions.
For the Helmholtz equation, these approximations for the solution \( p \) and the test function \( v \) involve basis functions given by \( \mathbf{N} \).
This approach leads to a system of linear equations, facilitating computation by breaking down the continuous problem into a discrete one.
The Galerkin method is crucial because it ensures that the errors for the governing equations' residuals are orthogonal to the chosen finite-dimensional subspace, which effectively minimizes the global error within that space.
Stiffness matrix
The stiffness matrix, denoted as \( K \), is a key component in the finite element method.
It arises from discretizing the weak form of the Helmholtz equation by applying the divergence theorem and integrating the basis functions: \[ K_{ij} = \int_{V} abla N_{i} \cdot abla N_{j} \, dV \]This matrix quantifies the relationship between different nodes in the element, representing the change in the system due to applied forces or constraints.
Computing \( K \) accurately is essential because it determines how the structure deforms and vibrates under load.
It arises from discretizing the weak form of the Helmholtz equation by applying the divergence theorem and integrating the basis functions: \[ K_{ij} = \int_{V} abla N_{i} \cdot abla N_{j} \, dV \]This matrix quantifies the relationship between different nodes in the element, representing the change in the system due to applied forces or constraints.
Computing \( K \) accurately is essential because it determines how the structure deforms and vibrates under load.
Mass matrix
The mass matrix, denoted as \( M \), accompanies the stiffness matrix in the finite element method.
It represents mass distribution across the element and is given by:\[ M_{ij} = \int_{V} N_{i}N_{j} \, dV \]The mass matrix is crucial because it relates to the inertia of the system, specifically how mass is exchanged between nodes and how it affects the system's dynamic response.
When combined with the stiffness matrix in the equation \((K + \frac{\omega^{2}}{c^{2}}M)p_{e} = 0\), the mass matrix aids in simulating the vibrational characteristics of the system, essential for accurate dynamic analysis.
It represents mass distribution across the element and is given by:\[ M_{ij} = \int_{V} N_{i}N_{j} \, dV \]The mass matrix is crucial because it relates to the inertia of the system, specifically how mass is exchanged between nodes and how it affects the system's dynamic response.
When combined with the stiffness matrix in the equation \((K + \frac{\omega^{2}}{c^{2}}M)p_{e} = 0\), the mass matrix aids in simulating the vibrational characteristics of the system, essential for accurate dynamic analysis.
Other exercises in this chapter
Problem 8
Consider a uniform, simply supported beam of length \(L\). For each of the following loadings, use the Galerkin method to determine constant \(a\) in the approx
View solution Problem 10
The equation of motion of a string is \(T w_{x x x}-\rho_{L} \ddot{w}=0\), where \(T=\) constant axial tension, \(w=\) small lateral displacement, \(x=\) axial
View solution Problem 19
The differential equation for wind-driven circulation in a shallow lake is $$ \psi_{, x x}+\psi_{1 y y}+A \psi_{x}+B \psi_{y}+C=0 $$ where \(\psi\) is the strea
View solution Problem 5
Consider the differential equation \(u_{3 x}+2 u-16 x=0\) with the boundary condition \(u=0\) at \(x=0 .\) The exact solution is \(u=4\left(e^{-2 x}-1\right)+8
View solution