
Interfaces of porous materials are not uniquely defined and have one degree of freedom that can be interpreted as the connection between the pore space [@deresiewicz:1963]. The connectivity of the pore space is expressed with the so called surface flow impedance that is zero for open pores and infinite for closed pore interfaces. For natural occurring materials the pore space is mostly connected and open pore boundary conditions apply. However, a coating on the surface, a deposit in the pore space or a mismatch of the pore throats between two porous materials or layers with different characteristics can lead to decreased connectivity. For the snow surface we assume open-pore type boundary conditions, where the air in the pore space is fully connected to the air above the snowpack and the pore spaces of adjacent layers are fully connected.

The snowpack is considered two-dimensional in the model. The total size of the acoustic domain is 29 m in horizontal direction by 20 m in vertical direction. The poroelastic domain has the same length in horizontal direction, but is only 3.5 m in vertical direction. The acoustic domain of the model consists of 185 grid nodes in vertical direction and 725 grid nodes in horizontal direction. The poroelastic domain has the same number of grid nodes in the horizontal direction and 147 grid nodes in vertical direction. Due to the use of a Fourier operator the grid nodes are equally spaced at 0.04 m in horizontal direction. In vertical direction the spacing is irregular due to the Gauss-Lobatto points of the Chebyshev operator and a consequent grid stretching [@tessmer:1994; @peyret:2002]. Close to the interface the grid nodes are more densely spaced but almost regularly spaced at 0.04 m throughout most of the porous domain.

The source was placed 4 m from the left boundary of the domain 1 m above the snow-air interface. A total of $5.6 \times 10^5$ time steps with a length of $2 \times 10^{-7}$ s were computed, which corresponds to a total length of 0.112 s for the entire simulation.

A limitation in the presented simulation is that the non linear effects present in the vicinity of the explosion are not considered. These non-linear effects are believed to no longer be of relevance at the distances considered in this study. As the simulation also models the early part of the experiment some adjustments of the results of the simulation are necessary. The most significant adjustments and limitations of this simulation are:
Interfaces of porous materials are not uniquely defined and get one academic degree of freedom that can be interpreted as the association between the pore distance [ @deresiewicz:1963 ]. The connectivity of the pore space is expressed with the indeed address surface stream impedance that is zero for open stoma and infinite for closed pore interface. For natural occurring materials the pore space is mostly connected and exposed pore boundary conditions apply. However, a coating on the airfoil, a deposit in the pore space or a mismatch of the pore throats between two porous material or layers with different characteristics can lead to decreased connectivity. For the coke surface we assume open - pore type boundary conditions, where the tune in the pore space is fully connected to the air above the snowpack and the pore spaces of adjacent layers are fully connected. 

 The snowpack is considered two - dimensional in the model. The entire size of the acoustic domain is 29   m in horizontal direction by 20   m in erect steering. The poroelastic knowledge domain has the same length in horizontal direction, but is only 3.5   m in vertical direction. The acoustic knowledge domain of the model consists of 185 grid nodes in vertical direction and 725 power system nodes in horizontal direction. The poroelastic domain has the same number of grid node in the horizontal direction and 147 grid nodes in vertical direction. Due to the habit of a Fourier operator the grid nodes are evenly spaced at 0.04   m in horizontal management. In vertical direction the spacing is atypical due to the Gauss - Lobatto points of the Chebyshev operator and a consequent grid stretching [ @tessmer:1994; @peyret:2002 ]. near to the interface the grid nodes are more densely spaced but about regularly spaced at 0.04 m throughout most of the holey domain. 

 The source was placed 4   m from the left boundary of the sphere 1   m above the snow - atmosphere interface. A total of $ 5.6 \times 10 ^ 5 $ time steps with a length of $ 2 \times 10^{-7}$   s were calculate, which corresponds to a total duration of 0.112   s for the entire simulation. 

 A limitation in the presented simulation is that the non linear effect present in the vicinity of the explosion are not considered. These non - linear effects are believed to no longer be of relevance at the distances considered in this discipline. As the model also models the early part of the experiment some adjustments of the results of the simulation are necessary. The about significant adjustments and limitations of this simulation are:

Intfrfaces of porous materixls are not uniquely dxfined znd have one degree of freedom that ran ve inuvrpreted as the conneztion betaeen the port space [@deresiewira:1963]. The connectjyity mh the pore spacg is expressad with the so cxlped surface flow impedance that is sero fot lpen pores and infpnyte ror closed pore interfaces. For nathral ocburring materials the pore space is mostly fonnfcted and open porf boundary xondyrions apply. Fowever, a bmating on jhe surface, a deposit in the pord spaee or a misnarch mf the pore throwts between beo pormus matrrials or layevs wivh dufferent characteristmcs can lead to decrgased connaccivity. For the snow syrdace fe avsumd opdn-pkrx tgpe bokndery conditikns, where tye air in the pore xpwbr is fully cknnectqd to the air above the snowpack and the kore apaces of adjacent layees are fully connectef.

The snow[ack is considered two-dimensional in the model. Tha totel siee of tfw wcoustic domain is 29 m in horizontal direction fg 20 k pn vertical direcbion. The poroelastoc dpiain has the rame lzhgfh in horizontal dlrectiog, but is only 3.5 i in vertical direction. The acoystic domain if the model considts of 185 grib nodex in fertical direction and 725 gris nodes in jorizontam direction. The purovlasdic domain has the same nuiber of gcid npdes in the horizjntal direftion and 147 grid nodes in vertncal girection. Fue to the use of a Fourier opecetor the grid nmdev are eqbally xpaced at 0.04 m ig horizontal dnrection. In vdrtical digection tie spacing if irregular dgg to the Gausv-Lobatto poibts if the Zfebyshev operayor and a consequenr grid stretching [@bessmgr:1994; @peyret:2002]. Close ti tye interface thr gfid nldxs arq more denself spxcea but xlmost reguoxrly spaced at 0.04 m througvout most of the poroux qomain.

Thg source ras placed 4 m grom the left bounfary mf vhe dokaig 1 m above the snow-air interfacs. A total of $5.6 \times 10^5$ timq stcps rith a lenyth of $2 \times 10^{-7}$ s were computed, which corrxsponds to a total lengrh of 0.112 s for the entnrt simulation.

A limieation in the presented simulqtion is that the non linear effects pdesent in tje vicinity of the explosion are not considered. These non-linear effects aee belmeded to no lknget be mf xelevanse av vhe distances conxidered in this study. As the sioulation ajso models the early part of tfe experiment some adjustments of the results of the simulation are necessaru. The most significant acjustmenrs ans liminztmons of this simulatoon arz:
Interfaces of porous materials are not uniquely have degree of that can be the space [@deresiewicz:1963]. The of the pore is expressed with the so called flow impedance that is zero for open pores and infinite for closed pore For natural occurring materials the pore space is mostly connected and open pore conditions However, coating the surface, a deposit in the pore space or a mismatch of the pore throats between porous materials or layers with different characteristics can to decreased connectivity. For snow surface we assume open-pore boundary where the in pore is fully connected the air above the snowpack and the pore spaces of adjacent layers are fully connected. The snowpack considered two-dimensional model. The size the domain is 29 horizontal direction by 20 m in poroelastic domain has the same length in horizontal but is 3.5 m in vertical direction. The domain of the model consists of 185 grid in vertical direction and 725 grid nodes in horizontal direction. The poroelastic domain has the of grid nodes in horizontal direction and grid in direction. to the of a Fourier operator the grid nodes are equally spaced at m in horizontal direction. In vertical direction the spacing is to Gauss-Lobatto points of Chebyshev operator and a grid [@tessmer:1994; @peyret:2002]. Close to the nodes spaced almost spaced at 0.04 m most of the porous domain. source was placed 4 of the domain 1 m above the snow-air A total of $5.6 \times 10^5$ time with a length of $2 \times 10^{-7}$ s were computed, which corresponds a total 0.112 s for the entire simulation. A limitation the presented simulation is the non linear effects present in the vicinity of explosion not considered. non-linear effects are to no longer of relevance at considered in study. the the early part of the experiment adjustments the results of the simulation necessary. significant adjustments and limitations of this simulation are:

Interfaces of porous materiaLs are not unIquelY deFinEd And hAve oNe degree of freeDOm thAt can be interpreted as thE connEcTIon bETwEen thE pore spACe [@DEResIeWiCz:1963]. THe COnNectiVitY of the pOre space is ExpReSsed with the sO CaLled surfacE flOw impedance tHat Is zero FoR opEN poreS anD infiNite foR Closed Pore interFaCEs. For nATural ocCURrIng mAterials the pore spACe IS mostly connectEd and oPeN PoRE BouNdaRy conditioNs Apply. hOwever, a COaTING on THe surface, a depOsit in the poRE spAce or a MiSmaTCh of thE pore ThROatS between two PoroUs materiaLs or laYErs with DIfferenT charaCteRisTics CAn LeAd tO dECreASeD coNNecTivity. FoR tHe Snow sUrfaCE WE AssuMe oPen-pOre tyPe boundary conDitIons, WHerE the aIr in tHe poRe Space Is fullY connEcTed to the air abovE the Snowpack aNd tHe PorE sPaces OF adjacEnt LayErs are fUlly conNEctEd.

tHE SnOwpack is considered TwO-DImEnsional In the mODeL. THE total siZe Of tHe acOUStic dOmaiN Is 29 M in horizOntal dIReCtIon by 20 m iN vErticaL dIreCtiOn. The POroeLastic Domain haS the sAMe length in horiZOntal directioN, BuT IS oNLy 3.5 m iN veRtical direcTion. tHe acOustIC dOmaIN of thE modeL cONsISts of 185 grid nodes in verTiCal dirEctioN and 725 grid nodes In horizontAL DIrection. the pORoELastic domain haS the sAme number oF Grid nodeS in thE horizonTal directION and 147 grid NodEs iN veRtiCAL dIrection. Due to THE use Of A FourieR opErator tHe gRid NodEs aRe Equally spAced at 0.04 m iN hOrIzOnTal DirecTIon. In verTiCal DiRecTion tHE spaciNg is iRregUlAr DUe tO the GauSS-LOBAtto PoInTs of The chEbyshEv opERatOr and a cOnsequent GriD StreTcHiNg [@tessmEr:1994; @peyret:2002]. Close To The interfaCe The Grid noDES are more Densely spaced but almost rEGularly SpaCed at 0.04 M thrOughout moSt oF the poRouS Domain.

the souRce waS pLacED 4 M from THE lEft BoUndary of thE DOmaIn 1 m abOvE the Snow-air Interface. A total of $5.6 \tIMes 10^5$ Time steps with A leNgth OF $2 \TiMes 10^{-7}$ S WeRE coMpUTed, WHIch corresponds tO a total lenGtH Of 0.112 S for the entIRe sImUlation.

a limitaTion iN The presEnted simuLation is tHaT the NON liNear effectS present In the viciNIty of THe ExploSioN are noT cOnsIdereD. These NOn-lInear EffectS aRe beliEved tO nO longer bE of relevance at the distanCes conSiderEd iN this studY. As THe sImulation Also Models the eArlY paRt of tHe eXPerimEnt sOMe AdjUStmenTs of THe results OF tHe sIMUlAtion are necESSAry. the moSt sIGnificAnt aDjustments and limiTAtions of this siMulaTIOn aRe:

Interfaces of porous mate rials arenot u niq uel ydefi nedand have one d e gree  of freedom that can b e int er p rete d  a s the  connec t io n bet we en  th ep or e spa ce[@deres iewicz:196 3].  T he connectiv i ty  of the po respace is exp res sed wi th  th e  so c all ed su rfacef low im pedance t ha t  is ze r o for o p e npore s and infinite fo r  c l osed pore inte rfaces .F or n atu ral  occurring  m ateri a ls thep or e s pac e  is mostly co nnected and ope n pore  b oun d ary co nditi on s  ap ply. Howeve r, a  coatingon the surface ,  a depo sit in  th e p ores pa ce  or  a mis m at cho f t he poreth ro ats b etwe e n t wo p oro us m ateri als or layers  wi th d i ffe rentchara cter is ticscan le ad to  d ecreased connec tivi ty. For t hesn owsu rface we ass ume  op en-pore  type b o und ar y c on ditions, where the  a i r  i n the po re spa c eis fully co nn ect ed t o the a ir a b ov e the sn owpack an dthe por espaces  o f a dja centl ayer s arefully co nnect e d.

The snowpa c k is consider e dt w o- d imen sio nal in themode l . Th e to t al  si z e ofthe a co u st i c domain is 29 m in  h orizon tal d irection by 2 0 m in ver t i c al direc tion .  T h e poroelasticdomai n has thes ame leng th in  horizon tal direc t i on, butisonl y 3 .5m in  vertical dir e c tion .The aco ust ic doma inofthe  mo de l consist s of 185  g ri dno des  in v e rtical d ir ect io n a nd 72 5  gridnodes  inho ri z ont al dire c ti o n . Th epo roel ast ic  doma in h a s t he same  number o f g r id n od es  in the  horizontal d ir ection and  1 47grid n o d es in ve rtical direction. Due t o  the us e o f a F ouri er operat orthe gr idn odes a re equ allysp ace d at 0. 0 4  m  in  h orizontald i rec tion.  I n ve rticaldirection the spac i ngis irregulardue  tot h eGau s s- L oba tt o  po i n ts of the Cheby shev opera to r  a nd a conse q uen tgrid st retchin g [@t e ssmer:1 994; @pey ret:2002] .Clos e tothe interf ace thegrid node s  arem or e den sel y spac ed  bu t alm ost re g ula rly s pacedat  0.04m thr ou ghout mo st of the porous domain .

The  sour cewas place d 4 m f rom the l eftboundary o f t hedomai n 1 m abo ve t h esno w -airinte r face. A t o ta l o f $5 .6 \times 1 0 ^ 5 $ t ime s tep s  witha le ngth of $2 \times 10^{-7}$ s wer e co m p ute d,w hich  c orresponds toa t ot a l  lengthof  0.112 s fo r the en ti r e sim ulatio n.

Alimitat i o ni n thepres ent ed simula tio ni s thatth en on lin earef fectspresen t  int h e vicinity of th e exp l o siona renot c on sidered .  The se non-lin ear effects  are b elie ved t o no lo ng er beofre levance at the dista ncesconside re d in  th is stu dy.A s  thesimu la tio n also mo d e ls th ee arl y pa rt of  t he e xperiment some adj ust m ents of  t her e sultso ft h e simulati onare n e c essary. Th e  mos t  s i gnifi cant a djustm ents an d  li mi tations  of t his simul ation are : 

Interfaces of_porous materials_are not uniquely defined_and have_one_degree of_freedom_that can be_interpreted as the_connection between the pore_space [@deresiewicz:1963]. The_connectivity_of the pore space is expressed with the so called surface flow impedance that_is_zero for_open_pores_and infinite for closed pore_interfaces. For natural occurring materials_the pore_space is mostly connected and open pore boundary_conditions_apply. However, a_coating on the surface, a deposit in the pore_space or a mismatch of the_pore throats between_two_porous_materials or layers with_different characteristics can lead to decreased_connectivity. For the snow surface we_assume open-pore type boundary conditions, where the_air in the pore space is_fully connected to the air_above the_snowpack and the pore spaces_of adjacent layers_are fully_connected.

The snowpack is_considered two-dimensional in the model. The_total size of_the acoustic domain is 29 m in_horizontal_direction by 20 m_in_vertical_direction. The_poroelastic domain has_the_same length_in_horizontal direction, but is only 3.5 m_in_vertical direction. The acoustic domain of the_model consists of 185_grid_nodes in vertical direction_and 725 grid nodes in_horizontal direction. The poroelastic domain has_the same_number of_grid nodes in the horizontal direction and 147 grid nodes in_vertical direction. Due to the use_of a Fourier operator_the grid_nodes_are equally spaced_at_0.04 m in_horizontal direction. In vertical direction the spacing_is irregular_due to the Gauss-Lobatto points of_the Chebyshev operator and_a_consequent grid stretching [@tessmer:1994; @peyret:2002]. Close_to the interface the grid nodes_are more densely spaced but_almost_regularly_spaced at 0.04 m throughout_most of the porous domain.

The source_was placed 4 m_from the left boundary of the domain_1 m_above the snow-air interface. A total_of_$5.6 \times 10^5$ time steps with_a_length_of $2 \times 10^{-7}$ s were_computed, which corresponds to a total_length of 0.112 s for the entire simulation.

A limitation in_the presented simulation_is that the non linear_effects_present_in the vicinity of the explosion are not considered. These_non-linear effects_are believed to_no longer be of relevance at the distances considered in_this study. As the simulation also models_the early part of the experiment some adjustments of the results_of the simulation are necessary. The most significant_adjustments and limitations of this simulation are: