mmolero / SimNDT

Ultrasonic NDT Simulator with engine core based on the Elastodynamic Finite Integration Technique (EFIT)
47 stars 13 forks source link

Acoustic ultrasonic wave propagation #15

Open Nicko222 opened 2 years ago

Nicko222 commented 2 years ago

Dear @mmolero, Thanks for your amazing ultrasound simulator!!! I would like to know if I use it in acoustic wave propagation how to get the pressure field. e.g., in the homogeous 2D acoustic wave equation:

wave.

In the equation, v represents velocity, u represents the pressure response in an acoustic medium (known as the “wavefield”). If the output of Txx and Tyy have the same meaning of wavefield u. Hope to receive your reply. Thanks a lot!

mmolero commented 1 year ago

Please see the issue #14.

Nicko222 commented 1 year ago

Please see the issue #14.

Dear @mmolero , Thanks for your answer. I am sure that simndt can be used in asonic ultrasound simulation. However, I am still confused that how to extract the pressure wavefield u in the whole computational domain not from the outputs of the receivers . e.g., if the numerical grid is 10001000, how to get the pressure output with size 10001000. I have met the issure for a long time. Very sorry to disturb you.

mmolero commented 1 year ago

You need to observe for example a set of points, imagine that you have a receiver that can be represented as x number of points so that the time response it can be seen as the average for all these points for each time step. In this way you can have a time transient signal coming from averaging the receiver line. Does it make sense to you?

Nicko222 commented 1 year ago

You need to observe for example a set of points, imagine that you have a receiver that can be represented as x number of points so that the time response it can be seen as the average for all these points for each time step. In this way you can have a time transient signal coming from averaging the receiver line. Does it make sense to you?

Dear @mmolero , Thanks for your explanation. If I catch your idea, if the numerical grid is 1000 by 1000, I could use 1000 receivers to save the signal, right??? In addition, I have attached a simulation figure. In the figure, the pressure wavefield is in the whole computational domain. It is the pressure I want to exact, right? If this, how to save the outputs in the whole compuatational domain? Thanks a lot! simndt1