Staggered-Grid Finite-Difference Acoustic Modeling With the Time-Domain Atmospheric Acoustic ...

107 Pages·2008·3.01 MB·English

SANDIA REPORT
SAND2006-2540
UnlimitedRelease
PrintedMay2005
Staggered-Grid Finite-Difference Acoustic
Modeling With the Time-Domain
Atmospheric Acoustic Propagation Suite
(TDAAPS)
NeillP.Symons,DavidF.Aldridge,SandiaNationalLaboratories;DavidH.Marlin,SandraL.
CollierU.S.ArmyResearchLaboratory;D.KeithWilson,U.S.ArmyColdRegionsResearch
&EngineeringLab.;VladimirE.Ostashev,NOAA/EnvironmentalTechnologyLaboratory
Preparedby
SandiaNationalLaboratories
Albuquerque,NewMexico87185andLivermore,California94550
SandiaisamultiprogramlaboratoryoperatedbySandiaCorporation,
aLockheedMartinCompany,fortheUnitedStatesDepartmentofEnergy’s
NationalNuclearSecurityAdministrationunderContractDE-AC04-94-AL85000.
Approvedforpublicrelease;furtherdisseminationunlimited.
IssuedbySandiaNationalLaboratories,operatedfortheUnitedStatesDepartmentof
EnergybySandiaCorporation.
NOTICE:Thisreportwaspreparedasanaccountofworksponsoredbyanagencyof
theUnitedStatesGovernment. NeithertheUnitedStatesGovernment,noranyagency
thereof,noranyoftheiremployees,noranyoftheircontractors,subcontractors,ortheir
employees,makeanywarranty,expressorimplied,orassumeanylegalliabilityorre-
sponsibilityfor the accuracy, completeness, or usefulness of any information, appara-
tus,product,orprocessdisclosed,orrepresentthatitsusewouldnotinfringeprivately
ownedrights. Reference hereintoanyspeci(cid:2)ccommercial product,process, orservice
bytradename,trademark,manufacturer,orotherwise,doesnotnecessarilyconstitute
or imply its endorsement, recommendation, or favoring by the United States Govern-
ment,anyagencythereof,oranyoftheircontractorsorsubcontractors. Theviewsand
opinionsexpressedhereindonotnecessarilystateorre(cid:3)ectthoseoftheUnitedStates
Government,anyagencythereof,oranyoftheircontractors.
PrintedintheUnitedStatesofAmerica. Thisreporthasbeenreproduceddirectlyfrom
thebestavailablecopy.
AvailabletoDOEandDOEcontractorsfrom
U.S.DepartmentofEnergy
Of(cid:2)ceofScienti(cid:2)candTechnicalInformation
P.O.Box62
OakRidge,TN37831
Telephone: [phone]
Facsimile: [phone]
E-Mail: [email protected]
Onlineordering: http://www.doe.gov/bridge
Availabletothepublicfrom
U.S.DepartmentofCommerce
NationalTechnicalInformationService
5285PortRoyalRd
Spring(cid:2)eld,VA22161
Telephone: [phone]
Facsimile: [phone]
E-Mail: [email protected]
Onlineordering: http://www.ntis.gov/ordering.htm
DEPARTMENTOFENERGY
•(cid:13) •(cid:13)
UNITEDSTATESOF AMERICA
2
SAND2006-2540
UnlimitedRelease
PrintedMay2005
Staggered-Grid Finite-Difference Acoustic
Modeling With the Time-Domain Atmospheric
Acoustic Propagation Suite (TDAAPS)
Neill P. Symons, David F. Aldridge
Sandia National Laboratories
Geophysics Department
Albuquerque, NM 87185-0750
[email protected]
David H. Marlin, Sandra L. Collier
U.S. Army Research Laboratory
D. Keith Wilson
U.S. Army Cold Regions Research
Engineering Lab.
Vladimir E. Ostashev
NOAA/Environmental Technology Laboratory
3
Abstract
This document is intended to serve as a users guide for the time-domain atmospheric
acoustic propagation suite (TDAAPS) program developed as part of the Department of
Defense High-Performance Modernization Of(cid:2)ce (HPCMP) Common High-Performance
ComputingScalableSoftwareInitiative(CHSSI). TDAAPSperformsstaggered-grid(cid:2)nite-
differencemodelingoftheacousticvelocity-pressuresystemwiththeincorporationofspa-
tially inhomogeneous winds. Wherever practical the control structure of the codes are
writteninC++ usingan objectoriented design. Sectionsof code where a large numberof
calculationsarerequiredarewritteninCorF77inordertoenablebettercompileroptimiza-
tion of these sections. The TDAAPS program conforms to a UNIX style calling interface.
Most of the actions of the codes are controlled by adding (cid:3)ags to the invoking command
line. Thisdocumentpresentsa large numberof examplesand providesnewusers withthe
necessarybackgroundtoperformacousticmodelingwithTDAAPS.
4
Acknowledgment
FundingforthisworkprovidedbytheDepartmentofDefenseHigh-PerformanceModern-
izationOf(cid:2)ce(HPCMP)CommonHigh-PerformanceComputingScalableSoftwareInitia-
tive (CHSSI), project CEA-11. Thanks to the TDAAPS Alpha tester Micheal White from
ERDC-CERL and Beta tester Rodney Whitaker from Los Alamos National Laboratory,
who also provided many helpful comments on this document. This document was vastly
improved after thoughtful and complete reviews by Mattew M. Haney and Sandford Bal-
lard.
Forcodedistribution,contactDr. DavidMarlin([email protected]).
5
Contents
Nomenclature 11
1 Introduction 13
1.1 DocumentPurpose..................... ........................ 13
1.2 ExportControl ........................ ........................ 13
1.3 DesignChoices ....................... ........................ 14
1.3.1 C++ andC .................... ........................ 14
1.3.2 NetCDF Files .................. ........................ 14
1.3.3 MessagePassingProtocol ......... ........................ 15
1.4 ConcurrentVersionSystem(CVS) ........ ........................ 15
1.4.1 CVS Logging .................. ........................ 17
1.5 Background .......................... ........................ 19
1.5.1 Moving-MediaAcousticEquations . ........................ 19
1.5.2 Finite-Difference . ............... ........................ 19
1.5.3 ParallelImplementation .......... ........................ 20
2 Running the TDAAPSAlgorithm 23
2.1 BasicsofRunningTDAAPS.............. ........................ 23
2.1.1 Checkpoints.................... ........................ 24
2.2 DispersionandStability................. ........................ 24
6
2.3 InputFiles ........................... ........................ 25
2.3.1 Units ......................... ........................ 25
2.3.2 ModelConstructionandMatlabTM .. ........................ 25
2.3.3 TextInputFiles ................. ........................ 30
2.4 OutputFileFormats .................... ........................ 30
2.5 BoundaryConditions ................... ........................ 32
2.5.1 Zwikker-Kosten(Mass-Resistance)BoundaryCondition......... 33
2.5.2 Rock Property(IrregularSurface) BoundaryCondition .......... 34
2.6 Quasi-Wavelets ....................... ........................ 34
2.7 UsageandDe(cid:2)nitionofFlags ............ ........................ 35
2.7.1 BNF De(cid:2)nitionofTDAAPScall .... ........................ 37
2.7.2 DescriptionofTDAAPS(cid:3)ags ...... ........................ 37
3 Model Generation 41
3.1 Introduction .......................... ........................ 41
3.2 ModelBuildingwithMatlabTM ........... ........................ 42
3.3 ModelBuildingwithbuildSgfdModel ...... ........................ 48
3.3.1 UsageandDe(cid:2)nitionofFlags...... ........................ 48
3.3.2 BNF De(cid:2)nitionofbuildSgfdModelcall ...................... 50
3.3.3 DescriptionofbuildSgfdModel(cid:3)ags. ........................ 50
4 Examples 51
4.1 TransmissionLosswithVerticalWindGradient ...................... 51
4.2 ExtinctionandCoherence ............... ........................ 54
4.3 AlphaTestScalability .................. ........................ 59
4.4 LongRangeUltra-LowFrequency ........ ........................ 62
7
Bibliography 71
Appendix
A TDAAPSCallingFlagsQuick Reference 73
B Examplesfrom the BetaTest 75
B.1 HillTest ............................. ........................ 75
B.2 ExtinctionandCoherence ............... ........................ 83
B.3 Zwikker-KostenPartiallyAbsorbingBoundaryCondition .............. 90
C Large Eddy Simulation(LES) Example 95
8
List of Figures
1.1 DoubleTime-StepTemporalUpdate ....... ........................ 21
1.2 DomainDecomposition ................. ........................ 22
2.1 ComparisonofSlicesfrom3Dand1DModel ....................... 29
2.2 QWVelocityDistributions............... ........................ 36
3.1 ThreeCrossSectionsThroughtheHillModel ....................... 49
4.1 100Hzmono-frequencysourcewavelet..... ........................ 53
4.2 TraceoutputfromTransmissionLossExample....................... 55
4.3 TraceoutputfromExtinctionandCoherence Example ................. 58
4.4 3DWindModelforLowFrequencyExample........................ 67
4.5 Time-SlicefromLowFrequencyExample .. ........................ 68
9
List of Tables
1.1 CVSRevisionNumbers: src/acoustic sgfd .. ........................ 16
1.2 CVSRevisionNumbers: src/sgfd ......... ........................ 16
1.3 CVSRevisionNumbers: src/utils ......... ........................ 17
10