A global method for coupling transport with chemistry in heterogeneous porous media

Modeling reactive transport in porous media, using a local chemical equilibrium assumption, leads to a system of advection–diffusion PDEs coupled with algebraic equations. When solving this coupled system, the algebraic equations have to be solved at each grid point for each chemical species and at each time step. This leads to a coupled non-linear system. In this paper, a global solution approach that enables to keep the software codes for transport and chemistry distinct is proposed. The method applies the Newton–Krylov framework to the formulation for reactive transport used in operator splitting. The method is formulated in terms of total mobile and total fixed concentrations and uses the chemical solver as a black box, as it only requires that one be able to solve chemical equilibrium problems (and compute derivatives) without having to know the solution method. An additional advantage of the Newton–Krylov method is that the Jacobian is only needed as an operator in a Jacobian matrix times vector product. The proposed method is tested on the MoMaS reactive transport benchmark.     

Experimental and Numerical Studies of Velocity and Turbulence Intensities for Mid-Channel Bar

Water Resources - Turbulent phenomenon in braided rivers is much more complex compared to the straight and meandering rivers. The turbulent flow structure at locations upstream and downstream of...     

Performance comparison between logistic and generalized surplus-production models applied to the Sillago sihama fishery in Pakistan

The catch and effort data of Sillago sihama fishery in Pakistani waters were used to investigate the performance of two closely related stock assessment models: logistic and generalized surplus-production models. Compared with the generalized production model, the logistic model produced more reasonable estimates for parameters such as maximum sustainable yield. The Akaike’s Information Criterion values estimated at 4.265 and -51.152 respectively by the logistic and generalized models. Simulation analyses of the S. sihama fishery showed that the estimated and observed abundance indices for the logistic model were closer than those for the generalized production model. Standardized residuals were distributed closer for logistic model, but exhibited a slightly increasing trend for the generalized model. Statistical outliers were seen in 1989 and 1993 for the logistic model, and in 1981 and 1999 for the generalized model. Simulated results revealed that the logistic estimates were close to the true value for low CVs (coefficients of variation) but widely dispersed for high CVs. In contrast, the generalized model estimates were loose for all CV levels. The estimated production model curve parameter was not reasonable at all the tested levels of white noise. With the increase in white noise R2 for the catch per unit effort decreased. Therefore, we conclude that the logistic model performs more reasonably than the generalized production model.     

Mapping dune dynamics by InSAR coherence

It is generally held that subtle changes in sandy environments are very difficult to detect in imagery. Nonetheless, this study demonstrates how synthetic aperture radar (SAR) interferometric decorrelation can be used to identify changes in individual sand dunes. The use of coherence maps over time facilitates the analysis of dune dynamics, both temporally and spatially. The Ashdod‐Nizzanim coastal dunes, along the southern coastal plain of Israel, were chosen as an illustrative example of the analysis of dune dynamics. High‐resolution TerraSAR‐X (TSX) radar images covering the entire research area were acquired for the period February to July 2012, together with meteorology data (wind and rain) for the area. The coherence results enabled the stability of individual dunes to be described as a function of time. It was found that the dune crests were more stable than the windward slopes and that the degree of stability was dependent on the distance of the dune from the sea. The results of this study show the potential of using interferometric synthetic aperture radar (InSAR) decorrelation for aeolian studies, even in areas characterized by low coherence. Copyright © 2017 John Wiley & Sons, Ltd.     

Application of Hamilton's law of varying action to the restricted three-body problem

The Law of Varying Action, originally published by Hamilton in 1834, was recently employed by Bailey in devising a technique for generating power series characterizing the motions of dynamical systems. Furthermore, Bailey's method permits one to construct these series in a simple and direct way, without using the associated differential equations of motion. In the present paper, Hamilton's law is developed in its most general form and is used to produce series solutions of the restricted three-body problem. Finally, for illustrative purposes, numerical results are presented for several symmetric periodic orbits.The authors are indebted to Professor Thomas R. Kane of Stanford University for this idea.     

Deriving snow‐cover depletion curves for different spatial scales from remote sensing and snow telemetry data

During the melting of a snowpack, snow water equivalent (SWE) can be correlated to snow‐covered area (SCA) once snow‐free areas appear, which is when SCA begins to decrease below 100%. This amount of SWE is called the threshold SWE. Daily SWE data from snow telemetry stations were related to SCA derived from moderate‐resolution imaging spectroradiometer images to produce snow‐cover depletion curves. The snow depletion curves were created for an 80 000 km² domain across southern Wyoming and northern Colorado encompassing 54 snow telemetry stations. Eight yearly snow depletion curves were compared, and it is shown that the slope of each is a function of the amount of snow received. Snow‐cover depletion curves were also derived for all the individual stations, for which the threshold SWE could be estimated from peak SWE and the topography around each station. A station's peak SWE was much more important than the main topographic variables that included location, elevation, slope, and modelled clear sky solar radiation. The threshold SWE mostly illustrated inter‐annual consistency. Copyright © 2015 John Wiley & Sons, Ltd.     

Nodal-based three-dimensional discontinuous deformation analysis (3-D DDA)

This paper presents a new numerical model that can add a finite element mesh into each block of the three-dimensional discontinuous deformation analysis (3-D DDA), originally developed by Gen-hua Shi. The main objectives of this research are to enhance DDA block’s deformability. Formulations of stiffness and force matrices in 3-D DDA with conventional Trilinear (8-node) and Serendipity (20-node) hexahedral isoparametric finite elements meshed block system due to elastic stress, initial stress, point load, body force, displacement constraints, inertia force, normal and shear contact forces are derived in detail for program coding. The program code for the Trilinear and Serendipity hexahedron elements have been developed, and it has been applied to some examples to show the advantages achieved when finite element is associated with 3-D DDA to handle problems under large displacements and deformations. Results calculated for the same models by use of the original 3-D DDA are far from the theoretical solutions while the results of new numerical model are quite good in agreement with theoretical solutions; however, for the Trilinear elements, more number of elements are needed.     

Cartography-Oriented Design of 3D Geospatial Information Visualization – Overview and Techniques

In economy, society and personal life map-based interactive geospatial visualization becomes a natural element of a growing number of applications and systems. The visualization of 3D geospatial information, however, raises the question how to represent the information in an effective way. Considerable research has been done in technology-driven directions in the fields of cartography and computer graphics (e.g., design principles, visualization techniques). Here, non-photorealistic rendering (NPR) represents a promising visualization category – situated between both fields – that offers a large number of degrees for the cartography-oriented visual design of complex 2D and 3D geospatial information for a given application context. Still today, however, specifications and techniques for mapping cartographic design principles to the state-of-the-art rendering pipeline of 3D computer graphics remain to be explored. This paper revisits cartographic design principles for 3D geospatial visualization and introduces an extended 3D semiotic model that complies with the general, interactive visualization pipeline. Based on this model, we propose NPR techniques to interactively synthesize cartographic renditions of basic feature types, such as terrain, water, and buildings. In particular, it includes a novel iconification concept to seamlessly interpolate between photorealistic and cartographic representations of 3D landmarks. Our work concludes with a discussion of open challenges in this field of research, including topics, such as user interaction and evaluation.     

Stability of rubble-mound breakwater using H50 wave height parameter

The prediction of rubble mound breakwaters' stability is one of the most important issues in coastal and maritime engineering. The stability of breakwaters strongly depends on the wave height. Therefore, selection of an appropriate wave height parameter is very vital in the prediction of stability number. In this study, H₅₀, the average of the 50 highest waves that reach the breakwater in its useful life, was used to predict the stability of the armor layer. First, H₅₀ was used instead of the significant wave height in the most recent stability formulas. It was found that this modification yields more accurate results. Then, for further improvement of the results, two formulas were developed using model tree.To develop the new formulas, two experimental data sets of irregular waves were used. Results indicated that the proposed formulas are more accurate than the previous ones for the prediction of the stability parameter. Finally, the proposed formulas were applied to regular waves and a wide range of damage levels and it was seen that the developed formulas are applicable in these cases as well.