Dispersivity Determination Through a Modeling Approach From a Tracer Test Based on Total Br Concentration in Soil Samples

Determination of reliable solute transport parameters is an essential aspect for the characterization of the mechanisms and processes involved in solute transport (e.g., pesticides, fertilizers, contaminants) through the unsaturated zone. A rapid inexpensive method to estimate the dispersivity parameter at the field scale is presented herein. It is based on the quantification by the X-ray fluorescence solid-state technique of total bromine in soil, along with an inverse numerical modeling approach. The results show that this methodology is a good alternative to the classic Br− determination in soil water by ion chromatography. A good agreement between the observed and simulated total soil Br is reported. The results highlight the potential applicability of both combined techniques to infer readily solute transport parameters under field conditions.

3D numerical modelling of acoustic horns using the method of fundamental solutions

In the present work, a three-dimensional (3D) formulation based on the method of fundamental solutions (MFS) is applied to the study of acoustic horns. The implemented model follows and extends previous works that only considered two-dimensional and axisymmetric horn configurations. The more realistic case of 3D acoustic horns with symmetry regarding two orthogonal planes is addressed. The use of the domain decomposition technique with two interconnected sub-regions along a continuity boundary is proposed, allowing for the computation of the sound pressure generated by an acoustic horn installed on a rigid screen. In order to reduce the model discretization requirements for these cases, Green’s functions derived with the image source methodology are adopted, automatically accounting for the presence of symmetry conditions. A strategy for the calculation of an optimal position of the virtual sources used by the MFS to define the solution is also used, leading to improved reliability and flexibility of the proposed method. The responses obtained by the developed model are compared to reference solutions, computed by well-established models based on the boundary element method. Additionally, numerically calculated acoustic parameters, such as directivity and beamwidth, are compared with those evaluated experimentally.

Fluid-Fluid Reactions: Study of the Surface Renewal Theory in Chemical Engineering Courses

The Surface Renewal Theory (SRT) is one of the most unfamiliar models in order to characterize fluid-fluid and fluid-fluid-solid reactions, which are of considerable industrial and academicals importance. In the present work, an approach to the resolution of the SRT model by numerical methods is presented, enabling the visualization of the influence of different variables which control the heterogeneous overall process. Its use in a classroom allowed the students to reach a great understanding of the process.

Implementation of non-local boundary layer schemes in the Regional Atmospheric Modeling System and its impact on simulated mesoscale circulations

This paper proposes the implementation of different non-local Planetary Boundary Layer schemes within the Regional Atmospheric Modeling System (RAMS) model. The two selected PBL parameterizations are the Medium-Range Forecast (MRF) PBL and its updated version, known as the Yonsei University (YSU) PBL. YSU is a first-order scheme that uses non-local eddy diffusivity coefficients to compute turbulent fluxes. It is based on the MRF, and improves it with an explicit treatment of the entrainment. With the aim of evaluating the RAMS results for these PBL parameterizations, a series of numerical simulations have been performed and contrasted with the results obtained using the Mellor and Yamada (MY) scheme, also widely used, and the standard PBL scheme in the RAMS model. The numerical study carried out here is focused on mesoscale circulation events during the summer, as these meteorological situations dominate this season of the year in the Western Mediterranean coast. In addition, the sensitivity of these PBL parameterizations to the initial soil moisture content is also evaluated. The results show a warmer and moister PBL for the YSU scheme compared to both MRF and MY. The model presents as well a tendency to overestimate the observed temperature and to underestimate the observed humidity, considering all PBL schemes and a low initial soil moisture content. In addition, the bias between the model and the observations is significantly reduced moistening the initial soil moisture of the corresponding run. Thus, varying this parameter has a positive effect and improves the simulated results in relation to the observations. However, there is still a significant overestimation of the wind speed over flatter terrain, independently of the PBL scheme and the initial soil moisture used, even though a different degree of accuracy is reproduced by RAMS taking into account the different sensitivity tests.

Thermodynamic, surface and structural properties of liquid Co-Si alloys

The knowledge of thermophysical properties of liquid Co-Si alloys is a key requirement for manufacturing of composite materials by infiltration method. Despite this need, the experimental and predicted property data of the Co-Si system are scarce and often inconsistent between the various sources. In the present work the mixing behaviour of Co-Si melts has been analysed through the study of the concentration dependence of various thermodynamic, surface (surface tension and surface composition) and structural properties (concentration fluctuations in the long-wavelength limit and chemical short-range order parameter) in the framework of the Compound Formation Model (CFM) and Quasi Chemical Approximation for regular solutions (QCA). In addition, the surface tension of the Co22·5Si77.5 (in at%) eutectic alloy, that is proposed to be used as the infiltrant, has been measured by the pendant drop method at temperatures ranging from 1593 to 1773 K. The results obtained were discussed with respect to both, temperature and concentration, and subsequently compared with the model predictions and literature data.

A Finite Element Numerical Algorithm for Modelling and Data Fitting in Complex Systems

Numerical modelling methodologies are important by their application to engineering and scientific problems, because there are processes where analytical mathematical expressions cannot be obtained to model them. When the only available information is a set of experimental values for the variables that determine the state of the system, the modelling problem is equivalent to determining the hyper-surface that best fits the data. This paper presents a methodology based on the Galerkin formulation of the finite elements method to obtain representations of relationships that are defined a priori, between a set of variables: y = z(x1, x2,...., xd). These representations are generated from the values of the variables in the experimental data. The approximation, piecewise, is an element of a Sobolev space and has derivatives defined in a general sense into this space. The using of this approach results in the need of inverting a linear system with a structure that allows a fast solver algorithm. The algorithm can be used in a variety of fields, being a multidisciplinary tool. The validity of the methodology is studied considering two real applications: a problem in hydrodynamics and a problem of engineering related to fluids, heat and transport in an energy generation plant. Also a test of the predictive capacity of the methodology is performed using a cross-validation method.

ROCK OUTCROPS IN THE CANADIAN SHIELD: AN INVESTIGATION OF CONTAMINANT TRANSPORT FROM SURFACE SOURCES IN FRACTURED ROCK AQUIFERS

Flow, recharge and transport dynamics in fractured rock aquifers with low lying rock outcrops is a largely unexplored area of study in hydrogeology. The purpose of this thesis is to examine these topics in an agricultural area in Eastern Ontario. The study consists of a regional scale groundwater quality study, an infiltration experiment that considers bacteria transport from the ground surface to a well, and a numerical modelling study that tests the parameters that affect surface infiltration of a tracer from a rock outcrop to a deeper horizontal fracture. In the water quality study, approximately 65% of the samples contained total coliform, 16% contained E. coli, and 1% contained nitrate-N at greater than 5 mg/L. Occurrence of E. coli increased when considering seasonality, where wells were drilled on rock outcrops, and for shallow well intervals. Nitrate-N did not occur above the Guidelines for Canadian Drinking Water Quality (Health Canada, 2012) of 10 mg/L. Rapid arrival times were observed in the infiltration study for both the microspheres (30 minutes) and a dye tracer (45 minutes) in a well approximately 6.0 m in horizontal and 2.8 m in vertical distance from the tracer source. Transport velocities were approximately 38.9 m/day for the dye tracer and 115.2 m/day for the colloidal tracer. Results of the model runs indicate that overburden can provide an effective protective layer to transport in fractures, that high groundwater velocities occur in larger fracture apertures and higher gradients dilute tracer concentrations, and that lower groundwater velocities occur with smaller fracture apertures and lower gradients result in elevated tracer concentrations. Lower rainfall rates, larger fracture apertures, early tracer time, larger gradients, and lower water levels maintained unsaturated conditions for longer time periods such that tracer transport was delayed until saturated conditions were attained. The overall heterogeneity of this aquifer environment creates a source water protection conundrum where the water quality is generally good, while transport can occur very quickly in proximity to rock outcrops and in areas with limited overburden.

Overcoming performance and convergence issues of discrete transform based modeling of crosslinking classical and reversible deactivated radical polymerizations

Population balances of polymer species in terms 'of discrete transforms with respect to counts of groups lead to tractable first order partial differential equations when ali rate constants are independent of chain length and loop formation is negligible [l]. Average molecular weights in the absence ofgelation are long known to be readily found through integration of an initial value problem. The extension to size distribution prediction is also feasible, but its performance is often lower to the one provided by methods based upon real chain length domain [2]. Moreover, the absence ofagood starting procedure and a higher numerical sensitivity hás decisively impaired its application to non-linear reversibly deactivated polymerizations, namely NMRP [3].

Numerical and experimental investigation of structural stiffness influence on ratcheting convection cell in granular soils under cyclic loading

Lateral cyclic loaded structures in granular soils can lead to an accumulation of irreversible strains by changing their mechanical response (densification) and forming a closed convective cell in the upper layer of the bedding. In the present thesis the convective cell dimension, formation and grain migration inside this closed volume have been studied and presented in relation to structural stiffness and different loads. This relation was experimentally investigated by applying a cyclic lateral force to a scaled flexible vertical element embedded in dry granular soil. The model was monitored with a camera in order to derive the displacement field by means of the PIV technique. Modelling large soil deformation turns out to be difficult, using mesh-based methods. Consequently, a mesh-free approach (DEM) was chosen in order to investigate the granular flow with the aim of extracting interesting micromechanical information. In both the numerical and experimental analyses the effect of different loading magnitudes and different dimensions of the vertical element were considered. The main results regarded the different development, shape and dimensions of the convection cell and the surface settlements. Moreover, the Discrete Element Method has proven to give satisfactory results in the modelling of large deformation phenomena such as the ratcheting convective cell.

Benchmarking numerical models of brittle thrust wedges

We report quantitative results from three brittle thrust wedge experiments, comparing numerical results directly with each other and with corresponding analogue results. We first test whether the participating codes reproduce predictions from analytical critical taper theory. Eleven codes pass the stable wedge test, showing negligible internal deformation and maintaining the initial surface slope upon horizontal translation over a frictional interface. Eight codes participated in the unstable wedge test that examines the evolution of a wedge by thrust formation from a subcritical state to the critical taper geometry. The critical taper is recovered, but the models show two deformation modes characterised by either mainly forward dipping thrusts or a series of thrust pop-ups. We speculate that the two modes are caused by differences in effective basal boundary friction related to different algorithms for modelling boundary friction. The third experiment examines stacking of forward thrusts that are translated upward along a backward thrust. The results of the seven codes that run this experiment show variability in deformation style, number of thrusts, thrust dip angles and surface slope. Overall, our experiments show that numerical models run with different numerical techniques can successfully simulate laboratory brittle thrust wedge models at the cm-scale. In more detail, however, we find that it is challenging to reproduce sandbox-type setups numerically, because of frictional boundary conditions and velocity discontinuities. We recommend that future numerical-analogue comparisons use simple boundary conditions and that the numerical Earth Science community defines a plasticity test to resolve the variability in model shear zones.

Benchmarking analogue models of brittle thrust wedges

We performed a quantitative comparison of brittle thrust wedge experiments to evaluate the variability among analogue models and to appraise the reproducibility and limits of model interpretation. Fifteen analogue modeling laboratories participated in this benchmark initiative. Each laboratory received a shipment of the same type of quartz and corundum sand and all laboratories adhered to a stringent model building protocol and used the same type of foil to cover base and sidewalls of the sandbox. Sieve structure, sifting height, filling rate, and details on off-scraping of excess sand followed prescribed procedures. Our analogue benchmark shows that even for simple plane-strain experiments with prescribed stringent model construction techniques, quantitative model results show variability, most notably for surface slope, thrust spacing and number of forward and backthrusts. One of the sources of the variability in model results is related to slight variations in how sand is deposited in the sandbox. Small changes in sifting height, sifting rate, and scraping will result in slightly heterogeneous material bulk densities, which will affect the mechanical properties of the sand, and will result in lateral and vertical differences in peak and boundary friction angles, as well as cohesion values once the model is constructed. Initial variations in basal friction are inferred to play the most important role in causing model variability. Our comparison shows that the human factor plays a decisive role, and even when one modeler repeats the same experiment, quantitative model results still show variability. Our observations highlight the limits of up-scaling quantitative analogue model results to nature or for making comparisons with numerical models. The frictional behavior of sand is highly sensitive to small variations in material state or experimental set-up, and hence, it will remain difficult to scale quantitative results such as number of thrusts, thrust spacing, and pop-up width from model to nature.