Geometric integration applied to moving mesh methods and degenerate Lagrangians

Moving mesh methods (also called r-adaptive methods) are space-adaptive strategies used for the numerical simulation of time-dependent partial differential equations. These methods keep the total number of mesh points fixed during the simulation, but redistribute them over time to follow the areas where a higher mesh point density is required. There are a very limited number of moving mesh methods designed for solving field-theoretic partial differential equations, and the numerical analysis of the resulting schemes is challenging. In this thesis we present two ways to construct r-adaptive variational and multisymplectic integrators for (1+1)-dimensional Lagrangian field theories. The first method uses a variational discretization of the physical equations and the mesh equations are then coupled in a way typical of the existing r-adaptive schemes. The second method treats the mesh points as pseudo-particles and incorporates their dynamics directly into the variational principle. A user-specified adaptation strategy is then enforced through Lagrange multipliers as a constraint on the dynamics of both the physical field and the mesh points. We discuss the advantages and limitations of our methods. The proposed methods are readily applicable to (weakly) non-degenerate field theories---numerical results for the Sine-Gordon equation are presented.

In an attempt to extend our approach to degenerate field theories, in the last part of this thesis we construct higher-order variational integrators for a class of degenerate systems described by Lagrangians that are linear in velocities. We analyze the geometry underlying such systems and develop the appropriate theory for variational integration. Our main observation is that the evolution takes place on the primary constraint and the 'Hamiltonian' equations of motion can be formulated as an index 1 differential-algebraic system. We then proceed to construct variational Runge-Kutta methods and analyze their properties. The general properties of Runge-Kutta methods depend on the 'velocity' part of the Lagrangian. If the 'velocity' part is also linear in the position coordinate, then we show that non-partitioned variational Runge-Kutta methods are equivalent to integration of the corresponding first-order Euler-Lagrange equations, which have the form of a Poisson system with a constant structure matrix, and the classical properties of the Runge-Kutta method are retained. If the 'velocity' part is nonlinear in the position coordinate, we observe a reduction of the order of convergence, which is typical of numerical integration of DAEs. We also apply our methods to several models and present the results of our numerical experiments.

The lakes of Michoacan (Mexico): A brief history and alternative point of view

This article is intended to open a discussion about the historical development of lakes Zirahuen, Patzcuaro and Cuitzeo in the state of Michoacan, and the postulated relationships between lake ecology and evolution. Dr Fernando De Buen was the first man dedicated to limnology in Mexico who came to the country in the 1930s. He was adviser at the Estacion Limnologica de Patzcuaro and wrote outstanding papers dealing with Mexican lakes. The lakes of Michoacan probably formed in the late Pliocene or Holocene, and were part of a tributary to the Lerma River, which became isolated by successive volanic barriers to form lake basins. Lake Zirahuen is a warm monomictic waterbody with unique water dynamics amongst the Michoacan lakes. Because it is relatively deep (max depth 40m), seasonal patterns of alternating circulation and thermal stratification develop in the lake, a feature not shared by the other two polymictic shallow lakes, Patzcuaro and Cuitzeo.

Construction of the refractive index profiles for few-mode planar optical waveguides

We present a nondestructive technique to predict the refractive index profiles of isotropic planar waveguides, on which a thin gold film is deposited to as the cladding. The negative dielectric constant of the metal results in significant differences of effective indices between TE and TM modes. The two polarized modes and a surface plasmon resonance (SPR) with abundant information of the surface index can be used to construct the refractive index profiles of single-mode and two-mode waveguides at a fixed wavelength. (c) 2005 Elsevier B.V. All rights reserved.

An Index of Abundance for Coastal Species of Juvenile Sharks from the Northeast Gulf of Mexico

A fishery-independent assessment of juvenile coastal shark populations in U.S. waters of the northeast Gulf of Mexico was conducted using two methods: gillnets and longlines. Surveys were conducted monthly during April–October in two fixed sampling areas from 1996 to 1998. The Atlantic sharpnose shark, Rhizoprionodon terraenovae, and the blacktip shark, Carcharhinus limbatus, were the most common species captured with either longlines or gillnets. An additional 14 shark species were captured, and juvenile indices of abundance were developed for 8 species with gillnets and 6 species of sharks with longlines. Trends in catch-per-unit-effort were found to vary depending on species. Length-frequency information revealed that the majority of sharks captured were juveniles. Given the direct relationship between stock and recruitment for sharks, continued monitoring of juvenile abundance will aid in determining the strength of the parental stock size and for predicting future population strength.

Multi scale shape index for 3D object recognition

We present Multi Scale Shape Index (MSSI), a novel feature for 3D object recognition. Inspired by the scale space filtering theory and Shape Index measure proposed by Koenderink & Van Doorn [6], this feature associates different forms of shape, such as umbilics, saddle regions, parabolic regions to a real valued index. This association is useful for representing an object based on its constituent shape forms. We derive closed form scale space equations which computes a characteristic scale at each 3D point in a point cloud without an explicit mesh structure. This characteristic scale is then used to estimate the Shape Index. We quantitatively evaluate the robustness and repeatability of the MSSI feature for varying object scales and changing point cloud density. We also quantify the performance of MSSI for object category recognition on a publicly available dataset. © 2013 Springer-Verlag.

生物质气化与废弃物焚烧联合发电的技术经济分析

介绍了生物质气化与废弃物焚烧联合发电技术项目，确定了该项目经济效益的评价指标，定量计算了项目的投资回收期、净现值和内部收益率。同时还对燃料费用、上网电价和固定资产变化引起的敏感性进行分析。结果表明，该联合发电技术的动态投资回收期为9.05a，净现值为2770万元，内部收益率为15.82％，三个经济指标均符合行业标准。从经济角度看是完全可行的。

Study on QSPR of alcohols with a novel edge connectivity index F-m

A novel edge degree f(i) for heteroatom and multiple bonds in molecular graph is derived on the basis of the edge degree delta(e(r)). A novel edge connectivity index F-m is introduced. The multiple linear regression by using the edge connectivity index F-m and alcohol-type parameter delta, alcohol-distance parameter L can provide high-quality QSPR models for the normal boiling points (BPs), molar volumes (MVs), molar refraction (MRs), water solubility(log(1/S)) and octanol/water partition (logP) of alcohols with up to 17 non-hydrogen atoms. The results imply that these physical properties may be expressed as a liner combination of the edge connectivity index and alcohol-type parameter, 6, alcohol-distance parameter, L. For the models of the five properties, the correlation coefficient r and the standard errors are 0.9969,3.022; 0.9993, 1.504; 0.9992, 0.446; 0.9924,0.129 and 0.9973,0.123 for BPs, MVs, MRs, log(1/S) and logP, respectively. The cross-validation by using the leave-one-out method demonstrates the models to be highly reliable from the point of view of statistics.

Monitoring the vegetation start of season (SOS) across the island of Ireland using the MERIS global vegetation index

The aim of this study was to develop a methodology, based on satellite remote sensing, to estimate the vegetation Start of Season (SOS) across the whole island of Ireland on an annual basis. This growing body of research is known as Land Surface Phenology (LSP) monitoring. The SOS was estimated for each year from a 7-year time series of 10-day composited, 1.2 km reduced resolution MERIS Global Vegetation Index (MGVI) data from 2003 to 2009, using the time series analysis software, TIMESAT. The selection of a 10-day composite period was guided by in-situ observations of leaf unfolding and cloud cover at representative point locations on the island. The MGVI time series was smoothed and the SOS metric extracted at a point corresponding to 20% of the seasonal MGVI amplitude. The SOS metric was extracted on a per pixel basis and gridded for national scale coverage. There were consistent spatial patterns in the SOS grids which were replicated on an annual basis and were qualitatively linked to variation in landcover. Analysis revealed that three statistically separable groups of CORINE Land Cover (CLC) classes could be derived from differences in the SOS, namely agricultural and forest land cover types, peat bogs, and natural and semi-natural vegetation types. These groups demonstrated that managed vegetation, e.g. pastures has a significantly earlier SOS than in unmanaged vegetation e.g. natural grasslands. There was also interannual spatio-temporal variability in the SOS. Such variability was highlighted in a series of anomaly grids showing variation from the 7-year mean SOS. An initial climate analysis indicated that an anomalously cold winter and spring in 2005/2006, linked to a negative North Atlantic Oscillation index value, delayed the 2006 SOS countrywide, while in other years the SOS anomalies showed more complex variation. A correlation study using air temperature as a climate variable revealed the spatial complexity of the air temperature-SOS relationship across the Republic of Ireland as the timing of maximum correlation varied from November to April depending on location. The SOS was found to occur earlier due to warmer winters in the Southeast while it was later with warmer winters in the Northwest. The inverse pattern emerged in the spatial patterns of the spring correlates. This contrasting pattern would appear to be linked to vegetation management as arable cropping is typically practiced in the southeast while there is mixed agriculture and mostly pastures to the west. Therefore, land use as well as air temperature appears to be an important determinant of national scale patterns in the SOS. The TIMESAT tool formed a crucial component of the estimation of SOS across the country in all seven years as it minimised the negative impact of noise and data dropouts in the MGVI time series by applying a smoothing algorithm. The extracted SOS metric was sensitive to temporal and spatial variation in land surface vegetation seasonality while the spatial patterns in the gridded SOS estimates aligned with those in landcover type. The methodology can be extended for a longer time series of FAPAR as MERIS will be replaced by the ESA Sentinel mission in 2013, while the availability of full resolution (300m) MERIS FAPAR and equivalent sensor products holds the possibility of monitoring finer scale seasonality variation. This study has shown the utility of the SOS metric as an indicator of spatiotemporal variability in vegetation phenology, as well as a correlate of other environmental variables such as air temperature. However, the satellite-based method is not seen as a replacement of ground-based observations, but rather as a complementary approach to studying vegetation phenology at the national scale. In future, the method can be extended to extract other metrics of the seasonal cycle in order to gain a more comprehensive view of seasonal vegetation development.

An EORTC phase I study of epirubicin in combination with fixed doses of cyclophosphamide and infusional 5-fu (CEF-infu) as primary treatment of large operable or locally advanced/inflammatory breast cancer.

PURPOSE: The association of continuous infusion 5-fluorouracil, epirubicin (50 mg/m2 q 3 weeks) and a platinum compound (cisplatin or carboplatin) was found to be very active in patients with either locally advanced/inflammatory (LA/I) [1, 2] or large operable (LO) breast cancer (BC) [3]. The same rate of activity in terms of response rate (RR) and response duration was observed in LA/I BC patients when cisplatin was replaced by cyclophosphamide [4]. The dose of epirubicin was either 50 mg/m2 [ 1, 2, 3] or 60 mg/m2/cycle [4]. The main objective of this study was to determine the maximum tolerated dose (MTD) of epirubicin when given in combination with fixed doses of cyclophosphamide and infusional 5-fluorouracil (CEF-infu) as neoadjuvant therapy in patients with LO or LA/I BC for a maximum of 6 cycles. PATIENTS AND METHODS: Eligible patients had LO or LA/I BC, a performance status 0-1, adequate organ function and were <65 years old. Cyclophosphamide was administered at the dose of 400 mg/m2 day 1 and 8, q 4 weeks and infusional 5-fluorouracil 200 mg/m2/day was given day 1-28, q 4 weeks. Epirubicin was escalated from 30 to 45 and to 60 mg/m2 day 1 and 8; dose escalation was permitted if 0/3 or 1/6 patients experienced dose limiting toxicity (DLT) during the first 2 cycles of therapy. DLT for epirubicin was defined as febrile neutropenia, grade 4 neutropenia lasting for >7 days, grade 4 thrombocytopenia, or any non-haematological toxicity of CTC grade > or =3, excluding alopecia and plantar-palmar erythrodysesthesia (this toxicity was attributable to infusional 5-fluorouracil and was not considered a DLT of epirubicin). RESULTS: A total of 21 patients, median age 44 years (range 29-63) have been treated. 107 courses have been delivered, with a median number of 5 cycles per patient (range 4-6). DLTs on cycles I and 2 on level 1, 2, 3: grade 3 (G3) mucositis occurred in 1/10 patients treated at the third dose level. An interim analysis showed that G3 PPE occurred in 5/16 pts treated with the 28-day infusional 5-FU schedule at the 3 dose levels. The protocol was subsequently amended to limit the duration of infusional 5-fluorouracil infusion from 4 to 3 weeks. No G3 PPE was detected in 5 patients treated with this new schedule. CONCLUSIONS: This study establishes that epirubicin 60mg/m2 day 1 and 8, cyclophosphamide 400mg/m2 day 1 and 8 and infusional 5-fluorouracil 200 mg/m2/day day 1-21. q 4 weeks is the recommended dose level. Given the encouraging activity of this regimen (15/21 clinical responses) we have replaced infusional 5-fluorouracil by oral capecitabine in a recently activated study.

A cross-site, comparative effectiveness study of an integrated HIV and substance use treatment program.

Co-occurrence of HIV and substance abuse is associated with poor outcomes for HIV-related health and substance use. Integration of substance use and medical care holds promise for HIV patients, yet few integrated treatment models have been reported. Most of the reported models lack data on treatment outcomes in diverse settings. This study examined the substance use outcomes of an integrated treatment model for patients with both HIV and substance use at three different clinics. Sites differed by type and degree of integration, with one integrated academic medical center, one co-located academic medical center, and one co-located community health center. Participants (n=286) received integrated substance use and HIV treatment for 12 months and were interviewed at 6-month intervals. We used linear generalized estimating equation regression analysis to examine changes in Addiction Severity Index (ASI) alcohol and drug severity scores. To test whether our treatment was differentially effective across sites, we compared a full model including site by time point interaction terms to a reduced model including only site fixed effects. Alcohol severity scores decreased significantly at 6 and 12 months. Drug severity scores decreased significantly at 12 months. Once baseline severity variation was incorporated into the model, there was no evidence of variation in alcohol or drug score changes by site. Substance use outcomes did not differ by age, gender, income, or race. This integrated treatment model offers an option for treating diverse patients with HIV and substance use in a variety of clinic settings. Studies with control groups are needed to confirm these findings.

The Sports Overload Injuries of the Lower Extremities from the Biomechanical Point of View

This work studies the major sports overload injuries of the lower extremities from the biomechanical point of view. At the same time, the main paradigms of podiatric biomechanics and the application of new biomechanical theories in the study of these lesions are reviewed. With current legislation, clinical gait biomechanical studies should be carried out in health centres and the only health professionals who can perform them are podiatrists and doctors (because they both can diagnose). Graduates in physical education can carry out studies in the field or in the sports court for the sole purpose of improving athletic performance, but never intended to treat a pathology overload.

Using idiothetic cues to swim a path with a fixed trajectory and distance: Necessary involvement of the hippocampus but not the retrosplenial cortex.

Rats rapidly learned to find a submerged platform in a water maze at a constant distance and angle from the start point, which changed on every trial. The rats performed accurately in the light and dark, but prior rotation disrupted the latter condition. The rats were then retested after receiving cytotoxic hippocampal or retrosplenial cortex lesions. Retrosplenial lesions had no apparent effect in either the light or dark. Hippocampal lesions impaired performance in both conditions but spared the ability to locate a platform placed in the center of the pool. A hippocampal deficit emerged when this pool-center task was run in the dark. The spatial effects of hippocampal damage extend beyond allocentric tasks to include aspects of idiothetic guidance.

Long-term occlusal and soft-tissue profile outcomes after treatment of Class II Division 1 malocclusion with fixed appliances

Introduction: Our objective was to determine which factors were predictive of good long-term outcomes after fixed appliance treatment of Class II Division 1 malocclusion. Methods: Two hundred seven patients with Class II Division 1 malocclusion were examined in early adulthood at a mean of 4.6 years after treatment with fixed appliances. The peer assessment rating index was used to evaluate dental alignment and occlusal relationships. The soft-tissue profile was assessed with the Holdaway angle. Results: Logistic regression identified 3 pretreatment variables that were predictive of a good facial profile (Holdaway angle) at recall: the lower lip to E-plane distance (P

Developing an integrated hydrograph separation and lumped modelling approach to quantifying hydrological pathways in Irish river catchments

An appreciation of the quantity of streamflow derived from the main hydrological pathways involved in transporting diffuse contaminants is critical when addressing a wide range of water resource management issues. In order to assess hydrological pathway contributions to streams, it is necessary to provide feasible upper and lower bounds for flows in each pathway. An important first step in this process is to provide reliable estimates of the slower responding groundwater pathways and subsequently the quicker overland and interflow pathways. This paper investigates the effectiveness of a multi-faceted approach applying different hydrograph separation techniques, supplemented by lumped hydrological modelling, for calculating the Baseflow Index (BFI), for the development of an integrated approach to hydrograph separation. A semi-distributed, lumped and deterministic rainfall runoff model known as NAM has been applied to ten catchments (ranging from 5 to 699 km2). While this modelling approach is useful as a validation method, NAM itself is also an important tool for investigation. These separation techniques provide a large variation in BFI, a difference of 0.741 predicted for BFI in a catchment with the less reliable fixed and sliding interval methods and local minima turning point methods included. This variation is reduced to 0.167 with these methods omitted. The Boughton and Eckhardt algorithms, while quite subjective in their use, provide quick and easily implemented approaches for obtaining physically realistic hydrograph separations. It is observed that while the different separation techniques give varying BFI values for each of the catchments, a recharge coefficient approach developed in Ireland, when applied in conjunction with the Master recession Curve Tabulation method, predict estimates in agreement with those obtained using the NAM model, and these estimates are also consistent with the study catchments’ geology. These two separation methods, in conjunction with the NAM model, were selected to form an integrated approach to assessing BFI in catchments.

Analysis of Fixed-Speed Wind Farm Low-Frequency Power Pulsations using a Wavelet-Prony Method

The increasing penetration of wind generation on the Island of Ireland has been accompanied by close investigation of low-frequency periodic pulsations contained within the active power flow from different wind farms. A primary concern is excitation of existing low-frequency oscillation modes already present on the system, particularly the 0.75 Hz mode as a consequence of the interconnected Northern and Southern power system networks. Recently grid code requirements on the Northern Ireland power system have been updated stipulating that wind farms connected after 2005 must be able to control the magnitude of oscillations in the range of 0.25 - 1.75 Hz to within 1% of the wind farm's registered output. In order to determine whether wind farm low-frequency oscillations have a negative effect (excite other modes) or possibly a positive impact (damping of existing modes) on the power system, the oscillations at the point of connection must be measured and characterised. Using time - frequency methods, research presented in this paper has been conducted to extract signal features from measured low-frequency active power pulsations produced by wind farms to determine the effective composition of possible oscillatory modes which may have a detrimental effect on system dynamic stability. The paper proposes a combined wavelet-Prony method to extract modal components and determine damping factors. The method is exemplified using real data obtained from wind farm measurements.