Insights into the Design and Characteristics of the Sedimentation Process at Multi-Barrel Culverts in Iowa, TR-596, 2010

The present study is an integral part of a broader study focused on the design and implementation of self-cleaning culverts, i.e., configurations that prevent the formation of sediment deposits after culvert construction or cleaning. Sediment deposition at culverts is influenced by many factors, including the size and characteristics of material of which the channel is composed, the hydraulic characteristics generated under different hydrology events, the culvert geometry design, channel transition design, and the vegetation around the channel. The multitude of combinations produced by this set of variables makes the investigation of practical situations a complex undertaking. In addition to the considerations above, the field and analytical observations have revealed flow complexities affecting the flow and sediment transport through culverts that further increase the dimensions of the investigation. The flow complexities investigated in this study entail: flow non-uniformity in the areas of transition to and from the culvert, flow unsteadiness due to the flood wave propagation through the channel, and the asynchronous correlation between the flow and sediment hydrographs resulting from storm events. To date, the literature contains no systematic studies on sediment transport through multi-box culverts or investigations on the adverse effects of sediment deposition at culverts. Moreover, there is limited knowledge about the non-uniform, unsteady sediment transport in channels of variable geometry. Furthermore, there are few readily useable (inexpensive and practical) numerical models that can reliably simulate flow and sediment transport in such complex situations. Given the current state of knowledge, the main goal of the present study is to investigate the above flow complexities in order to provide the needed insights for a series of ongoing culvert studies. The research was phased so that field observations were conducted first to understand the culvert behavior in Iowa landscape. Modeling through complementary hydraulic model and numerical experiments was subsequently carried out to gain the practical knowledge for the development of the self-cleaning culvert designs.

Simulation of surface waves in porous media

We present a novel numerical algorithm for the simulation of seismic wave propagation in porous media, which is particularly suitable for the accurate modelling of surface wave-type phenomena. The differential equations of motion are based on Biot's theory of poro-elasticity and solved with a pseudospectral approach using Fourier and Chebyshev methods to compute the spatial derivatives along the horizontal and vertical directions, respectively. The time solver is a splitting algorithm that accounts for the stiffness of the differential equations. Due to the Chebyshev operator the grid spacing in the vertical direction is non-uniform and characterized by a denser spatial sampling in the vicinity of interfaces, which allows for a numerically stable and accurate evaluation of higher order surface wave modes. We stretch the grid in the vertical direction to increase the minimum grid spacing and reduce the computational cost. The free-surface boundary conditions are implemented with a characteristics approach, where the characteristic variables are evaluated at zero viscosity. The same procedure is used to model seismic wave propagation at the interface between a fluid and porous medium. In this case, each medium is represented by a different grid and the two grids are combined through a domain-decomposition method. This wavefield decomposition method accounts for the discontinuity of variables and is crucial for an accurate interface treatment. We simulate seismic wave propagation with open-pore and sealed-pore boundary conditions and verify the validity and accuracy of the algorithm by comparing the numerical simulations to analytical solutions based on zero viscosity obtained with the Cagniard-de Hoop method. Finally, we illustrate the suitability of our algorithm for more complex models of porous media involving viscous pore fluids and strongly heterogeneous distributions of the elastic and hydraulic material properties.

New horizons for black holes and branes

We initiate a systematic scan of the landscape of black holes in any spacetime dimension using the recently proposed blackfold effective worldvolume theory. We focus primarily on asymptotically flat stationary vacuum solutions, where we uncover large classes of new black holes. These include helical black strings and black rings, black odd-spheres, for which the horizon is a product of a large and a small sphere, and non-uniform black cylinders. More exotic possibilities are also outlined. The blackfold description recovers correctly the ultraspinning Myers-Perry black holes as ellipsoidal even-ball configurations where the velocity field approaches the speed of light at the boundary of the ball. Helical black ring solutions provide the first instance of asymptotically flat black holes in more than four dimensions with a single spatial U(1) isometry. They also imply infinite rational non-uniqueness in ultraspinning regimes, where they maximize the entropy among all stationary single-horizon solutions. Moreover, static blackfolds are possible with the geometry of minimal surfaces. The absence of compact embedded minimal surfaces in Euclidean space is consistent with the uniqueness theorem of static black holes

Transonic flow past on enlarged missile nose parts

Numerical computation of a viscid heat-conducting transonic flow over a generic commercial rocket profile with symmetric oversized nose part was carried out. It has been shown that at zero angle of attack for some free-streamvelocity value flow pattern loses its symmetry. This results in non-uniform pressure distribution on rocket surface in angle direction which may yield in additional oscillating stress on the rocket. Also it has been found that obtained non-symmetric flow patterns are stable for small velocity perturbations.

Experimental and Numerical Analysis of Different Volutes in a Centrifugal Compressor

In a centrifugal compressor the flow around the diffuser is collected and led to the pipe system by a spiral-shaped volute. In this study a single-stage centrifugal compressor with three different volutes is investigated. The compressorwas first equipped with the original volute, the cross-section of which was a combination of a rectangle and semi-circle. Next a new volute with a fully circular cross-section was designed and manufactured. Finally, the circular volute wasmodified by rounding the tongue and smoothing the tongue area. The overall performance of the compressor as well as the static pressure distribution after the impeller and on the volute surface were measured. The flow entering the volute was measured using a three-hole Cobra-probe, and flow visualisations were carriedout in the exit cone of the volute. In addition, the radial force acting on theimpeller was measured using magnetic bearings. The complete compressor with thecircular volute (inlet pipe, full impeller, diffuser, volute and outlet pipe) was also modelled using computational fluid dynamics (CFD). A fully 3-D viscous flow was solved using a Navier-Stokes solver, Finflo, developed at Helsinki University of Technology. Chien's k-e model was used to take account of the turbulence. The differences observed in the performance of the different volutes were quite small. The biggest differences were at low speeds and high volume flows,i.e. when the flow entered the volute most radially. In this operating regime the efficiency of the compressor with the modified circular volute was about two percentage points higher than with the other volutes. Also, according to the Cobra-probe measurements and flow visualisations, the modified circular volute performed better than the other volutes in this operating area. The circumferential static pressure distribution in the volute showed increases at low flow, constant distribution at the design flow and decrease at high flow. The non-uniform static pressure distribution of the volute was transmitted backwards across the vaneless diffuser and observed at the impeller exit. At low volume flow a strong two-wave pattern developed into the static pressure distribution at the impeller exit due to the response of the impeller to the non-uniformity of pressure. The radial force of the impeller was the greatest at the choke limit, the smallest atthe design flow, and moderate at low flow. At low flow the force increase was quite mild, whereas the increase at high flow was rapid. Thus, the non-uniformityof pressure and the force related to it are strong especially at high flow. Theforce caused by the modified circular volute was weaker at choke and more symmetric as a function of the volume flow than the force caused by the other volutes.

Spontaneous Emergence of Chirality in the Limited Enantioselectivity Model: Autocatalytic Cycle Driven by an External Reagent

The emergence of chirality in enantioselective autocatalysis for compounds unable to transform according to the Frank-like reaction network is discussed with respect to the controversial limited enantioselectivity (LES) model composed of coupled enantioselective and non-enantioselective autocatalyses. The LES model cannot lead to spontaneous mirror symmetry breaking (SMSB) either in closed systems with a homogeneous temperature distribution or in closed systems with a stationary non-uniform temperature distribution. However, simulations of chemical kinetics in a two-compartment model demonstrate that SMSB may occur if both autocatalytic reactions are spatially separated at different temperatures in different compartments but coupled under the action of a continuous internal flow. In such conditions, the system can evolve, for certain reaction and system parameters, toward a chiral stationary state; that is, the system is able to reach a bifurcation point leading to SMSB. Numerical simulations in which reasonable chemical parameters have been used suggest that an adequate scenario for such a SMSB would be that of abyssal hydrothermal vents, by virtue of the typical temperature gradients found there and the role of inorganic solids mediating chemical reactions in an enzyme-like role.

Correcting surface coil excitation inhomogeneities in single-shot SPEN MRI.

Given their high sensitivity and ability to limit the field of view (FOV), surface coils are often used in magnetic resonance spectroscopy (MRS) and imaging (MRI). A major downside of surface coils is their inherent radiofrequency (RF) B1 heterogeneity across the FOV, decreasing with increasing distance from the coil and giving rise to image distortions due to non-uniform spatial responses. A robust way to compensate for B1 inhomogeneities is to employ adiabatic inversion pulses, yet these are not well adapted to all imaging sequences - including to single-shot approaches like echo planar imaging (EPI). Hybrid spatiotemporal encoding (SPEN) sequences relying on frequency-swept pulses provide another ultrafast MRI alternative, that could help solve this problem thanks to their built-in heterogeneous spatial manipulations. This study explores how this intrinsic SPEN-based spatial discrimination, could be used to compensate for the B1 inhomogeneities inherent to surface coils. Experiments carried out in both phantoms and in vivo rat brains demonstrate that, by suitably modulating the amplitude of a SPEN chirp pulse that progressively excites the spins in a direction normal to the coil, it is possible to compensate for the RF transmit inhomogeneities and thus improve sensitivity and image fidelity.

From Regional Landslide Detection to Site-Specific Slope Deformation Monitoring and Modelling Based on Active Remote Sensors

Landslide processes can have direct and indirect consequences affecting human lives and activities. In order to improve landslide risk management procedures, this PhD thesis aims to investigate capabilities of active LiDAR and RaDAR sensors for landslides detection and characterization at regional scales, spatial risk assessment over large areas and slope instabilities monitoring and modelling at site-specific scales. At regional scales, we first demonstrated recent boat-based mobile LiDAR capabilities to model topography of the Normand coastal cliffs. By comparing annual acquisitions, we validated as well our approach to detect surface changes and thus map rock collapses, landslides and toe erosions affecting the shoreline at a county scale. Then, we applied a spaceborne InSAR approach to detect large slope instabilities in Argentina. Based on both phase and amplitude RaDAR signals, we extracted decisive information to detect, characterize and monitor two unknown extremely slow landslides, and to quantify water level variations of an involved close dam reservoir. Finally, advanced investigations on fragmental rockfall risk assessment were conducted along roads of the Val de Bagnes, by improving approaches of the Slope Angle Distribution and the FlowR software. Therefore, both rock-mass-failure susceptibilities and relative frequencies of block propagations were assessed and rockfall hazard and risk maps could be established at the valley scale. At slope-specific scales, in the Swiss Alps, we first integrated ground-based InSAR and terrestrial LiDAR acquisitions to map, monitor and model the Perraire rock slope deformation. By interpreting both methods individually and originally integrated as well, we therefore delimited the rockslide borders, computed volumes and highlighted non-uniform translational displacements along a wedge failure surface. Finally, we studied specific requirements and practical issues experimented on early warning systems of some of the most studied landslides worldwide. As a result, we highlighted valuable key recommendations to design new reliable systems; in addition, we also underlined conceptual issues that must be solved to improve current procedures. To sum up, the diversity of experimented situations brought an extensive experience that revealed the potential and limitations of both methods and highlighted as well the necessity of their complementary and integrated uses.

3D Surface Modeling From Point Clouds

The goal of this thesis is to implement software for creating 3D models from point clouds. Point clouds are acquired with stereo cameras, monocular systems or laser scanners. The created 3D models are triangular models or NURBS (Non-Uniform Rational B-Splines) models. Triangular models are constructed from selected areas from the point clouds and resulted triangular models are translated into a set of quads. The quads are further translated into an estimated grid structure and used for NURBS surface approximation. Finally, we have a set of NURBS surfaces which represent the whole model. The problem wasn’t so easy to solve. The selected triangular surface reconstruction algorithm did not deal well with noise in point clouds. To handle this problem, a clustering method is introduced for simplificating the model and removing noise. As we had better results with the smaller point clouds produced by clustering, we used points in clusters to better estimate the grids for NURBS models. The overall results were good when the point cloud did not have much noise. The point clouds with small amount of error had good results as the triangular model was solid. NURBS surface reconstruction performed well on solid models.

Deformation of Gabor systems

We introduce a new notion for the deformation of Gabor systems. Such deformations are in general nonlinear and, in particular, include the standard jitter error and linear deformations of phase space. With this new notion we prove a strong deformation result for Gabor frames and Gabor Riesz sequences that covers the known perturbation and deformation results. Our proof of the deformation theorem requires a new characterization of Gabor frames and Gabor Riesz sequences. It is in the style of Beurling's characterization of sets of sampling for bandlimited functions and extends significantly the known characterization of Gabor frames 'without inequalities' from lattices to non-uniform sets.

Construção de câmara de luz ultravioleta para fotopolimerização de fases estacionárias monolíticas

The assembly of a photochemical reactor with six fluorescent lamps, used for photopolymerizations is described. This chamber presents a mobile support, allowing the placement of samples at different heights and a safety lock that interrupts the radiation, if it is opened during operation. The mirrored internal walls avoid the dispersion and non-uniform distribution of light. There is no high heating because the own character of the used lamps. All parts could be purchased in commerce with less than U$ 150,00. This reactor was successfully used for monolithic stationary phase photopolymerization.

Computational Evaluation of Print Unevenness According to Human Vision

Print quality and the printability of paper are very important attributes when modern printing applications are considered. In prints containing images, high print quality is a basic requirement. Tone unevenness and non uniform glossiness of printed products are the most disturbing factors influencing overall print quality. These defects are caused by non ideal interactions of paper, ink and printing devices in high speed printing processes. Since print quality is a perceptive characteristic, the measurement of unevenness according to human vision is a significant problem. In this thesis, the mottling phenomenon is studied. Mottling is a printing defect characterized by a spotty, non uniform appearance in solid printed areas. Print mottle is usually the result of uneven ink lay down or non uniform ink absorption across the paper surface, especially visible in mid tone imagery or areas of uniform color, such as solids and continuous tone screen builds. By using existing knowledge on visual perception and known methods to quantify print tone variation, a new method for print unevenness evaluation is introduced. The method is compared to previous results in the field and is supported by psychometric experiments. Pilot studies are made to estimate the effect of optical paper characteristics prior to printing, on the unevenness of the printed area after printing. Instrumental methods for print unevenness evaluation have been compared and the results of the comparison indicate that the proposed method produces better results in terms of visual evaluation correspondence. The method has been successfully implemented as ail industrial application and is proved to be a reliable substitute to visual expertise.

Application of Computational Fluid Dynamics in Modeling Blood Flow in Human Thoracic Aorta

The aim of this study was to simulate blood flow in thoracic human aorta and understand the role of flow dynamics in the initialization and localization of atherosclerotic plaque in human thoracic aorta. The blood flow dynamics in idealized and realistic models of human thoracic aorta were numerically simulated in three idealized and two realistic thoracic aorta models. The idealized models of thoracic aorta were reconstructed with measurements available from literature, and the realistic models of thoracic aorta were constructed by image processing Computed Tomographic (CT) images. The CT images were made available by South Karelia Central Hospital in Lappeenranta. The reconstruction of thoracic aorta consisted of operations, such as contrast adjustment, image segmentations, and 3D surface rendering. Additional design operations were performed to make the aorta model compatible for the numerical method based computer code. The image processing and design operations were performed with specialized medical image processing software. Pulsatile pressure and velocity boundary conditions were deployed as inlet boundary conditions. The blood flow was assumed homogeneous and incompressible. The blood was assumed to be a Newtonian fluid. The simulations with idealized models of thoracic aorta were carried out with Finite Element Method based computer code, while the simulations with realistic models of thoracic aorta were carried out with Finite Volume Method based computer code. Simulations were carried out for four cardiac cycles. The distribution of flow, pressure and Wall Shear Stress (WSS) observed during the fourth cardiac cycle were extensively analyzed. The aim of carrying out the simulations with idealized model was to get an estimate of flow dynamics in a realistic aorta model. The motive behind the choice of three aorta models with distinct features was to understand the dependence of flow dynamics on aorta anatomy. Highly disturbed and nonuniform distribution of velocity and WSS was observed in aortic arch, near brachiocephalic, left common artery, and left subclavian artery. On the other hand, the WSS profiles at the roots of branches show significant differences with geometry variation of aorta and branches. The comparison of instantaneous WSS profiles revealed that the model with straight branching arteries had relatively lower WSS compared to that in the aorta model with curved branches. In addition to this, significant differences were observed in the spatial and temporal profiles of WSS, flow, and pressure. The study with idealized model was extended to study blood flow in thoracic aorta under the effects of hypertension and hypotension. One of the idealized aorta models was modified along with the boundary conditions to mimic the thoracic aorta under the effects of hypertension and hypotension. The results of simulations with realistic models extracted from CT scans demonstrated more realistic flow dynamics than that in the idealized models. During systole, the velocity in ascending aorta was skewed towards the outer wall of aortic arch. The flow develops secondary flow patterns as it moves downstream towards aortic arch. Unlike idealized models, the distribution of flow was nonplanar and heavily guided by the artery anatomy. Flow cavitation was observed in the aorta model which was imaged giving longer branches. This could not be properly observed in the model with imaging containing a shorter length for aortic branches. The flow circulation was also observed in the inner wall of the aortic arch. However, during the diastole, the flow profiles were almost flat and regular due the acceleration of flow at the inlet. The flow profiles were weakly turbulent during the flow reversal. The complex flow patterns caused a non-uniform distribution of WSS. High WSS was distributed at the junction of branches and aortic arch. Low WSS was distributed at the proximal part of the junction, while intermedium WSS was distributed in the distal part of the junction. The pulsatile nature of the inflow caused oscillating WSS at the branch entry region and inner curvature of aortic arch. Based on the WSS distribution in the realistic model, one of the aorta models was altered to induce artificial atherosclerotic plaque at the branch entry region and inner curvature of aortic arch. Atherosclerotic plaque causing 50% blockage of lumen was introduced in brachiocephalic artery, common carotid artery, left subclavian artery, and aortic arch. The aim of this part of the study was first to study the effect of stenosis on flow and WSS distribution, understand the effect of shape of atherosclerotic plaque on flow and WSS distribution, and finally to investigate the effect of lumen blockage severity on flow and WSS distributions. The results revealed that the distribution of WSS is significantly affected by plaque with mere 50% stenosis. The asymmetric shape of stenosis causes higher WSS in branching arteries than in the cases with symmetric plaque. The flow dynamics within thoracic aorta models has been extensively studied and reported here. The effects of pressure and arterial anatomy on the flow dynamic were investigated. The distribution of complex flow and WSS is correlated with the localization of atherosclerosis. With the available results we can conclude that the thoracic aorta, with complex anatomy is the most vulnerable artery for the localization and development of atherosclerosis. The flow dynamics and arterial anatomy play a role in the localization of atherosclerosis. The patient specific image based models can be used to diagnose the locations in the aorta vulnerable to the development of arterial diseases such as atherosclerosis.

Performance evaluation of 3D computer vision techniques

This work presents the implementation and comparison of three different techniques of three-dimensional computer vision as follows: • Stereo vision - correlation between two 2D images • Sensorial fusion - use of different sensors: camera 2D + ultrasound sensor (1D); • Structured light The computer vision techniques herein presented took into consideration the following characteristics: • Computational effort ( elapsed time for obtain the 3D information); • Influence of environmental conditions (noise due to a non uniform lighting, overlighting and shades); • The cost of the infrastructure for each technique; • Analysis of uncertainties, precision and accuracy. The option of using the Matlab software, version 5.1, for algorithm implementation of the three techniques was due to the simplicity of their commands, programming and debugging. Besides, this software is well known and used by the academic community, allowing the results of this work to be obtained and verified. Examples of three-dimensional vision applied to robotic assembling tasks ("pick-and-place") are presented.

Weed control via intercropping with gliricidia: I. cotton crop

The majority of cotton grown commercially in the world has white lint, but recently, there has been a growing interest in colored lint cotton in several countries, including Brazil. The use of naturally-colored fiber reduces chemical pollution. The objective of this paper was to evaluate cotton cultivar fiber yield in response to weed control via intercropping with gliricídia. Cultivars BRS-Verde (greenish fibers), BRS-Rubi (reddish brown fibers), BRS-Safira (brown fibers), and BRS-187 8H (white fibers) were submitted to the following treatments: no hoeing, two hoeings (at 20 and 40 days after transplanting), and cotton intercropped with gliricídia. In the intercropped treatment, gliricídia was planted between rows of cotton plants, using one seedling pit-1, in pits spaced 50.0 cm apart. Twelve weed species predominated in the experiment, many of them belonging to the Poaceae family. Weeds occurred at different frequencies and in a non-uniform manner in the experimental area. Cultivars did not influence weed dry matter. Intercropping with gliricídia reduced weed dry matter but did not prevent reductions in cotton fiber and seed cotton yield, which were higher in hoed plots. Cultivar BRS Safira had the highest fiber yield, but no differences were observed between cultivars regarding to seed cotton yield.