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
Surface Wave Analysis
-
- Site Admin
- Posts: 4832
- Joined: 22 Aug 2009, 11:57
- Antispam: Yes
- Location: Espoo, Finland
- Contact:
Re: Surface Wave Analysis
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
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
-
- Posts: 6
- Joined: 08 Jan 2015, 08:20
- Antispam: Yes
Re: Surface Wave Analysis
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
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
Re: Surface Wave Analysis
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:
So, in the context of WaveEq.f90
In Ubuntu, we convert the filename.f90 file into filename.so file by the following command:
and, in context of WaveEq.f90,
For detailed information, you can refer the Elmer Programmers Manual.
Yours
Anil Kunwar
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
Code: Select all
>elmerf90 -o WaveEq.dll WaveEq.f90
Code: Select all
$elmerf90 -o filename.so filename.f90
Code: Select all
$elmerf90 -o WaveEq.so WaveEq.f90
Yours
Anil Kunwar
Anil Kunwar
Faculty of Mechanical Engineering, Silesian University of Technology, Gliwice
Faculty of Mechanical Engineering, Silesian University of Technology, Gliwice
-
- Posts: 6
- Joined: 08 Jan 2015, 08:20
- Antispam: Yes
Re: Surface Wave Analysis
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
Sincerely,
Sangkyun Kim
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
Sincerely,
Sangkyun Kim
-
- Posts: 6
- Joined: 08 Jan 2015, 08:20
- Antispam: Yes
Re: Surface Wave Analysis
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
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
-
- Posts: 6
- Joined: 08 Jan 2015, 08:20
- Antispam: Yes
Re: Surface Wave Analysis
I add a more attachment, a .gif file of the simulation.
how can I solve this?... do you have any idea?..
Sincerely,
Sangkyun Kim
how can I solve this?... do you have any idea?..
Sincerely,
Sangkyun Kim
- Attachments
-
- Wave Propagation.gif (942.94 KiB) Viewed 8159 times
Re: Surface Wave Analysis
Hi Sangkyun Kim,
could you help me please?, i have a propbleme with the propagation in elastic plate.
Fouzi
could you help me please?, i have a propbleme with the propagation in elastic plate.
Fouzi