Mercurial > hg > propeller-model

annotate freqCalcCyl7HiFi.pd @ 0:cef74e8afc39

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Propeller Blog Patch Added
author Rod Selfridge <r.selfridge@qmul.ac.uk>
date Wed, 05 Apr 2017 08:31:12 +0100
parents
children
rev   line source
r@0 1 #N canvas 72 88 655 381 10;
r@0 2 #X obj 52 169 *~;
r@0 3 #X text 16 281 Calculates the frequency;
r@0 4 #X text 17 293 based on the physics of;
r@0 5 #X text 280 299 Calculates the gain of the dipole source;
r@0 6 #X text 392 216 Azimuth;
r@0 7 #X text 371 145 Elevation;
r@0 8 #X text 354 114 Distance;
r@0 9 #X text 318 54 Diameter;
r@0 10 #N canvas 322 154 491 423 Dipoles 0;
r@0 11 #X obj 210 217 *~;
r@0 12 #X obj 249 315 throw~ l;
r@0 13 #X obj 276 292 throw~ r;
r@0 14 #X obj 249 271 fcpan;
r@0 15 #X text 101 35 Generates the sound associated with the lift dipole.
r@0 16 ;
r@0 17 #X obj 210 108 vcf~ 10;
r@0 18 #X obj 210 131 vcf~ 10;
r@0 19 #X obj 284 95 r \$0-Qval;
r@0 20 #X floatatom 294 270 5 0 0 0 - - -;
r@0 21 #X obj 362 158 r myMet;
r@0 22 #X obj 362 178 snapshot~;
r@0 23 #X floatatom 362 199 0 0 0 0 - - -;
r@0 24 #X obj 210 67 inlet~;
r@0 25 #X obj 261 155 r broBoost;
r@0 26 #X obj 210 176 *~ 0.0001;
r@0 27 #X obj 225 196 inlet~;
r@0 28 #X obj 131 67 inlet~;
r@0 29 #N canvas 304 98 651 476 rangeLimit 0;
r@0 30 #X obj 285 58 inlet~;
r@0 31 #X obj 128 282 outlet~;
r@0 32 #X obj 103 60 sig~ 1;
r@0 33 #X obj 117 145 exp~;
r@0 34 #X obj 87 186 /~;
r@0 35 #X obj 102 166 +~;
r@0 36 #X obj 118 104 -~ 1;
r@0 37 #X obj 117 125 *~ 100;
r@0 38 #X obj 87 207 *~ -1;
r@0 39 #X obj 87 228 +~ 1;
r@0 40 #X obj 163 87 sig~ 1;
r@0 41 #X obj 173 173 exp~;
r@0 42 #X obj 143 214 /~;
r@0 43 #X obj 158 194 +~;
r@0 44 #X obj 173 130 -~ 1;
r@0 45 #X obj 173 151 *~ 1000;
r@0 46 #X obj 128 259 *~;
r@0 47 #X obj 173 108 /~ 15000;
r@0 48 #X obj 372 168 r myMet;
r@0 49 #X obj 372 188 snapshot~;
r@0 50 #X floatatom 372 209 0 0 0 0 - - -;
r@0 51 #X obj 118 81 /~ 10;
r@0 52 #X connect 0 0 17 0;
r@0 53 #X connect 0 0 21 0;
r@0 54 #X connect 2 0 4 0;
r@0 55 #X connect 2 0 5 0;
r@0 56 #X connect 3 0 5 1;
r@0 57 #X connect 4 0 8 0;
r@0 58 #X connect 5 0 4 1;
r@0 59 #X connect 6 0 7 0;
r@0 60 #X connect 7 0 3 0;
r@0 61 #X connect 8 0 9 0;
r@0 62 #X connect 9 0 16 0;
r@0 63 #X connect 10 0 12 0;
r@0 64 #X connect 10 0 13 0;
r@0 65 #X connect 11 0 13 1;
r@0 66 #X connect 12 0 16 1;
r@0 67 #X connect 13 0 12 1;
r@0 68 #X connect 14 0 15 0;
r@0 69 #X connect 15 0 11 0;
r@0 70 #X connect 16 0 1 0;
r@0 71 #X connect 17 0 14 0;
r@0 72 #X connect 18 0 19 0;
r@0 73 #X connect 19 0 20 0;
r@0 74 #X connect 21 0 6 0;
r@0 75 #X restore 99 156 pd rangeLimit;
r@0 76 #X obj 133 205 *~;
r@0 77 #X obj 276 249 r centrePan;
r@0 78 #X obj 249 227 inlet;
r@0 79 #X obj 210 248 spigot~;
r@0 80 #X connect 0 0 21 0;
r@0 81 #X connect 3 0 1 0;
r@0 82 #X connect 3 1 2 0;
r@0 83 #X connect 5 0 6 0;
r@0 84 #X connect 6 0 14 0;
r@0 85 #X connect 7 0 5 2;
r@0 86 #X connect 7 0 6 2;
r@0 87 #X connect 9 0 10 0;
r@0 88 #X connect 10 0 11 0;
r@0 89 #X connect 12 0 5 0;
r@0 90 #X connect 13 0 14 1;
r@0 91 #X connect 14 0 0 0;
r@0 92 #X connect 15 0 18 1;
r@0 93 #X connect 15 0 0 1;
r@0 94 #X connect 16 0 5 1;
r@0 95 #X connect 16 0 6 1;
r@0 96 #X connect 16 0 17 0;
r@0 97 #X connect 17 0 18 0;
r@0 98 #X connect 19 0 3 1;
r@0 99 #X connect 19 0 8 0;
r@0 100 #X connect 20 0 21 1;
r@0 101 #X connect 21 1 3 0;
r@0 102 #X connect 21 1 10 0;
r@0 103 #X restore 220 333 pd Dipoles;
r@0 104 #X obj 239 307 noise~;
r@0 105 #X text 17 306 the object and environment;
r@0 106 #X text 279 8 Source Speed;
r@0 107 #X obj 310 120 inlet~;
r@0 108 #X obj 258 35 inlet~;
r@0 109 #X obj 292 94 /~ 343.56;
r@0 110 #N canvas 110 127 1340 718 GoldfreqGain 0;
r@0 111 #N canvas 759 125 172 246 pow6 0;
r@0 112 #X obj 47 67 *~;
r@0 113 #X obj 76 67 *~;
r@0 114 #X obj 76 35 inlet~;
r@0 115 #X obj 47 163 outlet~;
r@0 116 #X obj 106 67 *~;
r@0 117 #X obj 47 93 *~;
r@0 118 #X obj 47 118 *~;
r@0 119 #X connect 0 0 5 0;
r@0 120 #X connect 1 0 5 1;
r@0 121 #X connect 2 0 0 0;
r@0 122 #X connect 2 0 0 1;
r@0 123 #X connect 2 0 1 1;
r@0 124 #X connect 2 0 4 0;
r@0 125 #X connect 2 0 4 1;
r@0 126 #X connect 2 0 1 0;
r@0 127 #X connect 4 0 6 1;
r@0 128 #X connect 5 0 6 0;
r@0 129 #X connect 6 0 3 0;
r@0 130 #X restore 218 82 pd pow6;
r@0 131 #X obj 218 123 *~;
r@0 132 #X text 587 13 Distance;
r@0 133 #X text 48 531 Calculates Lift Dipole gain;
r@0 134 #X obj 218 60 inlet~;
r@0 135 #X text 292 166 Diameter;
r@0 136 #X obj 219 618 outlet~;
r@0 137 #X text 298 252 Length;
r@0 138 #X obj 234 469 /~ 360;
r@0 139 #X obj 234 518 *~;
r@0 140 #X obj 219 541 *~;
r@0 141 #X text 74 100 Strouhal number;
r@0 142 #X text 216 46 Wind speed;
r@0 143 #X text 281 146 Corr Length;
r@0 144 #X text -86 235 Mass Density of Air;
r@0 145 #X text 108 242 Sqrt(2xpi);
r@0 146 #X obj 332 470 /~ 360;
r@0 147 #X obj 332 516 *~;
r@0 148 #X obj 332 492 cos~;
r@0 149 #X obj 219 567 *~;
r@0 150 #X obj 442 257 cos~;
r@0 151 #X obj 442 300 *~;
r@0 152 #X obj 427 196 sig~ 1;
r@0 153 #X obj 427 320 -~;
r@0 154 #N canvas 739 248 286 289 pow4 0;
r@0 155 #X obj 47 67 *~;
r@0 156 #X obj 76 67 *~;
r@0 157 #X obj 76 35 inlet~;
r@0 158 #X obj 47 163 outlet~;
r@0 159 #X obj 47 93 *~;
r@0 160 #X connect 0 0 4 0;
r@0 161 #X connect 1 0 4 1;
r@0 162 #X connect 2 0 0 0;
r@0 163 #X connect 2 0 0 1;
r@0 164 #X connect 2 0 1 1;
r@0 165 #X connect 2 0 1 0;
r@0 166 #X connect 4 0 3 0;
r@0 167 #X restore 427 340 pd pow4;
r@0 168 #X obj 412 360 *~;
r@0 169 #X obj 219 385 /~;
r@0 170 #X text 763 12 Azimuth;
r@0 171 #X text 662 14 Elevation;
r@0 172 #X text 512 235 Remember - Goldstein moves cylinder;
r@0 173 #X obj 442 238 /~ 360;
r@0 174 #X obj 219 210 *~;
r@0 175 #X obj 234 495 sin~;
r@0 176 #X obj 567 29 inlet~;
r@0 177 #X obj 665 32 inlet~;
r@0 178 #X obj 767 30 inlet~;
r@0 179 #X obj 411 155 *~;
r@0 180 #X text -214 358 If adding additional signal at 0 degrees then you
r@0 181 will have to add it before you divid by distance - for all sources!
r@0 182 ;
r@0 183 #X obj 219 597 *~ 100000;
r@0 184 #X obj 234 146 sig~ 8;
r@0 185 #X text -154 219 Scale factor to bring in line with CFD calcs (divided
r@0 186 by 400);
r@0 187 #X obj 219 233 *~ 0.0076766;
r@0 188 #X text 340 105 Speed of Sound to power 3 x 32;
r@0 189 #X obj 411 175 *~ 1.29764e+09;
r@0 190 #X obj 233 102 sig~ 0.85;
r@0 191 #X obj 457 278 inlet~;
r@0 192 #X obj 234 186 *~ 1;
r@0 193 #X obj 249 165 inlet;
r@0 194 #X obj 234 254 r length;
r@0 195 #X obj 219 276 *~ 1;
r@0 196 #X obj 556 271 r myMet;
r@0 197 #X obj 556 291 snapshot~;
r@0 198 #X floatatom 556 311 0 0 0 0 - - -;
r@0 199 #X connect 0 0 1 0;
r@0 200 #X connect 1 0 31 0;
r@0 201 #X connect 4 0 0 0;
r@0 202 #X connect 8 0 32 0;
r@0 203 #X connect 9 0 10 1;
r@0 204 #X connect 10 0 19 0;
r@0 205 #X connect 16 0 18 0;
r@0 206 #X connect 17 0 19 1;
r@0 207 #X connect 18 0 17 0;
r@0 208 #X connect 18 0 17 1;
r@0 209 #X connect 19 0 38 0;
r@0 210 #X connect 20 0 21 0;
r@0 211 #X connect 21 0 23 1;
r@0 212 #X connect 22 0 23 0;
r@0 213 #X connect 23 0 24 0;
r@0 214 #X connect 24 0 25 1;
r@0 215 #X connect 25 0 26 1;
r@0 216 #X connect 26 0 10 0;
r@0 217 #X connect 30 0 20 0;
r@0 218 #X connect 31 0 41 0;
r@0 219 #X connect 32 0 9 0;
r@0 220 #X connect 32 0 9 1;
r@0 221 #X connect 33 0 36 0;
r@0 222 #X connect 33 0 36 1;
r@0 223 #X connect 34 0 30 0;
r@0 224 #X connect 34 0 8 0;
r@0 225 #X connect 35 0 16 0;
r@0 226 #X connect 36 0 43 0;
r@0 227 #X connect 38 0 6 0;
r@0 228 #X connect 39 0 46 0;
r@0 229 #X connect 41 0 49 0;
r@0 230 #X connect 43 0 25 0;
r@0 231 #X connect 44 0 1 1;
r@0 232 #X connect 45 0 21 1;
r@0 233 #X connect 46 0 31 1;
r@0 234 #X connect 47 0 46 1;
r@0 235 #X connect 48 0 49 1;
r@0 236 #X connect 49 0 26 0;
r@0 237 #X connect 50 0 51 0;
r@0 238 #X connect 51 0 52 0;
r@0 239 #X restore 258 281 pd GoldfreqGain;
r@0 240 #X obj 52 123 sig~ 0.85;
r@0 241 #N canvas 161 91 608 545 doppler 0;
r@0 242 #X obj 91 27 inlet~;
r@0 243 #X obj 91 266 outlet~;
r@0 244 #X obj 121 47 sig~ 323.56;
r@0 245 #X obj 121 190 +~;
r@0 246 #X obj 106 212 /~;
r@0 247 #X obj 91 234 *~;
r@0 248 #X obj 136 90 delta~;
r@0 249 #X obj 136 111 *~ 44100;
r@0 250 #X obj 333 191 r myMet;
r@0 251 #X obj 333 211 snapshot~;
r@0 252 #X floatatom 333 232 0 0 0 0 - - -;
r@0 253 #X obj 136 69 inlet~;
r@0 254 #X obj 136 141 lop~ 1;
r@0 255 #X obj 136 163 lop~ 1;
r@0 256 #X connect 0 0 5 0;
r@0 257 #X connect 2 0 3 0;
r@0 258 #X connect 2 0 4 0;
r@0 259 #X connect 3 0 4 1;
r@0 260 #X connect 4 0 5 1;
r@0 261 #X connect 5 0 1 0;
r@0 262 #X connect 6 0 7 0;
r@0 263 #X connect 7 0 12 0;
r@0 264 #X connect 8 0 9 0;
r@0 265 #X connect 9 0 10 0;
r@0 266 #X connect 11 0 6 0;
r@0 267 #X connect 12 0 13 0;
r@0 268 #X connect 13 0 3 1;
r@0 269 #X restore 52 228 pd doppler;
r@0 270 #N canvas 290 200 365 347 scatter 0;
r@0 271 #X obj 138 20 inlet~;
r@0 272 #X obj 125 272 outlet~;
r@0 273 #X obj 177 53 snapshot~;
r@0 274 #X obj 177 73 >= 0;
r@0 275 #X obj 138 94 spigot~;
r@0 276 #X obj 178 141 spigot~;
r@0 277 #X obj 125 163 spigot~;
r@0 278 #X obj 182 21 r myMet;
r@0 279 #X text -5 301 Limits the value of the elevation to 25 to replicate
r@0 280 the scatter that would occur in the propeller when on a plane.;
r@0 281 #X obj 125 140 min~ -5;
r@0 282 #X obj 178 115 max~ 5;
r@0 283 #X connect 0 0 2 0;
r@0 284 #X connect 0 0 4 0;
r@0 285 #X connect 2 0 3 0;
r@0 286 #X connect 3 0 4 1;
r@0 287 #X connect 3 0 5 1;
r@0 288 #X connect 3 0 6 1;
r@0 289 #X connect 4 0 9 0;
r@0 290 #X connect 4 1 10 0;
r@0 291 #X connect 5 1 1 0;
r@0 292 #X connect 6 0 1 0;
r@0 293 #X connect 7 0 2 0;
r@0 294 #X connect 9 0 6 0;
r@0 295 #X connect 10 0 5 0;
r@0 296 #X restore 338 170 pd scatter;
r@0 297 #X obj 476 72 *~ 1.225;
r@0 298 #X obj 476 118 /~ 1.81e-05;
r@0 299 #X text 496 137 Reynolds Number;
r@0 300 #X text 550 113 Viscosity;
r@0 301 #N canvas 308 58 1027 635 StBandwidth 0;
r@0 302 #X obj 460 365 swap 1;
r@0 303 #X obj 460 391 /;
r@0 304 #X obj 461 434 outlet;
r@0 305 #X text 437 453 Q value;
r@0 306 #X obj 477 -36 r myMet;
r@0 307 #X obj 477 -14 snapshot~;
r@0 308 #X obj 203 164 spigot;
r@0 309 #X obj 219 187 spigot;
r@0 310 #X obj 236 140 > 0;
r@0 311 #X obj 219 233 + 0.9797;
r@0 312 #X obj 219 208 * 4.624e-05;
r@0 313 #X obj 602 150 spigot;
r@0 314 #X obj 618 173 spigot;
r@0 315 #X obj 651 149 < 1e+06;
r@0 316 #X obj 618 216 * 1.227e-10;
r@0 317 #X obj 618 196 *;
r@0 318 #X obj 690 215 * -8.553e-05;
r@0 319 #X obj 618 261 + 16.5;
r@0 320 #X obj 618 238 +;
r@0 321 #X text -84 298 First section of the Reynolds number - DeltaF/F relationship
r@0 322 is linear and did not appear to be reliant on current TU values;
r@0 323 #X text 565 329 Third section of the Reynolds number - DeltaF/F relationship
r@0 324 is linear and did not appear to be reliant on current TU values;
r@0 325 #X text 114 -16 Middle section of Reynolds number - DeltaF/F relationship
r@0 326 is also linear but depends on TU Value;
r@0 327 #X obj 461 412 max 0.005;
r@0 328 #X obj 635 126 > 193260;
r@0 329 #X obj 252 163 < 193260;
r@0 330 #X floatatom 352 446 0 0 0 0 - - -;
r@0 331 #X floatatom 588 23 0 0 0 0 - - -;
r@0 332 #X obj 460 342 / 100;
r@0 333 #X obj 531 -36 inlet~;
r@0 334 #X msg 593 195 10;
r@0 335 #X connect 0 0 1 0;
r@0 336 #X connect 0 1 1 1;
r@0 337 #X connect 1 0 22 0;
r@0 338 #X connect 4 0 5 0;
r@0 339 #X connect 5 0 13 0;
r@0 340 #X connect 5 0 11 0;
r@0 341 #X connect 5 0 8 0;
r@0 342 #X connect 5 0 6 0;
r@0 343 #X connect 5 0 23 0;
r@0 344 #X connect 5 0 24 0;
r@0 345 #X connect 5 0 26 0;
r@0 346 #X connect 6 0 7 0;
r@0 347 #X connect 7 0 10 0;
r@0 348 #X connect 8 0 6 1;
r@0 349 #X connect 9 0 27 0;
r@0 350 #X connect 10 0 9 0;
r@0 351 #X connect 11 0 12 0;
r@0 352 #X connect 11 0 29 0;
r@0 353 #X connect 12 0 16 0;
r@0 354 #X connect 12 0 15 1;
r@0 355 #X connect 12 0 15 0;
r@0 356 #X connect 13 0 12 1;
r@0 357 #X connect 14 0 18 0;
r@0 358 #X connect 15 0 14 0;
r@0 359 #X connect 16 0 18 1;
r@0 360 #X connect 18 0 17 0;
r@0 361 #X connect 22 0 2 0;
r@0 362 #X connect 22 0 25 0;
r@0 363 #X connect 23 0 11 1;
r@0 364 #X connect 24 0 7 1;
r@0 365 #X connect 27 0 0 0;
r@0 366 #X connect 28 0 5 0;
r@0 367 #X connect 29 0 27 0;
r@0 368 #X restore 476 154 pd StBandwidth;
r@0 369 #X obj 476 177 s \$0-Qval;
r@0 370 #X floatatom 551 177 0 0 0 0 - - -;
r@0 371 #X obj 327 146 inlet~;
r@0 372 #X obj 345 213 inlet~;
r@0 373 #N canvas 290 200 365 347 scatter 0;
r@0 374 #X obj 138 20 inlet~;
r@0 375 #X obj 125 272 outlet~;
r@0 376 #X obj 177 53 snapshot~;
r@0 377 #X obj 177 73 >= 0;
r@0 378 #X obj 138 94 spigot~;
r@0 379 #X obj 178 141 spigot~;
r@0 380 #X obj 125 163 spigot~;
r@0 381 #X obj 182 21 r myMet;
r@0 382 #X text -5 301 Limits the value of the elevation to 25 to replicate
r@0 383 the scatter that would occur in the propeller when on a plane.;
r@0 384 #X obj 125 140 min~ -5;
r@0 385 #X obj 178 115 max~ 5;
r@0 386 #X connect 0 0 2 0;
r@0 387 #X connect 0 0 4 0;
r@0 388 #X connect 2 0 3 0;
r@0 389 #X connect 3 0 4 1;
r@0 390 #X connect 3 0 5 1;
r@0 391 #X connect 3 0 6 1;
r@0 392 #X connect 4 0 9 0;
r@0 393 #X connect 4 1 10 0;
r@0 394 #X connect 5 1 1 0;
r@0 395 #X connect 6 0 1 0;
r@0 396 #X connect 7 0 2 0;
r@0 397 #X connect 9 0 6 0;
r@0 398 #X connect 10 0 5 0;
r@0 399 #X restore 362 259 pd scatter;
r@0 400 #X obj 546 261 r myMet;
r@0 401 #X obj 546 281 snapshot~;
r@0 402 #X floatatom 546 301 0 0 0 0 - - -;
r@0 403 #X obj 275 65 inlet;
r@0 404 #X obj 476 96 *~ 1;
r@0 405 #X obj 52 193 /~ 1;
r@0 406 #X obj 464 233 inlet;
r@0 407 #X connect 0 0 34 0;
r@0 408 #X connect 9 0 8 1;
r@0 409 #X connect 12 0 15 3;
r@0 410 #X connect 12 0 17 1;
r@0 411 #X connect 13 0 14 0;
r@0 412 #X connect 13 0 15 0;
r@0 413 #X connect 13 0 0 1;
r@0 414 #X connect 13 0 19 0;
r@0 415 #X connect 14 0 15 2;
r@0 416 #X connect 15 0 8 2;
r@0 417 #X connect 16 0 0 0;
r@0 418 #X connect 17 0 8 0;
r@0 419 #X connect 19 0 33 0;
r@0 420 #X connect 20 0 23 0;
r@0 421 #X connect 23 0 24 0;
r@0 422 #X connect 23 0 25 0;
r@0 423 #X connect 26 0 18 0;
r@0 424 #X connect 26 0 15 4;
r@0 425 #X connect 27 0 15 5;
r@0 426 #X connect 29 0 30 0;
r@0 427 #X connect 30 0 31 0;
r@0 428 #X connect 32 0 33 1;
r@0 429 #X connect 32 0 15 1;
r@0 430 #X connect 32 0 34 1;
r@0 431 #X connect 33 0 20 0;
r@0 432 #X connect 34 0 17 0;
r@0 433 #X connect 34 0 30 0;
r@0 434 #X connect 35 0 8 3;