Theoretical Study on Nematic Liquid Crystalline Ordering in
Thin Confined Systems
Muniriding YASEN†(Department of Physics Engineering, Division of Material Science)
Keywords: Nematic Phase, Maier-Saupe Model, Phase Transition, Effective Field,
Anchoring Wall, Uniaxial and Biaxial Order, Crossover

1.

Introduction
An order of nematic is generally soft and reacts sensitively to an external field and also is affected to acondition of boundary walls which confine the liquid crystal.

The direction of the order is changed by anelectric field under suitable anchoring condition due to walls, and competition between both agents is applied toliquid crystal displays.

Therefore, it is important to understand how the ordering and direction are controlled

by anchoring conditions.
In this thesis, nematic crstalline ordering in thin confined systems under various type of anchoring conditions
are investigated in the framework of the Maier-Saupe model based on statistical mechanism method.

In orderto clarify the mechanism of continuous change of phase, non-uniformity due to boundary effect is described interms of effective filed which is conjugate to the order parameters, and behavior of the system is analyzed byobserving loci of the effective fields on the phase diagram of the bulk.
2.

Namatic ordering in a homeotropic cell
The nematic ordering in a thin system under the homeotropic anchoring condition is studied.

As the systemthat correspondent to the homeotropic system, the surface molecule orientation that fixed into a wall surface
normal was introduced. The first order transition disappears and continuous change of phase occurs in the systemwith thickness smaller than a certain critical thickness. The mechanism of this continuous change occurring at
the thickness just smaller than a critical thickness is disclosed; a metasteble high-temperature phase changes to ametastable low-temperature one continuously via an unstable phase between them.

On the other hand, at theextremely thin system, the effective fields do not cross the transition line nor the bound lines of the metastablestates and move in an analytic domain, and the system undergoes the usual continuous change.
† Product Development Center, Optical Business Headquarters, Nitto Denko Corporation

3. Nematic ordering under external fields
Another type of anchoring wall, that is the state in which the molecule aligns to a substance in parallel isstudied. The bulk systems under external field corresponding to this are nematic phase having negative magnetic
anisotropy.

So, nematic phase transition with both

uniaxial and biaxial order parameters is studied in thetwo kind of external fields where are conjugate to theorder parameters, respectively.

Figure 1 shows aglobal phase diagram on the fields versus temperature
space is obtained first in the mean field theory, which issimilar topologically to the phase diagram of thethree-state Potts in three dimensions.

Based on thisdiagram, the ordering phenomenon in the system withplaner anchoring walls is predicted, where the firstorder phase transition changes to the second order one
as the thickness of the system becomes small.

Fig. 1 Phase diagram on the oblique hx - hy - hz
plane. TCP is the tricritical point.

4.

Effects of Biaxial Anchoring Wall on Nematic Ordering
Nematic Ordering in the system with biaxial anchoring walls are investigated as the correspondent to a bulk

nematic system exposed two kinds of external fields. The analyses are carried out in the framework of themean field theory.

The wall condition is given by a probability density function of the long axis orientation ofmolecule on the wall and by a parameter denoting strength of interaction between a nematic molecule and amolecule on the wall.
simply.

And Landau unfolding method was used in this study to do a numerical calculation

It is shown that the biaxial order suppresses the uniaxial order and the crossover between the

homeotropically and homogeneously anchored phases is proved to occur under the condition which is strictly
similar to that at bulk in the fields.

Also, when the interaction of molecule on wall was strong, it is shown thata wetting phenomenon is caused even under biaxial walls.

Author’s publications[1] M.Yasen, M. Torikai, M.Yamashita: J. Phys. Soc. Jpn., 73, (2004) 2453-2457.
[2] M. Yasen, M. Torikai and M. Yamashita: Mol. Cryst. Liq. Cryst. 438 (2005) 77.
[3] M. Yasen, M. Torikai and M. Yamashita: 2nd Japanese-Italian Workshop on Liquid Crystals (Tsu (Japan),
2004.9.5-8), Book of Abstracts, P5.

A Study of Positive Active Materials for Lithium Ion
Rechargeable Battery
Kenji Shizuka †(Energy Conversion Chemistry, Division of Materials Science)
Keywords: Lithium-ion battery, Positive active materials, Cathode
1. The objective of this study
The objective of this study is to characterize the positive active materials for lithium ion rechargeable battery and to identify
the relationship between the properties and the electrochemical performances, and to bring these results to develop the excellent
materials.
2. Summary【Spinel Li-Mn-O materials】Li(Mn2-xLix)O4-δ have been extensively studied for cathode material in large capacity
rechargeable Li-ion batteries for EV. However, it has large capacity fading upon charge and discharge cycling. It is important toclarify the mechanism of capacity loss in order to improve the cycle life. Oxygen nonstoichiometry of LiMn2O4-δ has beenreported, however, only qualitative relationship between the cycle life and oxygen deficiency ‘δ’ has been reported so far.
Therefore, the quantitative details of relationship between the cycle life and the oxygen deficiency were studied. δwas
precisely determined by chemical analysis. The capacity retention is linearly decreased with increasingδ. The changes inlattice parameter at the end of discharged state after the cycle depended onδ. They were in good accordance with thedegradation in the low voltage region of potential curve, which was strongly affected by oxygen loss and rapidly saturated
within about the 100th cycle. These results implied the influence of Jahn-Teller ion Mn3+ introduced with oxygen deficiency.
In addition, the mechanisms of manganese spinels dissolution and capacity fade at high temperature were studied. Tests on
the stability of spinels stored at different states of charge were performed. Mn dissolution took place irreversibly from thecharged state with formation of MnF2, ramsdellite-Li0.5MnO2. In the discharged state, the Mn at the surface of the LiMn2O4 wasin equilibrium with a soluble species, which leads to an increase in the cathode resistance. The amount of the dissolved Mn
increased with increasing LiPF6 concentration and that 2 equiv. of PO2F2- were generated. Concurrent formation of carbonate
decomposition products suggests the mechanism of Mn dissolution from the charged state was composed of three steps asshown in Fig. 2.
EtO

1. Decarboxylation Catalyzed by λMnO2

EtO
O

DEC
LiPF6

λMnO2

EtOH etc.

λMnO2

[Mn]
LiPF6

2. Reduction of λMnO2 by EtOH
2EtOH + 2λMnO2

MnO2

AcOEt, H2O
LiPF6

EtOH + CO2 + C2H4

AcOEt + 2MnO + 2H2O

3. Mn Dissolution
4MnO + 3λMnO2 + 8LiPF6 + 6H2O

Mn(Ⅱ), LiPO2F2, HF
etc.

2Mn(PF6)2 + MnF2 + 4Li0.5MnO2+ 4LiPO2F2 + 2LiFHF + 10HF

Fig. 1. The relationship between oxygen deficiency δ and capacity

Fig. 2. Proposed mechanism of Mn dissolution from λMnO2.

retention rate after the 100th cycle.

【Layered Li-Ni-O materials】The effect of CO2 on layered Li1+zNi1-x-yCoxMyO2 (M = Al, Mn) materials was studied.
Li1+zNi(1-x)/2CoxMn(1-x)/2O2 (Ni/Mn = 1) singularly exhibit high storage stability. On the other hand, Li1+zNi0.80Co0.15Al0.05O2 werevery unstable. The surface of CO2-treated sample is overspread with passivation which might be Li2CO3. The relationshipbetween degree of carbonation and cell performance revealed the discharge capacity was not decreased and the capacity
retention was improved during the consumption of excess lithium ‘z’, and then deteriorated rapidly above ‘z’ inLi1+zNi0.80Co0.15Al0.05O2.
† Mitsubishi Chemical Group Science and Technology Research Center, Inc.

【Layered Li-Ni-Mn-Co-O materials】Li1+yNixCo1-2xMnxO2 materials were characterized by various methods. The cell
performance depended on the composition. The rate performance deteriorated as Co content decreased, and improved as excess
Li content increased. The change of cell performance associated with composition variations was closely related to both thecation mixing (structural disorder) and the volume resistivity (electrical conductivity). The degree of cation mixing depended
greatly on Co content, and the volume resistivity depended greatly on excess Li content.
1.0E+08

1.00

y=0.050.98

Volume resistivity / Ω･cm

Occupation parameters g

0.99

y=0.10

0.97
0.96
0.95
0.94
0.93
0.92

y=0.05y=0.101.0E+07

1.0E+06

1.0E+05

1.0E+040.91
1.0E+03

0.90
0.3

0.35

0.4

0.45

0.5

x in Li1+yNixCo1-2xMnxO 2

0.3

0.35

0.4

0.45

0.5

x in Li1+y NixCo 1-2xMnxO2

Fig. 3. Relationship between composition and occupation parameter

Fig. 4. Relationship between composition and volume resistivity forg(Li) for Li1+yNixCo1-2xMnxO2 (y = 0.05 and 0.10).

Li1+yNixCo1-2xMnxO2 (0.33 ≦ x ≦ 0.475, y = 0.05 and 0.10).

X-ray absorption fine structure on layered Li(Ni, Mn, Co)O2 materials was studied. This study revealed that Li contributes tothe rate capability through the concentration of Ni3+ in the transition metal layers, which is closely related to the volume
resistivity of the material. On the other hand, Co contributes to the Li-ion conductivity through the less cation mixing rather than
to the electronic conductivity of the pristine powder.
The electrochemical analysis of Li1+yNi5/12Co1/6Mn5/12O2 (y = 0.04, 0.13) materials was performed to elucidate how thechange of the volume resistivity of the pristine powder contributes to the cell performance. Electrochemical impedance
measurements indicated the charge transfer resistance decreased with increased excess Li content ‘y’. The equivalent circuits’
model considering the volume resistivity was proposed.
【Li-Cu-Ni-O materials】The solid solution materials between Li2CuO2 and Li2NiO2 were studied with respect of theelectrochemical performance and the crystal structure. Pure Li2CuO2 shows only 100 mAh/g of capacity, but nickel substituted
compositions delivered reversible capacity of 250 mAh/g corresponding one equivalent of lithium. The enhancement ofreversible capacity is attributed to the stability of Ni3+ state compared to Cu3+. Li2CuO2 undergoes phase change to monoclinic
during lithium deintercalation, while the nickel-rich compositions show a new phase different from either original orthorhombic
or monoclinic structure. The latter phase shows two-step voltage profile at 2 and 4 V. The middle composition in the solid
solutions, e.g. Li2Cu0.6Ni0.4O2, indicates high capacity and high resistance against the structural degradation under high ratelithium intercalation.
3. Concluding remarks
The present study has identified the relationship between the properties and the electrochemical performances about thepositive active materials, and brought these results to develop the excellent materials.
Author’s publications[1] Y. Chida, H. Wada, K. Shizuka, Relationship between cycle life of Li(Mn2-xLix)O4-δ and oxygen deficiency δ, J. Power
Sources, 81-82 (1999) 454.
[2] T. Aoshima, K. Okahara, C. Kiyohara, K. Shizuka, Mechanisms of manganese spinels dissolution and capacity fade at hightemperature, J. Power Sources, 97-98 (2001) 377.
[3] K. Shizuka, C. Kiyohara, K. Shima, Y. Takeda, Effect of CO2 on layered Li1+zNi1-x-yCoxMyO2 (M = Al, Mn) cathodematerials for lithium ion batteries, J. Power Sources, accepted for publication.
[4] K. Shizuka, T. Kobayashi, K. Okahara, K. Okamoto, S. Kanzaki, R. Kanno, Characterization of Li1+yNixCo1-2xMnxO2
positive active materials for lithium ion batteries, J. Power Sources, 146 (2005) 589[5] K. Okamoto, K. Shizuka, T. Akai, Y. Tamaki, K. Okahara, M. Nomura, X-ray absorption fine structure study on layered
LiMO2 (M=Ni, Mn, Co) cathode materials, J. Electrochem. Soc., 153 (2006) A1120.
[6] K. Shizuka, S. Takano, K. Okahara, Relationship between physical properties and electrochemical performances for
Li1+yNixCo1-2xMnxO2 positive active materials, J. Power Sources, submitted for publication.
[7] N. Imanishi, K. Shizuka, T. Matsumura, A. Hirano, Y. Takeda, Preparation and electrochemical properties of
Li2CuO2-LiNiO2 solid solutions as lithium intercalation electrode, Solid State Ionics, 177 (2006) 1341.

Estimation, Parameter Learning and Prediction of Time-Varying
Communication Channels
Gagik Mkrtchyan †Keywords: Digital communication systems, time-varying signal propagation channels, inter-carrier interference cancellation, channel estimation,
dynamic parameter learning, model-based prediction.
1. The purpose
The aim of this research is feasibility study of high data rate communication in time varying signal propagation channels. The targets ofresearch are the effective inter-carrier interference cancellation method and provision of the estimation, dynamic parameter learning and predictionof time varying channel response.
2. The related prior arts
There is the known bound for communication systems performance, but that bound is mostly unattainable because of the channel response
knowledge error and the interference in time varying channels. The time varying channel response estimation is required for many wirelesscommunication systems because of the user mobility. In some mobile wireless applications, such MC-CDMA, the precise channel estimation is notsatisfactory and each user must predict the uplink channel response in the time varying fading channel. The main problem of linear regression based
prediction methods is the requirement of the channel’s statistics knowledge, or otherwise the large number of known data transmission (pilot
symbols or pilot tones) must be done for the precise calculation of the statistical characteristics, which leads to the significant loss of thecommunication efficiency. The inter-carrier interference (ICI) can become the strong degradation factor for channel estimation and data detection,
due to the carrier frequency offset and Doppler frequency spread. While the problem of the carrier frequency offset is comparatively unsophisticatedand has multiple suitable methods for the offset correction; the ICI caused by Doppler spread has very complex representation and there is no fully
satisfying solution. The ICI occurs because when using IDFT, the spectral density function of n-th subcarrier is zero on the frequencies of all othersubcarriers, but it is not zero even very near to them.
3. Proposed methods abstract
This research provides the detail representation of the temporal and spatial characteristics of the time varying channel models. The
employed in this thesis dynamically quantized parameter based models describe the time varying channel by limited number of parameter set,
which are assumed to be constant during some quantized period of the time, but has random behavior from one time step to another transition. The
transition between sets of the parameters depends only from the previous state of the parameters by the Markov chain rule, and the statistics of thetime varying channel model is given by the transition probabilities. The conventional preamble symbol based channel estimation methods calculate
the channel response separately for the each of preambles and then use linear and nonlinear interpolation for channel response calculation during thedata symbols, transmitted between them. In this research is employed a different approach, based on the quantized limited number of the channel
parameters, so the received information from the multiple preamble symbols is directly combined for maximum likelihood (ML) channel
estimation results calculation.
The conventional pilot tones based channel estimation, widely used in the orthogonal frequency division multiplexing (OFDM) and the multicarried code division multiple access (MC-CDMA), employs the orthogonality of subcarriers for separate sub channels responses estimation andthen linearly or nonlinearly interpolates the response through the whole communication channel. This method has two drawbacks for time varying
channel: it suffers from introduced inter carrier interference (ICI) and it has a loss of the estimation precision during the interpolation. For solution ofthese problems in thesis is proposed the time domain correction of the subcarriers to reduce effect of ICI and frequency-time domain estimation forimproving the___________________________________________† Division of Systems Engineering, Department of lectrical and Electronic Engineering, Mie University

precision. In many of the wireless communication systems it is not only required the precise channel estimation but also it’s near future prediction.
For example the orthogonality of the MC-CDMA uplink transmission can be ensured only by means of the pre-equalization, and for the timevarying channels the pre-equalization coefficients must be predicted using the estimation results of the downlink. In this thesis the quantized model
based prediction approach is provided, which doesn’t use the knowledge of channel statistics. The conventional methods based on the Wiener,
Kalman or other filtering techniques need the channel response’s autocorrelation knowledge, which is not known for the most of the practicalpropagation channel conditions..
The final part of the thesis is dedicated to the time varying channel’s parameters dynamic learning. Due to the strong tendency to theunification of the mobile communication systems, the software-defined radio (SDR) with reconfigurable signal processing devices becomes the
new standard for future generation of the wireless communication systems. SDR requires the learning of the multiple unknown parameters to realize
the potential capabilities of the cognitive communication. The dynamic learning of the time varying channel parameters is valuable for SDR because
it doesn’t employ any knowledge about temporal and spatial characteristics of the propagation paths distribution. The dynamic parameter learning
can be executed by means of the evolutional algorithms in combination with the complex weighted neural networks. Evolutional algorithms are
optimal multi dimensional search techniques, and are well known as the effective methods for the parameter learning. Evolutional algorithms can beconditionally divided into the three categories: evolution strategies, genetic algorithms and evolutionary programming. While the neural networks
and genetic algorithms were widely used for the time invariant channel response learning, they were not applied to the time varying channel’s
parameters learning. The problem actually was that both pure neural networks and genetic algorithms couldn’t be directly applied to the timevarying conditions because the convergence time is too long and small fluctuation can create the chaotic behavior. In this thesis is proposed theusage of the evolutional strategy based method instead of the genetic algorithm. The proposed method controls the fit result of the nonlinearly
dependent parameters by means of the choosing the principal components from both of the parents during the “crossover” and keeps the “elitism” of
the “population set” by means of introducing the minimal “learning and crossover ages”. The complex weighted neural networks, which areemployed for learning of the propagation paths amplitudes, are simplified to the single layer perceptrons with the linear activation functions. This
allows employing of the direct error gradient based learning of linear parameters, which decreases the number of the required learning steps and ismore robust to the fluctuations due to the noise. The convergence problem is solved by the separation of the dynamic learning process for the
nonlinearly and linearly dependent parameters. The first is realized by means of the evolutional algorithm, while the second is realized by means of
the complex weighted neural network.
4. Concluding remarks
Unfortunately the degradation of data communication performance in time varying channels is inevitable, and all the channels are timevarying to a greater or less extent. But in the most of the actual cases it is possible to strongly minimize the degradation effect if precisely to know
channel response dependence in time. This thesis proposed the novel methods of ICI cancellation, dynamic learning and prediction of channel.
1) The proposed time domain signal form correction method can effectively compensate the ICI without significant increase of the signal processing
complexity and is applicable up to the very high vehicle mobility.
2) Next contribution of the thesis is in the time varying channel’s parameters learning by means of the adaptive quantization and the model based
prediction of time varying channel, which allows compensation of the MUI in the uplink transmission of TDD MC-CDMA systems.
3) The finally proposed dynamic learning method combines the complex weighted neural networks with evolutionary learning algorithm forsimultaneous learning and prediction of channel response, and it provides good accuracy for nearly ideal compensation of the degradation due to thetime varying channel.
5. Publications[1] Gagik Mkrtchyan, Kazuo Mori and Hideo Kobayashi, “Correction of OFDM Signal Form in Time Domain to Reduce ICI Due to the Doppler
Spread and Carrier Frequency Offset,” IEICE Transactions on Communications, Vol.E88-B, no.1 pp.122-133, Jan. 2005.
[2] Gagik Mkrtchyan, Katsuhiro Naito, Kazuo Mori and Hideo Kobayashi, “Doppler Spread Estimation Method for OFDM Signal Using Mean
Square of Channel Impulse Response’s Time Derivative,” IEICE Transactions on Communications, vol.E89-B, No.10, pp2961-2966, Oct. 2006.

[3] Gagik Mkrtchyan, Katsuhiro Naito, Kazuo Mori and Hideo Kobayashi, “Evolutional Algorithm Based Learning ofTime-Varying Multipath Fading Channels for Software Defined Radio,” IEICE Transactions on Communications,
Vol.E89-B, No.12 pp.3269-3273, Dec. 2006.

Study on Improvement and Application of Titanium
Dioxide - Based Photocatalysts
Naomi Nishikawa †( Advanced Inorganic Chemistry, Division of Materials Science )Keywords ： Photocatalysts, TiO2, sol-gel method, visible-light-driven

1. Introduction
Semiconductor-based photocatalysts which create electrons and holes by absorbing light have been
applied to the design of environmental photocatalytic systems and dye-sensitized solar cells (DSSC). TiO2
has been extremely focused on its variety of applications, because it shows a number of attractive
characteristics such as high photoreactivity, low cost, non-toxicity and chemical stability. Recently, many
attempts have been made to improve and apply of titanium dioxide-based photocatalysts.
In this research, it was aimed to prepare an improved photocatalyst based on the titanium dioxide byusing the sol-gel method under the environment-friendly condition, and to evaluate their photocatalytic
activity. In addition, visible-light-driven TiO2 thin films were applied to a dye-sensitized solar cell system.

2. Adsorption and Photocatalytic Characteristics of Activated Carbon Made from TiO2-Coated Woody
Waste
Composite powder materials consisting of activated carbon and TiO2 were prepared through thecarbonization of woody waste powders coated with TiO2 gel by using the sol-gel method. The effect of theconditions for preparing TiO2 sols and the carbonization process on the properties of the resulting composites
was examined. It was found that the addition of diethanolamine (DEA) to the TiO2 sol was effective forachieving a smooth coating on woody waste-derived activated carbon. The composites obtained byheat-treating the wood powders coated with DEA and polyethylene glycol (PEG)-added TiO2 gel with at 530oC and 700 oC exhibited both the adsorption and the photocatalytic decomposition of methylene blue in water.

3. Preparation of Visible-Light-Driven Photocatalysts by Sol-Gel Method and Their Photocatalytic
Activity
The Pr-doped photocatalytic thin films were prepared by the sol-gel method. The photocatalytic activityof the films under visible light irradiation was evaluated.
(1) Doping of Pr was very effective to shift the absorption edge of TiO2 thin film toward the longwavelength. This means that the absorption edge shift toward visible wavelength side in the Pr-doped TiO2
thin film, probably through the formation of an impurity level due to Pr3+ or Pr4+ in the band gap. For
degradation of methylene blue solution, the Pr-doped TiO2 thin film showed excellent photocatalytic activity
under visible light (>400nm) irradiation, especially in the range of 400-500nm, while TiO2 thin film showedno activity under the same condition (Fig.1). In other words, it was found that Pr-doped TiO2 thin films
acted as visible-light-driven photocatalyst.
( 2 )Pr-doped titanophosphate glass thin films were prepared by the sol-gel method. Under full arc
irradiation of Xe lamp, the photocatalytic activity was observed for all samples. But, under visible lightirradiation of Xe lamp, it was only seen for Pr-doped samples.

† Mie Science and Technology Promotion Center, Industrial Research Division, Mie Prefectural Government

100: TiO2: Pr0.01

80% Degraded% Degraded

80

100

60
40

: TiO2: Pr0.01

60
40
20

20
0
0

100
200
Time/min

300

0
0

100
200
Time/min300

Fig.1 Photocatalytic activity of Pr-TiO2 film under visible light ( l >400 nm ) irradiation of Xe lamp.
4. Reduction of Mite Allergen Activities by Photocatalytic Treatment with Titanium Dioxide for a
Cleaner Indoor Environment
The reduction of mite allergen activities by photocatalytic treatment with titanium dioxide (TiO2) wasinvestigated for a cleaner indoor environment.
Degussa P25 TiO2 and recombinant Derf 1 protein wereused as the photocatalyst and the mite allergen, respectively.
With increasing the amounts ofphotocatalyst, the rate for the reduction of allergen activities increased gradually.
At the irradiation timeof 60 min, the reduction efficiencies of allergen activities were more than 96.9%, when the photocatalyst
TiO2 of more than 1 mg was used.
As a consequence, the photocatalytic treatment with TiO2 was very
effective for the reduction of allergen activities.
Since aspartic acid and glycine concentrations increased
during the initial degradation, the destruction of amino acid may become a main reason for the reduction ofallergen activities by the photocatalytic treatment.
The photocatalytic treatment for the reduction ofallergen activities is simple, easy handling and low cost.
5. Application of Visible-Light-Driven Pr-Doped TiO2 Thin Films to a Dye-Sensitized Solar Cell
SystemSince it has been expected that absorption in a visible region by the Pr-doped TiO2 layer was effective
to enhance the conversion efficiency of dye-sensitized solar cell (DSSC), the photovoltalic properties of the
Pr-doped TiO2 thin film were evaluated for the DSSC system under the simulated sun light.
As a consequence, the light conversion efficiency (η) of DSSC using the Pr-doped TiO2 thin film was
1.07%, although the efficiency was 0.71% for DSSC using the pure TiO2 thin film.
6. Conclusion
In the present study, an improved TiO2 photocatalyst were prepared by using the sol-gel method
under the environment-friendly condition.

Author's Publications[ 1 ] Naomi Nishikawa ， Masaki Murayama ， Kanichi Kamiya ：“ Adsorption and Photocatalytic
Characteristics of Activated Carbon Made from TiO2-coated Woody Waste” ， Mokuzai Gakkaishi,Vol.52,
No.1, pp.50-54,2006[ 2 ] Naomi Nishikawa, Satoshi Kaneco, Hideyuki Katsumata, Tohru Suzuki, Kiyohisa Ohta, Yasuo Takeda,
Masaki Murayama, Eiji Yamazaki, Noritsugu Hashimoto, Kazuaki Masuyama and Masashi Shoyama ：“Fabrication of Visible Light Response Praseodymium-doped TiO2 Photocatalysts and Its Application to
Dye-sensitized Solar Cells” ， Photo/ Electrochemistry & Photobiology in the Environment, Energy &Fuel,pp.311-319,2006[ 3 ] Naomi Nishikawa, Satoshi Kaneco, Hideyuki Katsumata, Tohru Suzuki, Kiyohisa Ohta, Eiji Yamazaki,

Kazuaki Masuyama, Tadanori Hashimoto and Kanichi Kamiya ： ”Photocatalytic Degradation of allergen inwater with Titanium Dioxide”, Fresenius Environmental Bulletin (FEB), accepted,

Investigation of novel asymmetric reactions and synthetic studyof a gastrin receptor antagonist, AG-041R
Takashi Emura†(Organic chemistry, Chemistry for materials)
Keywords: Gastrin, Antagonist, Asymmetric synthesis, Diastereoselective, Menthol1.

The purpose
The aim of the research is to establish an efficient synthetic route for the novel anti gastrin agent AG-041R
1, which has a tetra-substituted chiral center in the structure. The synthesis has to meet the requirement for the
manufacturing cost and scalability, therefore efficient construction of the chiral center was investigated.
2.

The substantiation of the synthesis
The approach taken for the retrosynthetic analysis
SCHEME 1. Retrosynthetic Analysis(Scheme 1) revealed that a strategy based on the
Me
Me
O
O
alkylation of the oxindole derivatives would be the
HN
N
HN
N
most straightforward in order to obtain the basic
H H
H
N
HH
N
framework of our target compound. From the outset it
+ Br
N
O O
N
O
Me
O
was realized that one of the synthetic challenges
Me
OEt
OEt
would be the asymmetric alkylation of oxindole
OEt
OEt
AG-041R 1
2
derivative 2. Therefore, the focus was on identifying
an efficient method for the asymmetric alkylation
reaction.
O
Me
Asymmetric alkylation using a phase transfer
OEt
Br
+
+
catalyst is one of the choices for construction of the
OEt
O
N
OCN
H
tetra-substituted chiral center. However, this method
still requires a relatively large amount of chiral catalyst, which has a complicated structure and is difficult toobtain.
Due to availability and ease of control of the reaction process, it was decided to investigate thestoichiometric alkylation of oxindole enolates with haloacetic acid esters using a commercially available chiralalcohol such as l-menthol as an auxiliary. Reaction of this type has rarely been studied to date, even though a
substantial amount of effort in using l-menthol as a chiral auxiliary has been demonstrated. To the best of ourknowledge, only a single report has been published in which poor stereoselectivity was shown in the alkylationof oxindole derivatives with l-menthyl chloroacetate as a chiral electrophile.
3.

The results of the asymmetric reaction
The alkylation reaction of urea 2 proceeded
efficiently in aprotic polar solvents such as DMF andDMSO and the racemic AG-041 was synthesized fromthealkylation
of
urea
2
with
2-bromo-N-p-tolylacetamide, and so we first tried thealkylation reaction with l-menthyl bromoacetate in
DMF. The resulting diastereoselectivity of the reaction
was low (Table 1, entry 1). Even at lower temperaturesin DMF, only a slight degree of improvement in thestereoselectivity was observed (Table 1, entries 2 and
3). By using solvent THF, the selectivity showed aslight increase with t-BuOK as the base (Table 1, entry4). When the reaction was performed in a less polar
solvent with a lithium cation base, the†Chugai Pharmaceutical CO., LTD.

TABLE 1. Diastereoselective Alkylation of the Urea 2O

MeNHCONHp-Tol

Base
HN
N
H H BrCH2CO2l-MenthylN2

Entrya
1
2
3
4
5b
6b, c
8
9
10d
11

O
OEt

Solvent

OEt

Solvent
DMF
DMF
DMF – Toluene
THF
THF
Dioxane
THF
THF
THF
THF

NHCONHp-Tol

O
N
3

Base
t-BuOK
t-BuOK
t-BuOK
t-BuOK
LiHMDS
LiHMDS
LiHMDS
LiHMDS
LiHMDS
LiHMDS

O
O
OEt
OEt

Equiv of base
1.1
1.1
1.1
1.1
1.0
1.0
1.0
1.1
1.2
1.1

O
+
N
4

Temp. (oC)
-10
-60
-78
-10
20
20
0
0
0
-15

O
O
OEt
OEt

Ratio (9:10)
60:40e75:25e80:20e75:25e90:10e
92:8e
92:8f
93:7f
94:6f
94:6f

a
Unless otherwise noted, the reaction was carried out with 1.2 equiv of l-menthyl
bromoacetate. b1.1 equiv of l-menthyl bromoacetate was used. cLiHMDS in hexane was used.
d
1.3 equiv of l-menthyl bromoacetate was used. eDetermined by 1H-NMR analysis of thecrude reaction mixture. fDetermined by HPLC analysis of the crude reaction mixture.

diastereoselectivity was greatly increased even at rt (Table 1, entry 5). Although dioxane showed the best resultsin terms of diastereoselectivity for the reaction (Table 1, entry 6), it requires special handling because of its
possible carcinogenic activity. Therefore, THF was selected for further study. The diastereoselectivity increased
when the reaction was performed at lower temperature (Table 1, entry 8) and was improved by the use of aslight excess of base (Table 1, entries 9 and 10).
Finally, we selected the reaction conditions shown above (Table 1, entry 9) for further scale-up. We havesuccessfully performed the reaction on a scale of 100 g with a 91.8:8.2 diastereoselectivity ratio.
Recrystallization of the crude product
SCHEME 2.
from a methanol/water solution gave a
NH
NOMe
OEt
Pd/C
Br
H NOMe HCl
O
NOMe
K HPO
H
OEt
NCO
diastereomerically pure product inN
ONO
KOt-Bu
N
O
N
O
OEt
OEt
CH CN
55% yield.
80%
55% in 2 steps
H
H
DMA
100%
8
5
6
7
OEt
OEt
110 C, 4h
The transformation from theMeOmenthyl ester 3 to the final compound
NHCONHp-Tol
NHCONHp-Tol
NHCONHp-Tol
H
p-TolNH
LiHMDS
HN
N
KOH
OH
N
O
1 was performed by the hydrolysis and
EDC, HOBt
H H BrCH CO l-Menthyl
O
O
O
EtOH
H
O
THF
N
O
amidation by using EDC-HCl as a
N
O
N
O
N
O
Toluenereflux
OEt
OEt
OEt
OEt
coupling agent, the efficient total
OEt
2 OEt
3 OEt
9
1 OEt
synthesis was achieved (Scheme 2).
2
2

2

2

4

3

o

2

2

2

2

FIGURE 1. l-Menthyl Bromoacetate

4.

Stereochemistry of the reaction
H
From the results shown in Table 1 and because the use of a moreH
O
coordinatable lithium base gave the best stereoselectivity, there is possible
Br H
coordination of l-menthyl bromoacetate with the lithium enolates at its
O
carbonyl group in the transition state. In the most stable conformer of thel-menthyl bromoacetate the carbonyl group should be directed downward
as depicted in Figure 1, which suggests that the oxindole FIGURE 2. Plausible Favored and Unfavored Transition States
enolates should approach from the down side of theconformer coordinating with the lithium cation, and thusO
HN
the approach of the enolates from the upper side would
HN
N
HN
H
O
O
be eliminated. And the results with the other type of
Li
N
O
N
O O
oxindole derivatives indicate that the carbonyl group in
O
R
R
Br
major
the nitrogen protective group of the oxindole enolatesOshould participate in the coordination with the lithium
Transition state A
favoredcation. Based on these results, the two transition states
shown in Figure 2 are possibilities. Of the two, transitionBrstate A is preferable due to the steric interaction betweenOO
the isopropyl group in the l-menthyl group and the
HN
HN
N
H
HN
O
nitrogen protective group in transition state B.
O
O
N

5.

Conclusion

R
Transition state B

O

Li

N
R

O

O

minor

unfavored

A highly efficient asymmetric alkylation of oxindole derivatives has been demonstrated by using
commercially inexpensive l-menthol as a chiral auxiliary. The use of lithium cation is essential for the highdiastereoselection. A six-step chromatography-free synthesis of the optically pure oxindole derivative 1 wassuccessfully achieved by the finding of the diastereoselective reaction in overall yield of 26%.
Author’s publication
T. Emura, T. Esaki, K. Tachibana, M. Shimizu. Efficient Asymmetric Synthesis of Novel Gastrin Receptor
Antagonist AG-041R via Highly Stereoselective Alkylation of Oxindole Enolates. J. Org. Chem. 2006, 71,
8559.

Fundamental Research on Mechanical Evaluation of
Cardiac Function Using Numerical Simulation
Masakazu Tsutsumi(Material Science, Graduate School of Engineering)
Keywords: Biomechanics, Biomechanics, Human left ventricle, Myocardial wall motion,
Finite element method, Mechanical evaluation, Cardiac function
1. Introduction
It is extremely important to estimate quantitatively the mechanical functions of the left ventricle froma viewpoint of the medical diagnosis of heart disease. In order to estimate the ventricular functions, theauthors have constructed a fundamental numerical simulation system based on a finite element ventricle
model connected with both an electric stimulus transmission model and a blood circulation system model.
2. Mathematical model of left ventricle and finite element simulation
The mathematical models of the left ventricle are composed of the mechanical model of the
myocardial fiber, the electric stimulus transmission model and the blood circulation system model. The
myocardium consists of numerous contractile muscle fiber elements called “sarcomere” which produce anactive tensile force and contract by themselves after receiving the electric stimulus sent from the so-called
pacemaker of the heart. The magnitude of active tensile force produced by the sarcomeres may beapproximated through a set of simple relations proposed by Beyar and Sideman. The relationships of themagnitude of active tensile force against time and length are indicated in Fig.1a and 1b, respectively. The
installation of the fiber orientation into the finite element model is realized by transforming the standard
element in the local coordinate system to the real element in the global coordinate system. As thecirculation system model, a simplified electric circuit analogy model is employed in the present numericalsimulation. The 3-dimensional geometry of the left ventricle is here assumed to be a prorate spheroid forsimplicity. The isoparametric parallelepiped finite elements (198 elements in total) are arranged along theventricular wall. The mechanical properties of the myocardial fiber and the other various material
constants are not determined for a specific individual, but chosen on the basis of studies up to the present
as well as the experience and knowledge of medical doctors.
To confirm the reliability of proposed simulation system, the numerical results were compared with themeasurements results obtained by the corresponding MR-tagging technique. Fig.2 shows three strain
components in myocardium at end systole for normal heart. As recognized from this figure, the simulated
results may reproduce well the corresponding results obtained by MR-tagging technique.
3. Numerical results and discussion
The dynamics changes and biomechanical properties are analyzed through a complete cardiac cycle.
The reliability of the obtained numerical results is verified through the comparison with the measurement
results obtained from the corresponding MR-tagging technique.
In the research work, the analysis of the cardiac functions supposed the following three heart diseases
were performed using the proposed simulation system.
(1) HHD (Hypertensive Heart Disease)(2) DCM (Dilated Cardiomyopathy) accompanied with LBBB (Left Bundle Branch Block)(3) MI (Myocardial Infarction)
From simulated results (HHD), it is considered that the amount of contraction of the myocardium
decreases when the heart becomes the elevated blood pressure condition. It is also considered that thehypertrophy of the myocardial wall in HHD is one of acclimatization developments to compensate areduction of the stroke volume.
From simulated results (DCM accompanied with LBBB), it is predicted that deformation of the
myocardial wall during early systole is affected by the change of electric stimulus transmission pathway
due to LBBB, and that during late systole is affected by the decrease of active force due to DCM.
From simulated results (MI), it is assumed that that the blood pressure is kept to some degree when theamount of contraction of the myocardium decreases 50% compared to normal heart.
4. Conclusion
The numerical system to analyze the mechanical properties and functions of the left ventricle has been
constructed in this study by combining the mechanical model of the left ventricle with the circulatorysystem model. The proposed simulation system can reproduce the generally well-accepted properties andfunctions of the left ventricle, and also makes it possible to estimate, for example, the stress distributions
which are quite difficult to measure, while all of the basic elements of the system are very simple andfundamental. By substituting higher-grade models for the present fundamental models, the more
complicated performance of the left ventricle could be reproduced in a more precise and reliable manner.
Hence, in the future we can expect objective and quantitative diagnoses for heart diseases by improving
the present numerical system.
There are still many problems to solve before we can realize the ideal simulation system. One of them isthat there is a shortage of reliable invivo information about the mechanical properties of the bio-tissues
and organ. And even if we could obtain such information, the characteristic complexity, the timedependence and the individual differences of organisms are turned to considerable barriers for thenumerical simulation. However, it is expected that numerical simulators like those presented here will pay
a major role in overcoming these difficulties.

0.5
0.4

radialcircumferential
meridional

Strain

0.3
0.2
0.1
0
-0.1
-0.2
-0.3

Measurement(28yrs, male)

Simulation

Author’s publications[1] Tadashi INABA, Takahiro NAKANO, Masakazu TSUTSUMI, Shingo KAWASAKI, Yasutomi
KINOSADA, Masataka TOKUDA, Quantitative Evaluation of Left Ventricular Wall Motion in Patientwith Coronary Artery Bypass Grafting Using Magnetic Resonance Tagging Technique, JSMEInternational Journal, Series A, Vol.49, No.4, pp. 597-603, 2006.
[2] Masakazu TSUTSUMI, Akiyoshi SHIOMI, Tadashi INABA, Yutaka SAWAKI, Masataka TOKUDA,
Finite Element Analysis of Cardiac Function in Hypertensive Heart Disease, Proceedings of 4th World
Congress of Biomechanics, Proceedings CD, 2002.
[3] Masakazu TSUTSUMI, Kazuyuki YANAGISAWA, Tadashi INABA, Masataka TOKUDA,
Mechanical Evaluation of Cardiac Function in Heart with Disease by Numerical Simulation, Proceedingsof 5th World Congress of Biomechanics, pp.563-566, 2006.

Experimental Research on Mechanical Properties ofZr based Bulk Metallic Glass
Takamasa Yoshikawa(Mechanical Properties of Materials, Division of Materials Engineering)
Keywords: Zirconium based bulk metallic glass, Amorphous, Mechanical properties,
Thermoplastic deformation, Crystallization

1. Introduction
Bulk metallic glasses (BMGs) are alloys that have various features that are not observed in conventional
metallic materials. BMGs consist of metallic elements with amorphous structures. It has the thickness and sizeenough to be applied to a structural material, i.e. bulk size, and exhibits the specific glass transition
phenomenon different from conventional amorphous alloys. The BMGs show the excellent mechanicalproperties, high corrosion resistance and exhibit a low magnetic coercive force. In particular, Zirconium-based
BMGs have not only high mechanical properties but also high glass forming ability. Due to these remarkableproperties, it is expected that Zr-based BMGs will be applied to new structural materials. However, BMGs
must be worked and formed at high temperatures into the desired shape for mechanical applications, because it
is difficult to plastically deform BMGs at room temperature (RT).
In this study, the influence of thermoplastic deformation on the mechanical properties of BMGs at RT isinvestigated. In particular, the effects of temperature, strain, and strain rate in thermoplastic deformation on thefracture stress of BMGs at RT is focused on. A comprehension of the dependence of these parameters under hot
conditions on the material strength at RT will expedite the industrial use of BMGs.
2. Experimental method
The strength properties of Zr55Al10Cu30Ni5 and Zr60Al10Cu25Ni5 BMGs were investigated at RT (293 K). To
study the influence of heat alone on the mechanical and atomic structural properties of these BMGs at RT,
specimens were heated in a muffle furnace. Subsequently, mechanical loading tests and analysis of atomic
structure by using a XRD were performed at RT after cooling. Moreover the plastic deformation behavior ofthese BMGs at high temperature was studied by using a hot loading machine to which an electrical furnace is
attached. The tensile tests were performed at various temperature conditions and strain rates. All thethermoplastic deformation experiments were stopped before the fracture of the specimens. After thermoplastic
deformation, the strength properties were studied at RT, and compared with the as-cast materials. To investigate
the thermal properties of the Zr55Al10Cu30Ni5 and Zr60Al10Cu25Ni5 BMGs, a DSC was used.
3. Result and discussion
The fracture stresses of as-cast Zr55Al10Cu30Ni5 and Zr60Al10Cu25Ni5 BMGs at RT were more than 1500 MPa.
These materials exhibited a large elastic strain, extending up to 2.0%. On the other hand, it was difficult toplastically deform these BMGs at RT. The stress-strain curves of these BMGs are similar to those of so-called
brittle materials, particularly for the tensile tests. However, several experimental results indicated that the
BMGs have mechanical properties similar to those of ductile materials such as steel, although no plastic strain
appears in their stress-strain curves at RT.

Crystallization Temperature [K]Fracture Stress (RT) [ MPa ]2000
When the Zr55Al10Cu30Ni5 BMGs were heated at a temperature
673[K]As Castgreater than 685 K in the muffle furnace, this material was1500
crystallized. Furthermore, the fracture stress of crystallizedmaterial drastically decreased at RT after cooling. The material
1000
673[K]
undergoing heating above 685 K is embrittled and loses theductile properties that the as-cast materials possess. However, the500
673[K]
Zr55Al10Cu30Ni5 BMG after being heated below 685 K was not
693[K]
crystallized and its mechanical strength can be maintained at0levels similar to that of the as-cast samples. On the other hand, the
0
20
40
60
80
100
Thermoplastic Strain (%)
Zr60Al10Cu25Ni5 BMG also maintains the as-cast strength
properties after being heated at a temperature below 710 K.
5.0E-05 [1/s] 673[K] 2.0E-04 [1/s] 673[K]3.0E-03 [1/s] 673[K] 2.0E-04 [1/s] 693[K]Therefore, these temperatures can be considered as the upper limit
temperature in order to maintain the as-cast strength after a
Fig. 1 Influence of thermoplastic strain atheating process.
673 K corresponding to strain rate andThese materials exhibited plastic strain at a high temperature temperature on fracture stress ofalthough no plastic deformation appeared at RT. The temperature to Zr60Al10Cu25Ni5 BMG at RT
allow the plastic deformation lied over not only a higher but also a673[K]As Castlower temperature than the upper limit temperature. The property of762
plastic deformation of these BMGs at a high temperature depended
757
673[K]
strongly on not only the temperature but also the strain rate.
752
The three parameters—strain, temperature, and strain rate of the673[K]
thermoplastic deformation—can be considered to have an effect on
747
the structure of the BMG, and thus on the strength after the693[K]
742
deformation, even if the environmental temperature is lower than the
0
20
40
60
80
100
upper limit temperature. Therefore, the influence of each parameter
Thermoplastic Strain (%)
on the mechanical properties of the Zr60Al10Cu25Ni5 BMG at RT was
5.0E-05[1/s] 673[K] 2.0E-04[1/s] 673[K]studied. The experimental results showed that the strength of this
3.0E-03[1/s] 673[K] 2.0E-03[1/s] 693[K]BMG at RT after the thermoplastic deformation depends strongly on
Fig. 2 Influence of thermoplastic
the amount of thermoplastic strain, the temperature and the strain rate
deformation parameters on crystallization
during deformation. When the strain rate is comparatively low, the temperature of Zr Al Cu Ni BMG at
60 10
25
5
strength of this material at RT gradually decreases with increasing the RT
thermoplastic strain. On the other hand, a higher strain rate prevents
the strength after thermoplastic deformation from decreasing. Furthermore, the decrease of the strength of thisBMG at RT is significant at a higher temperature even if the temperature is lower than the upper limit
temperature.
From the thermal analysis by using a DSC, the crystallization temperature of the specimens deformed undervarious hot conditions changed with the thermoplastic strain. The tendency of the change coincided with thoseof the strength at RT after the deformation depending on the three parameters. The decrease of crystallization
temperature is considered to imply that the material changes into the state to be more easily crystallized
because the crystalline phase is precipitated and then a local free volume increases inside the amorphous phase.

4. Conclusion
From the experimental results of Zr-based BMGs, the influence of the thermoplastic deformation on themechanical properties at RT was clarified. It is important and useful from the industrial viewpoint that these
BMGs can be freely deformed and the excellent strength properties can be maintained after thermoplastic
deformation.
Author’s publications[1] T. Yoshikawa, M. Tokuda and T. Inaba, Key Engineering Materials, Vols. 340-341, (2007) pp.113-118[2] T. Yoshikawa, M. Tokuda and T. Inaba, Journal of The Society of Materials Science Japan, Vol. 56, No.2,(2007) pp.171-177

A Study on the Modification of the Engineering Plastics
to the Advanced High-Performance Materials
Hatsuhiko Harashina†(Materials Chemistry, Division of Materials Science)
Keywords: Engineering plastics, Poly(oxymethylene), Poly(ethylene terephthalate), Poly(ester amide),
Ethanolamine, Flame retardancy, Thermal properties, Thermal analysis, Molecular orbital calculation

1. Introduction
Engineering plastics, including Poly(oxymethylene)s (POMs) and poly(ethylene terephthalate) (PET),
have superior physical, mechanical, thermal, electrical and chemical properties, and they have been used inmany industrial fields such as the automotive, electrical/electronic, construction and consumer markets.
Moreover, in order to apply POM and PET to various industrial fields as the further advanced materials, theprogressive and drastic improvement of the fundamental issues having POM and PET needs to be performed.
As a design methodology for the advanced high-performance engineering plastics, the formulation design bypackaging of the optimized functional additives (ex. flame retardants) or the structural modification byintroducing of hydrogen-bondable amide groups to the polymer main chain can be expected.
The present paper describes the results of a study on the modification of the engineering plastics to theadvanced high-performance materials, including POM and PET, the objective of which was to produce newadvanced plastic compositions and polymers, characterize them, and obtain the fundamental design
methodology to create the advanced high-performance materials.
2. Results and Discussion
2.1. A Study on Design Methodology for High-Performance (flame retardancy) of POM [1]New flame retardant combination system for POM was investigated. It was found that POM issynergistically flame-retarded by the combination of red phosphorus with novolac as a phenolic resin andmelamine as an aminotriazine compound. Furthermore, the detailed combustion behavior of the flame
retarded POM has been studied to confirm its flame retarding mechanism by using thermal and FTIR
analyses. The results of cone calorimetry, thermogravimetry and FTIR analysis suggested that the flame
retarding mechanism is the intumescent char formation in the condensed phase. Novolac having a phenolichydroxyl group is miscible with POM, and in the flaming process, red phosphorus yields phosphine and its
acidic product such as phosphoric acid due to hydrolysis and oxidation reactions.
In addition, all of novolac, melamine and phosphine are able to readily react with formaldehyde
generating from POM during burning to give the reinforced and cross-linked char network through thepolyaddition and polycondensation reactions. Therefore, the red phosphorus/novolac/melamine ternary
combination system could synergistically promote the high flame retardancy of POM without the flaming
drips.
2.2. A Study on Design Methodology for High-performance of semiaromatic polyesters [2, 3]
Semiaromatic poly(ester amide)s (PEAs) were synthesized by the melt polycondensation of
ethanolamine (EA)-derivatives with dimethyl terephthalate (DMT) and ethylene glycol (EG) in the presenceof tetrabutyl titanate as a catalyst, and their crystallization and thermal properties were investigated.
†Polyplastics Co., Ltd.

Introduction of an amide group into the semiaromatic polyesters such as poly(ethylene terephthalate)(PET) gave the PEAs (EA-modified PET polymers) having an increase of melting point. Interestingly,
these PEAs were found to decompose at a lower temperature than PET on the basis of the TGA analysis.
Moreover, the direct pyrolysis/mass spectrometry suggested that an initial step of the thermal decomposition
was a b-CH hydrogen transfer reaction via a six-member ring transition state at the ester-ethylene-amide unit,
at which the carbon-oxygen bond scission takes place to yield carboxyl and N-vinylamide end groups(Scheme 1). Furthermore, molecular orbital calculations using trimer model molecules supported strongly
that the b-CH hydrogen transfer reaction in the thermal decomposition of PEAs occur more easily at the
methylene group next to the amide group in an ester-ethylene-amide unit rather than at the methylene groupnext to the ester group in an ester-ethylene-ester one (Figure 1).

Scheme 1.
The b-CH hydrogen-transfer reactions of (a) PET, (b) poly(ethylene terephthalamide), and (c) PEAs.

Figure 1.
Frontier electron densities on ethylene protons for trimer models of PET and PEAs.

3. Conclusionsi) The red phosphorus/novolac/melamine ternary combination system could promote the high flame
retardancy of POM without the flaming drips. The flame retarding mechanism was the intumescent charformation to give the reinforced and cross-linked char network through the polyaddition andpolycondensation reaction in the condensed phase.
ii) Semiaromatic poly(ester amide)s (PEAs) were synthesized by the melt polycondensation of ethanolamine
derivatives with DMT and EG. The primary thermal decomposition of PEAs was a b-CH hydrogen transfer
reaction via six-member ring transition state in the ester-ethylene-amide unit to form carboxyl andN-vinylamide end groups. The thermal decomposition behavior of PEAs was explained well by
semiempirical molecular orbital calculations on the basis of timer model molecules.
Author’s Publications[1] H. Harashina, Y. Tajima, T. Itoh, “Synergistic effect of red phosphorus, novolac and melamine ternary
combination on flame retardancy of poly(oxymethylene)”, Polymer Degradation and Stability, vol.91,
1996-2002 (2006).
[2] H. Harashina, T. Nakane, T. Itoh, “Synthesis of Poly(ester amide)s by the Melt Polycondensation of
Semiaromatic Polyesters with Ethanolamine and Their Characterization”, Journal of Polymer Science Part
A: Polymer Chemistry, vol.45, No.11, 2184-2193 (2007).
[3] H. Harashina, T. Nakane, T. Itoh, “Synthesis and Properties of Semi-aromatic Polyesteramides Using
Ethanolamine”, Polymer Preprints, Japan vol.53, No.2, 4119 (2004).

A Study of the Reception Space Planning of Homes for the Elderly
Hideki FUJIEDA(Architecture and planning,Division of Systems Engineering)
Keywords:Visitor,Residents,Change of receptions,State of receptions,Setting,Nursing home
1. The significance of setting reception space of a nursing home as a research theme
The act of receiving guests in a homelike environment has been considered to be an important act that ties people inside
the house and those outside, and various types of reception space have been created in the private homes. However, in spite ofits importance, the reception space of nursing homes has never been discussed in detail, nor even the actual situation is grasped
yet. Current national policy for the preparation of the nursing homes is to attain the continual homelike environment before andafter entering them. Thus, more and more attention is to be paid to the reception space. In this sense, it is necessary to grasp thecurrent conditions of the reception space of the nursing homes and to propose a guideline to them.
2. 2 Purpose of the Study
A study of the reception space planning is a new perspective. Therefore, there are almost no previous studies of this matter.
This study aims to set up some guidelines for the nursing homes to become a place where both the residents and the visitors
feel at home. This is done by (1) overviewing the current trends in order to extract their characteristics, and (2) analyzing thereception space in detail.
3. Definition of the reception at a nursing home
It is defined as residents and visitors "sharing time comfortably."
4. Status of visitors
There are two types of visitors to nursing homes: one is those who come to see the residents individually; the other is thoseadministrating the programs of the nursing homes.
5. The institutions investigatedIn order to grasp the current situation of various types of homes for the elderly, the institutions investigated are those
covered by the public nursing care insurance set by the Ministry of Health, Labor and Welfare. The author sent out aquestionnaire to a total of 55 institutions: 19 special nursing homes, 9 group homes, 10 care houses, 2 private nursing homes,
13 geriatric health care institutions, 2 geriatric medical care hospitals. In addition, the author visited 11 institutions: 4 special
nursing homes, 3 geriatric health care institutions, 1 geriatric medical care hospital, 1 group home, 2 care houses.
6. The method of this study and its results
This study is carried out in five phases through "interviews and questionnaires to people involved in the act of reception" and"the observations of the receptions space of the institutions."[Issues of the Reception at the Homes for the Elderly]
Processes and problems involving the act of reception (chapter 2) aim to grasp the actual situation of the reception at theinstitutions, to organize the significance of reception and the problems of reception at institutions, and to propose the topics tobe discussed in the following chapters. What is clarified here is that the act of reception is necessary for residents of theinstitutions since it facilitates psychologically affirmative effects and that the institutional programs are necessary especially at
the institutions for the aged since it leads to heartwarming mingling at various places in the site of the institutions. Following
matters are to be considered for individual visits: frequencies of the visits, location of the institutions, equipment of the rooms,
reception at the common space, and the awareness of the people involved in the act of reception. As to the institutional
programs, environments that generate conversations and the awareness of the people involved in the act of reception are to beconsidered.
As to the changes of reception after moving into institutions (chapter 5), some specific characteristics of the life in the
NAITO Architects Co.,Ltd,Dr.Eng.

institutions are extracted in contrast to those found in the private homes. What is clarified to be of importance as to the locationof the institutions is the accessibility by private cars. This is because many of the visitors visit the institutions by car regardlessof the distance from their residence. Next, reception at the common space and the attitudes of the inmates and the guests areinvestigated. The act of reception is found to be performed not only in the private rooms but also at various places inside andoutside the institutions. The private rooms are preferred because in there both the guests and the inmates are able to enjoy
private conversations in an intimate atmosphere. On the other hand, the guests have to care a lot about the health condition of
the inmates and the inmates are sometimes reluctant to show their private space to the guests. That is why they go outside theprivate rooms for a walk or for a change of mood. In order for the inmates to enjoy the reception outside the private rooms"without feeling any constraints," following places are preferred: half-private living rooms, dining rooms, corridors, entrance
halls, meeting rooms, bathrooms and outside institutions. Of course, different institutions have different facilities (i.e., presenceor absence of certain equipment) and these also affect the details and frequencies of the act of reception.
[Extraction of the Conditions of the Reception Space Planning]
As to the equipment of the reception space of the institutional program (chapter 3), observation and analysis are done inorder to grasp the conditions in which the act of reception through the program prepared by each institution becomes
successful by the visitors who take part in the program and to clarify the condition for the preparation of the reception space.
As a result, six types of institutional programs (i.e., viewing-type, circle-type, service-type, luncheon-type, conversation-type,and festival-type) are proved to be used differently in the care space, in the common space outside the care units in the building,and outside the building. Also, necessary conditions for each space division are clarified.
As to the case study of how the change of the reception space affects the atmosphere of the reception space (chapter 4), in
order to clarify the factors that may improve the quality of life (QOL), observation and analysis from the point of view of therelationship between the improvement of the living conditions and the reception circumstances are done through the cases ofreconstruction of the hospitals for long-stay inpatients. As a result, necessary conditions for the preparation of the space in theprivate rooms, necessity of the places for the visitors to have a rest and exchange information and their requirements are
extracted.
As to the equipment of the reception space of the institutions (chapter 6), in order to clarify the factors that bring about thereception, observation and analysis of the equipment are done by comparing the equipment of the places where reception isdone. As a result, conditions for the staff members' ease of assistance and environmental conditions of private rooms, dining
rooms, living rooms, and the space for residents and visitors to move together are clarified.
Lastly, a reception space planning checklist is proposed by summing up the results obtained by analyzing the situation andfactors of reception, the conditions necessary for the preparation of the reception space, and the program conditions necessary
for the preparation of the reception space. Furthermore, the checklist is applied to the two institutions (new KT hospital for
long-stay inmates, and Y special nursing home) of which the author involved in designing and the designing concept is "theactualization of the homelike atmosphere." The problems of the design concept are in conformity with the ones drawn fromthe reception space planning checklist. This demonstrates that the reception space planning checklist obtained as a result of thisstudy is effective in the actual design of the construction.
The reception space planning checklist proposed in this study is useful not only in designing an institution that aims tofulfill homelike atmosphere but also in administrating an institution since it gives us beneficial information on the
improvement of the residents' facial expressions and physical conditions. The author, therefore, expects it to be widely
accepted.
Although this study focuses on the reception space planning for the institutions of the elderly, its concept is applicable tothe institutions where the mentally and/or physically malfunctioning people live or receive medical care as well as to theinstitutions for the elderly who require nursing care since the act of reception is an indispensable activity for a human to live.
Author's publications(1)H.FUJIEDA,S.IMAI,S.KINOSHITA:A Study of the Spatial Conditions of the Reception in Nursing homes, Journal ofArchitecture and Planning,No616,pp39-46,June.,2007(2)H.FUJIEDA,S.IMAI,S.KINOSHITA:A Study on Scenes of Reception in Nursing home, Regional Community Facilities
Planning and Design,AIJ,Vol.23,pp251-258,Jul.,2005(3)N.OCHI,H.FUJIEDA,S.IMAI,S.KINOSHITA:Programed Reception to Neighborhood in Nursing homes, Regional
Community Facilities Planning and Design,AIJ,Vol.23,pp259-264,Jul.,2005

A quantitative study on the camera calibration error
Takashi Fujimoto*(Machine vision，Division of system engineering)
Keyword: camera parameter calibration, camera parameter, formulization, optimization, quantification

1.Purpose
The object’s position on the image taken by the camera gives the information on the direction of avisual line cast by the camera on the object. Like a triangular surveying where the number of observing
points is two, by using the information on the direction of a visual line from two or more observing points,
3-D position for the target can be presumed by fusing these. In order to perform this with high accuracy,
we require the directions of visual lines with sufficient accuracy. Here, cameras involve some importantcamera parameters such as the principal point that is the intersection of the optical axis of camera with theimage plane, the principal distance that is the distance between the center of lens and the image plane, andthe length of 1-pixel on the image plane. So, these parameters are used for transforming image
coordinates into visual lines with respect to the camera, not only measuring the image coordinates ofobject’s feature points, but also we have to calibrate the camera parameters. So, camera parameter
calibration is the basic and practical technique in the field of 3-D measurement based on visual
information. So, in this research, we focus on the 2-plane calibration method is one of the simplest
calibration methods, in the field of calibration technique, that use the reference such as the points andlines whose configuration are known, and conduct a study of following 3 perspectives.
1. Raising the precision of camera parameter calibration
2. Quantification of the calibration error to the true value of camera parameter
3. Evaluation of 3-D estimation error due to the calibration error
The outline of each viewpoint is as follows.
1. The accuracy of 2-plane calibration method calibrates the camera parameters by taking a picture of thecalibration fiducial chart moving in the depth direction is affected by following factors.
･Amount of movement in depth direction of the calibration fiducial chart･Imaging condition while photographing the calibration fiducial chart
So, we propose a new imaging condition from the perspective improving the accuracy of calibration, andoptimize the amount of movement in the depth direction of the calibration fiducial chart.
2. New imaging condition and traditional imaging condition, in addition, we suppose the imaging
condition that is possible in practical use, and formulized the accuracy of camera parameter calibration
into the amount of movement in depth direction of the calibration fiducial chart.
3. So, the essential information obtained from 2-D image is visual angle that is the angle in whichobserving the object from camera, and is 3-D information, it is important to evaluate the accuracy ofcamera parameter calibration by using visual angle. In this research, we suppose the method to evaluate
the error of visual angle accurately, and ascertain the characteristics between camera parameter
calibration error and visual angle.
1. The candidate in this research
As before, in this research, we focus on the 2-plane calibration method. And, the parameters to becalibrated are following camera parameters.
･Principal point : the intersection of the optical axis of camera with the image plane･Principal distance : the distance between the center of lens and the image plane･The length of 1-pixel on the image plane
Suzuka Fuji Xerox Co.,Ltd *And, we are premised on pinhole model as camera model. In addition, the 2-plane calibration method isthe calibration method calibrating the camera parameters by taking a picture of the calibration fiducial
chart moving in the depth direction.
2. Calculational procedure
In this research, We assume that the error is caused from the observation error of image coordinates ofthe fiducial points arisen in which we observe the fiducial points arranged on the calibration fiducial chart.
Therefore, we calculate the error from camera to object calculated based on value of camera parameter
calibration by using law of error propagation.
3. The influence of the imaging condition on the accuracy of calibration
The accuracy of 2-plane calibration method is affected by the imaging conditon in which we take aphoto of the calibration fiducial chart. So, we optimize the distance of the calibration fiducial chart, andformulize the calibation error of principal distance under the condition that the number of fiducial points
are photed at the position of the calibration fiducial chart is near and far from camera is same (imaging
condition.1). In addition, we formulize the calibration error of principal distance and the length of 1-pixelon the image plane under the condition that the additional fiducial points are photed at the position of thecalibration fiducial chart is far position(imaging condition.2). By taking a photo of more fiducial points,
we can advance the accuracy of calibration in 2-plane calibration method.
4. The influence of the setting error of the calibration fiducial chart on the accuracy of calibration
In actual use, there is a possibility that the position and the angle of the calibration fiducial chart. So,
we formulize the accuracy of calibration if the position and the angle of the calibration fiducial chart
involve errors.
5. The evaluation of the camera parameter calibration error by using visual angle
The essential information obtained from 2-D image is visual angle that is 3-D information and is theangle in which observing the object from camera. So, it is important to evaluate the accuracy of camera
parameter calibration by using visual angle. In this research, we suppose the method to evaluate the errorof visual angle accurately, and ascertain the characteristics between camera parameter calibration error
and visual angle. And finally, we formulize the error of visual angle.
6. Conclusion
In this research, we focused on the 2-plane calibration method is one of the simplest calibration
methods in the field of calibration technique. And, we examined the accuracy of calibration under theimaging condition.1 and imaging condition.2. Then, we optimized the amount of movement in depth
direction of the calibration fiducial chart, and ascertained accuracy of calibration is improved under theimaging condition.2, and formulized the calibration error of camera parameters under each imaging
condition. In addition, in practical use, there is a possibility that the position and the angle of thecalibration fiducial chart, so that we formulize the accuracy of calibration if the position and the angle ofthe calibration fiducial chart involve errors. And finally, so the essential information obtained from 2-D
image is visual angle that is 3-D information and is the angle in which observing the object from camera,
we ascertained the characteristics between camera parameter calibration error and visual angle, andformulize the error of visual angle.
Author’s publications[1] Y.NOMURA, T.FUJIMOTO, and D.ZHANG, “A Formulization of Camera-Parameter Calibration Errors andResultant Visual-Line Estimation Errors” Transactions of the JSME(Series C) Vol.71, No.706(2005) pp.147~ 155[2] T.FUJIMOTO, Y.NOMURA, “ Relationships between the camera calibration accuracy and the imaging
condition” Transactions of the JSME(Series C) Vol.72, No.717(2006)[3] Takashi Fujimoto, Yoshihiko Nomura, “A formulation of camera parameter calibration error due to reference
point setup error”, PROCEEDINGS of Asia International Symposium on Mechatronics (AISM) 2006, December
12-15, (2006), CD-ROM[4] T.FUJIMOTO, Y.NOMURA, and D.ZHANG, “Theoretical Error Analysis with Camera Parameter Calibration”
PROCEEDINGS OF SPIE in Machine Vision and its Optomechatronic Applications, Vol.5603(2004), pp.182~
190.
[5] T.FUJIMOTO, Y.NOMURA, and D.ZHANG, “A formulation of camera parameter calibration considering all
errors about reference points ” PROCEEDINGS OF JSME in Mechanical Engineering Congress ,Japan(2006),
Vol.Ⅶ, pp.327~ 328.

Study on Systematization of Porous Concretewith Small to Large Size Aggregates
Akihiro Maegawa †( Architecture, Division of Systems Engineering )Keywords ： Porous concrete, Aggregate, Particle size, Environment, Geometrical model

1. INTRODUCTION
Porous concrete is environmentally friendly material that has continuous voids inside to be used widelyin various fields. The performance of porous concrete is largely determined by the size and quantity of theinternal pores, and pore size is controlled and adjusted by the aggregate grain diameter used as well as thevolume of bonding material. However, according to the past papers, most of studies on the porous concrete ，the aggregate of the particle size 2 - 20 mm has been used. Therefore, in this study, to expand the usage ofporous concrete, manufacturing method and fundamental properties of porous concrete with small to large size
aggregates have been examined. With small size aggregate ( particle size: 0.6 - 2.5 mm ) , it was defined as
small particle size porous concrete, while with large one ( particle size: 40 - 400 mm ) , it was defined aslarge particle size porous concrete. Table 1 shows the placement of this study in porous concrete. In addition,
to analyze the theory of compressive strength of porous concrete, a geometric model with ideally sphericalaggregates was proposed.
2. EXPERIMENTAL STUDY ON FUNDAMENTAL PROPERTIES OF SMALL PARTICLE SIZEPOROUS CONCRETE
Basic properties of small particle size porous concrete ( hereafter referred to as “ SPOC ” ) were
examined. According to the test results, SPOC was superior in water retention and pumping compared tonormal particle size porous concrete ( particle size: 5 - 13 mm ) . Fig. 1 shows that the smaller the aggregate
grain, the higher the water retention and pumping rate. It should be noted, however, that permeability of
SPOC is rather low.
3. EXPERIMENT AND MODELING ON RELATIONSHIP BETWEEN COMPRESSIVE STRENGTH
AND VOID RATIO OF POROUS CONCRETE
There was little theoretical examination on the relationship between void ratio and compressive strength
Table 1 The placement of this study in porous concrete
The object of this research

Area investigated in the past

The object of this researchClassification
Small particle size porous concrete( Porous mortar )Aggregate

Particle size
Kinds of material

0.6～ 1.2mm

1.2～ 2.5mm

Porous concrete
2.5～5mm 5～ 13mm 13～20mm 20～40mm

Crushed fine stone, Molten fine slag

Crushed stone, Recycle aggregate

Mixing ｂｙ mixeｒ

Mixing ｂｙ mixeｒ

Tamping vibration / Spraying

Tamping vibration

Construction method

Large particle size porous
concrete
40～ 400mm
Crushed stone, Concrete rubble
Integration of coarse aggregateby spraying binder

Water retentive and pumping
pavement
Application field

Water penetration pavement
Base for greenery(G rasses)(Trees)(Fishing reef)† Mie Science and Technology Promotion Center, Industrial Research Division, Mie Prefectural Government

6. CONCLUSIONSIn the present study, a new performance of porous
concrete was recognized. Moreover, a geometrical model
concerning compressive strength of porous concrete wasproposed. In the future, it is expected that results of this study
contribute to the development for porous concrete.
Author's Publications

Water retention rate (％)(aggregate
1size)◇：5-13mm△：1.2-2.5mm●：0.6-1.2mm

40
30
20
10
0

10

20

30

40

Void ratio (％)( a ) Water retention rate
Water pumping rate (％)

5. APPLICABILITY AS FISHING REEF OF LARGEPARTICLE SIZE POROUS CONCRETE
Experiments have been carried out to confirm
applicability of LPOC as the fishing reef. Firstly, to confirm
stability of fishing reef placed in a sea, LPOC was sunk in thelarge-scale waterway where waves and currents can begenerated. It has been found that void of fishing reef
decreases the power of current and wave, and the effect onwave is especially large. Next, to confirm applicability as afishing reef for marine creatures, 70 LPOC blocks ( size: 1 m
× 1 m × 0.7 m ) were sunk into Ise bay of Mie Prefecture.
After six months, lobsters, abalones, sea cucumbers etc. were
found to live in LPOC. Judging from the result, LPOC isuseful to fishing reeves.

50

0

50

(aggregate
1size)◇：5-13mm△：1.2-2.5mm●：0.6-1.2mm

40
30
20
10
0

0

10

20

30

40

Void ratio (％)( b ) Water pumping rate
Fig. 1 One example of test results of small
particle size porous concrete
Derection 2

Derection 1

F

4. FUNDAMENTAL STUDY ON MANUFACTURING
METHOD OF LARGE PARTICLE SIZE POROUS
CONCRETE
Manufacturing method of large particle size porous
concrete ( hereafter, referred to as “ LPOC ” ) using concreterubble was examined. As the aggregate grain diameter waslarge, it was difficult to mix LPOC by means of a normal
mixer. Therefore, in stead of manufacturing with a mixer, abinder was splayed on the upper surface of the concreterubbles of every layer. By the experiment, it was confirmedthat LPOC was producible by the spraying method.

Compressive strength ratioof porous concrete. Therefore, a simple geometrical model
concerning compressive strength was proposed ， and thecalculated values were compared with the experimental ones.
From the result, it is possible to predict the experimental
results by means of the proposed model. Especially, in thecase of direction 1, orthorhombic lattice array showed an
excellent presumption accuracy ( Fig. 2 ) .

Theoretical formula F =π（2tr＋t2）
2 √3 r2●：Theory ( Direction 1 )◇：Theory ( Direction 2 )△：Theory ( Direction 3 )○：Experiment（aggregate size: 0.6 - 400 mm）

Derection 3

1.2
1.0
0.8
0.6
0.4
0.2

F

0.0
0

10

20

30

40

Void ratio (％)

Fig. 2 Comparison between theoretical value[ 1 ] Akihiro Maegawa, Shigemitsu Hatanaka, Naoki Mishima andand experimental value
Yukihisa Yuasa, Fundamental study on manufacturing method( In the case of orthorhombic lattice array )of large particle porous concrete using concrete rubble andapplicability as fishing bank, Journal of Structural and Construction Engineering, AIJ, No. 589, 43-48, Mar., 2005
[ 2 ] Akihiro Maegawa, Shigemitsu Hatanaka, Naoki Mishima and Yukihisa Yuasa, Fundamental study on manufacture oflarge particle porous concrete using concrete rubble, Proceeding of Japan Concrete Institute, Vol. 26, No. 1,
1455-1460, 2004
[ 3 ] Akihiro Maegawa, Shigemitsu Hatanaka, Naoki Mishima and Moe Kuroda, Fundamental study on bending fracture
properties of large particle porous concrete, Proceeding of Japan Concrete Institute, Vol. 27, No. 1, 1273-1278, 2005
[ 4 ] Akihiro Maegawa, Toshihiko Shakouchi, Yukihisa Yuasa, Naoki Mishima and Shigemitsu Hatanaka, Research on
stability of fishing bank made by large particle porous concrete placed under sea, Journal of Technology and Design,
AIJ, No. 22, 53-58, Dec., 2005
[ 5 ] Akihiro Maegawa, Akira Yamamoto, Naoki Mishima and Shigemitsu Hatanaka, Experimental study on manufacturingmethod of small particle size porous concrete, Proceeding of Japan Concrete Institute, Vol. 28, No. 1, 1397-1402,
2006[ 6 ] Akihiro Maegawa, Shigemitsu Hatanaka, Naoki Mishima and Akira Yamamoto, Experimental study on fundamental
properties of small particle size porous concrete, Cement Science and Concrete Technology, No. 60, 264-270, 2006

Flow Analysis of Particle Laden Two-Phase Jet Flow from
Rectangular Nozzle and Its Application to Micro-Blasting Process
Masaki Sugimoto*(Fluid Engineering, Division of Systems Engineering)

Keywords: Rectangular jet flow, Gas-solid two phase flow, Micro-blasting
1. Introduction
Micro-blasting process, MBP by particle laden impinging air jet flow with the order of 10mm particles has
been increasing in order to process brittle materials, such as glass or ceramic, as the fine process for highaccuracy and productivity. In MBP, a circular nozzle is used commonly, yet it takes a large number of scanningtimes to process a wide area, and then a long machining time. Recently, use of a two-dimensional rectangularnozzle has been considered to process a large area with high efficiency. However, it is difficult to obtain uniform
and high accuracy for cutting depth over the nozzle width.
In this study, in order to improve the cutting performance of an MBP nozzle, a new rectangular nozzle witha large aspect ratio is proposed and the flow characteristics and cutting performance are examined precisely.
2. Experimental and results
Figure 1 shows the nozzles used. The first nozzle is considered in order to change a circular cross section toa rectangular one. The second nozzle has three bars in it to diffuse or spread the flow. The third and forth nozzles
are proposed to improve the cutting performance. Particle laden impinging jet is used to process the brittle
material, glass plate. Particle used is silicon carbide, SiC of 50% diameter dp50 is 25mm.
The particle behavior from each rectangular nozzle was measured by PIV. The cutting performance of
nozzle was examined by evaluating the shape of cutting cross section of grass plate. The cutting efficiency wasevaluated by comparison of cutting removals per supplied power with conventional circular nozzle.
Figure 2 shows the visualized flow pattern in the half right had the particle velocity vectors and in the half
left had the particle concentration for the free jet of each nozzle, respectively. In the case of 1st nozzle, particle
does not disperse to y direction, it concentrated at center of nozzle, y/h=0. However setting the vane in nozzle,
particle dispersed to y direction, 2nd nozzle shows uniform concentration except edge of nozzle. The
concentration of 3rd nozzle has uniform all over the nozzle width, then the velocity vector becomes uniform, too.
The cutting performance of each nozzle is shown in Fig.3. All of the case, cutting depth correspond withparticle concentration, dense region becomes deep and thin region becomes shallow. Center of nozzle and near
the both side edge in y direction the cutting depths of 1st and 2nd nozzle are considerably. On the other hand, the
3rd nozzle has a high cutting performance, the accuracy of cutting depth of it is within 10% in the wide range at–70<y/h<70 since particle flow is uniform.
Figure 4 shows the cutting efficiency of 3rd nozzle with conventional nozzle. The 3rd nozzle has the largest(12 ﾟ )

100

10

60

(38 ﾟ )

1(a) 1st nozzle(b) 2nd nozzle
Fig. 1 Details of new rectangular nozzle*Sintobrator, Ltd.

100
(53 ﾟ)

160(c) 3rd nozzle

60

100

3-φ5

100

φ32

-50

0

y/h

50

-50

50

0

y/h

-80

50

50

50

100

100

100

150

150
x/h

(a) 1st nozzle

x/h

-30

-60

:50m/sec

16(a) 1st nozzle(b) 2nd nozzle

-60
0
-30

y/h

(b) 2nd nozzle(c) 3rd nozzle
Figure 2 Particle flow image and particle velocity vector

Removal weight, W [g/min]Cutting depth, d [mm]

-60
0

80

150
x/h

0
-30

0

(c) 3rd nozzle

3rd noz zle
p ressure ty p e nozz le(conventional)
suction ty pe noz zle(conventional)

12

8

4

0
-120

-60

0

60

y /h

Figure 3 Cutting profile of work-piece

120

0

5
10
15
Supplied power, Wa [kW ]20

Figure 4 Cutting efficiency

cutting efficiency under the same supplied power. The cutting rate of the 3rd nozzle is almost 1.2 times of theconventional circular nozzle, and moreover it has a larger processing speed. It results from reduction of particleparticle interaction and improvement of particle discharging from impinging spot at work-piece.
3. Conclusion
Main results are as follows,(1) The circular vane is effective to uniform the particle flow. Also, uniform particle flow can be obtained byexpanding the nozzle exit.
(2) The cutting accuracy shows a high quality of within 10% at －70 < y/h < 70 for the 3rd nozzle.
(3) The 3rd nozzle takes the largest removal weight compared to the conventional circular nozzles. For
example, cutting efficiency of the 3rd nozzle improves by 20% of the pressure type nozzle at the supplied
power of W = 10kW.
Author’s publications[1] M. Sugimoto, T. Shakouchi, K, Hayakawa, M. Izawa, Gas-Particle Two-Phase Jet Flow from Slot Nozzle
and Micro-Blasting Process, Trans. of Jpn. Soc. Mech. Eng., Series B, Vol.71, No.710 (2005), pp. 24502457.
[2] M. Sugimoto, T. Shakouchi, K, Hayakawa, M. Izawa, A study on Machining Characteristics of TwoDimensional Rectangular Abrasive Jet Machining Nozzle, Journal of Jpn. Soc. Precis. Eng., Vol.72, No.1
(2006), pp. 121-126.
[3] M. Sugimoto, T. Shakouchi, K. Hayakawa and M. Izawa, Gas-Particle Two-Phase Jet Flow from Slot
Nozzle and Micro-Blasting Process，JSME Int. Journal, Series B, Vol.49, No.3 (2006), pp. 705-713.
[4] M Sugimoto, T. Shakouchi, M. Izawa, Behaviors of Particle Laden Jet Flow and Machining Characteristics
of Micro-Blasting Process, Journal of Jpn. Soc. Abras. Technol., Vol.51, No.1 (2007), pp. 46-51.
[5] M. Sugimoto, T. Shakouchi, K. Hayakawa, M. Okazaki and M. Izawa, Particle Laden Impinging Jet Flowfrom Rectangular Nozzle and Abrasive Jet Machining，Proc. of Int. Conf. on Jets, Wakes and Separated
Flows, ICJWSF-2005, JSME No.05-201 (2005), pp. 325-330.
[6] M. Sugimoto and T. Shakouchi, Development of Rectangular Microblasting Nozzle, Proc. of Int. Conf. on
Precision Engineering, ICPE 11th, JSPE Publication Series No.5 (2006), pp. 181-186.

Investigations on a Vertical Upward Gas-Liquid Two-Phase Flowthrough Sudden Contraction or Expansion Pipe
Voutsinas Alexandros+(Fluid and Thermo Engineering Laboratory, Division of Systems Engineering)
Keywords: Bubbly two-phase flow, Sudden expansion-contraction pipe, Fluctuating flow, Flow control
1. The purpose
The aim of this research is to study the pipe flow characteristics when a gas-liquid two-phase flow passes
through a vertical abrupt contraction or expansion area. The effect of several parameters, such as Reynolds
number and volumetric gas flow rate, are considered. Drag and fluctuating phenomena are investigated andsimple plus economic flow control methods are proposed and their results are presented.

without ring
Lr/H= 1=2

3.E+03PSD

2. Related prior arts
Two-phase gas-liquid flow systems are frequently seen in
the industrial field. Several studies have been conducted sofar, such as Aloui et al. (1999), conducting experiments forthe singularity pressure drop, void fraction, wall shear stress,
bubble velocity and size, or Kondo et al. (2003) giving
valuable information concerning flow regime difference,
local void fraction flow observation on bubble deformation.
However, studies dealing with the fluctuating phenomena
and drag of the expansion /contraction area, and methods offlow control are scarcely seen in the bibliography.

Reduction
58%

2.E+03

1.E+030.E+000

5

10f [Hz]

15

20

D R%

3. Results and discussion
Fig.1 FFT analysis results
3.1 Sudden contraction pipe flow
(A =0.18, Re=2×104, αv=5%)
The present work is divided into two main parts. Part one,is the case of a sudden contraction and part two, theexpansion case. In both cases flow parameters were varied
A=0.18A=0.33and the effect on the flow pattern was investigated.
20
For the case of a sudden contraction, prior to and after thecontraction area, vortices generate, affecting significantly thedownstream flow. The flowing conditions affect the size and10length of it and fluctuating phenomena were observed. Whenthe flow control method was applied, it showed a remarkable
suppression of the vortex generation downstream, which is0closely related to drag and fluctuation phenomena. In the
0
5
10
a v%
following, fig.1 shows the comparative results from FFTanalysis conducted between the plain contraction and when aring shaped obstacle (proposed method) is attached in two
Fig.2 Drag reduction rate, Lr /H =1
different upstream positions. It is clear that the PSD value,(Re=1×104)
which represents the intensity of the fluctuation, is largely
suppressed when the flow control method is applied. This verifies that flow control is been achieved.
Figure 2 shows the results from the drag reduction measurements. Pressure distribution along the vertical axis
was measured and the pressure loss difference due to the contraction effect was investigated.
3.2 Sudden expansion pipe flow
After the positive results from the proposed flow control method, investigations on the sudden expansion casewere conducted. Similar to the contraction, on the downstream, vortex area is generated. In addition, due to thebuoyancy force of the bubbles a gas-free region exists just downstream the expanded area. The area of thisregion is strongly affected by the liquid velocity. In addition, even though an axisymmetric sudden expansion,4. Concluding remarks

f [Hz]

10
Re=10000=200008=30000

6
4
2
0
0

10

20

av %

30

Fig.3 Second peak frequency variation undervarious Re number and gas fraction

62%

plain expansi on
ri ng, Lr/H=3
step-ri ng

69%

57%

60

71%

80

PSD

small changes in the force balance within the vortex region
cause instabilities and fluctuating phenomena occur.
Two flow control methods were investigated. One, mounting
the ring shaped obstacle at several downstream positions,dividing the vortex region into two parts and reducing its
intensity. The second, installing a step-ring just after the
expansion in order to reduce the recirculation area size.
It was found that two dominant frequency peaks exist. In
addition, the second peak is strongly affected by the gas
fraction as shown in fig.3. Increasing the volumetric gas
fraction, it shows a tendency to shift its frequency to lower
values even under different Reynolds number.
In fig.4, the FFT result between the plain expansion and thetwo suggested methods is illustrated. The fluctuation amplitude
is considerably suppressed for both dominant peaks.
The drag reduction occurring due to the flow control method
was also calculated and results are presented in fig.5. Reduction
under all flow conditions was achieved, with values reaching up
to 40% at its maximum.

40

20

DR

1) From visualization results in both cases it was made apparent
0
0
10
20
that the vortex region is strongly affected by the Reynolds
f [Hz] 30
number and volumetric gas fraction.
Fig.4 FFT analysis results2) Fluctuation phenomena were investigated downstream the(Re =1×104, av=5%)
contraction/expansion. The dominant fluctuation frequency isaffected significantly by gas fraction. Two dominant peaks were40
observed up to av=20% for the sudden expansion, while justone was seen for higher values showing a tendency to shift its
30
frequency to lower values in all cases.
3) When a ring is mounted, the frequency value did not change.
20
This denotes that the geometric structure does not affect it. The
most effective ring position shifts downstream for all Re
10number when increasing av.
ring, Lr/H=4
step-ring4) Applying the flow control method, manages to suppress the0vortex generation and reduces the fluctuation intensity, reaching
0
10
20
30
a %
values up to 70% depending on the flow conditions5) The losses occurring due to sudden changes were
Fig.5 Drag reduction rate, DR, (Re=2×104)
investigated and by the proposed methods, drag reduction wasachieved. Positive results could be retrieved at all conditions with values reaching up to 40%.
v

Author’s Publications
1] A. Voutsinas., T. Shakouchi., k. Tsujimoto., T. Ando., “Flow Characteristics and Drag Reduction of Vertical Upward
Gas-Liquid Two-Phase Flow through Sudden Contraction Pipe”(In Japanese), Transactions of JSME, Series B, Vol.72,
No.720, (2006), pp.1888-1894.
2] A. Voutsinas., T. Shakouchi., J. Takamura., K. Tsujimoto and T. Ando., “Flow Characteristics and Control of Vertical
Upward Gas-Liquid Two-Phase Flow Through a Sudden Contraction Pipe” ,JSME-International Journal, Series B1,
November, (2006), No.06-4043.
3] A. Voutsinas., T. Shakouchi., J. Takamura., K. Tsujimoto., T. Ando., “Flow and Control of Vertical Upward Gas-Liquid
Two-Phase Flow through Sudden Contraction Pipe”, Proceedings of the ICJWSF, International Conference on Jets, Wakes
and Separated Flows, CD-ROM (2005), pp.307-312, Toba, Mie, Japan.
4] A. Voutsinas., T. Shakouchi., T. Tsujimoto., T. Ando., ”Gas-Liquid Two-Phase Oscillating, Fluctuating Flow throughVertical Sudden Contraction Pipe”, Transactions of JSME, Series B, Vol.72, No.717, (2006), pp.1131-1136.
5] T. Shakouchi., T. Nakamura., J. Kawaguchi., A. Voutsinas., K. Tsujimoto and T. Ando., “PIV Measurement of Flow
Characteristics Around Two-Tandem Square Cylinders in Gas-Liquid Bubbly Flow”, Proceeding of 4th International
Symposium on Measurement Techniques for Multiphase Flows, Sept. 10-12, (2003), Hangzhou, China.

A study of efficient downlink packet transmissions technique forCDMA cellular mobile communications
Abubaker Faraj Khumsi †Keywords: Mobile communications, CDMA, Downlink, Transmit power control, Soft handoff, Transmission
control, IP packets, Geographical unfairness, TPWC, MLPT, Downlink capacity
1. The purpose
The aim of this study is to improve the system capacity in CDMA cellular systems, especially in downlink
communications. The focus on the applications which do not have stringent delay requirements such as web
browsing and file transfer, the real time transmissions are not always necessary and time delay is tolerable
such as IP packets. Since the CDMA cellular system is the major system and succeeded in the commercial
application, the focus in this thesis is on this technique.
2. The related prior arts
The rapid growth in demand for mobile communications has led into intense research and development
effort towards a new generation of cellular systems. The new system must be able to provide quality ofservice (QoS), support a wide range of services and improve the system capacity. In the cellular systems, themobility causes dynamic variations in link quality and interference levels, sometimes requiring that a
particular user change its serving base stations. This change is known as a handoff. Usually, continuous
service is achieved by supporting handoff from one cell to another. The availability of cheap and reliablewireless Internet access will shift the service base traditionally found in mobile cellular networks toward
emerging wireless Internet service. This will result in significant demand being placed on both existing andnext generation cellular and IP networks. The key challenge in realizing an IP wireless network is how to
support soft handoff in order to minimize the packet losses and enhance the system performance. Although
the soft handoff is an effective way to increase channel capacity, reliability, and coverage of CDMA (CodeDivision Multiple Access) systems, it has some disadvantages especially in downlink channels. One of themain problems is that the simultaneous multiple base station transmissions will cause an increase in theinterference that affects other radio links, and consequently limits the downlink capacity. The conventional
scheme, Site Selection Diversity Transmission Power Control (SSDT), which realizes site selection
transmission diversity instead of the full site transmission diversity during soft handoff mode has also someshortcomings. The disadvantage of this scheme is that in bad channel conditions the larger transmitted power
is required to compensate the fading. As a result of this, the downlink interference increases.
3. Proposed methods abstract
This thesis focuses on the applications which do not have stringent delay requirements such as web browsing
and file transfer, the real time transmissions are not always necessary and time delay is tolerable. The aim of
this thesis is to enhance the system performance for the downlink IP packet transmission by improving thetransmission delay and geographical unfairness in transmission quality.
This thesis divided into three parts; in the first part we propose a transmission control scheme to support IP
packets transmission on CDMA wireless networks. The objective of this proposal is to improve the system
capacity with maintaining the minimum service disruption during handoffs of the mobile stations. In this
proposed scheme, if the channel condition is bad due to fading fluctuation, the base station delays the packet
transmission until better channel conditions are available. In other words, the base station sends the packets
only when the instantaneous attenuation from the transmitting base station is less than a threshold value. The
performance assessment of the simulation results proved that by applying this proposed scheme the systemthroughput would be efficiently improved.
_________________________________________† Division of Systems Engineering, Department ofElectrical and Electronic Engineering, Mie University

In the second part of this thesis; the location-based transmission power and window control scheme(location-based TPWC scheme) is proposed. The objective is to improve the transmission delay performanceand the geographical unfairness in transmission quality. In the conventional TPWC and multi-link packet
transmission (MLPT) scheme, the packets are sent only during a transmit time window according to the
propagation conditions. Besides, MLPT scheme reduces the delay caused by TPWC. However, by applying
this conventional scheme, there will be a considerable increase in the transmission delay under somecircumstances owing to the delay packet transmission. To avoid these shortcomings, the proposed scheme
offers the same transmission window size for all mobile stations in order to improve the delay and unfairness.
By using the soft handoff status as an indicator of the mobile locations, the base station individually adjust thetransmission window for each mobile station. By maintaining almost the same size of transmission window
for all the mobile stations, which is attained by adjusting the required threshold value for each mobile station,
the transmission delay performance and the geographical unfairness can be improved.
In the third part of this thesis; an adaptive transmission window control scheme based on soft handoff status
and traffic load is proposed. This scheme is proposed to avoid the shortcomings of our previous proposed
scheme (a location-based TPWC scheme). The location-based TPWC scheme improved the delay
performance and geographical unfairness, but the selected threshold values which constrain the transmission
window are not optimum for all traffic loads. In this proposed scheme the threshold values are controlled bytaking into account not only the locations of the mobile stations but also the traffic load. The advantage of thisscheme is that the optimum threshold values are chosen for each traffic load. The throughput and transmission
delay performances for the proposed scheme were evaluated by computer simulations. Based on thesimulation results, the proposed scheme can achieve the same throughput performance and improve the delay
performance compared with the conventional schemes. Furthermore, the performance comparison concludes
that the proposed scheme is better than the conventional ones for maintaining fairness of service to all mobile
stations at different traffic loads. Accordingly, the proposed scheme is effective for packet transmissions inCDMA cellular packet communication systems.
4. Concluding remarks
This thesis investigated the packet transmission control in downlink CDMA cellular systems. Several
techniques are proposed to cope with the recent trend in wireless communications. This trend is to provide themobile users with ubiquitous access to the Internet, mobility and handoff. Mobility management is an
important issue in the area of mobile communications. These proposed schemes reduced the service disruption
during soft handoff phase and improved the system capacity. The techniques proposed in this thesis contributeto the development of the future mobile systems.
5. Publications[1] Abubaker Khumsi, Kazuo Mori, Hideo KOBAYASHI, “Downlink Packet Transmission Control in Soft
Handoff Status on CDMA Wireless IP Networks,” ECTI-EEC Transactions on Electrical Eng., Electronics,
and Communications, Vol. 3, No.2, pp125-131, August 2005.
[2] Abubaker Khumsi, Kazuo MORI, Hideo KOBAYASHI, “Downlink Packet Transmission Control Basedon Soft Handoff Status in CDMA Cellular Packet Networks,” IEICE Transactions on Communications, Vol.
E90-B, No.5, to be published in May 2007.
[3] Abubaker Khumsi, Kazuo Mori, Hideo KOBAYASHI, “Soft Handoff Technique for Downlink IP Packet
Transmissions in Cellular Environments,” Proc. of ISCIT 2004, pp.1138-1143, Sapporo, Japan, Oct. 2004.
[4] Kazuo Mori, Khomsi ABUBAKER, Katsuhiro NAITO, Hideo KOBAYASHI, Hamid AGHVAMI,
Downlink Packet Transmission Control Based on Handoff Status in CDMA Cellular Packet Networks,” Proc.
of IEEE International Zurich seminar on communications, pp.62-65, Feb. 2006.
[5] Abubaker KHUMSI, Kazuo MORI, Katsuhiro Naito, Hideo KOBAYASHI, “Adaptive Transmission
Window Control for Downlink Packet Transmission in CDMA Cellular Networks,” Proc. of 2006
International Symposium on Intelligent Signal Processing and Communication Systems (ISPACS), pp. 743746, Tottori, Japan, Dec. 2006.

Study on Mitigation Methods of Non-Linear Distortion for OFDM
Based Wireless Communications Systems
Pisit Boonsrimuang †Keywords: OFDM, PAPR, Non-Linear channel, Inter-modulation noise, WLAN channel, Satellite channel
Clipping, DAR, IDAR and Dummy sub-carrier insertion.
1. The purpose
The purpose of this research achieves the high transmission data rate in the non-linear channel for
Orthogonal Frequency Division Multiplexing (OFDM) based wireless communications systems. The targetsof research are to achieve the best bit error rate (BER) performance with efficient usage of non-linear power
amplifier at the transmitter. In this research proposes the several solutions of non-linear distortion, which
could provide various practical solutions for next generation of multimedia wireless communications systems
employing the OFDM technique.
2. The related prior arts
One of the limitations of using OFDM technique is the larger Peak to Averaged Power Ratio (PAPR) of
its time domain signal. The larger PAPR signal would cause the severe degradation of BER performance andthe undesirable frequency spectrum re-growth both due to the non-linear distortion occurring in the non-linear
amplifier which is usually required at the transmitter in the wireless communications systems. The simple
solution to overcome this problem is to operate the non-linear amplifier at the linear region with taking theenough larger input back-off (IBO). However, this approach leads the inefficient usage of non-linear power
amplifier, and would lead a serious problem on battery consumption especially for the cases of mobile
terminal and portable wireless LAN terminal. In order to maximize the usage of power efficiency, the nonlinear amplifier is usually forced to operate at the near its saturation region. However this approach will leadto the severe degradation of BER performance and undesirable frequency spectrum re-growth for the larger
PAPR signal due to the occurring of inevitably higher non-linear distortion.
3. Proposed methods abstract
The PAPR problem is currently recognized as one of the essential research topics when employing
OFDM technique in the non-linear channel. Up to today, various kinds of methods were proposed in order tosolve the PAPR problem. These proposed methods could be mainly categorized into two types of methods; 1)
to reduce the PAPR to extend possible at the transmission side, and 2) to compensate the non-linear distortion
at the receiver side.
As for the first type of methods, there are three major proposals which can improve the PAPR
performance at the transmission side, that are the selected mapping method (SLM), partial transmit sequence
method (PTS) and dummy sequence insertion method (DSI). All of these methods can provide the better
PAPR performance by controlling the phase of data or dummy sub-carriers. The SLM and PTS methods
control the phase of data sub-carrier and the DSI method controls the phase of dummy sub-carriers at thetransmission side. First two methods are required to inform the phase information controlled for the data subcarriers to the receiver as the side information. On the other hand, the DSI method requests no sideinformation, because of using dummy sub-carriers, which could be discarded at the receiver. From this reason,
the DSI method could be realized with almost the same system efficiency and less overall system complexity
as compared with the PTS and SLM methods. Although the PAPR performance could be improved asincreasing the number of predetermined discrete phases, the computation complexity would increase due tothe necessity of larger number of iterations in the flipping algorithm.
From this fact, the conventional DSI method has a limitation on PAPR improvement when assuming thesmall number of predetermined discrete phases so as to keep the acceptable computational complexity.
To solve the above problem on the conventional DSI method, this research proposes two novel PAPR
reduction methods both of which can achieve the better PAPR performance with less computation complexity
than the conventional DSI method. The feature of the first proposed DSI method is to employ the timefrequency domain swapping algorithm in the determination of frequency data for dummy sub-carriers. In the
proposed method, the larger amplitude levels are detected in the time domain and these amplitude levels areconverted to the frequency domain to determine the frequency data for the dummy sub-carriers. These
procedures are repeated up to reaching either of the predetermined PAPR performance or the number of___________________________________________† Division of Systems Engineering,
Graduate School of Engineering, Mie University

predetermined iterations. The feature of second proposed method is to enable the determination of frequency
data for each dummy sub-carrier theoretically by using the continuous phase. From this fact, the proposed
method can achieve the better PAPR performance with less computation complexity than that for theconventional DSI method. From the computer simulation results, it is confirmed that these two proposed
method can achieve much better PAPR performance and BER performance than the conventional methodwith less computation complexity in the non-linear channel.
As for the second type of methods, the clipping method in conjunction with the Decision Aided
Reconstruction (DAR) method was proposed to compensate the clipping noise at the receiver side. The
clipping method is well known simple method to enable the reduction of PAPR performance by clipping thelarger amplitude level of transmission OFDM signal at the transmission side. Although the clipping method
can achieve the better PAPR performance, this method produces newly the clipping noise which causes thedegradation of BER performance. The DAR method was proposed to compensate the clipping noise in which
the clipping noise is reconstructed at the receiver side by using the decision data information. However the
DAR method can only mitigate the clipping noise and not for the non-linear distortion occurring in the nonlinear amplifier.
To solve the above problem on the conventional DAR method, this research proposes the Improved DAR(IDAR) method, which can mitigate both the clipping noise and inter-modulation noise. In the proposed
IDAR method, the characteristics of non-linear amplifiers are required to be known at the receiver for
mitigating the inter-modulation noise. This research also proposes the estimation method for AM-AM and
AM-PM conversions characteristics of non-linear amplifiers by using low PAPR (Peak to Averaged Power
Ratio) preamble symbols. This research presents two example systems which employ the proposed IDAR
method in the wireless LAN system and the broad band satellite communication systems. From the computer
simulation results, it is concluded that the proposed IDAR method can achieve the higher transmission data
rate and higher efficient usage of non-linear power amplifier with keeping the better BER performance evenin the non-linear channel.
4. Concluding remarks
This research proposed the solutions for PAPR problems. These proposed methods could be mainly
categorized into two types of methods; 1) to reduce the PAPR to extend possible at the transmission side, and2) to compensate the non-linear distortion at the receiver side.
First type, the common name can be called by PAPR reduction methods. The PAPR performance isrequired to reduce as more as possible before inputted to the HPA at the transmitter side. All of PAPR
reduction methods can improve both for the spectral growth and BER performance for OFDM technique inthe wireless communications. This research proposes two PAPR reduction methods by using dummy subcarriers. The proposed methods can achieve better PAPR and BER performances with less complexity thanconventional method. These methods suitable for low cost receiver because it have no required additional
circuits. The conventional OFDM receiver can work with these methods.
Second type, this research proposes the IDAR method that can work well both for WLAN and satellite
channel. In the simulation results show that the proposed method can improve BER performance and efficient
usage of HPA. The proposed method can achieve higher transmission data rate than single carrier modulationin the satellite channel by using multi-level QAM. To realize in the practical system, research also proposes
the non-linear characteristic estimation method that the non-linear amplifier characteristics in the IDAR
method is required to be known at the receiver.
The multi-level QAM have capability to improve the transmission data rate, however have very
sensitivity to the non-linearity of HPA such as 64QAM. The last chapter proposed combination between
PAPR reduction and IDAR methods. The proposed method can employ multi-level QAM in the satellite
channel. The simulation results show that the proposed methods show the higher data transmission rate andbetter BER performance than single carrier modulation.
As a conclusion of researches in this research, the proposed PAPR reduction methods and mitigation
methods of non-linear distortion noise could provide various practical solutions for next generation ofmultimedia wireless communications systems employing the OFDM technique.
5. Publications[1] Pisit Boonsrimuang, Kazuo Mori, Tawil Paungma and Hideo Kobayashi, “Proposal of Clipping and InterModulation Noise Mitigation Method for OFDM Signal in Non-linear Channel,” IEICE Transactions on
Communications, Vol.E88-B, No.2, pp.427-435, February 2005.
[2] Pisit Boonsrimuang, Pornpawit Boonsrimuang, Kazuo Mori, Tawil Paungma and Hideo Kobayashi,“Mitigation of Non-linear Distortion using PTS and IDAR Method for Multi-Level QAM-OFDM
system,” ECTI Transactions on Computer and Information Technology, Vol.1, No.2, pp 84-90, Nov 2005.
[3] Pisit Boonsrimuang, Katsuhiro Naito, Kazuo Mori, Tawil Paungma and Hideo Kobayashi, “PAPR
Reduction of OFDM Signal by Use of DSI Method with Time-Frequency Domain Swapping Algorithm,”
IEICE Transactions on Communication, Vol.E90-B, No.4, pp.1001-1006, April 2007.

Ergonomic Study on VDT Workstation Designfor Decreasing Workload
Shin Saito†(System Design, Division of System Engineering)
Keywords: VDT, Workstation, System design, Multi display, Display height, Workload

1. Introduction
In this study, the height and the number of the liquid crystal display (LCD) were discussed to propose theergonomic design for visual display terminal (VDT) workstation that is able to decrease the workload of VDToperators. There were two kinds of studies in the doctor’s thesis. The first study was to decide the appropriate
height of LCD, and the second study was to discuss the number of the LCD, with using multi display type VDT.
In order to design the ideal VDT workstation, these VDT operations were analyzed by work efficiency,
work posture, electromyogram (EMG) activities and subjective assessment in ergonomic aspects.
2. Appropriate height of LCD in VDT work
The appropriate heights of LCD, including a notebook personal computer (NPC), were discussed byanalyzing of work posture and EMG activities of neck muscles of VDT operators. The comparative experiments
were carried out in each experimental condition.
Results were as follows. It is necessary to separate the display of NPC from the keyboard and central
processing unit (CPU). The position of the NPC display prevents the head from inclining forward, while
maintaining an adequate viewing distance. The height of the display should be adjusted to each user with 1200mm as an upper limit. The appropriate height can be estimated by summing of the sitting height and the chairheight of each user. Moreover, a manuscript stand should be attached to the display in data entry tasks to prevent
poor posture.
3. Characteristics of multi LCD type VDTThe characteristics of multi display VDT was evaluated by the analyses of work performance, subjective
assessment, work posture, EMG activities, critical fusion frequency (CFF), NASA Task Load Index(NASA-TLX) and subjective symptoms of fatigue.
Results were as follows. Work efficiency of the single LCD was inferior compared to multi LCD and large
one. The multi display type VDT workstation was evaluated highly, in visibility and operability, than the single
display. A head moved wider for the system with multi LCD than those with a single LCD and large one. EMG
activities did not differ significantly in all experiment conditions. Fatigue and workload of multi LCD did not
differ significantly than single LCD.
4. Conclusion
The appropriate height of LCD and the multi LCD were examined to design the comfortable VDTworkstation for decreasing workload. These design requirements are going to obtain ergonomic data applicableto the guidelines and recommendation for VDT operations.
†: Mie Prefectural College of Nursing, JAPANDevelopment of A Flexible Structural Force Sensor
and Application to A Robot Hand
Naoki Saito(System Engineering, Division of Mechanical Engineering)
Keywords: Flexible structure, Force sensor, Robot hand,
1. The purpose
The aim of this research is to develop a flexible structural force sensor which is mounted on a surface of afingertip of a robot hand and to examine feasibility of new dexterous tasks by the sensor. The sensor’s flexibility
is modeling on a human finger cushion. This research describes a structure and fundamental characteristics ofmeasurement of the flexible sensor and several applications about object grasping by a robot hand which
installed the sensor on its finger.
2. Outline of the flexible structural force sensor
In this research, we propose a new type of force
sensor. This sensor can measure a contact location fC andfour-axis contact forces which are consisted of a verticalforce FZ, two-axis shearing forces FX, FY and a rotational
moment around the contact point MC. This sensor
structure is shown in Figure 1. This sensor consists of twoparts: a contact part and a measurement part. The contact
part is an arch beam that is made of silicon rubber. Siliconrubber is an elastic material, so the contact part is easily
deformed.
On the other hand, the measurement part consists offour beams and a support part. When a contact force istransmitted from the contact part to the measurement part,
Figure 1: Structure of the flexible force.
these beams are deformed depending on the amount ofcontact force or the direction of the contact force.
Four-axis contact forces are calculated in terms of the amount of these deformations which are measured bystrain gages eX1, eX3, eZ1 and eZ3 on the beam of the measurement part.
3. Characteristics of measurement
A contact location and four-axis contact forces
are determined with the use of some equations which
are derived based on Castigliano’s theorem and simple
analysis model of the contact part which is included
all of dimensions and material characteristics of thesensor. Validity of equations for determines theseinformation and measurement accuracy are confirmed
experimentally. One of the experimental results isshown as Figure 2. This result shows that desired
values and experimental results are advisableagreement.
In addition, this sensor can also calculate
four-axis deformations of the flexible contact part ofthe sensor. Deformations of the contact part are
determined with the use of measurement result of eachcontact force.

Exp. condition
fC=0deg, FX=FY=0N, FZ=2N MC=-13.5Nmm.

Figure 2: Experimental results of a contact locationand four-axis contact force measurement.

4. Applications
4-1. Object grasping against translational and rotational slippage
The issue of translational and rotational slippage that occurs when a robot hand grasps an object is
discussed. To grasp an object against each directional slippage is important for a dexterous task realization. In
this application, a sufficient conditional expression to grasp an object against translational and rotational slippageis presented. This expression is derived based on the consideration of a contact area between a flexible contact
part and object’s surface such as Figure 3 which occurs the sensor is pressed to the object. An effectiveness ofthe expression is confirmed experimentally with the use of a translational gripper where the flexible contact
sensor is mounted. Figure 4 shows experimental result of object grasping. From this result, it is confirmed thatgrasping force is increasing with increasing translational force and rotational moment. On the other hand,
object’s position does not change approximately. Therefore, it is confirmed that the robot hand can grasp against
translational and rotational slippage.

Figure 3: A contact area between the flexible
contact part of the sensor and an object.

Figure 4: Experimental results of object grasping
against translational and rotational slippage.

4-2. Object placing with a low impulsive force
A placing motion of a grasped object by a robot hand equippedwith the sensor on the finger is discussed. The aim of this application isto realize a quick placing motion with the impulsive force between theobject and the floor small. To derive a dynamic model of the motion, we
consider the deformation of the flexible sensor occurring when the handgrasps the object. The dynamic model represents the relationshipbetween the impulsive force and an approaching trajectory of the robot
hand to the floor. From this model, we can obtain the trajectory of thehand which ensures that the impulsive force is less than the object’s own
weight.
Figure 5 and Figure 6 show that the experimental results of placing
motion. From Figure 5, it is confirmed that the robot hand position isapproximately same as a desired trajectory which is obtained by themodel, and the sensor deformation is decreased. From Figure 6, it isconfirmed that a reaction force is less than the object’s weight.
Therefore, it is confirmed that the hand equipped the flexible sensor can
put down the grasping object with low impulsive force.

Figure 5: Experimental result of robot
position and sensor deformation.

Figure 6: Experimental result ofreaction force from the floor.

5. Conclusion
In this research, we proposed a new type of force sensor which has a flexible contact mechanism. Thissensor can measure a contact location and four-axis contact forces.
An effectiveness of this flexible structural force sensor demonstrated through two applications. These
results of application suggest that the robot hand which has a flexible contact part on its finger surface canperform dexterous tasks.
Author’s publications (in part)[1]N. Saito, T. Satoh, S. Kajikawa, H. Okano, Development of Four-Axis Contact Force Sensor with a Flexible
Contact Part, IEEJ Trans. SM, Vol.125, No.1, pp.7-14, 2005.
[2]N. Saito, T. Satoh, H. Okano, Grasping Force Control in Consideration of Translational and Rotational
Slippage of an Object by a Flexible Contact Sensor, IEEJ Trans. EIS, Vol.126, No.8, pp.990-996, 2006.
[3]N. Saito, T. Satoh, Y. Suzuki, H. Okano, Placing Motion of an Object by a Robot Hand with a Flexible Sensor,
Trans. of JSME C, Vol.73, No.725, pp.206-213, 2007.

Estimation on The Amount of Movement Lines of Material
Transportation with Systematization of The Information and Material
Management- A study on the facility management of hospital architecture Shinsuke Kawai*(Arichitectural Design and Planning, Division of Systems Engineering)
Keywords: hospital architecture, facility management, connection with departments, ordering system,
material management and supply systems, model, amount of movement lines,
material transportation, estimation, nursing activities locations, gradual composition
This paper aims at the hospital architecture offering essential medical services. The foremost tasks of
hospital architecture are to contribute to the sound management and to maintain the high functional levels in thelong run. To settle these issues, facility management is thought to be quite effective.
Medical services by doctors and nurses include not only treatment directly for patients but also consequential
chores. Increasing the latter means decreasing the quality of medical treatment. Under the financial pressure,hospital managers must save life cycle cost. When hospital architecture is planned, flows of people, articles andinformation are considered to be indicators. These flows need to be investigated closely, and to be improved sothat medical staff can take more time in the direct medical treatment. They have to reduce the burden of movingaround, to diminish the stress of indirect consequential chores, and to save and simplify their labors. It can bealso successful in financial management. Targeting material transportation, the medical information system andthe material management and supply system, to systematize the movement of materials and information enables
hospitals to control effectively. This means achieving greater efficiencies in management and medical treatment.
The purpose of this study is as follows; one is to organize a hospital supporting system and to suggest whatto evaluate; another is to estimate the amount of movement lines, which have been made more efficient by ahospital supporting system and a department arrangement. Some previous work have reported the flows ofmaterials, the movement of the transporter, the amount of transportation, the department arrangements, and theestimation of plan type and so on. Generally speaking, the contents of the most previous work can be divided
into hospital management and architectural space. In this study, the medical information system and material
management and supply system are treated systematically. Considering the movement lines as an indicator,
hospital management and architectural space are treated comprehensively in the Chapter 1.
In the first half of the second chapter and in the third chapter, the first purpose of this study is achieved. One
example of hospital supporting system, i. e. ordering system has been investigated in the survey questionnaires.
The results of the survey show that the interlock of the doctor's processing order and accounting procedure isdistinguishing. Two difficulties are still left; how to make it a standard and how to abolish many kinds of tickets
and vouchers. The evolution items concerning what to introduce first are as follows; first, where many systems
are adopted-[B: prescription and checkup], [C' : B plus medical radiation], [C: C' plus injection], and [E: C plus
interlocking accounting]; second, how many tickets and vouchers are eliminated-[I: many], [II: few], and [III:none]. A pattern of [E-III] is the most systematic. Also interviews are conducted at four hospitals. As a result, it
is revealed that plural systems are introduced. The types of the material management and supply system aredepend on the combination of supply methods, management methods, and planar forms of supplying section.
From the actual conditions of those four hospitals, three typical models are derived. Making an allowance for the
introduction of the POS system and the material management at consumption of each model are simulated.
The latter half of the second chapter and the Chapter 4 treat the connection with departments. The patterns ofdepartment arrangement have developed with systemization of the medical information system. In the pattern of
E-III, characteristically the flow of job cards has replaced that of web orders. The common materials are pick up
from a flow chart of supplying operation and presented as a flow modeling. Examined qualitatively, the flows ofarticles and supplies in three models are estimated.
In the Chapter 5, the amount of transportation between each department is figured out. Another research ontransportation is conducted at the hospitals for three days. The results shows that the conveyance during the* Redeveloping Project Team on Univ.Hospital and School of Medicine, Mie Univ.

daytime on weekdays occupies 70 % of the whole, and that increase in scheduled transport makes the number oftransporting operations decrease. The rate of transportation load by the types of transportation is demonstrated.
The amount of movement lines of materials are calculated in the Chapter 6. That of each hospital is different,
but by establishing the medical information system, it can be reduced by about 25 to 50 %, and when thematerial management and supply systems are added, it can be cut by more 39 %, and when the POS system isadded, by 50 %. It is concluded that the material management and supply systems are effective in providing theeasy access, and that three models are successful in helping the medical staff to shorten the walking distance and
the POS system in helping the administrative staff. More two points are clarified; an intensive supplying
sections needs ingenuity of department arrangement in order to lessen the movement lines; and dispersive
supplying section should systematize the material management service.
The seventh chapter deals with the flow of staff and patients in hospital ward. Another inquiry has done at theforegoing hospitals. According to the results of the survey, it is determined that nurses walk around more than 4kilometer during the working hours. Most of the destinations are nurse station to exchange information,
document the nursing record, and so on. At the other destinations, they give patients treatment and talk with theirfamilies. And it has revealed that the locations of nursing activities are in a gradual composition. The
combinations of the materials and facilities are examined and the results show that the hospital supporting
system can decentralize working locations and shorten the movement lines. Moreover, even the existing
hospitals can make improvements.
The position in facility management (FM) of this research product is considered in the next chapter andestimated through the viewpoints of space and systems, which is one of the methods for evaluation of FM
productivity. One characteristic of hospital architecture is that a hospital has plural systems. Systematization not
only in each department but also between departments can make the flows and transportation easier. In eachdepartment, nursing service and the hospital supporting system are connected closely and the system should be
typical and customized. In the hospitals used for a long term, which needs more flexible systems, the model 3
and the POS system are promising. The space arrangement under the influence of systematization is the
positional relation between demanding spaces and supplying spaces. Scattering the supplying spaces can cut
down the amount of movement lines. On the other hand, concentrating them make it easier to upgrade systems
and renovate. In the case of the POS system, because out-of-hospital carriers are involved in the transporting
service, hospitals needs high security and the simple and quick access to demanding places from the outside.
Another example of FM in hospital architecture is functional estimation. In order to maintain the level ofsystems, periodical and uncomplicated checking is significant. When the material arrangement is not considered
for nursing service in a unit of a hospital ward, it might not be consistent with its hospital supporting system
and need to reassemble the whole system.
Three points are left for the future investigation; first, exploration in the United States about the developmentof the POS system; second, simulation of the changing service in a nursing ward by grading the hospital
supporting system; third, establishment of screening self-diagnosis in a practical way.
Author's publications[1]S.Kawai and S.Imai, The Developing Pattern of Ordering System's Introduction - A study on theinformational connection with departments of Hospital Architecture, Hospital Administration, Journal of theJapanese Society on Hospital Administration, Vol.40 No.2, pp.5-12, 2003.4[2]S.Kawai and S.Imai, A Study of Stock Management and Supply Systems through Modeling and Estimationon the Changes with Advancements in Information Technology, Hospital Administration, Journal of theJapanese Society on Hospital Administration, Vol.43 No.1, pp.23-34, 2006.1[3]S.Kawai and S.Imai, Evaluation of The Material Management and Supply Systems' Modeling through The
Material Flow - A study on the connection with departments of hospital vol.1, Journal of Architecture andplanning, Transactions of Architectural Institute of Japan, No.605, pp.55-62, 2006.7[4]S.Kawai and S.Imai, The Amount of Material Transportation through The Material Management and Supply
Systems - A study on the connection with departments of hospital vol.2, Journal of Architecture and planning,
Transactions of Architectural Institute of Japan, No.613, pp.45-52, 2007.3[5]S.Kawai, K.Katsuno and S.Imai, A Basic Study on Gradual Composition of Locations for Nursing Activities
in Wards in Relations to The Locations of Materials, Journal of Architecture and planning, Transactions ofArchitectural Institute of Japan, No.614, 2007.4

Motion tracking ofleft ventricular myocardium using ultrasoundWataru OHYAMA ∗Keyword: Ultrasonic diagnosis system, Myocardial motion tracking

1 Introduction
The echocardiogram is an essential tool for diagnosis and treatment of cardiovascular diseases. The useof echocardiograms for clinical diagnosis serves clinitians and doctors with multifold advantages including
noninvasiveness of diagnosis and real–time imaging. Recent technical innovations of ultrasonic diagnosis
equipment reclaim their new availabilities in clinics. These novel modalities require to track a speciﬁc
tissue of myocardium to keep their focus on a relevant tissue to evaluate. In this research, four novel
algorithms, that aim to contribute for these new availabilities of ultrasonic in diagnosis, are investigated.
2 Automatic Extraction of Left Ventricular Endocardium
Automatic extraction of left ventricular endocardium in echocardiograms is required for quantitativeevaluation of the functional performance of the left ventricle. In this section, a novel automatic extraction
method based on double thresholding for echocardiograms is proposed, and its eﬀectiveness and theaccuracy are evaluated. B-mode echocardiograms are ﬁrst binarized with a threshold determined by thediscriminant analysis for the gray level histogram. Then the binary images are contracted n times toremove small regions and to disconnect the region of cardiac cavity from the other false regions. Amongthe obtained regions which corresponds to the cardiac cavity is selected and dilated 2n times to createa mask which restricts the region of the second thresholding operation. The size and the location of the
cardiac cavity in the preceding frame are utilized to select the corresponding region. The masked image ofeach frame is binarized in the restricted area in the same way as in the ﬁrst thresholding operation. The
evaluation test is carried out using the scatter diagram of radius of contours extracted by two observers
and automatic extraction method. These results showed that the accuracy of the extracted contours wasfavorably compared to the accuracy of manually traced contours.
3 Automatic Tracking for Regional Myocardial Motion
In this section, an automatic tracking method for 2-D motion of Regional myocardial motion in highframe rate echocardiography by using the correlation method with connective multiple ROIs is proposed.
The proposed method is a combination method of maximization for correlation of brightness distribution
and optimization the geometric location of ROIs. A comparison test of tracking accuracies between the∗ Mie University

conventional method and proposed one is carried out. Results of these tests shows that proposed method
is able to derive more accurate trajectories from echocardiography than conventional method.
4 Local Myocardial Motion Tracking Using Ultrasonic Doppler Signal
In this section, a new method for automatically tracking the motion of local region in left ventricularmyocardium by means of ultrasonic pulsed Doppler signal is proposed. This method consists of a velocity
detection procedure based on correlation weighted mean instantaneous phase diﬀerence and a motion
tracking procedure employing a elastic model. Most of ultrasonic pulsed Doppler signals contain considerable amount of speckled noise, which causes detection error of velocity. The procedure of correlation
weighted mean velocity is aimed to reduce the velocity detection error, and the elastic model is used to
avoid the accumulation of the error to keep track of the motion of the myocardium in reasonable accuracy.
5 Local Myocardial Motion Tracking Using Ultrasonic RF Signal
This section describes a new algorithm to track identical myocardium accurately through a few cardiac
cycles. Hierarchical correlation is used with wide correlation window as ﬁrst and narrow one as second.
Then, to increase the reliability of the estimated position of myocardial motion, neighboring points areweighted by correlation coeﬃcients of ultrasonic waveforms between sequential pulses and are summed
up. This method was applied to normal and diseased hearts. Tracking error during one cardiac cycle
decreased to one-ﬁfth of simple method and the stability increased substantially. Regional myocardial
strains in the wall were estimated and color-coded using this tracking method. Result: this strain image
revealed clearly the diﬀerences of myocardial function between endocardium and epicardium, and thespeciﬁc characters of myocardial diseases. Moreover, the cyclic variation of integrated backscatters is
automatically obtained using this tracking method. This cyclic variation also indicated the myocardial
tissue characters of diseased heart in terms of ultrasonic backscatter intensity.
6 Conclusions
In this research, new modalities to track the myocardial motion from echocardiogram. The proposed
approaches are able to increase the performance of motion tracking. These methods contribute to provide
new applications of ultrasonic for diagnosing cardiovascular diseases.
Author’s publications
1. W.Ohyama, T.Wakabayashi, F.Kimura, S.Tsuruoka, and K.Sekioka: ”Automatic Extraction of Left
Ventricular Endocardium in Echocardiograms Using Double Thresholding Method (in Japanese)”:The Transactions of The Institute of Electrical Engineers of Japan: Vol. 121-C, No. 9, pp. 14481456
2. W. Ohyama, N. Ismail, T. Wakabayashi, F. Kimura, S. Tsuruoka, and K. Sekioka: ”Local Myocardial Motion Tracking Based on Correlation Phase Diﬀerence Method (In Japanese)”:IEICE
Transactions Vol. J86-A No.9 pp.917-928
3. W.Ohyama, M.Inami, T.Wakabayashi, F.Kimura, S.Tsuruoka, and K.Sekioka:

Automatic Track-ing for Regional Myocardial Motion by Correlation Method with Connecting Multiple ROIs (inJapanese) : The Transactions of The Institute of Electrical Engineers of Japan: Vol. 124-C, No.
10, pp. 2079-2086