3D acoustic wave modeling with a time-space domain dispersion-relation-based Finite-difference scheme
|
|
- Abel Leonard
- 6 years ago
- Views:
Transcription
1 P-8 3D acoustic wave odeling with a tie-space doain dispersion-relation-based Finite-difference schee Yang Liu * and rinal K. Sen State Key Laboratory of Petroleu Resource and Prospecting (China University of Petroleu, Beijing), China The Institute for Geophysics, John A. and Katherine G. Jackson School of Geosciences, The University of Texas at Austin, U.S.A. Suary Spatial finite-difference (FD) stencils designed in the space doain are usually eployed in wave equation odeling. In this paper, we adopt a spatial FD stencil which is devised in the tie-space doain based on a dispersion relation to iprove the accuracy of FD odeling for 3D acoustic wave equation. In addition, we eploy a new hybrid absorbing boundary condition (ABC) to eliinate edge reflection due to finite coputational doain. Dispersion analysis and odeling results deonstrate that the tie-space doain dispersion-relation-based spatial FD stencils can indeed iprove the odeling accuracy. The odeling accuracy can be iproved further by using a truncated FD ethod. The new hybrid ABC can effectively absorb boundaries reflections. The odeling schee presented in this paper is ore accurate and efficient than the conventional one and can be used routinely. Introduction The finite-difference (FD) ethod is an iportant ethod for the nuerical solution of partial differential equations and has been widely utilized in seisic odeling and igration (e.g., Kelly et. al, 976; Dablain, 986; Liu and Sen, 9a). A high-order spatial FD is a coon approach to increase odeling accuracy (e.g., Dablain, 986; Fornberg, 987; Etgen and O Brien, 7). A loworder FD algorith uses a shorter operator but needs ore grid points for discretization. A high-order FD algorith uses a longer operator but needs fewer grid points. Since a conventional explicit high-order teporal FD is usually unstable in the wave equation odeling, spatial derivatives are used to replace high-order teporal derivatives (e.g., Dablain, 986) to increase the accuracy of teporal derivatives with additional coputational cost. Generally, ost FD ethods deterine the FD stencils for spatial derivatives only in the space doain. However, the seisic wave propagation calculation is done both in space and tie doains. To iprove the accuracy of conventional FD ethods, a new tie-space doain FD ethod was proposed to derive the spatial FD coefficients in the joint tie-space doain (Finkelstein and Kastner, 7). The key idea of the ethod is that the dispersion relation is copletely satisfied at designated frequencies and thus the spatial FD coefficients are frequency dependent. This ethod was developed further for D, D and 3D acoustic wave odeling using a plane wave theory and the Taylor series expansion (Liu and Sen, 7b). These new spatial FD coefficients, dependent on the Courant nuber and space point nuber, are frequency independent though they lead to a frequency dependent nuerical solution. In this paper, we report on the nuerical results in 3D using FD that replaces the conventional 3D spatial FD coefficients with these new coefficients to iprove the odeling accuracy without increasing the calculation aount. A truncated FD ethod is also used to enhance the odeling accuracy further. A hybrid absorbing boundary condition (ABC) is also developed for 3D acoustic odeling. Conventional finite-difference schee for the 3D acoustic wave equation We start with the 3D acoustic wave equation given by p p p p + + =, () x y z v t * China University of Petroleu, State Key Laboratory of Petroleu Resource and Prospecting, Changping, Beijing, 49, China. Eail: wliuyang@vip.sina.co.
2 3D Acoustic odeling with Dispersion-relation-based FD ethod Where, p = p( x, y, z, t) is a scalar wave field, and v is velocity. The following nd-order FD is usually used for the tie derivatives, p δ p = p,, + ( p,, + p,, ) t δt τ, () where, n p, l, j = p( x + h, y + lh, z + jh, t + nτ ), (3) τ is tie step, h is grid size. Generally, the odeling accuracy is iproved by high-order FD on the space derivatives given by p δ p p,, a ( p,, p,, ) x δ x h + +, (4a) p δ p p,, a ( p,, p,, ) y δ y h + +, (4b) p δ p p,, a ( p,, p,, ) z δ z h + +. (4c) A Taylor series expansion is generally used to deterine the FD coefficients (Fornberg, 987; Dablain, 986). The FD coefficients in Equations (4a), (4b) and (4c) are as follows (Liu and Sen, 9a) a a ( ) + n = ( =,,..., ), (5a) n < n, n a =. (5b) It can be proved that when a -order FD space stencil and a nd-order FD tie stencil are used to solve the 3D acoustic wave equation, nuerical odeling has nd-order accuracy. It is noteworthy that increasing ay help reduce the agnitude of the conventional FD error without increasing the accuracy order. Tie-space doain dispersion-relation-based finitedifference schee for the 3D acoustic wave equation For the tie-space doain dispersion-relation-based FD schee, spatial FD coefficients can be derived by using the following steps (Liu and Sen, 9b). First, derive dispersion relation of FD odeling by using plane wave theory and Equations ()-(4). Second, apply the Taylor series expansion for trigonoetric functions in the dispersion relation equation. Then, copare coefficients of power series and obtain the following equations ( cos θ cos φ cos θ sin φ sin θ ) + + a = r j j j j j j j ( j =,,..., ), (6) where, r = vτ / h, θ is the plane wave propagation angle easured fro the horizontal plane perpendicular to z axis, φ is the aziuth of the plane wave. Note that these equations indicate that the coefficients a are a function of θ and φ. We solve Equation (6) to obtain a by using an optial direction with θ = and φ = π /8 (Liu and Sen, 9b), which akes FD odeling attain the highest ()th-order accuracy along 48 directions. Absorbing boundary condition Here, we adopt the new hybrid schee developed by Liu and Sen (9c) to absorb reflections fro the odel boundaries in nuerical solutions of wave equations (WEs). This schee divides the coputational doain into boundary, transition and inner areas. The wavefields within the inner and the boundary areas are coputed by the WE and the one-way wave equation (OWWE) respectively. The wavefields within the transition area are deterined by a weighted cobination of the wavefields coputed by WE and the OWWE to obtain a sooth variation fro the inner area to the boundary via the transition zone. We develop this hybrid ABC for 3D acoustic wave odeling. For 3D odeling, boundaries include six sides, twelve edges and eight corners. Second-order OWWEs are adopted for these sides, and first-order OWWEs for these edges and corners (Clayton and Engquist, 977). Ten grid points are used to define the transition area to absorb boundary reflections. Dispersion analysis The dispersion of FD for 3D odeling is described by the following expression (Liu and Sen, 9b), ( kh θ ) ( ) ( kh θ φ ), + sin cos cos / (7) sin sin / vfd sin sin δ = = r a + kh cos θ sin φ / v rkh
3 3D Acoustic odeling with Dispersion-relation-based FD ethod Where, k is the wavenuber. If δ equals, there is no dispersion. If δ is far fro, a large dispersion will occur. Since kh is equal to π at the Nyquist frequency, when calculating δ, kh only ranges fro to π. Figure illustrates the variation of dispersion paraeter δ ( φ, θ ) with kh along nine directions, which deonstrates that the accuracy of the tie-space ethod is greater than that of the conventional ethod. v FD /v v FD /v kh (a) Conventional ethod kh Exact (, ) (, π/8) (, π/8) (π/8, π/8) (π/8, π/8) (π/8, 3π/8) (π/8, π/8) (π/8, π/8) (π/8, 3π/8) Exact (, ) (, π/8) (, π/8) (π/8, π/8) (π/8, π/8) (π/8, 3π/8) (π/8, π/8) (π/8, π/8) (π/8, 3π/8) (b) Tie-space ethod Figure : Plot of dispersion curves of the conventional and the tie-space FD ethods for 3D acoustic wave equation odeling. =, v = 3/s, τ =.s, h =. Nuerical odeling exaples First, both the conventional and the tie-space FD ethods are used to siulate 3D acoustic wave propagation in a hoogeneous acoustic ediu under the sae discretization. The odel and siulation paraeters are listed in the caption of Figure. Coputed snapshots respectively by the conventional and the tie-space FD ethods are shown in Figures (a) and (b). Coparing these two figures, we can see that the tie-space ethod aintains the wavefor better than the conventional ethod and thus has greater precision. To iprove the odeling accuracy further, we introduce the truncated FD ethod into the tie-space FD ethod. It is known that with the increase of the nuber of grid points involved in FD discretization, the accuracy increases but the coputational cost also increases. Liu and Sen (9d) found that there exist soe very sall coefficients for highorder FD coefficients and with the increase of order the nuber of these sall coefficients increases but their values decrease sharply. They also deonstrated that oitting these sall coefficients can aintain approxiately the sae level of accuracy of FD but reduce coputational cost significantly. Figure (c) shows snapshots by the truncated tie-space ethod, it follows that the dispersion decreases further copared with Figure (b). Next, we adopt the truncated tie-space FD ethod to perfor nuerical odeling. Figure 3 displays the seisogras coputed for a horizontally layered acoustic odel with the Clayton-Engquist ABC (Clayton and Engquist, 977) and the new hybrid ABC. The odel and siulation paraeters are given in the figure caption. The figure suggests that the boundary reflections are still strong copared with the reflections fro the true reflectors for the Clayton-Engquist ABC. The new hybrid ABC achieves nearly perfect absorption. Conclusions We have developed a tie-space doain dispersionrelation-based FD schee cobined with the truncated FD ethod and the hybrid ABC for 3D acoustic wave equation odeling. Dispersion analysis and nuerical odeling results deonstrate that this schee has greater accuracy and can effectively suppress dispersion and boundary reflections. Acknowledgents This research is partially supported by NSFC under contract No and the National 863 Progra of China under contract No. 7AA6Z8. 3
4 3D Acoustic odeling with Dispersion-relation-based FD ethod (a) Conventional ethod (a) Snapshots at t=5s by the Clayton-Engquist ABC (left) and (b) Tie-space ethod (b) Snapshots at t=7s by the Clayton-Engquist ABC (left) and (c) Snapshots at t=s by the Clayton-Engquist ABC (left) and 7 3 Trace no Tie (s) (c) Truncated tie-space ethod Figure : Snapshots coputed by finite-difference odeling for a 3D hoogeneous acoustic odel respectively with conventional, tie-space and truncated tie-space FD ethods. Each figure includes panels with snapshots at 5s and 5s fro left to right. The odel velocity is 3/s, grid size is, grid diensions are, grid point coordinates range fro (,,) to (99,99,99), tie step is s, =. The sae discretization is involved in all the three ethods, but their spatial FD coefficients are different. A source pulse of 5Hz sine function with one period length is located at the center of the odel. The free surface condition and ABC are not included here. Top, front and right surfaces of snapshots are recorded at z=, y=5 and x=99 respectively (d) Seisogras with the Clayton-Engquist ABC 4
5 3D Acoustic odeling with Dispersion-relation-based FD ethod Tie (s) Trace no (e) Seisogras with the hybrid ABC Figure 3: Snapshots and seisogras coputed by truncated tiespace FD odeling for 3D horizontally layered acoustic odel respectively with the Clayton-Engquist ABC and the hybrid ABC. Each figure of seisogras includes 3 receiver lines. The odel has 6 layers, whose velocities are velocity are 5/s, 3/s, 35/s, 3/s, 36/s and 4/s fro shallow to deep, 5 interfaces depth are 3, 5, 8, and 4. Grid size is, grid diensions are, grid coordinates range fro (,,) to (99,99,99), tie step is s. =. A source pulse of Hz sine function with one period length is located (99, 99, ). Receivers are located on the surface, and 3 receiver lines shown in this figure are located at y=99, 49 and 99 fro left to right. Trace no. of seisogra for each receiver line ranges fro to. Traces whose no. varies fro to 9 are shown since grids of width are used for the hybrid ABC. The free surface condition is included here. Top, front and right surfaces of snapshots are recorded at z=, y= and x=99 respectively. References Clayton, R. W., and B. Engquist, 977, absorbing boundary conditions for acoustic and elastic wave equations: Bulletin of the Seisological Society of Aerica, 6, Dablain,. A., 986, the application of high-order differencing to the scalar wave equation: Geophysics, 5, Etgen, J. T., and. J. O Brien, 7, Coputational ethods for large-scale 3D acoustic finite-difference odeling: A tutorial: Geophysics, 7, S3 S3. Fornberg B., 987, the pseudospectral ethod - coparisons with finite differences for the elastic wave equation: Geophysics, 5, Finkelstein B., and R. Kastner, 7, Finite difference tie doain dispersion reduction schees: Journal of Coputational Physics,, Kelly, K. R., R. Ward, W. S. Treitel, and R.. Alford, 976, Synthetic seisogras: A finite-difference approach: Geophysics, 4, 7. Liu Y., and. K. Sen, 9a, A practical iplicit finitedifference ethod: exaples fro seisic odeling: Journal of Geophysics and Engineering, 6, Liu Y., and. K. Sen, 9b, A new tie-space doain high-order finite-difference ethod for the acoustic wave equation, Journal of Coputational Physics, doi:.6/j.jcp Liu Y., and. K. Sen, 9c, A hybrid schee for absorbing edge reflections in nuerical odeling of wave propagation: Geophysics, in revision. Liu Y., and. K. Sen, 9d, Nuerical odeling of wave equation by a truncated high-order finite-difference ethod: Earthquake Science, ():
P032 3D Seismic Diffraction Modeling in Multilayered Media in Terms of Surface Integrals
P032 3D Seisic Diffraction Modeling in Multilayered Media in Ters of Surface Integrals A.M. Aizenberg (Institute of Geophysics SB RAS, M. Ayzenberg* (Norwegian University of Science & Technology, H.B.
More informationData-Driven Imaging in Anisotropic Media
18 th World Conference on Non destructive Testing, 16- April 1, Durban, South Africa Data-Driven Iaging in Anisotropic Media Arno VOLKER 1 and Alan HUNTER 1 TNO Stieltjesweg 1, 6 AD, Delft, The Netherlands
More informationNumerical Solution of the MRLW Equation Using Finite Difference Method. 1 Introduction
ISSN 1749-3889 print, 1749-3897 online International Journal of Nonlinear Science Vol.1401 No.3,pp.355-361 Nuerical Solution of the MRLW Equation Using Finite Difference Method Pınar Keskin, Dursun Irk
More informationINTRODUCTION. Residual migration has proved to be a useful tool in imaging and in velocity analysis.
Stanford Exploration Project, Report, June 3, 999, pages 5 59 Short Note On Stolt prestack residual igration Paul Sava keywords: Stolt, residual igration INTRODUCTION Residual igration has proved to be
More informationDESIGN OF THE DIE PROFILE FOR THE INCREMENTAL RADIAL FORGING PROCESS *
IJST, Transactions of Mechanical Engineering, Vol. 39, No. M1, pp 89-100 Printed in The Islaic Republic of Iran, 2015 Shira University DESIGN OF THE DIE PROFILE FOR THE INCREMENTAL RADIAL FORGING PROCESS
More informationExperimental Design For Model Discrimination And Precise Parameter Estimation In WDS Analysis
City University of New York (CUNY) CUNY Acadeic Works International Conference on Hydroinforatics 8-1-2014 Experiental Design For Model Discriination And Precise Paraeter Estiation In WDS Analysis Giovanna
More informationlecture 37: Linear Multistep Methods: Absolute Stability, Part I lecture 38: Linear Multistep Methods: Absolute Stability, Part II
lecture 37: Linear Multistep Methods: Absolute Stability, Part I lecture 3: Linear Multistep Methods: Absolute Stability, Part II 5.7 Linear ultistep ethods: absolute stability At this point, it ay well
More informationA Nonlinear Sparsity Promoting Formulation and Algorithm for Full Waveform Inversion
A Nonlinear Sparsity Prooting Forulation and Algorith for Full Wavefor Inversion Aleksandr Aravkin, Tristan van Leeuwen, Jaes V. Burke 2 and Felix Herrann Dept. of Earth and Ocean sciences University of
More informationProbabilistic micro-earthquake location for reservoir monitoring and characterization
CWP-63 Probabilistic icro-earthquake location for reservoir onitoring and characterization Ran Xuan and Paul Sava Center for Wave Phenoena, Colorado School of Mines ABSTRACT Micro-seisicity is used to
More informationComparison of Stability of Selected Numerical Methods for Solving Stiff Semi- Linear Differential Equations
International Journal of Applied Science and Technology Vol. 7, No. 3, Septeber 217 Coparison of Stability of Selected Nuerical Methods for Solving Stiff Sei- Linear Differential Equations Kwaku Darkwah
More informationNUMERICAL MODELLING OF THE TYRE/ROAD CONTACT
NUMERICAL MODELLING OF THE TYRE/ROAD CONTACT PACS REFERENCE: 43.5.LJ Krister Larsson Departent of Applied Acoustics Chalers University of Technology SE-412 96 Sweden Tel: +46 ()31 772 22 Fax: +46 ()31
More informationNumerical Studies of a Nonlinear Heat Equation with Square Root Reaction Term
Nuerical Studies of a Nonlinear Heat Equation with Square Root Reaction Ter Ron Bucire, 1 Karl McMurtry, 1 Ronald E. Micens 2 1 Matheatics Departent, Occidental College, Los Angeles, California 90041 2
More informatione-companion ONLY AVAILABLE IN ELECTRONIC FORM
OPERATIONS RESEARCH doi 10.1287/opre.1070.0427ec pp. ec1 ec5 e-copanion ONLY AVAILABLE IN ELECTRONIC FORM infors 07 INFORMS Electronic Copanion A Learning Approach for Interactive Marketing to a Custoer
More informationFluid Substitution Model to Generate Synthetic Seismic Attributes: FluidSub.exe
Fluid Substitution Model to Generate Synthetic Seisic Attributes: FluidSub.exe Sahyun Hong and Clayton V. Deutsch Geostatistical data integration is a ature research field and any of algoriths have been
More informationInspection; structural health monitoring; reliability; Bayesian analysis; updating; decision analysis; value of information
Cite as: Straub D. (2014). Value of inforation analysis with structural reliability ethods. Structural Safety, 49: 75-86. Value of Inforation Analysis with Structural Reliability Methods Daniel Straub
More informationBernoulli Wavelet Based Numerical Method for Solving Fredholm Integral Equations of the Second Kind
ISSN 746-7659, England, UK Journal of Inforation and Coputing Science Vol., No., 6, pp.-9 Bernoulli Wavelet Based Nuerical Method for Solving Fredhol Integral Equations of the Second Kind S. C. Shiralashetti*,
More informationThe Wilson Model of Cortical Neurons Richard B. Wells
The Wilson Model of Cortical Neurons Richard B. Wells I. Refineents on the odgkin-uxley Model The years since odgkin s and uxley s pioneering work have produced a nuber of derivative odgkin-uxley-like
More informationSPATIAL AND TEMPORAL DISTRIBUTIONS OF MIGRATION IN BIO-RETENTION SYSTEMS
Li, Y.-Q., et al.: Spatial and Teporal Distributions of Moisture THERMAL SCIENCE, Year 04, Vol. 8, No. 5, pp. 557-56 557 SPATIAL AND TEMPORAL DISTRIBUTIONS OF MIGRATION IN BIO-RETENTION SYSTEMS by Yu-Qi
More informationComparison of Charged Particle Tracking Methods for Non-Uniform Magnetic Fields. Hann-Shin Mao and Richard E. Wirz
42nd AIAA Plasadynaics and Lasers Conferencein conjunction with the8th Internati 27-30 June 20, Honolulu, Hawaii AIAA 20-3739 Coparison of Charged Particle Tracking Methods for Non-Unifor Magnetic
More informationCOS 424: Interacting with Data. Written Exercises
COS 424: Interacting with Data Hoework #4 Spring 2007 Regression Due: Wednesday, April 18 Written Exercises See the course website for iportant inforation about collaboration and late policies, as well
More informationAnalytical solution for the electric potential in arbitrary anisotropic layered media applying the set of Hankel transforms of integer order
Geophysical Prospecting, 2006, 54, 651 661 Analytical solution for the electric potential in arbitrary anisotropic layered edia applying the set of Hankel transfors of integer order E. Pervago 1,A.Mousatov
More informationP016 Toward Gauss-Newton and Exact Newton Optimization for Full Waveform Inversion
P016 Toward Gauss-Newton and Exact Newton Optiization for Full Wavefor Inversion L. Métivier* ISTerre, R. Brossier ISTerre, J. Virieux ISTerre & S. Operto Géoazur SUMMARY Full Wavefor Inversion FWI applications
More informationInterferometry in dissipative media: Addressing the shallow sea problem for Seabed Logging applications
Evert Slob and Kees Wapenaar, Departent of Geotechnology, Delft University of Technology Roel Snieder, Center for Wave Phenoena and Departent of Geophysics, Colorado School of Mines SUMMARY We derive interferoetric
More informationProc. of the IEEE/OES Seventh Working Conference on Current Measurement Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES
Proc. of the IEEE/OES Seventh Working Conference on Current Measureent Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES Belinda Lipa Codar Ocean Sensors 15 La Sandra Way, Portola Valley, CA 98 blipa@pogo.co
More informationUsing a De-Convolution Window for Operating Modal Analysis
Using a De-Convolution Window for Operating Modal Analysis Brian Schwarz Vibrant Technology, Inc. Scotts Valley, CA Mark Richardson Vibrant Technology, Inc. Scotts Valley, CA Abstract Operating Modal Analysis
More informationMulti-Scale/Multi-Resolution: Wavelet Transform
Multi-Scale/Multi-Resolution: Wavelet Transfor Proble with Fourier Fourier analysis -- breaks down a signal into constituent sinusoids of different frequencies. A serious drawback in transforing to the
More informationSEISMIC FRAGILITY ANALYSIS
9 th ASCE Specialty Conference on Probabilistic Mechanics and Structural Reliability PMC24 SEISMIC FRAGILITY ANALYSIS C. Kafali, Student M. ASCE Cornell University, Ithaca, NY 483 ck22@cornell.edu M. Grigoriu,
More informationNon-Parametric Non-Line-of-Sight Identification 1
Non-Paraetric Non-Line-of-Sight Identification Sinan Gezici, Hisashi Kobayashi and H. Vincent Poor Departent of Electrical Engineering School of Engineering and Applied Science Princeton University, Princeton,
More informationREPRODUCING KERNEL PARTICLE METHOD FOR CONDUCTION RADIATION INTERACTION IN 3D PARTICIPATING MEDIA
Proceedings of the Asian Conference on Theral Sciences 2017, 1st ACTS March 26-30, 2017, Jeu Island, Korea ACTS-P00337 REPRODUCING KERNEL PARTICLE METHOD FOR CONDUCTION RADIATION INTERACTION IN 3D PARTICIPATING
More informationNew Slack-Monotonic Schedulability Analysis of Real-Time Tasks on Multiprocessors
New Slack-Monotonic Schedulability Analysis of Real-Tie Tasks on Multiprocessors Risat Mahud Pathan and Jan Jonsson Chalers University of Technology SE-41 96, Göteborg, Sweden {risat, janjo}@chalers.se
More informationABSTRACT INTRODUCTION
Wave Resistance Prediction of a Cataaran by Linearised Theory M.INSEL Faculty of Naval Architecture and Ocean Engineering, Istanbul Technical University, TURKEY A.F.MOLLAND, J.F.WELLICOME Departent of
More informationANALYSIS ON RESPONSE OF DYNAMIC SYSTEMS TO PULSE SEQUENCES EXCITATION
The 4 th World Conference on Earthquake Engineering October -7, 8, Beijing, China ANALYSIS ON RESPONSE OF DYNAMIC SYSTEMS TO PULSE SEQUENCES EXCITATION S. Li C.H. Zhai L.L. Xie Ph. D. Student, School of
More informationAN EXPLICIT METHOD FOR NUMERICAL SIMULATION OF WAVE EQUATIONS
The 4 th World Conference on Earthquake Engineering October -7, 8, Beiing, China AN EXPLICIT ETHOD FOR NUERICAL SIULATION OF WAVE EQUATIONS Liu Heng and Liao Zheneng Doctoral Candidate, Det. of Structural
More informationStudy of radiation in spherical media using moment method
Aerican Journal of Physics and Applications 03; (): 5-3 Published online July 0, 03 (http://www.sciencepublishinggroup.co/j/ajpa) doi: 0.648/j.ajpa.0300.5 Study of radiation in spherical edia using oent
More informationDynamic Response of SDOF System: A Comparative Study Between Newmark s Method and IS1893 Response Spectra
Dynaic Response of SDOF Syste: A Coparative Study Between Newark s Method and IS893 Response Spectra [] P. Manoj Reddy [2] Dr. A. Viala [] Post Graduate Student, [2] Associate Professor [][2] Chaitanya
More information16.333: Lecture # 7. Approximate Longitudinal Dynamics Models
16.333: Lecture # 7 Approxiate Longitudinal Dynaics Models A couple ore stability derivatives Given ode shapes found identify sipler odels that capture the ain responses Fall 24 16.333 6 1 More Stability
More informationCh 12: Variations on Backpropagation
Ch 2: Variations on Backpropagation The basic backpropagation algorith is too slow for ost practical applications. It ay take days or weeks of coputer tie. We deonstrate why the backpropagation algorith
More informationCHAPTER 19: Single-Loop IMC Control
When I coplete this chapter, I want to be able to do the following. Recognize that other feedback algoriths are possible Understand the IMC structure and how it provides the essential control features
More informationExtension of CSRSM for the Parametric Study of the Face Stability of Pressurized Tunnels
Extension of CSRSM for the Paraetric Study of the Face Stability of Pressurized Tunnels Guilhe Mollon 1, Daniel Dias 2, and Abdul-Haid Soubra 3, M.ASCE 1 LGCIE, INSA Lyon, Université de Lyon, Doaine scientifique
More informationPh 20.3 Numerical Solution of Ordinary Differential Equations
Ph 20.3 Nuerical Solution of Ordinary Differential Equations Due: Week 5 -v20170314- This Assignent So far, your assignents have tried to failiarize you with the hardware and software in the Physics Coputing
More informationREDUCTION OF FINITE ELEMENT MODELS BY PARAMETER IDENTIFICATION
ISSN 139 14X INFORMATION TECHNOLOGY AND CONTROL, 008, Vol.37, No.3 REDUCTION OF FINITE ELEMENT MODELS BY PARAMETER IDENTIFICATION Riantas Barauskas, Vidantas Riavičius Departent of Syste Analysis, Kaunas
More informationAnisotropic reference media and the possible linearized approximations for phase velocities of qs waves in weakly anisotropic media
INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF PHYSICS D: APPLIED PHYSICS J. Phys. D: Appl. Phys. 5 00 007 04 PII: S00-770867-6 Anisotropic reference edia and the possible linearized approxiations for phase
More informationTheory and Applications of the Indoor Noise Module (VDI 3760)
Fullerton, CA 9835 USA Theory and Applications of the Indoor Noise Module (VDI 3760) The following Technical Note suarizes the atheatical concepts used in the SoundPLAN Indoor Factory Noise Module. 1.
More informationANALYTICAL INVESTIGATION AND PARAMETRIC STUDY OF LATERAL IMPACT BEHAVIOR OF PRESSURIZED PIPELINES AND INFLUENCE OF INTERNAL PRESSURE
DRAFT Proceedings of the ASME 014 International Mechanical Engineering Congress & Exposition IMECE014 Noveber 14-0, 014, Montreal, Quebec, Canada IMECE014-36371 ANALYTICAL INVESTIGATION AND PARAMETRIC
More informationA remark on a success rate model for DPA and CPA
A reark on a success rate odel for DPA and CPA A. Wieers, BSI Version 0.5 andreas.wieers@bsi.bund.de Septeber 5, 2018 Abstract The success rate is the ost coon evaluation etric for easuring the perforance
More informationA Self-Organizing Model for Logical Regression Jerry Farlow 1 University of Maine. (1900 words)
1 A Self-Organizing Model for Logical Regression Jerry Farlow 1 University of Maine (1900 words) Contact: Jerry Farlow Dept of Matheatics Univeristy of Maine Orono, ME 04469 Tel (07) 866-3540 Eail: farlow@ath.uaine.edu
More informationQ ESTIMATION WITHIN A FORMATION PROGRAM q_estimation
Foration Attributes Progra q_estiation Q ESTIMATION WITHIN A FOMATION POGAM q_estiation Estiating Q between stratal slices Progra q_estiation estiate seisic attenuation (1/Q) on coplex stratal slices using
More informationNumerical Study of Tsunami Propagation in Mentawai Islands West Sumatra
Universal Journal of Geoscience 5(4): 112-116, 201 DOI: 10.13189/ujg.201.050404 http://www.hrpub.org Nuerical Study of Tsunai Propagation in Mentawai Islands West Suatra Noverina Alfiany *, Kazuhiro Yaaoto,
More informationAn Alternative Discretization and Solution Procedure for the Dual Phase-Lag Equation
An Alternative Discretization and Solution Procedure for the Dual Phase-Lag Equation J. M. McDonough, Departents of Mechanical Engineering and Matheatics, University of Kentucky, Lexington, Kentucky 40506-0503,
More informationSome Perspective. Forces and Newton s Laws
Soe Perspective The language of Kineatics provides us with an efficient ethod for describing the otion of aterial objects, and we ll continue to ake refineents to it as we introduce additional types of
More informationAnalyzing Simulation Results
Analyzing Siulation Results Dr. John Mellor-Cruey Departent of Coputer Science Rice University johnc@cs.rice.edu COMP 528 Lecture 20 31 March 2005 Topics for Today Model verification Model validation Transient
More informationA Boundary Element Model for Underwater Acoustics in Shallow Water
A Boundary Eleent Model for Underwater Acoustics in hallow Water J.A.F. antiago,l.c.wrobel Abstract: This work presents a boundary eleent forulation for two-diensional acoustic wave propagation in shallow
More informationPOD-DEIM MODEL ORDER REDUCTION FOR THE MONODOMAIN REACTION-DIFFUSION EQUATION IN NEURO-MUSCULAR SYSTEM
6th European Conference on Coputational Mechanics (ECCM 6) 7th European Conference on Coputational Fluid Dynaics (ECFD 7) 1115 June 2018, Glasgow, UK POD-DEIM MODEL ORDER REDUCTION FOR THE MONODOMAIN REACTION-DIFFUSION
More informationMAKE SURE TA & TI STAMPS EVERY PAGE BEFORE YOU START
Laboratory Section: Last Revised on Deceber 15, 2014 Partners Naes: Grade: EXPERIMENT 8 Electron Beas 0. Pre-Laboratory Work [2 pts] 1. Nae the 2 forces that are equated in order to derive the charge to
More informationModel Fitting. CURM Background Material, Fall 2014 Dr. Doreen De Leon
Model Fitting CURM Background Material, Fall 014 Dr. Doreen De Leon 1 Introduction Given a set of data points, we often want to fit a selected odel or type to the data (e.g., we suspect an exponential
More informationResearch Article Rapidly-Converging Series Representations of a Mutual-Information Integral
International Scholarly Research Network ISRN Counications and Networking Volue 11, Article ID 5465, 6 pages doi:1.54/11/5465 Research Article Rapidly-Converging Series Representations of a Mutual-Inforation
More informationA Simplified Analytical Approach for Efficiency Evaluation of the Weaving Machines with Automatic Filling Repair
Proceedings of the 6th SEAS International Conference on Siulation, Modelling and Optiization, Lisbon, Portugal, Septeber -4, 006 0 A Siplified Analytical Approach for Efficiency Evaluation of the eaving
More informationANALYSIS OF HALL-EFFECT THRUSTERS AND ION ENGINES FOR EARTH-TO-MOON TRANSFER
IEPC 003-0034 ANALYSIS OF HALL-EFFECT THRUSTERS AND ION ENGINES FOR EARTH-TO-MOON TRANSFER A. Bober, M. Guelan Asher Space Research Institute, Technion-Israel Institute of Technology, 3000 Haifa, Israel
More informationProjectile Motion with Air Resistance (Numerical Modeling, Euler s Method)
Projectile Motion with Air Resistance (Nuerical Modeling, Euler s Method) Theory Euler s ethod is a siple way to approxiate the solution of ordinary differential equations (ode s) nuerically. Specifically,
More informationRecovering Data from Underdetermined Quadratic Measurements (CS 229a Project: Final Writeup)
Recovering Data fro Underdeterined Quadratic Measureents (CS 229a Project: Final Writeup) Mahdi Soltanolkotabi Deceber 16, 2011 1 Introduction Data that arises fro engineering applications often contains
More informationIntelligent Systems: Reasoning and Recognition. Perceptrons and Support Vector Machines
Intelligent Systes: Reasoning and Recognition Jaes L. Crowley osig 1 Winter Seester 2018 Lesson 6 27 February 2018 Outline Perceptrons and Support Vector achines Notation...2 Linear odels...3 Lines, Planes
More informationGrid performance models using Design of Experiments (DoE) methods
Grid perforance odels using Design of Experients (DoE) ethods IEPC 25-139 Presented at the 29 th International Electric Propulsion Conference, Princeton University, October 31 Noveber, 25 Michael Tartz
More informationAn improved self-adaptive harmony search algorithm for joint replenishment problems
An iproved self-adaptive harony search algorith for joint replenishent probles Lin Wang School of Manageent, Huazhong University of Science & Technology zhoulearner@gail.co Xiaojian Zhou School of Manageent,
More informationA Better Algorithm For an Ancient Scheduling Problem. David R. Karger Steven J. Phillips Eric Torng. Department of Computer Science
A Better Algorith For an Ancient Scheduling Proble David R. Karger Steven J. Phillips Eric Torng Departent of Coputer Science Stanford University Stanford, CA 9435-4 Abstract One of the oldest and siplest
More informationDepartment of Electronic and Optical Engineering, Ordnance Engineering College, Shijiazhuang, , China
6th International Conference on Machinery, Materials, Environent, Biotechnology and Coputer (MMEBC 06) Solving Multi-Sensor Multi-Target Assignent Proble Based on Copositive Cobat Efficiency and QPSO Algorith
More informationMulti-Window Post-Doppler Dimensionality Reduction for Multi-Waveform STAP
Multi-Window Post-Doppler Diensionality Reduction for Multi-Wavefor STAP Luuba arnett 1, Patrick M. McCorick 1, Shannon D. Blunt 1, Justin G. Metcalf 2 1 Radar Systes Lab, University of Kansas, Lawrence,
More informationŞtefan ŞTEFĂNESCU * is the minimum global value for the function h (x)
7Applying Nelder Mead s Optiization Algorith APPLYING NELDER MEAD S OPTIMIZATION ALGORITHM FOR MULTIPLE GLOBAL MINIMA Abstract Ştefan ŞTEFĂNESCU * The iterative deterinistic optiization ethod could not
More informationSuppress Parameter Cross-talk for Elastic Full-waveform Inversion: Parameterization and Acquisition Geometry
Suress Paraeter Cross-talk for Elastic Full-wavefor Inversion: Paraeterization and Acquisition Geoetry Wenyong Pan and Kris Innanen CREWES Project, Deartent of Geoscience, University of Calgary Suary Full-wavefor
More informationQualitative Modelling of Time Series Using Self-Organizing Maps: Application to Animal Science
Proceedings of the 6th WSEAS International Conference on Applied Coputer Science, Tenerife, Canary Islands, Spain, Deceber 16-18, 2006 183 Qualitative Modelling of Tie Series Using Self-Organizing Maps:
More informationN-Point. DFTs of Two Length-N Real Sequences
Coputation of the DFT of In ost practical applications, sequences of interest are real In such cases, the syetry properties of the DFT given in Table 5. can be exploited to ake the DFT coputations ore
More informationStochastic Subgradient Methods
Stochastic Subgradient Methods Lingjie Weng Yutian Chen Bren School of Inforation and Coputer Science University of California, Irvine {wengl, yutianc}@ics.uci.edu Abstract Stochastic subgradient ethods
More informationOcean 420 Physical Processes in the Ocean Project 1: Hydrostatic Balance, Advection and Diffusion Answers
Ocean 40 Physical Processes in the Ocean Project 1: Hydrostatic Balance, Advection and Diffusion Answers 1. Hydrostatic Balance a) Set all of the levels on one of the coluns to the lowest possible density.
More informationAdapting the Pheromone Evaporation Rate in Dynamic Routing Problems
Adapting the Pheroone Evaporation Rate in Dynaic Routing Probles Michalis Mavrovouniotis and Shengxiang Yang School of Coputer Science and Inforatics, De Montfort University The Gateway, Leicester LE1
More informationSupport Vector Machine Classification of Uncertain and Imbalanced data using Robust Optimization
Recent Researches in Coputer Science Support Vector Machine Classification of Uncertain and Ibalanced data using Robust Optiization RAGHAV PAT, THEODORE B. TRAFALIS, KASH BARKER School of Industrial Engineering
More informationA DESIGN GUIDE OF DOUBLE-LAYER CELLULAR CLADDINGS FOR BLAST ALLEVIATION
International Journal of Aerospace and Lightweight Structures Vol. 3, No. 1 (2013) 109 133 c Research Publishing Services DOI: 10.3850/S201042862013000550 A DESIGN GUIDE OF DOUBLE-LAYER CELLULAR CLADDINGS
More informationAnton Bourdine. 1. Introduction. and approximate propagation constants by following simple ratio (Equation (32.22) in [1]):
Matheatical Probles in Engineering olue 5, Article ID 843, pages http://dx.doi.org/.55/5/843 Research Article Fast and Siple Method for Evaluation of Polarization Correction to Propagation Constant of
More informationFeature Extraction Techniques
Feature Extraction Techniques Unsupervised Learning II Feature Extraction Unsupervised ethods can also be used to find features which can be useful for categorization. There are unsupervised ethods that
More informationAnalysis of ground vibration transmission in high precision equipment by Frequency Based Substructuring
Analysis of ground vibration transission in high precision equipent by Frequency Based Substructuring G. van Schothorst 1, M.A. Boogaard 2, G.W. van der Poel 1, D.J. Rixen 2 1 Philips Innovation Services,
More informationAn Improved Particle Filter with Applications in Ballistic Target Tracking
Sensors & ransducers Vol. 72 Issue 6 June 204 pp. 96-20 Sensors & ransducers 204 by IFSA Publishing S. L. http://www.sensorsportal.co An Iproved Particle Filter with Applications in Ballistic arget racing
More informationStability Analysis of the Matrix-Free Linearly Implicit 2 Euler Method 3 UNCORRECTED PROOF
1 Stability Analysis of the Matrix-Free Linearly Iplicit 2 Euler Method 3 Adrian Sandu 1 andaikst-cyr 2 4 1 Coputational Science Laboratory, Departent of Coputer Science, Virginia 5 Polytechnic Institute,
More informationThis model assumes that the probability of a gap has size i is proportional to 1/i. i.e., i log m e. j=1. E[gap size] = i P r(i) = N f t.
CS 493: Algoriths for Massive Data Sets Feb 2, 2002 Local Models, Bloo Filter Scribe: Qin Lv Local Models In global odels, every inverted file entry is copressed with the sae odel. This work wells when
More informationNumerically repeated support splitting and merging phenomena in a porous media equation with strong absorption. Kenji Tomoeda
Journal of Math-for-Industry, Vol. 3 (C-), pp. Nuerically repeated support splitting and erging phenoena in a porous edia equation with strong absorption To the eory of y friend Professor Nakaki. Kenji
More informationSPECTRUM sensing is a core concept of cognitive radio
World Acadey of Science, Engineering and Technology International Journal of Electronics and Counication Engineering Vol:6, o:2, 202 Efficient Detection Using Sequential Probability Ratio Test in Mobile
More informationEfficient dynamic events discrimination technique for fiber distributed Brillouin sensors
Efficient dynaic events discriination technique for fiber distributed rillouin sensors Carlos A. Galindez,* Francisco J. Madruga, and Jose M. Lopez-Higuera Photonics Engineering Group, Universidad de Cantabria,Edif.
More informationCelal S. Konor Release 1.1 (identical to 1.0) 3/21/08. 1-Hybrid isentropic-sigma vertical coordinate and governing equations in the free atmosphere
Celal S. Konor Release. (identical to.0) 3/2/08 -Hybrid isentropic-siga vertical coordinate governing equations in the free atosphere This section describes the equations in the free atosphere of the odel.
More informationSpine Fin Efficiency A Three Sided Pyramidal Fin of Equilateral Triangular Cross-Sectional Area
Proceedings of the 006 WSEAS/IASME International Conference on Heat and Mass Transfer, Miai, Florida, USA, January 18-0, 006 (pp13-18) Spine Fin Efficiency A Three Sided Pyraidal Fin of Equilateral Triangular
More informationRole of Control-Structure Interaction in Protective System Design
Published in the ASCE Journal of Engineering Mechanics, ol. 121, No. 2, Feb. 1995 pp. 322 338. Role of Control-Structure Interaction in Protective Syste Design S. J. Dyke 1, Student Meber ASCE, B. F. Spencer,
More informationAn Approximate Model for the Theoretical Prediction of the Velocity Increase in the Intermediate Ballistics Period
An Approxiate Model for the Theoretical Prediction of the Velocity... 77 Central European Journal of Energetic Materials, 205, 2(), 77-88 ISSN 2353-843 An Approxiate Model for the Theoretical Prediction
More informationSensorless Control of Induction Motor Drive Using SVPWM - MRAS Speed Observer
Journal of Eerging Trends in Engineering and Applied Sciences (JETEAS) 2 (3): 509-513 Journal Scholarlink of Eerging Research Trends Institute in Engineering Journals, 2011 and Applied (ISSN: 2141-7016)
More informationDirect numerical simulation of matrix diffusion across fracture/matrix interface
Water Science and Engineering, 2013, 6(4): 365-379 doi:10.3882/j.issn.1674-2370.2013.04.001 http://www.waterjournal.cn e-ail: wse2008@vip.163.co Direct nuerical siulation of atrix diffusion across fracture/atrix
More informationPattern Recognition and Machine Learning. Learning and Evaluation for Pattern Recognition
Pattern Recognition and Machine Learning Jaes L. Crowley ENSIMAG 3 - MMIS Fall Seester 2017 Lesson 1 4 October 2017 Outline Learning and Evaluation for Pattern Recognition Notation...2 1. The Pattern Recognition
More informationOptimization of ripple filter for pencil beam scanning
Nuclear Science and Techniques 24 (2013) 060404 Optiization of ripple filter for pencil bea scanning YANG Zhaoxia ZHANG Manzhou LI Deing * Shanghai Institute of Applied Physics, Chinese Acadey of Sciences,
More informationBootstrapping Dependent Data
Bootstrapping Dependent Data One of the key issues confronting bootstrap resapling approxiations is how to deal with dependent data. Consider a sequence fx t g n t= of dependent rando variables. Clearly
More informationA NEW ROBUST AND EFFICIENT ESTIMATOR FOR ILL-CONDITIONED LINEAR INVERSE PROBLEMS WITH OUTLIERS
A NEW ROBUST AND EFFICIENT ESTIMATOR FOR ILL-CONDITIONED LINEAR INVERSE PROBLEMS WITH OUTLIERS Marta Martinez-Caara 1, Michael Mua 2, Abdelhak M. Zoubir 2, Martin Vetterli 1 1 School of Coputer and Counication
More informationEffective joint probabilistic data association using maximum a posteriori estimates of target states
Effective joint probabilistic data association using axiu a posteriori estiates of target states 1 Viji Paul Panakkal, 2 Rajbabu Velurugan 1 Central Research Laboratory, Bharat Electronics Ltd., Bangalore,
More informationGREY FORECASTING AND NEURAL NETWORK MODEL OF SPORT PERFORMANCE
Journal of heoretical and Applied Inforation echnology 3 st March 03 Vol 49 No3 005-03 JAI & LLS All rights reserved ISSN: 99-8645 wwwatitorg E-ISSN: 87-395 GREY FORECASING AND NEURAL NEWORK MODEL OF SPOR
More informationOn Poset Merging. 1 Introduction. Peter Chen Guoli Ding Steve Seiden. Keywords: Merging, Partial Order, Lower Bounds. AMS Classification: 68W40
On Poset Merging Peter Chen Guoli Ding Steve Seiden Abstract We consider the follow poset erging proble: Let X and Y be two subsets of a partially ordered set S. Given coplete inforation about the ordering
More informationA model reduction approach to numerical inversion for a parabolic partial differential equation
Inverse Probles Inverse Probles 30 (204) 250 (33pp) doi:0.088/0266-56/30/2/250 A odel reduction approach to nuerical inversion for a parabolic partial differential equation Liliana Borcea, Vladiir Drusin
More informationHyperbolic Horn Helical Mass Spectrometer (3HMS) James G. Hagerman Hagerman Technology LLC & Pacific Environmental Technologies April 2005
Hyperbolic Horn Helical Mass Spectroeter (3HMS) Jaes G Hageran Hageran Technology LLC & Pacific Environental Technologies April 5 ABSTRACT This paper describes a new type of ass filter based on the REFIMS
More informationANALYSIS OF REFLECTOR AND HORN ANTENNAS USING MULTILEVEL FAST MULTIPOLE ALGORITHM
European Congress on Coputational Methods in Applied Sciences and Engineering ECCOMAS 2 Barcelona, 11-14 Septeber 2 ECCOMAS ANALYSIS OF REFLECTOR AND HORN ANTENNAS USING MULTILEVEL FAST MULTIPOLE ALGORITHM
More information