Numerical Homogenization of the Rigidity Tensor in Hooke's Law Using the Node-Based Finite Element Method

Combining a geological model with a geomechanical model, it generally turns out that the geomechanical model is built from units that are at least a 100 times larger in volume than the units of the geological model. To counter this mismatch in scales, the geological data model's heterogeneous fine-scale Young's moduli and Poisson's ratios have to be “upscaled” to one “equivalent homogeneous” coarse-scale rigidity. This coarse-scale rigidity relates the volume-averaged displacement, strain, stress, and energy to each other, in such a way that the equilibrium equation, Hooke's law, and the energy equation preserve their fine-scale form on the coarse scale. Under the simplifying assumption of spatial periodicity of the heterogeneous fine-scale rigidity, homogenization theory can be applied. However, even then the spatial variability is generally so complex that exact solutions cannot be found. Therefore, numerical approximation methods have to be applied. Here the node-based finite element method for the displacement as primary variable has been used. Three numerical examples showing the upper bound character of this finite element method are presented.     

The identifiability of parameters in a water quality model of the Biebrza River, Poland

The identifiability of model parameters of a steady state water quality model of the Biebrza River and the resulting variation in model results was examined by applying the Monte Carlo method which combines calibration, identifiability analysis, uncertainty analysis, and sensitivity analysis. The water quality model simulates the steady state concentration profiles of chloride, phosphate, ammonium, and nitrate as a function of distance along a river. The water quality model with the best combination of parameter values simulates the observed concentrations very well. However, the range of possible modelled concentrations obtained for other more or less equally eligible combinations of parameter values is rather wide. This range in model outcomes reflects possible errors in the model parameters. Discrepancies between the range in model outcomes and the validation data set are only caused by errors in model structure, or (measurement) errors in boundary conditions or input variables. In this sense the validation procedure is a test of model capability, where the effects of calibration errors are filtered out. It is concluded that, despite some slight deviations between model outcome and observations, the model is successful in simulating the spatial pattern of nutrient concentrations in the Biebrza River.     

Turbulent Pressure Statistics in the Atmospheric Boundary Layer from Large-Eddy Simulation

We use large-eddy simulation (LES) to study the turbulent pressure field in atmospheric boundary layers with free convection, forced convection, and stable stratification. We use the Poisson equation for pressure to represent the pressure field as the sum of mean-shear, turbulence–turbulence, subfilter-scale, Coriolis, and buoyancy contributions. We isolate these contributions and study them separately. We find that in the energy-containing range in the free-convection case the turbulence–turbulence pressure dominates over the entire boundary layer. That part dominates also up to midlayer in the forced-convection case; above that the mean-shear pressure dominates. In the stable case the mean-shear pressure dominates over the entire boundary layer.We find evidence of an inertial subrange in the pressure spectrum in the free and forced-convection cases; it is dominated by the turbulence–turbulence pressure and has a three-dimensional spectral constant of about 4.0. This agrees well with quasi-Gaussian predictions but is a factor of 2 less than recent results from direct numerical simulations at moderate Reynolds numbers. Measurements of the inertial subrange pressure spectral constant at high Reynolds numbers, which might now be possible, would be most useful.     

Preliminary Characterisation of New Glass Reference Materials (GSA-1G, GSC-1G, GSD-1G and GSE-1G) by Laser Ablation-Inductively Coupled Plasma-Mass Spectrometry Using 193 nm, 213 nm and 266 nm Wavelengths

New glass reference materials GSA-1G, GSC-1G, GSD-1G and GSE-1G have been characterised using a prototype solid state laser ablation system capable of producing wavelengths of 193 nm, 213 nm and 266 nm. This system allowed comparison of the effects of different laser wavelengths under nearly identical ablation and ICP operating conditions. The wavelengths 213 nm and 266 nm were also used at higher energy densities to evaluate the influence of energy density on quantitative analysis. In addition, the glass reference materials were analysed using commercially available 266 nm Nd:YAG and 193 nm ArF excimer lasers. Laser ablation analysis was carried out using both single spot and scanning mode ablation. Using laser ablation ICP-MS, concentrations of fifty-eight elements were determined with external calibration to the NIST SRM 610 glass reference material. Instead of applying the more common internal standardisation procedure, the total concentration of all element oxide concentrations was normalised to 100%. Major element concentrations were compared with those determined by electron microprobe. In addition to NIST SRM 610 for external calibration, USGS BCR-2G was used as a more closely matrix-matched reference material in order to compare the effect of matrix-matched and non matrix-matched calibration on quantitative analysis. The results show that the various laser wavelengths and energy densities applied produced similar results, with the exception of scanning mode ablation at 266 nm without matrix-matched calibration where deviations up to 60% from the average were found. However, results acquired using a scanning mode with a matrix-matched calibration agreed with results obtained by spot analysis. The increased abundance of large particles produced when using a scanning ablation mode with NIST SRM 610, is responsible for elemental fractionation effects caused by incomplete vaporisation of large particles in the ICP.     

Development and Extension of An Aggregated Scale Model: Part 1 – Background to ASMITA

Whilst much attention has been given to models that describe wave, tide and sediment transport processes in sufficient detail to determine the local changes in bed level over a relatively detailed representation of the bathymetry, far less attention has been given to models that consider the problem at a much larger scale (e.g. that of geomorphological elements such as a tidal flat and tidal channel). Such aggregated or lumped models tend not to represent the processes in detail but rather capture the behaviour at the scale of interest. One such model developed using the concept of an equilibrium concentration is the Aggregated Scale Morphological Interaction between Tidal basin and Adjacent coast (ASMITA). In this paper we provide some new insights into the concepts of equilibrium, and horizontal and vertical exchange that are key components of this modelling approach. In a companion paper, we summarise a range of developments that have been undertaken to extend the original model concept, to illustrate the flexibility and power of the conceptual framework. However, adding detail progressively moves the model in the direction of the more detailed process-based models and we give some consideration to the boundary between the two.     

An ecosystem-based composite spatial model for floodplain vulnerability assessment: a case study of Artigas,Uruguay

Floods are natural processes that constitute a hazard to society when associated to improper land use. Anthropic activities in floodplains are a factor of vulnerability that converts a natural hazard into a threat factor, eventually leading to disaster. Nowadays, natural and social complex processes demand integrated assessments in order to improve their understanding, helping decision making over sustainable use of territory, as well as integrating society’s activity in ecosystems and potentials, restrictions and benefits that society obtain from them. In this context, the objective of this work was to build a composite vulnerability model for a floodplain under urban influence, using an integrated assessment approach. This model was based on three dimensions; threat, fragility and an ecosystem services provision. These dimensions were calculated using both primary and secondary information, and weights by specialists. Main results show that the area presents high vulnerability with an increasing gradient towards high and urbanized areas, associated with an important number and relevant ecosystem services. Also, a spatial heterogeneity of the three dimensions emerged, making evident this area’s complexity and the need of integrated assessments to approach it. The composite vulnerability model proposed presents an elevated potential for natural and social processes analysis in floodplains, which is crucial for these territory management. Moreover, these integrated dimensions could contribute to decision making in different levels, as well as generating important supplies for environmental management and land planning.

     

Effective fetch and non-linear four-wave interactions during wave growth in slanting fetch conditions

Wave growth in slanting fetch (with wind blowing obliquely off a coast) is investigated with 7 years worth of routine wave measurements in Lake IJssel in The Netherlands and with the SWAN wave model. Two aspects are considered in particular for this case: the validity of the concept of effective fetch and the role of the non-linear four-wave interactions. For slanting and parallel fetch conditions, we found some significant deviations from the effective fetch assumption, leading to 20–35% mismatch in either the peak period T_p or the significant wave height H_m0 respectively. However, the effect of discrepancies between various widely accepted wave growth formulas turned out to be even more important. The wave directions during slanting fetch are significantly ‘steered’ by the coastline, especially in the first kilometre(s) off the coast. The role of the non-linear four-wave interactions is investigated by running the SWAN (version 40.41) wave model with three different quadruplet formulations. Exact quadruplet methods (Xnl) yielded relatively strong wave steering, despite the four-wave interactions being relatively weak. Application of Xnl did not lead to better overall agreement with measurements — improvements for the mean wave period T_m01 were offset by some deterioration for the wave height H_m0.     

Uncertainty in the application of the parabolic bay shape equation: Part 2

This paper attempts to highlight issues that are relevant in the application of the Parabolic Bay Shape Equation (PBSE) to a non-equilibrium bay. For this case, the bay of Imbituba in southern Brazil was chosen. The construction of a breakwater to shelter the port of Imbituba in the south of the bay was accompanied by an increase in sedimentation to the port area from an eroding downdrift beach. Superimposed plots of the coastline of the Bay of Imbituba from different years confirm a general trend of accretion of the southern part of the bay accompanied with a retreat of the coastline in the northern part. After the application of the PBSE it became clear that the breakwater caused a change in the equilibrium state of the bay. Between 1947 and 2001 the Bay of Imbituba has changed from a dynamic equilibrium to a natural beach reshaping or self-reshaping. The tendency of the sedimentation of the southern part of the bay can be explained by the SEP associated with the new updrift diffraction point (tip of the breakwater): The seaward position of the SEP predicts a need for additional sediment in order to achieve a stable plan form.