How to Implement FEM with a Discontinuity in the Exact Solution?

  • Thread starter Thread starter ageralo
  • Start date Start date
  • Tags Tags
    Discontinuity Fem
Click For Summary
SUMMARY

The discussion focuses on implementing Finite Element Methods (FEM) to solve a problem with a discontinuous exact solution. The exact solution is defined as u = x^5/20 - x/20 for x in the interval [0, 0.6) and u = sin(x) for x in (0.6, 1]. The user has successfully set up the matrices Au = f but is uncertain how to handle the discontinuity in the code. The problem emphasizes the need for a robust approach to manage discontinuities in FEM implementations.

PREREQUISITES
  • Understanding of Finite Element Methods (FEM)
  • Familiarity with numerical methods for solving differential equations
  • Knowledge of matrix operations and formulations in FEM
  • Basic understanding of piecewise functions and discontinuities
NEXT STEPS
  • Research techniques for handling discontinuities in FEM, such as mesh refinement or special boundary conditions
  • Learn about piecewise linear basis functions in FEM
  • Explore the implementation of jump conditions in numerical simulations
  • Study the application of the Galerkin method in FEM for discontinuous solutions
USEFUL FOR

Students and professionals in computational mechanics, engineers implementing FEM for simulations, and researchers dealing with discontinuous solutions in numerical analysis.

ageralo
Messages
4
Reaction score
0

Homework Statement



Use FEM to solve this problem. The difficulty lies in the fact that the exact solution has a discontinuity in it. From x=[0,0.6) the exact solution u is x5/20 - x/20 and from x =(0.6,1] u is sin(x). The problem I'm having is I'm not sure what to do at the jump in my code. I have already set up my matrices Au = f, and it solves it fine given no discontinuity.

uxx = f

Homework Equations



u(0) = u0
u(1) = u1
 
Physics news on Phys.org
ageralo said:

Homework Statement



Use FEM to solve this problem. The difficulty lies in the fact that the exact solution has a discontinuity in it. From x=[0,0.6) the exact solution u is x5/20 - x/20 and from x =(0.6,1] u is sin(x). The problem I'm having is I'm not sure what to do at the jump in my code. I have already set up my matrices Au = f, and it solves it fine given no discontinuity.

uxx = f

Homework Equations



u(0) = u0
u(1) = u1

What is FEM?

RGV
 

Similar threads

Replies
2
Views
2K
Replies
4
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 6 ·
Replies
6
Views
6K
Replies
9
Views
3K
  • · Replies 3 ·
Replies
3
Views
4K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 4 ·
Replies
4
Views
4K
  • · Replies 5 ·
Replies
5
Views
3K
Replies
3
Views
17K