A5056180795- earthquake and tectonic studies
- Advanced Numerical Methods in Computational Mathematics
- Computational Fluid Dynamics and Aerodynamics
- Numerical methods for differential equations
- Meteorological Phenomena and Simulations
- Seismology and Earthquake Studies
- Fluid Dynamics and Turbulent Flows
- Electromagnetic Simulation and Numerical Methods
- Seismic Waves and Analysis
- Geotechnical and Geomechanical Engineering
- High-pressure geophysics and materials
- Seismic Imaging and Inversion Techniques
- Climate variability and models
- Groundwater flow and contamination studies
- Earthquake Detection and Analysis
- Reservoir Engineering and Simulation Methods
- Landslides and related hazards
- Differential Equations and Numerical Methods
- Lattice Boltzmann Simulation Studies
- Numerical methods in engineering
- Hydraulic Fracturing and Reservoir Analysis
- Atmospheric and Environmental Gas Dynamics
- Earthquake and Tsunami Effects
- Geophysics and Sensor Technology
- Matrix Theory and Algorithms
Naval Postgraduate School
2014-2024
Stanford University
2007-2013
Grantmakers for Effective Organizations
2012
University of California, Santa Cruz
2004
University of Washington
2004
We study dynamic rupture propagation on flat faults using 2D plane strain models featuring strongly rate-weakening fault friction (in a rate-and-state framework) and off-fault Drucker-Prager viscoplasticity.Plastic deformation bounds stresses near the front limits slip velocities to ∼10 m=s, bound expected be independent of earthquake magnitude.As originally shown for ruptures in an elastic medium (Zheng Rice, 1998), consequence is existence critical background stress level at which...
Research Article| April 04, 2018 A Suite of Exercises for Verifying Dynamic Earthquake Rupture Codes Ruth A. Harris; Harris aU.S. Geological Survey, 345 Middlefield Road, MS 977, Menlo Park, California 94025 U.S.A., harris@usgs.gov Search other works by this author on: GSW Google Scholar Michael Barall; Barall bInvisible Software Inc., P.O. Box 6541, San Jose, 95150 U.S.A. Brad Aagaard; Aagaard Shuo Ma; Ma cDepartment Sciences, Diego State University, Diego, 92182 Daniel Roten; Roten Kim...
Research Article| May 01, 2013 Rupture to the Trench: Dynamic Simulations of 11 March 2011 Tohoku Earthquake Jeremy E. Kozdon; Kozdon Department Geophysics, Stanford University, 397 Panama Mall, Stanford, California 94305jekozdon@nps.edu *Now at Applied Mathematics, Naval Postgraduate School, 833 Dyer Road, Monterey, 93943. Search for other works by this author on: GSW Google Scholar Eric M. Dunham Geophysics and Institute Computational Mathematical Engineering, 94305edunham@stanford.edu...
Abstract Numerical simulations of sequences earthquakes and aseismic slip (SEAS) have made great progress over past decades to address important questions in earthquake physics. However, significant challenges SEAS modeling remain resolving multiscale interactions between nucleation, dynamic rupture, slip, understanding physical factors controlling observables such as seismicity ground deformation. The increasing complexity calls for extensive efforts verify codes advance these with rigor,...
ABSTRACT Numerical modeling of earthquake dynamics and derived insight for seismic hazard relies on credible, reproducible model results. The sequences earthquakes aseismic slip (SEAS) initiative has set out to facilitate community code comparisons, verify advance the next generation physics-based models that reproduce all phases cycle. With goal advancing SEAS robustly incorporate physical geometrical complexities, here we present comparison results from two new benchmark problems: BP1-FD...
A methodology for handling block-to-block coupling of nonconforming, multiblock summation-by-parts finite difference methods is proposed. The based on the construction projection operators that move a grid solution along an interface to space piecewise defined functions; we specifically consider discontinuous, polynomial functions. constructed are compatible with underlying energy norm. Using linear wave equation in two dimensions as model problem, stability coupled numerical method proven...
Abstract. We introduce ClimateMachine, a new open-source atmosphere modeling framework which uses the Julia language and is designed to be scalable on central processing units (CPUs) graphics (GPUs). ClimateMachine common both for coarser-resolution global simulations high-resolution, limited-area large-eddy (LESs). Here, we demonstrate LES configuration of model in canonical benchmark cases atmospheric flows using total energy-conserving nodal discontinuous Galerkin (DG) discretization...
Abstract Dynamical cores used to study the circulation of atmosphere employ various numerical methods ranging from finite‐volume, spectral element, global spectral, and hybrid methods. In this work, we explore use Flux‐Differencing Discontinuous Galerkin (FDDG) simulate a fully compressible dry at resolutions. We show that method offers judicious compromise between high‐order accuracy stability for large‐eddy simulations atmospheric general circulation. particular, filters, divergence...
Abstract We couple a node-centered finite volume method to high order difference simulate dynamic earthquake ruptures along nonplanar faults in two dimensions. The is implemented on an unstructured mesh, providing the ability handle complex geometries. geometric complexities are limited small portion of overall domain and elsewhere used, enhancing efficiency. Both methods summation-by-parts form. Interface conditions coupling numerical solution across physical interfaces like faults,...
Abstract We present an efficient numerical method for earthquake sequences in 2D antiplane shear that incorporates wave propagation. A vertical strike‐slip fault governed by rate‐and‐state friction is embedded a heterogeneous elastic half‐space discretized using high‐order accurate Summation‐by‐Parts finite difference method. Adaptive time‐stepping applied during the interseismic periods; coseismic rupture we apply non‐stiff method, enabling variety of explicit time stepping methods....
Modern analysis of X-ray absorption fine structure (XAFS) is usually based on a traditional least-squares fitting procedure. Here an alternative Bayes–Turchin method discussed which has number advantages. In particular the takes advantage priori estimates model parameters and their uncertainties avoids restriction size parameter space or necessity for Fourier filtering. Thus permits full spectra (XAS), including both XAFS near-edge (XANES). The approach leads to set linear equations...
Summary Multidimensional transport for reservoir simulation is typically solved by applying 1D numerical methods in each spatial-coordinate direction. This approach simple, but the disadvantage that errors become highly correlated with underlying computational grid. In many real-field applications, this can result strong sensitivity to grid design not only computed saturation/composition fields also critical integrated data such as breakthrough times. Therefore, increase robustness of...
Numerical modeling of earthquake dynamics and derived insight for seismic hazard relies on credible, reproducible model results. The SEAS (Sequences Earthquakes Aseismic Slip) initiative has set out to facilitate community code comparisons, verify advance the next generation physics-based models that reproduce all phases cycle. With goal advancing robustly incorporate physical geometrical complexities, here we present comparison results from two new benchmark problems: BP1-FD considers full...