So, I am running a simple acoustic simulation with the Helmholtz solver. Here my sif file:
Code: Select all
Header
CHECK KEYWORDS Warn
Mesh DB "." "."
Include Path ""
Results Directory ""
! Simulation Parameters:
! Frequency [Hz]:
$ f = 3583
! Velocity [m / s]
$ v = 0.000001
End
Simulation
Max Output Level = 5
Coordinate System = Cartesian
Coordinate Mapping(3) = 1 2 3
Simulation Type = Steady state
Steady State Max Iterations = 1
Output Intervals = 1
Timestepping Method = BDF
BDF Order = 1
Solver Input File = case.sif
Post File = case.ep
Coordinate Scaling = 0.001
End
Constants
Gravity(4) = 0 -1 0 9.82
Stefan Boltzmann = 5.67e-08
Permittivity of Vacuum = 8.8542e-12
Boltzmann Constant = 1.3807e-23
Unit Charge = 1.602e-19
End
Body 1
Target Bodies(1) = 1
Name = "Body 1"
Equation = 1
Material = 1
Body Force = 1
End
Solver 1
Equation = Helmholtz Equation
Procedure = "HelmholtzSolve" "HelmholtzSolver"
Variable = -dofs 2 Pressure Wave
Exec Solver = Always
Stabilize = True
Bubbles = False
Lumped Mass Matrix = False
Optimize Bandwidth = True
Steady State Convergence Tolerance = 1.0e-5
Nonlinear System Convergence Tolerance = 1.0e-7
Nonlinear System Max Iterations = 20
Nonlinear System Newton After Iterations = 3
Nonlinear System Newton After Tolerance = 1.0e-3
Nonlinear System Relaxation Factor = 1
Linear System Solver = Iterative
Linear System Iterative Method = BiCGStab
Linear System Max Iterations = 500
Linear System Convergence Tolerance = 1.0e-10
BiCGstabl polynomial degree = 2
Linear System Preconditioning = Diagonal
Linear System ILUT Tolerance = 1.0e-3
Linear System Abort Not Converged = False
Linear System Residual Output = 1
Linear System Precondition Recompute = 1
End
Solver 2
Equation = Result Output
Procedure = "ResultOutputSolve" "ResultOutputSolver"
Save Geometry Ids = True
Output File Name = case
Output Format = Vtu
Save Linear Elements = True
Exec Solver = After Simulation
End
Equation 1
Name = "Helmholtz EQ"
Angular Frequency = Real MATC "2 * pi * f"
Active Solvers(1) = 1
End
Equation 2
Name = "Result Output EQ"
Active Solvers(1) = 2
End
Material 1
Name = "Air (room temperature)"
Viscosity = 1.983e-5
Heat expansion Coefficient = 3.43e-3
Heat Conductivity = 0.0257
Relative Permittivity = 1.00059
Sound speed = 343.0
Heat Capacity = 1005.0
Density = 1.205
End
Body Force 1
Name = "SPL"
Press_dB = Variable Pressure Wave 1, Pressure Wave 2
Real MATC "20*log(((sqrt(tx(0)^2+tx(1)^2))/sqrt(2)))"
End
Boundary Condition 1
Target Boundaries(9) = 1 2 3 4 5 6 7 8 9
Name = "Hard Wall"
Wave Flux 1 = 0
Wave Flux 2 = 0
End
Boundary Condition 2
Target Boundaries(1) = 10
Name = "Radiator"
Wave Flux 1 = 0
Wave Flux 2 = Real MATC "2 * pi * f * 1.205 * v"
End
By the way, I am running Elmer on Mac, installed with this.
Edit: Wrong sif file uploaded
All seems to be fine on the log:
Code: Select all
ResultOutputSolver: -------------------------------------
ResultOutputSolve: Saving with prefix: case
ResultOutputSolver: Dimension of mesh: 3
ResultOutputSolver: Creating list for saving - if not present
CreateListForSaving: Field Variables for Saving
CreateListForSaving: Scalar Field 1: pressure wave 1
CreateListForSaving: Scalar Field 2: pressure wave 2
ResultOutputSolver: Saving in unstructured VTK XML (.vtu) format
VtuOutputSolver: Saving results in VTK XML format with prefix: case
VtuOutputSolver: Saving number of partitions: 1
ResultOutputSolver: -------------------------------------
Any of you able to suggest how can I export the dB values defined in Body Force 1? Should I actually define it somewhere else?