Knowledge (XXG)

Computational aeroacoustics

Source đź“ť

172:
integral methods. The acoustic sources have to be known from some different source (e.g. a Finite Element simulation of a moving mechanical system or a fluid dynamic CFD simulation of the sources in a moving medium). The integral is taken over all sources at the retarded time (source time), which is the time at that the source is sent out the signal, which arrives now at a given observer position. Common to all integral methods is, that they cannot account for changes in the speed of sound or the average flow speed between source and observer position as they use a theoretical solution of the wave equation. When applying Lighthill's theory to the Navier Stokes equations of Fluid mechanics, one obtains volumetric sources, whereas the other two analogies provide the far field information based on a surface integral. Acoustic analogies can be very efficient and fast, as the known solution of the wave equation is used. One far away observer takes as long as one very close observer. Common for the application of all analogies is the integration over a large number of contributions, which can lead to additional numerical problems (addition/subtraction of many large numbers with result close to zero.) Furthermore, when applying an integral method, usually the source domain is limited somehow. While in theory the sources outside have to be zero, the application can not always fulfill this condition. Especially in connection with CFD simulations, this leads to large cut-off errors. By damping the source gradually to zero at the exit of the domain or adding some additional terms to correct this end-effect, these cut-off errors can be minimized.
212:
applies, can be produced as a superposition of monopoles and dipoles on the surface. The theory follows directly from the wave equation. The source strength of monopoles and dipoles on the surface can be calculated if the normal velocity (for monopoles) and the pressure (for dipoles) on the surface are known respectively. A modification of the method allows even to calculate the pressure on the surface based on the normal velocity only. The normal velocity could be given by a FE-simulation of a moving structure for instance. However, the modification to avoid the acoustic pressure on the surface to be known leads to problems, when considering an enclosed volume at its resonant frequencies, which is a major issue of the implementations of their method. The Kirchhoff integral method finds for instance application in
188:
it. Part of the acoustic wave is removed by one source and a new wave is radiated to fix the different wave speed. This often leads very large volumes with strong sources. Several modifications to Lighthill's original theory have been proposed to account for the sound-flow interaction or other effects. To improve Lighthill's analogy different quantities inside the wave operator as well as different wave operators are considered by following analogies. All of them obtain modified source terms, which sometimes allow a more clear sight on the "real" sources. The acoustic analogies of Lilley, Pierce, Howe and Möhring are only some examples for aeroacoustic analogies based on Lighthill's ideas. All acoustic analogies require a volume integration over a source term.
184:'. To obtain Lighthill's aeroacoustic analogy the governing Navier-Stokes equations are rearranged. The left hand side is a wave operator, which is applied to the density perturbation or pressure perturbation respectively. The right hand side is identified as the acoustic sources in a fluid flow, then. As Lighthill's analogy follows directly from the Navier-Stokes equations without simplification, all sources are present. Some of the sources are then identified as turbulent or laminar noise. The far-field sound pressure is then given in terms of a volume integral over the domain containing the sound source. The source term always includes physical sources and such sources, which describe the propagation in an inhomogeneous medium. 104:. Later in the beginning 1990s the growing CAA community picked up the term and extensively used it for any kind of numerical method describing the noise radiation from an aeroacoustic source or the propagation of sound waves in an inhomogeneous flow field. Such numerical methods can be far field integration methods (e.g. FW-H) as well as direct numerical methods optimized for the solutions (e.g.) of a mathematical model describing the aerodynamic noise generation and/or propagation. With the rapid development of the computational resources this field has undergone spectacular progress during the last three decades. 216:(BEM). A non-zero flow velocity is accounted by considering a moving frame of reference with the outer flow speed, in which the acoustic wave propagation takes place. Repetitive applications of the method can account for obstacles. First the sound field on the surface of the obstacle is calculated and then the obstacle is introduced by adding sources on its surface to cancel the normal velocity on the surface of the obstacle. Variations of the average flow field (speed of sound, density and velocity) can be taken into account by a similar method (e.g. dual reciprocity BEM). 135:(CFD) tool and secondly an acoustic solver. The flow field is then used to calculate the acoustical sources. Both steady state (RANS, SNGR (Stochastic Noise Generation and Radiation), ...) and transient (DNS, LES, DES, URANS, ...) fluid field solutions can be used. These acoustical sources are provided to the second solver which calculates the acoustical propagation. Acoustic propagation can be calculated using one of the following methods: 739: 229:
Different from the Kirchhoff method, these sources follow directly from the Navier-Stokes equations through Lighthill's analogy. Sources outside the FW-H surface can be accounted by an additional volume integral over quadrupole sources following from the Lighthill Tensor. However, when considering the same assumptions as Kirchhoffs linear theory, the FW-H method equals the Kirchhoff method.
22: 428: 416: 191:
The major difficulty with the acoustic analogy, however, is that the sound source is not compact in supersonic flow. Errors could be encountered in calculating the sound field, unless the computational domain could be extended in the downstream direction beyond the location where the sound source has
187:
The wave operator of Lighthill's analogy is limited to constant flow conditions outside the source zone. No variation of density, speed of sound and Mach number is allowed. Different mean flow conditions are identified as strong sources with opposite sign by the analogy, once an acoustic wave passes
121:
describes both the flow field, and the aerodynamically generated acoustic field. Thus both may be solved for directly. This requires very high numerical resolution due to the large differences in the length scale present between the acoustic variables and the flow variables. It is computationally
971:
A Fourier pseudospectral time-domain method can be applied to wave propagation problems pertinent to computational aeroacoustics. The original algorithm of the Fourier pseudo spectral time domain method works for periodical problems without the interaction with physical boundaries. A slip wall
211:
showed, that the radiation of sound from a limited source region can be described by enclosing this source region by a control surface - the so-called Kirchhoff surface. Then the sound field inside or outside the surface, where no sources are allowed and the wave operator on the left hand side
171:
There are multiple methods, which are based on a known solution of the acoustic wave equation to compute the acoustic far field of a sound source. Because a general solution for wave propagation in the free space can be written as an integral over all sources, these solutions are summarized as
228:
and Hawkings is based on Lighthill's acoustic analogy. However, by some mathematical modifications under the assumption of a limited source region, which is enclosed by a control surface (FW-H surface), the volume integral is avoided. Surface integrals over monopole and dipole sources remain.
192:
completely decayed. Furthermore, an accurate account of the retarded time-effect requires keeping a long record of the time-history of the converged solutions of the sound source, which again represents a storage problem. For realistic problems, the required storage can reach the order of 1
734:{\displaystyle \mathbf {U} ={\begin{bmatrix}\rho \\u\\v\\p\\\end{bmatrix}}\ ,\ \mathbf {F} ={\begin{bmatrix}\rho _{0}u+\rho u_{0}\\u_{0}u+p/\rho _{0}\\u_{0}v\\u_{0}p+\gamma p_{0}u\\\end{bmatrix}}\ ,\ \mathbf {G} ={\begin{bmatrix}\rho _{0}v\\0\\p/\rho _{0}\\\gamma p_{0}v\\\end{bmatrix}},} 92:
the field of computational fluid mechanics has been advancing rapidly in the past few years and now offers the hope that "computational aeroacoustics," where noise is computed directly from a first principles determination of continuous velocity and vorticity fields, might be possible,
130:
In this approach the computational domain is split into different regions, such that the governing acoustic or flow field can be solved with different equations and numerical techniques. This would involve using two different numerical solvers, first a dedicated
1032:
Hardin, J. C. and Lamkin, S. L., "Computational aeroacoustics - Present status and future promise," IN: Aero- and hydro-acoustics; Proceedings of the Symposium, Ecully, France, July 3–6, 1985 (A87-13585 03-71). Berlin and New York, Springer-Verlag, 1986, p.
99:
Later in a publication 1986 the same authors introduced the abbreviation CAA. The term was initially used for a low Mach number approach (Expansion of the acoustic perturbation field about an incompressible flow) as it is described under
972:
boundary condition, combined with buffer zone technique to solve some non-periodical aeroacoustic problems has been proposed. Compared to other computational methods, pseudospectral method is preferred for its high-order accuracy.
324: 1152:
X. Huang and X. Zhang, "A Fourier Pseudospectral Method for Some Computational Aeroacoustics Problems", International Journal of Aeroacoustics, Vol 5, No 3, 2006. pp.279-294.
932: 954: 884: 262: 842: 762: 316: 289: 1125:
Howe, M. S., "Contributions to the theory of aerodynamic sound, with application to excess jet noise and the theory of the flute", J. Fluid Mech., 71:625-673, 1975
86:
The origin of computational aeroacoustics can only very likely be dated back to the middle of the 1980s, with a publication of Hardin and Lamkin who claimed, that
822: 802: 782: 411:{\displaystyle {\frac {\partial \mathbf {U} }{\partial t}}+{\frac {\partial \mathbf {F} }{\partial x}}+{\frac {\partial \mathbf {G} }{\partial y}}=\mathbf {S} } 963:
flows in compressible regimes, the acoustic propagation may be influenced by non-linearities and the LEE may no longer be the appropriate mathematical model.
43: 30: 1268: 1248: 1143:
X. X. Chen, X. Huang and X. Zhang, "Sound Radiation from a Bypass Duct with Bifurcations", AIAA Journal, Vol. 47, No. 2, 2009. pp.429-436.
1233: 1213: 991:
Refer to the paper "Acoustic Perturbation Equations Based on Flow Decomposition via Source Filtering" by R.Ewert and W.Schroder.
1116:
Pierce, A. D., "Wave equation for the sound in fluids with unsteady inhomogeneous flow", J. Acoust. Soc. Am., 87:2292-2299, 1990
1228: 1162:
Ewert, R.; Schröder, W. (July 2003). "Acoustic perturbation equations based on flow decomposition via source filtering".
956:
on the right-side represents distributed unsteady sources. The application of LEE can be found in engine noise studies.
132: 1068:
C. K. W. Tam, and J. C. Webb, "Dispersion-Relation-Preserving Finite Difference Schemes for Computational Acoustics",
101: 1134:
Mohring, W. A well posed acoustic analogy based on a moving acoustic medium. 2010, arXiv preprint arXiv:1009.3766.
35: 1023:
Hardin, J.C. and Lamkin, S. L., "Aeroacoustic Computation of Cylinder Wake Flow," AIAA Journal, 22(1):51-57, 1984
1055:
Ffowcs Williams, J. E., and Hawkings, D. L., "Sound Generated by Turbulence and Surfaces in Arbitrary Motion",
213: 118: 1258: 208: 1171: 889: 225: 937: 1263: 1253: 847: 240: 827: 1179: 204: 181: 747: 294: 267: 1005: 1175: 807: 787: 767: 1183: 1242: 1000: 67: 1199:
Lighthill, M. J., "A General Introduction to Aeroacoustics and Atmospheric Sounds",
1223: 960: 237:
Considering small disturbances superimposed on a uniform mean flow of density
75: 1218: 21: 193: 1042:
Ffowcs Williams, "The Noise from Turbulence Convected at High Speed",
1107:
Lilley, G. M., "On the noise from air jets", AGARD CP 131, 13.1-13.12
318:, the Euler equations for a two dimensional model is presented as: 71: 15: 1201:
ICASE Report 92-52, NASA Langley Research Centre, Hampton, VA
1229:
Computational Aeroacoustics at Technische Universität Berlin
1094:
Lighthill, M. J., "On Sound Generated Aerodynamically, ii",
1081:
Lighthill, M. J., "On Sound Generated Aerodynamically, i",
1219:
Computational Aeroacoustics at the Ecole Centrale de Lyon
1224:
Computational Aeroacoustics at the University of Leuven
655: 502: 445: 122:
very demanding and unsuitable for any commercial use.
1234:
A CAA lecture script of Technische Universität Berlin
940: 892: 850: 830: 810: 790: 770: 750: 431: 327: 297: 270: 243: 948: 926: 878: 836: 816: 796: 776: 756: 733: 410: 310: 283: 256: 113:Direct numerical simulation (DNS) approach to CAA 1057:Philosophical Transactions of the Royal Society 1044:Philosophical Transactions of the Royal Society 8: 941: 939: 912: 903: 897: 891: 870: 861: 855: 849: 829: 809: 789: 769: 749: 711: 694: 685: 662: 650: 642: 616: 597: 580: 566: 557: 542: 528: 509: 497: 489: 440: 432: 430: 403: 384: 378: 359: 353: 334: 328: 326: 302: 296: 275: 269: 248: 242: 46:of all important aspects of the article. 1016: 70:that aims to analyze the generation of 42:Please consider expanding the lead to 7: 980:Expansion about Incompressible Flow 1214:Examples in Aeroacoustics from NASA 824:are the acoustic field variables, 391: 381: 366: 356: 341: 331: 14: 78:flows through numerical methods. 1164:Journal of Computational Physics 1070:Journal of Computational Physics 988:Acoustic Perturbation Equations 942: 643: 490: 433: 404: 385: 360: 335: 20: 927:{\displaystyle c_{p}/c_{v}=1.4} 34:may be too short to adequately 1059:, Vol. A264, 1969, pp. 321-342 1046:, Vol. A255, 1963, pp. 496-503 44:provide an accessible overview 1: 1269:Computational fields of study 1184:10.1016/S0021-9991(03)00168-2 1072:, Vol. 107, 1993, pp. 262-281 1249:Computational fluid dynamics 1085:, Vol. 211, 1952, pp 564-587 949:{\displaystyle \mathbf {S} } 844:the ratio of specific heats 133:Computational fluid dynamics 879:{\displaystyle c_{p}/c_{v}} 64:Computational aeroacoustics 1285: 233:Linearized Euler Equations 224:The integration method of 1098:, Vol. 222, 1954, pp 1-32 257:{\displaystyle \rho _{0}} 886:, for air at 20 Â°C 214:Boundary element methods 837:{\displaystyle \gamma } 291:and velocity on x-axis 175: 950: 934:, and the source term 928: 880: 838: 818: 798: 778: 758: 735: 412: 312: 285: 258: 119:Navier-Stokes equation 97: 951: 929: 881: 839: 819: 799: 779: 759: 757:{\displaystyle \rho } 736: 413: 313: 311:{\displaystyle u_{0}} 286: 284:{\displaystyle p_{0}} 259: 88: 938: 890: 848: 828: 808: 788: 768: 748: 429: 325: 295: 268: 241: 1176:2003JCoPh.188..365E 176:Lighthill's analogy 142:Lighthill's analogy 946: 924: 876: 834: 814: 794: 774: 754: 731: 722: 627: 474: 408: 308: 281: 254: 200:Kirchhoff integral 145:Kirchhoff integral 1096:Proc. Roy. Soc. A 1083:Proc. Roy. Soc. A 817:{\displaystyle p} 797:{\displaystyle v} 777:{\displaystyle u} 641: 635: 488: 482: 398: 373: 348: 139:Integral methods 117:The compressible 61: 60: 1276: 1188: 1187: 1159: 1153: 1150: 1144: 1141: 1135: 1132: 1126: 1123: 1117: 1114: 1108: 1105: 1099: 1092: 1086: 1079: 1073: 1066: 1060: 1053: 1047: 1040: 1034: 1030: 1024: 1021: 955: 953: 952: 947: 945: 933: 931: 930: 925: 917: 916: 907: 902: 901: 885: 883: 882: 877: 875: 874: 865: 860: 859: 843: 841: 840: 835: 823: 821: 820: 815: 803: 801: 800: 795: 783: 781: 780: 775: 763: 761: 760: 755: 740: 738: 737: 732: 727: 726: 716: 715: 699: 698: 689: 667: 666: 646: 639: 633: 632: 631: 621: 620: 602: 601: 585: 584: 571: 570: 561: 547: 546: 533: 532: 514: 513: 493: 486: 480: 479: 478: 436: 417: 415: 414: 409: 407: 399: 397: 389: 388: 379: 374: 372: 364: 363: 354: 349: 347: 339: 338: 329: 317: 315: 314: 309: 307: 306: 290: 288: 287: 282: 280: 279: 263: 261: 260: 255: 253: 252: 182:Acoustic Analogy 167:Integral methods 56: 53: 47: 24: 16: 1284: 1283: 1279: 1278: 1277: 1275: 1274: 1273: 1239: 1238: 1210: 1196: 1191: 1161: 1160: 1156: 1151: 1147: 1142: 1138: 1133: 1129: 1124: 1120: 1115: 1111: 1106: 1102: 1093: 1089: 1080: 1076: 1067: 1063: 1054: 1050: 1041: 1037: 1031: 1027: 1022: 1018: 1014: 1006:Acoustic theory 997: 986: 978: 969: 936: 935: 908: 893: 888: 887: 866: 851: 846: 845: 826: 825: 806: 805: 786: 785: 766: 765: 746: 745: 721: 720: 707: 701: 700: 690: 679: 678: 672: 671: 658: 651: 626: 625: 612: 593: 590: 589: 576: 573: 572: 562: 538: 535: 534: 524: 505: 498: 473: 472: 466: 465: 459: 458: 452: 451: 441: 427: 426: 390: 380: 365: 355: 340: 330: 323: 322: 298: 293: 292: 271: 266: 265: 244: 239: 238: 235: 226:Ffowcs Williams 222: 202: 178: 169: 128: 126:Hybrid approach 115: 110: 84: 66:is a branch of 57: 51: 48: 41: 29:This article's 25: 12: 11: 5: 1282: 1280: 1272: 1271: 1266: 1261: 1256: 1251: 1241: 1240: 1237: 1236: 1231: 1226: 1221: 1216: 1209: 1208:External links 1206: 1205: 1204: 1195: 1192: 1190: 1189: 1170:(2): 365–398. 1154: 1145: 1136: 1127: 1118: 1109: 1100: 1087: 1074: 1061: 1048: 1035: 1025: 1015: 1013: 1010: 1009: 1008: 1003: 996: 993: 985: 982: 977: 974: 968: 967:Pseudospectral 965: 944: 923: 920: 915: 911: 906: 900: 896: 873: 869: 864: 858: 854: 833: 813: 793: 773: 753: 742: 741: 730: 725: 719: 714: 710: 706: 703: 702: 697: 693: 688: 684: 681: 680: 677: 674: 673: 670: 665: 661: 657: 656: 654: 649: 645: 638: 630: 624: 619: 615: 611: 608: 605: 600: 596: 592: 591: 588: 583: 579: 575: 574: 569: 565: 560: 556: 553: 550: 545: 541: 537: 536: 531: 527: 523: 520: 517: 512: 508: 504: 503: 501: 496: 492: 485: 477: 471: 468: 467: 464: 461: 460: 457: 454: 453: 450: 447: 446: 444: 439: 435: 420: 419: 406: 402: 396: 393: 387: 383: 377: 371: 368: 362: 358: 352: 346: 343: 337: 333: 305: 301: 278: 274: 251: 247: 234: 231: 221: 218: 201: 198: 177: 174: 168: 165: 164: 163: 160: 157: 156:Pseudospectral 154: 151: 150: 149: 146: 143: 127: 124: 114: 111: 109: 106: 83: 80: 59: 58: 38:the key points 28: 26: 19: 13: 10: 9: 6: 4: 3: 2: 1281: 1270: 1267: 1265: 1262: 1260: 1257: 1255: 1252: 1250: 1247: 1246: 1244: 1235: 1232: 1230: 1227: 1225: 1222: 1220: 1217: 1215: 1212: 1211: 1207: 1202: 1198: 1197: 1193: 1185: 1181: 1177: 1173: 1169: 1165: 1158: 1155: 1149: 1146: 1140: 1137: 1131: 1128: 1122: 1119: 1113: 1110: 1104: 1101: 1097: 1091: 1088: 1084: 1078: 1075: 1071: 1065: 1062: 1058: 1052: 1049: 1045: 1039: 1036: 1029: 1026: 1020: 1017: 1011: 1007: 1004: 1002: 1001:Aeroacoustics 999: 998: 994: 992: 989: 983: 981: 975: 973: 966: 964: 962: 957: 921: 918: 913: 909: 904: 898: 894: 871: 867: 862: 856: 852: 831: 811: 791: 771: 751: 728: 723: 717: 712: 708: 704: 695: 691: 686: 682: 675: 668: 663: 659: 652: 647: 636: 628: 622: 617: 613: 609: 606: 603: 598: 594: 586: 581: 577: 567: 563: 558: 554: 551: 548: 543: 539: 529: 525: 521: 518: 515: 510: 506: 499: 494: 483: 475: 469: 462: 455: 448: 442: 437: 425: 424: 423: 400: 394: 375: 369: 350: 344: 321: 320: 319: 303: 299: 276: 272: 249: 245: 232: 230: 227: 219: 217: 215: 210: 206: 199: 197: 195: 189: 185: 183: 180:Also called ' 173: 166: 161: 158: 155: 152: 147: 144: 141: 140: 138: 137: 136: 134: 125: 123: 120: 112: 107: 105: 103: 96: 94: 87: 81: 79: 77: 73: 69: 68:aeroacoustics 65: 55: 52:February 2013 45: 39: 37: 32: 27: 23: 18: 17: 1259:Aerodynamics 1200: 1167: 1163: 1157: 1148: 1139: 1130: 1121: 1112: 1103: 1095: 1090: 1082: 1077: 1069: 1064: 1056: 1051: 1043: 1038: 1028: 1019: 990: 987: 979: 970: 958: 743: 421: 236: 223: 203: 190: 186: 179: 170: 129: 116: 98: 91: 89: 85: 63: 62: 49: 33: 31:lead section 961:Mach number 264:, pressure 1243:Categories 1012:References 1264:Mechanics 1254:Acoustics 959:For high 832:γ 752:ρ 705:γ 692:ρ 660:ρ 610:γ 564:ρ 522:ρ 507:ρ 449:ρ 392:∂ 382:∂ 367:∂ 357:∂ 342:∂ 332:∂ 246:ρ 209:Helmholtz 205:Kirchhoff 196:of data. 76:turbulent 36:summarize 1033:253-259. 995:See also 194:terabyte 1194:Sources 1172:Bibcode 108:Methods 82:History 1203:, 1992 744:where 640:  634:  487:  481:  422:where 72:noise 804:and 220:FW-H 207:and 148:FW-H 1180:doi 1168:188 984:APE 976:EIF 922:1.4 162:APE 159:EIF 153:LEE 102:EIF 74:by 1245:: 1178:. 1166:. 784:, 764:, 95:" 1186:. 1182:: 1174:: 943:S 919:= 914:v 910:c 905:/ 899:p 895:c 872:v 868:c 863:/ 857:p 853:c 812:p 792:v 772:u 729:, 724:] 718:v 713:0 709:p 696:0 687:/ 683:p 676:0 669:v 664:0 653:[ 648:= 644:G 637:, 629:] 623:u 618:0 614:p 607:+ 604:p 599:0 595:u 587:v 582:0 578:u 568:0 559:/ 555:p 552:+ 549:u 544:0 540:u 530:0 526:u 519:+ 516:u 511:0 500:[ 495:= 491:F 484:, 476:] 470:p 463:v 456:u 443:[ 438:= 434:U 418:, 405:S 401:= 395:y 386:G 376:+ 370:x 361:F 351:+ 345:t 336:U 304:0 300:u 277:0 273:p 250:0 90:" 54:) 50:( 40:.

Index


lead section
summarize
provide an accessible overview
aeroacoustics
noise
turbulent
EIF
Navier-Stokes equation
Computational fluid dynamics
Acoustic Analogy
terabyte
Kirchhoff
Helmholtz
Boundary element methods
Ffowcs Williams
Mach number
Aeroacoustics
Acoustic theory
Bibcode
2003JCoPh.188..365E
doi
10.1016/S0021-9991(03)00168-2
Examples in Aeroacoustics from NASA
Computational Aeroacoustics at the Ecole Centrale de Lyon
Computational Aeroacoustics at the University of Leuven
Computational Aeroacoustics at Technische Universität Berlin
A CAA lecture script of Technische Universität Berlin
Categories
Computational fluid dynamics

Text is available under the Creative Commons Attribution-ShareAlike License. Additional terms may apply.

↑