Vibroacoustic Simulation. Alexander Peiffer

Чтение книги онлайн.

Читать онлайн книгу Vibroacoustic Simulation - Alexander Peiffer страница 33

Vibroacoustic Simulation - Alexander Peiffer

Скачать книгу

t EndFraction equals minus StartFraction partial-differential left-parenthesis rho u Subscript x Superscript 2 Baseline right-parenthesis Over partial-differential x EndFraction minus StartFraction partial-differential upper P Over partial-differential x EndFraction plus f Subscript x"/> (2.5)

      Here, fx=Fx/(Adx) is the volume force density (force per volume). Using the chain law the partial derivatives of the first and second term lead to

      The term in brackets is the homogeneous continuity Equation (2.1), and Equation (2.6) simplifies to

       rho StartFraction partial-differential v Subscript normal x Baseline Over partial-differential t EndFraction plus rho v Subscript normal x Baseline StartFraction partial-differential v Subscript normal x Baseline Over partial-differential x EndFraction plus StartFraction partial-differential upper P Over partial-differential x EndFraction equals f Subscript x (2.7)

      As with the conservation of mass, this can be extended to three dimensions:

      This equation is the non-linear, inviscid momentum equation called the Euler equation.

      2.2.3 Equation of State

      The above equations relate pressure, velocity and density. For further reducing this set we need a third equation. The easiest way would be to introduce the . Here we start with the first law of thermodynamics in order to show the difference between isotropic (or adiabatic) equation of state and other relationships.

       StartLayout 1st Row 1st Column d q 2nd Column period period period 3rd Column specific heat q equals q left-parenthesis upper T comma rho right-parenthesis 2nd Row 1st Column d v 2nd Column period period period 3rd Column specific volume v equals upper V slash upper M 3rd Row 1st Column upper P d v 2nd Column period period period 3rd Column specific expansion work 4th Row 1st Column d r 2nd Column period period period 3rd Column specific friction losses EndLayout (2.10)

      With the specific entropy ds=dq+drT we get:

       StartLayout 1st Row 1st Column d s equals left-parenthesis StartFraction partial-differential u Over partial-differential upper T EndFraction right-parenthesis Subscript upper T Baseline StartFraction d upper T Over upper T EndFraction 2nd Column plus 3rd Column StartFraction upper P Over upper T EndFraction d upper V 2nd Row 1st Column equals StartFraction c Subscript normal v Baseline Over upper T EndFraction d upper T 2nd Column plus 3rd Column StartFraction upper P Over upper T EndFraction d upper V 3rd Row 1st Column equals StartFraction c Subscript normal v Baseline Over upper T EndFraction d upper T 2nd Column plus 3rd Column StartFraction upper P Over upper T EndFraction d left-parenthesis StartFraction 1 Over rho EndFraction right-parenthesis 4th Row 1st Column equals StartFraction c Subscript normal v Baseline Over upper T EndFraction d upper T 2nd Column plus 3rd Column StartFraction upper P Over upper T rho squared EndFraction d rho EndLayout (2.10)

      The relation dv=d(1/ρ) comes from the fact that v is a mass specific value and therefore the reciprocal of the density ρ=1/v. For an ideal gas we have

      cp and cv are the specific thermal heat capacities for constant pressure and volume, respectively. That is the ratio of temperature change ∂T per increase of heat ∂q. From the total differential

      we can derive

       StartFraction d upper T Over upper T EndFraction equals StartFraction d upper P Over upper P EndFraction minus StartFraction d rho Over rho EndFraction (2.13)

      Using all above relations the change in density dρ is:

       d rho equals StartFraction rho Over kappa upper P EndFraction d upper P minus StartFraction rho Over c Subscript normal p Baseline EndFraction d s (2.14)

      with κ=cv/cp. In most acoustic cases the process is isotropic: i.e. time scales are too short for heat exchange in a free gas; thus ds=0, and the change of pressure per density is

      In case of constant temperature (isothermal) dT=0 we get with (2.12) and the ideal gas law (2.11):

       upper 
						<noindex><p style= Скачать книгу