Vehicle Routing Problems that Minimize the Completion Time: Heuristics, Worst-Case Analyses, and Computational Results

In the standard Vehicle Routing Problem (VRP), we route a fleet of vehicles to deliver the demands of all customers such that the total distance traveled by the fleet is minimized. In this dissertation, we study variants of the VRP that minimize the completion time, i.e., we minimize the distance of the longest route. We call it the min-max objective function. In applications such as disaster relief efforts and military operations, the objective is often to finish the delivery or the task as soon as possible, not to plan routes with the minimum total distance. Even in commercial package delivery nowadays, companies are investing in new technologies to speed up delivery instead of focusing merely on the min-sum objective. In this dissertation, we compare the min-max and the standard (min-sum) objective functions in a worst-case analysis to show that the optimal solution with respect to one objective function can be very poor with respect to the other. The results motivate the design of algorithms specifically for the min-max objective. We study variants of min-max VRPs including one problem from the literature (the min-max Multi-Depot VRP) and two new problems (the min-max Split Delivery Multi-Depot VRP with Minimum Service Requirement and the min-max Close-Enough VRP). We develop heuristics to solve these three problems. We compare the results produced by our heuristics to the best-known solutions in the literature and find that our algorithms are effective. In the case where benchmark instances are not available, we generate instances whose near-optimal solutions can be estimated based on geometry. We formulate the Vehicle Routing Problem with Drones and carry out a theoretical analysis to show the maximum benefit from using drones in addition to trucks to reduce delivery time. The speed-up ratio depends on the number of drones loaded onto one truck and the speed of the drone relative to the speed of the truck.

Artificial fertilization of oocytes and sperm activation in pacu: effects of the spermatozoa:oocyte ratio, water volume, and in natura semen preservation

Objetivou-se foi avaliar a fertilização artificial e a duração da motilidade espermática em pacus com diferentes doses inseminantes, volumes de água e preservação do sêmen in natura. Foram realizados quatro experimentos para avaliação do efeito de doses inseminantes (7x10³, 7x10(4), 7x10(5), 7x10(6) e 7x10(7) espermatozoides ovócito-1) sobre a fertilização artificial dos ovócitos; do efeito do volume de água (0,5; 15,0; 30,0; 45,0 e 60,0 mL de água mL-1 de ovócitos) com doses inseminantes de 105.481 e 210.963 espermatozoides ovócito-1; do efeito de diluição do sêmen (0,005; 0,05; 0,5 e 5,0 µL de sêmen mL-1 de água) sobre a duração da motilidade espermática; e do efeito do armazenamento a 15 ºC por 9 h sobre a duração da motilidade espermática e o índice de sobrevivência espermática. Os maiores resultados obtidos foram: doses inseminantes entre 7x10³ e 7x10(7) espermatozoides ovócito-1; 15 a 60 mL de água mL-1 de ovócitos; diluição de 0.005 µL sêmen mL-1 de água e 98,65% de sobrevivência espermática até o tempo de preservação de 2h45min36s. A preservação a 15ºC por 9 horas não influencia a duração da motilidade espermática. As maiores taxas de fertilização podem ser observadas no emprego de 0,27 a 270 µL de sêmen mL-1 de ovócitos, com 15 a 60 mL de água para ativação.

Vortex Shedding Dynamics in Long Aspect-Ratio Aerodynamics Bodies

The focus of the current dissertation is to study qualitatively the underlying physics of vortex-shedding and wake dynamics in long aspect-ratio aerodynamics in incompressible viscous flow through the use of the KLE method. We carried out a long series of numerical experiments in the cases of flow around the cylinder at low Reynolds numbers. The study of flow at low Reynolds numbers provides an insight in the fluid physics and also plays a critical role when applying to stalled turbine rotors. Many of the conclusions about the qualitative nature of the physical mechanisms characterizing vortex formation, shedding and further interaction analyzed here at low Re could be extended to other Re regimes and help to understand the separation of the boundary layers in airfoils and other aerodynamic surfaces. In the long run, it aims to provide a better understanding of the complex multi-physics problems involving fluid-structure-control interaction through improved mathematical computational models of the multi-physics process. Besides the scientific conclusions produced, the research work on streamlined and bluff-body condition will also serve as a valuable guide for the future design of blade aerodynamics and the placement of wind turbines and hydrakinetic turbines, increasing the efficiency in the use of expensive workforce, supplies, and infrastructure. After the introductory section describing the main fields of application of wind power and hydrokinetic turbines, we describe the main features and theoretical background of the numerical method used here. Then, we present the analysis of the numerical experimentation results for the oscillatory regime right before the onset of vortex shedding for circular cylinders. We verified the wake length of the closed near-wake behind the cylinder and analysed the decay of the wake at the wake formation region, and then studied the St-Re relationship at the Reynolds numbers before the wake sheds compared to the experimental data. We found a theoretical model that describes the time evolution of the amplitude of fluctuations in the vorticity field on the twin vortex wake, which accurately matches the numerical results in terms of the frequency of the oscillation and rate of decay. We also proposed a model based on an analog circuit that is able to interpret the concerning flow by reducing the number of degrees of freedom. It follows the idea of the non-linear oscillator and resembles the dynamics mechanism of the closed near-wake with a common configured sine wave oscillator. This low-dimensional circuital model may also help to understand the underlying physical mechanisms, related to vorticity transport, that give origin to those oscillations.

Vortex evolution in the wake of accelerating low-aspect-ratio plates and three-dimensional bodies of revolution

To find examples of effecient locomotion and manoeuvrability, one need only turn to the elegant solutions natural flyers and swimmers have converged upon. This dissertation is specifically motivated by processes of evolutionary convergence, which have led to the propulsors and body shapes in nature that exhibit strong geometric collapse over diverse scales. These body features are abstracted in the studies presented herein using low-aspect-ratio at plates and a three-dimensional body of revolution (a sphere). The highly-separated vortical wakes that develop during accelerations are systematically characterized as a function of planform shape, aspect ratio, Reynolds number, and initial boundary conditions. To this end, force measurements and time-resolved (planar) particle image velocimetry have been used throughout to quantify the instantaneous forces and vortex evolution in the wake of the bluff bodies. During rectilinear motions, the wake development for the flat plates is primarily dependent on plate aspect ratio, with edge discontinuities and curvature playing only a secondary role. Furthermore, the axisymmetric case, i.e. the circular plate, shows strong sensitivity to Reynolds number, while this sensitivity quickly diminishes with increasing aspect ratio. For rotational motions, global insensitivity to plate aspect ratio has been observed. For the sphere, it has been shown that accelerations play an important role in the mitigation of flow separation. These results - expounded upon in this dissertation - have begun to shed light on the specific vortex dynamics that may be coopted by flying and swimming species of all shapes and sizes towards efficient locomotion.

Variability in the international normalised ratio (INR) in patients with antiphospholipid syndrome and positive lupus anticoagulant: should the INR targets be higher?

El síndrome antifosfolípido es un desorden autoinmune caracterizado por hipercoagulabilidad que requiere terapia anticoagulante como pilar fundamental, siendo la warfarina el tratamiento de elección en los casos que requieren manejo por largos periodos. Sin embargo, los pacientes con anticoagulante lúpico positivo representan un reto porque tienen mayor riesgo de presentar eventos trombóticos, sumado a que el seguimiento con el International Normalized Ratio (INR) no es confiable, ya que estos anticuerpos generan interferencia con las pruebas de laboratorio basadas en fosfolípidos, como es el caso del tiempo de protrombina (PT) con INR basal prolongado, incluso antes del inicio de la terapia anticoagulante. Por tal razón, se ilustra el caso de una paciente con síndrome antifosfolípido primario y anticoagulante lúpico positivo quien ha presentado múltiples episodios trombóticos, a pesar de recibir terapia anticoagulante. Además se hace una revisión de la literatura disponible y se postulan nuevas metas de INR en estos pacientes diferentes de las que se plantean actualmente.

Characterization of the volatile profile of Brazilian Merlot wines through comprehensive two dimensional gas chromatography time-of-flight mass spectrometric detection.

Wine aroma is an important characteristic and may be related to certain specific parameters, such as raw material and production process. The complexity of Merlot wine aroma was considered suitable for comprehensive two-dimensional gas chromatography (GCGC), as this technique offers superior performance when compared to one-dimensional gas chromatography (1D-GC). The profile of volatile compounds of Merlot wine was, for the first time, qualitatively analyzed by HS-SPME-GCxGC with a time-of-flight mass spectrometric detector (TOFMS), resulting in 179 compounds tentatively identified by comparison of experimental GCxGC retention indices and mass spectra with literature 1D-GC data and 155 compounds tentatively identified only by mass spectra comparison. A set of GCGC experimental retention indices was also, for the first time, presented for a specific inverse set of columns. Esters were present in higher number (94), followed by alcohols (80), ketones (29), acids (29), aldehydes (23), terpenes (23), lactones (16), furans (14), sulfur compounds (9), phenols (7), pyrroles (5), C13-norisoprenoids (3), and pyrans (2). GCxGC/TOFMS parameters were improved and optimal conditions were: a polar (polyethylene glycol)/medium polar (50% phenyl 50% dimethyl arylene siloxane) column set, oven temperature offset of 10ºC, 7 s as modulation period and 1.4 s of hot pulse duration. Co-elutions came up to 138 compounds in 1D and some of them were resolved in 2D. Among the coeluted compounds, thirty-three volatiles co-eluted in both 1D and 2D and their tentative identification was possible only due to spectral deconvolution. Some compounds that might have important contribution to aroma notes were included in these superimposed peaks. Structurally organized distribution of compounds in the 2D space was observed for esters, aldehydes and ketones, alcohols, thiols, lactones, acids and also inside subgroups, as occurred with esters and alcohols. The Fischer Ratio was useful for establishing the analytes responsible for the main differences between Merlot and non-Merlot wines. Differentiation among Merlot wines and wines of other grape varieties were mainly perceived through the following components: ethyl dodecanoate, 1-hexanol, ethyl nonanoate, ethyl hexanoate, ethyl decanoate, dehydro-2-methyl-3(2H)thiophenone, 3-methyl butanoic acid, ethyl tetradecanoate, methyl octanoate, 1,4 butanediol, and 6-methyloctan-1-ol.

Detailed analysis of a conservative difference approximation for the time fractional diffusion equation

Diffusion equations that use time fractional derivatives are attractive because they describe a wealth of problems involving non-Markovian Random walks. The time fractional diffusion equation (TFDE) is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α ∈ (0, 1). Developing numerical methods for solving fractional partial differential equations is a new research field and the theoretical analysis of the numerical methods associated with them is not fully developed. In this paper an explicit conservative difference approximation (ECDA) for TFDE is proposed. We give a detailed analysis for this ECDA and generate discrete models of random walk suitable for simulating random variables whose spatial probability density evolves in time according to this fractional diffusion equation. The stability and convergence of the ECDA for TFDE in a bounded domain are discussed. Finally, some numerical examples are presented to show the application of the present technique.

Physiological and symptomatic responses to cycling and walking in intermittent claudication

To shed light on the potential efficacy of cycling as a testing modality in the treatment of intermittent claudication (IC), this study compared physiological and symptomatic responses to graded walking and cycling tests in claudicants. Sixteen subjects with peripheral arterial disease (resting ankle: brachial index (ABI) < 0.9) and IC completed a maximal graded treadmill walking (T) and cycle (C) test after three familiarization tests on each mode. During each test, symptoms, oxygen uptake (VO2), minute ventilation (VE), respiratory exchange ratio (RER) and heart rate (HR) were measured, and for 10 min after each test the brachial and ankle systolic pressures were recorded. All but one subject experienced calf pain as the primary limiting symptom during T; whereas the symptoms were more varied during C and included thigh pain, calf pain and dyspnoea. Although maximal exercise time was significantly longer on C than T (690 +/- 67 vs. 495 +/- 57 s), peak VO2, peak VE and peak heart rate during C and T were not different; whereas peak RER was higher during C. These responses during C and T were also positively correlated (P < 0.05) with each other, with the exception of RER. The postexercise systolic pressures were also not different between C and T. However, the peak decline in ankle pressures from resting values after C and T were not correlated with each other. These data demonstrate that cycling and walking induce a similar level of metabolic and cardiovascular strain, but that the primary limiting symptoms and haemodynamic response in an individual's extremity, measured after exercise, can differ substantially between these two modes.

Determination of preventive maintenance lead time using hybrid analysis

The time for conducting Preventive Maintenance (PM) on an asset is often determined using a predefined alarm limit based on trends of a hazard function. In this paper, the authors propose using both hazard and reliability functions to improve the accuracy of the prediction particularly when the failure characteristic of the asset whole life is modelled using different failure distributions for the different stages of the life of the asset. The proposed method is validated using simulations and case studies.

Effective Cell-Centred Time-Domain Maxwell's Equations Numerical Solvers

This research work analyses techniques for implementing a cell-centred finite-volume time-domain (ccFV-TD) computational methodology for the purpose of studying microwave heating. Various state-of-the-art spatial and temporal discretisation methods employed to solve Maxwell's equations on multidimensional structured grid networks are investigated, and the dispersive and dissipative errors inherent in those techniques examined. Both staggered and unstaggered grid approaches are considered. Upwind schemes using a Riemann solver and intensity vector splitting are studied and evaluated. Staggered and unstaggered Leapfrog and Runge-Kutta time integration methods are analysed in terms of phase and amplitude error to identify which method is the most accurate and efficient for simulating microwave heating processes. The implementation and migration of typical electromagnetic boundary conditions. from staggered in space to cell-centred approaches also is deliberated. In particular, an existing perfectly matched layer absorbing boundary methodology is adapted to formulate a new cell-centred boundary implementation for the ccFV-TD solvers. Finally for microwave heating purposes, a comparison of analytical and numerical results for standard case studies in rectangular waveguides allows the accuracy of the developed methods to be assessed.