Estimation of dispersion in estuaries by a self-consistent method

Longitudinal and transverse dispersion in estuaries may be determined in principle by measuring the change in time of the concentration field of a water soluble dye. In practice, measurements may take the form of several sampling runs crossing the dye patch. The coefficients of dispersion are usually estimated from the second moments of the concentration field; the use of second moments, in particular, introduces errors because measurements at long times or distance are unduly emphasized.A self-consistent method, free of the disadvantages inherent in the use of higher moments, is described in which the distribution of concentration in space is related to the zeroth spatial moment of concentration. In this formulation, based upon sampling runs across the dye patch, the concentration field in the estuary may be directly and simply related to the dispersion coefficients in both the intercept and gradient of a linear graph.The methods are illustrated by the characterization of a series of experiments in the Inner Estuary of Milford Haven.     

Numerical calculations in the orbital determination of an artificial satellite for a long arc

Numerical methods are usually used for the computation of ephemerides with perturbations for the precise orbital determination of an artificial satellite. But their numerical stability will be encountered in a long arc. In this case the use the improved Encke special perturbation methods has been suggested. The results of this paper show that Encke's method does indeed have a certain effectiveness, but cannot yet completely resolve the numerical stability, and the more efficient method is to use the energy integral or its variational relation to control the growth of the along-track error in general numerical calculations so that the aim of stabilization can be achieved.     

A numerical method for differential equations of second order

A new algorithm of order five is presented for the solution of the initial value problem where the system of ordinary differential equations is of second order and does not contain the first derivative.     

Numerical modelling of the impact crater depth-diameter dependence in an acoustically fluidized target

2D numerical modelling of impact cratering has been utilized to quantify an important depth-diameter relationship for different crater morphologies, simple and complex. It is generally accepted that the final crater shape is the result of a gravity-driven collapse of the transient crater, which is formed immediately after the impact. Numerical models allow a quantification of the formation of simple craters, which are bowl-shaped depressions with a lens of rock debris inside, and complex craters, which are characterized by a structural uplift. The computation of the cratering process starts with the first contact of the impactor and the planetary surface and ends with the morphology of the final crater. Using different rheological models for the sub-crater rocks, we quantify the influence on crater mechanics. To explain the formation of complex craters in accordance to the threshold diameter between simple and complex craters, we utilize the Acoustic Fluidization model. We carried out a series of simulations over a broad parameter range with the goal to fit the observed depth/diameter relationships as well as the observed threshold diameters on the Moon, Earth and Venus.     

Three-dimensional discontinuous deformation analysis (3-D DDA) using a new contact resolution algorithm

This paper presents a new contact calculating algorithm for contacts between two polyhedra with planar boundaries in the three-dimensional discontinuous deformation analysis (3-D DDA). In this algorithm, all six type contacts in 3-D (vertex-to-face, vertex-to-edge, vertex-to-vertex, face-to-face, edge-to-edge, and edge-to-face) are simply transformed into the form of point-to-face contacts. The presented algorithm is a simple and efficient method and it can be easily coded into a computer program. In this paper, formulations of normal contact, shear contact and frictional force submatrices based on the new method are derived and the algorithm has been programmed in VC⁺⁺. Examples are provided to demonstrate the new contact rule between two blocks.     

On CO<Subscript>2</Subscript> fluid flow and heat transfer behavior in the subsurface,following leakage from a geologic storage reservoir

Geologic storage of CO₂ is expected to produce plumes of large areal extent, and some leakage may occur along fractures, fault zones, or improperly plugged pre-existing wellbores. A review of physical and chemical processes accompanying leakage suggests a potential for self-enhancement. The numerical simulations presented here confirm this expectation, but reveal self-limiting features as well. It seems unlikely that CO₂ leakage could trigger a high-energy run-away discharge, a so-called “pneumatic eruption,” but present understanding is insufficient to rule out this possibility. The most promising avenue for increasing understanding of CO₂ leakage behavior is the study of natural analogues.     

大鹏湾潮流数值计算

应用分步杂交方法建立了大鹏湾边界潮流数值模型，并对其进行了模拟计算，重现了该区域潮波及潮流的分布规律，并根据实测资料对其进行了验证，从而可系统地了解这一海区的潮波系统及潮流场等海洋要素的分布情况。     

Generation of random waves in time-dependent extended mild-slope equations using a source function method

We develop techniques of numerical wave generation in the time-dependent extended mild-slope equations of Suh et al. [1997. Time-dependent equations for wave propagation on rapidly varying topography. Coastal Engineering 32, 91–117] and Lee et al. [2003. Extended mild-slope equation for random waves. Coastal Engineering 48, 277–287] for random waves using a source function method. Numerical results for both regular and irregular waves in one and two horizontal dimensions show that the wave heights and the frequency spectra are properly reproduced. The waves that pass through the wave generation region do not cause any numerical disturbances, showing usefulness of the source function method in avoiding re-reflection problems at the offshore boundary.     

Freely floating-body simulation by a 2D fully nonlinear numerical wave tank

Nonlinear interactions between large waves and freely floating bodies are investigated by a 2D fully nonlinear numerical wave tank (NWT). The fully nonlinear 2D NWT is developed based on the potential theory, MEL/material-node time-marching approach, and boundary element method (BEM). A robust and stable 4th-order Runge–Kutta fully updated time-integration scheme is used with regriding (every time step) and smoothing (every five steps). A special φ_n-η type numerical beach on the free surface is developed to minimize wave reflection from end-wall and wave maker. The acceleration-potential formulation and direct mode-decomposition method are used for calculating the time derivative of velocity potential. The indirect mode-decomposition method is also independently developed for cross-checking. The present fully nonlinear simulations for a 2D freely floating barge are compared with the corresponding linear results, Nojiri and Murayama’s (Trans. West-Jpn. Soc. Nav. Archit. 51 (1975)) experimental results, and Tanizawa and Minami’s (Abstract for the 6th Symposium on Nonlinear and Free-surface Flow, 1998) fully nonlinear simulation results. It is shown that the fully nonlinear results converge to the corresponding linear results as incident wave heights decrease. A noticeable discrepancy between linear and fully nonlinear simulations is observed near the resonance area, where the second and third harmonic sway forces are even bigger than the first harmonic component causing highly nonlinear features in sway time series. The surprisingly large second harmonic heave forces in short waves are also successfully reproduced. The fully updated time-marching scheme is found to be much more robust than the frozen-coefficient method in fully nonlinear simulations with floating bodies. To compare the role of free-surface and body-surface nonlinearities, the body-nonlinear-only case with linearized free-surface condition was separately developed and simulated.