Time-Harmonic Elastic Wave Problem

More
5 years 4 months ago #1701 by Amad
Time-Harmonic Elastic Wave Problem was created by Amad
Hello everyone!

I am trying to solve the time-harmonic elastic wave problem using the Galerkin method.
But the results are not correct.

I am attaching a jupyter-notebook with the case that I tried. Could you tell me what I am doing wrong?
In addition, I would like to include a pointwise source in this problem, could you also help me with this?

Finally, how can I use the complex conjugate to measure L2-norm error?
Code:
gfu = GridFunction(fes) u_exact = CoefficientFunction(r_source) u = CoefficientFunction(gfu) #Error analysis print ("Displacement L2-error:", sqrt (Integrate ( (u-u_exact)*(u-u_exact), mesh)))

Thank you in advance for your time.

Kind Regards,

Alan
Attachments:
More
5 years 4 months ago #1709 by joachim
Replied by joachim on topic Time-Harmonic Elastic Wave Problem
Hi Alan,

it looks like you wanted to specify Dirichlet boundary conditions, but you gave Neumann ...

For the beginning, replace the point-source by a small Gaussian peak,

Joachim
The following user(s) said Thank You: Amad
More
5 years 4 months ago #1725 by Amad
Replied by Amad on topic Time-Harmonic Elastic Wave Problem
Dear Joachim,

Thank you for your answer. It's working now! :)
How to use L^2 norm for a complex number?

I tried to include a complex conjugate in the code below, but it didn't work.
Code:
#Error analysis print ("Displacement L2-error:", sqrt (Integrate ( (u-u_exact)*(u-u_exact), mesh)))

Thank you in advance.

Best Regards,

Alan
More
5 years 3 months ago #1730 by christopher
Replied by christopher on topic Time-Harmonic Elastic Wave Problem
Hi Amad,
you can either use the Norm CoefficientFunction or the Conj (conjugate).]
Best
Christopher
The following user(s) said Thank You: Amad
Time to create page: 0.101 seconds