Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “International Journal for Numerical Methods in Fluids”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

7 recordsLinked to original sources
USGS2005

Semi-implicit finite difference methods for three-dimensional shallow water flow

A semi-implicit finite difference method for the numerical solution of three-dimensional shallow water flows is presented and discussed. The governing equations are the primitive three-dimensional turbulent mean flow equations where the pressure distribution in the vertical has been assumed to be hydrostatic. In the method of solution a minimal degree of implicitness has been adopted in such a fashion that the resulting algorithm is stable and gives a maximal computational efficiency at a minimal computational cost. At each time step the numerical method requires the solution of one large linear system which can be formally decomposed into a set of small three-diagonal systems coupled with one five-diagonal system. All these linear systems are symmetric and positive definite. Thus the existence and uniquencess of the numerical solution are assured. When only one vertical layer is specified, this method reduces as a special case to a semi-implicit scheme for solving the corresponding two-dimensional shallow water equations. The resulting two- and three-dimensional algorithm has been shown to be fast, accurate and mass-conservative and can also be applied to simulate flooding and drying of tidal mud-flats in conjunction with three-dimensional flows. Furthermore, the resulting algorithm is fully vectorizable for an efficient implementation on modern vector computers.

International Journal for Numerical Methods in Flu↗
USGS2011

Adaptive finite volume methods with well-balanced Riemann solvers for modeling floods in rugged terrain: Application to the Malpasset dam-break flood (France, 1959)

The simulation of advancing flood waves over rugged topography, by solving the shallow-water equations with well-balanced high-resolution finite volume methods and block-structured dynamic adaptive mesh refinement (AMR), is described and validated in this paper. The efficiency of block-structured AMR makes large-scale problems tractable, and allows the use of accurate and stable methods developed for solving general hyperbolic problems on quadrilateral grids. Features indicative of flooding in rugged terrain, such as advancing wet-dry fronts and non-stationary steady states due to balanced source terms from variable topography, present unique challenges and require modifications such as special Riemann solvers. A well-balanced Riemann solver for inundation and general (non-stationary) flow over topography is tested in this context. The difficulties of modeling floods in rugged terrain, and the rationale for and efficacy of using AMR and well-balanced methods, are presented. The algorithms are validated by simulating the Malpasset dam-break flood (France, 1959), which has served as a benchmark problem previously. Historical field data, laboratory model data and other numerical simulation results (computed on static fitted meshes) are shown for comparison. The methods are implemented in GEOCLAW, a subset of the open-source CLAWPACK software. All the software is freely available at. Published in 2010 by John Wiley & Sons, Ltd.

International Journal for Numerical Methods in Flu↗
USGS2006

On the interpretation of energy and energy fluxes of nonlinear internal waves: An example from Massachusetts Bay

A self-consistent formalism to estimate baroclinic energy densities and fluxes resulting from the propagation of internal waves of arbitrary amplitude is derived using the concept of available potential energy. The method can be applied to numerical, laboratory or field data. The total energy flux is shown to be the sum of the linear energy flux ??? u??? p??? dz (primes denote baroclinic quantities), plus contributions from the non-hydrostatic pressure anomaly and the self-advection of kinetic and available potential energy. Using highly resolved observations in Massachusetts Bay, it is shown that due to the presence of nonlinear internal waves periodically propagating in the area, ??? u??? p??? dz accounts for only half of the total flux. The same data show that equipartition of available potential and kinetic energy can be violated, especially when the nonlinear waves begin to interact with the bottom. ?? 2006 Cambridge University Press.

Massachusetts Bay↗
USGS1981

A transient laboratory method for determining the hydraulic properties of 'tight' rocks-I. Theory

Transient pulse testing has been employed increasingly in the laboratory to measure the hydraulic properties of rock samples with low permeability. Several investigators have proposed a mathematical model in terms of an initial-boundary value problem to describe fluid flow in a transient pulse test. However, the solution of this problem has not been available. In analyzing data from the transient pulse test, previous investigators have either employed analytical solutions that are derived with the use of additional, restrictive assumptions, or have resorted to numerical methods. In Part I of this paper, a general, analytical solution for the transient pulse test is presented. This solution is graphically illustrated by plots of dimensionless variables for several cases of interest. The solution is shown to contain, as limiting cases, the more restrictive analytical solutions that the previous investigators have derived. A method of computing both the permeability and specific storage of the test sample from experimental data will be presented in Part II.

International Journal of Rock Mechanics and Mining↗
USGS2001

Flow of variably fluidized granular masses across three-dimensional terrain I. Coulomb mixture theory

Rock avalanches, debris flows, and related phenomena consist of grain-fluid mixtures that move across three-dimensional terrain. In all these phenomena the same basic forces, govern motion, but differing mixture compositions, initial conditions, and boundary conditions yield varied dynamics and deposits. To predict motion of diverse grain-fluid masses from initiation to deposition, we develop a depth-averaged, threedimensional mathematical model that accounts explicitly for solid- and fluid-phase forces and interactions. Model input consists of initial conditions, path topography, basal and internal friction angles of solid grains, viscosity of pore fluid, mixture density, and a mixture diffusivity that controls pore pressure dissipation. Because these properties are constrained by independent measurements, the model requires little or no calibration and yields readily testable predictions. In the limit of vanishing Coulomb friction due to persistent high fluid pressure the model equations describe motion of viscous floods, and in the limit of vanishing fluid stress they describe one-phase granular avalanches. Analysis of intermediate phenomena such as debris flows and pyroclastic flows requires use of the full mixture equations, which can simulate interaction of high-friction surge fronts with more-fluid debris that follows. Special numerical methods (described in the companion paper) are necessary to solve the full equations, but exact analytical solutions of simplified equations provide critical insight. An analytical solution for translational motion of a Coulomb mixture accelerating from rest and descending a uniform slope demonstrates that steady flow can occur only asymptotically. A solution for the asymptotic limit of steady flow in a rectangular channel explains why shear may be concentrated in narrow marginal bands that border a plug of translating debris. Solutions for static equilibrium of source areas describe conditions of incipient slope instability, and other static solutions show that nonuniform distributions of pore fluid pressure produce bluntly tapered vertical profiles at the margins of deposits. Simplified equations and solutions may apply in additional situations identified by a scaling analysis. Assessment of dimensionless scaling parameters also reveals that miniature laboratory experiments poorly simulate the dynamics of full-scale flows in which fluid effects are significant. Therefore large geophysical flows can exhibit dynamics not evident at laboratory scales.

Journal of Geophysical Research B: Solid Earth↗
USGS2007

Spectral element modelling of fault-plane reflections arising from fluid pressure distributions

The presence of fault-plane reflections in seismic images, besides indicating the locations of faults, offers a possible source of information on the properties of these poorly understood zones. To better understand the physical mechanism giving rise to fault-plane reflections in compacting sedimentary basins, we numerically model the full elastic wavefield via the spectral element method (SEM) for several different fault models. Using well log data from the South Eugene Island field, offshore Louisiana, we derive empirical relationships between the elastic parameters (e.g. P-wave velocity and density) and the effective-stress along both normal compaction and unloading paths. These empirical relationships guide the numerical modelling and allow the investigation of how differences in fluid pressure modify the elastic wavefield. We choose to simulate the elastic wave equation via SEM since irregular model geometries can be accommodated and slip boundary conditions at an interface, such as a fault or fracture, are implemented naturally. The method we employ for including a slip interface retains the desirable qualities of SEM in that it is explicit in time and, therefore, does not require the inversion of a large matrix. We performa complete numerical study by forward modelling seismic shot gathers over a faulted earth model using SEM followed by seismic processing of the simulated data. With this procedure, we construct post-stack time-migrated images of the kind that are routinely interpreted in the seismic exploration industry. We dip filter the seismic images to highlight the fault-plane reflections prior to making amplitude maps along the fault plane. With these amplitude maps, we compare the reflectivity from the different fault models to diagnose which physical mechanism contributes most to observed fault reflectivity. To lend physical meaning to the properties of a locally weak fault zone characterized as a slip interface, we propose an equivalent-layer model under the assumption of weak scattering. This allows us to use the empirical relationships between density, velocity and effective stress from the South Eugene Island field to relate a slip interface to an amount of excess pore-pressure in a fault zone. ?? 2007 The Authors Journal compilation ?? 2007 RAS.

Geophysical Journal International↗
USGS2003

A numerical simulation of magma motion, crustal deformation, and seismic radiation associated with volcanic eruptions

The finite difference method is used to calculate the magma dynamics, seismic radiation, and crustal deformation associated with a volcanic eruption. The model geometry consists of a cylindrical reservoir and narrow cylindrical conduit embedded in a homogeneous crust. We consider two models of eruption. In the first model, a lid caps the vent and the magma is overpressurized prior to the eruption. The eruption is triggered by the instantaneous removal of the lid, at which point the exit pressure becomes equal to the atmospheric pressure. In the second model, a plug at the reservoir outlet allows pressurization of only the magmatic fluid in the reservoir before the eruption. Magma transfer between the reservoir and conduit is triggered by the instantaneous removal of the plug, and the eruption occurs when the pressure at the conduit orifice exceeds the material strength of the lid capping the vent. In both models, magma dynamics are expressed by the equations of mass and momentum conservation in a compressible fluid, in which fluid expansion associated with depressurization is accounted for by a constitutive law relating pressure and density. Crustal motions are calculated from the equations of elastodynamics. The fluid and solid are dynamically coupled by applying the continuity of wall velocities and normal stresses across the conduit and reservoir boundaries. Free slip is allowed at the fluid-solid boundary. Both models predict the gradual depletion of the magma reservoir, which causes crustal deformation observed as a long-duration dilatational signal. Superimposed on this very-long-period (VLP) signal generated by mass transport are long-period (LP) oscillations of the magma reservoir and conduit excited by the acoustic resonance of the reservoir-conduit system during the eruption. The volume of the reservoir, vent size, and magma properties control the duration of VLP waves and dominant periods of LP oscillations. The second model predicts that when the magmatic fluid reaches the vent, a high-pressure pulse occurs at this location in accordance with the basic theory of compressible fluid dynamics. This abrupt pressure increase just beneath the vent is consistent with observed seismograms in which pulse-like Rayleigh waves excited by a shallow source are dominant. The strength of the lid plays an important role in the character of the seismograms and in defining the type of eruption observed.

Geophysical Journal International↗