notebook.community

Edit and run 





In [2]:



    


import sys
sys.path.append('/home/fjargsto/AcousticBEM/Python')

from HelmholtzSolver3D import *
from ExampleBoundaries import Sphere

# Test Problem 1
# Dirichlet boundary condition with phi chose to be the same
# as a point-source at the center of the square.
#
frequency = 100.0 # frequency [Hz]
k = frequencyToWavenumber(frequency)

solver = HelmholtzSolver3D(*(Sphere()))

boundaryCondition = BoundaryCondition(solver.aElement.shape[0])
boundaryCondition.alpha.fill(1.0)
boundaryCondition.beta.fill(0.0)

for i in range(solver.aCenters.shape[0]):
    R = norm(solver.aCenters[i, :])
    boundaryCondition.f[i] = np.exp(1j * k * R) / R

boundaryIncidence = BoundaryIncidence(solver.aElement.shape[0])
boundaryIncidence.phi.fill(0.0)
boundaryIncidence.v.fill(0.0)

exteriorPoints = np.array([[0.0000, 0.0000,  2.0000],
                           [0.0500, 0.0000,  4.0000],
                           [0.1000, 0.0000,  8.0000],
                           [0.0500, 0.1000, -2.0000]], dtype=np.float32)

exteriorIncidentPhi = np.zeros(exteriorPoints.shape[0], dtype=np.complex64)

boundarySolution = solver.solveExteriorBoundary(k, boundaryCondition, boundaryIncidence)
exteriorPhi = solver.solveExterior(boundarySolution, exteriorIncidentPhi, exteriorPoints)
print "\n\nTest Problem 1"
print "==============\n"
print boundarySolution
printInteriorSolution(boundarySolution, exteriorPhi)

# Test Problem 2
# Neumann boundary condition with v chosen to be the same
# as a point-source at the center of the square.
#
boundaryCondition.alpha.fill(0.0)
boundaryCondition.beta.fill(1.0)

for i in range(solver.aCenters.shape[0]):
    R = norm(solver.aCenters[i, :])
    a = solver.aVertex[solver.aElement[i, 0], :]
    b = solver.aVertex[solver.aElement[i, 1], :]
    c = solver.aVertex[solver.aElement[i, 2], :]
    normal = Normal3D(a, b, c)
    dot = np.dot(solver.aCenters[i, :], normal)
    boundaryCondition.f[i] = np.exp(1j * k * R) * dot * (1j * k * R - 1.0) / R**3

boundarySolution = solver.solveExteriorBoundary(k, boundaryCondition, boundaryIncidence)
exteriorPhi = solver.solveExterior(boundarySolution, exteriorIncidentPhi, exteriorPoints)
print "\n\nTest Problem 2"
print "==============\n"
print boundarySolution
printInteriorSolution(boundarySolution, exteriorPhi)

# Test Problem 3
# Neumann boundary condition with v=0 on the full boundary.
# The 400Hz field in this case is generated by a point source
# at (0.5, 0.25)
boundaryCondition.alpha.fill(1.0)
boundaryCondition.beta.fill(0.0)

sourcePoint = np.array([0.0, 0.0, 3.0], dtype=np.float32)
for i in range(solver.aCenters.shape[0]):
    x = solver.aCenters[i, :]
    r = sourcePoint - x
    R = norm(r)
    boundaryIncidence.phi[i] = np.exp(1j * k * R) / (4.0 * np.pi * R)
    boundaryCondition.f[i] = boundaryIncidence.phi[i]
    a = solver.aVertex[solver.aElement[i, 0], :]
    b = solver.aVertex[solver.aElement[i, 1], :]
    c = solver.aVertex[solver.aElement[i, 2], :]
    normal = Normal3D(a, b, c)
    dot = -np.dot(r, normal)
    boundaryIncidence.v[i] = np.exp(1j * k * R) * dot * (1j * k * R - 1.0) / (4.0 * np.pi * R**3)

boundarySolution = solver.solveExteriorBoundary(k, boundaryCondition, boundaryIncidence)

for i in range(exteriorPoints.shape[0]):
    R = norm(exteriorPoints[i, :] - sourcePoint)
    exteriorIncidentPhi[i] = np.exp(1j * k * R) / (4.0 * np.pi * R)
    
exteriorPhi = solver.solveExterior(boundarySolution, exteriorIncidentPhi, exteriorPoints)
print "\n\nTest Problem 3"
print "==============\n"
print boundarySolution
printInteriorSolution(boundarySolution, exteriorPhi)

# Test Problem 4
# Neumann boundary condition with v=0 on the full boundary.
# The 400Hz field in this case is generated by a point source
# at (0.5, 0.25)
boundaryCondition.alpha.fill(0.0)
boundaryCondition.beta.fill(1.0)
boundaryCondition.f.fill(0.0)

sourcePoint = np.array([0.0, 0.0, 3.0], dtype=np.float32)
for i in range(solver.aCenters.shape[0]):
    x = solver.aCenters[i, :]
    r = sourcePoint - x
    R = norm(r)
    boundaryIncidence.phi[i] = np.exp(1j * k * R) / (4.0 * np.pi * R)  
    a = solver.aVertex[solver.aElement[i, 0], :]
    b = solver.aVertex[solver.aElement[i, 1], :]
    c = solver.aVertex[solver.aElement[i, 2], :]
    normal = Normal3D(a, b, c)
    dot = -np.dot(r, normal)
    boundaryIncidence.v[i] = np.exp(1j * k * R) * dot * (1j * k * R - 1.0) / (4.0 * np.pi * R**3)

boundarySolution = solver.solveExteriorBoundary(k, boundaryCondition, boundaryIncidence)

for i in range(exteriorPoints.shape[0]):
    R = norm(exteriorPoints[i, :] - sourcePoint)
    exteriorIncidentPhi[i] = np.exp(1j * k * R) / (4.0 * np.pi * R)
    
exteriorPhi = solver.solveExterior(boundarySolution, exteriorIncidentPhi, exteriorPoints)
print "\n\nTest Problem 4"
print "==============\n"
print boundarySolution
printInteriorSolution(boundarySolution, exteriorPhi)



    


















    







Test Problem 1
==============

Density of medium:      1.205 kg/m^3
Speed of sound:         344.0 m/s
Wavenumber (Frequency): 1.82650735674 (100.0 Hz)

index          Potential                   Pressure                    Velocity              Intensity
    1  -5.9354e-02+  1.1234e+00i   -8.5057e+02+ -4.4938e+01i   -1.8609e+00+ -1.1710e+00i     8.1774e+02
    2  -5.9402e-02+  1.1234e+00i   -8.5054e+02+ -4.4974e+01i   -1.8610e+00+ -1.1710e+00i     8.1777e+02
    3  -5.9354e-02+  1.1234e+00i   -8.5057e+02+ -4.4938e+01i   -1.8609e+00+ -1.1711e+00i     8.1774e+02
    4  -5.9354e-02+  1.1234e+00i   -8.5057e+02+ -4.4938e+01i   -1.8609e+00+ -1.1710e+00i     8.1774e+02
    5  -5.9402e-02+  1.1234e+00i   -8.5054e+02+ -4.4974e+01i   -1.8610e+00+ -1.1710e+00i     8.1777e+02
    6  -5.9354e-02+  1.1234e+00i   -8.5057e+02+ -4.4938e+01i   -1.8609e+00+ -1.1711e+00i     8.1774e+02
    7  -4.8280e-02+  1.1310e+00i   -8.5632e+02+ -3.6554e+01i   -1.8703e+00+ -1.1178e+00i     8.2122e+02
    8   1.0774e-02+  1.1696e+00i   -8.8553e+02+  8.1574e+00i   -1.9635e+00+ -1.1314e+00i     8.6474e+02
    9  -4.8379e-02+  1.1309e+00i   -8.5627e+02+ -3.6629e+01i   -1.8702e+00+ -1.1176e+00i     8.2118e+02
   10   1.0774e-02+  1.1696e+00i   -8.8553e+02+  8.1574e+00i   -1.9635e+00+ -1.1314e+00i     8.6474e+02
   11  -4.8280e-02+  1.1310e+00i   -8.5632e+02+ -3.6554e+01i   -1.8703e+00+ -1.1178e+00i     8.2122e+02
   12   1.0515e-02+  1.1694e+00i   -8.8541e+02+  7.9609e+00i   -1.9629e+00+ -1.1313e+00i     8.6448e+02
   13  -4.8280e-02+  1.1310e+00i   -8.5632e+02+ -3.6554e+01i   -1.8703e+00+ -1.1178e+00i     8.2122e+02
   14   1.0774e-02+  1.1696e+00i   -8.8553e+02+  8.1574e+00i   -1.9635e+00+ -1.1314e+00i     8.6474e+02
   15  -4.8379e-02+  1.1309e+00i   -8.5627e+02+ -3.6629e+01i   -1.8702e+00+ -1.1176e+00i     8.2118e+02
   16   1.0774e-02+  1.1696e+00i   -8.8553e+02+  8.1574e+00i   -1.9635e+00+ -1.1314e+00i     8.6474e+02
   17  -4.8280e-02+  1.1310e+00i   -8.5632e+02+ -3.6554e+01i   -1.8703e+00+ -1.1178e+00i     8.2122e+02
   18   1.0515e-02+  1.1694e+00i   -8.8541e+02+  7.9609e+00i   -1.9629e+00+ -1.1313e+00i     8.6448e+02
   19  -4.8280e-02+  1.1310e+00i   -8.5632e+02+ -3.6554e+01i   -1.8703e+00+ -1.1178e+00i     8.2122e+02
   20   1.0774e-02+  1.1696e+00i   -8.8553e+02+  8.1574e+00i   -1.9635e+00+ -1.1314e+00i     8.6474e+02
   21  -4.8379e-02+  1.1309e+00i   -8.5627e+02+ -3.6629e+01i   -1.8702e+00+ -1.1176e+00i     8.2118e+02
   22   1.0774e-02+  1.1696e+00i   -8.8553e+02+  8.1574e+00i   -1.9635e+00+ -1.1314e+00i     8.6474e+02
   23  -4.8280e-02+  1.1310e+00i   -8.5632e+02+ -3.6554e+01i   -1.8703e+00+ -1.1178e+00i     8.2122e+02
   24   1.0515e-02+  1.1694e+00i   -8.8541e+02+  7.9609e+00i   -1.9629e+00+ -1.1313e+00i     8.6448e+02
   25  -4.8280e-02+  1.1310e+00i   -8.5632e+02+ -3.6554e+01i   -1.8703e+00+ -1.1178e+00i     8.2122e+02
   26   1.0774e-02+  1.1696e+00i   -8.8553e+02+  8.1574e+00i   -1.9635e+00+ -1.1314e+00i     8.6474e+02
   27  -4.8379e-02+  1.1309e+00i   -8.5627e+02+ -3.6629e+01i   -1.8702e+00+ -1.1176e+00i     8.2118e+02
   28   1.0774e-02+  1.1696e+00i   -8.8553e+02+  8.1574e+00i   -1.9635e+00+ -1.1314e+00i     8.6474e+02
   29  -4.8280e-02+  1.1310e+00i   -8.5632e+02+ -3.6554e+01i   -1.8703e+00+ -1.1178e+00i     8.2122e+02
   30   1.0515e-02+  1.1694e+00i   -8.8541e+02+  7.9609e+00i   -1.9629e+00+ -1.1313e+00i     8.6448e+02
   31  -5.9354e-02+  1.1234e+00i   -8.5057e+02+ -4.4938e+01i   -1.8609e+00+ -1.1710e+00i     8.1774e+02
   32  -5.9402e-02+  1.1234e+00i   -8.5054e+02+ -4.4974e+01i   -1.8610e+00+ -1.1710e+00i     8.1777e+02
   33  -5.9354e-02+  1.1234e+00i   -8.5057e+02+ -4.4938e+01i   -1.8609e+00+ -1.1711e+00i     8.1774e+02
   34  -5.9354e-02+  1.1234e+00i   -8.5057e+02+ -4.4938e+01i   -1.8609e+00+ -1.1710e+00i     8.1774e+02
   35  -5.9402e-02+  1.1234e+00i   -8.5054e+02+ -4.4974e+01i   -1.8610e+00+ -1.1710e+00i     8.1777e+02
   36  -5.9354e-02+  1.1234e+00i   -8.5057e+02+ -4.4938e+01i   -1.8609e+00+ -1.1711e+00i     8.1774e+02


Sound pressure at the sample points

index          Potential                    Pressure               Magnitude         Phase

    1  -4.6280e-01+ -1.8973e-01i    1.4365e+02+ -3.5040e+02i     1.4555e+02 dB       -1.1817
    2   1.5555e-01+  1.9605e-01i   -1.4843e+02+  1.1777e+02i     1.3953e+02 dB       2.4709
    3  -4.3217e-02+  1.1747e-01i   -8.8943e+01+ -3.2721e+01i     1.3351e+02 dB       -2.7891
    4  -4.6099e-01+ -1.9206e-01i    1.4541e+02+ -3.4903e+02i     1.4553e+02 dB       -1.1760


Test Problem 2
==============

Density of medium:      1.205 kg/m^3
Speed of sound:         344.0 m/s
Wavenumber (Frequency): 1.82650735674 (100.0 Hz)

index          Potential                   Pressure                    Velocity              Intensity
    1  -1.3097e-01+  1.2302e+00i   -9.3138e+02+ -9.9161e+01i   -1.9844e+00+ -1.3717e+00i     9.9211e+02
    2  -1.3108e-01+  1.2301e+00i   -9.3136e+02+ -9.9247e+01i   -1.9843e+00+ -1.3717e+00i     9.9209e+02
    3  -1.3097e-01+  1.2302e+00i   -9.3138e+02+ -9.9161e+01i   -1.9844e+00+ -1.3717e+00i     9.9211e+02
    4  -1.3097e-01+  1.2302e+00i   -9.3138e+02+ -9.9161e+01i   -1.9844e+00+ -1.3717e+00i     9.9211e+02
    5  -1.3108e-01+  1.2301e+00i   -9.3136e+02+ -9.9247e+01i   -1.9843e+00+ -1.3717e+00i     9.9209e+02
    6  -1.3097e-01+  1.2302e+00i   -9.3138e+02+ -9.9161e+01i   -1.9844e+00+ -1.3717e+00i     9.9211e+02
    7  -1.1266e-01+  1.2612e+00i   -9.5489e+02+ -8.5299e+01i   -2.0097e+00+ -1.3675e+00i     1.0178e+03
    8  -3.9075e-02+  1.2997e+00i   -9.8403e+02+ -2.9584e+01i   -2.1415e+00+ -1.3437e+00i     1.0735e+03
    9  -1.1283e-01+  1.2612e+00i   -9.5488e+02+ -8.5428e+01i   -2.0094e+00+ -1.3676e+00i     1.0178e+03
   10  -3.9074e-02+  1.2997e+00i   -9.8403e+02+ -2.9584e+01i   -2.1415e+00+ -1.3437e+00i     1.0735e+03
   11  -1.1266e-01+  1.2612e+00i   -9.5489e+02+ -8.5299e+01i   -2.0097e+00+ -1.3675e+00i     1.0178e+03
   12  -3.9372e-02+  1.2995e+00i   -9.8392e+02+ -2.9809e+01i   -2.1408e+00+ -1.3437e+00i     1.0732e+03
   13  -1.1266e-01+  1.2612e+00i   -9.5489e+02+ -8.5299e+01i   -2.0097e+00+ -1.3675e+00i     1.0178e+03
   14  -3.9075e-02+  1.2997e+00i   -9.8403e+02+ -2.9584e+01i   -2.1415e+00+ -1.3437e+00i     1.0735e+03
   15  -1.1283e-01+  1.2612e+00i   -9.5488e+02+ -8.5428e+01i   -2.0094e+00+ -1.3676e+00i     1.0178e+03
   16  -3.9074e-02+  1.2997e+00i   -9.8403e+02+ -2.9584e+01i   -2.1415e+00+ -1.3437e+00i     1.0735e+03
   17  -1.1266e-01+  1.2612e+00i   -9.5489e+02+ -8.5299e+01i   -2.0097e+00+ -1.3675e+00i     1.0178e+03
   18  -3.9372e-02+  1.2995e+00i   -9.8392e+02+ -2.9809e+01i   -2.1408e+00+ -1.3437e+00i     1.0732e+03
   19  -1.1266e-01+  1.2612e+00i   -9.5489e+02+ -8.5299e+01i   -2.0097e+00+ -1.3675e+00i     1.0178e+03
   20  -3.9075e-02+  1.2997e+00i   -9.8403e+02+ -2.9584e+01i   -2.1415e+00+ -1.3437e+00i     1.0735e+03
   21  -1.1283e-01+  1.2612e+00i   -9.5488e+02+ -8.5428e+01i   -2.0094e+00+ -1.3676e+00i     1.0178e+03
   22  -3.9075e-02+  1.2997e+00i   -9.8403e+02+ -2.9584e+01i   -2.1415e+00+ -1.3437e+00i     1.0735e+03
   23  -1.1266e-01+  1.2612e+00i   -9.5489e+02+ -8.5299e+01i   -2.0097e+00+ -1.3675e+00i     1.0178e+03
   24  -3.9372e-02+  1.2995e+00i   -9.8392e+02+ -2.9809e+01i   -2.1408e+00+ -1.3437e+00i     1.0732e+03
   25  -1.1266e-01+  1.2612e+00i   -9.5489e+02+ -8.5299e+01i   -2.0097e+00+ -1.3675e+00i     1.0178e+03
   26  -3.9075e-02+  1.2997e+00i   -9.8403e+02+ -2.9584e+01i   -2.1415e+00+ -1.3437e+00i     1.0735e+03
   27  -1.1283e-01+  1.2612e+00i   -9.5488e+02+ -8.5428e+01i   -2.0094e+00+ -1.3676e+00i     1.0178e+03
   28  -3.9075e-02+  1.2997e+00i   -9.8403e+02+ -2.9584e+01i   -2.1415e+00+ -1.3437e+00i     1.0735e+03
   29  -1.1266e-01+  1.2612e+00i   -9.5489e+02+ -8.5299e+01i   -2.0097e+00+ -1.3675e+00i     1.0178e+03
   30  -3.9372e-02+  1.2995e+00i   -9.8392e+02+ -2.9809e+01i   -2.1408e+00+ -1.3437e+00i     1.0732e+03
   31  -1.3097e-01+  1.2302e+00i   -9.3138e+02+ -9.9161e+01i   -1.9844e+00+ -1.3717e+00i     9.9211e+02
   32  -1.3108e-01+  1.2301e+00i   -9.3136e+02+ -9.9247e+01i   -1.9843e+00+ -1.3717e+00i     9.9209e+02
   33  -1.3097e-01+  1.2302e+00i   -9.3138e+02+ -9.9161e+01i   -1.9844e+00+ -1.3717e+00i     9.9211e+02
   34  -1.3097e-01+  1.2302e+00i   -9.3138e+02+ -9.9161e+01i   -1.9844e+00+ -1.3717e+00i     9.9211e+02
   35  -1.3108e-01+  1.2301e+00i   -9.3136e+02+ -9.9247e+01i   -1.9843e+00+ -1.3717e+00i     9.9209e+02
   36  -1.3097e-01+  1.2302e+00i   -9.3138e+02+ -9.9161e+01i   -1.9844e+00+ -1.3717e+00i     9.9211e+02


Sound pressure at the sample points

index          Potential                    Pressure               Magnitude         Phase

    1  -5.0106e-01+ -2.3891e-01i    1.8088e+02+ -3.7937e+02i     1.4645e+02 dB       -1.1259
    2   1.6130e-01+  2.2745e-01i   -1.7221e+02+  1.2212e+02i     1.4047e+02 dB       2.5247
    3  -5.5039e-02+  1.2836e-01i   -9.7183e+01+ -4.1672e+01i     1.3446e+02 dB       -2.7365
    4  -4.9894e-01+ -2.4135e-01i    1.8273e+02+ -3.7776e+02i     1.4644e+02 dB       -1.1203


Test Problem 3
==============

Density of medium:      1.205 kg/m^3
Speed of sound:         344.0 m/s
Wavenumber (Frequency): 1.82650735674 (100.0 Hz)

index          Potential                   Pressure                    Velocity              Intensity
    1  -1.9228e-02+ -2.9438e-02i    2.2288e+01+ -1.4558e+01i   -4.9138e-02+  1.9805e-02i    -6.9175e-01
    2  -1.9228e-02+ -2.9438e-02i    2.2288e+01+ -1.4558e+01i   -4.9147e-02+  1.9801e-02i    -6.9183e-01
    3  -1.9228e-02+ -2.9438e-02i    2.2288e+01+ -1.4558e+01i   -4.9138e-02+  1.9805e-02i    -6.9175e-01
    4  -1.9228e-02+ -2.9438e-02i    2.2288e+01+ -1.4558e+01i   -4.9138e-02+  1.9805e-02i    -6.9175e-01
    5  -1.9228e-02+ -2.9438e-02i    2.2288e+01+ -1.4558e+01i   -4.9147e-02+  1.9801e-02i    -6.9183e-01
    6  -1.9228e-02+ -2.9438e-02i    2.2288e+01+ -1.4558e+01i   -4.9138e-02+  1.9805e-02i    -6.9175e-01
    7   4.7159e-03+ -2.9464e-02i    2.2308e+01+  3.5705e+00i   -8.1902e-03+ -5.5131e-03i    -1.0120e-01
    8   1.5058e-02+ -2.2965e-02i    1.7388e+01+  1.1401e+01i    6.7827e-03+  2.8962e-03i     7.5476e-02
    9   4.7168e-03+ -2.9464e-02i    2.2308e+01+  3.5712e+00i   -8.1721e-03+ -5.4868e-03i    -1.0095e-01
   10   1.5058e-02+ -2.2965e-02i    1.7388e+01+  1.1401e+01i    6.7827e-03+  2.8962e-03i     7.5476e-02
   11   4.7159e-03+ -2.9464e-02i    2.2308e+01+  3.5705e+00i   -8.1902e-03+ -5.5131e-03i    -1.0120e-01
   12   1.5060e-02+ -2.2964e-02i    1.7387e+01+  1.1402e+01i    6.8121e-03+  2.9187e-03i     7.5858e-02
   13   4.7159e-03+ -2.9464e-02i    2.2308e+01+  3.5705e+00i   -8.1902e-03+ -5.5131e-03i    -1.0120e-01
   14   1.5058e-02+ -2.2965e-02i    1.7388e+01+  1.1401e+01i    6.7827e-03+  2.8962e-03i     7.5476e-02
   15   4.7168e-03+ -2.9464e-02i    2.2308e+01+  3.5712e+00i   -8.1721e-03+ -5.4868e-03i    -1.0095e-01
   16   1.5058e-02+ -2.2965e-02i    1.7388e+01+  1.1401e+01i    6.7827e-03+  2.8962e-03i     7.5476e-02
   17   4.7159e-03+ -2.9464e-02i    2.2308e+01+  3.5705e+00i   -8.1902e-03+ -5.5131e-03i    -1.0120e-01
   18   1.5060e-02+ -2.2964e-02i    1.7387e+01+  1.1402e+01i    6.8121e-03+  2.9187e-03i     7.5858e-02
   19   2.2262e-02+  3.6185e-03i   -2.7397e+00+  1.6855e+01i   -4.2742e-03+  2.6234e-02i     2.2694e-01
   20   2.3410e-02+ -4.9282e-03i    3.7312e+00+  1.7724e+01i    3.7485e-03+  1.5965e-02i     1.4847e-01
   21   2.2262e-02+  3.6190e-03i   -2.7400e+00+  1.6855e+01i   -4.2573e-03+  2.6224e-02i     2.2684e-01
   22   2.3410e-02+ -4.9282e-03i    3.7312e+00+  1.7724e+01i    3.7485e-03+  1.5965e-02i     1.4847e-01
   23   2.2262e-02+  3.6185e-03i   -2.7397e+00+  1.6855e+01i   -4.2742e-03+  2.6234e-02i     2.2694e-01
   24   2.3410e-02+ -4.9268e-03i    3.7302e+00+  1.7724e+01i    3.7479e-03+  1.5939e-02i     1.4825e-01
   25   2.2262e-02+  3.6185e-03i   -2.7397e+00+  1.6855e+01i   -4.2742e-03+  2.6234e-02i     2.2694e-01
   26   2.3410e-02+ -4.9282e-03i    3.7312e+00+  1.7724e+01i    3.7485e-03+  1.5965e-02i     1.4847e-01
   27   2.2262e-02+  3.6190e-03i   -2.7400e+00+  1.6855e+01i   -4.2573e-03+  2.6224e-02i     2.2684e-01
   28   2.3410e-02+ -4.9282e-03i    3.7312e+00+  1.7724e+01i    3.7485e-03+  1.5965e-02i     1.4847e-01
   29   2.2262e-02+  3.6185e-03i   -2.7397e+00+  1.6855e+01i   -4.2742e-03+  2.6234e-02i     2.2694e-01
   30   2.3410e-02+ -4.9268e-03i    3.7302e+00+  1.7724e+01i    3.7479e-03+  1.5939e-02i     1.4825e-01
   31   1.6509e-02+  1.2863e-02i   -9.7388e+00+  1.2499e+01i   -2.4259e-02+  2.6930e-02i     2.8643e-01
   32   1.6509e-02+  1.2863e-02i   -9.7389e+00+  1.2499e+01i   -2.4259e-02+  2.6934e-02i     2.8646e-01
   33   1.6509e-02+  1.2863e-02i   -9.7388e+00+  1.2499e+01i   -2.4259e-02+  2.6930e-02i     2.8643e-01
   34   1.6509e-02+  1.2863e-02i   -9.7388e+00+  1.2499e+01i   -2.4259e-02+  2.6930e-02i     2.8643e-01
   35   1.6509e-02+  1.2863e-02i   -9.7389e+00+  1.2499e+01i   -2.4259e-02+  2.6934e-02i     2.8646e-01
   36   1.6509e-02+  1.2863e-02i   -9.7388e+00+  1.2499e+01i   -2.4259e-02+  2.6930e-02i     2.8643e-01


Sound pressure at the sample points

index          Potential                    Pressure               Magnitude         Phase

    1  -1.9481e-02+  7.7344e-02i   -5.8559e+01+ -1.4749e+01i     1.2960e+02 dB       -2.8949
    2  -2.0465e-02+  7.6555e-02i   -5.7962e+01+ -1.5494e+01i     1.2954e+02 dB       -2.8804
    3  -1.5173e-02+  4.3979e-03i   -3.3297e+00+ -1.1488e+01i     1.1553e+02 dB       -1.8529
    4  -1.4195e-02+  4.5823e-03i   -3.4694e+00+ -1.0748e+01i     1.1504e+02 dB       -1.8830


Test Problem 4
==============

Density of medium:      1.205 kg/m^3
Speed of sound:         344.0 m/s
Wavenumber (Frequency): 1.82650735674 (100.0 Hz)

index          Potential                   Pressure                    Velocity              Intensity
    1  -4.3048e-02+ -3.6742e-02i    2.7818e+01+ -3.2593e+01i    0.0000e+00+  0.0000e+00i     0.0000e+00
    2  -4.3049e-02+ -3.6743e-02i    2.7819e+01+ -3.2593e+01i    0.0000e+00+  0.0000e+00i     0.0000e+00
    3  -4.3048e-02+ -3.6742e-02i    2.7818e+01+ -3.2593e+01i    0.0000e+00+  0.0000e+00i     0.0000e+00
    4  -4.3048e-02+ -3.6742e-02i    2.7818e+01+ -3.2593e+01i    0.0000e+00+  0.0000e+00i     0.0000e+00
    5  -4.3049e-02+ -3.6743e-02i    2.7819e+01+ -3.2593e+01i    0.0000e+00+  0.0000e+00i     0.0000e+00
    6  -4.3048e-02+ -3.6742e-02i    2.7818e+01+ -3.2593e+01i    0.0000e+00+  0.0000e+00i     0.0000e+00
    7   1.0180e-03+ -3.8552e-02i    2.9189e+01+  7.7076e-01i    0.0000e+00+  0.0000e+00i     0.0000e+00
    8   1.6749e-02+ -2.7232e-02i    2.0618e+01+  1.2681e+01i    0.0000e+00+  0.0000e+00i     0.0000e+00
    9   1.0245e-03+ -3.8540e-02i    2.9180e+01+  7.7566e-01i    0.0000e+00+  0.0000e+00i     0.0000e+00
   10   1.6749e-02+ -2.7232e-02i    2.0618e+01+  1.2681e+01i    0.0000e+00+  0.0000e+00i     0.0000e+00
   11   1.0180e-03+ -3.8552e-02i    2.9189e+01+  7.7076e-01i    0.0000e+00+  0.0000e+00i     0.0000e+00
   12   1.6762e-02+ -2.7223e-02i    2.0611e+01+  1.2691e+01i    0.0000e+00+  0.0000e+00i     0.0000e+00
   13   1.0180e-03+ -3.8552e-02i    2.9189e+01+  7.7076e-01i    0.0000e+00+  0.0000e+00i     0.0000e+00
   14   1.6749e-02+ -2.7232e-02i    2.0618e+01+  1.2681e+01i    0.0000e+00+  0.0000e+00i     0.0000e+00
   15   1.0245e-03+ -3.8540e-02i    2.9180e+01+  7.7566e-01i    0.0000e+00+  0.0000e+00i     0.0000e+00
   16   1.6749e-02+ -2.7232e-02i    2.0618e+01+  1.2681e+01i    0.0000e+00+  0.0000e+00i     0.0000e+00
   17   1.0180e-03+ -3.8552e-02i    2.9189e+01+  7.7076e-01i    0.0000e+00+  0.0000e+00i     0.0000e+00
   18   1.6762e-02+ -2.7223e-02i    2.0611e+01+  1.2691e+01i    0.0000e+00+  0.0000e+00i     0.0000e+00
   19   1.4134e-02+  9.3995e-03i   -7.1166e+00+  1.0701e+01i    0.0000e+00+  0.0000e+00i     0.0000e+00
   20   2.1619e-02+ -2.1880e-03i    1.6566e+00+  1.6368e+01i    0.0000e+00+  0.0000e+00i     0.0000e+00
   21   1.4143e-02+  9.3982e-03i   -7.1156e+00+  1.0708e+01i    0.0000e+00+  0.0000e+00i     0.0000e+00
   22   2.1619e-02+ -2.1880e-03i    1.6566e+00+  1.6368e+01i    0.0000e+00+  0.0000e+00i     0.0000e+00
   23   1.4134e-02+  9.3995e-03i   -7.1166e+00+  1.0701e+01i    0.0000e+00+  0.0000e+00i     0.0000e+00
   24   2.1624e-02+ -2.1961e-03i    1.6627e+00+  1.6372e+01i    0.0000e+00+  0.0000e+00i     0.0000e+00
   25   1.4134e-02+  9.3995e-03i   -7.1166e+00+  1.0701e+01i    0.0000e+00+  0.0000e+00i     0.0000e+00
   26   2.1619e-02+ -2.1880e-03i    1.6566e+00+  1.6368e+01i    0.0000e+00+  0.0000e+00i     0.0000e+00
   27   1.4143e-02+  9.3982e-03i   -7.1156e+00+  1.0708e+01i    0.0000e+00+  0.0000e+00i     0.0000e+00
   28   2.1619e-02+ -2.1880e-03i    1.6566e+00+  1.6368e+01i    0.0000e+00+  0.0000e+00i     0.0000e+00
   29   1.4134e-02+  9.3995e-03i   -7.1166e+00+  1.0701e+01i    0.0000e+00+  0.0000e+00i     0.0000e+00
   30   2.1624e-02+ -2.1961e-03i    1.6627e+00+  1.6372e+01i    0.0000e+00+  0.0000e+00i     0.0000e+00
   31  -6.1759e-03+  1.7138e-02i   -1.2975e+01+ -4.6759e+00i    0.0000e+00+  0.0000e+00i    -0.0000e+00
   32  -6.1747e-03+  1.7138e-02i   -1.2975e+01+ -4.6750e+00i    0.0000e+00+  0.0000e+00i    -0.0000e+00
   33  -6.1759e-03+  1.7138e-02i   -1.2975e+01+ -4.6759e+00i    0.0000e+00+  0.0000e+00i    -0.0000e+00
   34  -6.1759e-03+  1.7138e-02i   -1.2975e+01+ -4.6759e+00i    0.0000e+00+  0.0000e+00i    -0.0000e+00
   35  -6.1747e-03+  1.7138e-02i   -1.2975e+01+ -4.6750e+00i    0.0000e+00+  0.0000e+00i    -0.0000e+00
   36  -6.1759e-03+  1.7138e-02i   -1.2975e+01+ -4.6759e+00i    0.0000e+00+  0.0000e+00i    -0.0000e+00


Sound pressure at the sample points

index          Potential                    Pressure               Magnitude         Phase

    1  -1.7040e-02+  6.9910e-02i   -5.2930e+01+ -1.2902e+01i     1.2870e+02 dB       -2.9025
    2  -2.2546e-02+  7.8827e-02i   -5.9682e+01+ -1.7070e+01i     1.2984e+02 dB       -2.8630
    3  -1.6519e-02+  4.1656e-03i   -3.1539e+00+ -1.2507e+01i     1.1619e+02 dB       -1.8178
    4  -1.7032e-02+ -4.8191e-03i    3.6487e+00+ -1.2896e+01i     1.1652e+02 dB       -1.2951

Copyright (C) 2017 Frank Jargstorff

This file is part of the AcousticBEM library.

AcousticBEM is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.

AcousticBEM is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.

You should have received a copy of the GNU General Public License along with AcousticBEM. If not, see http://www.gnu.org/licenses/.

Content source: Gjacquenot/AcousticBEM

Similar notebooks:

notebook.community | gallery | about