First-principles analysis of electron correlation, spin ordering and phonons in the normal state of FeSe1-x

We present first-principles density-functional-theory-based calculations to determine the effects of the strength of on-site electron correlation, magnetic ordering, pressure and Se vacancies on phonon frequencies and electronic structure of FeSe1-x. The theoretical equilibrium structure (lattice parameters) of FeSe depends sensitively on the value of the Hubbard parameter U of on-site correlation and magnetic ordering. Our results suggest that there is a competition between different antiferromagnetic states due to comparable magnetic exchange couplings between first- and second-neighbor Fe sites. As a result, a short range order of stripe antiferromagnetic type is shown to be relevant to the normal state of FeSe at low temperature. We show that there is a strong spin-phonon coupling in FeSe (comparable to its superconducting transition temperature) as reflected in large changes in the frequencies of certain phonons with different magnetic ordering, which is used to explain the observed hardening of a Raman-active phonon at temperatures (similar to 100 K) where magnetic ordering sets in. The symmetry of the stripe antiferromagnetic phase permits an induced stress with orthorhombic symmetry, leading to orthorhombic strain as a secondary order parameter at the temperature of magnetic ordering. The presence of Se vacancies in FeSe gives rise to a large peak in the density of states near the Fermi energy, which could enhance the superconducting transition temperature within the BCS-like picture.

Correlation of Oxygen Storage Capacity and Structural Distortion in Transition-Metal-, Noble-Metal-, and Rare-Earth-Ion-Substituted CeO2 from First Principles Calculation

Oxygen storage/release (OSC) capacity is an important feature common to all three-way catalysts to combat harmful exhaust emissions. To understand the mechanism of improved OSC for doped CeO2, we undertook the structural investigation by X-ray diffraction (XRD), X-ray photoelectron spectroscopy (XPS), H-2-TPR (temperature-programmed hydrogen reduction) and density functional theoretical (DFT) calculations of transition-metal-, noble-metal-, and rare-earth (RE)-ion-substituted ceria. In this report, we present the relationship between the OSC and structural changes induced by the dopant ion in CeO2. Transition metal and noble metal ion substitution in ceria greatly enhances the reducibility of Ce1-xMxO2-delta (M = Mn, Fe, Co, Ni, Cu, Pd, Pt, Ru), whereas rare-earth-ion-substituted Ce(1-x)A(x)O(2-delta) (A = La, Y) have very little effect in improving the OSC. Our simulated optimized structure shows deviation in cation oxygen bond length from ideal bond length of 2.34 angstrom (for CeO2). For example, our theoretical calculation for Ce28Mn4O62 structure shows that Mn-O bonds are in 4 + 2 coordination with average bond lengths of 2.0 and 3.06 angstrom respectively. Although the four short Mn-O bond lengths spans the bond distance region of Mn2O3, the other two Mn-O bonds are moved to longer distances. The dopant transition and noble metal ions also affects Ce coordination shell and results in the formation of longer Ce-O bonds as well. Thus longer cation oxygen bonds for both dopant and host ions results in enhanced synergistic reduction of the solid solution. With Pd ion substitution in Ce1-xMxO2-delta (M = Mn, Fe, Co, Ni, Cu) further enhancement in OSC is observed in H-2-TPR. This effect is reflected in our model calculations by the presence of still longer bonds compared to the model without Pd ion doping. The synergistic effect is therefore due to enhanced reducibility of both dopant and host ion induced due to structural distortion of fluorite lattice in presence of dopant ion. For RE ions (RE = Y, La), our calculations show very little deviation of bonds lengths from ideal fluorite structure. The absence of longer Y-O/La-O and Ce-O bonds make the structure much less susceptible to reduction.

Struggles Over Legitimacy in Global Organizational Restructuring: A Rhetorical Perspective on Legitimation Strategies and Dynamics in a Shutdown Case

Critical organization scholars have focused increasing attention on industrial and organizational restructurings such as shutdown decisions. However, little is known about the rhetorical strategies used to legitimate or resist plant closures in organizational negotiations. In this article, we draw from New Rhetoric to analyze rhetorical struggles, strategies and dynamics in unfolding organizational negotiations. We focus on the shutdown of the bus body unit of the Sweden-based Volvo Bus Corporation in Finland. We distinguish five types of rhetorical legitimation strategies and dynamics. These include the three classical dynamics of logos (rational arguments), pathos (emotional moral arguments), and ethos (authority-based arguments), but also autopoiesis (autopoietic narratives), and cosmos (cosmological constructions). Our analysis contributes to previous studies on organizational restructuring by providing a more nuanced understanding of how contemporary industrial closures are legitimated and resisted in organizational negotiations. This study also increases theoretical understanding of the role of rhetoric in legitimation more generally.

Antarctic ice shelf melting and its impact on the global sea ice-ocean system

Earth s ice shelves are mainly located in Antarctica. They cover about 44% of the Antarctic coastline and are a salient feature of the continent. Antarctic ice shelf melting (AISM) removes heat from and inputs freshwater into the adjacent Southern Ocean. Although playing an important role in the global climate, AISM is one of the most important components currently absent in the IPCC climate model. In this study, AISM is introduced into a global sea ice-ocean climate model ORCA2-LIM, following the approach of Beckmann and Goosse (2003; BG03) for the thermodynamic interaction between the ice shelf and ocean. This forms the model ORCA2-LIM-ISP (ISP: ice shelf parameterization), in which not only all the major Antarctic ice shelves but also a number of minor ice shelves are included. Using these two models, ORCA2-LIM and ORCA2-LIM-ISP, the impact of addition of AISM and increasing AISM have been investigated. Using the ORCA2-LIM model, numerical experiments are performed to investigate the sensitivity of the polar sea ice cover and the Antarctic Circumpolar Current (ACC) transport through Drake Passage (DP) to the variations of three sea ice parameters, namely the thickness of newly formed ice in leads (h0), the compressive strength of ice (P*), and the turning angle in the oceanic boundary layer beneath sea ice (θ). It is found that the magnitudes of h0 and P* have little impact on the seasonal sea ice extent, but lead to large changes in the seasonal sea ice volume. The variation in turning angle has little impact on the sea ice extent and volume in the Arctic but tends to reduce them in the Antarctica when ignored. The magnitude of P* has the least impact on the DP transport, while the other two parameters have much larger influences. Numerical results from ORCA2-LIM and ORCA2-LIM-ISP are analyzed to investigate how the inclusion of AISM affects the representation of the Southern Ocean hydrography. Comparisons with data from the World Ocean Circulation Experiment (WOCE) show that the addition of AISM significantly improves the simulated hydrography. It not only warms and freshens the originally too cold and too saline bottom water (AABW), but also warms and enriches the salinity of the originally too cold and too fresh warm deep water (WDW). Addition of AISM also improves the simulated stratification. The close agreement between the simulation with AISM and the observations suggests that the applied parameterization is an adequate way to include the effect of AISM in a global sea ice-ocean climate model. We also investigate the models capability to represent the sea ice-ocean system in the North Atlantic Ocean and the Arctic regions. Our study shows both models (with and without AISM) can successfully reproduce the main features of the sea ice-ocean system. However, both tend to overestimate the ice flux through the Nares Strait, produce a lower temperature and salinity in the Hudson Bay, Baffin Bay and Davis Strait, and miss the deep convection in the Labrador Sea. These deficiencies are mainly attributed to the artificial enlargement of the Nares Strait in the model. In this study, the impact of increasing AISM on the global sea ice-ocean system is thoroughly investigated. This provides a first idea regarding changes induced by increasing AISM. It is shown that the impact of increasing AISM is global and most significant in the Southern Ocean. There, increasing AISM tends to freshen the surface water, to warm the intermediate and deep waters, and to freshen and warm the bottom water. In addition, increasing AISM also leads to changes in the mixed layer depths (MLD) in the deep convection sites in the Southern Ocean, deepening in the Antarctic continental shelf while shoaling in the ACC region. Furthermore, increasing AISM influences the current system in the Southern Ocean. It tends to weaken the ACC, and strengthen the Antarctic coastal current (ACoC) as well as the Weddell Gyre and the Ross Gyre. In addition to the ocean system, increasing AISM also has a notable impact on the Antarctic sea ice cover. Due to the cooling of seawater, sea ice concentration and thickness generally become higher. In austral winter, noticeable increases in sea ice concentration mainly take place near the ice edge. In regards with sea ice thickness, large increases are mainly found along the coast of the Weddell Sea, the Bellingshausen and Amundsen Seas, and the Ross Sea. The overall thickening of sea ice leads to a larger volume of sea ice in Antarctica. In the North Atlantic, increasing AISM leads to remarkable changes in temperature, salinity and density. The water generally becomes warmer, more saline and denser. The most significant warming occurs in the subsurface layer. In contrast, the maximum salinity increase is found at the surface. In addition, the MLD becomes larger along the Greenland-Scotland-Iceland ridge. Global teleconnections due to AISM are studied. The AISM signal is transported with the surface current: the additional freshwater from AISM tends to enhance the northward spreading of the surface water. As a result, more warm and saline water is transported from the tropical region to the North Atlantic Ocean, resulting in warming and salt enrichment there. It would take about 30 40 years to establish a systematic noticeable change in temperature, salinity and MLD in the North Atlantic Ocean according to this study. The changes in hydrography due to increasing AISM are compared with observations. Consistency suggests that increasing AISM is highly likely a major contributor to the recent observed changes in the Southern Ocean. In addition, the AISM might contribute to the salinity contrast between the North Atlantic and North Pacific, which is important for the global thermohaline circulation.

Global Processes, Local Strategies : migration of educated Dominicans.

The aim of this thesis is to examine migration of educated Dominicans in light of global processes. Current global developments have resulted in increasingly global movements of people, yet people tend to come from certain places in large numbers rather than others. At the same time, international migration is increasingly selective, which shows in the disproportional number of educated migrants. This study discovers individual and societal motivations that explain why young educated Dominicans decide to migrate and return. The theoretical framework of this thesis underlines that migration is a dynamic process rooted in other global developments. Migratory movements should be seen as a result of interacting macro- and microstructures, which are linked by a number of intermediate mechanisms, meso-structures. The way individuals perceive opportunity structures concretises the way global developments mediate to the micro-level. The case of the Dominican Republic shows that there is a diversity of local responses to the world system, as Dominicans have produced their own unique historical responses to global changes. The thesis explains that Dominican migration is importantly conditioned by socioeconomic and educational background. Migration is more accessible for the educated middle class, because of the availability of better resources. Educated migrants also seem less likely to rely on networks to organize their migrations. The role of networks in migration differs by socioeconomic background on the one hand, and by the specific connections each individual has to current and previous migrants on the other hand. The personal and cultural values of the migrant are also pivotal. The central argument of this thesis is that a veritable culture of migration has evolved in the Dominican Republic. The actual economic, political and social circumstances have led many Dominicans to believe that there are better opportunities elsewhere. The globalisation of certain expectations on the one hand, and the development of the specifically Dominican feeling of ‘externalism’ on the other, have for their part given rise to the Dominican culture of migration. The study also suggests that the current Dominican development model encourages migration. Besides global structures, local structures are found to ve pivotal in determining how global processes are materialised in a specific place. The research for this thesis was conducted by using qualitative methodology. The focus of this thesis was on thematic interviews that reveal the subject’s point of view and give a fuller understanding of migration and mobility of the educated. The data was mainly collected during a field research phase in Santo Domingo, the Dominican Republic in December 2009 and January 2010. The principal material consists of ten thematic interviews held with educated Dominican current or former migrants. Four expert interviews, relevant empirical data, theoretical literature and newspaper articles were also comprehensively used.

Effects of Zr and Ti doping on the dielectric response of CeO2: A comparative first-principles study

Zr doping in ceria (CeO2) results in enhanced static dielectric response compared to pure ceria. On the other hand, Ti doping in ceria keeps its dielectric constant unchanged. We use first-principles density functional theory calculations based on pseudopotentials and a plane wave basis to determine electronic properties and dielectric response of Zr/Ti-doped and oxygen-vacancy-introduced ceria. Softening of phonon modes is responsible for the enhancement in dielectric response of Zr-doped ceria compared to that of pure ceria. The ceria-zirconia mixed oxides should have potential use as high-k materials in the semiconductor industry. (c) 2010 Elsevier Ltd. All rights reserved.

On the principles of Japanese forest policy since 1950.

XVIII IUFRO World Congress, Ljubljana 1986.

Ühis- ning erijooni Soome ja Eesti Kaitseliitude puhkpillimuusikas aastatel 1925-1934.

This study focuses on the similarities and differences between the Estonian Defence League and the Finnish Civil Guard brass bands during the period 1925-1934. By 1934 this paramilitary volunteer state defence organisation had reached stability in its development, such that social, cultural and patriotic education of the people - with the help of brass band music among other means- had acquired a significant role, in addition to prioritised military and sports activities. The study begins with introductory paragraphs I and II, which describes the founding of the organisations, their participation in the Wars of Independence and their subsequent peace time activities as well as their representation in the media at the time. The thesis also briefly introduces military music in Finland and Estonia, as well as describes the influence of military music on the Defence League brass bands. The period under review includes the global economic crisis, which undoubtedly concerned the Defence League/Civil Guard and the Lapua and War of Independence movements, which greatly affected the apolitical principles of the organisations. The main emphasis of the thesis is the Defence League/Civil Guard brass band`s musical activities in two counties - Etelä-Pohjanmaa and Pärnumaa, while also including a general overview of the Estonian Defence League brass bands´ activities. One of the most important benefits of the thesis is its introduction of the brass band repertoire in use at the time, which was played by both professional and as well as amateur orchestras the latter of which also included the brass bands of the Defence League/Civil Guard and the Fire Services. Brass band music held a secondary, yet significant position in the Defence League/Civil Guard, where the orchestra as a musical grouping was obliged to perform not only at inner-organisational and national celebrations but also at any event requiring brass band music, such as song festivals, singing days, and other local cultural events. The professional preparation of the band conductors at the beginning of the period under review was not well specialised, but the training of the Defence League/Civil Defence brass band conductors was carried out regularly in both republic according to the opportunities and dedicated training programmes available. The musicians of the Defence League/Civil Defence brass bands were at the same time members of the military organisations as well as amateur musicians, which placed upon them extra demands - they were under close public supervision in all situations. Based on the principle of chronology it appeared that both Finnish and Estonian respective organisations´ brass bands used the gradually improving economic situation for purchasing musical instruments, obtaining repertoire and training musicians/conductors. Despite the fact that brass band music in the Defence League/Civil Guard was considered an amateur activity and a hobby, the more far-reaching objective of the organisation was to resemble the Defence Forces´ orchestras as closely as possible in all aspects. The Defence League/Civil Guard brass band music definitely had a significant influence on forming, developing and enriching music life in both republics. The reviewed nine-year period introduced the musical activities of the Defence League/Civil Guard against the background of the everyday life of the organisation and the need for brass band music and its continuity in the voluntary state defence institutions of both republics.

The Liminality of Cosmopolitanism from the Global South – Transnational Activists in Costa Rica between Colonial Legacies and Decolonial Attitudes

The Master’s thesis examines whether and how decolonial cosmopolitanism is empirically traceable in the attitudes and practices of Costa Rican activists working in transnational advocacy organizations. Decolonial cosmopolitanism is defined as a form of cosmopolitanism from below that aims to propose ways of imagining – and putting into practice – a truly globe-encompassing civic community not based on relations of domination but on horizontal dialogue. This concept has been developed by and shares its basic presumptions with the theory on coloniality that the modernity/coloniality/decoloniality research group is putting forward. It is analyzed whether and how the workings of coloniality as underlying ontological assumption of decolonial cosmopolitanism and broadly subsumable under the three logics of race, capitalism, and knowledge, are traceable in intermediate postcolonial transnational advocacy in Costa Rica. The method of analysis chosen to approach these questions is content analysis, which is used for the analysis of qualitative semi-structured in-depth interviews with Costa Rican activists working in advocacy organizations with transnational ties. Costa Rica was chosen as it – while unquestionably a Latin American postcolonial country and thus within the geo-political context in which the concept was developed – introduces a complex setting of socio-cultural and political factors that put the explanatory potential of the concept to the test. The research group applies the term ‘coloniality’ to describe how the social, political, economic, and epistemic relations developed during the colonization of the Americas order global relations and sustain Western domination still today through what is called the logic of coloniality. It also takes these processes as point of departure for imagining how counter-hegemonic contestations can be achieved through the linking of local struggles to a global community that is based on pluriversality. The issues that have been chosen as most relevant expressions of the logic of coloniality in the context of Costa Rican transnational advocacy and that are thus empirically scrutinized are national identity as ‘white’ exceptional nation with gender equality (racism), the neoliberalization of advocacy in the Global South (capitalism), and finally Eurocentrism, but also transnational civil society networks as first step in decolonizing civic activism (epistemic domination). The findings of this thesis show that the various ways in which activists adopt practices and outlooks stemming from the center in order to empower themselves and their constituencies, but also how their particular geo-political position affects their work, cannot be reduced to one single logic of coloniality. Nonetheless, the aspects of race, gender, capitalism and epistemic hegemony do undeniably affect activist cosmopolitan attitudes and transnational practices. While the premisses on which the concept of decolonial cosmopolitanism is based suffer from some analytical drawbacks, its importance is seen in its ability to take as point of departure the concrete spaces in which situated social relations develop. It thus allows for perceiving the increasing interconnectedness between different levels of social and political organizing as contributing to cosmopolitan visions combining local situatedness with global community as normative horizon that have not only influenced academic debate, but also political projects.

Geometric Characterizations for Patterson-Sullivan Measures of Geometrically Finite Kleinian Groups

The main results of this thesis show that a Patterson-Sullivan measure of a non-elementary geometrically finite Kleinian group can always be characterized using geometric covering and packing constructions. This means that if the standard covering and packing constructions are modified in a suitable way, one can use either one of them to construct a geometric measure which is identical to the Patterson-Sullivan measure. The main results generalize and modify results of D. Sullivan which show that one can sometimes use the standard covering construction to construct a suitable geometric measure and sometimes the standard packing construction. Sullivan has shown also that neither or both of the standard constructions can be used to construct the geometric measure in some situations. The main modifications of the standard constructions are based on certain geometric properties of limit sets of Kleinian groups studied first by P. Tukia. These geometric properties describe how closely the limit set of a given Kleinian group resembles euclidean planes or spheres of varying dimension on small scales. The main idea is to express these geometric properties in a quantitative form which can be incorporated into the gauge functions used in the modified covering and packing constructions. Certain estimation results for general conformal measures of Kleinian groups play a crucial role in the proofs of the main results. These estimation results are generalizations and modifications of similar results considered, among others, by B. Stratmann, D. Sullivan, P. Tukia and S. Velani. The modified constructions are in general defined without reference to Kleinian groups, so they or their variants may prove useful in some other contexts in addition to that of Kleinian groups.

Global, Local and Vector-Valued Tb Theorems on Non-Homogeneous Metric Spaces

Various Tb theorems play a key role in the modern harmonic analysis. They provide characterizations for the boundedness of Calderón-Zygmund type singular integral operators. The general philosophy is that to conclude the boundedness of an operator T on some function space, one needs only to test it on some suitable function b. The main object of this dissertation is to prove very general Tb theorems. The dissertation consists of four research articles and an introductory part. The framework is general with respect to the domain (a metric space), the measure (an upper doubling measure) and the range (a UMD Banach space). Moreover, the used testing conditions are weak. In the first article a (global) Tb theorem on non-homogeneous metric spaces is proved. One of the main technical components is the construction of a randomization procedure for the metric dyadic cubes. The difficulty lies in the fact that metric spaces do not, in general, have a translation group. Also, the measures considered are more general than in the existing literature. This generality is genuinely important for some applications, including the result of Volberg and Wick concerning the characterization of measures for which the analytic Besov-Sobolev space embeds continuously into the space of square integrable functions. In the second article a vector-valued extension of the main result of the first article is considered. This theorem is a new contribution to the vector-valued literature, since previously such general domains and measures were not allowed. The third article deals with local Tb theorems both in the homogeneous and non-homogeneous situations. A modified version of the general non-homogeneous proof technique of Nazarov, Treil and Volberg is extended to cover the case of upper doubling measures. This technique is also used in the homogeneous setting to prove local Tb theorems with weak testing conditions introduced by Auscher, Hofmann, Muscalu, Tao and Thiele. This gives a completely new and direct proof of such results utilizing the full force of non-homogeneous analysis. The final article has to do with sharp weighted theory for maximal truncations of Calderón-Zygmund operators. This includes a reduction to certain Sawyer-type testing conditions, which are in the spirit of Tb theorems and thus of the dissertation. The article extends the sharp bounds previously known only for untruncated operators, and also proves sharp weak type results, which are new even for untruncated operators. New techniques are introduced to overcome the difficulties introduced by the non-linearity of maximal truncations.

Size dependence of the bulk modulus of semiconductor nanocrystals from first-principles calculations

The variation in the bulk modulus of semiconductor nanoparticles has been studied within first-principles electronic-structure calculations using the local density approximation (LDA) for the exchange correlation. Quantum Monte Carlo calculations carried out for a silicon nanocrystal Si87H76 provided reasonable agreement with the LDA results. An enhancement was observed in the bulk modulus as the size of the nanoparticle was decreased, with modest enhancements being predicted for the largest nanoparticles studied here, a size just accessible in experiments. To access larger sizes, we fit our calculated bulk moduli to the same empirical law for all materials, the asymptote of which is the bulk value of the modulus. This was found to be within 2-10% of the independently calculated value. The origin of the enhancement has been discussed in terms of Cohen's empirical law M.L. Cohen, Phys. Rev. B 32, 7988 (1985)] as well as other possible scenarios.

The top 100 questions of importance to the future of global agriculture

Despite a significant growth in food production over the past half-century, one of the most important challenges facing society today is how to feed an expected population of some nine billion by the middle of the 20th century. To meet the expected demand for food without significant increases in prices, it has been estimated that we need to produce 70-100 per cent more food, in light of the growing impacts of climate change, concerns over energy security, regional dietary shifts and the Millennium Development target of halving world poverty and hunger by 2015. The goal for the agricultural sector is no longer simply to maximize productivity, but to optimize across a far more complex landscape of production, rural development, environmental, social justice and food consumption outcomes. However, there remain significant challenges to developing national and international policies that support the wide emergence of more sustainable forms of land use and efficient agricultural production. The lack of information flow between scientists, practitioners and policy makers is known to exacerbate the difficulties, despite increased emphasis upon evidence-based policy. In this paper, we seek to improve dialogue and understanding between agricultural research and policy by identifying the 100 most important questions for global agriculture. These have been compiled using a horizon-scanning approach with leading experts and representatives of major agricultural organizations worldwide. The aim is to use sound scientific evidence to inform decision making and guide policy makers in the future direction of agricultural research priorities and policy support. If addressed, we anticipate that these questions will have a significant impact on global agricultural practices worldwide, while improving the synergy between agricultural policy, practice and research. This research forms part of the UK Government's Foresight Global Food and Farming Futures project.

Challenges in Earth sciences: the 21st century

The 21st century poses many challenges for global sustainability. Among them, most importantly, the human race will encounter scarcity of raw materials and conventional energy resources. And, India may have to take the brunt of these problems as it is going to be the most populated region of the world with concomitant increase in energy demand and requirement of other resources. India will be the testing ground for introducing newer ways of green technology and innovative principles of resource management and utilization. With the vagaries of potential climate change gathering clouds in the background, Earth sciences will have a special and predominant role in guiding the society in prioritizing our resource discovery, utilization and their consumption and the upkeep of environment. On the fundamental level, Earth sciences are going through a most exciting phase of development as a born-again science. Technological breakthroughs including the satellite-based observations augur well for gaining new insights into Earth processes. A set of exciting fundamental problems that are globally identified will set the stage for an exhilarating period of new discoveries. Improvements in numerical and computer-based techniques will assist in modelling of Earth processes to unprecedented levels. India will have to take special effort in improving the existing experimentation facilities in the Earth science departments of the country, and also the general level of Earth science education to meet the global standards. This article presents an Earth science vision for the 21st century in an Indian context.

Global solutions describing the collapse of a spherical or cylindrical cavity

The collapse of a spherical (cylindrical) cavity in air is studied analytically. The global solution for the entire domain between the sound front, separating the undisturbed and the disturbed gas, and the vacuum front is constructed in the form of infinite series in time with coefficients depending on an ldquoappropriaterdquo similarity variable. At timet=0+, the exact planar solution for a uniformly moving cavity is assumed to hold. The global analytic solution of this initial boundary value problem is found until the collapse time (=(gamma–1)/2) for gamma le 1+(2/(1+v)), wherev=1 for cylindrical geometry, andv=2 for spherical geometry. For higher values of gamma, the solution series diverge at timet — 2(beta–1)/ (v(1+beta)+(1–beta)2) where beta=2/(gamma–1). A close agreement is found in the prediction of qualitative features of analytic solution and numerical results of Thomaset al. [1].