AcousticPDEComponent

AcousticPDEComponent[vars,pars]

yields an acoustic PDE term component with variables vars and parameters pars.

Details

  • AcousticPDEComponent returns a sum of differential operators to be used as a part of partial differential equations:
  • AcousticPDEComponent models the propagation of sound in isotropic media in both the time and frequency domain by mechanisms such as diffusion.
  • AcousticPDEComponent models acoustic phenomena in fluids with dependent variable pressure [TemplateBox[{InterpretationBox[, 1], "Pa", pascals, "Pascals"}, QuantityTF]], independent variables [TemplateBox[{InterpretationBox[, 1], "m", meters, "Meters"}, QuantityTF]] and time variable [TemplateBox[{InterpretationBox[, 1], "s", seconds, "Seconds"}, QuantityTF]] or frequency variable [TemplateBox[{InterpretationBox[, 1], {"rad", , "/", , "s"}, radians per second, {{(, "Radians", )}, /, {(, "Seconds", )}}}, QuantityTF]].
  • Time-dependent variables vars are vars={p[t,x1,,xn],t,{x1,,xn}}.
  • Frequency-dependent variables vars are vars={p[x1,,xn],ω,{x1,,xn}}.
  • The time domain acoustics model AcousticPDEComponent is based on a wave equation with time variable , density , sound speed and sound sources and :
  • The frequency domain acoustics model AcousticPDEComponent is based on a Helmholtz equation with angular frequency :
  • The units of the acoustic PDE terms are in units of [1/TemplateBox[{InterpretationBox[, 1], {{"s", ^, 2}}, seconds squared, {"Seconds", ^, 2}}, QuantityTF]].
  • The following parameters pars can be given:
  • parameterdefaultsymbol
    "DipoleSource"{0,}, dipole source [TemplateBox[{InterpretationBox[, 1], {"N", , "/", , {"m", ^, 3}}, newtons per meter cubed, {{(, "Newtons", )}, /, {(, {"Meters", ^, 3}, )}}}, QuantityTF]]
    "MassDensity"1, density of media [TemplateBox[{InterpretationBox[, 1], {"kg", , "/", , {"m", ^, 3}}, kilograms per meter cubed, {{(, "Kilograms", )}, /, {(, {"Meters", ^, 3}, )}}}, QuantityTF]]
    "Material"Automatic
    "MonopoleSource"0, monopole source [1/TemplateBox[{InterpretationBox[, 1], {{"s", ^, 2}}, seconds squared, {"Seconds", ^, 2}}, QuantityTF]]
    "RegionSymmetry"None
    "SoundSpeed"1, speed of sound [TemplateBox[{InterpretationBox[, 1], {"m", , "/", , "s"}, meters per second, {{(, "Meters", )}, /, {(, "Seconds", )}}}, QuantityTF]]
  • All parameters may depend on any of , and as well as other dependent variables with the exception of , resulting in a nonlinear eigenvalue problem.
  • AcousticPDEComponent allows for sources in the time domain and sources in the frequency domain.
  • Monopole source,
    Dipole source,
  • A monopole source models a point source that radiates sound isotropically.
  • A dipole source models a two-point source that radiates sound anisotropically.
  • The number of independent variables in specifies the length of .
  • If no parameters are specified, the default time domain acoustics PDE is
  • If no parameters are specified, the default frequency domain acoustics PDE is
  • A possible choice for the parameter "RegionSymmetry" is "Axisymmetric".
  • "Axisymmetric" region symmetry represents a truncated cylindrical coordinate system where the cylindrical coordinates are reduced by removing the angle variable as follows:
  • dimensionreductionequation
    1D
    2D
  • If the AcousticPDEComponent depends on parameters that are specified in the association pars as ,keypi,pivi,], the parameters are replaced with .

Examples

open allclose all

Basic Examples  (4)

Define a time domain acoustic PDE term:

Define a frequency domain acoustic model:

Define model variables vars for a transient acoustic pressure field with model parameters pars:

Define initial conditions ics of a right-going sound wave :

Set up the equation with a sound hard boundary at the right end:

Solve the PDE:

Visualize the sound field in the time domain:

Define model variables vars for a frequency domain acoustic pressure field with model parameters pars:

Set up the equation with a radiation boundary at the left end:

Solve the PDE:

Visualize the sound field in the frequency domain at various frequencies :

Scope  (21)

Define a time- or frequency-independent acoustic model:

Define a time- or frequency-independent acoustic axisymmetric model:

Time Domain  (7)

Define model variables vars for a transient acoustic pressure field with model parameters pars:

Set up initial conditions ics of a right-going plane wave :

Set up the equation with an acoustic absorbing boundary at the right end for a plane wave:

Solve the PDE:

Visualize the solution:

Define model variables vars for a transient acoustic pressure field with model parameters pars:

Define initial conditions ics of a right-going plane wave :

Set up the equation with an acoustic impedance boundary at the right and an impedance of :

Solve the PDE:

Visualize the solution:

Define model variables vars for a transient acoustic pressure field with model parameters pars:

Define silent initial conditions ics:

Set up the equation with an acoustic normal velocity boundary with the sound particle velocity v of at the left end:

Solve the PDE on a refined mesh:

Visualize the solution:

Define model variables vars for a transient acoustic pressure field with model parameters pars:

Define silent initial conditions ics:

Set up the equation with an acoustic pressure boundary and a pressure source of at the left end:

Solve the PDE on a refined mesh:

Visualize the solution:

Define model variables vars for a frequency domain acoustic pressure field with model parameters pars:

Define silent initial conditions ics:

Set up the equation with an acoustic radiation boundary at the left end, a pressure source of and a radiation angle of :

Solve the PDE:

Visualize the solution:

Define model variables vars for a transient acoustic pressure field with model parameters pars:

Define initial conditions ics of a right-going plane wave :

Set up the equation with an acoustic sound hard boundary at the right end:

Solve the PDE:

Visualize the solution:

Define model variables vars for a transient acoustic pressure field with model parameters pars:

Define initial conditions of a right-going plane wave :

Set up the equation with an acoustic sound soft boundary at the right end:

Solve the PDE:

Visualize the solution:

Frequency Domain  (10)

Define a frequency domain acoustic model with particular sound speed and mass density:

Define a frequency domain acoustic model for a particular material:

Define a frequency domain acoustic model for a particular material:

Define model variables vars for a frequency domain acoustic pressure field with model parameters pars:

Set up the equation with a radiation boundary at the left end and an acoustic absorbing boundary at the right end:

Solve the PDE:

Visualize the solution in the frequency domain at various frequencies :

Convert the solution to the time domain:

Define model variables vars for a frequency domain acoustic pressure field with model parameters pars:

Set up the equation with a radiation boundary at the left, an acoustic impedance boundary at the right and an impedance of :

Solve the PDE:

Visualize the solution in the frequency domain at various frequencies :

Convert the solution to the time domain:

Define model variables vars for a frequency domain acoustic pressure field with model parameters pars:

Set up the equation with an acoustic normal velocity boundary at the left, the sound particle velocity of and an acoustic absorbing boundary at the right:

Solve the PDE:

Visualize the solution in the frequency domain at various frequencies :

Convert the solution to the time domain:

Define model variables vars for a frequency domain acoustic pressure field with model parameters pars:

Set up the equation with an acoustic pressure boundary at the left, a pressure source of and an acoustic absorbing boundary at the right:

Solve the PDE:

Visualize the solution in the frequency domain at various frequencies :

Convert the solution to the time domain:

Define model variables vars for a frequency domain acoustic pressure field with model parameters pars:

Set up the equation with an acoustic radiation boundary at the left end, a pressure source of and a radiation angle of :

Solve the PDE:

Visualize the solution in the frequency domain at various frequencies :

Convert the solution to the time domain:

Define model variables vars for a frequency domain acoustic pressure field with model parameters pars:

Set up the equation with an acoustic radiation boundary at the left, a pressure source of and an acoustic sound hard boundary at the right:

Solve the PDE:

Visualize the solution in the frequency domain at various frequencies :

Convert the solution to the time domain:

Define model variables vars for a frequency domain acoustic pressure field with model parameters pars:

Set up the equation with a radiation boundary at the left end and an acoustic absorbing boundary at the right end:

Solve the PDE:

Visualize the solution in the frequency domain at various frequencies :

Convert the solution to the time domain:

Units  (2)

Set up an acoustic time domain equation for xenon:

Set up an acoustic frequency domain model for xenon by specifying material parameters:

Applications  (1)

The following acoustic model describes an open pipe, wherein a vibrating piston is placed inside one end of the pipe while the other end of the pipe opens into an infinite domain. In this case, an impedance boundary condition is placed on one end to model the infinite domain. The pipe that will be modeled is a flanged circular pipe, as shown in the following figure:

Since the geometry of the pipe and the boundary conditions are rotationally symmetric about the axis, an axisymmetric model can be used. The governing equation for describing the sound wave propagation is the axisymmetric Helmholtz equation.

Set up the variables and parameters:

The axisymmetric geometry can be approximated by a 2D rectangle, which represents a cross-section of the pipe in the plane:

Set up the rectangle region with as the radius of the tube and as the length of the tube:

In the model, there are two boundary conditions. One is a NeumannValue, which expresses the acceleration of the piston with :

The second boundary condition is an AcousticImpedanceValue with impedance . The impedance is given by the following approximation, where is the wavenumber:

Set up the equation:

Solve the PDE with with a MaxCellMeasure defined by and with a resolution of 12 to get an accurate result:

Visualize the pressure distribution in the full 3D region:

Wolfram Research (2020), AcousticPDEComponent, Wolfram Language function, https://reference.wolfram.com/language/ref/AcousticPDEComponent.html (updated 2023).

Text

Wolfram Research (2020), AcousticPDEComponent, Wolfram Language function, https://reference.wolfram.com/language/ref/AcousticPDEComponent.html (updated 2023).

CMS

Wolfram Language. 2020. "AcousticPDEComponent." Wolfram Language & System Documentation Center. Wolfram Research. Last Modified 2023. https://reference.wolfram.com/language/ref/AcousticPDEComponent.html.

APA

Wolfram Language. (2020). AcousticPDEComponent. Wolfram Language & System Documentation Center. Retrieved from https://reference.wolfram.com/language/ref/AcousticPDEComponent.html

BibTeX

@misc{reference.wolfram_2024_acousticpdecomponent, author="Wolfram Research", title="{AcousticPDEComponent}", year="2023", howpublished="\url{https://reference.wolfram.com/language/ref/AcousticPDEComponent.html}", note=[Accessed: 12-June-2024 ]}

BibLaTeX

@online{reference.wolfram_2024_acousticpdecomponent, organization={Wolfram Research}, title={AcousticPDEComponent}, year={2023}, url={https://reference.wolfram.com/language/ref/AcousticPDEComponent.html}, note=[Accessed: 12-June-2024 ]}