helmholtz's People
helmholtz's Issues
Milestone 1: Write the initial version of the manual/readme.md file
Milestone 1: Update readme.md for the new frequency-dependent stuff
The readme.md file has a section on verification that we compute the correct answer. This still uses the scheme from before #9, where density and wave speed are provided in the .prm file for all frequencies. The description of the inputs for these verification cases needs to be updated now that #9 is done.
Milestone 2: Fix the output format of the P and U matrices
From Jason's write-up:
Eventually for a system with 2 ports, Ares is going to need a matrix that looks like this:
| m11 m12 m13 m14 | | P1 | = |b1|
| m21 m22 m23 m24 | | U1 | |b2|
| P2 |
| U2 |
where the velocities U1
and U2
are defined to be into the node (or into the volume).
For a system with 3 ports, I'm going to need a matrix that looks like this:
| m11 m12 m13 m14 m15 m16| | P1 | = |b1|
| m21 m22 m23 m24 m25 m26| | U1 | |b2|
| m31 m32 m33 m34 m35 m36| | P2 | |b3|
| U2 |
| P3 |
| U3 |
etc.
For a volume with N ports, I need N linear equations relating the 2*N DOF's (P1-PN and U1-UN) to each other.
(Of course, you need 2*N equations to solve for the full system. The other N equations come from how the volume's ports are connected to the other components in the equivalent circuit model in Ares.)
For your problem, the b coefficients will all be 0 (i.e. there are no sources in your volume).
Also, we don't need all of the m coefficients. Because this is a linear system, we know that the 2N DOF's can be related to each other through a simple NxN system of equations that can look like:
|Q11 Q12 |* |P1| = |U1|
|Q21 Q22 | |P2| |U2|
which can be rewritten as
|Q11 -1 Q12 0| * |P1| = |0|
|Q21 0 Q22 -1| |U1| |0|
|P2|
|U2|
The form for the N=3 system will be
|Q11 -1 Q12 0 Q13 0| * |P1| = |0|
|Q21 0 Q22 -1 Q23 0| |U1| |0|
|Q31 0 Q32 0 Q33 -1| |P2| |0|
|U2|
|P3|
|U3|
So half of our mij coefficients are either -1 or 0.
To find the N*N Qij coefficients, I solve the system N times:
once with P1=1, P2-PN=0, and then I solve for U1-UN
once with P2=1, P1=P3-PN=0, and again I solve for U1-UN
once with P3=1, P1=P2=P4-PN=0, and again I solve for U1-UN
. . .
once with PN=1, P1-P(N-1)=0, and again I solve for U1-UN
Then from these N solutions, I'm able to solve for the N*N Qij coefficients. I'll go through the math of this for the N=2 case.
Let's rewrite
|Q11 Q12 | |P1| = |U1|
|Q21 Q22| |P2| |U2|
as
|P1 P1 0 0 | * |Q11| = |U1|
| 0 0 P1 P2 | |Q12| |U2|
|Q21|
|Q22|
So we've got N=2 equations and N*N unknown.
Now suppose we know P1, U1, P2, and U2 for N separate unique conditions. Let's call these 2 unique conditions:
P1_1, U1_1, P2_1, U2_1 for condition "_1"
P1_2, U1_2, P2_2, U2_2 for condition "_2"
|P1_1 P2_1 0 0 | * |Q11| = |U1_1|
| 0 0 P1_1 P2_1| |Q12| |U2_1|
|P1_2 P2_2 0 0 | |Q21| |U1_2|
| 0 0 P1_2 P2_2| |Q22| |U2_2|
So now we have N unknowns and N equations.
We can solve for Qij's and from there assign the mij coefficients.
For a N=3 system, we have NxN=3*3 = 9 unknowns, but we get 3 equations for every solution and the system to be solved is
|P1_1 P2_1 P2_1 0 0 0 0 0 0 | * |Q11| = |U1_1|
| 0 0 0 P1_1 P2_1 P3_1 0 0 0 | |Q12| |U2_1|
| 0 0 0 0 0 0 P1_1 P2_1 P3_1 | |Q13| |U3_1|
|P1_2 P2_2 P2_2 0 0 0 0 0 0 | |Q21| |U1_2|
| 0 0 0 P1_2 P2_2 P3_2 0 0 0 | |Q22| |U2_2|
| 0 0 0 0 0 0 P1_2 P2_2 P3_2 | |Q23| |U3_2|
|P1_3 P2_3 P2_3 0 0 0 0 0 0 | |Q31| |U1_3|
| 0 0 0 P1_3 P2_3 P3_3 0 0 0 | |Q32| |U2_3|
| 0 0 0 0 0 0 P1_3 P2_3 P3_3 | |Q33| |U3_3|
Also write output data into a CSV file
Milestone 2: Implement the mesh refinement criterion
As stated in the readme.md, we want to do this:
On the other hand, if the number of mesh refinement steps is negative,
then the mesh is adapted using the following algorithm. For a given
wave speed, we can compute the wavelength of oscillations for a
given frequency _w_ via
lambda = c/f
= c/(w/(2pi))
If `lambda` is larger than the diameter of the domain, then we replace
it by the diameter of the domain since that is the largest wavelength
at which the solution can vary.
We can then determine a frequency-adapted mesh size in the following
way: We want that there are at least _N_ grid points per wave
length. Since there are as many grid points per cell as the polynomial
degree _p_, this equations to requiring that the mesh size _Delta x_
satisfies _Delta x / p <= lambda/N_
or equivalently: _Delta x <= lambda/N * p_. The `Mesh refinement steps`
parameter is then interpreted as _N_ if it is negative.
In the code, we just need to compute lambda correctly -- right now it is simply set to 1000.
Milestone 1: Name output files by the number of the frequency
As $title says. This requires moving the number of the frequency through from the place where we call the function.
Milestone 1: Correct the definition of the velocity on ports during postprocessing
Concatenate all .csv files into one
Ask Jason whether we need this.
Finish documenting point values in the example section
Reverse the order of indices in the U, P, and M matrices
Columns should be source port, rows destination port. Currently, it's the other way around.
Milestone 1: Verify correctness of the code
Use a simple pipe with two ports, for which I can write down the exact solution of the problem.
Visualize boundary indicators
...so that we can actually see what geometry we are given.
Milestone 2: Read COMSOL files
Write a reader for the COMSOL .mphtxt file format:
https://www.comsol.com/fileformats
See also Jason's email from 11/2/2021.
Signal failure
<output prefix>solver_failure_signal.txt
This file is created by dealII_3D and indicates that there was an error in solving the problem passed to it. Information about the error will be placed into the output_log.txt file. ealII_3D will terminate after this file is created.
Move surface output into visualization/ folder
Related to #19.
Flag that determines that we only want diagnostic output, no solution
Document CSV output file
Output boundary indicators and areas
Could you report the surface areas and names of the interface surfaces?
It would be nice if this could be a separate command call to helmholtz.exe where helmholtz.exe loads the mesh file and just reports some information about the mesh. Currently the number, names and area of the interface areas are needed.
Milestone 1: Correct the boundary conditions on the output ports
Right now, we prescribe unit pressure at the source port and zero pressure at the other ports. For all other parts of the boundary, we prescribe Neumann values
nabla p = 0
I suspect, however, that setting the pressure to zero at the outlet ports is not what is really wanted. Need to discuss with Jason.
Milestone 1: Implement frequency-dependent density and bulk modulus
For example, for air:
% frequency(Hz) density(kg/m3) bulk modulus(Pa)
10 1.205728 142090.344491053
60.2010050251256 1.205728 142090.344491053
110.402010050251 1.205728 142090.344491053
160.603015075377 1.205728 142090.344491053
210.804020100503 1.205728 142090.344491053
261.005025125628 1.205728 142090.344491053
311.206030150754 1.205728 142090.344491053
361.407035175879 1.205728 142090.344491053
411.608040201005 1.205728 142090.344491053
461.809045226131 1.205728 142090.344491053
512.010050251256 1.205728 142090.344491053
562.211055276382 1.205728 142090.344491053
612.412060301507 1.205728 142090.344491053
662.613065326633 1.205728 142090.344491053
712.814070351759 1.205728 142090.344491053
763.015075376884 1.205728 142090.344491053
813.21608040201 1.205728 142090.344491053
863.417085427136 1.205728 142090.344491053
913.618090452261 1.205728 142090.344491053
963.819095477387 1.205728 142090.344491053
1014.02010050251 1.205728 142090.344491053
1064.22110552764 1.205728 142090.344491053
1114.42211055276 1.205728 142090.344491053
1164.62311557789 1.205728 142090.344491053
1214.82412060302 1.205728 142090.344491053
1265.02512562814 1.205728 142090.344491053
1315.22613065327 1.205728 142090.344491053
1365.42713567839 1.205728 142090.344491053
1415.62814070352 1.205728 142090.344491053
1465.82914572864 1.205728 142090.344491053
1516.03015075377 1.205728 142090.344491053
1566.23115577889 1.205728 142090.344491053
1616.43216080402 1.205728 142090.344491053
1666.63316582915 1.205728 142090.344491053
1716.83417085427 1.205728 142090.344491053
1767.0351758794 1.205728 142090.344491053
1817.23618090452 1.205728 142090.344491053
1867.43718592965 1.205728 142090.344491053
1917.63819095477 1.205728 142090.344491053
1967.8391959799 1.205728 142090.344491053
2018.04020100503 1.205728 142090.344491053
2068.24120603015 1.205728 142090.344491053
2118.44221105528 1.205728 142090.344491053
2168.6432160804 1.205728 142090.344491053
2218.84422110553 1.205728 142090.344491053
2269.04522613065 1.205728 142090.344491053
2319.24623115578 1.205728 142090.344491053
2369.4472361809 1.205728 142090.344491053
2419.64824120603 1.205728 142090.344491053
2469.84924623116 1.205728 142090.344491053
2520.05025125628 1.205728 142090.344491053
2570.25125628141 1.205728 142090.344491053
2620.45226130653 1.205728 142090.344491053
2670.65326633166 1.205728 142090.344491053
2720.85427135678 1.205728 142090.344491053
2771.05527638191 1.205728 142090.344491053
2821.25628140704 1.205728 142090.344491053
2871.45728643216 1.205728 142090.344491053
2921.65829145729 1.205728 142090.344491053
2971.85929648241 1.205728 142090.344491053
3022.06030150754 1.205728 142090.344491053
3072.26130653266 1.205728 142090.344491053
3122.46231155779 1.205728 142090.344491053
3172.66331658291 1.205728 142090.344491053
3222.86432160804 1.205728 142090.344491053
3273.06532663317 1.205728 142090.344491053
3323.26633165829 1.205728 142090.344491053
3373.46733668342 1.205728 142090.344491053
3423.66834170854 1.205728 142090.344491053
3473.86934673367 1.205728 142090.344491053
3524.07035175879 1.205728 142090.344491053
3574.27135678392 1.205728 142090.344491053
3624.47236180905 1.205728 142090.344491053
3674.67336683417 1.205728 142090.344491053
3724.8743718593 1.205728 142090.344491053
3775.07537688442 1.205728 142090.344491053
3825.27638190955 1.205728 142090.344491053
3875.47738693467 1.205728 142090.344491053
3925.6783919598 1.205728 142090.344491053
3975.87939698492 1.205728 142090.344491053
4026.08040201005 1.205728 142090.344491053
4076.28140703518 1.205728 142090.344491053
4126.4824120603 1.205728 142090.344491053
4176.68341708543 1.205728 142090.344491053
4226.88442211055 1.205728 142090.344491053
4277.08542713568 1.205728 142090.344491053
4327.2864321608 1.205728 142090.344491053
4377.48743718593 1.205728 142090.344491053
4427.68844221106 1.205728 142090.344491053
4477.88944723618 1.205728 142090.344491053
4528.09045226131 1.205728 142090.344491053
4578.29145728643 1.205728 142090.344491053
4628.49246231156 1.205728 142090.344491053
4678.69346733668 1.205728 142090.344491053
4728.89447236181 1.205728 142090.344491053
4779.09547738693 1.205728 142090.344491053
4829.29648241206 1.205728 142090.344491053
4879.49748743719 1.205728 142090.344491053
4929.69849246231 1.205728 142090.344491053
4979.89949748744 1.205728 142090.344491053
5030.10050251256 1.205728 142090.344491053
5080.30150753769 1.205728 142090.344491053
5130.50251256281 1.205728 142090.344491053
5180.70351758794 1.205728 142090.344491053
5230.90452261307 1.205728 142090.344491053
5281.10552763819 1.205728 142090.344491053
5331.30653266332 1.205728 142090.344491053
5381.50753768844 1.205728 142090.344491053
5431.70854271357 1.205728 142090.344491053
5481.90954773869 1.205728 142090.344491053
5532.11055276382 1.205728 142090.344491053
5582.31155778894 1.205728 142090.344491053
5632.51256281407 1.205728 142090.344491053
5682.7135678392 1.205728 142090.344491053
5732.91457286432 1.205728 142090.344491053
5783.11557788945 1.205728 142090.344491053
5833.31658291457 1.205728 142090.344491053
5883.5175879397 1.205728 142090.344491053
5933.71859296482 1.205728 142090.344491053
5983.91959798995 1.205728 142090.344491053
6034.12060301508 1.205728 142090.344491053
6084.3216080402 1.205728 142090.344491053
6134.52261306533 1.205728 142090.344491053
6184.72361809045 1.205728 142090.344491053
6234.92462311558 1.205728 142090.344491053
6285.1256281407 1.205728 142090.344491053
6335.32663316583 1.205728 142090.344491053
6385.52763819095 1.205728 142090.344491053
6435.72864321608 1.205728 142090.344491053
6485.92964824121 1.205728 142090.344491053
6536.13065326633 1.205728 142090.344491053
6586.33165829146 1.205728 142090.344491053
6636.53266331658 1.205728 142090.344491053
6686.73366834171 1.205728 142090.344491053
6736.93467336683 1.205728 142090.344491053
6787.13567839196 1.205728 142090.344491053
6837.33668341709 1.205728 142090.344491053
6887.53768844221 1.205728 142090.344491053
6937.73869346734 1.205728 142090.344491053
6987.93969849246 1.205728 142090.344491053
7038.14070351759 1.205728 142090.344491053
7088.34170854271 1.205728 142090.344491053
7138.54271356784 1.205728 142090.344491053
7188.74371859296 1.205728 142090.344491053
7238.94472361809 1.205728 142090.344491053
7289.14572864322 1.205728 142090.344491053
7339.34673366834 1.205728 142090.344491053
7389.54773869347 1.205728 142090.344491053
7439.74874371859 1.205728 142090.344491053
7489.94974874372 1.205728 142090.344491053
7540.15075376884 1.205728 142090.344491053
7590.35175879397 1.205728 142090.344491053
7640.5527638191 1.205728 142090.344491053
7690.75376884422 1.205728 142090.344491053
7740.95477386935 1.205728 142090.344491053
7791.15577889447 1.205728 142090.344491053
7841.3567839196 1.205728 142090.344491053
7891.55778894472 1.205728 142090.344491053
7941.75879396985 1.205728 142090.344491053
7991.95979899498 1.205728 142090.344491053
8042.1608040201 1.205728 142090.344491053
8092.36180904523 1.205728 142090.344491053
8142.56281407035 1.205728 142090.344491053
8192.76381909548 1.205728 142090.344491053
8242.9648241206 1.205728 142090.344491053
8293.16582914573 1.205728 142090.344491053
8343.36683417085 1.205728 142090.344491053
8393.56783919598 1.205728 142090.344491053
8443.76884422111 1.205728 142090.344491053
8493.96984924623 1.205728 142090.344491053
8544.17085427136 1.205728 142090.344491053
8594.37185929648 1.205728 142090.344491053
8644.57286432161 1.205728 142090.344491053
8694.77386934673 1.205728 142090.344491053
8744.97487437186 1.205728 142090.344491053
8795.17587939699 1.205728 142090.344491053
8845.37688442211 1.205728 142090.344491053
8895.57788944724 1.205728 142090.344491053
8945.77889447236 1.205728 142090.344491053
8995.97989949749 1.205728 142090.344491053
9046.18090452261 1.205728 142090.344491053
9096.38190954774 1.205728 142090.344491053
9146.58291457286 1.205728 142090.344491053
9196.78391959799 1.205728 142090.344491053
9246.98492462312 1.205728 142090.344491053
9297.18592964824 1.205728 142090.344491053
9347.38693467337 1.205728 142090.344491053
9397.58793969849 1.205728 142090.344491053
9447.78894472362 1.205728 142090.344491053
9497.98994974874 1.205728 142090.344491053
9548.19095477387 1.205728 142090.344491053
9598.391959799 1.205728 142090.344491053
9648.59296482412 1.205728 142090.344491053
9698.79396984925 1.205728 142090.344491053
9748.99497487437 1.205728 142090.344491053
9799.1959798995 1.205728 142090.344491053
9849.39698492462 1.205728 142090.344491053
9899.59798994975 1.205728 142090.344491053
9949.79899497487 1.205728 142090.344491053
10000 1.205728 142090.344491053
and for a porous metarial:
%frequency(Hz) density(kg/m3) bulk modulus(Pa)
10 1.7507 +i*-251.40 112615.491+i*119.669
60.201 1.7507 +i*-41.763 112631.842+i*720.021
110.402 1.7507 +i*-22.776 112671.498+i*1318.65
160.603 1.7506 +i*-15.66 112734.259+i*1914.15
210.804 1.7505 +i*-11.934 112819.809+i*2505.12
261.005 1.7504 +i*-9.6422 112927.724+i*3090.22
311.206 1.7502 +i*-8.0903 113057.471+i*3668.13
361.407 1.75 +i*-6.97 113208.422+i*4237.61
411.608 1.7498 +i*-6.1235 113379.853+i*4797.49
461.809 1.7496 +i*-5.4615 113570.962+i*5346.66
512.0101 1.7493 +i*-4.9296 113780.870+i*5884.10
562.2111 1.7491 +i*-4.493 114008.635+i*6408.88
612.4121 1.7487 +i*-4.1283 114253.261+i*6920.16
662.6131 1.7484 +i*-3.8192 114513.708+i*7417.19
712.8141 1.7481 +i*-3.5538 114788.905+i*7899.32
763.0151 1.7477 +i*-3.3236 115077.755+i*8366.00
813.2161 1.7473 +i*-3.1221 115379.149+i*8816.76
863.4171 1.7468 +i*-2.9441 115691.971+i*9251.25
913.6181 1.7464 +i*-2.786 116015.112+i*9669.18
963.8191 1.7459 +i*-2.6445 116347.473+i*10070.3
1014.0201 1.7454 +i*-2.5171 116687.976+i*10454.7
1064.2211 1.7449 +i*-2.4019 117035.568+i*10822.1
1114.4221 1.7443 +i*-2.2973 117389.226+i*11172.7
1164.6231 1.7437 +i*-2.2018 117747.966+i*11506.7
1214.8241 1.7432 +i*-2.1144 118110.842+i*11824.0
1265.0251 1.7425 +i*-2.034 118476.952+i*12125.1
1315.2261 1.7419 +i*-1.9598 118845.440+i*12410.1
1365.4271 1.7412 +i*-1.8913 119215.500+i*12679.4
1415.6281 1.7406 +i*-1.8277 119586.372+i*12933.3
1465.8291 1.7399 +i*-1.7685 119957.347+i*13172.3
1516.0302 1.7392 +i*-1.7134 120327.767+i*13396.7
1566.2312 1.7384 +i*-1.6619 120697.022+i*13607.0
1616.4322 1.7377 +i*-1.6137 121064.551+i*13803.7
1666.6332 1.7369 +i*-1.5685 121429.842+i*13987.3
1716.8342 1.7361 +i*-1.526 121792.427+i*14158.1
1767.0352 1.7353 +i*-1.4859 122151.885+i*14316.7
1817.2362 1.7345 +i*-1.4482 122507.835+i*14463.7
1867.4372 1.7336 +i*-1.4126 122859.939+i*14599.4
1917.6382 1.7328 +i*-1.3789 123207.898+i*14724.4
1967.8392 1.7319 +i*-1.3469 123551.447+i*14839.2
2018.0402 1.731 +i*-1.3166 123890.358+i*14944.2
2068.2412 1.7301 +i*-1.2879 124224.435+i*15040.0
2118.4422 1.7292 +i*-1.2605 124553.510+i*15126.9
2168.6432 1.7283 +i*-1.2345 124877.445+i*15205.4
2218.8442 1.7273 +i*-1.2097 125196.127+i*15275.9
2269.0452 1.7264 +i*-1.186 125509.466+i*15338.9
2319.2462 1.7254 +i*-1.1635 125817.397+i*15394.8
2369.4472 1.7244 +i*-1.1419 126119.870+i*15444.0
2419.6482 1.7234 +i*-1.1212 126416.859+i*15486.8
2469.8492 1.7224 +i*-1.1014 126708.349+i*15523.7
2520.0503 1.7214 +i*-1.0824 126994.343+i*15554.9
2570.2513 1.7204 +i*-1.0643 127274.858+i*15580.7
2620.4523 1.7194 +i*-1.0468 127549.922+i*15601.6
2670.6533 1.7183 +i*-1.03 127819.573+i*15617.9
2720.8543 1.7173 +i*-1.0139 128083.859+i*15629.7
2771.0553 1.7162 +i*-0.9983 128342.838+i*15637.4
2821.2563 1.7152 +i*-0.9834 128596.573+i*15641.3
2871.4573 1.7141 +i*-0.9690 128845.135+i*15641.6
2921.6583 1.713 +i*-0.9551 129088.599+i*15638.6
2971.8593 1.7119 +i*-0.9417 129327.048+i*15632.4
3022.0603 1.7108 +i*-0.9288 129560.564+i*15623.3
3072.2613 1.7097 +i*-0.9163 129789.237+i*15611.4
3122.4623 1.7086 +i*-0.9042 130013.157+i*15597.1
3172.6633 1.7075 +i*-0.8925 130232.416+i*15580.4
3222.8643 1.7064 +i*-0.8813 130447.11 +i*15561.5
3273.0653 1.7053 +i*-0.8703 130657.333+i*15540.7
3323.2663 1.7042 +i*-0.8597 130863.183+i*15518.0
3373.4673 1.7031 +i*-0.8495 131064.755+i*15493.5
3423.6683 1.702 +i*-0.8395 131262.148+i*15467.5
3473.8693 1.7008 +i*-0.8299 131455.457+i*15440.0
3524.0704 1.6997 +i*-0.8205 131644.779+i*15411.3
3574.2714 1.6986 +i*-0.8114 131830.209+i*15381.2
3624.4724 1.6974 +i*-0.8026 132011.842+i*15350.1
3674.6734 1.6963 +i*-0.7940 132189.771+i*15318.0
3724.8744 1.6952 +i*-0.7856 132364.090+i*15284.9
3775.0754 1.694 +i*-0.7775 132534.888+i*15251.0
3825.2764 1.6929 +i*-0.7696 132702.256+i*15216.4
3875.4774 1.6918 +i*-0.7619 132866.281+i*15181.0
3925.6784 1.6906 +i*-0.7544 133027.051+i*15145.1
3975.8794 1.6895 +i*-0.7471 133184.649+i*15108.6
4026.0804 1.6884 +i*-0.7400 133339.159+i*15071.7
4076.2814 1.6872 +i*-0.7331 133490.662+i*15034.3
4126.4824 1.6861 +i*-0.7263 133639.237+i*14996.6
4176.6834 1.685 +i*-0.7197 133784.963+i*14958.5
4226.8844 1.6838 +i*-0.7133 133927.914+i*14920.2
4277.0854 1.6827 +i*-0.7070 134068.166+i*14881.7
4327.2864 1.6816 +i*-0.7009 134205.789+i*14843.0
4377.4874 1.6805 +i*-0.6949 134340.855+i*14804.1
4427.6884 1.6793 +i*-0.6890 134473.431+i*14765.2
4477.8894 1.6782 +i*-0.6833 134603.585+i*14726.2
4528.0905 1.6771 +i*-0.6777 134731.382+i*14687.1
4578.2915 1.676 +i*-0.6722 134856.884+i*14648.0
4628.4925 1.6749 +i*-0.6669 134980.152+i*14608.9
4678.6935 1.6738 +i*-0.6616 135101.248+i*14569.9
4728.8945 1.6727 +i*-0.6565 135220.228+i*14530.9
4779.0955 1.6716 +i*-0.6515 135337.15 +i*14492.0
4829.2965 1.6705 +i*-0.6465 135452.067+i*14453.2
4879.4975 1.6694 +i*-0.6417 135565.033+i*14414.6
4929.6985 1.6683 +i*-0.6370 135676.100+i*14376.0
4979.8995 1.6672 +i*-0.6324 135785.318+i*14337.6
5030.1005 1.6661 +i*-0.6278 135892.735+i*14299.4
5080.3015 1.665 +i*-0.6234 135998.399+i*14261.3
5130.5025 1.664 +i*-0.6190 136102.356+i*14223.5
5180.7035 1.6629 +i*-0.6147 136204.650+i*14185.8
5230.9045 1.6618 +i*-0.6105 136305.324+i*14148.3
5281.1055 1.6608 +i*-0.6064 136404.420+i*14111.0
5331.3065 1.6597 +i*-0.6023 136501.979+i*14074.0
5381.5075 1.6587 +i*-0.5983 136598.04 +i*14037.2
5431.7085 1.6576 +i*-0.5944 136692.641+i*14000.6
5481.9095 1.6566 +i*-0.5906 136785.819+i*13964.2
5532.1106 1.6556 +i*-0.5868 136877.610+i*13928.1
5582.3116 1.6545 +i*-0.5831 136968.050+i*13892.3
5632.5126 1.6535 +i*-0.5794 137057.171+i*13856.7
5682.7136 1.6525 +i*-0.5759 137145.007+i*13821.3
5732.9146 1.6515 +i*-0.5723 137231.588+i*13786.2
5783.1156 1.6505 +i*-0.5689 137316.947+i*13751.4
5833.3166 1.6495 +i*-0.5655 137401.113+i*13716.8
5883.5176 1.6485 +i*-0.5621 137484.115+i*13682.5
5933.7186 1.6475 +i*-0.5588 137565.981+i*13648.5
5983.9196 1.6465 +i*-0.5555 137646.737+i*13614.7
6034.1206 1.6455 +i*-0.5523 137726.412+i*13581.2
6084.3216 1.6445 +i*-0.5492 137805.029+i*13547.9
6134.5226 1.6436 +i*-0.5461 137882.614+i*13514.9
6184.7236 1.6426 +i*-0.5430 137959.192+i*13482.2
6234.9246 1.6416 +i*-0.5400 138034.785+i*13449.7
6285.1256 1.6407 +i*-0.5371 138109.416+i*13417.5
6335.3266 1.6397 +i*-0.5341 138183.107+i*13385.5
6385.5276 1.6388 +i*-0.5313 138255.879+i*13353.8
6435.7286 1.6379 +i*-0.5284 138327.754+i*13322.4
6485.9296 1.6369 +i*-0.5256 138398.751+i*13291.2
6536.1307 1.636 +i*-0.5229 138468.89 +i*13260.3
6586.3317 1.6351 +i*-0.5201 138538.189+i*13229.6
6636.5327 1.6342 +i*-0.5175 138606.667+i*13199.2
6686.7337 1.6333 +i*-0.5148 138674.342+i*13169.0
6736.9347 1.6323 +i*-0.5122 138741.231+i*13139.1
6787.1357 1.6315 +i*-0.5096 138807.351+i*13109.4
6837.3367 1.6306 +i*-0.5071 138872.718+i*13080.0
6887.5377 1.6297 +i*-0.5046 138937.348+i*13050.8
6937.7387 1.6288 +i*-0.5021 139001.257+i*13021.8
6987.9397 1.6279 +i*-0.4997 139064.459+i*12993.1
7038.1407 1.627 +i*-0.4973 139126.968+i*12964.6
7088.3417 1.6262 +i*-0.4949 139188.800+i*12936.3
7138.5427 1.6253 +i*-0.4925 139249.968+i*12908.3
7188.7437 1.6245 +i*-0.4902 139310.484+i*12880.5
7238.9447 1.6236 +i*-0.4879 139370.362+i*12852.9
7289.1457 1.6228 +i*-0.4857 139429.615+i*12825.5
7339.3467 1.6219 +i*-0.4835 139488.255+i*12798.3
7389.5477 1.6211 +i*-0.4812 139546.293+i*12771.4
7439.7487 1.6203 +i*-0.4791 139603.741+i*12744.7
7489.9497 1.6195 +i*-0.4769 139660.610+i*12718.2
7540.1508 1.6186 +i*-0.4748 139716.912+i*12691.9
7590.3518 1.6178 +i*-0.4727 139772.657+i*12665.7
7640.5528 1.617 +i*-0.4706 139827.855+i*12639.8
7690.7538 1.6162 +i*-0.4686 139882.516+i*12614.1
7740.9548 1.6154 +i*-0.4665 139936.651+i*12588.6
7791.1558 1.6146 +i*-0.4645 139990.268+i*12563.3
7841.3568 1.6138 +i*-0.4625 140043.377+i*12538.2
7891.5578 1.6131 +i*-0.4606 140095.987+i*12513.3
7941.7588 1.6123 +i*-0.4587 140148.107+i*12488.5
7991.9598 1.6115 +i*-0.4567 140199.745+i*12464.0
8042.1608 1.6108 +i*-0.4548 140250.91 +i*12439.6
8092.3618 1.61 +i*-0.4530 140301.609+i*12415.4
8142.5628 1.6092 +i*-0.4511 140351.851+i*12391.4
8192.7638 1.6085 +i*-0.4493 140401.644+i*12367.5
8242.9648 1.6078 +i*-0.4475 140450.995+i*12343.8
8293.1658 1.607 +i*-0.4457 140499.911+i*12320.3
8343.3668 1.6063 +i*-0.4439 140548.399+i*12297.0
8393.5678 1.6055 +i*-0.4421 140596.467+i*12273.8
8443.7688 1.6048 +i*-0.4404 140644.121+i*12250.8
8493.9698 1.6041 +i*-0.4387 140691.368+i*12228.0
8544.1709 1.6034 +i*-0.4370 140738.214+i*12205.3
8594.3719 1.6027 +i*-0.4353 140784.666+i*12182.8
8644.5729 1.602 +i*-0.4336 140830.729+i*12160.4
8694.7739 1.6013 +i*-0.4320 140876.410+i*12138.2
8744.9749 1.6006 +i*-0.4303 140921.715+i*12116.1
8795.1759 1.5999 +i*-0.4287 140966.649+i*12094.2
8845.3769 1.5992 +i*-0.4271 141011.218+i*12072.4
8895.5779 1.5985 +i*-0.4255 141055.428+i*12050.8
8945.7789 1.5978 +i*-0.4240 141099.283+i*12029.3
8995.9799 1.5971 +i*-0.4224 141142.789+i*12008.0
9046.1809 1.5965 +i*-0.4209 141185.952+i*11986.8
9096.3819 1.5958 +i*-0.4194 141228.775+i*11965.7
9146.5829 1.5951 +i*-0.4178 141271.264+i*11944.8
9196.7839 1.5945 +i*-0.4163 141313.424+i*11924.0
9246.9849 1.5938 +i*-0.4149 141355.259+i*11903.3
9297.1859 1.5932 +i*-0.4134 141396.773+i*11882.8
9347.3869 1.5925 +i*-0.4119 141437.972+i*11862.4
9397.5879 1.5919 +i*-0.4105 141478.860+i*11842.1
9447.7889 1.5912 +i*-0.4091 141519.441+i*11822.0
9497.9899 1.5906 +i*-0.4077 141559.718+i*11802.0
9548.191 1.59 +i*-0.4063 141599.697+i*11782.1
9598.392 1.5894 +i*-0.4049 141639.380+i*11762.3
9648.593 1.5887 +i*-0.4035 141678.773+i*11742.6
9698.794 1.5881 +i*-0.4021 141717.878+i*11723.1
9748.995 1.5875 +i*-0.4008 141756.700+i*11703.7
9799.196 1.5869 +i*-0.3995 141795.243+i*11684.4
9849.397 1.5863 +i*-0.3981 141833.509+i*11665.2
9899.598 1.5857 +i*-0.3968 141871.502+i*11646.1
9949.799 1.5851 +i*-0.3955 141909.226+i*11627.1
10000 1.5845 +i*-0.3942 141946.684+i*11608.3
Complete point value output for the with-attenuation case.
Document this in readme.md.
Milestone 2: Standardize and document error reporting
One thing I didn't put into my interface document is a way for deal II to report errors. !!!!
We're going to need a way to report them in a way that gives the user a decent chance to find out what's wrong with ow he setup the mesh and interface areas.
Can you add this to the interface document? It can't just be "error while reading file"... It has to give them enough information to help figure out what's wrong.
My general rule for how software should work is : what would YOU want the software to do/say if you were the user and not the programmer??? That's the way the software should be written.
Milestone 1: Create success file
<output prefix>success_signal.txt
This file is created by dealII_3D and indicates that all of the frequencies specified in <cmd_filename> have been successfully solved for. dealII_3D will terminate after this file is created.
Milestone 1: Implement a way to scale STL geometry to SI units
there is some confusion with units because the STL file does not use
SI units but uses mm for lengths, whereas internally I use SI units)
When I send an STL file to someone, I usually tell them the units. Usually it's mm, sometimes it's inches. Then when the people import it into whatever CAD software they're using, they apply a scaling factor to represent it in whatever units they use in their CAD file.
So I suspect we'll need to have them provide the units used in the STL file, or perhaps a simple multiplicative value that will convert the units to meters.
Fix check_for_termination
In both the Helmholtz and Biharmonic programs, check_for_termination_signal()
exits the program when requested, but also returns a value. Furthermore, that value is ignored almost everywhere. Fix this in both program: if it exits, we don't need a return value.
Milestone 1: Output pressure and velocity at "listening points"
From Jason:
Also, I wanted the ability to output information on "listening" points inside of the volume. I don't see this mentioned anywhere in your SOW that I have. Along with the port interface areas, I want the user to be able to specify points in the volume where the pressure and particle velocity vector can be computed for. This will consist of a similar transfer matrix dependent on the inputs at the ports.
Validate the correctness of results for the 5cm cylinder
See an earlier email by Jason with the specifics of the testcase.
Figure out whether the solution for the test case should/should not be symmetric
Make sure we write to the common output files under a mutex.
Document the surface output in the readme.md file
Milestone 2: Write each frequency results into its own files
From Jason:
One thing that I had to do with biharmonic was to have it write each frequency result to it's own file. I tried to use frequency_response.txt, but I'd get read/write collision errors when your solver and Ares were both trying to access it at the same time.
Finish documenting the plain text output files
Add a table of contents to the readme.md file
Using this style:
# Table of contents
1. [Introduction](#introduction)
2. [Some paragraph](#paragraph1)
1. [Sub paragraph](#subparagraph1)
3. [Another paragraph](#paragraph2)
## This is the introduction <a name="introduction"></a>
Some introduction text, formatted in heading 2 style
## Some paragraph <a name="paragraph1"></a>
The first paragraph text
### Sub paragraph <a name="subparagraph1"></a>
This is a sub paragraph, formatted in heading 3 style
## Another paragraph <a name="paragraph2"></a>
The second paragraph text
Milestone 2: Use file prefix
Related to #12:
Also, I want all of the files to start with their own unique prefix so I can run multiple instances of your solver and still use the same directory. So each instance will have it's own unique prefix for the interface files.
Standardize writing complex numbers as a+bj
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.