|
[Sponsors] |
November 22, 2010, 15:15 |
Aeroacoustics in OpenFOAM...
|
#1 |
New Member
Bertrand Delvaux
Join Date: Apr 2010
Posts: 29
Rep Power: 16 |
Hello,
I am pretty new to the great OF Community. For my PhD, I need to study the acoustics coming out of a shaped tube (like blowing into an instrument). My geometry is the following : inlet, shaped tube, outlet consisting in a hemisphere (open atmosphere). I settled periodic boundary condition for velocity as the inlet and "open atmosphere" for the outlet . I solve the problem with rhoPisoFoam, LES. Up to now, I've found pretty nice results, using, among others this great tool audioToProbe, to check up with the acoustic signal at different places, but unfortunately, my signal gets noisier with the distance increasing from the source, right when it goes into the hemisphere. I suspect the outlet boundary condition. Up to now, I tried p: outlet { type totalPressure; p0 uniform 1e5; U U; phi phi; rho rho; psi none; gamma 1; value uniform 1e5; } and U: outlet { type inletOutlet; inletValue uniform (0 0 0); value uniform (0 0 0); } Or p with "waveTransmissive" But I don't know how to deal with this issue about "non-reflective boundary condition" when dealing with aeroacoustics. It would be great if someone could give me some hints on these issues : 1. How to set this "open atmosphere" boundary condition? 2. How to take into account the non-reflective boundary condition? 3. How to adapt U outlet condition in accordingly? 4. Is that a good solver for Aeroacoustics? Thank you very much for your answers... Bertrand |
|
November 23, 2010, 00:33 |
|
#2 |
Super Moderator
Takuya OSHIMA
Join Date: Mar 2009
Location: Niigata City, Japan
Posts: 518
Blog Entries: 1
Rep Power: 20 |
Hi Bertrand,
I can't say much about compressible aeroacoustics itself, but one thing to note about probeToAudio is that it normalizes the input acoustic waveform with its maximum amplitude to avoid saturation and to obtain maximum S/N ratio. However there could be several potential problems by this feature under your scenario, two of which I can think of are:
Takuya |
|
November 25, 2010, 05:15 |
|
#3 |
New Member
Bertrand Delvaux
Join Date: Apr 2010
Posts: 29
Rep Power: 16 |
Hi Takuya,
Thank you very much for your answer! Your insights gave me some ideas... 1. How can I (simply?) modify probeToAudio in order not to normalize the sampled pressure with the maximum amplitude? 2. Concerning your suggestion 2, I have an input of a sinus (or half a sinus) for the velocity, with a fundamental frequency of 100 Hz. So there is no problem about being out of the audible range. Whenever I am "in the tube" it is a nice warm sound of 100 Hz, but as soon as I am out of the tube, in the open-atmosphere, the sound is still 100 Hz fundamental frequency (hopefully!) but gets noisier. 3. Therefore, when I add the visualization of what happens about p and U in this outlet hemisphere, I am very suspicious about the outlet boundary condition for both p and U : I really wonder how to implement such an "open atmospehre"... Would you have any idea? Regards, Bertrand |
|
November 27, 2010, 09:20 |
|
#4 |
Super Moderator
Takuya OSHIMA
Join Date: Mar 2009
Location: Niigata City, Japan
Posts: 518
Blog Entries: 1
Rep Power: 20 |
Hi Bertrand,
@1., @2: I haven't tried by myself but you should be able to modify the following bits to replace max(bias - minValue, maxValue - bias) with a fixed value. Code:
// scale and convert the data to be fit into the int16_t range scalar scaleFactor = max(bias - minValue, maxValue - bias); Takuya |
|
December 3, 2010, 06:57 |
|
#5 |
New Member
Bertrand Delvaux
Join Date: Apr 2010
Posts: 29
Rep Power: 16 |
Hi Takuya,
Thank you for your answer, I am now using Audacity, as you suggested months ago, in order to plot the spectrum : really helpfull in my case. 1. I have been modifying the code for audioToProbe, so that it doesn't normalize it, but how shall I "compile" this new code, (and from which location), in order to use it as an alternative audioToProbe utility? Sorry if my question may sound like a newbie one, but I am pretty new to OF, and still didn't implement my own functions, but I am looking forward to it. 2. I am still busy trying to figure out a way of implementing a far-field condition such as it would be an open atmosphere, non-reflective boundary condition, and wave-transmissive. As I saw in other posts, it looks like you are busy with this as well, so if you have more hints on that, I would definitely enjoy if you could share... 3. I've downloaded your last contribution to OF, still have to give it a look, but the principle looks amazing! Congrats' Regards, Bertrand |
|
December 3, 2010, 23:32 |
|
#6 |
Super Moderator
Takuya OSHIMA
Join Date: Mar 2009
Location: Niigata City, Japan
Posts: 518
Blog Entries: 1
Rep Power: 20 |
Hi Bertrand,
@Compiling: take a look at the User Guide first, http://www.openfoam.com/docs/user/co...p#x10-690003.2 @farfield condition: I know a convective outlet BC works well for low-Mach incompressible cases (which I studied with my own Fortran code) but not sure about compressible cases thus can't be of much help. As far as I've read in this forum the waveTransmissive BC should work? Takuya |
|
December 7, 2010, 06:59 |
|
#7 |
New Member
Bertrand Delvaux
Join Date: Apr 2010
Posts: 29
Rep Power: 16 |
Hi Takuya,
@Compiling : it didn't work, because there were some files missing in the official link to probeToAudio, but I could use the one from potentialWaveFoam, so no problem. @probeToAudio : I have a question : what is the meaning of 32767.0 in line 209? @convectiveOutlet : I am gonna try now, and do comparison cases with waveTransmissive. Bertrand |
|
Tags |
aeroacoustics, open atmosphere, rhopisofoam, wavetransmissive |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Superlinear speedup in OpenFOAM 13 | msrinath80 | OpenFOAM Running, Solving & CFD | 18 | March 3, 2015 06:36 |
Cross-compiling OpenFOAM 1.7.0 on Linux for Windows 32 and 64bits with Mingw-w64 | wyldckat | OpenFOAM Announcements from Other Sources | 3 | September 8, 2010 07:25 |
Modified OpenFOAM Forum Structure and New Mailing-List | pete | Site News & Announcements | 0 | June 29, 2009 06:56 |
64bitrhel5 OF installation instructions | mirko | OpenFOAM Installation | 2 | August 12, 2008 19:07 |
OpenFOAM Debian packaging current status problems and TODOs | oseen | OpenFOAM Installation | 9 | August 26, 2007 14:50 |