Surface Wave Analysis

The graphical user interface of Elmer
Post Reply
InsertCoin
Posts: 6
Joined: 08 Jan 2015, 08:20
Antispam: Yes

Surface Wave Analysis

Post by InsertCoin »

Hi :)
I have used Elmer for analyzing stress and strain usually.
And now, I have a chance to try it for solving transient wave equation.
I just wish to see the propagation of wave in time domain by using 3D square bar model.

I saw a video from youtube that shows "surface wave initiated by random drops", and i thought it might be possible to analyze the transient wave by using Elmer.

So I received some test examples through svn, however, I couldn't find any example for transient wave equation.

Anyhow, without any examples, I've tried to solve it by using helmholtz eqn, N-S eqn, but failed and failed and screwed.. :(

Here's my questions.
1. Is it possible to analyze transient wave equation by using up-to-date Elmer?
2. If it is possible, would you please recommend me equations that I have to use?
3. If you don't mind, would you please give some hints for the wave model?
4. Is there any way that I could get the project file of "surface wave initiated by random drops" for instance?

Sincerely,

-Sangkyun Kim
raback
Site Admin
Posts: 4812
Joined: 22 Aug 2009, 11:57
Antispam: Yes
Location: Espoo, Finland
Contact:

Re: Surface Wave Analysis

Post by raback »

Hi

There is a test case "WaveEq" that includes a wave solver. There is somewhat modified version of that with Sommerfeld BCs which is not in the distro. I might put it there though. Have to polish it a little time permitting.

-Peter
InsertCoin
Posts: 6
Joined: 08 Jan 2015, 08:20
Antispam: Yes

Re: Surface Wave Analysis

Post by InsertCoin »

Dear Peter,

Thank you for your kind reply.

I am using windows and Elmer GUI now.
When I tried to perform the test example, it seemed I need WaveEq.dll file to solve.

How can I implement solver .dll file in the example, which is "WaveEq.f90"?
I've tried to make it as .dll but I got plenty of errors..

Would you please give some notice for me how to run this example?
or
Is there any ways to get the required dll files?

Sincerely

Sangkyun Kim
annier
Posts: 1168
Joined: 27 Aug 2013, 13:51
Antispam: Yes

Re: Surface Wave Analysis

Post by annier »

Hi Kim Sangkyun,
In windows, we convert the filename.f90 file into filename.dll file by the following command after jumping into the directory containing the filename.f90 file:

Code: Select all

>elmerf90 -o filename.dll filename.f90
So, in the context of WaveEq.f90

Code: Select all

>elmerf90 -o WaveEq.dll WaveEq.f90
In Ubuntu, we convert the filename.f90 file into filename.so file by the following command:

Code: Select all

$elmerf90 -o filename.so filename.f90
and, in context of WaveEq.f90,

Code: Select all

$elmerf90 -o WaveEq.so WaveEq.f90
For detailed information, you can refer the Elmer Programmers Manual.

Yours
Anil Kunwar
Anil Kunwar
Faculty of Mechanical Engineering, Silesian University of Technology, Gliwice
InsertCoin
Posts: 6
Joined: 08 Jan 2015, 08:20
Antispam: Yes

Re: Surface Wave Analysis

Post by InsertCoin »

Dear Anil Kunwar,

Thank you very much for your help!!
It is nicely worked :)
This would be good starting point for surface wave analysis in 3D bar model :lol:

Sincerely,

Sangkyun Kim
InsertCoin
Posts: 6
Joined: 08 Jan 2015, 08:20
Antispam: Yes

Re: Surface Wave Analysis

Post by InsertCoin »

Hi :)

With your graceful help, I've proceeded one step further.

However, unfortunately, I meet a problem as I attach figure file.

Wave propagation direction is opposite to my intended direction.

My intention was
1. wave is generated at the left side edge
2. wave is propagated to the right side

But the analysis result was
1. position of wave generation is moving from left to right side (thus, red bar is generated more)
2. wave itself is moving to left side

it seems like there is a vacuum cleaner at the left side to sucking, and therefore, wave generated.

Here's my SIF file.

Check Keywords "Warn"

Header
Mesh DB "."
End

Simulation
Coordinate System = "Cartesian 2D"
Simulation Type = Transient
Timestep Sizes = 0.01
Timestep Intervals = 1000
Output Intervals(1) = 1
Steady State Max Iterations = 1
Max Output Level = 3
End

Body 1
Equation = 1
Initial Condition = 1
Material = 1
End

Initial Condition 1
Potential = Real 0
Potential Velocity = Real 0
End

Equation 1
Active Solvers(2) = 1 2
End

Material 1
Name = "Material1"
Density = 1
Viscosity = 1.0
Permittivity = Real 1.0
Reference Temperature = Real 293.0
End

Solver 1
Equation = "Poisson"
Variable = "Potential"
Variable DOFs = 1
Element = "p:4"
Procedure = "WaveEq" "WaveSolver"
Linear System Solver = "Direct"
Linear System Direct Method = Umfpack
Steady State Convergence Tolerance = 1e-09
End

Solver 2
Equation = Result Output
Output Format = Vtu
Output File Name = case
Procedure = "ResultOutputSolve" "ResultOutputSolver"
Exec Solver = Always
End

Boundary Condition 1
Target Boundaries(1) = 3
Potential = Variable Time
Real MATC "10*sin(2*3.1415*tx(0))"
Potential Velocity = Variable Time
Real MATC "-10*sin(2*3.1415*tx(0))"
End
Attachments
1.png
(42.42 KiB) Not downloaded yet
InsertCoin
Posts: 6
Joined: 08 Jan 2015, 08:20
Antispam: Yes

Re: Surface Wave Analysis

Post by InsertCoin »

I add a more attachment, a .gif file of the simulation.

how can I solve this?... do you have any idea?..

Sincerely,

Sangkyun Kim
Attachments
Wave Propagation.gif
Wave Propagation.gif (942.94 KiB) Viewed 7990 times
Fouzi
Posts: 55
Joined: 14 Mar 2017, 18:03
Antispam: Yes

Re: Surface Wave Analysis

Post by Fouzi »

Hi Sangkyun Kim,

could you help me please?, i have a propbleme with the propagation in elastic plate.


Fouzi
Post Reply