Scribd
Upload
4 views3 pages

Adaptive FEM Detailed Notes

The document discusses the implementation of Adaptive Finite Element Method (FEM) with explicit time integration for solving the wave equation, detailing its strong and weak formulations, finite element discretization, and the leapfrog scheme for time-stepping. It emphasizes the importance of the CFL condition, local time-stepping, and a posteriori error estimation for adaptive algorithms. Additionally, it highlights potential extensions to nonlinear wave equations and applications in seismic studies.

Uploaded by

mritunjayddu
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
4 views3 pages

Adaptive FEM Detailed Notes

The document discusses the implementation of Adaptive Finite Element Method (FEM) with explicit time integration for solving the wave equation, detailing its strong and weak formulations, finite element discretization, and the leapfrog scheme for time-stepping. It emphasizes the importance of the CFL condition, local time-stepping, and a posteriori error estimation for adaptive algorithms. Additionally, it highlights potential extensions to nonlinear wave equations and applications in seismic studies.

Uploaded by

mritunjayddu
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

Research Notes: Adaptive FEM with Explicit Time

Integration for Wave Equation

1. Strong Form of Wave Equation

u_tt - div(c^2 grad u) = f(x,t) with initial and boundary conditions.

2. Weak Formulation

Multiply by test function v in H1 and integrate: (u_tt, v) + ∫ c^2 ∇u·∇v dx = (f, v). Define bilinear form
a(u,v) = ∫ c^2 ∇u·∇v dx.

3. Finite Element Discretization

Approximate solution in finite-dimensional space Vh. u_h(x,t) = Σ U_i(t) φ_i(x). Leads to system: M
U'' + K U = F.

4. FEM Matrices

Mass matrix: M_ij = ∫ φ_i φ_j dx. Stiffness matrix: K_ij = ∫ c^2 ∇φ_i·∇φ_j dx.

5. Leapfrog Scheme

Second-order explicit scheme: U^{n+1} = 2U^n - U^{n-1} + ∆t^2 M^{-1}(F^n - K U^n).

6. CFL Condition

∆t ≤ C h_min. Small mesh size → very small time-step.

7. Local Time-Stepping

Use smaller time-step only in refined regions to overcome CFL restriction.

8. A Posteriori Error Estimation

Error bounded by sum of indicators: 1. Spatial residual 2. Time discretization 3. Mesh change 4.
LTS error 5. Data error.

9. Residual-Based Error Indicator


Element residual: h_K^2 ||f + div(c^2 grad u_h)||. Edge jump terms included.

10. Mesh Change Error

Error from projection between meshes controlled using interpolation estimates.

11. Adaptive Algorithm

Compute → estimate error → refine/coarsen → repeat.

12. Functional Analysis Framework

Uses Hilbert spaces H1, L2 and energy norms.


13. Research Insights

Combines explicit schemes with adaptivity efficiently. Important for wave propagation and
large-scale simulations.

14. Possible Extensions

Nonlinear wave equations, heterogeneous media, high-order FEM, seismic applications.

You might also like