ISSN 1880-8468
Technical Report ofInternational Development Engineering࿖㓙㐿⊒Ꮏቇႎ๔TRIDE-2007- 01

February 1, 2007

Abstracts of Master Theses
Presented in February 2007

Department of International Development Engineering,
Graduate School of Science and Engineering,Tokyo Institute of Technology
http://www.ide.titech.ac.jp/TR

Technical Report of International Development Engineering
TRIDE-2007-01
Table of Contents
Effect of preparation condition applied by Chevrel phase sulfides synthesis on
the hydrodesulfurization activity of new developing catalysts. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Masakazu KURATA

1

Introduction of ICT (information and communication technology) to the world
heritage site in developing countries – Case study of Luang Prabang
Lao PDR –. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Nobuo OOKA

5

Study of text classification with statistical pattern recognition method. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Yusuke WATADA

9Extraction of liquid organic multi-component mixture with ionic liquids assolvents. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Kenya AKAISHI

13

Moving picture coding with wavelet transform. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Takashi OHTA

17

Creation and evaluation of e-Learning materials for introduction of engineering. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Masanori KANEKO

21

Study on applicability of steel slag hydrated matrix to steel reinforced members
under marine environment. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Tomonari KIMURA

25

A study on gas-liquid two-phase flow patterns in rectangular micro-channels. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .Shintaro KOBAYASHI

29

In-situ soil remediation by sedimentation method – Heavy metal ions’ adsorption characteristics on soil and adsorbents –. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Ryoichi SUGA

33

The psychological effects of behavioral change caused by public works on place
attachment. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Haruna SUZUKI

i37

A procedure of determining parameters by adapting to the results calculated
by MEAM’ 92. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .Kento TOKUMARU

41

Long term prediction of strength deterioration due to Ca leaching from cement
treated soil
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Katsufumi HASHIMOTO

45

Analysis of turbulent organized structures at the urban outdoor scalemodel(COSMO). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Satoshi HIROOKA

49

Efficiency of the implicit finite element method of geotechnical engineering. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Takahiro FURUMURA

53

Synthesis of zeolites from lake sludge. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Wei QU

57

Influence of creep phenomenon on manipulation by shape control of rubberelastic membrane. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Xuanxuan CHANG

61

Development of a stored channel model for UWB link level system design andevaluation. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .Huynh Thi thanh TRIEU

65

Numerical simulation on temperature field for Newtonian fluid enclosed inconcentric sphere. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Bian HU

69

A Study on corrosion of paint-coated steel with defects in marine environment. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Aung Kyaw MIN

73

Study on antennas embedded in mobile phone housing. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Jeonghoon HAN

77ii

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

Effect of preparation condition applied by Chevrel phase sulfides synthesis on
the hydrodesulfurization activity of new developing catalysts
Student number:04M18104 Name: Masakazu KURATA Supervisor: Hirofumi HINODE

シェブレル相硫化物合成法を応用した新規水素化脱硫触媒の研究

倉田 雅一
本研究ではシェブレル相硫化物合成法を応用して調製したサンプルの水素化脱硫活性を評価した。750ºC で水素化処理を行ったサンプルは水素化時間の増加と共に活性が減少したが、450ºC
で水素化処理を行ったサンプルは水素化時間の増加と共に活性が増大し、水素化 24 時間以上では市販の工業触媒を上回る活性を示した。XRD、BET、TEM の結果から水素化時間 750ºC のサンプ
ルでは熱劣化及び結晶化が起こり、水素化時間 450ºC のサンプルでは活性相である CoMoS 相の量が促進されたことが確認できた。

1.

the hydrogenation of supported CPS needed to bedone at high temperature for a few days.
In this research, the effect of preparation
conditions, in particular, hydrogenating time andtemperature, on HDS activity and thecharacterization of supported Chevrel phasesulfides have been investigated.

Introduction

Sulfur removal in gasoline and diesel oil is
now strongly desirable because sulfur may
adversely affect catalytic performance for exhaustemission treatment. So worldwide refiners aretrying to provide clean transportation fuels tomeet increasingly demanding environmental
regulations. Hydrodesulfurization (HDS) is sulfur
compounds removal process of petroleum
feedstocks in petroleum refiners. Now numerous
attempts have already been made to develop more
active HDS catalyst.
In the last half centuries, alumina supported
Co-Mo, Ni-Mo, Ni-W sulfides catalysts were
widely used as HDS catalyst. Active site ofconventional CoMo sulfides catalyst is theso-called CoMoS phase, in which Co is locatedon the edges of highly dispersed MoS 2 slabs1).
MoS2 edge plane has two planes, rim plane andedge plane. Rim plane has both HDS andhydrogenation (HYD) ability, while edge plane
has only HDS ability2). There are two types of
CoMoS phase, Type I and Type II. Type I isrelated to highly dispersed single MoS2 slab
maintaining their interactions with the support
thus is less activity. Type II is related to lessdispersed multi MoS2 slabs with less interaction
with support, resulting in high activity.
In this research, samples were prepared using
methods applied by Chevrel phase sulfides (CPS)synthesis3). CPS are ternary molybdenum sulfides
and have the general formula MxMo6S8-y, where
M represents ternary metal, one of about forty
metal4). Previous study showed supported CPS
had a high activity5,6), and we presumed that CPS
and M-Mo-S phase of active site for commercial
catalysis is similar structure. On the other hand,

Figure 1 Structure of Chevrel phase sulfides

2.

Experimental

Supported CoCPS catalysts were prepared byimpregnation method. Precursors of the catalyst
were Co(NO3)2·6H2O and (NH4)6Mo7O24·4H2O(Wako Pure Chemical Industries, Ltd.).
JRC-ALO2 (Al2O3, Reference Catalyst, Catalysis
Society of Japan) was used as a support. The Mo
loading level was fixed at 15wt%. Co/Mo ratio
was 0.28. In the first step, JRC-ALO2 was
impregnated with precursors aqueous solution for
24h. Then, the impregnated sample was dried up,
and calcined at 500ºC for 4h. The sample was1Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

pelletized, crushed, and sieved to 0.71~1.00mm
in order to reduce pressure drop In the second
step, sulfurization was carried out at 400ºC for 2hunder 10% H2S/H2 flow in order to prepare
sulfides. The hydrogenation was then carried out
under H2 flow at 450ºC or 750ºC for 0-72h.
Commercial catalyst (Nippon Ketjen Co., Ltd.,KF 752-1.5E) was pelletized, crushed, and sieved
to 0.71~1.00mm, and then sulfurized at 400ºC for
2h under 10% H2S/H2 flow.
HDS reaction was carried out in a fixed-bed
flow reactor under atmospheric pressure4). 0.2g ofthe sample was packed in the reactor and the total
flow rate of the feed mixture gas was fixed at
20ml/min. Thiophene was used as a target
material of sulfur compounds. The concentrationof thiophene in the feed stream was maintained at
2.1mol%. The reaction temperature was selected
at 300 and 350ºC. Thiophene was quantitatively
analyzed by GC-MS (Shimadzu Corp.,
QP5050A), and HDS activity and selectivity of
products were calculated.
The samples were characterized by XRD(Rigaku Co., MutiFlex), BET specific surface
area measurement (COULTER Co., SA3100) andTEM (JEOL Ltd., JEM-2010F).

3.

Figure 3 shows thiophene conversion of
CoCPS/ALO2 hydrogenating at 450ºC withdifferent
hydrogenatingtime.
Thiophene
conversion of the samples increased with anincrease in hydrogenating time, and sample ofhydrogenated for 48h showed thiophene
conversion at 300 and 350ºC of 96.4% and 99.4%,
respectively. There was, however, a slightdecrease in the thiophene conversion of samplehydrogenated for 72h compared to that of 48h.

Thiophene conversion [%]

100

80
70
60

300ºC
350ºC

50
0

24

48

Figure3 Thiophene conversion of CoCPS/ALO2hydrogenated at 450ºC with differenthydrogenating time

Figure 2 shows thiophene conversion of
CoCPS/ALO2 hydrogenating at 750ºC withdifferent
hydrogenating time.
Thiophene
conversions of samples decreased with anincrease in hydrogenating time. Thiophene
conversions at 300ºC and 350ºC of nothydrogenated sample were 79.63% and 96.34%,
respectively, and those of hydrogenated for 24h
were 53.05% and 76.82%, respectively.

Figure 4 shows thiophene conversion ofcommercial CoMoS catalyst samples withdifferent hydrogenation temperature and time. H0means not hydrogenated sample, and H450 and
H750 mean hydrogenation for 24h at 450 and
750ºC, respectively. Thiophene conversions ofH450 sample were higher than those ofcommercial CoMoS catalyst.
100

100

72

Hydrogenating time [h]Result and Discussion

300ºC
350ºC

Thiophene conversion [%]Thiophene conversion [%]

90

90
80
70
60

300ºC
350ºC

90
80
70
60
50

50
0

8
16
Hydrogenating time [h]

comm.

24

H0

H450

H750

Figure 4 Thiophene conversions of commercial
CoMoS catalyst and H0, H450, and H750

Figure2 Thiophene conversion of CoCPS/ALO2hydrogenated at 750ºC with differenthydrogenating time2

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

Figure 5 shows XRD patterns of H0, H450 and
H750 samples. XRD pattern of H750 sample
showed clearly peaks of MoS2 at 2θ =14.38º(002),
32.68º(100),
33.51º(101),
58.33º(110) or/and CoCPS at 2θ = 13.78(101),
33.58(122), and the peaks of JRC-ALO2 (2θ =38º, 46º, and 67º). According to XRD patterns ofH450 and H0 samples showed small peaks at 2θ
= 33º and 58º. From the XRD pattern, MoS2
or/and CoCPS well crystallized in H750 sample.
Figure 6 shows TEM image of H750 sample. 2
or 3 slabs and large MoS2 crystals were observedin this image.

of other samples. Due to hydrogenating at hightemperature, specific surface area of H750
sample considerably decreased.
Table 1 Specific surface area of H0, H450 and H750

SampleSpecific surface areaH0

202.18m2/g

H450

197.72m2/g

H750

154.31m2/g

JRC-ALO2

280.10m2/g

For these reason, low activity of H750 sample
was attributed to increased crystallinity anddecreased specific surface area.

Figure 5 XRD patterns of H0, H450, and H750 samples
and JRC-ALO2

Figure 7 TEM image of H0 sample

Figure 6 TEM image of H750 sample

Table 1 shows specific surface area of H0, H450
and H750 samples. Specific surface areas of H0and H450 were almost the same, while that of
H750 sample was significantly smaller than those

Figure 8 TEM image of H450 sample

3Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

Differences between H0 and H450 samples
were not identified by XRD patterns and specificsurface areas. On the other hand, TEM images ofboth samples, H0 in Figure 7 and H450 in Figure
8 show remarkable difference. H0 sample in
Figure 7 had only 2-4 slabs, so CoMoS phase of
H0 sample was mainly Type I and had lessactivity. As shown in Figure 8, H450 sample
had 3-8 slabs, so CoMoS phase of H450 sample
was mainly Type II and had high activity. For thisreason, hydrogenating at 450ºC promoted theformation of CoMoS Type II phase.

4.

Conclusion

Thiophene conversion of hydrogenated at
750ºC samples decreased with an increase inhydrogenating time, due to increased crystallinityand decreased specific surface area of the sample.
On the other hand, thiophene conversion ofhydrogenated at 450ºC samples increased with anincrease in hydrogenating time, due to the formationof CoMoS Type II. Hydrogenated at 450ºC samples
had higher HDS activity than commercial CoMoS
catalyst. H450 sample has a high HYD activity
because of the formation of rim plane mainly.

Reference[1] H.Topsøe et al., Bull.Soc.Chim.Belg., 90, 1189
(1981)[2] M.Daage et al., J.Catal., 149, 414 (1994)[3] R.Chevrel et al., J.Solid.State Chem., 3, 519(1971)[4] G.J.Dudly et al., J.Solid.State Chem., 32, 259(1980)[5] V.Harel-Michaud et al., J.Alloys.Comp., 317-318,
195 (2001)[6] S.Ooi et al., React.Kinet.Catal.Lett., 82, 89(2004)[7] A.Nishijima et al., Sekiyugakkaishi, 32, 35(1989)[8] H.Zhang. Master’s thesis, Tokyo Institute oftechnology (2006)
Figure 9 TEM image of commercial CoMoS catalyst8)

In order to measure the HYD activity of thesamples, selectivity of the products was measured.
1-butene, trans-2-butene, cis-2-butene andn-butane were detected. Table 2 shows selectivity
for H0 and H450 samples. H450 sample washigher HDS activity than H0 sample, and also
H450 sample was higher HYD activity than H0
sample. These can be explained according to theparticle size of MoS2 slabs, since MoS2 particle
size is correlated with the amount of rim and edge
planes7). H450 sample in Figure 8 showed small
particle size of MoS2 slabs compared tocommercial CoMoS catalyst shown in Figure 98)
indicating formation of rim plane mainly.
Therefore, H450 sample had a high HYD activity.
Table 2 Products selectivities of H0 and H450 samples

SampleButaneButenesH071.80%

28.20%

H450

75.36%

24.64%

4

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

Introduction of ICT (Information and Communication Technology) to the World Heritage
Site in developing countries - Case study of Luang Prabang, Lao PDR Student Number: 04M51248

Name: Nobuo OOKA Supervisors: Jun-ichi TAKADA, Shinobu YAMAGUCHI

開発途上国の世界文化遺産地域における情報通信の導入-ラオス、ルアンパバーンにおけるケーススタディ大岡信夫
本論文では、ルアンパバーンにおいて現地政府・UNESCO との協働のもと行ってきたニーズアセ
スメントおよび、それをもとに実施したデータベース開発・現地人材育成・ICT センターの設置運営などの活動ついて報告する。また、実施を通して実現した、開発途上国の世界文化遺産を核にした持続可能な地域開発への ICT 分野での貢献を評価するとともに、その経験から持続可能性に対する留意事項を抽出した。

field, and chapter 6 touches upon the further issues. Finally,
chapter 7 discusses on the sustainability of the development
activities through the experiences of the involvement in thefield activities.

1 Introduction
Many development projects have been implemented in thedeveloping countries, but in 1980s, many assistance projects
faced the difficulties on the sustainability of theirimplementation. Then, sustainable development has been
receiving a lot of attention. In other word, experts, government,
and international organizations are very much concerned how
development activities can be continued as a locally owned
implementation. Through the experiences on the actually
implemented activities in Luang Prabang, Lao PDR, this thesispresents the case study as one of the applicable example tocontribute to the sustainable development in developing
countries.

2 Background2.1

Regional background

Lao PDR is a landlocked country in the Indochinese Peninsula
and regarded as one of the underdeveloped countries with low
GDP per capita (PPP) of $1,969US1 and low HDI (Human
Development Indicator) of 0.5532 which is 133rd out of 177
countries. Luang Prabang is the provincial capital of Luang
Prabang province which is located northern Lao PDR. In 1995,
the town of Luang Prabang was inscribed into UNESCOWorld Heritage site, due to its well-preserved townscape andnatural environment. After the inscription, the number oftourists has been remarkably increasing; an average of
50,000–60,000 visitors in 1995 and 1996 per year has
increased to 196,106 3 in 2003. Since the town of Luang
Prabang is the second largest city in Lao PDR and is facing
rapidly growing tourism, regional development including
surrounding area, putting the heritage at the core, has been
expected. Generally, it is often the case that the World Heritage
inscription brings not only the positive impacts, such asincreased income from tourist, more employment
opportunities, improved infrastructure, promotion ofcommunity identity, and international cooperation and subsidy,
but also negative impacts such as rapid and unplanned
construction, environmental destruction, and increase ofincome gap. Therefore, in order to maximize the positive
impacts with reducing negative impacts, it has been concernedto promote comprehensive approach integrating the various
development components, such as economic development,

The town of Luang Prabang, inscribed into the UNESCOWorld Heritage Site in 1995, has experienced regional
economic growth due to the increasing tourism. In order topromote sustainable regional development, it has been
discussed how different sectors can be well-integrated for theregional development in Luang Prabang, such astransportation and infrastructure and communication
technologies, education and health, commerce, agriculture andforestry and so on.
In the meanwhile, ICT development in developing countries
has been considered as one of the important factors in theregional development since 1980s [1]. Based on the requestfrom Lao government, since 2003 Tokyo Tech team has been
taking a role on technical assistance to identify and provide
appropriate and applicable ICT application in Luang Prabang.
While the author has been involved in the needs assessment
and its implementation at the local site, it has been focused
how ICT can contribute to the heritage conservation anddevelopment.
In this thesis, chapter 2 covers the background information,and chapter 3 describes the current difficulties and needs in
ICT area identified through the field survey. Then chapter 4
highlights actually implemented activities. According to theimplementation, chapter 5 shows the observed impacts in the1

Source: CIA, “The World Factbook – Laos-”,

https://www.cia.gov/cia/publications/factbook/geos/la.html
2

Source: UNDP, “Human Development Report 2006”

Sources: ADUC, SCOT, November 2004; Luang Prabang,
Provincial Tourism Office, 2004

3

5

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

environmental protection, human development, social
development and so on.
2.2

participation of local communities have been regarded as animportant factor to promote the sustainable development. Also,
having no official web site, Luang Prabang would have missed
the potential interest of the global public including tourists,
researchers, fund agency, business people and so on. Therefore,
effective means to disseminate the public information to thevarious target groups were expected.

Collaboration framework on the activities

In order to promote comprehensive approach in Luang
Prabang, the Inter-institutional committee, so-called Local
Heritage Committee (LHC) is organized. La Maison duPatrimoine (MdP) is one of the focal institutions in LHC,
created in 1996 according to the UNESCO recommendation,
mainly taking the role on rehabilitation and management of thecultural heritage property and natural heritage, development ofmarket-place, and coordination of international projects.
Meanwhile, the tripartite collaboration among Laogovernment, Tokyo Tech, and UNESCO has been held on
ICT development area since late 2003, when Tokyo Techparticipated in the inter-government symposium in Luang
Prabang according to the request of UNESCO. The role ofTokyo Tech is to identify and provide the appropriate andlocally applicable ICT application to contribute to thesustainable regional development in Luang Prabang.
Especially in the field activities, Tokyo Tech has been incooperation with MdP.

Active reuse of restored architecture: There have been anumber of restored heritage architectures in Luang Prabang,but some of them have been under utilized. Not only
conservation but also innovative reuse of heritage architectures
was vital to promote benefits for the local communities andpeople.
3.2

Identified component for ICT development
According to the needs mentioned above, the following fivecomponents for ICT development shown in Figure 1 were
proposed.

1. Development of database application: With
well-structured database of digitalized data to enable long-term
storage of data, and effective interfaces to fully utilize storeddata.

3 Needs analysis
3.1

2. Establishment of network infrastructure: To encourage theinformation sharing among institutions.

Needs analysis

3. Creation of Web sites: To give the effective means todisseminate information to global public.

Based on the request from the local government, Tokyo Techteam conducted needs assessment in Luang Prabang in 2004.
Through the field survey, the following five needs for ICT
development were identified.

4. Development and management of ICT centers: To
provide information to the public in Luang Prabang, and topromote reuse of the restored architectures.

Long-term storage of data: There was abundant information
accumulated through the heritage conservation anddevelopment, which was not only given the worth as heritage
data, but also regarded as the valuable resource of theexperiences in heritage management and urban development.
However, most of information was not kept in good condition,
being stored in paper-based. The appropriate condition to
enable long-term storage was needed as the fundamental for
effective information management.

5. Human resource development: Human resource is one of
the essential factors to realize sustainable development [2].
Since there were no ICT professionals in Luang Prabang,
continuous human resource development was required to
implement and sustain the activities.

Utilization of stored information: The stored paper-based
data has had disadvantages in searching, updating, and sharing
data. Recently, the effort to integrate DBMS (Database
Management System) has been made by the local stuffs, but it
was not so actively used as they expected because of lack ofstrategy and human resource. The effective tools which could
enhance the utilization of stored data were needed.
Information sharing: In Luang Prabang, network
infrastructure and information sharing framework wereinsufficient. In order to promote the comprehensive approach,the means to encourage intercommunication among
institutions were required.

Figure 1 Components of ICT development [4]4 Implementation of the pilot activities
In the implementation phase, the author has been mainly
involved in local human resource development, development
of database application, and development of ICT center out of
the five components mentioned above.

Information provision to the public: There was no facility
where the population and visitors in Luang Prabang can accessto the reliable information of the town, although awareness and

4.1
6

Development of prototype database

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

to public in the first pilot ICT center established in July 2006.
For the sustainable operation, self-sustainable financial cycle
has been required by reducing the initial/running cost andpromoting income generating activities. Also, involvement of
the local communities was appreciated such as employment
the local staff, encourage the communities’ participation in thediscussion, and the promotion of the business by local people
in the site.

application
Since the development of database application was the firsttrial in Luang Prabang, the prototype database application,so-called “HeritageDB”, was created in Jan 2005, in order toshow the example and clarify the requirements at the local site.
The prototype contains data on the heritage architectures,
wetlands, roads, drawings, and maps. Some data had beenpreviously managed using Microsoft Access, but there weresome problems, such as lack of regularization in data structure,
buggy macros left due to the leave of the developer, anddifficulty in sharing via the network. Therefore, the prototype
was implemented to solve the problems of the existing
database with some extension.

The center was established in the site of the representative
heritage of Lao wooden architecture. It has been equipped withthe five PCs to access to those two prototype databases andweb site of MdP, browsing space of the books, CDROM,
brochures, posters and goods of Luang Prabang, of which parts
are also for sale in order to earn the operation cost of the ICT
center.

The prototype system was composed with Debian
GNU/Linux as OS (Operating System), PostgreSQL as
DBMS (Database Management System) and Zope as Web
application server, having the following three features;web-based interface to promote sharing and dissemination ofdata, use of FOSS (Free and Open Source Software) to keep
the legality and cost efficiency, and multi-language support ofthe interface for the local population and various target groups
who use different languages.
4.2

In the first sixteen weeks, 1,500 people have visited the center,of which 73% were tourists from foreign countries. According
to the feedbacks from the visitors, most of visitors regard the
database application as the attracting resource of information,and they were satisfied with the well-equipped center. On theother hand, also quite some comments for the required
improvements were given. For the technical aspects, theenhancement of multi-language support, the moreuser-friendly graphical interface, and accessibility to the
database application from Internet were suggested. Also for thecontents, more information about the Luang Prabang history,
and frequent update were requested.

Local human development

According to the requirement for local human resource, thelocal ICT team was organized with six Lao staffs from thegovernment institutions, who have been equipped with insighton heritage conservation and urban development owing to theprevious works.

5 Impacts from the activities

The training workshops have been conducted with thefollowing two features; one is needs based practical training
which provides not only ICT skills but also the productions
which are actually desired from the institutions. The other isthe usage of FOSS in order to promote understanding on legal
use of software and to save the cost.

Through the implementation, the following impacts wereobserved at the local site.

Further development of database application

1. Trained local ICT team members have become
well-recognized ICT engineers in Luang Prabang. For instance,
the governor’s office requests some technical support by the
ICT team and it has increased motivation for the local ICTteam to contribute to Luang Prabang

After the training was started, further development of databaseapplication has been conducted with the local ICT team basedon the skills learned.

2. Introduction and training of FOSS raised awareness on legal
use of the software, giving the options besides illegal copy ofthe commercial software.

The database of the photos collected in MdP on the culture andon the natural environment in Luang Prabang, so-called“PhotoDB”, was created with thousands of photos which havebeen stored as the films or prints. The photoDB has beenequipped with photo management interface to enable thesearch by various criteria from number of photos, of which
detailed profiles are input and updated.

3. Motivation of the organizations to participate in databasecreation has been increased. Some institutions have been
motivated to provide their information to the ICT team fordatabase creation, and some other institutions have interest to
create their database application for their work.

4.3

4. The local communities have become more aware of public
information, by visiting ICT center and participating in giving
suggestion and discussion on the contents as public
information.

Also, further improvement of the user-interface of HeritageDB
was implemented utilizing the graphical documents such as
maps and drawings, in order to promote usability both forinternal and external use.
4.4

5. The activities have gathered attention from Internationalorganizations and Media. For instance, the activity on the ICT

Setting up the pilot ICT center

The database applications created in 4.1 and 4.3 ware opened

7

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

center was featured by the booklet, “Stay Another Day4”.

activities are operated vertically. In order to promote thesustainability of ICT development, close communication
among institutions is necessary. Especially, in the casemulti-institutional project, it is necessary to continuously
promote consensus by sharing information and mind,
since different institutions might have different interest
based on different scopes and culture.

6 Further issues
6.1

Lack of Internet environment

Since the technical documents, technical supports, and updated
programs of FOSS are available mainly in Internet, lack ofInternet has made difficulties on human resource development
and operation and maintenance of the system. Rapidly
improving Internet connectivity in Luang Prabang hopefully
will solve this problem.
6.2

Local ownership: It is often the case that the ownershipsof the activities belong to experts, government, or
agencies. But it is inevitable to enhance local ownershipto effectively continue the development activities.

Difficulty in multi-language support

8 Conclusion

Lao language support of the system is currently insufficient toobtain the maximum benefits from the database application.
Input and display of Lao characters are possible on Microsoft
Windows XP with additional software and fonts. Also, due tothe lack of human resource for translation among languages,
the multi-lingual contents development has been in the difficult
situation.
6.3

According to the needs assessment based on field survey, theproposed components for ICT development activities havebeen implemented, especially the author has been mainly
involved in the local human resource development,
development of database application and ICT center. Throughthe implementation, the activities have contributed to give thebasis and view for the further ICT development in Luang
Prabang. Especially, the local ICT team members who are
equipped both with ICT skills and insight of heritage
conservation and development can contribute to enable the
ICT development in the context of heritage centered regional
development.

Difficulty in self-sustainable financial cycle

The operation cost of ICT center was higher than expected. It
was revealed that not only the cost for expensive ICT
equipments, but also the running cost including employmentof the security guards to save the expensive equipments are
significant. Also, the number of visitors was not sufficient to
get the sufficient amount income by selling goods and by theadvertisement.

Furthermore, five important issues for sustainable
development were identified through the experience of thiscase study in particular attention to ICT field. Those issues
hopefully contribute to promote the sustainable developmentin further activities in Luang Prabang and in other places in theworld.

7 Discussion on sustainability
Through the experiences of this international
development project, the author has tried to clarify
important issues in the context of sustainable
development. The identified issues are listed as follows:

References[1] The Independent Commission for Worldwide
Telecommunications Development, (1984) “Missing

Integration: Sustainability cannot be realized only byfocusing on the area of expertise. When introducing
appropriate ICT, it is necessary to consider the various
related aspects, such as human resource, economics,
logistics, politics and laws, technology, and culture.

Link”, ITU, CH-1211 Geneva 20, Switzerland.
[2] N. Ooka, S. Yamaguchi, J. Takada, (2004) "Construction of
telecommunication network for rural development: Feasibility
study of ICT application to the World Heritage Site of Luang
Prabang", 3rd Great Asian Street Symposium, Singapore.

Adoption to change: Since the local situation is rapidlychanging, especially in developing counties it is also
required to continuously observe and adapt to thesechanges. Technology in not an exception, and available
and affordable hardware/software and infrastructure in
locality should be followed up.

[3] S. Yamaguchi, J. Takada, N. Ooka, J. Abe, (2005) “ICT
contribution to promote sustainable development, Case study in
the World Heritage site of Luang Prabang”, 8th World
Symposium of the Organization of World Heritage Cities, Cusco,
Peru.
[4] S. Yamaguchi, M. Mochizuki, N. Ooka, (2005) The
presentation slides of Tokyo Tech Team in Ministry ofInformation and Culture, Vientiane, Lao PDR

Needs based development activities: Any development
project can be truly effective only when it is based on thelocal needs of communities and government. Therefore,
careful needs assessment is a pre-requisite to thedevelopment project.
Information sharing: It is true that many development

4

Stay Another Day: http://www.stay-another-day.org/

8Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)
Study about Text Classification with Statistical Pattern Recognition
MethodStudent Number: 04M51254 Name: Yusuke WatadaSupervisor: Yukihiko Yamashita

統計的パターン認識手法を用いた文書分類に関する研究
和多田 佑介「文書分類タスク」は従来, 自然言語処理の分野で研究されてきたが, 近年では統計的パターン認識による手法も用いられている. String Subsequence Kernel(SSK) を用いた文書分類もその一つである. ところが, この手法の計算コストは非常に高いため, 近似的な手法を用いる必要がある. そのために, 特徴選択に工夫が必要となる. 本研究では, 特徴選択の方法を改良し性能向上を試みる. また SSK の特徴である言語に依存しないという点を利用し, スパムメール分類においても有効であることを実験で示す.

1 Introduction
The growth of the World Wide Web brings us the information society. A huge amount of information is generated, and we should process it every day. However,today we get so huge amount text data such as e-mails,
web pages, and so on that we cannot process them manually. Therefore we require a good method to process
them more quickly and automatically.
Natural Language Processing has been studied somany years to solve such a problem that we can process
natural language automatically with computers. For example, asking computers of some questions by speaking natural language or interpreting a text written inforeign language into another automatically. However,generally it is very difficult for computers to solve likethese kinds of problems, because natural language has
so much variety and requires a lot of vocabraly.
“Text classification” has traditionally been studiedin Natural Language Processing field. Recently it has
also been studied in Statistical Pattern Recognition field.
String Subsequence Kernel(SSK) is one of such methods, but it takes too high calculation cost. Although we
suppose approximated methods, another problem occured that we should select appropriate features. In thispaper, we modify the feature selection method to improve the performance. And we adapt SSK to the Spammail filter it never depends on any languages.
This paper consists of 6 Sections. In Section 1 we
explain intruduction of study. In Section 2 we explain
Text classification. In Section 3 we explain Support
Vector Machine. In Section 4 we explain SSK and provide a new feature selection method. In Section 5 weshow computer experiments and its results. In Section
6 we provide conclusion.

2 Text Classificationtext

categoryearnchampion products inc said its
board of directors approved a
two-for-one stock split of its
common shares for shareholdersof record as of april 1, 1987.
the company also said its board
voted to recommend toshareholders at the annual
meeting april 23 an increase in
the authorized capital stockfrom five mln to 25 mln shares.
reuterTopics of earning.

acq
Topics of acquisition of company.

corn
Topics of corn.

crude
Topics of crude oil.

Figure 1: Text classification“Text Classification” is a problem to map a certain
input text to a category which is defined in advance.
Traditionally, this kind of problem is solved by Natural Language Processing methods. At first, input texts
are split into a sequence of smaller units, such as words
or n-grams, and in each unit, we can calculate the probability of appearance from its appearance count. For
example, bayesian classifier is one of the famous textclassifiers. Naive bayes is the simplest bayesian classifier, and its expression is following:argmax P(Ci |S) =i=

P(Ci ∩ S)P(Ci )P(S)Π j P(w j )P(Ci ).
P(S)S is an input text, and Ci represents the classes for eachi. w j represents words which appeared in S for each j.
P(Ci |S) means the probability where we get S as input
text and S map to Ci . To determine the class which S
map to, we calculate P(Ci |S) for each i, and we selectthe most probable class of i.

93 Support Vector Machine

Then, we adapt TRIDE-2007-01,
the following kernel
Februayfunction:
22, 2006 (ISSN 1880-8468)

Technical Report of International Development Engineeringhφ (x), φ (y)i = K(x, y).

Class 1
Class 2

We can replace the kernel function from the innerproduct. If we use the kernel function, we can cut down
the calculation cost.
As kernel functions, generally used like followings.

Figure 2: Linear Classifier
Support Vector Machine(SVM) is a 2-class classifier which uses the feature vectors and a hyperplaneclassifier with a kernel function. In this kind of methods, we previously transform training/test patterns intofeature vectors, and using training vectors, learn the parameter in the following linear classification function:

K(x, y) =hx, yi: Identical.

K(x, y) =(hx, yi + c) p

: Polinomial.

K(x, y) =||x−y||2exp(− 2σ 2 ): Gaussian.

4 String Subsequence Kernel

In kernel method, we only have to calculate a kernel
function but the inner product itself. Then, we suppose
a new kernel function what we call String Subsequenced(x) = hw, xi + b.
Kernel(SSK). This kernel provides the similarity of two
x represents the input vector, w is a weight vector input texts by retrieving the feature quantity based on
constant, and b is a scalar constant. If this function re- following expression.
turns positive, then x map to class 1, and negative, class
2. If we may devide feature vectors into two groups by
Kn (s,t) = ∑ φu (s)φu (t)hyperplane, we have many possible hyperplanes. So weu∈∑nwould like to determine the most proper one by using= ∑ ∑ λ l(i) ∑ λ l( j)training vectors. This means we determine w and b.
=u∈∑n i:u=s[i]j:u=t[ j]∑n ∑∑ λ l(i)+l( j) .

u∈∑ i:u=s[i] j:u=t[ j]Lefthand means SSK of length-n between text s and textt. ∑ is a set of all alphabets, and ∑n means a set of all
substring of which length is n.i and j are substring index vector. For example, if substring s consists of first,
third and fourth letter from original string, i = (1, 3, 4)
Figure 3: Margin Maximization Principle
and s[i] means substring consists of i. If we get the twostring s and t, we retrieve all capable substrings from s
One of the most proper strategies for classification and t, and compare with each other. If the substrings areis to determin constants so as to maximize its margin. same, we calculate and add the feature quantity. λ is a
Margin is the minimum distance between hyperplane to decaydance parameter, and l(i), l( j) are length of eachthe most neighbor samples for each class. This strategy substring. This means the longer the substring is, theis called Optimal Hyperplane Classifier(OHC) and it is lower the feature quantity.
good at classifying unknown samples. However someproblems are still remained. Primary, if we have no hyc-a c-t a-t b-a b-t c-r a-r b-r
perplane which can devide the samples, we cannot use
φ (cat) λ 2 λ 3 λ 2
0
0
0
0
0
this kind of method. Secondary, when we calculate clasφ (car) λ 2
0
0
0
0
λ3 λ2
0
sification function, we have to obtain the inner productφ (bat)00 λ2 λ2 λ3000of an input vectors. Generally, the dimension of featureφ (bar)
0
0
0
λ2
0
0
λ2 λ3
vectors is very high, and calculation cost of inner product is also high. Thus, we cannot treat high dimensional
For example, this table describes the results of mapping of 4 strings, cat, car, bat, and bar, For each string
vector spaces in this method.

10the feature
valuesof length-2-substrings,
c-a, c-t, a-t, bTechnical
Report ofInternational
Development Engineeringa, b-t, c-r, a-r and b-r are written in the table elements.
In cat, c and a are first and second letter, so feature value
c-a is λ 2−1+1 = λ 2 . All the same to bar, b and r are first
and third letter, so feature value b-r is λ 3−1+1 = λ 3 .
Then, we can calculate kernel value of SSK like followings:• K2 (cat, car) = λ 2 × λ 2 = λ 4 .
• K2 (bat, bar) = λ 2 × λ 2 = λ 4 .

4.2 Feature Selection
TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)
When we select features, we should consider what feature is proper to use. Feature selection effects the classification directly, so we should use a good feature select strategy. In existing studies, features are selectedwith its appearance count as substring. We propose amethod that takes feature’s inner class variance into account. That is we use the features which is large feature
values in one class but small in others. This criterior is
implemented by the following expression.

• K2 (cat, bar) = 0.
• K2 (car, bat) = 0.
Cat and car have common substring c-a, so its kernel value is λ 4 . Bat and bar are similar. However, cat
and bar or car and bat have no common substrings, so
their kernel values are zero.
Furthermore, by improving the calculation form withdynamic program, we can much more reduce the calculation cost.
K0′ (s,t) = 1, for all s,t.
Ki′ (s,t) = 0, if min(|s|, |t|) < i.
Ki (s,t) = 0, if min(|s|, |t|) < i.
′
Ki (sx,t) = λ Ki′ (s,t)+ ∑ Ki−1(s,t[1 : j − 1])λ |t− j+2| ).
j:t j =xi = 1, ..., n − 1.

v i j = f i j − mi j +s( fk j − mi j )2 
.
L−1
k6=i

∑

(1)

fi j is a mean of jth feature in ith class, mi j is the averaged value of fi j where i 6= j, and L is a number ofclasses. For each i, we calculate vi j and enumerate themby descendant order of vi j . Finally we choose the mostlargest n features of the class i.

5 Experiments
5.1 Text Classification
In this problem, we use the 3000 training samples and1000 test samples from Reuters-21578 data set. There
are 4 categories, earn, acq, crude and corn in it.

′
Kn (sx,t) = Kn (s,t) + ∑ Kn−1(s,t[1 : j − 1])λ 2.

class
earn
acq
corn
crude

j:t j =xPractically we uses the following normalized SSK
because of not depending on input text length.
K(s,t)
K̂n (s,t) = p.
K(s, s)K(t,t)

data size
1759
905
104
232

existing
96.42
89.47
56.61
78.37

proposed
96.35
91.95
72.35
81.69

The table is the comparrison between the existing
feture selection strategy and the proposed one. The results are expressed by F1 value. F1 value is calculated
4.1 Approximated SSK
by the following expression where p is precision and r
As we mentioned above, SSK takes too high calcula- is recall.
tion cost. Then we use the approximated SSK. In this2prmethod, previously we select the substrings and for each
.
F1 =
p
+r
input text, calclate the SSK with them. Generally, only
a few features contribute the performance in classificaFrom this result, we can say that we get good pertion but others are not so much. Therefore we can take formance in 3 data sets, acq, corn, crude.
approximate strategy in SSK. Each SSK value is the approximated feature value of SSK. Following expression 5.2 Spam Filter
is a function of approxmated SSK.
Secondly, we adapt the approximate SSK to mail classification. Recently spammers make complex tricks on
Kn (x, y) = ∑ K(x, s)K(y, s).
s∈S
junk mails and we cannot classify them by simple spamn
S represents a set of substring as SSK such as S ⊂ ∑ In filters which are based on word features. However, we
this strategy, we should determin a proper substring set consider that SSK will achive a good performance inS.
junk mail recognition, because SSK can treat not only

11

continuous
also noncontinuous
Technical
Report ofsubstrings
Internationalbut
Development
Engineering substrings.

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

In this experiment, we compare the performance among
bsfilter and SSKs with the existing and proposed feature selection strategies. The bsfilter is the opensource
bayesian filter software.
methodbsfilter
Existing
Proposed

accuracy
93.31± 2.66
90.51± 7.26
97.12 ± 1.21

From this result, we can say proposed method achive
better performance than bsfilter.

6 Conclusion
In this paper, we explained about Text Classification,
SVM SSK. SSK needs too high calculation cost to use
practically. Then we use the approximation of SSK. In
this strategy, we use a few features for classification.
However, it brings another problem that we should select proper features from all. Existing work do this task
by only selecting a set of n-grams depending on its appearance count. The new method we proposed is a newfeature selecting method. And we showed its advantages by experimental result of text classification andspam filtering problems. As future tasks, we should verify the results more precisely, introduce orthogonalization methods such as GSO, KL-extraction, and reduce
calculation cost.

References[1] Text Classification using String Kernel. Huma
Lodhi, Journal of Machine Learning Research2(2002)[2] SpamCop: A Spam Classification & Organization
Program, Patrick Pantel(1998)[3] David Haussler Convolution Kernels on Discrete
Structures(1999)

12

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

Extraction of liquid organic multi-component mixture with ionic liquids
as solvents
Student ID: 05M18018 Name: Kenya AKAISHI

Supervisor: Ryuichi Egashira & Sachio Hirose

イオン液体を溶媒とした液状有機系多成分混合物の抽出
赤石 憲哉
吸収油中の同素環式化合物（ナフタレンなど）に比較して含窒素複素環式化合物（インドールなど）の，改質ガ
ソリン中の飽和鎖式化合物（ヘキサン）に比較して芳香族化合物（ベンゼンなど）が選択的にイオン液体中に抽出された．従来用いられている抽出溶媒に比較して，イオン液体を溶媒とした場合においては収率が高く，分離の選択度は低かった．どちらの原料混合物においても，抽出相中からの抽出成分の分離に関して，イオン液体の分解点が抽出成分の沸点範囲より高いことから，蒸留を用いた分離操作が可能であった．
1. Introduction
Ionic liquids are organic salts that are liquid under 373
K. Ionic liquids have quite low volatility andflammability[4], they are inactive against water and theoxygen.
It has been already reported that the nitrogencompounds and the aromatic hydrocarbons wereseparated with ionic liquids as solvents from
nitrogen/homocyclic compound model mixtures[5],aromatic/aliphatic hydrocarbon model mixtures[2]respectively. However there is no report of the separationfrom an actual multi-component mixture.
In this work, the ionic liquids were applied to theseparation of the actual liquid organic multi-component
mixtures as extraction solvents. This would be done by
the confirmation of the thermal stability of the ionicliquids and the equilibrium extraction of the organic
multi-component oil mixtures. And the separation
process in which the extracted components wereseparated from the extract phase was considered.

Hydrophilic Ionic Liquids:

BF4
BF4
1-butyl-3-methylimidazolium
tetrafluoroborate

N,N-Diethyl-N-methyl-N-(2-methoxyethyl)ammonium
tetrafluoroborate

Hydrophobic Ionic Liquids:

PF6
1-butyl-3-methylimidazolium
hexafluorophosphate

N,N-Diethyl-N-methyl-N-(2-methoxyethyl)ammoniumbis(trifluoromethanesulfonyl)imide

Figure1: The structure of ionic liquids used in this research

2.3. Equilibrium Extraction
All organic multi-component mixtures were brought
into contact with each ionic liquid in a 3.5 x 10-6 m3 vial.
They were agitated intensely. After settlement, themiscibility between organic multi-component mixtures
and the ionic liquids were observed. The immiscible
combinations of organic multi-component mixtures andionic liquids were employed in extraction experiments.
The experimental conditions are shown in Table 1. The
samples were stirred for 120 hours at room temperature
to attain the equilibrium.

2. Experimental
2.1. Material System
Four types of ionic liquids were selected as extraction
solvents. 1-butyl-3-methyl-imidazolium tetrafluoroborate,[Bmim][BF4]; 1-butyl-3-methylimidazolium hexafluorophosphate, [Bmim][PF6], were purchased from Wako
Chemical Co. Ltd.
N,N-Diethyl-N-methyl-N-(2-methoxyethyl)ammonium
tetrafluoroborate,[Et2MeMeON][BF4];
and
N,N-Diethyl-N-methyl-N-(2-methoxyethyl)ammoniumbis(trifluoromethanesulfonyl)imide,[Et2MeMeON][Tf2N], were provided by Nisshinbo Ind. Inc. (Tokyo,
Japan). The structure of these ionic liquids is shown in
Figure 1. The actual organic multi-component mixtures
were coal tar absorption oil (AO), reformate gasoline(reformate) and light cycle oil (LCO).

Table 1: Experimental Conditions
Feed, R0[kg] 0.0010 - 0.0017
Solvent, E0 [kg] 0.0015 - 0.0023
E0/R0
0.5 - 1.0
[-]
Time
120
[h]
Temperature [K]296 ± 1

The compositions in the liquid phase were determinedby gas chromatography (GC17A, Shimazu) (GC) and
Karl-Fisher Titration (Metrohm. Ltd.) (KF). Those in theextract phase were also determined by GC and KF after
the reverse extraction was carried out with organic
solvents.

2.2. Thermal stability of ionic liquids
The thermal stabilities of ionic liquids were measured
by thermogravimetry (TG8120, Rigaku Denki Co. Ltd.)(TG).

3. Results and Discussion
3.1. Thermal stability of ionic liquids
Figure 2 shows the results of the TG analysis. The ionicliquids began to decompose around 600 K. All ionicliquids were thermally stable upto around 600 K.
13

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

Intensity [mV]

Weight loss [-]20

0[Et2MeMeON][BF4][Bmim][PF6]

3
2
1

15

5

6

1 Hexane
2 Benzene
3 Toluene
4 m-Xylene
5 Decane
6 Tridecane

4
10

5

[Et2MeMeON][Tf2N]
-1

0
0.0

[Bmim][BF4]

300

400

500

600

700

3.3. Equilibrium Extraction
The results of the miscibility were listed in Table 3. The
hydrophilicionicliquids;[Bmim][BF4]and[Et2MeMeON][BF4] were miscible with water, while thehydrophobicionicliquids;[Bmim][PF6]and[Et2MeMeON][Tf2N] were immiscible with water. All
the ILs were immiscible with reformate and LCO.

Intensity [mV]1

Table 3: The miscibility Yi of ionic liquids
with water, absorption oil, reformate and light cycle oil(O-miscible, X-immiscible)

5

AO

reformate

LC O

MeOH

EtOH

AcetonenC 6Toluene[BMIM] [BF 4 ]

O

X

X

X

O

O

O

X

X

[BMIM] [PF 6 ]

X

O

X

X

O

X

O

X

X

[Et2 MeMeON] [BF 4 ]

O

X

X

X

O

△

O

X

N.A.

[Et2 MeMeON] [Tf2 N]

X

O

X

X

O

O

O

X

N.A.

1 Toluene
2 Decane
3 Naphthalene
4 Quinoline
5 Isoquinoline
6 Indole
7 1-Methylnaphthalene
8 2-Methylnaphthalene
9 Biphenyl
10 Dibenzofuran

8
2
7
9

34

traditional so lvents
H2 O

0.6

(b)
Figure 3: the gas chromatograms of feed oils(a) AO; (b) reformate

Figure 2: TG diagram of ionic liquids; the ordinate is theweight loss of samples and the transverse is the temperature.

anion

0.4

Time [h]Temperature [K]catio n

0.2

800

5

0
0.0

10

6
0.5

1.0

Time [h](a)The compositions of feed oils are shown in Table 2, andthe gas chromatograms of feed oils are shown in Figure 3.
The feed oils contain a lot of components. Among those,the components listed in Table 2 were selected to be
identified and studied. And the mass fractions of water
were 0.001 and 0.0005 in AO and reformate respectively.

20

Intensity [mV]1

Table 2: The compositions of feed oils; AO and reformate
Absorption Oilconstituent abbrev.
Quinoline
Q
Isoquinoline
IQ
IndoleI
NaphthaleneN
1-Methylnaphthalene
1MN
2-Methylnaphthalene
2MN
BiphenylBP
Dibenzofuran
DBF

Intensity [mV]1
5

mass fra.
0.08
0.02
0.05
0.03
0.11
0.25
0.07
0.12

Reformateconstituent abbrev.
BenzeneB
TolueneTX
m -XyleneHn -Hexane

mass fra.
0.05
0.23
0.19
0.04

2

4

9

0
0.0

5

2

6

10

5

0.2

0.4

0.6

(b)
Figure 4: the gas chromatograms of the raffinate(a) AO; (b) reformate

The gas chromatograms of the raffinate are shown in
Figure 4. Compared with the gas chromatograms of feedoils, the peaks of nitrogen compounds in AO and thoseof aromatic hydrocarbons in reformate became smaller,
hence the nitrogen compounds and aromatichydrocarbons preferentially transferred into the extract
phase.
It was difficult to analyze the extract phase by means ofreverse extraction method, because of the polarity, thestructures of ionic liquids and so on, so that the materialbalance relationship during the runs was not sound. The
discussion was based on the composition in the raffinate
hereafter. Figures 5 and 6 show the yields of each

10

6
0.5

4

1 Hexane
2 Benzene
3 Toluene
4 m-Xylene
5 Decane
6 Tridecane

Time [h]

7

3

5

15

0
0.0

1 Toluene
2 Decane
3 Naphthalene
4 Quinoline
5 Isoquinoline
6 Indole
7 1-Methylnaphthalene
8 2-Methylnaphthalene
9 Biphenyl
10 Dibenzofuran

8

3

1.0

Time [h](a)14Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

component, Yi, plotted against solvent component in thecase of AO and reformate respectively. The yield was
defined by equation 1.
Rxi
Yi = 1 R0 x0
Equation 1

The definition of the separation selectivity is shown inequation 2. The separation selectivities of eachcompounds in AO and reformate are shown in Figure 7
and 8 respectively.

In case of AO, the yields with the solvent of ionicliquids were higher than those with the aqueous
methanolic solution[1]. The yields of indole with bothionic liquids reached up to 90 percents. In case ofreformate, the yields of hexane using ionic liquids weremuch higher than using sulfolane[3]. The yields ofaromatic hydrocarbons with ionic liquids were in therange from 30 to 50 percents.

Since the separation selectivities of nitrogen compounds
and aromatic hydrocarbons were higher than the unity incase of AO and reformate respectively, the nitrogencompounds and aromatic hydrocarbons could be
separated selectively from AO and reformate
respectively. The selectivities with ionic liquids werelower than those with conventional solvents. In AO,
indole was extracted more preferentially into the extract
phase than the other nitrogen compounds. The separation
selectivities of the nitrogen compounds with the ionicliquids constituted a reversal unexpectedly from thosewith the aqueous methanolic solution, which means theindole had a higher selectivity than the others using theionic liquids, on the other, hand it had a lower using theaqueous methanolic solution. In reformate, theselectivities with hydrophobic ionic liquids were higher
than those with hydrophilic ionic liquids.

βi , j =

[B
m

im

[B
m
][
B

F4

im

]

[B
][
P

F6

F4

]

[T
f

2N

]

Me

OH

]

Solvent

Figure 5: the yields, Yi of AO compositions○-Quinoline, ▽-Isoquinoline, □-Indole,●-Naphthalene, ▼- 1-Methylnaphthalene,■-2-Methylnaphthalene, ◆-Biphenyl,▲-Dibenzofuran

Equation 2

101

100

Miscible[Et2MeMeON]
0

Yj

Miscible

Selectivity of component i b i,2MN [-]Miscible

Yield Yi [-]Miscible

1

Yi[Et2MeMeON]
[B
m

im

[B
][
BF

m

im

4]

[B
][
P

F6

F4

[T

]

f2

N]

Me

OH

]

Solvent

Figure 7: the separation selectivity in case of AO.
○-Quinoline, ▽-Isoquinoline, □-Indole

1

Yield Yi [-]Selectivity of component i b i,H [-][Et2 MeMeON]

0

[B

m

im

[B
][

BF

m

4]

im

[B
][

PF

F4

]

[T

f2

N]

SU

LF

6]

Solvent

Figure 6: the yields, Yi of reformate compositions○-Benzene, ▽-Toluene, □-m-Xylene, ●-Hexane

101

100

[Et2MeMeON]
[B
m

im

[B
][B

F4

m
]

im

[B
][
P

F6

F4

]

[T
f2

N]

SU

LF

]

Solvent
Figure 7: the separation selectivity in case of reformate○-Benzene, ▽-Toluene, □-m-Xylene
15Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

7. Literature Cited[1] C. Salim
Doctoral thesis, Tokyo Tech, 2005
[2] G. Wytze Meindersma et. al.
Fuel Processing Technology, 87, 59-70, 2000
[3] H. Habaki
Master thesis, Tokyo Tech, 1995
[4] Jonathan G. Huddleston et. al.
Chemical Communication, 1765-1766, 1998
[5] M. Matsumoto et. al.
Journal of Japan Petroleum Institute, 49, (5), 256-261,
2006

Raffinate

Distillation

Feed

Extraction

3.3. Extraction Process
The separation of the extracted components from theextract phase is required after the extraction. The simple
scheme of the separation process is shown in Figure 8.
The According to TG results, ionic liquids wereremarkably thermally stable upto around 600 K. And
there is the difference between the decomposition
temperature of the ionic liquids and the boiling range ofextracted components. The extracted components inextract phase could be separated from extract phase andthe ionic liquid could be recovered by the distillation.

Extracted
Components

Solvent

Figure 8: the simple scheme of the separation process

4. Conclusions
The nitrogen compounds and the aromatichydrocarbons were separated from actual organic
multi-component mixtures. The ionic liquids used in thisresearch have the thermal stability upto around 600 K.
The solvent extraction could be operated owing to theimmiscibility of the ionic liquids. The nitrogencompounds and the aromatic hydrocarbons wereseparated preferentially from AO and reformate
respectively. The yields with ionic liquids were higher
than with the conventional solvents. The separation
selectivities with ionic liquids were inferior to those withconventional
solvents.
The
selectivities
withhydrophobic ionic liquids were higher than withhydrophilic ionic liquids in case of reformate. After the
extraction, it could be possible to separate thecomponents in extract and to recover ionic liquids by thedistillation
5. Acknowledgment
The author wish to thank Nisshinbo Co., Ltd., JFE
Chemical Co., Ltd., and Idemitsu Kosan Co., Ltd. forproviding the materials.
6. Nomenclatures
R: mass of raffinate
R0: mass of feed
xi: mass fraction of component i in raffinate phasex0,i: mass fraction of component i in feed
Yi: yield of component i<subscripts>
i: component i
j: component j
w: water<greek>β: separation selectivity

[kg]
[kg]
[-]
[-]
[-]

[-]

16

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

Moving picture coding with wavelet transform
Student Number：05M18024

Name：Takashi OHTA

Super visor：Yukihiko YAMASHITA

ウェーブレット変換を組み合わせた動画像符号化に関する研究
太田 昂志
MPEG に代表されるブロック単位の動画像圧縮方式には、復号した画像にブロック歪が現われる問題が生じる。この問題を解決するために「ウェーブレット変換を組み合わせた動画像符号化」を提案する。この手法は復号画像のブロック歪みを減らすだけでなく，予測画像にデブロッキングフィルタを適用することにより，差分画像におけるブロック歪みも減らすことで, より効率的な符号化が期待できる。

1

Introduction

cess per block.

In recent years, the information communication technology called IT (Information Technology) is splendidly
developed. Not only a document but also a picture, asound, etc. are widely used for communication. How-

2 MPEG
MPEG is a standard of the multimedia coding for ac-

ever, the amount of information of the multimedia data cumulation media, broadcast, communication, etc. It
which generally contains digitized pictures and sounds mainly consists of three regulations, such as the regulaare huge. Therefore, in order to treat the information, tion on the coding method of a video signal, the reguthe transmission way and the mass storage medium of
a broadband are needed. Then, the research on reduc-lation on the coding method of an audio signal, and forboth integrated method.

tion of the information by the data compression of pictures or sounds so-called compression coding came toVideo data is realized by the set of the still picture located in a line on the time-axis generally called a frame.
MPEG performs compression coding by reducing thosebe popular for the purpose of efﬁcient use. Although
the broadband communication can be used in wide area,
narrowband communication is also used. Then, efﬁ-

spatial and time redundancies. Reduction of spatial redundancy is called the coding in a frame, and performs

DCT (discrete cosine transform), quantization, and coding every 8×8 [pixel] block. Moreover, reduction ofducing the redundancy included in video data. There time redundancy is called the coding between frames,
are two kinds of redundancy of video data. One is the and is performed using the technique of motion comspatial redundancy and the other is time redundancy. pensation prediction. This is extracting and treating theciency of image coding has to be increased.
Video coding is the method of compressing by re-motion information on a certain domain in a picture in
two near frames in time. Generally, block matchingtional standard system of still picture coding is a tech- is performed by every 16×16 [pixel] block. And thenique used very widely. And the algorithm called mo- motion vectors which express motion information is extion compensation prediction is used in the latter re- tracted. The general procedure of motion compensation
duction. This is used in MPEG (Moving Picture Ex- prediction is shown below.
The former is also used in still picture coding. JPEG
(JointPhotographic Experts Group) which is an interna-

perts Group) which is the international-standard system

1. Extraction of moving vector by comparing theframe for encoding with the reference frameof video coding. It became possible to code and compress video data at high efﬁciency by these two techniques compared with only the former. However, suchalgorithm had the problem in a decoding picture. Visualdegradation called block distortion is produced from pro-17

2. Generation of the prediction picture by the motion vector and the reference frameTechnical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

3. Generation of the picture of difference by the substracting of the prediction picture from the frame
for coding

surroundings of the pixel which a target of the ﬁlter
whether the area is an edge or not. And, in the predicted
image we calculate the difference of the pixel values inwhich the boudary between blocks. By the result, we
4. Coding of the motion vector and the picture ofchange the direction and the strength of the ﬁlter. The
difference
algorithm of motion compensation predition with deblocking ﬁlter is as fellows.
5. Execution of 1. to the following two frames
The main coding parts of MPEG are realized with the
combination of the coding between frames and inside

1. A target block is made to match with a reference
frame using the block of 16×16[pixel].

of a frame. First of all as a basic procedure of MPEG,
the frame of the beginning of video or the frame used

2. the predicted image is created by the motion vector.

as a starting point is coded in a frame. These frames
are called I-picture. I-picture which had conversionquantization performed here is reverse converged by the

3. In the predicted image the deblocking ﬁlter to eachpixels on the boundary of the block is applyed.

local decoder, and is temporarily memorized by the frame
memory. Next, frames other than I-picture are coded

Except the boundary between two blocks that havethe same motion vector.

between frames which uses motion compensation preFor example, show the calcuration for p0,0 in the Fig.1.
diction. There are P-picture and B-picture in these frames. Calcuration ﬂow of deblocking ﬁlter for p0,0 is as felThe obtained prediction picture which is made based on lows.
a vector and it, and the picture of the difference of a target frame are coded. This reduces time redundancy and1. Start condition, state= 0, a, b, c, d = const
2. i f (|p0,−1 − p0,0 | < a, and |p0,0 − p0,1 | < b),state+1

more efﬁcient compression is performed.

3

3. i f (|p−1,−1 − p0,0 | < a, and |p0,0 − p1,1 | < b),state+2

Moving picture coding withwavelet transformIn the coding using the block matching method repre-

4. i f (|p1,−1 − p0,0 | < a, and |p0,0 − p−1,1 | < b),state+4

sented by MPEG, the problem of producing block distortion which is visual lattice-like degradation aooears
in the reconstructed video. On the other hand, in the

5. i f state= 1, 7
p0,0 = (p0,−2 + p0,−1 × 2 + p0,0 × 2 + p0,1 × 2 +p0,2 + 4) ÷ 8
p0,−1 = (p0,−2 + p0,−1 + p0,0 + p0,1 + 2) ÷ 4

ﬁeld of still picture coding, the coding which reduces
block distortion using wavelet transform and subband
conversion as a technique has been advocated. Then,in order to control block distortion produced in video,
we try to apply the coding (wavelet coding) which uses

6. i f state= 2p0,0 = (p−1,−1 + p0,0 × 2 + p1,1 + 2) ÷ 4
7. i f state= 3

wavelet transform. However, the conventional motion
compensation prediction which uses block matching is

p0,0 = (p0,−2 + p−1,−1 + p0,−1 + p0,0 × 2 + p0,1
+ p1,1 + p0,2 + 4) ÷ 8
p0,−1 = (p0,−2 + p0,−1 + p0,0 + p0,1 + 2) ÷ 4

not efﬁcient. Then, we propose appling the deblocking
ﬁlter to the motioncompensation prediction.

3.1

8. i f state= 4
p0,0 = (p1,−1 + p0,0 × 2 + p−1,1 + 2) ÷ 4

Deblocking ﬁlter

9. i f state= 5
At ﬁrst, we use the conventional block matching method
to predict the latter frame. Next, we use the deblocking ﬁlter to the predicted image. Then we examine the18

p0,0 = (p0,−2 + p1,−1 + p0,−1 + p0,0 × 2 + p0,1
+ p−1,1 + p0,2 + 4) ÷ 8
p0,−1 = (p0,−2 + p0,−1 + p0,0 + p0,1 + 2) ÷ 4

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

Fig. 1: Example of calucuration

Fig. 3: Rotating disk

Fig. 2: Diagram of proposed encoding method
10. i f state= 0, 6 Don’t apply ﬁlter.

3.2

Wavelet coding

Fig. 4: PSNR for different bit rates

Wavelet coding is used for the compression codingin the case of mainly inclining information for every MPEG2 used as codec of Hi-Vision or standard qualityfrequency band. When especially energy inclines to- of-image television. And quality-of-image degradation
ward the low cycle, it is known that very efﬁciencient of a decoding picture was evaluated numerically. PSNR
coding is possible. There is the coding method called is used for numerical assessment. PSNR is given by the
SPIHT (Set Partitioning In Hierarchical Trees) which is following formulas here.
specialized in the tree structure of wavelet. This coding method is used for the coding part of the proposal
technique. Correlation will become high supposing the
2-dimensional motion vector obtained by motion compensation prediction of a pixel unit is extracting the mo-

MSE =

255
PSNR(dB) = 20 log
MSE√−1 H−102∑W
x=0 ∑y=0 ( f (x, y) − f (x, y))
WH

(1)

(2)

In coding of I-picture, the block distortion is protion correctly. Therefore, motion information inclines
toward lower frequency domain. Moreover, the abso- duced by MPEG which performs block processing inlute value of the pixel of the picture of difference be- a decoding picture. However, block distortion is notcomes very small. Therefore, both information is com- produced by the proposal technique. By comparisonpressible at high efﬁciency using the above-mentioned by PSNR, it turned out in the amount of marks withthe same proposal technique that a picture better thanwave let coding.
MPEG is obtained. This result is shown in Fig.4.
In coding of P-picture, we compare the decoded im-

3.3

Image coding experimentages of MPEG2 and proposed. Fig.5 and Fig.6 are deIn an experiment, two frames of the standard video coded images which are zoomed by four times. In thesequence for assessment are used for I-picture and P- former image block distortion is seen, and in the latter
picture, respectively. And coding and decoding were is not. Deblocking ﬁlter disturbed that an unnecessary
performed by the proposal technique. The example of a high frequency elements included in the difference image. And, Fig.4 shows the good result at proposed P

frame is shown in Fig.3. Moreover, it is compared with19

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

References

picture. Especially in the low bit rate condition.

[1] Yoshinori Sakai，Toshiyuki Yoshida：“Image information encoding”, Ohmsha, 2001.
[2] Sadayasu Ono，Junji Suzuki：“Achievement
method of comprehensible JPEG/MPEG2”，
Ohmsha, 1995.
[3] Susumu Sakakibara：“Wavelet beginner’s guide”，
Tokyo Electrical Engineering College Publications
Service，1995.
[4] Amir Said, William A.Pearlman：“A New Fast andEfﬁcient Image Codec Based on Set Partitioningin Hierarchical Trees”，IEEE Trans. Circuits andSystems for Video Technology, vol.6, pp.243-250,
June. 1996.

Fig. 5: Decoded P picture of MPEG2

Fig. 6: Decoded P picture of proposed method

4Conclusions
In this paper, the video coding method using deblock-ing ﬁlter and waverlet tranceform was proposed, thecomputer experiment was conducted, and the good result was shown with the experiment especially in the
low bit rate condition. For future work, we have to develop the block noise detection method that provides
more highly efﬁcients, and a more effcient coding method
for the difference image.

20

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

Creation and Evaluation of e-Learning materials
for Introduction of Engineering
Student Number: 05M18030

Name: Masanori Kaneko Supervisor: Jun-ichi TAKADA

工学導入教育のための e-Learning 教材の作成とその評価
金子 昌永
最先端の科学技術研究は学問の深化・細分化を生んでおり, 低学年の学部学生にとって, 理数系基礎科目との関連性を理解するのは難しい. また, 日常生活に浸透した最先端技術はブラックボックス化されており, 使用者にとっての科学的興味の対象ではなくなっている. よって, 高校の理数系教育から工学教育への橋渡しを行う「工学導入教育」を提案し, 教材を作成した. また, 教材の形態は web コンテンツとした. これは, 主として学生の興味を喚起するリッチコンテンツを扱うためである. 本学の初年度学生を対象として, 作成した教材の試験評価を行った. 教材の消化時間およびその理解度を解析し, 学習前後の興味・意欲の変化をアンケート集計した.

1

Introduction

Since the latter part of the 20th century, there was aconsiderable growth in the various fields of engineering. As the latest scientific technologies are making remarkable progress and those corresponding academic
fields are becoming more advanced and complicated,
the gaps between the curriculum in senior high school
and the curriculum for engineering course in university are expanding. Moreover, another gap between
university education and practical business in company is also broadening. Furthermore, the scientificcourses in university are at much higher levels thanthose of high schools to produce superior scientists orengineers, and it is considered an important problem.
For this reason, it is difficult for first year bachelor
students to understand the relation between the current curriculum and the latest technologies. The effectis that first year bachelor students will find it diffucult
to be motivated in their study and to be interested inscience and technology.
Therefore, new learning materials are needed for students to help them to become motivated and interested. The materials also aim in filling the gap between the current curriculum and the latest technology.
Creating web content as an auxiliary learning material is one of the best solutions. It is not a traditional teacher-to-student lecture style, but students
can learn by themselves, anywhere, and at any time.
For this reason, the contents should be more attractive than a class lecture. Sounds and moving images
which are available in the web are persuasive and interactive contents which can make students actively
participate can be used.
In this paper, our interest is to create effective eLearning contents, and to analyze the effectivity. The
outline of this paper is as follows. Section 2 describes
the selected topics about this theme and introduce
the supposed e-Learning contents in detail. Section 3
represents data fitting model, which shows the resultsof user test. And finally, conclusion and future works
are made in section 5 and 6.

2

Selection of Study Topics

Requirements of topics are listed as below.
• The topics are familiar to students.
• The topics contains the subjects for 2nd year orlater.
• A prerequisite is an understanding of math andscience at the senior high school level.
For these reasons, the topic “Voice communication incellular phone” has been chosen in this study. Since
cellular phone is very familiar for students. It also contains the elements of the subjects for 2nd year or later.
(Signal processing, digital radio communication). It
is possible to describe this topic without any difficult
mathematical explanation as well.
The e-Learning material of “Voice communication incellular phone” has 9 sections.

2.1

Detail of Contents(1)content type: text, PNG image
This section explains the difference between analog and digital values (Figure 1).

Figure 1: Samples of analog clock and digital clock21

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)(2)content type: text, GIF animation
This section displays A/D conversion, sampling,and quantization. Students learn what is digital
signal (Figure 2).

Figure 4: Samples on information quantity calculator(7)content type: text, GIF animation
This section describes the non-linear quantization. The non-linear quantization for digital soundis one of information compression which can reproduce the characteristic of human ears (Figure
5).

Figure 2: Sampling simulation of sine waves(3)content type: text, GIF animation
This section explains that the quality of voice depends on sampling frequency and quantized bits.
It also mentions about the Nyquist-Shannon sampling theorem, but not mathematically.
(4)content type: text, wave sound
Related to the previous section, this section shows
several kinds of wave sounds (.wav), with eachof them having different sampling frequency andquantized bits (Figure 3).

Figure 5: Non-linear quantization(8)content type: text, Png image, soundThis section shows the waveform of a piano song,and makes its sound available (Figs. 6and7).
Students learn that all the waves consist of various sine waves.

Figure 3: Samples on the difference of sound quality(5)content type: text, Javascript
This section explains what is information quantity. Students can select sampling frequency andquantized bits and calculate the information quantity per a second by Javascript. These sampling
frequencies and quantized bits take the realistic
values for CD-Audio, DVD-Audio, MD-Audio,
and telephone, to make this content more familiar (Figure 4).
(6)content type: text

Figure 6: Time domain waveform(9)content type: Java applet
This section has a Java applet which is persuasive for explaining Fourier expansion. It shows awaveform which consists of five sine waves. Amplitude and frequency can of each sine wave beedited independently.

This section explains an easy example of information compression. Students learn the redundancy of information data.

22

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

Figure 7: Frequency domain waveform
Figure 9: the difference of elapsed time CDF F (t)After obtaining α and β of each contents, they arenormalized to get parameters α∗ and β ∗ . We can findcontents characters by plotting α∗ and β ∗ in a 2Dgraph. In Table 1, the contents characteristics by α∗
and β ∗ are given.
α∗ < 0, β ∗ > 0
simple, difficultα∗ < 0, β ∗ < 0simple, easy
Figure 8: summed up sine waveformα∗ > 0, β ∗ > 0
complicated, difficult
α∗ > 0, β ∗ < 0complicated, easy

Table 1: Contents characteristics by α∗ and β ∗3Analysis of Elapsed Time for 3.2 Data Analysis
Individual Contents

3.1

Data Fitting Model

Test UsersOS
Web browser
Sound playerJava VM

To analyze the test results, we introduce the Ueno’s
method [2]. According to Ueno’s method, the appearance of elapsed time for studying a certain content
follows the gamma distribution as follows (1).
f (t) =tα−1
β α (α − 1)!

texp(− )
β(1)

8 members of 1st year bachelor
student in Tokyo Tech.
MacOS X Leopard
Safari 2.0.1
iTunes 7
J2SE 5.0 Runtime

Table 2: Environment of user test

Table 2 shows the environment of user test. All the
Where α is a plain thinkin g process, and β is acontents do not depend on the platform.
required time for each process.
The parameter α shows whether the content requires
section# αβτ (sec) α∗β∗simple thinking or complicated thinking. The param1
9
8.56
77
0
-0.37
eter β shows the difficulty of the content. Therefore,2
10 12.25 122.5
0.1 -0.09
it can be shown as the equation (2)
3
14
7.1
99.38
0.5 -0.49
4
4 19.6378.5-0.5 0.48
αβ = τ
(2)
5
13 7.1492.88
0.4 -0.48
where τ is given as the mean of elapsed time for each
6
2
23.5
47
-0.7 0.78
student.
7
19 3.51
66.63
1
-0.76
To fit the collected datas to the model, the CDF of
8
6 12.48 74.88 -0.3 -0.07
elapsed time F (t) is
9
4 26.34
105
-0.5
1
½
0
(t < t0 )Rt
F (t) =
(3)
f
(t)dt
(t
≥ t0 )
Table 3: Result of elapsed time analysis
0
Where t0 shows the theoretical minimum elapsed
Table 3 shows the analyzed results of contents. As
time of the target learning content.
Figure 9 illustrates the dependence on parameter α. shown in Table 3,• High α appeared in section(3,5,7). These contents require more complicated thinking thanother contents.

Assuming the elapsed time as PDF f (t), we can estimate α and β which are best fit to F (t).

23Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)• High β is also seen in section(4,6,9). These con- [2] Maomi Ueno, Keizo Nagaoka: “Online analysis of elapsed time for e-Learning material by ustents requires more diffucult thinking than othering gamma distribution”, Japan Society for Edicacontents.
tional Technology Journal 29(2) 2005.
1
9:applet
0.8
6:txt
0.6

0.4

4:wav

β

*

0.2

0
8:txt+img+mp3−0.2

2:txt+gif−0.4
1:txt+img−0.6

5:txt+js
3:txt+gif

7:txt+gif

0.4

0.8

−0.8

−1
−1

−0.8

−0.6

−0.4

−0.2

0
α*

0.2

0.6

1

Figure 10: α∗ -β ∗ chart
Figure 10 shows the α∗ -β ∗ chart of all the contents.
As shown in figure 10, there are no relation between
content type and the parameter α∗ and β ∗ .

4Summary

The e-Learning web contents which have multimedia
and interactivity are created. And through the user
test, these contents are analyzed and evaluated by using Ueno’s method. The result of the user test revealed the difficulty and elapsed time for each contents.

5

Future Works
More topics and contents
The chosen topic “Voice communication in cellular phone” is an alternative plan of the topics,
and further subjects can be added (e.g. electro
magnetic wave, introduction of antenna, etc).
More test users
Because the number of test user is only 8, thedata is considered not fully reliable. According
to Ueno[2]’s user test, 78 is enough to have the
gamma distribution CDF F (t)(3), but 15 is notenough.
Online analysis system
As the number of contents and test users becomelarge, it requires much more time for analysis.
But since these contents are on the web. It isable to calculate and analyze the elapse time on
the server automatically in real time.

References[1] Walter Dick, Lou Carey, James O. Carey: “The
Systematic Design of Instruction”

24Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

Study on applicability of Steel Slag Hydrated Matrix
to steel reinforced members under marine environment
Student Number: 05M18047 Name: Tomonari KIMURA Supervisor: Nobuaki OTSUKI
海洋環境下における鉄鋼スラグ水和固化体の有筋部材への適用可能性に関する研究木村 智成
鉄鋼生産において副産物として生じる製鋼スラグと高炉スラグ微粉末を主要原料とした、環境低負荷の鉄鋼
スラグ水和固化体が土木建設材料として開発された。しかしながらその適用領域は港湾土木構造物の無筋部材に限定されている。本研究では、この鉄鋼スラグ水和固化体を有筋部材に適用する為に、物質透過に対する抵抗性及び鉄筋腐食に対する抵抗性に関する検討を行った。その結果、現在使用されている港湾土木構造物の鉄筋コンクリート部材と同等の、①塩化物イオン浸透に対する抵抗性、②酸素透過に対する抵抗性、③固化体内部の鉄筋腐食に対する抵抗性が認められ、有筋部材としての適用可能性が示された。
aggressive substances and corrosion are currently
limited. Therefore, to clarify the applicability of SSHM
in steel reinforced materials under marine
environment, this research study presents thefollowing objectives; ①to investigate the penetrationof aggressive substances into SSHM, ②to investigate
the corrosion resistance of reinforcement inside
SSHM, and ③to evaluate the durability of reinforced
SSHM.

1. Introduction
Recently, Steel Slag Hydrated Matrix (SSHM) has
been developed as a construction material for
reducing environmental problems. Its main
ingredients are pre-treatment slag and blast-furnace
slag powder which are by-products of steel making
process. In SSHM, the corresponding substitutematerial for cement is the mix of blast-furnace slag
powder, fly ash and slaked lime while thecorresponding substitute material for the fine andcoarse aggregate is pre-treatment slag. The processused in manufacturing SSHM (Fig 1.1) is the same asthat for normal concrete.
Steel making slag
(Pre-treatment slag)
Blast furnace slag
fine powder
Fly ash2. Outline of experiments
In this section, the outline of the experiments for thisstudy is presented.

Aging
Mix, Set, Cure

2.1 Materials used
In this research, two kinds of pre-treatment slag (K &F) were used for SSHM. The designation names camefrom the place of origin. Furthermore, in this study,
SSHM was compared with ordinary concrete. Table
2.1 shows the experimental cases and Table 2.2 shows
the specified mix proportions of SSHM and normal
concrete used.

SSHM

Water

Slaked lime

Fig1.1 Manufacturing process of SSHM
The application of this material is still limited
but it possesses promising features such as ①
manufacturing is possible using only by-products, ②
good viscosity and high segregation resistance during
fresh conditions, ③ high density, and ④ it isenvironmental friendly. When SSHM is applied toreinforced concrete, there should be a concern about
the decrease in durability due to the corrosion of thereinforcement. This is because corrosive materials in
SSHM move in a different way and also, the pH of theminute pore solution and the strength are rather low
when compared to concrete. However, studies on the
durability behavior of SSHM against penetration of

2.2 Items for measuring the infiltration of substances(Cl-, O2, CO2) into SSHM
Firstly, for both SSHM and normal concrete, thefollowing initial measurements were done; ①
penetration depth and ② chloride ion content (JIS A
1154) which are needed to determined the resistanceto chloride ion penetration. All these measurements
were done after the specimens were subjected toaccelerated test by repeated cycles of drying(temperature 40℃, relative humidity 35%) for 60

25

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

3.1 Results of chloride ion penetration into SSHM
Fig 3.1 shows the results of the chloride ion
penetration test while Fig 3.2 shows the results of thetest for total chloride ion contents. The calculated
chloride ion diffusion coefficients from the migration
test are shown in Fig 3.3.

hours and wetting (soaking in 3% NaCl solution) for
24 hours. Moreover, to compare SSHM with normal
concrete, the chloride ion diffusion coefficient was
examinedusing
themigration
test(JSCE-G571-2003).
Second, to compare the rate of ingress ofoxygen in SSHM to that of normal concrete, thelimiting current density was measured and theoxygen permeability was calculated.
Lastly, to compare the CO2 resistance of
SSHM with normal concrete, the carbonation depths(JIS A 1153) were measured.
Table 2.1 Experimental cases
Binder

SSHM(Sign : K40)

Blast furnace slag
fine powder
Fly ashSlaked lime

Normal concrete(Sign : NK40)
Normal concrete(Sign : NF40)

40

Gravels/a(%)

Pre-treatment
Slag (K)45

Pre-treatment
Slag (F)

67.2

Sand

15

10

5

0
0

Sand

Gravel67.2

Casting & Curing conditions① temp of 20℃ curing in water for 4 weeks② temp of 30℃ curing in water for 4 weeks③ temp of 20℃ curing in air for 3weeks

Table 2.2 Specified mix proportions of SSHM and normal concrete
Steel Slag Hydrated Matrix (Sign : K40)
Specified mix proportion (kg/m3)
40%s/aPre-treatment
slag
Blast furnace
Fly ashWater
Slaked lime（K）slag fine powder
5-20mm
53
297
100
182
1068
W/P

45%
Pre-treatment
slag
Admixture
（K）
0-5mm
874
4.046

Normal concrete (Sign : NF40)
Specified mix proportion (kg/m3)
W/C
40%
s/a
GravelOriginal Portland cement
Water
5-20mm
450
182
609

50

55

60

22
20
18
16
14
12
10
8
6
4
2
0 K40

K40
2～3cm

F40
0～1cm

F40
2～3cm

NK40
0～1cm

NK40
2～3cm

NF40
0～1cm

NF40
2～3cm

Distance from exposure surface

6
Chloride ion diffusion coefficient-82(×10 cm /sec)

Admixture
4.046

51 week
4 weeks
8 weeks
13 weeks

4
3
2
1
0

67.2%
Sand
Admixture
0-5mm
1247

15 20 25 30 35 40 45
Aacceleration period (cycles)

Fig 3.2 The measurement result of total chloride ions content

45%
Sand
0-5mm
874

10

0～1cm

Steel Slag Hydrated Matrix (Sign : F40)
Specified mix proportion (kg/m3)
W/P
40%
s/a
67.2%
Pre-treatment slag
Blast furnace
Fly ashWater（F）
Admixture
Slaked lime
slag fine powder
0-20mm
53
297
100
182
1867
6.391
Normal concrete (Sign : NK40)
Specified mix proportion (kg/m3)
W/C
40%
s/a
GravelOriginal Portland cement
Water
5-20mm
450
182
1068

5

Fig 3.1 The measurement result of the chloride ion penetration depth
45

OPC

K40
F40
NK40
NF40

20

Total chloride ions content(kg/m3)

SSHM(Sign : F40)

W/B (%)

Chloride ion penetration depth (mm)

Case

25

K40

F40

NK40

NF40

Fig 3.3 The calculation result of the chloride ion diffusion coefficient

In each of these examination results, the ingress ofthe chloride ion for F and NF in which the
sand-aggregate ratio is high (67.20%) was large.
However for the K and NK in which a lot of coarse
aggregates existed and only little amount of mortar
was present, the ingress was suppressed. Because it isconsidered that the chloride ion diffusivity ofpre-treatment slag is lower than that of mortar.
The chloride ion content of SSHM is generally
higher than concrete when measured near the surface(0-1cm) as influenced by the presence of fine particlesof blast-furnace slag powder. However, SSHM
exhibited improvement in water-tightness and longterm strength at its inner portions (2-3cm). Also, the

5.395

2.3 Items for measurement of corrosion of steel bar in
SSHM
In order to examine the corrosion resistance of steel
bar in SSHM, the anodic polarization curve, half cellpotential and polarization resistance were measured
using electrochemical techniques. Also, the corrosion
area (JCI-SC1) and mass loss in the steel bars
(JCI-SC1) were determined using
destructive
methods.
3. Examination results on the penetration of
aggressive substances (Cl-, O2, CO2) into SSHM

26

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

nondestructive test
Fig4.1 shows the measurement result of the anodic
polarization curve of SSHM (K40) and normal
concrete (NK40) after 55 cycles of acceleration test.
The state of passivity was judged by the grade ofpassivity according to past research [3] in terms of theanodic polarization curve. The higher grade indicatesthe better passivity.

chloride ion content at the interior of SSHM was lessthan that of concrete.
3.2 Results of oxygen permeability test in SSHM
Fig 3.4 shows the calculation result of the oxygen
permeability test while undergoing the accelerated
exposure.
5 cycles55 cycles1

5

Grade24

0.5

3Potential (V)Amount of oxygen permeability
×10 -11(mol/cm2/sec)

6

2
1
0

Grade20Grade3

K40
3cm

K40
5cm

F40
3cm

F40
5cm

NK40
3cm

NK40
5cm

NF40
3cm

NF40
5cm

K40 distance from exposure surface 3cm
K40 distance from exposure surface 5cm
NK40 distance from exposure surface 3cm
NK40 distance from exposure surface 5cm

Distance from exposure surface-1

Fig 3.4 The calculation result of the amount of the oxygen permeability

0.1

The values of oxygen permeability calculated after 55
cycles are much lower compared to those obtained
after 5 cycles. This is caused by the decrease in the
micro voids inside the SSHM as time passed by. This
improved behavior in oxygen permeability can
influence the corrosion of steel bar in SSHM as will be
discussed in the following sections.

100

Fig 4.1 The measurement result of the anodic polarization curve80
K40

Carbonation depth (mm)1
10Current density (μA/cm2)

The results indicated a better passive state existing in
SSHM when compared with concrete. This can be
explained by the following ① control of cathodereaction in steel bar by consumption of dissolved
oxygen that infiltrated from surface of SSHM by the
corrosion reaction of pre-treatment slag, ② control ofdissolution reaction in steel bars by pre-treatment
slag binding with chloride ions, ③ improvement ofsteel bars protection performance by demonstrating
good viscosity and good material separation resistanceof fresh SSHM.
On the other hand, an accurate value was notdetermined through the half cell potential andpolarization resistance measurements for SSHM.
This is because of the early corrosion of pre-treatment
slag itself that occurred at the surface part which hadsome effects in the measurements.

3.3 Results of carbonation test in SSHM
Fig 3.5 shows the depth of carbonation versus elapsed
time during the acceleration test. Resistance against
carbonation for SSHM was usually lower than that ofconcrete. This phenomenon is due to the consumptionof calcium hydroxide (supplied by slaked lime) by theblast-furnace slag powder and fly ash during the
pozzolanic reaction. Therefore, it is necessary to takecare of carbonation of SSHM especially inatmospheric zone.
70

Grade3

-0.5

F40
NK40

60

4.2 Results of steel bar corrosion by physicaldestructive test
Fig 4.2 shows the mass reduction ratio of steel bar by
the acceleration test after 25 cycles has passed.

NF40
50
40
30
20

Mass reduction ratio (%)

10
0
0

5

10
Acceleration period （√day）15

20

Fig 3.5 The measurement result of carbonation depth
4. Examination results of resistance to corrosion of steel
in SSHM

1 K1 : K40 F1 : F40
0.9 K2 : K40 casting & curing 30℃
F2 : F40 casting & curing 30℃
0.8 K3 : K40 curing in air condition
0.7 F3 : F40 curing in air condition
NK1 : NK40 NF1 : NF40
0.6 NK2 : NK40 casting & curing 30℃
NF2 : NF40 casting & curing 30℃
0.5 NK3 : NK40 curing in air condition
0.4 NF3 : NF40 curing in air condition
0.3
0.2
0.1
0 K1 K1 F1 F1 K2 K2 F2 F2 K3
3

4.1 Results of steel bar corrosion by electrochemical

5

3

5

3

5

3

5

3

K3
5

F3
3

F3 NK1 NK1 NF1 NF1 NK2 NK2 NF2 NF2 NK3 NK3 NF3 NF3
5
3
5
3
5
3
5
3
5
3
5
3
5

Distance from exposure surface (cm)Fig 4.2 The calculation result of the mass reduction rate27

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

The specimen used for this experiment is arectangular prism with horizontal steel bars. The
specimens were dosed with 10kg/m3 chloride ion
before mixing. This was done to further accelerate theattack of chloride ion and shorten the infiltration time.
The results show that after 25 cycles, the progress ofcorrosion was well controlled in SSHM specimens.
The reason for this is similar to those alreadyexplained above.

with reinforcing steel bar under marine environment,
all examinations and tests in the future should be
performed on specimens exposed to actual marine
environments for longer periods of time. This canallow for a more quantitative evaluation of the longterm durability of SSHM compared with existing
concrete structures. And, also it is necessary to study
the durability using other deterioration mechanisms
besides chloride attack.

5. Durability evaluation of SSHM as material withsteel reinforcement under marine environment
In this chapter, the prediction of deterioration due tochloride attack when SSHM is used as material withreinforcing steel is discussed.
In order to predict the deterioration of steel
reinforced members due to chloride attack, it is basic
to determine the length of time for each deterioration
stage: the incubation period, the propagation period,
the acceleration period, and the deterioration period.
Therefore this research tried to determine the
incubation and propagation period using the results
obtained above. The incubation period was calculatedusing the chloride ion diffusion coefficient of thespecimens cured for four weeks (Fig. 3.3), thegeneralized equation of the Fick’s second law and thechloride concentration limit of 1.2kg/m3 around thesteel bar (JSCE). While the propagation period was
calculated by using results of mass reduction due tocorrosion, the equation from JCI-SC1 for mass loss in
mg/dm2/day and the corrosion mass limit of 10mg/cm2
which is generally used as the quantity of corroded
mass of steel at the occurrence of crack(JSCE). Table
5.1 shows the calculation result of the average ratio of
deterioration periods between SSHM and normal
concrete．

6. Conclusions
The conclusions of this research are shown as follows.
1. The resistance to chloride ion and oxygen
penetration of SSHM is equal or even better thannormal cement based materials having the same mix
proportion and curing condition under marine
environments. This result is clearer in the interior of
the specimen. On the other hand, it is necessary toconsider appropriate measures for the use of SSHM
against carbonation, because the carbonation depth of
SSHM is larger than normal concrete materials.
2. The state of passive film of the steel bars that
influences corrosion in SSHM under marine
environment is at least equal to or even better than
the normal concrete materials having the same mix
proportion and curing condition.
3. The duration of the incubation period and the
propagation period of corrosion of steel in SSHM
under chloride attack is over 3 times longer thannormal concrete having the same mix proportion andcuring condition. Therefore it can be selected for use
as a substitute material for existing reinforced
structures in regions where its supply is highly
available and under marine environment.
References[1] Coastal Development Institute of Technology,“Steel Slag Hydrated Matrix technical manual – Useof Effective Technology For Pre-treatment Slag -”,Coastal Development Institute of Technology, 2003.
(in Japanese)[2] The Japan port and harbor association, “Standards
and explanations for the technology of facilities inharbors”, The Japan port and harbor association,
1999. (in Japanese)[3] N. Otsuki, “Research on corrosion of reinforced
concrete under marine environment”, Tokyo Instituteof Technology degree thesis, 1986. (in Japanese)[4] K. Kawata, S. Sakashita, Y. Hamazaki, and K.
Sugimoto, “Corrosion prevention of reinforcementsteel in mortar through small metal piece additions”,
Concrete Journal, Vol.38, No.2, pp. 29～33, Feb. 2000.
(in Japanese)

Table 5.1 The comparison of deterioration periods between
SSHM and normal concrete

Incubation period
Propagation period

SSHM/Normal concrete
3.24
4.92

The results show that SSHM had longer
periods of incubation and propagation compared withconcrete. This shows that the use of SSHM as
substitute to concrete in steel-reinforced structures
under chloride attack is quite feasible and also theservice life of the structure possibly becomes longer.
Therefore SSHM can be a substitute material inexisting reinforced structures such as caissons,cellular blocks, quays and mooring-posts which
structures that are highly used in marine
environments.
To further examine the application of SSHM

28

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

A study on Gas -liquid two - phase flow patterns in rectangular micro -channels
Student Number: 05M18053

Name: Shintaro KOBAYASHI

Supervisor: Shiro YOSHIKAWA

矩形微小流路における気液二相流の流動状態に関する研究小林慎太郎

近年、微小流路を有するマイクロ化学プロセスにおける反応・分離に関する研究が盛んに行われている。流路
内の流動状態を制御することはこれらの装置の設計上、重要である。本研究では、矩形微小流路に気体と液体を同時に流し、その流動状態を観察し、分類を行った。流動状態に関係するパラメーターとして表面張力、粘度、流量、圧力を挙げ、流動状態を決定する重要な因子となる無次元数に基づいた流動状態図の作成を行った。
さらに、操作条件と流動状態の間の関係について新たな知見を得た。

1. Introduction

important dimensionless numbers for the transitions of the flow

There are many kinds of equipment, where gas and liquid

patterns and to get the flow pattern map of Gas-liquidcontact with each other, for instance a chemical reactor, atwo-phase flow using dimensionless numbers as coordinates.

cooling tower, ink-jet systems, spray systems and so on. In the

2. Experiments

equipment, transport phenomena between gas and liquid phase

Figure 1 shows the experimental apparatus schematically.

should be important. The phenomena should be discussed on

Nitrogen gas was fed by a gas cylinder. Flow rate was

the basis of the pattern of the Gas-liquid two-phase flow.

controlled by a regulator and a needle valve. The gas piping

Recently, various kinds of micro chemical process with micro

was made of stainless steel. On the other hand, liquid was fed

channel have been studied. In a micro channel reactor,

with a fixed quantity by micro-feeder. The water feeding pipe

temperature and residence time are easily controlled. As a

was made of PTFE. The liquid was deionized water or 10wt%

result, the selectivity of the reaction is expected to be high.

aqueous ethanol.

Because of the high selectivity, the micro chemical processes

Figures 2 and 3 show the designs of the test sections. They are

would be utilized in the pharmaceutical industry and themade of glass. The channel height of each test section was lesscosmetic industry even if the amount of the products might notthan mm order. The channel length of each test section was

be so large.

enough long from the inlet length, so it is thought that the flow

Gas-liquid two-phase flow reactors would be important also inwere fully developed flow. Test section was put horizontally, so

the micro chemical process and some devices have been

it is thought that the influence of gravity was able to bestudied. As described above, it is important to get information

ignored.

on the flow patterns in the Gas-liquid two-phase reactor in

In the experiment, at first, the test channel was filled withorder to estimate the performance of the device. There is little

liquid. After that, nitrogen gas was fed to the test channel. The

information for controlling the flow patterns with theexperiments were carried out under the condition of various gas

operational conditions though it would be indispensable for theand liquid flow rate. The operation patterns are summarized inprocess very useful. The purpose of this study is to discuss the

Table 1.

29

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

Gas regulator

C

2

Needle valve

4mm or 2mm

N

Micro feeder

0.5mm

A B

65mm
Flow meter

CameraGas cylinder

Figure 3. Test section type β.

Test section
VesselFigure 1. Schematics of the test apparatus.

Table 1. Operation patterns.
LiquidGas

Test section

inlet

inlet

Experiment 1

type αA

B

Experiment 2

type αB

A

Experiment 3

type βA

B

Experiment 4

type βB

A

Figure 4. Schematics of the test sections.

3. Flow patterns and transition
In Experiments 1, 2, and 4, it was difficult to get the stable

There should be many factors, which influence the flow

flow pattern. In Experiment 3, three flow patterns were

patterns. But, it is impossible to take all of the factors into

observed.

account in making the flow pattern map. The most important

Figure 4(a) shows the Intermittent flow. It seems that there areones were discussed in the following from three view points.

large gas plugs in the liquid phase. Figure 4(b) shows the

3.1 Velocity and intervals of plugs

Separated flow. It seems that gas-liquid interface seems smooth.

Figure 5 shows transition from Intermittent flow to Separated

Figure 4(c) shows the Annular flow. It seems that gas-liquid

flow. Thinking about such a state and assuming that theinterface becomes turbulent.

intervals of gas plugs Ll and the length of the plugs Lg are

BCconstant, the time interval of the plug

0.5mm

A

30°

following equation.

T =4mm or 2mm

65mm

T is expressed by the

Lg + Llup(1)

u P shows the velocity of the plug. The superficial velocity ofgas j g and that of liquid jl are defined as follows.

Figure 2. Test section type α.

jg =

30

Qg
A

, jl =Ql (2)A

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

Q g , Ql show the volumetric flow rate and A is the area of

Eq. (8) means that µ , especially viscosity of liquid µ l , has

the cross section of the channel. A volume of a plug VB issignificant effects for the transition of flow patterns.

expressed as follows.

Flow

VB = QgT (3)On the assumption that the shape of a plug is half spheroid, thexvolume of a plug is expressed as follows.

VB =

π 2
x L g (4)12

Lg Ll=L
Distance between gas plugs is shorter.

u B shows the relative velocity between the plugs andsuperficial velocity of the two-phase mixture. So,u p is
Ll=L-dLDistance between gas plugs is far shorter.

written by

u P = jg + jl + u B (5)From Eqs. (1),(3),(4) and (5), L l is expressed as follows.

Ll=0 (Lg=∞)

π 2
x Lg ( j g + j l + u B ) - Lg j g A(6)
Ll = 12
jg A

Figure 5. Transition from Intermittent flow to Separated flow.

Flow

It is necessary for the transition from Intermittent flow toySeparated flow to be the condition Ll = 0 . When Ll = 0 ,x

Eq. (6) becomes as follows.

b

πx 2
0=
( jg + jl + u B ) - j g A
12(7)
jl u B 121Û ++=A j g j g πx 2
Eq. (7) means that

a

μg
μl

gas

vg

liquidvl

τyx

Figure 6. Flow of gas and liquid between two flat plates.

3.3 Surface tension

jl
has significant effects for the transitionjg

Figure 7 shows the growth of a plug. Thinking about such astate,

according to the Laplace-Yang equation, Pin and Poutare the pressures inside and outside of a plug.

of flow patterns not to speak of j g , j l .

Pin - Pout =3.2 Shear stress

2σ
r

(9)

It seems that the shear stress on the gas-liquid interface

The pressure P may be written by the characteristic length L

between gas and liquid has effects for the transition of flow

and the characteristic velocity U .

patterns. Figure 6 shows the flow of gas and liquid between

Pin = Pg =two flat plates. Thinking about such a state, the shear
stressτ yx on the interface is written by the pressure gradient P ,µ gU g2

L, Pout = Pl =µ lU lL2(10)

From Eqs. (9) and (10), the radius r of a plug isthe viscosity µ and the thickness of gas phase and liquidr=phasea , b .

2σL2
(11)
µ gU g - µ lU lEq. (11) means that σ may have effects for the transition ofæ 1 b2 µg - a 2 µl ö÷ (8)τ yx = Pçç 2 bµ g + aµ l ÷èøflow patterns.

31Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

flow seems to depend on the inertia force of gas phase. Thertransition from Intermittent flow to Annular flow seems todepend on the inertia force of gas phase and liquid phase,respectively. In the size of this study, inertia force is still

A plug grows big.

dominant. It seems that in the channel which is smaller than ther+drchannel of this study, viscous force is dominant. When viscous
force is dominant, Ca might be better thanWe .

Figure 7. Growth of a plug.

Flow pattern
Intermittent(N2−Water)

4. Flow pattern map

Separated(N2−Water)

4.1 Dimensionless numbers

Annular(N2−Water)

It seems that pressure P , density ρ , viscosity µ ,Intermittent(N2−Ethanol)

100

Separated(N2−Ethanol)

velocity U , surface tension σ and a characteristic length of theAnnular(N2−Ethanol)

channel D are the important parameters for the transitions of

10

Re l [-]

flow patterns. With these parameters, the radius of a plug R is
described as follows.

1æ æ Pg D ööçç÷÷÷ , (We ), (Ca ),çσç÷ (11)
è
ø
R
= fç
÷
D
æ Ul ö÷
ç æ Pl D ö÷÷ç÷, (We ), (Ca ), çç÷÷çç è σ øUè g øø
è

0.1
10

-5

0.0001

0.001

0.01

0.1

1

10

100

We [-]

We =Figure 8. Dimensionless flow pattern map.

µU (12)
ρU 2 D ,Ca =σσ

5.Conclusions

We shows the Weber number, and Ca shows the Capillary

The important factors in the transition of the patterns of thenumber. From Eq. (12), following equation is derived. There is

Gas-liquid two-phase flow in a micro rectangular channel were

the following relationship among We , Ca and Re.

discussed.

ρU D
2

We
=
Ca

µU

σ

Though two-phase flow has various parameters, in this study,σ = ρ UD = Re (13)µthe flow pattern map was able to be made using Re l andWe as the coordinates. As a result, the effect of inertia force is

Re shows the Reynolds number. By dimensional analysis, it isclear that Re, We and Ca should be used in the flow patternmore important then viscous force in the experimental

map.

conditions.

4.2 Dimensionless flow pattern map

6.References

Figure 8 shows the flow pattern map. The abscissa is We[1] 岡本秀穂, 橋爪新太,住友化学 2001−Ⅱマイクロリア

and the ordinate is liquid phase Reynolds number Re l . As a

クタの現状と展開 −合成実験、分析、製造への新規アプ

result of some trials in making maps, it appeared that We was

ローチ−,(2001), 32 – 45

better than Ca in correlating the experimental results. This[2] 植田辰洋，気液二相流 −流れと熱伝達− 第２版，養result suggests that the relationship between surface tension

賢堂(1989)and inertia force is more important than the effect of viscous[3] Tohru Fukano and Akira Kariyasaki, Characteristics offorce in the range of the experimental in this study. The solid

gas-liquid two-phase flow in a capillary tube, Nuclear

lines in the map show the boundary of the area of each flow

Engineering and Design, Volume 141, Issues 1-2, 2 Junepattern. The transitions from Intermittent flow to Separated

1993, Pages 59-68

32

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

IN-SITU SOIL REMEDIATION BY SEDIMENTATION METHOD~HEAVY METAL IONS’ ADSORPTION CHARACTERISTICS ON
SOIL AND ADSORBENTS~
汚染された砂質地盤の原位置攪拌洗浄Student Number: 05M18076 Name: Ryoichi SUGA Supervisor: Hideki OHTA and Pipatpongsa THIRAPONG

菅良一

土壌汚染対策法が施行されて以降，様々な浄化工法が提案されてきたものの，精度とコストの両面を満たす工法は数少ない．そこで，本研究では安価ですむ原位置地盤浄化工法の提案を目的としており，著者の研究はその基礎となるものである．まず人工的に重金属汚染させた土を，水，そして界面活性剤を加えた水とで洗浄することで，重金属の脱着特性を調べた．次い
で，洗浄水中の重金属を活性炭とゼオライトで回収することで，それらへの重金属の吸着特性を調べた．その他，ＳＳＥ，ＢＥＴなどの化学的なアプローチにより，これらの重金属の吸着・脱着特性の実験結果を理論的に実証することに成功した．
Key Words : adsorption, desorption, in-situ, heavy metal, biosurfactant, activated carbon, zeolite

fine soil particles, Desorption and Collection of HM
ions (See Fig.1).

1. INTRODUCTION
In July 2004, Tokyo Ohta Ward, unacceptable level ofcontamination (lead, oil, and trichloroethylene) was
detected from the ground owned by one of the factories
which was dealing with automobile manufacturing.
Ohta Ward is well known for the large number of small
factories and therefore there is a possibility that
contaminants are also detected from other areas in Ohta
Ward, which is undesirable for both land owners andalso residents near the sites. However not much effort
has been made for prevention and remediation of soil
contamination so far. What is worse, it is said that someof the land owners are hiding contaminated land andleave it even when they sell it to others. This is simply
because the cost of remediation is unaffordable for
them.
Therefore the purpose of this study is to propose a new
affordable soil remediation method. To keep the costlow, the method should be performed on site and ourtarget will be contaminated ground under small
factories, especially the ones dealing with plating.
Since the major contaminants detected from plating
factories are heavy metal, our target contaminants inthis paper will be limited in four Heavy Metals (HMs)(Copper, Lead, Zinc and Cadmium).

Adsorption

Desorption
Collection

fine particle(removed)

in watercoarse particle
in water with reagent
Adsorbent

contaminant
Biosurfactant

Fig.1 Process of the Sedimentation Method(2) Homogeneous Adsorption
Since HM ions generally exist only in the areaaround 3 to 5 meter from the ground surface (SeeFig.2), HM concentrations in soil in the casing are nothomogeneous. All the soil particles in the casing
should attach to the HM ions homogeneously so that theadsorbed HM ions can be removed efficiently as
explained in the folloing sections. In order that agitate
the suspension in the casing for a while.

P

P

pump
direction of water flow

2. CHARACTERISTICS OF THE
SEDIMENTATION METHOD

container
casing pipe

(1) General Description
To make the remediation cost-effective, in-situ soil
remediation method has been proposed. This
sedimentation method can be divided into 3 parts;
homogeneous adsorption, Sedimentation, Removal of

5mkind of a hammer grab
hole for water flow
contaminant

1m

Fig.2 Overall View of the Apparatus

133

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)(3) Removal of Fine Particles
Only the top layer of the sediment composed of fineparticles holding high percentage of contaminants in
the hole is taken out. According to the results ofsieving and sedimentation analysis that the author has
conducted (Results are as shown in Fig. 3), fineparticles have far larger surface area per gram thancoarse particles. See a red arrow in Fig. 3, representing
diameter of 0.074mm which is the boundary diameter
between fine and coarse particles in geotechnical fields.
If removing only a part of the soil, namely 10 percent ofsoil in weight settled at the top of the sediment, it will
be a great saving of time and cost (See Fig. 4)

2. EXPERIMENTAL RESULTAND DISCUSSION(1) Artificial Contamination of the Soil Sample
The uncontaminated soil sample was put into 20 millmolar (mM) solution of nitrate salts; Cu(NO3)2,
Cd(NO3)4, Pb(NO3)2 and Zn(NO3)2. According toReuss and Johnson, the nitrate ion is not significantly
adsorbed by soil fractures and is highly mobile. Hence,the presence of NO3 anion does not influence on the
behavior of HM ions. The soil-nitrate suspension wasaged in a reciprocating shaker for 16 days (See Fig. 5).
Grain Distribution Curve
The final concentrations of each HM ion after 16days
Percentage of Surface Area
exposure in the reciprocating shakers are Cu=2225ppm,
Pb=5102ppm, Zn=1101ppm and Cd=1393ppm.

100
Weight Percent

matter by holding them into the pores. There is no
guarantee, however, that AC is able to collect all sortsof heavy metal, and therefore another type ofadosorbent, zeolite, is also adapted.

80
60
40
20
0
0

0.01 0.1

1

10

100 1000
45

Diameter (mm)40

groundgroundPercentage Sorptionoo

smallerwater

Fig. 3 Grain Distritution Curve

35
30

Cu

25

Zn

20

Cd
Pb

15
10
5

particlesize

05

larger

0

10

15

20

Time (days)

Fig. 5 Daily Changes in the percentatige of HM adsorption
contaminated soil

10% in weight
90% in surface area

(2) Removal of the Upper Layer
Agitate and mix the suspension composed of artificiallycontaminated soil and deionized water for several
minutes and later leave the measuring cylinder, and thenthe double soil layer can be observed; fine particles on
top and coarse particles at bottom (See Fig. 6).

Fig. 4 Removal of fine pariticles in top layer(5) Desorption of HM Ions
it seems to be impossible to clear the environmental
standard only with those processes. Hence, in thissection, our focus is put on understanding of the HM
desorption behavior from soil surface and in the next
section, those released ions are going to be collected byusing adsorbents. It is well known that biosurfactant ishighly effective in removing oil from the surface ofparticular material as one can see in an advertisement ofdetergent. To enhance the desorption, therefore, we
adopt biosurfactant.

Wash Water
Fine material of soil

Coarse material of soil

Fig. 6 Two layers of sediment after agitation
14000

mg/kg (fine particles)(6) Collection of HM ions
In this section, collect all the ions released into washwater in a proper way. There are several types ofadsorbents suitable for this case. The most common
adsorbent of all is “Activated Carbon (AC)”, which is
the material mainly made of carbon and has a lot ofpores on its surface. AC is capable of collecting object

12000
10000
8000
6000
4000
2000
0

Cu

Zn

Pb

Cd

Fig. 7 Concentration of HM ions in fine soil

234

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

Concentration in Wash Water(mg/liter)(3) Water Wash and Biosurfactant Wash
In order to detach the HMs from soil surface, two
methods were proposed. The first one is just simply add
water and agitate the suspension with sufficient air
applied by air pump. From now on, we call this process
and water added “Water Wash” and “Wash Water”,
respectively. The second method is add biosurfactant(in our experiment, saponin) and enhance the efficiency
of desorption.
The following graphs show theconcentration in wash water.
Water Wash
800

800

600

600

Cd
Pb Cd
400

400

Cu Zn Pb
200

200
0
0

Cu
5
10
Elapsed Time (min)

15

0
0

found out which part of procedure can be and must be
improved. For HM existing in coarse particles, all thefine particles cannot be removed with one time washing
because in the sedimentation process, the degree ofseparation of coarse and fine particles is not perfect andnot all fine particles settle on top of the sediment.
Therefore, if the water wash is repeated and fineparticles are taken out continuously for several times, it
is expected that fine particles holding higher
concentrations of HM ions in bottom layer transfer to
Biosurfactant Wash
top layer progressively, which will be removed anddisposed. One more thing, for ions in wash water, isCd
Pb Cd
that zeolite, which is the adsorbent we adapted after AC,
Cu Zn Pb turns out to be ineffective in collecting HM ions at least
Cu
Cu, Zn, Pb, and Cd. Hence, another effective adsorbent
must be found.
5
10
15
Balance Sheet
(m g)
Cu
Zn
Pb
Cd

Elapsed Time (min)Fig. 8 HM concentration in wash water

1. Even after biosurfactant added, there is no distinct
difference found in the concentrations.
2. The initial part of inclination of the graph inbiosurfactant wash is steeper than that of water
wash, which indicates that desorption process withbiosurfactant can reach the equilibrium condition
faster than without it

(mg)
Cu
Zn
Pb
Cd

Concentration in Wash Water(mg/liter)(4) Collection Process

Initial Condition
total (mg)
44.5
22.0
102.0
27.9

After Biosurfactant
Coarse (mg)
Wash Water (mg)
0.02
25.2
5.3
10.33
48.2
0.07
1.8
17.67
After addition of Activated Carbon

Fig. 10 Balance of Heavy Metal transfer
Before Addition of Adsorbent

800

3. SOIL CHARACTERISTICS

After addition Adsorbent
600

Cu

400

Zn

200

Pb

0
0

(1) Carbon Content
Not only the surface area, but the presence of carbon
content also highly influences on the adsorption anddesorption characteristic. It is clear that fine particles
have more ability for adsorption than coarse particles.

Cd5
10
Elapsed Time (min)

15

%
H
C
N

Fig. 9 HM concentration in wash water after adsorbent wash

The figure indicates each HM ions concentration after
AC wash and how much AC is effective in thecollection of HM ions. However, AC turned out to beeffective only in the collection of Cu and Pb and any
distinct changes in the concentrations of Zn and Cdions were not found. It is possible to say that AC iseffective only in Cu and Pb ions, which is consistent
with the results from other studies, while only 10 to
20 % of the Zn and Cd ions were collected. Even after
the addition of zeolite, the result did not change
significantly.

Fine Particles
0.73
1.9
0.18

Coarse Particles
0.18
0.08
0.04

μ gramFine Particles
Sample Weight
8781
Content
H
64
C
167
N
16

Coarse Particles
12072
21
10
5

Fig. 11 Carbon content(2) BET Analysis
One of the main features of this remediation method is
to take advantage of large surface area of fine particles.
Hence, there is a need to determine how larger the areaof fine particles per 1g than that of coarse ones. The
result might be surprising. 1 gram of fine particles iscorresponding to around 28 m2 of surface area, while(5) Balance of HM transfer
At this moment, it seems that these laboratory resultsthat are still higher than environmental standard may
not live up to people’s expectations, however, thisresearch has just recently launched and at least we335

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

around 5m2 is for coarse particles. It is important tokeep in mind that the surface area of fine particles canbe much larger if the sepearation is more perfect.

fine particles, accounted for 90% in surface area. If all
the fine particles holding high concentrations of HM
ions are removed, in the last result, around 90% of HM
ions in target area can be removed.

Size
BET Surface Area
2
Fine Particles
27.831m /g
2
Coarse Particles
4.964m /g2
All (Fresh)
8.007m /g2
All (Contaminated)
9.804m /g

2. Adsorption and Desorption characteristics of 4
HM were revealed as follows
Cu: mostly exists in soil particle meaning Cu has large
adsorptive ability. However, Cu in coarse particles
needs to be transferred to fine particles, otherwise mostof Cu keep remain in coarse particles that we do nottake out from the casing. Repeat the sedimentation
process until perfect separation is achieved. As for AC
wash, Cu ions were dramatically adsorbed to AC.
Zn: spread out averagely to soil and wash water. Any
significant result was not obtained from water,
biosurfactant, and AC wash, and therefore, perfect
separation should be achieved and anothereffective adsorbent is needed.
Pb: The same tendency with Cu was observed.
Cd: It seems that Cd ions are easily and readily released
from soil particles from the release efficiency
obtained from water and biosurfactant wash, but
actually it does not. Cd ions are apparently
released, but in fact Cd ions never attached to soil
surface and just stay in wash water all the time, that
is, Cd does not have adsorptive ability to soil
surface nor AC.

Fig. 12 The result of BET Analysis

3. It is discovered that large amount of carbon content,
which is deemed to be highly influential to HM
adsorption, was detected in fine particle, moreover, 1gram of fine particle possess 28m2 of surface area.
Consequently, these results indicate that focusing on
fine particles for remediation is making a point and evenfabulous idea. The result from SEM analysis wasconsistent with the adsorption and desorption
characteristics mentioned above, and confirms thereasonableness of the laboratory experiment’s results.

4000
3500
3000

Cu

Zn

Pb

Cd
Exchangeable Fraction

2500

Carbonate Fraction

5. REFERENCE

2000

Fe-Mn Oxide Fraction

1)

BUDIANTA
Wawan,SALIM
Chris:LABORATORY EXPERIMENT ON THE
REMEDIATION OF METAL-CONTAMINATED
SOIL BY IN-SITU SOIL WASHING, 2006
2) Gillman, Sumpter: CEC Determination by the BaCl2
compulsive Exchange Method, 1986
3) Kyung-Jin Hong: Application of Plant-Derived
Biosurfactant to HM Removal from Fly Ash and Soil,
2000
4) Tessier et al: Sequential Extraction Protocol for
Analysis of HM Speciation in Soils and Sediments,
1979

Organic Matter Fraction

1500

residual fraction

1000
500
0

Co
nt
am
W ina
t
at
er ed
W S
Bi as oil
o
W h1
as 5
Co
h
nt
15
am
W ina
te
at
d
er
W S
Bi as oil
o
W h1
Co
as 5
nt
h
am
15
W ina
t
at
er ed
S
W
Bi a oil
o sh
W 1
Co
as 5
h
nt
15
am
W ina
t
at
er ed
W S
Bi as oil
o
W h1
as 5
h
15

Heavy Metal Concentration in Soil (mg/kg)(3) SEM Analysis
Sequential Extraction Method (SEM) is usuallyconducted to evaluate the speciation of particulate
metals. As a result of SEM analysis one can observe
how large or small the bonding strength is, how mobile
or immobile trace metals in soil are and how likely thetrace metals influence on the environment.
Fig. 13 indicates that Cd ions can be easily released into
water (soluble) or just never be adsorbed to soil surface.
For Zn, besides exchangeable fraction, carbonate
fraction and Fe-Mn matter fraction was also detected
and accounted for 30 to 50% of all. The second
fraction, namely carbonate fraction is still regarded ashighly soluble. Hence, it can be said that Zn are still
relatively soluble. For Pb, although exchangeable
fraction is still dominant of all fractions, carbonate
fraction and Fe-Mn Oxide fraction also constitute large
share of all fraction. Fe-Mn Oxide is regarded as
immobile part compared to other two above fractions.
Therefore, Pb is partially soluble but some of ions
possibly stick to soil surface. For Cu, Fe-Mn Oxide
fraction is dominant and it accounts for around 70 to
80% of all. Soil bound to Fe-Mn Oxide and Organic
matter is hardly released into water, that is, immobile.
Therefore, it can be concluded that Cu is less soluble
and hardly released into either water wash orbiosurfactant wash.

Fig. 13 The result of SEM analysis

4. CONCLUSION
1. According to the results from sieving andsedimentation analysis, it was discovered that 10% of
all object soil in weight, which is mainly composed of436Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

THE PSYCHOLOGICAL EFFECTS OF BEHAVIORAL CHANGE
CAUSED BY PUBLIC WORKS ON PLACE ATTACHMENT
Student Number: 05M18082 Name: Haruna SUZUKI Supervisor: Satoshi FUJII

土木事業による人々の行動の変化が地域への感情に与える影響に関する研究
鈴木春菜
人々がもつ地域への愛着の存在は，地域活動への協力を促すなど地域計画に大きな影響を及ぼすことが知られている．しかしながら，地域愛着の醸成過程，特に土木施策によってその過程にどのような変化がもたらされるかについては，十分に検討されていない．本研究では，心理学実験データに基づき，交通行動による「風土」への接触の程度の変化が長期的に地域愛着に与える影響，及び消費行動が地域風土との関わりや地域への感情に及ぼす影響について検証した．

2.2 Hypothesis of this study

1 Introduction

To examine the view of Hagihara [3], I proposed
hypothesis1.
Hypothesis 1: When individuals’ encounters withenvironment and neighbors in their residential area
increase, their place attachments to residential area
are enhanced.
While Hagihara [3] examined the effect of travel
behavior on place attachment, I go on to consider about theeffect of consumer behavior. Based on the concept ofHagihara [3], I proposed two hypotheses about consumer
behavior and place attachment.
In daily consumer behavior,
Hypothesis2: According to consumer behavior,encounters with environments and neighbors aredifferent.
Hypothesis3: The more encounters with environments andneighbors individuals have, the higher their place
attachment is.
I implemented surveys for residents in Japanese local
cities to examine hypotheses.

Many studies revealed that those who have positive
attitude toward the place tend to have sense of
responsibility to the place or the regional community andactively engaged in regional activities [1]. “Place
attachment”, one of the major concepts of human emotion
to the place, is studied in various fields. They showed thatplace attachment is affected by many factors, Individuality
like age, sex and religion, surroundings like public safety,
Individual behavior like daily activity[2].
Because Public works typically change people’ssurroundings, place attachment may be indirectly
influenced by public works. In turn, place attachment
would affect people’s public opinion and public behavior
that may have substantial effects on the content of public
works. Thus, place attachment is essential for deep
understandings of consequences and determinants of public
works. In spite of such importance, place attachment has
not yet well studied in planning of public works. With thisrecognition, this study focuses on place attachment.
Especially, the effect of public works on building of place
attachment was empirically investigated.

2.3 Definition of Place attachment
Researches on place attachment have been in the various
disciplines like a human geography, environmental
psychology, and social psychology. There are also various
definition and model of place attachment along each object.
In this study, I follow the definition by Hagihara [3] that is‘an affective bond or link between people and specificplaces’.
For the definition of the ‘place’ which people haveattachment to, Hagihara [3] defined as ‘school district ofelementary school or junior high school where resident livein’. Because this study also focuses on the role of place
attachment as a determinant of region-supporting behavior,
I use this definition of ‘place’.

2 Objects
2.1 Previous study
Hagihara & Fujii [3] focused on degree of encounterswith neighbors and environments, and its difference by theusage of transportation mode. They hypothesized thattravel behavior, especially car use, affects the degree of
encounters with environment and neighbors, and place
attachment. Their result showed the possibility of negative
impact of car use to encounters with social andenvironmental climate. Although they also hypothesizedthat the degree of encounters with them affect to the degreeof place attachment, they could not show it. They said that
because building of place attachment takes a long time, it
was difficult to explain by short term experiment.

3 The process of place attachment
building

37

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)(3) Amount of regional resources stock
To measure the amount of regional resources stock, 15
questions: “Are there (Park, convenience store,
supermarket, river and pond, shrine and temple, arable
land, historical landscape, forest, tourist spot, family
restaurant, arcade game, community centre, shopping
avenue, rail station, large shopping centre) your place?”
were used. Their questions consisted of three-point ranging,“nothing”, “a little” and “multiple”.

3.1 Method
In this study, I implemented a questionnaire survey tofollow up to the survey conducted in 2004 by Hagihara [3]to collect panel data.
3.1.1 Sampling and Data collection
Investigation areas are Hamamatsu city, Shizuoka pref.
and Toyohashi city, Aichi pref. Questionnaire was sent andbacked by mail. After data matching of two surveys, I got
117 samples in panel data. The response rate was 25.7% at
first survey and 38.8% at second survey. Participantsconsisted of 51.3% men. Their average age was 57.2 years
and average of resident years was 28.3 years (at 2nd
survey).
3.1.2 Measures
Questionnaire is consisted by contents about place
attachment, encounter with environment and neighbors,
stock of regional resources and Individual characteristic.
(1) Place attachment
To measure the place attachment, 13 questions shows in
TABLE 1 were used in the survey questionnaire. These
questions are made by Oya & Haga [4] and Hagihara [3].
Hagihara [3] classified them to three scales, Place
attachment: emotion (PAE), Place attachment: sustainment
desire (PAS), and Place attachment: preference (PAP).
Their questions consisted of five-point Likert-type ranging
from ‘I don’t think so’ (1) to ‘I think so’ (5). Cronbach’s αof each scale was statistically sufficient level as shown in
TABLE1.
(2) Encounters with environment and neighbors
To measure the degree of encounters with environment
and neighbors, 5 questions shown in TABLE 1 were used.
These are taken from Hagihara [3]. Their questions
consisted of five-point Likert-type ranging from ‘I seldom
do so’ (1) to ‘I frequently do so’ (5).

3.2 Analysis and Results
3.2.1 Correlation analysis
I used correlation analysis to examine the relation
between the change of encounters and change of place
attachment, shown in TABLE 2. TABLE2 show thatchange of encounters and change of PAP weresignificantly positive correlated. But change of encountersand other two scale of place attachment (PAS, PAE) were
not collated. PAP, PAS and PAE were significantly
correlated with each other.
TABLE2: Coefficients of correlation
Encounters
Encounters
Emotion
(PAE)

ｒ
p
n
Sustainment r
Placedesirepattachment(PAS)nPreference r
(PAP)
p
n

Place attachment
PAE

PAS

PAP

−

-

-

-

.019
.857
94
.009
.930
96

-

-

-

*

.179
.082
95

*

.304
.001
108

*

*

.507
<.001
106

.234
.015
107

* p < . 100

TABLE 1: Scales of place attachment and Encounters
Encounters with environment and neighbors （α=.84）

Positive effectError covariance

Each question asks ‘In your travel behavior, you….’
Hear insects and birds’ singing
Touch the outdoor air
Greet neighborhoods
Talk with neighborhoods

EncounterchangeGet a sniff of smell of nature, such as earth (soil) and flora.

PAEchangePAPchangePASchangePlace attachment: Preference (PAP)（α=.89）
I think that the place is livable.
I have the place which I love in the place.
I feel good when I walk in the place.
I prefer the ambiance and character of the place.
I prefer the place.
I feel relaxed in the place.
Place attachment: Emotion (PAE)（α=.83）

FIGURE 1: Estimated Structural Equation Model
TABLE 3 : Parameter Estimates fromStructural Equation
Modeling (B:standarized parameter)Bt-value
Encounter change⇒PAP change
0.12
2.02*
PAP change⇒PAE change 0.50
6.06**PAP change⇒PAS change 0.42
2.55*
PAE change⇔PAS change 0.13
2.36*
NFI: 0.982 CFI:1.000 RMSEA:0.000 n=116

I think the place cherished (or important).
I have an attachment to the place.
I feel that I have my own place in the place.
I feel the place is mine.
I want to live in the place forever.

Place attachment: Sustainment desire (PAS)（α=.91）
There are things what I don’t want to change.
There are things what I would be sorry if I lose it.

* p < .100

38

* * p < .050

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

during daily shopping, and consumer behavior.
(1) Place attachment
The same 13 questions as one used in 3.1.2(1) were also
used as scales of place attachment.
(2) Consumer behavior
The questionnaire includes questions asking about 10
daily-use shops, that is, shop name, shop size (small, large,
extra large), shop type (shopping street, supermarket,
convenience store, department store, others), shop location(urban area, suburb, near to station, near to home),transportation mode to shop (walk, bicycle, car, Train andBus (PT).), and frequency to use the shop (times/week),
average expenditure at the shop, attachment to the shop
(four-point ranging) and amount of conversation with store
staff (four-point ranging). I made two types of indicator for
consumer behavior from responses of these questions thatis, “modal share of each travel mode among all travel for
shopping” and “shopping frequency at each shop type(times/week)”. I consider 8 categories as shop type asfollows; 4 types (shopping street (SST) , small
supermarket(SSU), small [other] shop(SO), anythingelse(OT)) ×2 distances (near to home/ not near to home).
Where [other] denotes shops except convenience store,
shopping street, department store, and supermarket,
anything else means shops except SST, SSU, SO.
(3) Encounters with environment and neighbors during
daily shopping
In this survey, I made 5 questions to measure encounterswith environment and neighbors during daily shopping.
After principle component analysis, they were separated
into 2 components, i.e.”personal communication during
shopping” and “cognition and buying behavior of local
products” shown in TABLE4. There were some pertinent
questions at encounters with environment and neighbors inquestion item of consumer behavior, like distance from theshop. So, these are not only question to show the degree of
encounter. Their questions consisted of five-point
Likert-type ranging from ‘I seldom do so’ (1) to ‘I
frequently do so’ (5). Cronbach’s α of each scale wasstatistically sufficient level as also shown in TABLE4.

3.2.2 Structural equation modeling
Following the correlations, I estimated a structure equation
model shown in FIGURE 1. I presumed that it takes moretime to build PAE and PAS than PAP, and there isunobserved common factors between PAE and PAS as
shown by dot-arrow in FIGURE 1. The result shown in
TABLE 3 indicates that all of estimates were statisticallysignificant. This result support Hypothesis 1.
3.2.3 Factors which promote the place attachment
building
The effect of encounter change on PAP that we
empirically found may be dependent on regional attributesof the area. In order to investigate dependency on theregional attributes, dummy variables were created todenote attributes of the area. I used regression analysis totest interaction between encounters change and eachdummy variable to test the dependency. As a result, only
multiple existences of shrine and temple
hadsignificant positive effect on the PAP change. It was also
shown that increase of encounter to environment andneighbors have significant effect on place attachment if
there were some shrines and temples in the area, but such
effect was not significant otherwise.

4 The effects of consumer behavior on
Place attachment
4.1 Method
In this study, I implemented questionnaire surveys to
examine the effect of consumer behavior on place
attachment to test Hypothesis.
4.1.1 Sampling and Data collection
Investigation areas were Takamatsu city, Kagawa pref.,Toyohashi city, Aichi pref. and Kagoshima city,
Kagoshima Pref. A questionnaire was sent by postal mailto those who were randomly selected from city directories.
I requested a person who went for daily grocery shopping
to answer. An answered questionnaire was also sent backby postal mail. I got 507 samples from 3 cities. The total
response rate was 33.8%. Participants consisted of 83.2%
women. Their average age was 51.0 years and average ofresident years was 25.7 years.
4.1.2 Measures
Questionnaire was consisted of questions about place
attachment, encounters with environment and neighbors

4.2 Analysis and Results
4.2.1 Structural equation modeling
Following the concept by Hagihara [3], I presumed that
consumer behavior has effect on encounters and place
attachment. Based on this assumption, I estimate astructural equations model as shown in FIGURE2.

TABLE 4: Scales of Encounters with environment andneighbors during shopping
Communication during shopping (CS) （α=.83）

Encounter

Greet neighborhoods
Talk with neighborhoods
Talk with store staffs

Cognition and buying trend of local products(CBL)（α=.70）

Consumer behavior
CBL

I care about where products comes fromI often buy local products.
Note: Each question asks ‘In your shopping, you….’

CS

FIGURE 2: Estimated causal relation

39

PAP
PAS
PAE

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

Positive effect
Negative effect

CBL

Place attachment：
Emotion（PAE）

Error covariance

Frequency at near SST
Frequency at near SO

Communication during
Shopping (CS)Place attachment:
Preference (PAP)

Walking shareFrequency at
not near SO

Frequency at not near OT

Place attachment:
Sustainment desire(PAS)

PT share

FIGURE 3: Estimated model of consumer behavior and place attachment
Positive effect
Negative effectregional shops as an approach to understand place
attachment to an area. I presumed that shop type (dummy
variable denoting shopping street, convenience store,
department store, small size supermarket and others andextra large supermarket and others), location of shop
(dummy variance indicating urban area, suburb, near astation and near home) and frequency of visit to shop haveeffect on communication at shop and attachment to shop,communication at shop has also positive effect on
attachment to shop. Structural equations model was again
used for investigating structure of relations among thesevariables. The result indicates paths shown in FIGURE 4
were significant. As can be seen in FIGURE 4,
communications at shop were relatively high in small
supermarket and others, shopping street and department
store and relatively low in shops located at suburb. It wasalso shown that the more communication individual had ina shop, the more attachment to the shop he/she has.

Frequency
to visit

Type: shopping street
Type: department store
Type: small(supermarket + others)

Communication
at shop

Attachment
to shop

Location: suburb
Type: convenience store

FIGURE 4: Estimated model of attachment to shop
As shown in this figure, I consider CBL as an Index of
consumer behavior. In the result, relations which were
statistically significant were shown in FIGURE 3. FIGURE
3 show that all of 3variables of place attachment are
positively affected by CS. This suggest that the morecommunication during shopping people have, the moreplace attachment toward residential place they have.
FIGUEW3 also show; the more people go to shopping
street or small shops near to home, the more shopping
communication they have. These results support the
Hypothesis 2 and Hypothesis 3.
4.2.2 Relation between place attachment andattachment to shops
CS used in the analysis shown in TABLE 4 was not
limited to CS inside their residential place. But one may
presume that only CS inside their residential place has
direct effect on place attachment. For test this presumption,
I calculated correlations of them, in both the case of shops
near home (n=1266) and shops not-near home (n=1296).
The result indicates that attachment to the shops near home
had higher correlation with place attachment (PAP)(r=.207) than attachment to the shops not-near home had(r=.096).
4.2.3 Factors which promote the attachment to theshops
The analysis reported in 4.2.2 implies that factors of
attachment to regional shops in an area would be important
determinants of place attachment to the area. With thisrecognition, I tried to understand factors of attachment to

5 Conclusions
This study examined the hypotheses about process ofplace attachment building (Hypothesis1) and that theconsumer behavior affects the place attachment which was
mediated by the encounters with the environment andneighbors (Hypothesis2,3). To examine these hypotheses,
two surveys for residents were implemented. The results
imply that it takes a long time for place attachment to be
developed and consumer behavior had influences on
encounters and place attachment, as hypothesized in thisstudy.

References[1] Hidalgo, M. & Hernandez, B.: Place attachment: conceptual andempirical questions, Journal of Environmental Psychology 21,
pp.273-281, 2001.
[2] Brown, G., Brown, B. & Perkins, D.: New housing as neighborhood
revitalization –place attachment and confidence among residents-,
Environmental and behavior, vol.36 No.6, pp.749-775, 2004.
[3] Hagihara, G & Fujii, S.: The psychological effects of car use on place
attachment, the Proceeding of infrastructure planning, vol.33, 285,
2005 (in Japanese).
[4] OYA, H. & HAGA, S.: The effect of preferred daily mode oftransportation on elderly residents’ feelings towards theirneighborhood, Rikkyo Psychological Research, Vol. 45, pp.1-9,
2003 (in Japanese).

40Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

A procedure of determining parameters by adapting to
the results calculated by MEAM’92
Student Number: 05M18099 Name: Kento Tokumaru

Supervisor: Kunio Takahashi

MEAM’92 の計算結果に適合させた二量体参照修正埋め込み原子法のパラメータ決定法
得丸 建人
本論では，MEAM’92 で使用されているスクリーニング関数の不連続性が起きる原因を明らかにしている. そこで,スクリーニング関数をはずして計算できる, 参照構造を二量体とした二量体埋め込み原子法(Dimer Reference Modified Embedded Atom Method; DR-MEAM)のパラメータ決定手法を提案している．その決定手法は，スクリーニング関数の影響が及ばない範囲で MEAM’92 の計算結果に DR-MEAM パラメータを適合させるものである.
method of determining parameter withuniformity for many elements to expand theapplicability of MEAM to non-bulk systems.

1. Introduction
It is hoped that Modified Embedded Atom
Method (MEAM) will be applicable forengineering field, because MEAM cancalculate efficiently for large scaled simulation
with including physical implication to someextent. Recently, the applicability is getting to
expand to non-bulk systems, such as surface,cluster, however there are some issues. At first,
the parameters have not been determined formany elements consistently since MEAM’92
that is the base of MEAM series, and arbitrary
property is included during parameterization.
These causes are attributed to that it has adifficulty to fill up reliable physical properties.
Regarding arbitrary property, Yuan has solved
to some degree by utilizing physical property
of non-bulk systems. Secondly, “screening
function” that has been introduced since
MEAM’92 is considered to have discontinuity.
We demonstrate the application of particularcluster structure by MEAM’92, and clarify
how the screening function has influenced on
the discontinuity. Furthermore, we propose the2. Screening function
Screening function has been introduced forsimplification of calculation. Screening factor
is determined geometrically. We assume thecircumstance in Fig.2.1 to discuss discontinuityof the screening function. In case that the atom
j moves towards atom i, screening factor Sijk=0.
Adversely, in case that the atom j moves away
outside, screening factor Sijk=1.
We investigate the influence of the screening
factor geometrically by calculating cluster, as
shown in Fig2.2. The 7 atoms at first layer are
fixed. The 3 atoms at second layer are rotated,
and fixed relative atomic position, and released
only in the z direction. The discontinuity of thescreening factor appears at the angle 12˚ to 13˚
as shown in Fig. 2.3. Then it is obvious that therelationship of distance between d12 and d28
changes suddenly, as shown in Fig2.4.

41

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

3.

j

of

Dimer

Reference

Modified Embedded Atom Method

Rij(DR-MEAM)θi

Theoryk
RikTotal energy is defined as followings in the
MEAM theory.

Fig.2.1 Schematic in an example case of

Etotal = ∑ Ei

discontinuity(1)iEi = Fi ( ρ i ) +

5

z

6

x

9
10

y

8

where Ei is the energy of atom i, Fi is theembedding function, Rij is the distance between

3

atoms i and j, and φij expresses the pairinteraction between atoms i and j.
In former MEAM versions, in order todetermine the pair interaction for type-i atoms,a reference structure is used. The reference
structure is usually selected as equilibrium
crystal structure of type- i atoms. By
introducing cutoff or screening functions, thecalculation of the pair interactions can belimited to the nearest neighbor atoms only. On
the other hand, discontinuity of screening
function exists, as described previously. In thiswork, we apply diatomic structure as areference structure. Therefore, It needs not to
introduce the cutoff or screening function.
The interaction term can be obtained from a
universal energy function proposed by Rose etal. as the following.

Screening Factor, Si,k2
Fig.2.2 Atomic position to demonstrate

1.0
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0.0

S 2,8,S 8,2
S1,8
S8,1
S8,3

S 3,8
S4,8,S 8,4S 5,8,S8,5 ,S6,8,S 8,6,S7,8 ,S8,7
0

5

10
15
20
Angle,θ (deg.)

25

30

5.0

Distance of Two Atoms, d (m)Fig.2.3 Change of screening factor by rotating(×× 10 -10 )

4.5

d58

1
0
Eiu ( R ) = Fi ( ρ i ( R )) + φij ( R )2and then

d68

4.0

d 48

3.5

d 78
d 38

3.0

(2)

4

1

7

1
∑ φij ( Rij )
2 j ( ≠i )(3)

φij ( R ) = 2 Eiu ( R ) − 2 Fi ( ρ i 0 ( R ))
where

d 18

(4)

ρ i 0 is the background ground electron
u

2.5
2.0

d 28
0

5

density of dimer with distance of R. Ei is thed12,17,d 89,d810

10 15 20 25
Angle,θ (deg.)energy per atom for the dimer reference
structure as a function of the distance R, given
by a universal energy function.

30

Fig.2.4 Change of atomic distance by rotating

42

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

4. Determination of DR-MEAM
parameters+ a1 (

Eeq.DR. − Eeq.92 2
Req.DR. − Req.92 2)tetrag. + a1 ()tetrag.
Eeq.92
Req.92+ a1 (

Eeq.DR. − Eeq.92 2
Req.DR. − Req.92 2)tetrahed.+ a1 ()tetrahed.
Eeq.92
Req.92+ a2 (

There are 13 parameters in DR-MEAM, Z0,(5)→ min.

E0, R0, αi, Ai, βi(0)~(3), ti(0)~(3). Among thoseparameters, Z0 is the number of nearest
neighbors at the dimer (=1). E0, R0, and αi canbe obtained from the properties of the diatomic
molecule. ti(0) is a unit, just as MEAM’92.
The rests of 8 parameters, Ai, βi(0)~(3), ti(1)~(3)are determined by utilizing some results
calculated by MEAM’92. Ai is scaling factor

where Eeq.92 and Req.92 stand respectively forminimum binding energy and equilibrium
distance calculated by MEAM’92. Eeq.DR. andReq.DR stand for the same results calculated by
DR-MEAM. Additionally, a1 and a2 areweighting factors of clusters and bulk
respectively. In order to maintain reliability inbulk, the value of a2 needs to be higher thanthat of a1. The value of a1 is 1, the value of a2 is
4 in this work.

for the embedding energy, and βi(0)~(3) are
exponential decay factors for atomic densities,
and ti(1)~(3) are weighting factors for the partial
electron densities. Each optimized parameter issought by fitting to minimum binding energy
and equilibrium distance calculated by
MEAM’92. Table.1 shows each atomic
structure utilized in this work. We assume thatthe atomic position of the clusters, where thediscontinuity of the screening function does notoccur, is reliable.
Objective function is defined as,

5. Results and Discussion
Optimized parameters are listed in Table2,
and the values are the most minimal in Eq. (5).
It cannot be helped that there are some othercandidates, because some of them could be
locally minimal. However there is highpossibility that the optimized parameter exists
in the candidates. If the procedure in the samemanner is implemented for other elements, thevalues of parameters can be determined for thesame number of elements as MEAM’92.

Fobj{Ai , βi(0)～(3) , ti(1)～(3) }= a1 (

Eeq.DR. − Eeq.92 2
Req.DR. − Req.92 2)bulk + a2 ()bulk
Eeq.92
Req.92

Eeq.DR.− Eeq.92 2
Req.DR − Req.92 2)tri. + a1 ()tri.
Eeq.92
Req.92

Table.1. Figure of each regular atomic structure used for optimization of DR-MEAM
parameters, and the values of equilibrium distance Req.92 (Å) and binding energy Eeq.92(eV) calculated by MEAM92.
Triangle
RegularAtomic
Structure

Tetragon

Tetrahedron

Bulk
Req.92.bulk

Req.92.tri.

Req.92.tetrag.

Req.92.tetrahed.

Req.92 (Å)
Eeq.92 (eV)

2.55
5.45

2.39
8.01

2.64
7.19

43

2.56
3.54

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

Uncertainties of non-bulk systems, such as
the surface properties, are one of the main
reasons of the present study for are-determination of MEAM parameters. We
investigated the applicability of the DR-MEAM
to the surface energies using optimized Cu
parameters. The Cu surface energies of thethree low-index surfaces, (100), (110), and(111), compared with the experimental results
and the results calculated by 2NN-MEAM, arelisted in Table 3. The surface energies
calculated by this work are lower thanexperimental data. The results in this work
indicate as valid values as those of
2NN-MEAM, considering the experimental
results are extrapolated from adjacent melting
point. DR-MEAM has possibility to provide agood solution for surface properties, and othernon-bulk systems.

6. Conclusion
A determination method of the DR-MEAM
parameter values is proposed. The method
utilizes minimum binding energy andequilibrium distance calculated by MEAM’92.
The method can determine the parameter values
for many elements consistently as MEAM’92.
The surface properties of low-index surface
energies calculated by using the optimized
DR-MEAM parameters indicate reasonable
results just as 2NN-MEAM. The DR-MEAM
has the possibility to expand its applicability to
non-bulk systems.

Table 2. One of the optimized parameters for the DR-MEAM.
Ai β i( 0)
αi
Ei0
Ri0
Cu 1.01 2.22 4.42 0.64 4.77

β i(1)β i( 2 )β i( 3)

4.35

5.25

5.25

ti( 0 )
1

ti(1)
1.09

ti( 2 )
1.09

ti( 3)
1.11

Table 3. Comparison of surface properties by calculation methods. The parameters calculated by
DR-MEAM are optimized values by this work, as shown in Table2.

MethodExperiment[Ref.1]

2NN-MEAM[Ref.7]
DR-MEAM

Surface Energy, Esurf (mJ/m2)
1770
E(100)E (110)E(111)

1382
1404

1451
1489

1185
1341

References

50,152(1997)
5. T. Yamagishi, K. Takahashi, and T.Onzawa,
Surf.Sci. 445, 18(2000)
6. X. Yuan, K. Takahashi, Y. Yin, and T.
Onzawa, Modelling, Simul. Mater. Sci. Eng.
11(2003)
7. B. Lee, J. Shim, and M. I. Baskes, Phys.
Rev, B68.144112 (2003)

1. W.R. Tyson and W.A. Miller, Surf. Sci. 62,267(1977)
2. M. I. Baskes, Phys. Rev. B46, 2727(1992)
3. M. I. Baskes, J. E. Angelo and C. L. Bisson,
Modeling Simul. Mater. Sci. Eng. 2, 505(1994)
4. M. I. Baskes, Mater. Chemist. and Phys.

44

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

LONG TERM PREDICTION OF STRENGTH DETERIORATION DUE TO CA
LEACHING FROM CEMENT-TREATED SOIL
Student Number: 05M18113

Name: Katsufumi Hashimoto

Supervisor: Nobuaki Otsuki

セメント系改良地盤からの Ca 溶脱に伴う強度低下に関する長期予測橋本

勝文

本論文では，(1)粘土種類及び周辺環境中の Cl‑及び H+ が Ca 溶脱に及ぼす影響を定量的に把握すること，(2)(1)の結果を用いて，改良地盤からの Ca 溶脱に伴う化学的及び物理的変質に関する数値解析的長期予測手法を構築し，その妥当性を検証すること，(3)(2)で構築した長期予測手法を用いて，セメント系改良地盤の Ca 溶脱による化学的及び物理的変質に関する耐久性予測を行い，長期に亘る改良地盤の部材レベルにおける耐久性評価を目的とした．

1 Introduction

durability of cement-treated soil against Caleaching by using the above method.

Cement-treated soil is a composite material whosestrength and durability is improved by mixing soil
with cement as a hardening agent. Cement-treatedsoil constructed by Deep Mixing Method is
considered as highly durable material. However,previous researches indicate that there is a highpossibility of strength deterioration occurring in
cement-treated soil due to long term Ca leachingfrom Ca hydrates [1]. A characteristic of cementtreated soil is that soil particle itself has the
pozzolanic reactivity, comparing with othercementitious materials, such as mortar or concrete.
However, the influence of pozzolanic reactivity ofsoil on the deterioration due to Ca leaching is still
not clear. Additionally as cement-treated soil isconstructed by adding the hardening agent into thesoft ground at the site, the surrounding
environment is soft ground containing Cl- and H+.
However, the influence of Cl- and H+ insurrounding environment on the deterioration due
to Ca leaching is also still not clear. The above
influence on Ca leaching from cement-treated soil
is not considered in the long term prediction
method using numerical analysis in previos studies[2]. In addition, from the economic of view,
prediction result of the deterioration due to Caleaching should be considered during the design
stage in order to reduce unnecessary cost of
construction such as amount of cement or member
size. From the above background, the objectives ofthis paper are set as follows. Objective 1 is toinvestigate on the influences of pozzolanic
reactivity of soil, Cl- and H+ in surrounding
environment on the deterioration due to Caleaching quantitatively. Objective 2 is to propose
and confirm the validity of the long term
prediction method considering the above
influences. Objective 3 is to evaluate the long term

2 Outline of Experiment

In this section, outline of experiment issummarized as follows. Especially, experimental
cases, outline of specimen and measurement items
are described respectively.
2.1 Experimental cases
Experimental cases based on objective (1) aresummarized in Table-1. Experimental cases basedon objective (2) are summarized in Table-2. The
types of soil used in this paper are Daikoku,
Tokuyama, Kawasaki, Ariake, Kaolin soil. The
types of cement used in this paper are ordinary
portland cement (OPC), low heat portland cement(LPC) and blast furnace slag cement (BFSC).
Table-1 Experimental cases (1)3Soil type w (%) Ca(OH)2 content (kg/m )DaikokuTokuyama
100
30, 60, 90
Kawasaki
AriakeTable-2 Experimental cases (2)Soil type Cement type w(%)
OPC
Daikoku

87
BFSC
LPC
OPC

Tokuyama

132
BFSC
LPC
OPC

Kawasaki

110
BFSC
LPC
OPC

Ariake

100
BFSC
LPC

45

Cement content
3(kg/m )
42
84
125
84
84
42
84
125
84
84
42
84
125
84
84
42
84
125
84
84

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

exposure.
g) Estimated compressive strength
Estimated compressive strength is measured usingneedle penetration test and the relationshipbetween compressive strength and the slope ofneedle penetration as the parameter of physicaldeterioration after 13 and 26-week exposure.

Table-3 The amount of dissolvable silica and amorphousmaterial in soil
Daikoku Tokuyama Kawaskai
Amorphous
minerals (%)
Dissolvablesilica (%)

Ariake

17.9

17.4

9.8

8.5

7.5

9.4

3.8

4.2

2500

3 Influences of Soil Type, Cl- and H+ on
Deterioration Due to Ca Leaching

Compressive strength(kPa)Daikoku
2000

Tokuyama

High reactive

Kawasaki1500

In this section, specimens are exposed to soil andsolution respectively to investigate the influencesof soil type, Cl- and H+ on deterioration due to Caleaching.

Ariake1000
500

Low reactive

0
0

30
60
90
Ca(OH)2 content(kg/m3)

3.1 Influences of the amount of dissolvable silica
and amorphous material in soil on the
pozzolanic reactivity
The amount of dissolvable silica and amorphousmaterial in soil is shown in Table-3. The
relationship between soil type and compressive
strength is shown in Fig-1. From this figure, it isconfirmed that the larger amount of dissolvablesilica and amorphous material in soil, the higher
pozzolanic reactivity of soil. Therefore, it can besaid that Daikoku and Tokuyama soil is highreactive soil, Kawasaki and Ariake soil is low
reactive soil in this paper. The relationshipbetween distance form exposure surface andcompressive strength after 26-week exposure isshown in Fig-2. From this figure, the insidestrength after deterioration due to Ca leaching is
the highest in case of using LPC. The relationshipbetween exposure period and the deterioration
depth after 26-week exposure is shown in Fig-3.
The relationship between exposure period and thedecreasing ratio of compressive strength at 5mmfrom exposure surface after 26-week is shown in
Fig-4. From these figures, in the case of using
higher reactive soil, it is confirmed that thedeterioration more depthless. On the other hand, inthe case of using higher reactive soil, it isconfirmed that the decreasing ratio of compressive
strength at 5mm from exposure surface is larger..

120

Fig-1 The relationship between soil type and compressive
strength

2. 2 Outline of specimen
The specimens are cured for 24hours under sealedcondition with room temperature of 20oC. The
specimens are exposed to the soil or solution for
24hours after casting.
2.3 Measurement itemsa) Initial porosity
The initial porosity is measured and calculated
after each curing period. The initial porosity is
calculated from mass in condition of submerged,
oven dried and saturated surface dried.
b) Initial content of Ca(OH)2
The initial content of Ca(OH)2 is measured bydifferential thermal analysis after each curing
period. Raising rate of temperature is [40oC/min]and platinum cell is used in N2 gas from 40oC to1000oC.
c) Initial compressive strength
Initial compressive strength is measured bycompression testing machine. Rate of loading is
3[mm/min].
d) Ca leaching ratio
Concentration of Ca in outer solution is measured
using ion chromatography, and Ca leaching ratio
was calculated. Outer solution is exchanged every
7 days for 26-week exposure.
e) Ca content
Ca content is measured using energy-dispersible
X-rays as the parameter of chemical deterioration
after 13 and 26-week exposure.
f) Ca/Si mole ratio
Ca/Si mole ratio of cement matrix is measured
using energy-dispersible X-rays as the parameterof chemical deterioration after 13 and 26-week

3.2 Influences of Cl- and H+ on Ca leaching
The relationship between distance form exposure
surface and compressive strength after 26-week
exposure to NaCl and H2SO4 solution is shown in
Fig-5. From these figures, the higher theconcentration of Cl- and H+ in outer environment,
the faster the progress of deterioration due to Caleaching. In addition, it is confirmed that crack
occurs at the exposure surface in case of exposing
to H2SO4 of pH3-4.

46

Daikoku
1500
1000

OPC

500

BFSC
LPC

84kg/m3
0

0

10

20

30

40

50

2000

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

100

Tokuyama

Decreasing ratio of strength (%)

2000

Compressive strength (kPa)Compressive strength (kPa)

Technical Report of International Development Engineering

1500
1000

OPC
BFSC

500

LPC

84kg/m3
0
0

10

20

30

40

50

2000

Kawasaki
1500
1000

OPC

500

BFSC

84kg/m3
0
0

10

LPC

20

30

40

50

Compressive strength (kPa)Compressive strength (kPa)Distance from exposure surface (mm) Distance from exposure surface (mm)
2000

BFSC

10

LPC

20

30

40

50

Distance from exposure surface (mm) Distance from exposure surface (mm)Deterioration depth (mm)DaikokuTokuyama
Kawasaki
Ariake20
15

Low reactive

50High reactive
0

5
10 15 20 25
Exposure period (week)30

Fig-3 The relationship between exposure period anddeterioration depth
Decreasing ratio of compressive
strength @5mm (%)

120

Low reactive

100
80

LPC, 84kg/m3

60

DaikokuTokuyama
Kawasaki
Ariake402000High reactive

5
10 15 20 25
Exposure period (week)

30

600
500
400

NaCl=0%
NaCl=1%
NaCl=3%
NaCl=5%

300
200
100
0
0

10

20

30

40

50

Compressive strength (kPa)Compressive strength (kPa)Fig-4 The relationship between exposure period anddecreasing ratio at 5mm from exposure surface
700

700
600
500
400

pH3-4 (H2SO4)
pH4-5 (H2SO4)
pH5-6 (H2SO4)
Distilled water

300
200
100
0

0

10

20

30

40

b
+30%

20
40
60
80
Decreasing ratio of Ca (%)

100

4.1 Models for numerical prediction methoda) Ion migration model
Ion migration model, which is proposed inprevious study, Nernst-Plamk equation considering
Debye-Hückel theory and electrical neutrality
condition, is used.
b) Ca adsorption model
Ca adsorption model, which is based on Freundlich
adsorption equation, is used.
c) Ca dissolution model
Ca dissolution model of Ca(OH)2 and C-S-H
obtained from solubility test is used to consider the
influence of Cl- and H+.
d) Fixation model of chloride ion
Fixation model of chloride ion considering Ca
leaching, which is based on fixation model ofchloride ion proposed in previous study [3], is used.
e) Deposition and dissolution model of CaSO4
Deposition and dissolution model of CaSO4, which
considers the increase and decrease of pore volume,is use d.
f) Deposition and dissolution model of CaCO3
Deposition and dissolution model of CaCO3,
which considers the increase and decrease of porevolume, is used.
g) Strength deterioration model
Strength deterioration model, which is based on
the relationship between the decrease of Ca content
and that of strength obtained from exposure test, isused as shown in Fig-6.

LPC, 84kg/m3

10a

In this section, numerical prediction method for Caleaching is proposed. In addition, the validity ofthe method is confirmed by comparing to the resultof field investigation.

Fig-2 The relationship between distance form exposure
surface and compressive strength
25

-30%

20

a: trend line passing through point (100%,100%)b: line whose slope is +30% of ac: line whose slope is -30% of a
Fig-6 Strength deterioration model

OPC

0
0

c

40

0

1000

84kg/m3

60

0Ariake

1500

500

80

50

Distance from exposure surface (mm) Distance from exposure surface (mm)Fig-5 The relationship between distance form exposure
surface and compressive strength after 26-week exposure to
NaCl and H2SO4 solution

4.2 Confirmation of the validity of numericalprediction method
The result of field investigation of cement-treatedsoil, which is in Daikoku pier, after 20years ofservice life is used to confirm the validity of

4 Proposal of Numerical Prediction
Method for Ca Leaching47

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

investigation after 20years of service life.
Therefore, the durability against Ca leaching ofcolumnar cement-treated soil for 100years isinvestigated by the numerical prediction methodfrom the following view points. (1)In case ofusing the same cement content as cement-treatedsoil at Daikoku pier (156[kg/m3]). (2)In case ofusing the minimum cement content to fulfill thedesign strength (112[kg/m3]). The prediction
results of case (1) and (2) are shown in Fig-8.
From these figures, it is confirmed that the depth
where the strength is less than 2[MPa] is 10[cm] incase (1) and 23[cm] in case (2). In addition, theapparent strength is 3.15[MPa] (>2. 0[MPa]) incase (1) and 0.58 [MPa] (<2[MPa]) in case (2)after 100years if it is presumed that the strength indeteriorated area is 0 (Fig-9). Therefore, when
cement-treated soil is constructed, it is confirmedthat it is necessary to use amount of cement
content considering both design strength anddeterioration due to Ca leaching.

numerical prediction method proposed in thispaper [4]. The investigation result and prediction
result are shown in Fig-7. From this figure, it isconfirmed that the prediction result shows a good
agreement with the field investigation result.
Compressive strength (MPa)

12
10
8
6
No.1
No.2
Result
+30%
-30%

4
2
0
0

1 2 3 4 5 6 7 8 9 10Distance form exposure surface (cm)108Compressive strength (MPa)Compressive strength (MPa)Fig-7 The investigation result and prediction result(1)156kg/m3

6
4
2
0
0

5

10

15

20

25

Distance form exposure surface (cm)108(2)112kg/m3

6
4

6 Conclusions

2

As the results of this paper, the following
conclusions are obtained. (1)The influences of
pozzolanic reactivity of soil, Cl- and H+ insurrounding environment on the deterioration due
to Ca leaching are clarified. (2)The numericalprediction method for the deterioration due to Caleaching is proposed and the validity is confirmed.
(3)It is confirmed that it is necessary to set theamountofcementcontent
considering
deterioration due to Ca leaching as well as design
strength when cement-treated soil is constructed.

0
0

5

10

15

20

25

Distance form exposure surface (cm)Fig-8 The prediction results of case (1) and (2)
R
L

r

qd

qu

r
R

2

qd: Apparent compressive strength after deterioration
qu: Compressive strength before deterioration
L: deteriorated depth
Fig-9 Outline of deteriorated columnar cement-treated soil

7 References[1]Masaaki Terashi, Hiroyuki Tanaka, Tsukasa
Mitsumoto, Sadakichi Honma, Terumi Ohashi, The
Report of Port and Airport Research Institute,Vol.22, No.1, pp.69-95, 1983.
[2]Takahiro Nishida, Nobuaki Otsuki, Masakazu
Ikeda, CHAMPAPHANH Bouavieng, Masaaki
Terashi, Proceedings of JSCE, Journal ofConstruction Management and Engineering,
No.784, VI-66, pp.205-206, 2005.
[3]Tsuyoshi Maruya, Somnuk Tangtermsirikul,
Yasunori Matsuoka, Simulation of chloride
movement in hardened concrete Proceedings of
JSCE, Concrete Engineering and Pavements,
No.442, V-16, pp81-90, 1992.
[4]Masaki Ikegami, Ichiba Takehiro, Ohishi Kanta,
Terashi Masaaki, Long-Term strength change of
cement treated soil at Daikoku Pier, Proceedings ofSoft Ground Engineering in Coastal Areas, pp.241246, 2002.

5 Long Term Prediction of Durability of
Cement-Treated Soil against Ca Leaching
In this section, durability against Ca leaching ispredicted for cement-treated soil at Daikoku pier inservice life for 100 years. As mentioned before,
from the economic of view, prediction result of thedeterioration due to Ca leaching should beconsidered during the design stage in order toreduce unnecessary cost of construction such as
amount of cement or member size. On the other
hand, there is a possibility that the deterioration
due to Ca leaching can not be ignored in case ofusing less amount of cement or member size.
Actually cement content of 156[kg/m3] is used at
Daikoku pier, and it is confirmed that the strength
in sound area (8.2[MPa]) is too higher than thedesign strength (2[MPa]) from the result of field48

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

Analysis of turbulent organized structuresat the urban outdoor scale model (COSMO)
Student Number: 05M18120

Name: Satoshi HIROOKA

Supervisor: Manabu KANDA

屋外スケールモデル（COSMO）における乱流組織構造の解析
廣岡 智
本論文の目的は屋外スケールモデル COSMO における乱流の観測を行い、乱流組織構造の存在を確認し、その定量的把握をすることである。具体的には水平方向に熱電対と風速計を大規模に配置し、乱流観測を行った。そして、乱流構造の定量的解析に向け、OSIM(Organized Structure Identification Method)と呼ぶ乱流構造の抽出手法を開発した。これは乱流変動の連続点を構造として認識し、そのスケール、形状等を明らかにする方法である。この新しい方法で、風速の増加に伴い、より細長く筋立った構造になること、水平方向に大きな構造は時間方向にも大きい構造であること、正の温度変動をもたらす構造
は負の温度変動をもたらす構造より筋立った構造であることが分かった。

1. Introduction
horizontal extent with 1.5m (=H) cubic roughness arranged
uniformly(Fig2.1). The plan area index is 0.25.

Asian Mega-cities have been faced on overpopulation.
Overpopulated
citieshave
significant
negativeenvironmental effects on atmosphere, such as heat island
phenomena and localized heavy rain. To overcome these
problems, it is necessary to clarify the exchange process
between the urban canopy and the atmosphere.
Up to the present, the existence of the turbulent
organized structure (TOS, Fig1.1) on air flow is clarified byvarious field’s observations and wind tunnel experiments. It
is recognized that TOS account for major heat transport andthus is a key to solve the atmospheric environmental issues.
In real cities, due to the social restriction, only a few verticalobservations have been done. Using wind tunnel
experiments, Adrian proposed the figure of TOS (Fig.1.1),
over a flat wall [1], however, TOS over real cities are still
unknown.
In order to break through these gaps, I construct thescale model city called COSMO (Comprehensive Outdoor
Scale Model experiment for urban climate, Fig1.2).
At the COSMO, I conduct the turbulent observation byusing the 70 thermocouples and 30 sonic anemometers, inorder to detect the spatial distribution of TOS. Moreover, for
the quantitative analysis, I develop a new analytical way
called OSIM (Organized Structure Identification Method). I
applied it to COSMO data to investigate the TOS.

2.2 Instruments
Instruments are set as shown in Fig.2.1. To observe theturbulent properties, 30 sonic anemometers (YOUNG model
81000) and 60 thermocouples are used. Sonic anemometers
are arranged on upstream side and downstream side. Allsonic anemometers’ observation height are 3m (=2H) andhorizontal resolution are 3m. 60 thermocouples arehorizontally arranged on downstream side only.
The observation height is 3m(2H) and horizontal
resolution are 0.75m (0.5H).
2.3 Calculation
To detect the turbulent motion, I calculate thefluctuation of instantaneous wind velocity and temperaturefrom their 1 minute averaged values. The data is only
analyzed when the horizontal wind direction is less than 5
degree from the prevailing wind direction as shown on
Fig.2.1.

Model
81000

0.75 m3m

Mean wind3m(2H)

2. Observation setup
Upstream

Downstream

2.1 Scale Model site
The size of the urban scale model is 50 x 100m

d
win
an
Me

Height:3m (2H)
Resolution (Thermocouples): 0.75m (0.5H)
Resolution (Sonic anemometers): 3m (2H)
Sampling:
10 Hz

Fig1.1 TOS (by Adrian et al, 2000)

Fig2.1 Configuration of instruments

Fig1.2 COSMO

49

ThermocoupleTechnical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

44.25(m)(a)

044.25(m)(b)
044.25(m)(c)

044.25(m)(d)
044.25(m)(e)
0060(sec)Fig 2.2 Horizontal distributions((a): the fluctuation of stream wise velocity at upstream side, (b): the fluctuation of temperature at upstream side,(c): the fluctuation of stream wise velocity at downstream side, (d): the fluctuation of temperature at downstream side,(e): the fluctuation of temperature at downstream side by thermocouples)
2.4 Qualitative analysis
Fig. 2.2(a), (b), (c) and (d) show the horizontal
distribution (vertical-axis) and time series (horizontal-axis)of
the fluctuation of stream wise velocities on the upstream side,
temperature on the upstream side, stream wise velocities on
the downstream side and temperature on the downstreamside by sonic anemometers, respectively. And Fig 2.2(e) isthe horizontal distribution (vertical-axis) and time series(horizontal-axis) of the fluctuation of temperature bythermocouples.
In Fig.2.2 (c) and Fig.2.2 (d), the dotted circled areashows the area of relatively lower wind velocity (Fig.2.2(c))
and it corresponds to that of higher temperature (Fig.2.28(d)). It is due to the coherent motion of low speed and hightemperature air mass.
In Fig.2.2 (d) and Fig.2.2 (e), the horizontal
temperature distribution by thermocouples and that byanemometers are well matched.
Hereby, I consider that horizontal distribution of temperaturefluctuation by thermocouples at the downstream successfully
captures TOS.
In Fig.2.2 (b) and Fig.2.2 (d).The difference between
upstream observation data and downstream observation data
are clearly differed. At the upstream, the higher temperaturezone is relatively larger than its zone at downstream. If the
turbulent organized structure is not produced by cubicroughness at COSMO, these figures must be same.
Therefore, it is clear evidence that the downstream

For the quantitative analysis of temporal-spatial figureof TOS, OSIM has been developed. OSIM stands for
Organized Structure Identification Method. OSIM detects
the continuous grid points on time and space fields as
Turbulent Organized Structure.
As the first step, the fluctuation of temperature isdivided by standard deviations at each observation points(Dimensionless process). And I extract the grid points inwhich the temperature fluctuations are over the threshold(See 3.2). Secondly, these points are clustered if they arespatiotemporal continuous. As the final step, the cluster’s
characteristic values are output (See 3.3). Based on the
outputs value, the results are statistically processed.
3.2 Threshold set up<Concept of Threshold>I set threshold as a constant value. In addition, thethreshold value (α) is set positive and negative (±α).
Positive threshold value detect the low speed streaks(hereafter, I call it “the high temperature structure”), whereas,
negative threshold value detect the sink of the cooler air(hereafter, I call it “the low temperature structure”).
<Sensitivity analyses of α>
It seems that based on the threshold value, the result of
OSIM must be different, for example, the number of TOS
or the size of TOS detected by α might be different. In orderto decide the appropriate α value, I do two sensitivity
analyses of α. The first sensitivity analysis is focusing on thechange of the number of TOS per 1min and the second
sensitivity analysis is focusing on the change of temperature

3. New quantitative analysis (OSIM)
3.1 Algorism of OSIM

50

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

1

90
80
70
60
50
40
30
20
10
0

(a)

(b)

ratio of deviation squared
through all TOS(±α)

the number of TOS per 1 min.

Technical Report of International Development Engineering

(c)

0.8
0.6
0.4
0.2
0

0

1

2

3

4

0

1

2

α

3

4

α

Fig3.1 the sensitivity analysis of α in thenumber of TOS detected by α

Fig3.2 the sensitivity analysis of α
in the total energy occupation by TOS detected by αdeviation squared through TOS detected by ±α.
Firstly, Fig 3.1 is the results of the first sensitivity
analysis of α in the number of TOS detected by α through 81
hours. In Fig 3.1, horizontal axis is the range of α, andvertical axis is the number of TOS per 1 min. Here, I’d liketo divide this figure into three parts as shown. In zone (a),the number of TOS increases. This is due to the decouplingof TOS. And after reaching peaks (α=0.8-1.1, in zone (b)),
the number of TOS starts to decrease (in zone (c)). This isbecause too much strong threshold ignores some of TOS. I
consider that α should be set in peaks (α=0.8-1.1). In theseregions, it seems that TOS is properly divided and not somany TOS are ignored.
As the second sensitivity analysis, I do the sensitivity
analysis of α in deviation squared through all TOS detectedby ±α. Fig 3.2 shows the result of this sensitivity analysis.
The horizontal axis is same of Fig 3.1 and the vertical axis
shows the ratio of deviation squared through all TOS
detected by ±α (shown Eq.1). This value means that thesum of TOS detected by ±α assume how much heattransportation. Here if α is set as ±1.0, then, by all TOS
detected by ±α, energetically about 80% of heattransportation are expressed. Hereby, I set 1.0 as a threshold
and OSIM is operated.

( ratio of deviation squared through all TOS ( ±α ))all _ TOS

＝ ∑ (T ' )2all

∑ (T ' ) ・・・ ( Eq.1)
2

all _ TOS : area occupied by TOS ( ±α )all : total areasome values as follows.
ID number
All clusters have their own unique No. for identification,
called ID number.
Mean span wise widthThis is the average of the difference between maximum
value and minimum value of the span wise coordinate
which a TOS has.
Mean time duration
This is the difference between maximum value andminimum value of the time coordinate which a TOS has.
Mean time scale length
This is the value multiplied mean time duration by mean
wind velocity at downstream.
TOS areaThis is the value multiplied mean time duration by mean
span wise width.
Energy of TOS
This is the sum of temperature deviation squared througha TOS area.

3.3 Output values
For the determination of characteristic’s TOS, I define44.25(m)044.25(m)044.25(m)00

Fig 3.3 Examples of OSIM results(top: Raw data, middle: high temperature structures, bottom: low temperature structures
based on mean span wise width, TOS are colored. )51

60(sec)Technical Report of International Development Engineering

u=0-1(m/s)
u=1-2(m/s)
u=2-3(m/s)
u=5-6(m/s)

90
80

100
time scale lengh of TOS(m)100time scale length of TOS(m)TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

70
60
50
40
30
20
10

positive structure
negative structure

90
80
70
60
50
40
30
20
10

0

0
0

5

10

15

20

25

30

0

5

mean span wise width(m)1015mean span wise width(m)Fig4.1 mean span wise width
and time scale length

Fig4.2 high temperature structuresand low temperature structures

3.4 Output example
Let me show two examples of OSIM results. Fig 3.3(a)shows the raw horizontal distribution data. And Fig 3.3(b)shows the result of mean span wise width of hightemperature structures. Here, based on the mean span wise
width, the TOS is colored. As Fig 3.3(c), the result of mean
span wise width of low temperature structures is alsocolored.

and time scale length, and the pink dots show low
temperature structure’s relationship. For easy comparison,pink dots are just plotted at right side of blue dots. And therange of standard deviation of time scale length is shown asbars. This result means the high temperature structure must
be the low speed streak and it tends to be much longer timescale than low temperature structure. In other words, its
shape is streaky.

4. Result and discussion

5. Concluding remarks

In order to understand the behavior of TOS with windvelocity changes, 51days observation data is analyzed by
OSIM. Only ±5degree wind direction and unstable condition
data is used. The data is categorized by mean wind velocities(from 0-1(m/sec) to 5-6(m/sec) every 1(m/sec)). And
analysis follows these wind categories.
One of OSIM’s accents is to analysis with keeping
spatiotemporal information of TOS. So, all TOS iscategorized by their mean span wise width. Each category is
set every 0.75(m). Based on these span wise, the results are
statistically processed. Here, results are shown.
First of all, in Fig.4.1, the statistics results on each span
wise width categories are shown. Here, horizontal axis is themean span wise width and vertical axis is the average valueof the time scale length on each category. And the red line
shows the aspect ratio equal to 1.0, then, the TOS figure is
like a square. Four results are shown. In strong windcondition (u=5-6(m/s)), the dots are quite far from the redline. In contrast, in weak wind condition (u=0-1(m/s)), thendots are almost on the red line. This means, the strong wind
tends to produce much slender TOS. Moreover, until themean span wise width, 7m, the proportional relation is clear.
With the mean span wise width increased, the time scale
length of TOS also tends to become longer. Moreover, after
reaching to the mean span wise width (7m), the time scale
length has peaked. This means, there is the limitation of thetime scale length of TOS.
Finally, the high temperature structure tends to have
longer time scale length than the low temperature structure’s
time scale length. Fig.4.2 is the result analyzed by
5-6(m/sec) wind velocity data. The blue dots show the hightemperature structure’s relationship between mean span wise

Using OSIM, the spatiotemporal figure of TOS isclarified. Concluding remarks are follows;
1. Strong wind velocity produce slender TOS.
2. Mean span wise width and time scale length is in
proportion, until its mean span wise width (7m).
3. The time scale lengths have their limitation.
4. High temperature structure tends to have much longer
time scale length than low temperature structure have.
This means the high temperature structure’s shape isstreaky.
From now on, the observation instruments must be welldeveloped and much more accuracy observation will be done.
For such a situation, OSIM must be one of good ways toinvestigate TOS. Moreover, if I can define the quantificationof the complex shape, like meanderings, wheals and so on,
the method works better and better, and contribute to clarify
the behavior of TOS. And in future, if the quantification of
TOS is successfully done with various roughness
experiments such as changes of the plan area index or theheight of cubic blocks, environmental assessment will be
done.

6. Reference[1] Adrian, R. J., Meinhart, C. D., Tomkins, C. D. : Vortex
organization in the outer region of the turbulent boundary
layer, J. Fluid Mech., Vol. 422, pp1-54, 2000.

52

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

Efficiency of the Implicit Finite Element Method
on Geotechnical Engineering
Student Number: 05M18136, Name: FURUMURA, Takahiro, Supervisor: OHTA, Hideki

応力積分に陰的積分法を用いた弾塑性・静的有限要素法と、陽的積分を用いた弾塑性・静的有限要素法を、精度や計算コストの面から比較した。弾塑性構成モデルには修正カムクレイモデルを用い、通常の弾塑性変形解析のほかに、地盤中の間隙水の移動を考慮する擬似静的土／水連成解析についても検討した。

1.

Introduction

This paper discusses the applicability of the implicit
finite element (FE) method to geotechnical engineering.
The implicit FE method, proposed by Simo (e.g. [1]),has been considered to be more accurate and morestable than conventional explicit methods. Borja applied

Strain resolution

ε& p + ε& e = ε&Elastic constitutive modelσ& = Ce : ε& eConsistent requirement

f& = 0

Associated flow ruleε& p = γ& ∂f ∂σHardening rulepc = pc ε&vp MD

it to geomechanics with the modified Cam-Clay modelRate-formed stress-strain relationshipin [2] and [3], and the implicit FE method showed good

∂f
∂f


: Ce 
⊗
Ce :
 e
ep
∂
∂
σ
σ
σ& = C − : ε& = C : ε&
f
f
f
∂
∂
∂
e

:C :
+1: 
∂σ
∂σ
∂σ 


accuracy and stability.
Although the implicit FE methods have been
expected to be numerically accurate, this method
requires far more computational costs than the explicit

Increment-formed stress-strain relationshipmethod. Considering that the accuracy of the explicit∆σ = C|epn : ∆εmethod becomes more accurate as the increment
step-size decreases, the efficiency of the implicit FE
method is in question.

Fig.1 A formulation of the stiffness matrix for the explicit
methodSeveral researches have been

made on a comparison between the implicit FEM andthe explicit FE methods. For instance, Harewood et al.
and Sun et al. compared several implicit and explicit

Strain resolution∆ε p + ∆ε e = ∆εalgorithms on commercial FE software in [4] and [5].

Elastic constitutive model∆σ = D ( ∆ε e )However, the studies on geotechnical problems seem toConsistency requirement

f =0

be missing before now. In this paper, we compare these

Associated flow rule

∆ε p = ∆γ ∂f ∂σ

two methods on several practical examples about soil

Hardening rulepcn| +1 = pcn| exp(ε&vp MD)

structures.
Increment-formed stress-strain relationship2.

σ|n +1 = D ( ∆ε, σ|n )The implicit FE method on geotechnics

In elasto-plasticity, the stress-strain constitutive models
Fig.2 A formulation of the stiffness matrix for the implicit
methodare generally expressed with stress-rate and strain-rate.
Since numerical methods deal finite time-step size, this53

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

constitutive model must be integrated with respect tostress state, the stress update algorithm using thetime.

modified Cam-clay model gave a good stability andIn the explicit FE method, we firstly determine thevalid solutions on most cases.

rate-formed stiffness matrix, then integrate it with
explicit numerical integration as fig.1. The stiffness
matrix C

ep

3.

Efficiency improvement of the implicit FE
programis updated every step using the current

stress state, and the next stress state is calculated with it.

The standard Newton method would show a very

Since the forward Euler numerical integration is

rapid quadratic convergence, however, the calculationcommonly used, smaller increment step size is required

costs of the Jacobian matrix J itself and J-1 arefor reliable results. On the other hand, we can introduce

significant. To avoid the calculation of J and the inverse

analytical solutions or

numericalof J, various modification of the standard Newton

integration schemes if the implicit method is used. In

method have been proposed. A natural choice would bethis research, we used the analytical solution to integrate

the use of the modified Newton methods that use thethe elastic constitutive model and the Crank-Nicolson

same gradient matrix for some or all iterations. Thismethod for the associated flow rule.

results a large reduction of the CPU time concerning the

moreaccurate

The implicit FE method solves non-linear system of

matrixes, while this method requires far more number of

equilibrium equations, while the explicit FE method

iterations. Matthies and Straing showed that a moresolves linearized equilibrium equations. Therefore, thedrastic modification of the modified Newton method is

implicit method require iterative nonlinear solvers. The

possible ([7]), for instance, using elastic tangential

Newton method is commonly used. The gradient matrix

stiffness matrix instead of elasto-plastic consistentfor the Newton method, also called the consistent

tangential stiffness matrix. If we use the elastic

tangential matrix, can be obtained by solving thetangential matrix for over all calculations, the global

following linear system of equations.

tangential matrix K keeps constant. Therefore, we do  e ∂f ∂∆γ∂ f   ∂σ:+ ∆γ Ce : 2   := Ce
I +  C :
∂
σ
∂
σ
∂
σ
∂
ε

 

not need to solve the matrixes at each step and the LU

2[1]

factorization is required just once. Matthies et al. alsoThis matrix would be updated every iteration if theshowed that the difficulty would arise when we apply

standard Newton method is used.

the modified Newton methods including the line search

After solving the global matrix, we need tomethod to unloading problems [7].

determine the stress for updated displacement. The

A comparison between the standard Newton

updated stress can be obtained by solving the following

method and the modified Newton method is shown in

non-linear system of equations for the σn+1:

fig.3. The elastic stiffness matrix is used for the initial∂f∆σ n +1 = Ce :  ∆ε − ∆γ

′
∂
σ
+
|
1
n



tangential matrix of the modified Newton method.

[2]

While the standard Newton method showed a quadratic

It is not guaranteed that eq.[2] has a unique solution.

convergence, the convergence of the modified Newton

Even if eq.[2] has a unique solution, the Newton method

method is linear rather than quadratic. Considering themight fail to find it. Since the failure of the stress-update

implicit FE method requires high computational costs

algorithm lead the failure of over all calculation, somefor non-matrix calculation, such as the stress-update

techniques might be required according to the plasticity

algorithm, it can be say that the modified Newton

models. Despite these difficulties, as long as the initial

method is not efficient for elasto-plastic FE problems.

value of the Newton method was equal to the previous

Fig.4 shows how the tolerances affect the54

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

its convergence is quadratic. The appropriate magnitude+0

Normalized residual ||r|| / || r0||10

of matrix-solver tolerance seems to be 10-12.
4.

-5

10

Comparison with the explicit method

Fig.6 shows a 9m width surface foundation
subjected

standard Newton method
modified Newton method

-10

10

to

uniformly

loading.

Analyses

were

implemented with the implicit FE program and the
explicit FE program. Although we expected that lager

-15

10

5

10

15

20

25

30

increment step-size leaded less accurate results, no

Iteration

significant difference was observed between ∆P = 0.1

Fig.3 Convergence of the Newton-type methods

kN/m2 and ∆P = 1.0 kN/m2 in both the implicit FEM+0Normalized residual ||r|| / || r0||

10

-8

and the explicit FEM. Fig.7 indicates that the explicit

case 1: solver_tol = 10-12case 2: solver_tol = 10-16
case 3: solver_tol = 10-24case 4: solver_tol = 10

FEM might be slightly inaccurate at the step when some

10-5elements yield, however, the error is not very large.
From these results, we can say that the implicit FEM is

-10

10

completely not efficient for the problem on normallyconsolidated soil and has fewer advantages even in the-15

problem on over consolidated soil. The explicit FEM

10

5

10

15

20

25

Iteration

Case 1Case 2
Case 3
Case 4

tolerance of the
ILUBCG
10-8
10-12
10-16
10-24

would exhibit enough accuracy at a lower computational

tolerance of thestress-update algorithm
10-7
10-7
10-7
10-7

cost in this kind of practical problems.
The second example is an un-drained shear problem
on axial symmetric strain condition shown in fig.8. An

Fig.4 Normalized norm of the global residual vector

instable result was observed after q p over the criticalstate parameter M. From the result of the implicit FEM

convergence of the global equilibrium equations. The

are shown in fig.3.9, such instability was not observed

norm of residual vector normalized with the initial value

but the calculation was failed when q/p over the criticalis plotted. It seems that the tolerance of matrix-solverstate parameter. Cleary the stress state which satisfy q/p

has much effect on the convergence, while the accuracy> M cannot be accepted because the model assume thatof the stress-update algorithm dues not affects so much.

the soil material shall fail when q/p = M. It can be

The accuracy of the stress-update algorithm can reachconcluded that the explicit scheme gives some solutions

enough level even if the tolerance is very loose becauseeven in unrealistic situations and the solution is notalways reliable. In contrast, the implicit FEM wouldstop in such situation but it does not provide unrealistic

uniformly distributed load P

20m

fixed boundary ux = 0

fixed boundary ux = 0

solutions.

50m
fixed boundary ux = 0, uy = 0

Fig.6

Finite element mesh

55Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

However, even if we could reduce the number of02∆P = 1kN/m
2∆P = 0.1kN/m

settlement [m]-0.02

iterations, the computational cost of the implicit method-0.04

can not below that of the explicit method which need

-0.06

not use iterative methods. Considering that the explicit-0.08

method gave acceptable results even on lager step-sizes,

-0.10
0

2

4

6

8

it is doubtful that the implicit method is efficient in

10
2

terms of accuracy and computational costs.

uniformly distributed load [kN/m ]Fig.7.1 Settlement under the foundation (implicit)

Some issues related to the efficiency improvementof the implicit FE method were also discussed. Firstly,0∆P = 1kN/m2∆P = 0.1kN/m2

settlement [m]-0.02

the use of inexact Newton methods was discussed. The-0.04

line search algorithm required more than 30 iterations to-0.06

achieve the tolerance 10-5 while the standard Newton-0.08

method required 4. Considering the computational costs

-0.10
0

2

4

6

8

of non-solver processes, such as stress-update algorithm,

10
2

uniformly distributed load [kN/m ]are relatively high in implicit FE method. since inexact
Newtonmethods

reduces

the

CPU

time

for

1.2

1.2

matrix-solver but does not reduces the other processes,

1.0

1.0

it can be say that the standard Newton method should be

0.8

used. The appropriate selection of tolerances was also

Deviatoric stress

Deviatoric stressFig.7.2 Settlement under the foundation (explicit)

0.8
0.6
0.4
0.2

0.6

discussed and it is shown that the optimization of

0.4

tolerances can reduce the CPU time of analysis.

0.2

0

0.2

0.4

0.6

0.8

1.0

1.2

0

0.2

Effective mean stress

Fig.8

0.4

0.6

0.8

1.0

1.2

Effective mean stress6.

Un-drained shear

References[1] J.C. Simo and R.L. Taylor, “Consistent tangent operators for
rate-independent elastoplasticity”, Computer Methods in Applied

Mechanics and Engineering, vol.48, 101-118 (1985)

5.

Conclusion[2] R.I. Borja and S.R. Lee, Cam-Clay plasticity, part I: Implicit

A formulation of implicit elasto-plastic FE

integration of elasto-plastic constitutive relations, Computer Methods

in Applied Mechanics and Engineering, vol.78, 49-72, (1990)

method for geotechnical problems was shown and[3] R.I. Borja, Cam-Clay plasticity, part II: Implicit integration ofcompared with the explicit FE method. Although theelasto-plastic constitutive equation based on a nonlinear elastic stress

implicit FE methods had been expected to be morepredictor, Computer Methods in Applied Mechanics and Engineering,accurate when lager step-sizes are used, no significant

vol.78, 49-72, (1990)

difference was observed in our research. The implicit[4] F.J. Harewood and P.E. McHugh, Comparison of the implicit andexplicit finite element methods using crystal plasticity, Computational

method was slightly more accurate if the analysis

Materials Science, (2006)

includes yielding behavior, however, this advantage was[5] J.S. Sun, K.H. Lee, and H.P. Lee, Comparison of implicit andexplicit finite element methods for dynamic problems, Journal ofnot significant in practical analysis. In terms ofMaterials Processing Technology, 105, 110-118, (2000)

computational costs, the implicit method required far[6] A. Gens and D.M. Potts, Critical state models in computational

more CPU time relative to the explicit method. Since weGeomechanics, Engineering and Computers, 5, (1988)

used the standard Newton method and it worked well,[7] H. Matthies and G. Strang, “The solution of nonlinear finite

the convergence was quadratic. The required number of

element equations”. International Journal for Numerical Method inEngineering, vol.14, 1613-1626 (1979)

iterations was 3 or 4 for the tolerance we selected.

56

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

Synthesis of Zeolites from Lake Sludge
Student Number: 05M18142

Name: Wei QUSupervisor: Hirofumi HINODE

湖沼汚泥を原料としたゼオライトの合成
曲 巍
本研究では湖沼環境汚染の原因である湖沼汚泥の有効利用を目的とし、湖沼汚泥からアルカリ水
熱合成によってゼオライト結晶への転換の可能性を研究した。アルカリ濃度、反応温度、反応時間
および固液比の合成したゼオライトに対する影響について検討した。アルカリ水熱処理により、主に２つタイプのゼオライトが生成した(ゼオライト NaP1、Analcime-C)。合成したゼオライトの陽イオン交換容量(CEC)および重金属吸着実験を行い、吸着性能について検討した。その結果、天然
ゼオライトと同等の CEC 値を有することからゼオライトの用途である吸着剤また土壌改良剤としての利用の可能性を見出した。

1. Introduction

investigated.

Lake pollution is serious because lakes are closed

2. Experimental

water systems and contaminants can accumulate

2.1 Synthesis of zeolites

easily. Discharge from households and factory flow

(1)Materials

to lake. Contaminants in the discharge accumulate on

Lake sludges from Kasumigaura and Teganuma

the bottom of lake, and they become lake sludge.

were used in this study. After being crushed toLake sludge causes deterioration of water quality andparticle with size less than 150μm, the sludge wasproduces negative impacts on ecosystems anddried at 105℃ for 24h. The sieved sludge was stored

households using lake water. Therefore, it must be

in a plastic bottle at room temperature as raw

treated adequately.

material for this study.

Lake sludge is mainly treated by dredging method(2)Methodwhich remove sludge from the bottom of a lake

Lake sludge was mixed with sodium hydroxideusing a dredging ship or machine. This method

solution (NaOH) and was shaken for 24 h in a 23ml

generate large amount of lake sludge that still need to

Teflon reaction vessel. The vessel was then put

be treated. Some methods to treat dredged sludge areinside autoclave and incubated at a temperature of 90

landfill, incineration, but they have problems of theto 150℃ in an oven. After the reaction, the slurrylack of space, cost, energy efficiency. Previous study

was filtered and the solid phases were washedshowed that coal fly ash can be converted into

several times with distilled water to eliminate the[1]zeolites by hydrothermal treatment . Since the main

excess NaOH. The washed solid phase was air-driedcomponents of lake sludge are silica and alumina

and stored in the plastic bottle at room temperature

which is similar to the chemical components of coalas the products.

fly ash, it can be expected that lake sludge could also

2.2 Characterization

be used as raw materials for synthesizing zeolite .

The products obtained were characterized as

In this study, hydrothermal synthesis of zeolite

follows. The identification of crystalline materials infrom lake sludge in alkali solutions was tried, and thesludge and synthesized zeolitic products was carried

synthesis reaction conditions and the application ofout by X-ray diffraction (XRD). The specific surface

the zeolite synthesized from lake sludge werestructure was observed by scanning electron
157

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

microscopy (SEM). The surface area was measured

temperature of 120 to 150℃. In addition, zeoliteby using BET method.

NaX peak appeared on the diffractogram at 90℃.

2.3 Cation-exchange capacity(CEC) measurement
The cation-exchange capacity was determined byUU(f) 4.0M NaOH

U

QP

sodium acetate method[2]. The sample was saturated

UUwith NaOAc solution (pH 5.2), and then the Na+ wasU(e) 3.0M NaOH

extracted by NH4OAc solution. The amounts ofAreleased Na+ were measured using Inductively

A

P A

Coupled Plasma- Atomic Emission SpectrometryAP

P

equipment (ICP-AES).

P(d) 2.0M NaOH
P

A

A

A

A

A

(c) 1.0M NaOH

2.4 Heavy metal adsorption
The synthesized zeolites were evaluated on their(b) 0.5M NaOH
Q

capability to remove heavy metal in aqueous

Q
H D G

H

solutions. 50ml of solution containing different
5

concentrations of Pb2+ or Cd2+ were poured into a

10

15

20

N

25

30

H

M

35

40

Q

45

50

(a) sludge

55

60

[2θ]

glass bottle with 0.1 g sample and were shaken at

Figure 1. XRD pattern of the products obtained fromvarious NaOH concentrations (L/S=5ml/g, 120℃-24h)

150rpm for 24h. After equilibration period, solid andliquid were separated by centrifugation and filtration(0.45μm membrane filter). The concentration of Pb2+

A

or Cd2+ in solution were determined by ICP.

A
A

3.Results and Discussion
3.1 The synthesis of zeolites( XRD analysis)
The results from XRD analysis of obtained products(d) 150℃
A

with various NaOH concentrations are shown in

A

A

A

Figure 1. With NaOH concentration of 1M and below,
the zeolite NaP1(Si/Al=1.67)[1] was formed as shownPin the XRD patterns. Quartz phase was also present.

P

zeolite Analcime-C (Si/Al =2)

PP

At 2M NaOH concentration, it appeared that the[1][3](c) 120℃

PQ(b) 90℃

X

phase was formed.

(a) sludge

It is possible that zeolite type changes due to the
Si/Al ratio in the liquid phase.

5

At NaOH

10

15

20

25

30

35

40

45

50

55

60

[2θ]

concentration greater than 2M, zeolite NaP1 andFigure 2. XRD pattern of the products obtained from
Various reaction temperature(2M NaOH, L/S=5ml/g, 24h)

Analcime-C disappeared due to the formation of

Q: Quartz (SiO2)

Unnamed zeolite.

H: Halloysite-7A (Al2Si2O5(OH)4)

G: Greenalite-1M (Fe6Si4O10(OH)8) D: Dickite-2M1 (Al2Si2O5(OH)4)

Temperature is also an important factor affecting

M: Quintinite-3T (Mg4 Al2(OH)12CO3・3H2O)

the formation of zeolite (Figure 2). Formation of

P: Zeolite NaP1(Na6Al6Si10O32・12H2O)

zeolite NaP1 was observed at temperature of 90 to

U: Unnamed zeolite(Na6(AlSiO4)6・4H2O)

120℃, while zeolite Analcime-C was formed at

X: Zeolite NaX (Na2Al2Si2.5O9・6.2H2O)

A: Analcime-C(Na(AlSi2)O6・H2O)

2
58

N: Albite (NaAlSi3O8)

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

3.2 The synthesis of zeolites( Observations by SEM)
SEM

imagesofzeoliteNaP1obtainedbyhydrothermal treatment were shown in Figure 3. The
formation of zeolite NaP1 as a fine crystal with adiameter of approximately 5μm can be clearlyobserved.
SEM images of zeolite Analcime-C obtained byhydrothermal treatment were shown in Figure 4. At
the condition of 2M NaOH and 150℃, the formation

15kv ×500 50μm

of zeolite Analcime-C with a diameter of 10μm can be(a) 2M NaOH, L/S=4ml/g, 150℃-24h

clearly seen.

15kv ×2,00010μm(b) 2M NaOH, L/S=5ml/g, 150℃-24h

15kv ×5,000 5μm

Figure 4 SEM images of zeolite Analcime-C obtainedby hydrothermal treatment.

(a) 1M
NaOH,
L/S=6ml/g, 120℃-24h
1M
NaOH(SLA:1/6/12,24h-120℃)

3.3 The Cation-Exchange Capacity(CEC)
Figure 5 shows the value of CEC of products
obtained from different NaOH concentrations. The
CEC was greatly increased until product with 1M
NaOH, and then was decreased between 1M and 2M
NaOH, and were constant after 3M. Based on the
XRD results, the formation of zeolite NaP1 seemed
to be responsible for the increase of CEC, and the

15kv

×7,500 1μm(b) 1M
NaOH,(SLA:1/7/14,24h-120℃)
L/S=7ml/g, 120℃-24h
1M
NaOH

formation of zeolite analcime-C caused the decreaseof CEC. It is possible that changes in zeolite typelead to different CEC values. Figure 6 shows the

Figure 3. SEM images of zeolite NaP1 obtainedby hydrothermal treatment.

value of CEC of products obtained from different
temperatures. The CEC greatly decreased as thetemperature increased to 150℃. Based on XRD
359

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

results, it indicated that zeolite Analcime-C has

4. Conclusions

lower CEC value than zeolite NaP1.

Lake sludge was successfully converted into
zeolites. Zeolite NaP1 and Analcime-C was formed

200

CEC[meq/100g]

through the hydrothermal treatment in alkali
150

solutions. Zeolite NaP1 was formed at 1~2M of
NaOH solution and at 90~120℃. Analcime-C was

100

identified at 1~2M of NaOH solution and at50

treatment temperature between 120 and 150℃. The
cation exchange capacities of the products ranged

0
0

1

2

3

4

5

from 49 to 173 meq/100g. The zeolitized products

NaOH concentration[mol/L]

obtained

Figure 5. CEC of the products obtained from various
NaOH concentrations (L/S=5ml/g, 120℃-24h)

this

study

have

the

same

cation-exchange capacities as the commercial naturalzeolites (about 50 to 170 meq/100g). Obtained
zeolites also have adsorptive capability toward Pb2+

200

1M NaOH
2M NaOH

and Cd2+ ions.

150
CEC[meq/100g]

from

The zeolite synthesized from lake sludge in thisstudy could be applied for several applications such as

100

soil improvement or water purification.
50

References[1] X. Querol, N. Moreno, J. C. Umaña, A. Alastuey,

0
70

90

110

130

150

170

E. Hernández, A. López-Soler and F. Plana:’

Temperature[℃]

Synthesis of zeolites from coal fly ash: an

Figure 6. CEC of the products obtained from various
reaction temperature(L/S=5ml/g, 24h )overview ,’International Journal of Coal Geology,
Volume 50, Issues 1-4, May 2002, Pages 413-423

3.4 Heavy metal adsorption
Table 1 shows the adsorption capacity of Pb2+ andCd2+ forvarioussynthesis

conditions.

[2] United States Environmental Protection Agency,
http://www.epa.gov/.:Method9081,Cation-Excha

These

nge Capacity of Soils (Sodium Acetate)adsorption results demonstrated that the increase of[3] N. Shigemoto, S. Shirakami, S. Hirano and H.

BET surface area of treated lake sludge provided a

Hayashi,’ Preparation and characterization of

better adsorptive capability toward Pb and Cd

zeolites from coal fly ash. Nippon Kagaku Kaishi

removal.

5(1992),pp.484-492

Table 1. BET surface area, CEC, and adsorption capacity of Pb2+ and Cd2+ for various synthesis conditions.
Zeolite Product

SBET(m2/g)

CEC(meq/100g)

Xm(mg/g) for Pb

P1 (90℃)

93.522166

P2 (120℃ )59.407

P3 (150℃ )29.724

2+

2+

Xm(mg/g) for Cd

Major phase54.94539.841

NaP1+NaX(tr)173

42.19435.336

NaP1

58

35.54619.531

Analcime-C

460

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

Influence of creep phenomenon on manipulation by
shape control of rubber-elastic membrane
Student Number:05M18159 Name:CHANG Xuanxuan Supervisor:Kunio TAKAHASHI

ゴム弾性膜の形状制御によるマニピュレーションに及ぼすクリープ現象の影響
常軒軒

マニピュレーションを行う時に表面凝着力を調整できるなら簡単に操作が出来る。表面凝着力は、接触部分の等価曲率半径 R によって決められる。それ
で、プローブの先端を変形させて、操作を行う方法を学部論文時に提案した。
圧力によるゴム弾性膜の形状を制御する実験を行った。静電容量式変位計を用いて、弾性膜先端の変位を計測できた。しかし、圧力をかけると、クリープ現象が発生する。クリープ現象は、形状制御に対する影響を検討した。

1

Introduction

But, because impurities are contained a lot in rubber, it is difficult to obtain rubber with high purity.
So, the traceability is low. To solve it, we select 3
kinds of rubber, which are a natural rubber (250µm
thick), a natural rubber excluding protein (100µm
thick), and a nitrile rubber (100µm thick).

A tool such as a pair of tweezers that is designed touse with millimeter-sized or bigger objects, doesn’tfunction effectively in the micro-manipulation of
micro-scaled objects. In micro-manipulation, adhesion force becomes dominant rather than the gravitational force. In various manipulation methods, there
is a mechanical method for manipulation by using
the difference of adhesion forces between the object
and the tool, and the object and the substrate. When
we use the mechanical method, it is considered that
we can advantageously manipulate objects by changing adhesion force. According to JKR theory[1], adhesion force can be determined by the equivalent radius of curvature R and the work of adhesion ∆γ.
Based on this theory, we can control the adhesion
force in the micro-manipulation by controlling theequivalent radius of curvature. In earlier study, K.
J. Obata et al. have investigated the manipulation
using capillary force with a concave probe-tip[2]. In
these two methods, it can be considered that theshape at the contact part can influence the manipulation.
A thin elastic membrane seems easy to change its
shape (inflate or deflate). In this study, we focus on
how to control the shape of a thin elastic membranein order to use in the micro-manipulation by controlling the air pressure. An experiment was carried out
to study the influence of creep phenomenon on the
manipulation by shape control of the membrane.

2

Experiments

2.1

Sample preparation

Fig. 1: Probe

In the experiment, the rubber membrane is stuck
to one end of the cylindrical probe by bond, as shownin Fig.1. The membrane is then washed with theethanol for ten seconds, and by pure water using an
ultrasonic cleaner for three minutes.

2.2

The experiment system we designed to use in thisstudy, as shown in Fig.2, consists of two main parts.
The measuring part (Part A) is to measure the displacement of the center of the elastic membranecaused by the deformation, and the pressure-control
part (part B) is to control the pressure, and consequently control the deformation of the membrane.
Below, we are going into details of each part.

We selected rubber to be used as the material of themembrane, because rubber can have large deformation, and plastic deformation doesn’t happen easily.
161

Experimental system

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

to the glass ball and make the cantilever bend. The
corresponding displacement of the cantilever at distance l u(l) (also u(a)) will be measured by the displacement meter. The relation between F and u(l)can be express as,

F =

(1)

where I is the cross-sectional moment of inertia andE is the Young modulus. According to this, we canmeasure the position where the probe and the ballcome in contact. Therefore, the displacement of thecenter of the elastic membrane can be measured.
Since the force considered in this study is verysmall, it is necessary to design the measuring part tohave sufficiently high resolution to measure it. We
select a cantilever with the length l = 40mm, thewidth b = 2.65mm, and the thickness h = 0.15mm.
The resolution of the displacement meter is ∆u is
1 × 10−9 m. Although, the model of this study differsfrom those of the JKR theory, we attempt to roughly
estimate the adhesion force to be 47µN. Then, we
initially determine the resolution as 0.47µN, 1/100of that estimated by the JKR theory. When we
determine the distance a=30mm as the position to
stick the glass ball to satisfy the above-mentioned
condition. The force resolution of this system is
1.2 × 10−8 N.

Fig. 2: Experimental system

2.2.1

6EI
u(l)
a2 (3l − a)Measuring displacement

Fig.3 illustrates the configuration of the force sensor(part C in Fig.2). A capacitive displacement meter (Iwatsu ST-0536A) equipped with a probe (ST0536A) is used. This meter can measure the displacement of conductors, e.g. a cantilever for thisstudy, using a non-contact method, and outputs voltage with high precision in proportion to the displacement. A cantilever (SUS304) is set right below the probe of manipulation, on which a glass ball
(φ =2.1mm), to be used as the manipulated object,
is stuck at a distance a.
An elastic membrane (natural rubber), as the manipulator, is stuck to one end of the cylindrical probeof 4mm inner diameter, the other end of the manipulation probe is connected to the air-cylinder (partB) for which the shape of the membrane can be controlled by adding/reducing pressure. The membrane
is set below the glass ball and can be driven up-down
by a pulse motor controller (Suruga Seiki, Model
D100) with a distance resolution of 0.05µm.

2.2.2

Pressure control

To control the deformation (shape) of the elastic
membrane, we select an air cylinder (SMC, CDM30200) which the inner diameter is 30mm, and thestroke is 200mm. The air cylinder is connected to the
manipulation probe where the membrane is stuck to
via a pressure meter (KEYENCE, AP-C30) with aresolution of 0.1kPa, as shown in Fig.4.
Weight is loaded to the cylinder to compress or decompress the air, accordingly, the membrane will inflate or deflate. The pressure meter shows a positive
value for the compressed pressure, and a negativevalue for the decompressed pressure. The experiment
was carried out in atmosphere, hence the variationof pressure ranges from -0.1MPa to 0.1MPa.

2.3

The process of the experiment

The process of the experiment:
1. Drive the membrane upwards to the sample ball.
2. When they come in contact, the change of theforce will be measured.

Fig. 3: System of the cantilever

3. The probe’s displacement record, and this position are set up in 0.

By driving the membrane upward the glass ball,when they come in contact, a force F will be exerted

4. Drive the membrane down.
262

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

0-2 0 00kP a-4 0 k P a-5 5 k P a-6 0 k P a-6 5 k P a-7 0 k P a-8 0 k P a

F lex u r e ω ( µ m )-4 0 0-6 0 0-8 0 0-1 0 0 0-1 2 0 0-1 4 0 0

Fig. 4: Pressure control part

-1 6 0 0

0

100

200

300

400

500

T im e t (m in )5. Repeat step 1,2 after a time interval.

Fig. 5: Relationship between the pressure P and thestrain rate ε̇ (Natural rubber (250µm thick))

6. Record the probe’s displacement, and drive themembrane down.
7. Change the pressure, and repeat step 5,6.

the center of the membrane. The relation between
strain ε and flexure ω is (2).

By these processes, we can obtain the flexure of thecenter of the membrane.

1

3

3.1

−2
3R0
) −1
ε=(
4ω + 2R0

The creep phenomenon of
the deformation of the membrane

The result of the change of the strain is shown in
Fig.5.
Using the result, we can calculate the strain rateε̇. By the result of the experiments, we can knowthat pressure decreases, and the constant strain rate
decreases.
As shown in Fig.7, the constant strain rate ε̇ becomes almost 0, when pressure becomes below a certain value. The constant strain rate has the relation
when the material is metal[3].

Experiment result of deformation

The flexure ω of the center of each membrane wasmeasured. We use the probe of 4mm inner diameter
R0 , and the sample ball is a glass ball of 2.1mm diameter. According to the measurement result in Fig.5
The thinner the membrane, the larger the membrane
deforms.

3.2

ε̇ = kσ n

Influence of creep phenomenon(3)

σ is stress,ek, n are constants that depend at temperature. The constant strain rate approaches gradually 0, when the fitting is done by the least square
method like the metal. Because the metal’s n is 4～6, so we used n = 4, n = 6. It differs from the experiment result. So, we can be considered that the creep
phenomenon that is different from the metal occurs,
when the rubber membrane was deformed. And others material show the same tendency, as shown in
Fig.8.

In the experiment, we found that the membrane,
during a constant pressure is applied, keeps deforming with several speeds corresponding to the appliedpressure .e.g. Fig.5 shows that, when the elastic
membrane and the glass ball come in contact. The
membrane deforms at a constant strain rate whentime passes.
This can be considered as the creep
phenomenon[3][4].
The phenomenon was observed and examined throughout the experiment.
Because the thin membrane is used, we can suppose
the center of the membrane is uniform equibiaxial
tension[5]. Then, we use the natural rubber, so
we can suppose the rubber’s volume is constant[6].
By this assumption, we can calculate the strain at

4

Conclusion

It has been realized to be able to control the deformation of an elastic membrane that stuck on the
3

63

(2)

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)(× 1 0 -6 )

3

C

B

-40kP a-55kP a-60kP a-65kP a-70kP a-80kP a0.0 2.

S tr a in r a te εT h e c h a n g e o f str a in ∆ ε

2 .5
2
1 .5

A

1

0.0 1
0 .5
0

-8 0

-6 0

-4 0

-2 0

0

P r e ssu r e P (k P a )
0

100

200

300

400

Fig. 8: Relationship between the pressure P and theconstant strain rate ε̇, A: Natural rubber, B: Natural
rubber excluding protein, C: Nitrile rubber

T im e t (m in )Fig. 6: Relationship between the change of strain ∆εand time t (Natural rubber (250µm thick))

top of the manipulation probe by controlling pressure. It was confirmed that the creep phenomenon
occurs on the elastic membrane by pressure. By the
result of the experiments, the rubber membrane ofdifferent materials shows the same tendency that theconstant strain rate ε̇ becomes almost 0 when pressure becomes below a certain value. Because therubber membrane of different materials has the different characters, the values are also different. Whenthe deformation is controlled, the creep phenomenon
can be suppressed if the pressure is below a certain
value, which can be obtained from the experiment.

( × 1 0 -6 )1S tr a in r a te ε.

n = 6 , k = 4 .2 × 1 0 -1 8
n = 4 , k = 2 .4 × 1 0 -1 4

References[1] K. L. Johnson, K. Kendall, and A. D. Roberts:“ Surface energy and the contact of elastic solid, ”
Proceedings of the Royal Society of London,
A.324, pp.301-313, 1971

0 .5

0

-8 0

-6 0

-4 0

-2 0

[2] Shigeki Saito, Tomoyuki Motokado, Kenichi
J.Obata and Kunio Takahashi:“ Capillary force
with a concave probe-tip for micromanipulation,”
Applied Physics Letters 87,234103,2005

0

P r e ssu r e P ( k P a )[3] F.K.G. オドクヴィスト、J. ハルト 共著、村上澄男 訳:“ クリープ強さの理論 ”、培風館

Fig. 7: Relationship between the pressure P and
the constant strain rate ε̇ (Natural rubber (250µm
thick))[4] “ ゴム試験法 ”、日本ゴム協会[5] 渋谷寿一, 本間寛臣, 斎藤憲司“, 現代材料力学 ”, 朝
倉書店,2001[6] 久保亮五,“ ゴム弾性 [初版復刻版] ”, 裳華房,1996
464

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

Development of a Stored Channel Model for UWB
Link Level System Design and Evaluation
Student number: 05M18165

Name: Huynh Thi Thanh Trieu

Abstract
超広帯域（UWB）無線は低電力で高速なデータ伝送を実現できる技 術として近年注目を集めている．実環境下における UWB 通信の伝送特性評価および設計においては，ストア
ドチャネルシミュレーションが有効な方法として知られている．本研究は，まずストアドチャネルシミュレーションによる UWB 伝送特性評価の有効性を実機との比較により検証する．次に，伝搬路の特性を面的に網羅するストアドチャ
ネルモデルを構築し，伝搬損失や遅延広がりなど UWB 伝
送特性に影響する伝搬パラメタの分布を導出する。さらに，
インパルス無線型 UWB 伝送において用いられる Rake 合成法の性能を評価し，その最適設計について指針を与える．

Supervisor: Jun-ichi TAKADATx
UWB Signal
BPF
Generator
3.1-5.0 GHz
4 bit 10 GS/s DAC

Rx

LNA

BPF
3.1-5.0 GHz

Transmitter

Digital Sampling
Oscilloscope
12 bit 20 GS/s

ReceiverFig. 1: Block diagram of UWB testbed
Transmitted signal in the time domainx(t)Fourier transform
Transmitted signal in the frequency domain

X(f)Complex transfer functionof Tx - BPF G1(f)1. I NTRODUCTION

Complex transfer functionof stored channel G2(f)Ultrawideband (UWB) signals use a large bandwidth with low
power emissions. This can lead to new interesting possibilities for both communications and radar applications [1]. For
this reason, UWB systems have been investigated for many
years. UWB radio is adopted as a physical layer for highand low data rate wireless personal area networks (WPANs).
Standardization activities of WPAN have been performed in
the IEEE 802.15 Task Group 3a and 4a for high and low data
rate systems, respectively. In order to compare standardization
proposals, they have developed new standardized channel
models [2], [3]. These channel models are stochastic models
based on probability theories and thus are not suitable for
evaluating the transmission in realistic environments. On theother hand, stored channel simulations are often used to reflect
the properties of propagation channels in a real environment
for device testing [4]. However, whether the stored channel can
exactly simulate the real propagation environment has not beenverified. Besides, there have not been many stored channel
models [5], [6], [7]. To solve the problems stated above, ourresearch objectives are set to the following three-fold: (1) toverify the effectiveness of the stored channel simulation in
evaluating UWB system performance, (2) to develop a stored
channel model for an indoor office environment and (3) toevaluate the transmission performance using the developed
stored channel model.

Complex transfer functionof Rx - BPF and LNA G3(f)Received signal in the frequency domain S(f)Inverse Fourier transform
Received signal in the time domain s(t)Modeled noise n(t)Noise added received signal in the time domain y(t)Fig. 2: Flow chart of received waveform simulations using a stored channel

a UWB testbed [8] to the received waveform obtained fromstored channel simulations.
We conducted two experiments: transmission experimentand propagation experiment. Both experiments were performed in the same propagation environment with the sameantenna position in a rich multipath indoor office environment.
There were numerous clutters such as metal partitions, desks,
electronic equipments around the transmit and receive antennas.
1) Transmission experiment: The transmission experiment
was performed using the UWB testbed, whose diagram is
described in Fig. 1. On the transmitter side, the signal generator equipped with 10 GS/s, 4-bit DA Converter (DAC) wasused to directly generate the UWB signals up to 5 GHz.
The transmitted signal was the impulse response of a RootRaised Cosine (RRC) filter modulated by a sinusoidal wave
with frequency fc = 4.05 GHz [8]. These signals are thensent to the transmit antenna after limiting the bandwidth to
3.1 − 5.0 GHz using the BPF.

2. VALIDATING THE EFFECTIVENESS OF UWBTRANSMISSION SIMULATION USING A STORED CHANNEL
A. Experiment
In order to verify the effectiveness of stored channel simulations, we compared the actual received waveform measured by65

Technical Report of International Development Engineering

15

1

10
Voltage [mV]TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

23

5

4

0
−5
−10
−15

0

50
100
Relative time [ns](a) actual measurement

151Voltage [mV]

10

23

5

4
Fig. 4: Floor plan of the measurement environment

0−5

TABLE 1: S PECIFICATIONS OF THE MEASUREMENT
Bandwidth
Number of sweeping points
in the frequency domain
Transmit antenna
Height of transmite antenna
Receive antenna
Receive antenna position
Room size
Measured area
Number of array elements
Inter-element space
Inter-array interval
Number of spatial samples

−10
−15

0

50
100
Relative time [ns](b) stored channel simulations
Fig. 3: Received waveformOn the receiver side, the signal received from the receive
antenna passes through the BPF to reduce the out-of-band
interference. They are then amplified by LNA and finallytime domain waveform was detected at 12 bits, 20 GS/s highspeed DSO. The averaging function of DSO was used in orderto improve the signal-to-noise ratio (SNR) in the received
waveform.
2) Propagation experiment: In order to obtain the received
waveform via stored channel simulation to compare with thereceived waveform in the UWB testbed, we used a vector
network analyzer (VNA) to measure the transfer function
G2 (f ) of the propagation channels including both transmit
and receive antennas. The transfer functions G1 (f ) of BPF
on the transmitter side and G3 (f ) of BPF and LNA on thereceiver side of the UWB testbed were also measured with
the VNA. The measured transfer functions were then used toreconstruct the received waveform as described in Fig. 2. In thereconstruction of the received waveform, the same transmitted
waveform as the UWB testbed was considered.

3.1 − 10.6 GHz
1501 points (0 ≤ y ≤ 4600)
3001 points (5000 ≤ y ≤ 7600)
UWB monopole antenna
1300 mmUWB monopole antenna(−450, −20, 1300)
7015 × 9080 × 2700 mm3
plane of 5100 × 7600 mm2
5×5
25 mm
500 mm4205 points

B. Results and discussion
The waveform obtained from the transmission experimentand reconstructed waveform via stored channel simulations
are shown in Fig. 3. Comparison of the two waveforms
revealed that multipath signals appeared at the same delayed
time compared to the direct signal in both waveforms. In
addition, focusing on the four strongest waves, it was foundthat wave 1 was the direct wave and waves 2, 3, 4 were
the reflected waves with the delayed time of 4.3, 11.4, 21.7
[ns] respectively. Moreover, the peak amplitude values in
two received waveforms were almost the same. Finally, thecorrelation of this two waveforms was 0.87, which provedthat the received waveform obtained from actual measurement
was nearly identical to the one obtained from stored channel
simulations.

66

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

Cumulative distribution function

55

Pathloss [dB]

55
50

50
45

45

40
6000
4000
2000
y−coordinate [mm]
0

0

4000
2000
40
x−coordinate [mm]Fig. 5: Spatial distribution of pathloss

1
0.8

n=2

0.6

n=8
n = 20

0.4
0.2
0
0

2
4
6
8
Power improvement [dB]10

Fig. 7: Cumulative distribution of the improvement when SRake receiver isappliedDelay spread [ns]4040system can achieve. Figure 5 illustrates the spatial distrubutionof pathloss.
2) Delay spread: The root-mean-square delay spread,which is defined as the normalized second-order central moment of power delay profile (PDP), has a significant influencein the error probability [4]. The area distribution of the delay
spread is shown in Fig. 6.
It was found that both parameters were generally dependent
on the distance between TX and RX antennas. The pathloss
was about −40 dB when TX-RX distance was small. On theother hand, it became more when TX moved far from RX.
There was less spread when TX located near RX and more
spread when TX was far from RX.

30

30
20

20

10
0

10

4000
2000
x−coordinate [mm]

0

0

6000
4000
2000
0
y−coordinate [mm]Fig. 6: Spatial distribution of the delay spread

4. P ERFORMANCE EVALUATION OF A R AKE RECEIVER
BASED ON THE STORED CHANNEL MODELS

3. D EVELOPMENT OF A STORED CHANNEL MODEL FOR
THE DESIGN OF UWB COMMUNICATION SYSTEMS IN AN

In this chapter, a practical example in using the stored channel
models for link-level simulations is implemented. Particularly,
performance of the Rake receiver was investigated using thedeveloped stored channel models.

OFFICE ENVIRONMENT

Having verified its effectiveness in the previous section, we
next develop stored channel models for transmission simulations. Basic properties of propagation channels, such as
pathloss and delay spread, are also shown to facilitate discussion on the relationship between transmission performanceand propagation characteristics.

A. Rake receiver
The Rake receiver is designed to combat the effects ofmultipath fading. The receiver utilizes correlation between
the received signal and template waveform, and finds peaks
to collect the energy from different Multi Path Components(MPCs) so that receiving SNR improvement is attained. An
all Rake (ARake) receiver utilizes all the output from thecorrelators so that it can perfectly extract received energy.
However, because of power consumption and design complexity of the hardware, the implementation of ARake receivers are
impractical [5]. Thus, a selective Rake (SRake), a compromise
between improvement of SNR and hardware complexity byusing only the subset of taps with strong energy, is often usedinstead. It can achieve reasonable improvement of SNR withless hardware complexity compared to the ARake receiver. We
will consider SRake for the subsequent analyses.

A. Measurement campaigns
The spatial transfer functions were measured in a testroom
equipped with three-dimensional scanner and a VNA as shownin Fig. 4. There are desks, chairs, televisions and sofas in theroom. The walls were made of metal. In Fig. 4, the coordinate
system is also shown. The specifications of the measurement
are shown in Table 1.
B. Extraction of propagation parameters
1) Pathloss: The attenuation due to propagation effects
between transmitter and receiver is defined as pathloss. It
determines the average signal-to-noise ratio (SNR) that a67

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

Last, the performance of SRake receiver has been evaluatedbased on the stored channel model. The improvement ofreceiving signal-to-noise ratio was discussed. The numberof Rake fingers was optimized to be 8. It was also foundthat similarly to the propagation parameters, the distance
between TX and RX antennas is a factor that affects thesecharacteristics the most.

Power improvement [dB]

6
6

5

5
4

4

3
2

4000
2000

x−coordinate [mm]R EFERENCES

3

0

0

[1] A. F. Molisch, “Ultrawideband propagation channels-theory, measurement, and modeling,” IEEE Trans. on Vehic. Technol., Vol. 54, No. 5,
pp. 1528-1545, Sep. 2005.
[2] J. Foerster et al., “Channel modeling sub-committee report (final),”IEEE P802.15-02/490r1-SG3a, 2003.
[3] A. F. Molisch et al., “IEEE 802.15.4a channel model – final report,”
IEEE P802.15-04/662r0-SG4a, 2004.
[4] A. F. Molisch, Wireless Communications, 1-st edition, John Wiley &Sons, Chichester, 2005.
[5] M. Z. Win and R. A. Scholtz, “Characterization of ultra-wide bandwidth
wireless indoor channels: a communication-theoretic view,” IEEE J.
Selected Areas Commun., Vol. 20, No. 9, pp. 1613-1627, Dec. 2002.
[6] http://ultra.usc.edu/uwb database/[7] Y. Rikuta, S. Fujita, F. Ohkubo, H. Hosoya, K. Hamaguchi, J. Takada,
T. Kobayashi, “Indoor Channel Measurement of 26 GHz Band UWBCommunication System,” in Proc. IEEE International Conference on
Ultra-Wideband 2006 (ICUWB 2006), T2B-4, Waltham, USA, Sep.
2006.
[8] K. Takizawa, I. Nishiyama and Y. Rikuta, “Experimental evaluation ofvarious UWB signaling by using a UWB test bed,” in Proc. TechnicalCommittee on Wideband System, WBS2004-8, pp. 37-42, Shizuoka,
Japan, 2004 (In Japanese).

6000
4000
2000
2
y−coordinate [mm]Fig. 8: Spatial distribution of the minimum power improvement when a8-finger SRake receiver is applied

B. Results and discussion
It is important to choose an appropriate number of Rake fingers
to meet both performance and complexity requirements. First,
in Fig. 7, we showed the cumulative distribution of the
improved power in an array when the number of Rake fingers
varied from 2 to 20. It can be seen from the figure that the
improvement of the received power significantly changes until
8-finger SRake receiver were used, but the improvement was
gradually becoming less when the number of Rake fingers
exceeded 8. The improvement reaches 4 dB on average when
8-fingers were used. Thus, we can set the optimum number ofRake fingers to 8, and subsequent analyses were done with 8
fingers.
The improvement of the received power when using an 8finger Rake receiver is shown in Fig. 8. The farther the distance
between TX and RX antennas was, the higher the power could
be improved. This was due to the power composition between
the strongest MPC and scattered MPCs. When the RX antenna
located near the TX antenna, the strongest wave dominated thereceived power. On the other hand, power contribution fromthe strongest wave was less remarkable when the RX antenna
was far from the TX antenna. Therefore, the power did not
improve much when the strongest wave was much stronger
than the other waves, i.e., when the distance between antennas
was small.
5. C ONCLUSION
This paper first verified the effectiveness of stored channel
simulations. Thus, the stored channel simulations can be used
to replace the use of actual devices, which has a great contribution in the reduction of cost in developing new products.
Second, the store channel models have been developedin an indoor office environment. The important propagation
parameters have been extracted from this measured data andtheir distributions in the area have been investigated. We
noticed that all of these parameters are strongly dependent
on the distance between TX and RX antennas. The result canalso help to predict the outage performance inside this room.

68Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

Numerical simulation on Temperature Field for Newtonian
Fluid Enclosed in Concentric Sphere
05M18171

Bian Hu

Supervisor: Yoshihiro MOCHIMARU

同心二重球殻内の同心二重球殻内のニュートン流体温度場
ニュートン流体温度場シミュレーション
流体温度場シミュレーション
辺 虎
ニュートン流体でエネルギー式に対象を絞り、三次元空間場の数値解析を
スペクトル差分法によって行う。同心二重球殻内の流体内の熱対流を伴う温度場を研究対象とし、レイノルズ数,プラントル数等のパラメータの影響及び強制速度を与えたマクロ的対流の影響による温度場の差異を把握した。

INTRODUCTION

difference scheme. Effects of a Rayleighnumber, a Prandtl number, and a convection

In engineering fields a variety of heatflow on a temperature field are discussed.

convection phenomena are treated,e.g., in

In analysis, spherical harmonic functions

the field of air conditioning design in an

are introduced as a base for spectral

airplane, and of design of electronic device

decomposition to discretize the equation in acooling, or of design of light equipment[1].

finite difference scheme.

As an example, if the fluid enclosed in

ANALYSIS

concentric sphere possesses a non-uniform
temperature field, natural convection will be

Basic Equations
Equationsproduced.

It is assumed that the fluid is Newtonianand possessing constant thermal conductivity, constant viscosity, and a constant coeffici-

R1θ

ent of thermal expansion.No slip conditions

R0

ϕon the walls apply. Dissipation terms are
neglected.

Thus, the system of equations[2][3](continuity, momentum, energy) are:Fig.1 Model
Model

In the current study, attention is focused on
the energy equation itself for Newtonian
fluid

in

a

three-dimensional

sphere,analyzed with the aid of a spectral finite

169∇・u = 0

(1)

Du
1 2
Gr
= − ∇P +∇ u + er 2 Τ
Dτ
Re
Re

(2)

DΤ
1
=
∇ 2Τ
Dτ
PrRe

(3)

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

R − , θ − and ϕ − dir-

A Spherical Harmonicarmonic Function

components in the

A spherical harmonic function used is aections respectively.

function of two arguments, θconnected

toLegendreand ϕ ,

polynomials

BOUNDARY
BOUNDARY CONDITIONSCONDITIONS

mathematically. That is[4],Yn,m (θ , ϕ ) =2n +1 (n − m)! m
Pn (cosθ )eimϕ
4π (n + m)!

The following thermal boundary Dirichlet(4)

condition apply:

Orthogonality givers

2π

π

=0

=0

∫ϕ ∫θ sin θ P (cos θ ) P (cos θ ) dθ dϕ
m
n

0
=
4π (n + m )!

m
N

(9)

TH ( R = R1 )(10)

= sin θ × (sin θ + 0.8 cos ϕ + 1)

(5)[n ≠ N ](n − m )! (2n + 1) [n = N ]T0 ( R = R 0 ) = 0

"T0"

T
π
π
∫ dϕ ∫ Y (θ , ϕ )Y (θ , ϕ ) sin θ dθ = 1
2

*

0

0

n,m

n, m

"TH"

20.025

(6)

0.02
0.015
0.01
0.005

where the superscript * stands for a conju-

0gate. spectral decomposition is as follows
0∞n

T (R, θ , ϕ ) = ∑ ∑ Tn , m (R ) Yn , m (θ , ϕ )

0.5

1

1.5

2

2.5

3

3.50

1

2

3

4

5

6

7

ϕ

θ

(7)

n=0 m= −nFig.2.
Fig.2. Temperature distribution on the surface

Decomposition of the energy equation

For forced flow velocity

Substituting Eq.(7) into Eq.(3) with the aid

U = 0；

V = 0；(11)of Eqs.(5) and (6) gives

W = sin θ × (0.05 + 0.3 cos ϕ )∂Tn,m(R)
∂τ

−

(12)

1  1 ∂  2 ∂Tn,m(R)  n(n +1)
− 2 Tn,m(R)
 2 R∂R  R
PrRe R ∂R  (8) ∂Τ V ∂Τ W ∂Τ * Y n,m(θ,ϕ)sinθ dθ dϕ}= − ∫∫U ++ ∂R R ∂θ Rsinθ ∂ϕ 

RESULTS OF NUMERICAL ANALYSIS
In this study, effects of the forced flow on

where

U , V and, W stand for velocity

thermal fields are clarified.

2
70

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

τ = "t=0"
1
3
τ = "t=3"7
τ = "t=7"15
τ ="t=15"25
τ ="t=25"
τ "t=400"
= 400

T
1.8
1.6
1.4

Fig 7

"W=xsin(x)sin(x)"W = θ sin 2 (θ )"W=sin(x)(0.05+0.3cos(y))"
W = sin (θ )(0.05 + 0.3 cos(ϕ ))

(

T

1.6
1.4
1.2

1.2
1

1

0.8

0.8

0.6

0.6

0.4

0.4

0.2

0.2

0

0

0 0.5

1 1.5

2 2.5

1

3 3.5 0

2

5

4

3

6

7

3.5

3

2.5

2

1.5

1

ϕ

θ

)

sin (θ ) sin (θ ) + ϕ 2
W ="W=sin(x)(sin(x)+yy)"

1.8

0.5

0 10.45
10.75
10.7
10.65
10.6
10.55
10.5

θ

Fig3
Fig3. Temperature
Temperature development with time

Fig6
Fig6.

Temperature difference with flow
velocity at R = (R1 + R0 ) 2

Fig 5

T

"Pr=0.72 Re=1""Pr=7.0 Re=1"

T"Pr=0.72 Re=100"
1.8
1.6

T 1.4
1.2
1
0.8
0.6
0.4
0.2
0

0

0.5

1

1.5

2

2.5

3

3.5 10.45
10.75
10.7
10.65
10.6
10.55
10.5

W = θ sin 2 (θ )W = sin (θ )(0.05 + 0.3 cos(ϕ ))θ(W = sin (θ ) sin (θ ) + ϕ 2

)

θ
Fig.4
Fig.4 Temperature difference with

Fig.7
.7 Expansion of temperature field
Fig

Parameters at R = (R1 + R0 ) 2

cosθ

"q-t" using 1:2:3:4

T
1

200
0.5

0
-200

0

-400
-600

-0.5

-800
-1000

-1

-1200
-1400

-1
-0.8
-0.6
-0.4
-0.2

"[Pr=0.72(Re=1)]"

0
0.2
0.4
0.6
0.8

sin θ sin ϕ"Pr=7. Re=1]"

"

1

1

0.5

0

-0.5

-1

sin θ cos ϕ"[Pr=0.72 Re=100)]

θ
Fig.5
Fig.5 Expansion of temperature field

Fig.8
Fig
.8 Heat flux distribution on the outer surface

371Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)"q [Pr=0.72 Re=1]"

qH

qH"q [Pr=7 Re=1]""q [Pr=0.72 Re=100]""q [Pr=0.72 Re=1]"

0.5
0
-0.5
-1
-1.5
-2
-2.5

Fig 10

-3

0

0.5

1

1.5

2

2.5

3

3.50

1

2

3

4

5

6

7

ϕ

θ

θ
Fig.12 Expansion of temperature field

Fig.9
Fig.9 Heat flux difference with Pr

qH

CONCLUSIONSCONCLUSIONS"q [Pr=0.72 Re=1]"

It is found that a spectral finite difference
scheme can be applied at least to one case ofa three-dimensional thermal field with aprescribed convective flow for a variety of
combination of parameters.

REFERENCES

θ
Fig.10 Expansion of Heat flux

[1] Chung-Hyo JUNG, Takahiko TANAHASHI, “FEM Analysis of the Thermal Fluid"q [Pr=0.72 Re=1 k=100]"qH"q [Pr=0.72 Re=100 k=100]"Flows under Coriolis Force and Lorentz

50

Force” , 2003, Transactions of JSCES, pp.1-2

0-50[2] Tetsuya Kawamura, ”Application of sim-

-100
-150

ulation of a flow”, SANKAIDO, 2005.

-200
-250

Fig 11

-300

0

0.5

1

1.5

θ

2

2.5

3

3.50

1

2

3

4

5

6

[3]Frank M.White, ”VISCOUS FLUID FLOW”, 1991, McGRAW-HILL,Inc

7

ϕ[4] Kiyosato Okamoto, ”The prospects ofFourier analysis”,

Fig.11Fig.11 Heat flux difference with Re

4
72

Asakura Bookstore,1997

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

A Study on Corrosion of Paint-Coated Steel with Defects
in Marine Environment
Student Number: 05M18188 Name: Aung Kyaw Min Supervisor: Nobuaki OTSUKI
海洋環境下における塗装鋼材の欠陥の腐食メカニズムに関する研究
アウン・チョー・ミン
現在、港湾構造物においては塗装された鋼材が建設材料として使用される場合がある。このような構造物の耐久性を議論する場合、塩化物イオンによる鋼材の腐食が重要となる。特に、このような塗装鋼材の腐食は施工時あるいは供用時に生じる欠陥部を起点として進行する場合が多い。しかしながら、その速度やメカニズムは不明確である。そこで、本研究では、①鋼材表面におけるマクロセル腐食及びミクロセル腐食の評価手法の提案、②塗装した鋼材の欠陥部の腐食メカニズムの検討、③欠陥部の腐食に与える温度や暴露状況の影響に関する検討、④実構造物における腐食状況の検証を行った。
1. Introduction
The steel structures such as in bridges, expressways and port
structures are usually exposed to aggressive environments likemarine environment that promotes chloride attack. Therefore,
without provision of paint coating as protection, the design lifeof the steel structures will not reach. There are two main
corrosion protection mechanisms of paint coating such as (1)preventing the access of corrosion generative substances (Cl-,
H2O, O2) to the steel and (2) providing the high resistant media
between cathode and anode area in the formation of corrosion
cell on the steel surface. However, still severe corrosion in paintcoated steel structures is common in our daily life because
corrosion generative substances find ways to reach the steel
through the defects on the paint coating or through thedeteriorated paint coating. Once the corrosion generating
substances reaches the steel and the onset of corrosion process
occur. After the corrosion at the paint defects is generated, theproceeding process of corrosion is controlled or accelerated bysome influencing factors namely exposure situations andtemperature. Nowadays, the corrosion in those paint coatedsteel structures becomes a serious problem from the durabilityof infrastructures point of views. Unfortunately, there are only afew studies related to above matter.
From those backgrounds, the objectives of thisstudy became as follow: (1) to propose the divided steel plateevaluation method for detailed investigation of macrocell &microcell corrosion of the steel plate (2) to investigate corrosion
mechanism of paint-coated steel with defects exposed tomarine environment (3) to investigate the influences of
exposure situations and temperature on the corrosion and (4) topropose the estimation method for macrocell corrosion rate inexisting structures.

1cm

0.5cm

Solder

Epoxy

3cm

1c
m

2. Outline of experiment
In this section, the outline of the experiments for the study ispresented.
2.1 Materials used and Manufacturing of Specimens
In this study, mainly two types of specimens were used, divided
and undivided steel plates. The divided steel plate was made byconnecting 9 steel elements each having a dimension of 1cm x
1cm x 0.5cm with epoxy in the spacing of 0.1cm to form a 3cm
x 3cm divided steel plate as shown in Fig.1. Four lead wires
Steel element
3cm

Lead wires

0.1cm

Divided Steel
Plate
Fig.1 Outline of Divided Steel Plate
were connected separately to each steel element by soldering toallow the measurement of corrosion current between each steel
element [1]. The surface of divided steel plate was thenpolished to obtain a smooth and level surface. Moreover,undivided steel plates having the same size and make up of thesame material was prepared to do a comparative study between
divided and undivided steel plates. All the specimens werebrush painted with the paint (JIS-5621standards oil alkyd resinbase paint) for general corrosion except the specimens used
for proposing the divided steel plate. The paint condition was
set into two types good and poor conditions mainly by meansof thickness of paint. The paint coating thickness are (138~188μm) for good paint and (45~65μm) for poor paint.

73

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

2.2 Experimental Set Up(1) Proposal of Divided Steel Plate
In this experiment, the bare divided and undivided steel plates
were submerged in 3%NaCl solution under the same condition.
The specimens were submerged in 3%NaCl solution tosimulate the marine environment. The corrosion polarization
behavior of those two steel plates was compared every week
until 1-month.
(2) Investigation of the corrosion mechanism of paint-coatedsteel plate with defects
In order to investigate the corrosion mechanism of paint-coatedsteel plate with defects, the painted divided steel plates withdifferent types of defect (unpainted part) in middle steel
element were submerged in 3%NaCl solution. The details ofdefect types were as shown in Fig.2. Then, the corrosion
behavior of the specimens was investigated every month until
4-months.

corrosion behavior of the specimens was investigated every
month until 4-months
2.3 Investigation Items(1) Polarization Behavior
The factors expressing polarization behavior such as halfcell
potential and polarization resistance were measured by using
the Riken Denshi Co., Ltd made CT-7 corrosion monitor.
(2) Macrocell Corrosion Rate
The macrocell corrosion current was measured between the
middle steel element and the surrounding steel elements byusing zero resistance ammeters. Then the macrocell current
density can be calculated by using Equation (2-1). The
macrocell current density can be converted into the macrocell
corrosion rate by using conversion factor ( 100μA/cm2 = 1.16
mm/year) [1].
I mac =
where, Imac
Iin
Iout
S

Type-A

Type-E

Type-B

Type-CI in − I out
S

(2-1)

: macrocell corrosion current density (A/cm2): total current flowing in (A): total current flowing out (A): surface area of steel(2) Microcell corrosion rate
The microcell corrosion current density can be calculated bysubstituting polarization resistance in the equation (2-2). Finally,
microcell corrosion rate can be calculated by using theconversion factor (100μA/cm2 = 1.16 mm/year).

Type-D

Type-H
Type-F
Type-G(Poor Pt in Good Pt) (Poor Pt w.o Defect) (Good Pt w.o defect)

I micro =Fig.2 Painted Divided Steel Plates with defects

where, ImicK
Rp
S(3) Investigation of Influences of Exposure Situations
In this experiment, the type-B specimens were exposed todifferent exposure situations such as vacuum, submerged,
air-bubble (near), air-bubble (far), wet/dry and atmospheresituations. In vacuum (simulation of deep submerged situation)
situation, the dissolved oxygen (D.O) in 3%NaCl solution was
reduced under 1.0mg/l by using vacuum pump and thespecimens were exposed in low D.O content 3%NaCl solution.
In air-bubble near & far (simulation of direct and indirect wave
zone situation) situations, the air-bubble was supplied into thesolution near/far (from-15cm) to the specimens surface in therate of 10cm3/sec by using Suisaku Company made SSPP-7 air
pump. In wet/dry situation, the specimens were exposed in
3%NaCl solution same as the submerged situation for 1-day(wet situation) and removed the solution and put inenvironmental control chamber for 1-day (dry situation). For
the atmosphere situation, the specimens were sprayed with
3%NaCl everyday and put in environmental control chamber.
The temperature was controlled at 20℃ to reduce the effect oftemperature on the corrosion. Then, the corrosion behavior of
the specimens was investigated every month until 4-months.
(4) Investigation of Influences of Temperature
In this experiment, The solution temperature was controlled at(20, 30, 40℃) for the cases such as submerged, air-bubble near
& far and wet condition of wet/dry situation by using automaticheater. In the cases of atmospheric and dry condition of wet/dry
situation, the specimens were put in the environmental control
chamber and set the temperature at 20,30 and 40℃. Then the

K
Rp S

(2-2)

: microcell corrosion current density (A/cm2): constant (=0.0209V) [2]: polarization resistance (Ω): surface area of component (cm2)(3) Calculation of steel mass loss by corrosion
First, the corrosion rate measured at each month wastransformed to (g/cm2/year). Then, the steel mass loss at eachmonth by macrocell & microcell corrosion was calculated asflows:
Corrosion Rate (g/cm2/year) x Time (1/12yr) = Mass Loss (g/cm2)
Finally, the steel mass loss until 4-month was calculated bysummation of mass loss at each month.
3. Results and Discussion(1) Proposal of Divided Steel Plate
The average halfcell potential and polarization resistance at
different measurement points on the surface of divided andundivided steel plates were determined and compared until
4weeks as shown in Fig.5. From the results, it was confirmedthat the halfcell potential and polarization resistance of thedivided steel plate was almost the same as that of undivided
steel plate and indicated that the possibility and formation of
corrosion reaction on the divided steel plate was same as that ofundivided steel plate [3]. Therefore, it can be concluded that thedivided steel plate can be used to investigate the corrosion
behavior of undivided or normal steel plate and, the proposal of
divided steel plate evaluation method for macrocell corrosionof the steel plate was successfully proposed.

74

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

macrocell corrosion in 0-month &1-month. The same trend
was observed in all other specimens. Therefore, it can beconcluded that macrocell and microcell corrosion had a relationand the defect was the major cause of both corrosions of thepaint-coated steel plate.

Reference Electrode- Ag/AgCl

-700
-650
-600
-550
0week

1week

2week

3week

Microcell Corrosion
Current density (A/cm2)

Undivided
Divided
4week

Exposure Time
Undivided
Divided

4.50E+02

Ref Electrode- Ag/AgCl

3.00E+02
1.50E+02

Microcell Corrosion
Current density (A/cm2)

Polarization Resistance (Ω)

Halfcell Potential (mV)

Technical Report of International Development Engineering

0.00E+00
0week

1week

2week

3week

4week

Exposure Time

Defect Type-A

2.00E-05

DefectG-good paintP-poor paint

G P

G P

1.00E-05

Anodic3.00E-05

0.00E+00-1.00E-05

G P

G P

G P

-2.00E-05

Cathodic

Macrocell Corrosion
Current density (A/cm2)

Fig.5 Halfcell Potential & Polarization Resistance of Divided &Undivided Steel Plate(2) Investigation of the corrosion behavior of paint-coated steel
with defects
The macrocell corrosion current density between the defect andpaint part of the specimen type A & B are as shown in Fig.6.
From the results, it was confirmed that macrocell corrosion
occurred at the defect part and, the defect acted as anode andthe paint part acted as cathode. Then the paint part became
anode as the macrocell corrosion extended to the surrounding.
The other specimens had the same trend. In type-F specimen,
macrocell corrosion was also detected between good and poor
paint parts while poor paint part acted as anode and good paintas cathode. Almost no macrocell corrosion occurred in the caseof paint coating without defect (Type-G & H). Therefore, it canbe said that the defect part in the paint coating was the maincause of macorcell corrosion.
Moreover, the microcell corrosion current density of
type- A &B specimens is as shown in Fig.7. From the results, it

0-month 1-month 2-month 3-month 4-month

Defect2.00E-05

G-good paintP-poor paint1.00E-05

Anodic3.00E-05

0.00E+00

G P

G P

G P

G P

G P

-2.00E-05
0-month 1-month 2-month 3-month 4-month

Exposure Time

Fig.6 Macrocell Corrosion Current Density of the Specimens
was confirmed that microcell corrosion mainly occurred at
defect until 1-month and it spread to the surroundings after thesurrounding paint part degraded due to the cathode reaction of

Defect Type-A
DefectG-good paintP-poor paint2.00E-04

G P

G P

G P

G P

1.00E-04

0.00E+00
0-month 1-month 2-month 3-month 4-month

Exposure Time

3.00E-04

Defect Type-B

Defect2.00E-04

G P

G P

G P

G-good paintP-poor paintG P G P

1.00E-04

0.00E+00
0-month 1-month 2-month 3-month 4-month

Exposure Time

8.00E-02(Good Paint)(D.O-6.8mg/l)+
Bubble Near

6.00E-02
4.00E-02
2.00E-02(D.O-6.8mg/l)(D.O<1.0mg/l)(D.O-6.8mg/l)+
Bubble Far

0.00E+00
VacuumSubmerged

Bubble(Near)

Bubble(Far)Continuous Submerged Situations

75

G P

Fig.7 Microcell Corrosion Current Density of the Specimens(3) Investigation of Exposure Situations and Temperature
The calculated steel mass loss by macorcell and microcell of
the good painted type-B specimen are as shown in Fig.8. From
the results, it was confirmed that the calculated steel mass loss
in the specimens by macrocell and microcell corrosion was
influenced by exposure situations. The key factors for influenceof exposure situation were the dissolved oxygen content in thesolution and surface wetness condition. The vacuum situation
with the low D.O content gave the lowest calculated steel massloss by macrocell and microcell corrosion. The air-bubble(near) situation which had high D.O content (6.8 mg/l) in
addition with air-bubble near to the specimens gave the highest
steel mass loss by macrocell corrosion and second highest steel
mass loss by microcell corrosion. Although the wet/dry
situation showed the highest calculated steel mass loss by
microcell corrosion, its macrocell corrosion was very low when
it was dry. Therefore it can be concluded that the air-bubble(near) situation was the most severe exposure situation for
macrocell and microcell corrosion for paint-coated steel.
In order to investigate the influence of temperature,
the calculated mass loss by macrocell and microcell corrosion
under air-bubble (near) and wet/dry situations at different
temperature are compared as shown in Fig.9. From the results,
it was confirmed that the higher the temperature the higher the
calculated steel mass loss by macrocell and microcell corrosion.
In the case of the submerged situation, the higher thetemperature the lower the D.O content led to lower macrocell
Calculated Steel Mass Loss
(Macrocell) (g/cm2)-1.00E-05

Cathodic

Macrocell Corrosion
Current density (A/cm2)

Exposure Time
Defect Type-B

3.00E-04

Wet/DryATM

Wet/Dry cycle Situations

experiment results that the defect parts had the severe macrocell
and microcell corrosion rates. It can be said that the corrosion
rate of the existing structure was successfully estimated byusing the proposed method.

(Good Paint)

6.00E-01

45-cm

ATM

Wet/Dry cycle Situations

45cm

Fig.8 Calculated Steel Mass Loss by Macrocell &Microcell
Corrosion under Different Exposure Situations
and microcell corrosion. From the results, it can be said that the
air-bubble near situation at the high temperature was the most
severe situation for the corrosion of paint-coated steel withdefects.

9cm

Calculated Steel Mass
Loss (Microcell)(g/cm2)

8.00E-01

DefectPaint

Bubble (Near) Wet/Dry

-0.2
-0.4

C
B
A

DefectDefect(a) Macrocell Corrosion
Rate

6.00E-01
4.00E-01

S1

4
Defect2
0

Defect42cm-45cm

40℃

030cm-33cm

30℃

Temperature

0.2

C

6

36cm-39cm

20℃
Good Paint

0.4

AB

-Scale18cm-21cm

0.00E+00

0.1 Times

24cm-27cm

2.00E-02

C

6cm-9cm

Wet/Dry

4.00E-02AB

0.6

12cm-15cm

Defect Parts

Bubble (Near)

6.00E-02

Structure - H- Beam
Paint - Tar Epoxy

Defects

Good Paint
DefectPaint

Age - 30 Years

Exposure Situation - Wet/Dry

Corrosion rate (mm/yr)

Calculated Steel Mass
Loss (Macrocell)(g/cm2)

8.00E-02

9-cm

Corrosion rate (mm/yr)

Continuous Submerged Situations

Location - Yokosuka, Japan
Wet/Dry

0cm-3cmBubble(Far)

36cm-39cm

Bubble(Near)

42cm-45cm

Submerged

24cm-27cm

Vacuum30cm-33cm

0.00E+00

12cm-15cm

2.00E-01(D.O-6.8mg/l)+(D.O-6.8mg/l)
Bubble Near+
Bubble Far(D.O-6.8mg/l)(D.O<1.0mg/l)

18cm-21cm

4.00E-01

0cm-3cm

8.00E-01

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

6cm-9cmCalculated Steel Mass Loss
(Microcell) (g/cm2)

Technical Report of International Development Engineering

S1(b) Microcell Corrosion
Rate

Fig.11 Estimated Corrosion Rate of the existing Structure

2.00E-01

4. Conclusions(1) The divided steel plate evaluation method for macrocell
corrosion of the steel plate was successfully proposed.
(2) It was confirmed that the macrocell corrosion occurred at
the paint defect part, and the defect acted as anode and paintpart acted as cathode. Also, the macrocell corrosion was thecause of for spreading of microcell corrosion to the surrounding
paint part. The defect was the major cause of the corrosion ofpaint-coated steel.
(3) It was confirmed that the exposure situations andtemperature greatly influenced the corrosion of paint-coatedsteel with defects. The air-bubble (near) situations under higher
temperature was the most severe condition for the corrosion.
(4) The estimation method for macrocell corrosion rate of theexisting painted structure was successfully proposed. The
proposed method was successfully applied in existing structure.

0.00E+00
20℃

30℃

40℃

Temperature

Fig.9Calculated Steel Mass Loss by Macrocell &Microcell
Corrosion under Different Temperature(3) Proposal of Estimation Method for Macrocell Corrosion
Rate of Existing Paint-Coated Structure
In order to calculate the macrocell corrosion current between
individual parts of the steel plate, the circuit diagram wasconstructed on the steel plate as shown in Fig.10. Then thecircuit analysis was done base on the Ohm's law (V=IR) and
3

1

2
V1A

R R3
in

R2

R1
V1

Rout

Rin

V1B

V2

A

V1C

R out
V1B

R out

V2B

Rout
Rin
Rin

V3

R5 R out
V1C

Rout

R4

V5

V4

R7

V1C
R4

B

V3B

R6

Rout
R out

V3C

C

V2C
R9

V6
R8

V9

V8

V7

References[1] Miyazato, S., Otuski, N., Kimura, H., Estimation Method ofMacrocell Corrosion Rate of Rebar in Existing Concrete
Structures Using Non-destructive tests, East Asia-PacificConference (EASEC 8), 2, pp.531-542, 2001
[2] Tsuru, T., Maeda, R., Haruyama, S., Applying of AC.
method monitor to local corroion, Technique of Corrosion
Prevention 28, pp.638-664, 1979
[3] Fontana, M., Corrosion Engineering, McGraw-Hill, 1978

Fig.10 Circuit Diagram on the Surface of Divided Steel Plate
Kirchhoff's current law (summation of current at any node iszero) [1] and the simultaneous equations were solved by using
Mathematica 5.2 software. The input data for circuit analysis
were halfcell potential, paint/solution resistance andpolarization resistance. Finally the corrosion rate the existing
painted structure was estimated by using the proposed method.
The results of macrocell & microcell corrosion rate are asshown in Fig.11. The estimated results agreed with the76

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

Study on Antennas Embedded in Mobile Phone Housing
Student Number: 05M18194

Name: Jeonghoon HAN

Supervisor: Jun-ichi TAKADA

筐体組み込み型携帯電話用アンテナの研究
韓政勲

本研究では携帯電話内部にアンテナのためのスペースを確保する代わりに、薄型アンテナを筐体に組み込むことを提案する。一つの薄型アンテナにおいて２周波共振させることを第一目標とし、各帯域での帯域幅を確保することを第二目標とする。付加接地と切込みによりパッチ上の表面電流を調節することでアンテナの多共振化が図れるが、一方で切込みによる共振は帯域幅が狭く、これらの両立が難しいことが分かった。

1

Introduction

Table 1: Operationg frequencies of wireless communiFor last ten years we have noticed rapid development cation systems
of wireless communications, especially mobile communications. While a market for fixed telephones is ma- Systems Operating frequency
Overall bandwidth(Uplink and downlink:MHz) (MHz) [Relative BW]
ture, one for mobile phone has still been increasing in
250 [12.2%]
developed and developing countries. The requirements W-CDMA 1920 − 1980, 2110 − 2170
940 − 956, 810 − 826
146 [16.5%]
PDC
of customer and the market of mobile phone cause thatmobile phones are small sized and multifunction. Thesedays, mobile phones are handy and have various functions including, for example, an MP3 player, a digital
mobile phone can be also relieved.
camera, a mobile analog/digital TV and so on. Thus,
Table 1 shows the operating frequencies of some of thethe circuit of mobile phone has been required to be moremobile communication systems. Various mobile comdownsized and cheaper. This research is concerned withmunication systems are serviced and the bandwidths
antennas for mobile phones to approach to the solution.
for them have become wide. Preferably, a mobile phone
Half-wave dipole antennas are widely used for mobile
has the ability to adapt to multiple services. Therefore,
phones. They have advantages of easy manufacture, low
the multi-band operation with sufficient bandwidths ofprice, and good performance. However, they occupy
small antennas has become a very important design iswide space in the mobile phone and are protruded fromsue.
the mobile phone housing. These characteristics can be
restrictions in designing the mobile phones.
Low profile antennas, such as microstrip antennas (MSA),1.1 Objectives
can solve the problems. They are so thin that they canbe implemented in the mobile phones. These antennas In this research, a new antenna has been designed for a
The size of its upper
are usually placed at an edge of mobile phone and are specific model of mobile phones.
3
housing
is
80
×
50
×
12
mm
.
The
size
of its substrate is
implemented on dielectric materials with high relative
2
78
×
48
mm
.
The
heightbetweenthesubstrate to thepermittivity between 2.2 and 11. The usage of materials
with high permittivity can increase the cost of anten- housing is 4.5 mm. Therefore, the maximum3 size of thenas. They also occupy relatively inner space of mobile antennas is restricted to 78 × 48 × 4.5 mm . The detail will be described at Chapter 2 and 3. Further, thephones.
This paper proposes a new antenna for mobile phones. mobile phone was assumed to satisfy two wireless comThe antenna is proposed to be implemented in the mo- munication services in Japan; PDC (810-956MHz) andbile phone housing. In order to embed an antenna in W-CDMA (1920-2170MHz) bands. Preferably, a MSA
the mobile phone housing, the antenna should be very should resonate at dual frequencies (around 885MHz
thin. Thus, this research concentrates on MSA because and 2045MHz) and cover the PDC band (800MHz band)
MSA is the low profile antenna. A dielectric material and the W-CDMA band (2000MHz band). This refor antenna is not used; even though a dielectric ma- search focused on resonating a small sized MSA at dual
terial is not used, there is a permittivity between the frequency bands.
substrate and the antenna due to circuits interposed
between them, although the equivalent permittivity to
2 Microstrip Antennas (MSA)
the circuits is not considered in this research. Embedding an antenna into mobile phone housing and using
The MSA was first proposed by Deschanps [1] in 1953.
the air as the dielectric, it is expected that both the
The MSA is an attractive candidate for the use in mobile
inner space of mobile phone and cost of antenna can becommunications. The basic geometry of the MSA is
reduced. The complexity of manufacturing process forshown in Figure 1. The MSA provides advantages such

77

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

as low profile, light weight, low cost in manufacturing where WT is the peak energy stored, and Pt is the avand various shaped surface. However, the disadvantage erage power radiated.
is its inherent narrow bandwidth of a few percent due
Q-factor is inversely proportional to the bandwidth.
to a high quality factor.
Thus, lower Q indicates a broader bandwidth. Q-factor
can be approximately expressed in terms of the input
microstrip(patch)
impedance of the antenna as [5]feeddielectric materialha

Q = ω0

coaxial feed
ground plain
a

b)

a)

Figure 1: Microstrip Antenna a) top view b) side view
MSA can have various shapes. A rectangular MSA
is one of the most commonly used MSA. The transmission line mode in their original forms is applicable
to rectangular patches. It is impossible to adapt thetransmission line model to other shapes than rectangle.
In this thesis, an electromagnetic simulator is used. It ischaracterized by the patch area, its shape, the height ofpatch, permittivity of a dielectric material. A resonant
frequency (fr ) of a rectangular MSA can be obtained by
the following formula[2]:
fr =

c
√
2aeff εe·¸h (εe + 0.3){(a/h) + 0.262}
aeff = a 1 + 0.824,
a (εe − 0.258){(a/h) + 0.813}
(2)

(6)

where ω0 is the frequency at X(ω0 ) = 0, Z0 (ω) is thecomplex input impedance of the antenna, R0 (ω) is thereal part of Z0 (ω), X0 (ω0 ) is the imaginary part ofZ0 (ω), and Z00 (ω0 ) is the derivative of Z0 (ω0 ) with respect to ω.
The bandwidth can be determined by a frequency range
in which return loss is below a certain level. The returnloss is defined as the ratio of the reflected power to theincident power. The return loss is usually expressed indB as shown in the Eq. (7). Typically, a return lossof less than −10dB and −15dB are often need, whichcorrespond to VSWR=2 and 1.5, respectively. In thisresearch, the bandwidths of antennas are basically evaluated in accordance with the return loss of less than
−5dB. Thus, it is sufficient to evaluate the resonant frequencies and bandwidths of the antenna by the return
loss.

(1)

where

|Z00 (ω0 )|,
2R0 (ω0 )ReturnLoss[dB] = 10log10

2.2

ReflectedPower
IncidentPower(7)

Dualband Structures

In order to accommodate two wireless communication
services in a system, dualband antennas have been developed.
µ¶− 21
One of the dualband structure is a MSA with slots. By
εr + 1 εr − 1
h+
1 + 10(3) adding slots, it is possible to change the flows of current
εe =
2
2
a
on the patch. Dualband structures are introduced in
fr is a resonant frequency, c ia the velocity of light, aeff the L-slotted planar inverted-F antenna (PIFA) [6], thePIFA with H-type slit [7], and so on.
is an equivalent length,
Another dualband structure is a MSA having electriεe is an effective permittivity, and h is a height between
cal shorting pins. The shorting pin changes the input
the ground and the patch.
Roughly, a resonant frequency of a rectangular MSA is impedance of antenna. Thus, it is usually used forinversely proportional to the size of patch. The length impedance matching of antennas. However, adding thea of the patch is usually half-wavelength. Achievable shorting pins to antennas, it can also affect resonant frequencies of antennas. Generally, the structures having
relative bandwidth is known as a few percents.
the shorting pin is used together with slots.

2.1

Performance Evaluation of MSA

Simulation
The MSA is a resonant device. Thus, its impedance 3
varies more easily than its radiation pattern. The antenna impedance is related to the bandwidth of an an- Our objective is to develop a dualband antenna for thetenna, which determines an operating frequency range specific mobile phone. The requirement for an antennais that the maximum patch size of the antennas is 78 ×of the antenna.
48 × 4.5 mm3 as shown in Figure 2. An antenna is
The bandwidth(BW) can be defined by a voltage standingdesirable to resonate around 885MHz and 2045MHz.
wave ratio (VSWR) and a Q-factor of antenna as [4]:BW =VSWR − 1√
Q VSWR

(4)

3.1

(5)

The mobile phone under consideration has a sub-LCD
window on a side of the upper housing. The space
for the sub-LCD window should be prepared on thepatch. The slot for the sub-LCD results in the antenna

The Q-factor is defined as follow[3]
Q=ω

WT
,
Pt

78

Effects of a sub-LCD window

Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

with the size of the plate 15 × 4.5mm2 (shorting plate
B). Resonant frequencies at the 800-900MHz were selw
dom obtained due to the inherent structure of the and
24mm
48mm
tenna, i.e. the sub-LCD window and the fixed height.
The poor resonance at the lower band was obtained byadding a shorting plate. The antenna with the shorting
22mm
28mm
plate A resonates 820MHz and 2689MHz. The former
resonant frequency was closest to the lower design fre78mm
quency. However, its return loss at the frequency could
not satisfy the reference value. On the other hand, theFigure 2: The schematic layout of the upper housing
antenna with the shorting plate B resonates at 739MHz
and 2689MHz. And the antenna with the shorting plate
B has better impedance matching than the antenna withimpedance mismatching at ordinary feeding positions, the shorting plate A. These advantages of each shortingwhich are generally on the center line of the patch in a plate can be combined by the combination of the shorting plates. Using the combination of the shorting plates,
rectangular MSA.
the advantages of the shorting plates A and B could beobtained. The resonant frequencies are 864MHz and
3.2 Design of Dualband Structure MSA 1833MHz. Blue line in Figure 8 indicates the return
This section introduces the dual resonant structure of loss of the antenna with the combination of shorting
antennas. It is concentrated on the amendment of patch plates A and B. The resonant frequencies are closer toshape because other factors such as the thickness and the design frequencies (885MHz and 2045MHz).
permittivity of a dielectric material are restricted to the
fixed values.
In the design of a dualband MSA, the following threesteps were proposed. The first step is to find a preferable
feeding position by moving the feed point. The second
step is to match the impedance at a desired resonant frequency by adding electrical shorting plates. The lower
design frequency of 885MHz is prior to the upper design frequency of 2045MHz because it is more difficult
to resonate a small sized antenna at a low frequency.
Finally, the current paths are controlled by slotting a Figure 4: The currents on the patch at 864MHz in thepart of the MSA. It is found that the desired dual reso- antenna presented in Fig. 3 (length of the arrow correnant frequencies can be obtained by designing the MSA sponds to magnitude)
following the proposed design rule.
In the first step, the feed was moved on the patch. The
feeding position for impedance matching is obtained at
the right upper edge with d = 63mm and w = 6mm.
However, its resonant frequency is 2398MHz. This is
still high to the desired frequencies.
12mm

originshorting plate A

shorting plate B

0

d

w

12mm
24mm

48mm

feeding
position
22mm

Figure 5: The currents on the patch at 1833MHz in
the antenna presented in Fig. 3 (length of the arrowcorresponds to magnitude)

28mm
78mm

Figure 3: The positions of the shorting plates and the
The third step is to control current paths by slotting
feeda part of the MSA such that the desired dual resonant
frequencies are obtained. In the antenna with both of
the shorting plates, the current paths at 864MHz are
In the second step, the additional shorting plates are two ways as shown in Figure 4. One is from the feed
provided to resonate the antenna at a desired frequency. to the left upper edge. Another is from the feed to the
The shorting plates were added on the edges. The in- lower side, as turning at the sub-LCD window. The
teresting results were obtained in the case of two short- current does not reach the left lower edge. The curing plates. One plate was located at d = 50 mm and rent at 1833MHz started at the left upper edge and thew = 0 mm with its size 10 × 4.5mm2 (shorting plate A), feed in Figure 5. The current finally flows to the rightand another was placed at d = 63 mm and w = 0 mm lower edge. To raise the upper resonant frequency with-79Technical Report of International Development Engineering

TRIDE-2007-01, Februay 22, 2006 (ISSN 1880-8468)

out changing the lower resonant frequency, an edge ofpatch is cut out. Red line in Figure 8 indicates the return loss of the antenna whose the right lower edge is
cut out. And their current paths are shown in Figure
6 and 7. As expected, cutting the right lower edge affected the upper resonant frequency because the current
at the upper resonant frequency flows to the right loweredge;it does hardly effect the lower resonant frequency,
because the current at the lower resonant frequency seldom flow at the right lower edge. The upper resonant
frequency was shifted by 198MHz, whereas the lower
resonant frequency was shifted by 7MHz. Thus, dual
resonant frequencies at 871MHz and 2031MHz were obtained.

4

Conclusion

This paper proposed a new antenna for mobile phones,
an antenna embedded in mobile phone housing. The
resonant frequencies at the expected bands were obtained by designing the MSA following the establisheddesign rule. Further, antennas having various slots were
simulated. However, the bandwidths could not be satisfied. The restrictions imposed by the specificationsof the mobile phone housing, such as the ones for thedimension of the patch, and the thickness and permittivity of the dielectric material, are huge obstacles to enhance the bandwidths. In order to achieve the antenna
embedded in the mobile phone housing, research on enhancement of bandwidths should be performed withoutchanging the thickness and permittivity of the dielectricmaterial.

References[1] G.A. Deschanps,“Microstrip microwave antennas,”
3rd USAR Symposium on Antennas, 1953.
Figure 6: The currents on the patch at 871MHz of theantenna presented in Fig. 3 (length of the arrow corresponds to magnitude)[2] M. Haneishi, K. Hirasawa, and Y. Suzuki, Compact
and Plain Antenns, Tokyo,Japan, IEICE, 1996.
[3] N. Inagaki, N. Kikuma, and K. Isobe,“Size Reduction of an Inverted F Type Antenna,” Technical
Report of IEICE, A.P89-27, pp47-52, 1989.
[4] J.R. James and P.S. Hall, Handbook of Microstrip
Antennas, London,United Kingdom, Peter Peregrinus Ltd., 1989.
[5] D. Yaghjian and Steven R. Best, “Impedance,
Bandwidth, and Q of Antennas,” IEEE Trans. on
Antennas and Propagation, Vol.53, No. 4, pp12981324, April 2005.

Figure 7: The currents on the patch at 2031MHz of
the antenna presented in Fig. 3 (length of the arrowcorresponds to magnitude)[6] C. W. Chiu and F. L. Lin, “Compact dual-ban
PIFA with multi-resonators,” Electronics Letters,
Vol.39, No. 12, pp538-540, June 2002.
[7] Young-hun Lee, “Dual-band FIFA design for mobile phones using H-type slit,” IEEE, pp3111-3114,
March 2004.

Figure 8: The Return losses of antennas

80