〔主要な業績〕拡張された質点投射モデルの最適角とその近似解
坪井 一洋; 星野 貴弘, 筆頭著者, 日本応用数理学会
日本応用数理学会論文誌, 2023年09月25日, ［査読有り］

〔主要な業績〕硬式野球ボール型加速度センサにおける回転軸変化のシミュレーションによる再現
桧山 莉央,長岡 大志,坪井 一洋,野田 茂穂,姫野 龍太郎, 責任著者, 硬式野球ボール型加速度センサを用いた飛翔実験の結果から示唆された飛翔時に発生する回転軸の変化について考察した．センサユニット内の各センサが剛体的に配置された場合の各センサ位置での回転による加速度を定式化することでセンサ出力のシミュレータを作成した．ストレートとカーブの2球種に対して，センサ出力から得られた回転軸が時間的に変化する条件と回転軸をその初期値に固定して向きが変化しない条件でシミュレーションを行い，飛翔実験で得られたセンサ出力と比較した．その結果，前者の条件によってセンサから出力された加速度データが再現されることを確認し，ボール型センサの回転軸が飛翔時に変化することを示した．, 日本機械学会
日本機械学会論文集, 2021年07月25日, ［査読有り］

質点の投射角を再考する
坪井一洋, 日本物理学会
大学の物理教育, 2019年03月15日, ［査読有り］, ［招待有り］

加速度センサを内蔵したボールの飛翔状態での回転特性
長岡 大志; 桧山 莉央; 坪井 一洋; 野田 茂穂; 姫野 龍太郎,

Direct measurement of the rotational characteristics in flight of a baseball-like ball is performed by using a built-in acceleration sensor unit. This sensor unit consists of one three-axis sensor for low acceleration and three sensors for high acceleration. In the preliminary experiment, we determine the moment of inertia of the ball set up the sensor unit and it has different values in each rotation axis. At the same time, the accuracy of sensor outputs is confirmed. In order to investigate the variation of the rotation axis of the ball in flight, sensor outputs are simulated. As the result, measured outputs are reproduced successfully in the case of giving the time variation of the rotation axis obtained in the flight experiment, while the outputs are rather different under the condition of the fixed rotation axis. Furthermore, we study the temporal behavior of the translational acceleration of the ball and clarify the periodic variation of the direction of this acceleration component. These results suggest that the rotation axis of balls is varying in flight.

, 一般社団法人 日本機械学会
シンポジウム: スポーツ・アンド・ヒューマン・ダイナミクス講演論文集, 2019年

加速度センサを用いた飛翔中の野球ボールの回転特性計測
長岡 大志; 荻野 尚哉; 坪井 一洋; 野田 茂穂; 姫野 龍太郎,

We attempt to measure the rotational characteristics of a ball in flight by using a baseball-like ball with a built-in acceleration sensor unit which consists of one three-axis sensor for low G and three sensors for High G. Two approaches for obtaining the rotational speed and its axis of the ball from outputs of the sensors are formulated; one is based on rotational component in the acceleration outputs and the other is based on translational one. We employ in the present approaches the assumption that the time variation of rotation axis can be neglected because it is sufficiently slow. In flight experiments, we obtain nine time series data of acceleration for straight and curve ball trajectories, respectively. By using these experimental data, we estimate the rotational speed and rotation axis with the present two algorithms. As the results, it is shown that the rotation axis keeps almost constant during the flight in each algorithm's result and this means the validity of our assumption on rotation axis. The difference in the rotation axis obtained by the two algorithms ranges from 3.5 to 10.6 degrees in the straight ball data while it does from 3.9 to 18.8 degrees in the data of curve ball. Also, the averaged rotational speed in the nine trials is 12.56 rps in the straight data, and the difference between the algorithms indicates 0.10 rps (after pitching) and 0.17 rps (before catching). In the curve ball data, the averaged value and the difference between the algorithms are 19.22 rps, 0.10 rps (after pitching) and 0.14 rps (before catching), respectively. The present result shows the decrease of the rotational speed along flight trajectories, and it is probably the effect of aerodynamic torque.

, 一般社団法人 日本機械学会
日本機械学会論文集, 2019年, ［査読有り］

A Model for Merging Sections in Transportation Systems and Its Analysis
Takahiro Hoshino; Chuan Tian; Takafumi Komuro; Kazuhiro Tsuboi; Yoshio Hamamatsu, This paper proposes a model for merging sections in transportation systems and then analyzes traffic flow under this model. For this analysis, vehicles are represented by an idealized model based on a Markov chain model. Vehicles arrive at the merging section randomly in either of two inflow lines. When a vehicle arrives at the merging section, if no vehicle is present in the other line, then the vehicle must wait in order to connect with another vehicle on the other line before entering. In this merging section model, the queue length is restricted according to the travel time from the origin to destination. The mean queue length and mean queuing delay at the two inflow lines are explicitly derived here. The relation between the maximum queue length and the behavior of traffic is examined, and the obtained results clarify some characteristics of the proposed model., WILEY-BLACKWELL
ELECTRONICS AND COMMUNICATIONS IN JAPAN, 2016年09月, ［査読有り］

Merging Control Strategy with a Waiting-time Limit in a Personal Rapid Transit System　　　　　　　　　　　　　　　
Takafumi KOMURO; Kazuhiro TSUBOI; Takahiro HOSHINO; Yoshio HAMAMATSU, The Institute of Electrical Engineers of Japan
The International Conference on Electrical Engineering 2016, 2016年07月06日, ［査読有り］

A Metaheuristic Algorithm with Hybrid Insertion Procedure for the Traveling Salesman Problem
Takahiro Hoshino; Chuan Tian; Hisashi Kondo; Kazuhiro Tsuboi; Yoshio Hamamatsu
International Journal of Computer Science Issues, 2016年05月, ［査読有り］

角度に依存する初速度をもつ投射の最適角
坪井 一洋,星野 貴弘,浜松 芳夫, 筆頭著者, 一定速度で水平に移動する台車から質点が投射される力学モデルを考察した.この投射モデルは3個の変数をもち,スポーツ競技などの投射で見られる初速度の角度依存性を表すことが可能である.初速度と初速角の間に拘束条件をもつ本モデルによる初期条件下で飛距離が最大となる最適な投射角の理論的導出を試みた.その結果,ここで示した初等力学に基づく簡単な投射機構によって実際に観察される投射角が説明できることを示した., 一般社団法人 日本応用数理学会
日本応用数理学会論文誌, 2015年12月, ［査読有り］

個別輸送システム合流部の時間制約による制御戦略
小室孝文; 星野貴弘; 坪井一洋; 浜松芳夫, This paper deals with a merging section in a personal rapid transit system, and a new control strategy that allows a waiting time limit is proposed as a method to reduce the waiting queue. Vehicles arrive at a merging section from main and sub lines. Since vehicles on the main line have priority, a queue is formed on the sub line in order to avoid the collision of vehicles. Vehicles on the main line need to stop when the waiting time of the last vehicle on the sub line queue exceeds the limit. The relation between the waiting time limit and the behavior of traffic on the main and sub lines (characterized by the mean queue length and queuing delay) is derived explicitly by using a Markov chain model. Moreover, the control strategy for the queuing delay on each line is developed on the basis of the obtained relation. We compare the results obtained by using the present strategy with those obtained by using another strategy. The comparison shows the validity of the present strategy in cases where the probability of vehicle arrival on the sub line is low., Institute of Electrical Engineers of Japan
電気学会論文誌D（産業応用部門誌）, 2015年12月01日, ［査読有り］

輸送システム合流部のモデル化とその解析
星野 貴弘; 田 川; 小室 孝文; 坪井 一洋; 浜松 芳夫, This paper proposes a model for merging sections in transportation systems and then analyzes traffic flow under this model. For this analysis, vehicles are represented by an idealized model based on a Markov chain model. Vehicles arrive at the merging section randomly in either of two inflow lines. When a vehicle arrives at the merging section, if no vehicle is present in the other line, then the vehicle must wait in order to connect with another vehicle on the other line before entering. In this merging section model, the queue length is restricted according to the travel time from the origin to destination. The mean queue length and mean queuing delay at the two inflow lines are explicitly derived here. The relation between the maximum queue length and the behavior of traffic is examined, and the obtained results clarify some characteristics of the proposed model., Institute of Electrical Engineers of Japan
電気学会論文誌Ｄ（産業応用部門誌）, 2015年02月01日, ［査読有り］

軌道の特徴量を用いたボールの空力係数の推定
安田 海人,坪井 一洋,田中 晃平,宮嵜 武, 責任著者, スポーツで使用されるボールの空力係数を飛翔中の軌道形状データから求める方法を検証した．軌道形状の特徴量として，大域的形状を表わす量（飛距離や最大高さなど）と局所的に決まる量（曲率や捩率など）をそれぞれ考えた．前者を用いた方法では，特徴量の実測値と運動方程式による予測値の誤差を最小化することで空力係数を決定した．この方法でサッカーボールの空力係数の推定を行い妥当な結果が得られた．また，後者に基づく方法では空間曲線論を用いて局所的量と空力係数を関連付けた．この方法を卓球ボールの測定結果に適用することで，軌道に沿った空力係数の変化を評価することができた．, The Japan Society of Mechanical Engineers
日本機械学会論文集, 2014年06月25日, ［査読有り］

Aerodynamic Coefficients Estimation based on Local Geometry of Ball Trajectory　　　　　　　　　　　　　　　
Kazuhiro Tsuboi; Kaito Yasuda, 筆頭著者
Proceedings of the International Association of Computer Science in Sport Conference (IACSS2014), 2014年06月22日, ［査読有り］

Analysis of traffic congestion by considering merging on entry lines in transportation system
Takahiro Hoshino; Kazuhiro Tsuboi; Yoshio Hamamatsu, In a transportation system, a merging control strategy is necessary to avoid collisions between vehicles at a merging section with entry lines. In this study, under the assumption that each entry line is assigned time-independent priority, we analyze traffic congestion with a stochastic model. A quantitative estimation of the congestion is successfully obtained in terms of the average queue length and the average queuing delay. On the basis of analytical results, we propose a control strategy that changes the priority periodically. Using simulation, it is clarified that the proposed control strategy gives arbitrary values of the queue length and the queuing delay on each line. (c) 2012 Wiley Periodicals, Inc. Electron Comm Jpn, 95(8): 34-43, 2012; Published online in Wiley Online Library (). DOI 10.1002/ecj.11374, WILEY-BLACKWELL
ELECTRONICS AND COMMUNICATIONS IN JAPAN, 2012年08月, ［査読有り］

Application of Eulerian Approach to Simulate Hard Rime Accretion
Kazuhiro Tsuboi & Shigeo Kimura, 筆頭著者, We attempt to apply the Eulerian approach to simulation of ice adhesion on structures and to investigate macroscopic models for ice accretion processes. The methodology of the Eulerian approach in icing simulation is summarized and two technical difficulties are pointed out: the accuracy deterioration in representing the shape of accreted ice, and the stiffness of unsteady simulations. A method of adequately generating grid points overcomes the former difficulty while the latter remains an unsolved problem. Some computational examples obtained with the Eulerian simulation are presented. In the present simulations, we consider the hard rime accretion, i.e., an ice deposit formed from the impingement of super-cooled water droplets, and thus it can clarify the meanings of the macroscopic models for the local impingement efficiency and the growing direction of ice. The shapes of accreted hard rime on a stationary airfoil and the accretion rates are simulated, and they are compared with experimental results. The computational results show the importance of the velocity component of uniform flow direction in the incoming droplet velocity. Also, the local impingement efficiency along a cylinder oscillating perpendicular to a uniform air flow is obtained. The present simulations demonstrate the potential of the Eulerian approach in the prediction of ice adhesion phenomena. (C) Koninklijke Brill NV, Leiden, 2012, TAYLOR & FRANCIS LTD
Journal of Adhesion Science and Technology, 2012年07月24日, ［査読有り］

抵抗と揚力の働く投射体の最大投射角
坪井 一洋, 筆頭著者, 抵抗係数と揚力係数が一定の空気力が働く投射体の軌道を質点の運動で近似し，その最大投射角について考察した．速度に比例する空気力の働く質点の運動方程式から最大投射角を抵抗係数と揚力係数の摂動解として第２次項まで求めた．また，速度の自乗に比例する空気力に対する運動方程式から最大投射角の数値解を求めた．これらから抵抗と揚力が最大投射角に与える影響を評価した結果，空気力により最大投射角は減少する，減少の程度は抵抗より揚力の方が大きい，最大投射角は抵抗係数と揚力係数の2次式によって近似できる，などの結論が得られた, The Japan Society of Mechanical Engineers
日本機械学会論文集Ｃ編, 2012年06月25日, ［査読有り］

Application of Eulerian approach to simulate hard rime accretion
Kazuhiro Tsuboi; Shigeo Kimura, We attempt to apply the Eulerian approach to simulation of ice adhesion on structures and to investigate macroscopic models for ice accretion processes. The methodology of the Eulerian approach in icing simulation is summarized and two technical difficulties are pointed out: the accuracy deterioration in representing the shape of accreted ice, and the stiffness of unsteady simulations. A method of adequately generating grid points overcomes the former difficulty while the latter remains an unsolved problem. Some computational examples obtained with the Eulerian simulation are presented. In the present simulations, we consider the hard rime accretion, i.e., an ice deposit formed from the impingement of super-cooled water droplets, and thus it can clarify the meanings of the macroscopic models for the local impingement efficiency and the growing direction of ice. The shapes of accreted hard rime on a stationary airfoil and the accretion rates are simulated, and they are compared with experimental results. The computational results show the importance of the velocity component of uniform flow direction in the incoming droplet velocity. Also, the local impingement efficiency along a cylinder oscillating perpendicular to a uniform air flow is obtained. The present simulations demonstrate the potential of the Eulerian approach in the prediction of ice adhesion phenomena. © 2012 Copyright Taylor and Francis Group, LLC.
Journal of Adhesion Science and Technology, 2012年03月01日, ［査読有り］

Operation of Vehicles in AGV System and Method for Analyzing Operation
Takahiro Hoshino; Kazuhiro Tsuboi; Kazuo Yamanaka; Yoshio Hamamatsu, It is important to discuss the behavior of traffic at merging sections for the design and realization of an automated guided vehicle system (AGVS). In this study, we deal with a merging section of the AGVS by considering the time limit for the merging of vehicles. Near the merging section, the traffic flow in one direction forms a queue to avoid collisions between vehicles. We present an improvement of the control strategy proposed in previous research. However, since it is difficult to analyze the mathematical model of the merging section, analysis can be performed only for the case in which the time limit is set to the least values. We also propose a method for analyzing the model for an arbitrary time limit by using a linear recurrence relation. We discuss quantitatively the influence of the time limit on the average queue length and the detour rate. (C) 2011 Wiley Periodicals, Inc. Electr Eng Jpn, 178(1): 56-68, 2012; Published online in Wiley Online Library (wileyonlinelibrary.com). DOI 10.1002/eej.21190, WILEY-BLACKWELL
ELECTRICAL ENGINEERING IN JAPAN, 2012年01月, ［査読有り］

助走の効果を含む投射モデルのスポーツ競技への応用
坪井 一洋; 小長谷 康明, 筆頭著者, スポーツ競技で見られる助走を伴う投射に対する力学モデルとして３変数モデルを考察した．このモデルには投射時の水平速度，投射速度，投射角の3変数が含まれ，これらによって人体の体幹部と肢体の動きに付随する最小限の自由度を表現される．３変数投射モデルから初速度と初速角の間の関係式を導出し，それが助走を伴う投射に適した関数形であることを確認した．砲丸投げと走り幅跳びの実測データによる拘束条件と比較し，３変数モデルから最適初速角や最適初速度が推定できることを示した．このことから３変数モデルは初速度と初速角の間の拘束条件を含んだ投射モデルと考えることができ，助走を伴う投射モデルに対する力学的根拠を与える．, The Japan Society of Mechanical Engineers
日本機械学会論文集Ｃ編, 2011年10月25日, ［査読有り］

３自由度投射モデルによる砲丸投げ投法の特徴付け
坪井一洋, 筆頭著者, 助走の効果を含む３自由度投射モデルによって砲丸投げ投法の力学的差異について考察した．このモデルは砲丸の押し出し速度，押し出し角と助走の効果を表わす身体の移動速度からなる3個のパラメータで表現される．投射が最適な角度で行われたと仮定して，実測値からこれらのパラメータ値を推定した．その結果，押し出し速度と身体移動速度の間には投法によらず強い負の相関があり，これらの速度の間に緊密な関係があることを確認した．また，グライド投法と回転投法では身体移動速度の分布や最適押し出し角に違いがあることがわかった．このことは３自由度投射モデルによる砲丸投げ投法の力学的な特徴付けの可能性を示唆する結果と考えられる．, The Japan Society of Mechanical Engineers
日本機械学会論文集Ｃ編, 2011年03月, ［査読有り］

輸送システムにおける複数流入路での合流に伴う輻輳現象の解析
星野 貴弘; 坪井 一洋; 浜松 芳夫, In a transportation system, a merging control strategy is necessary to avoid collisions between vehicles at a merging section with entry lines. In this study, under the assumption that each entry line is assigned time-independent priority, we analyze traffic congestion with a stochastic model. A quantitative estimation of the congestion is successfully obtained in terms of the average queue length and the average queuing delay. On the basis of analytical results, we propose a control strategy that changes the priority periodically. Using simulation, it is clarified that the proposed control strategy gives an arbitrary values of the queue length and the queuing delay on each line. © 2011 The Institute of Electrical Engineers of Japan.
電気学会論文誌Ｄ（産業応用部門誌）, 2011年, ［査読有り］

A mathematical solution of the optimum takeoff angle in long jump
Kazuhiro Tsuboi, We successfully derive the optimum takeoff angle in long jump based on the maximization problem of the flight distance. The takeoff model proposed here includes three parameters: the horizontal speed of the centre of mass, the takeoff speed and the takeoff angle. The optimum takeoff angle is determined explicitly as a function of these parameters. The estimation of the optimum angle shows good agreement with the measured data of Mike Powell and Carl Lewis, which demonstrate the effectiveness of the present model and its solution. (C) 2009 Published by Elsevier Ltd., ELSEVIER SCIENCE BV
ENGINEERING OF SPORT 8: ENGINEERING EMOTION - 8TH CONFERENCE OF THE INTERNATIONAL SPORTS ENGINEERING ASSOCIATION (ISEA), 2010年, ［査読有り］

自動搬送システムの運行制御とその解析手法　　　　　　　　　　　　　　　
星野 貴弘,坪井 一洋,山中 一雄,浜松 芳夫
電気学会論文誌Ｄ, 2009年, ［査読有り］

着地時の重心低下を考慮した走り幅跳びの最適踏切角
坪井一洋, 筆頭著者, In long jump, it is possible to estimate a distance between takeoff and landing by assuming that the center of mass of a long jumper in aerial phase moves along a trajectory of projectile. Based on the assumption, we obtain an approximate trajectory of the jumper as a perturbed solution of flat aerial path, in which the effect of air drag is included as the linear form. Also, a maximization problem is formulated from the trajectory and the optimum angle of takeoff is theoretically derived. Then, we introduce two parameters, the horizontal velocity of the center of mass and the takeoff velocity into the present model of takeoff. In particular, the solution of the optimum angle includes the effect of vertical displacement at landing which affects significantly a distance of the aerial phase. The distance of the aerial phase and the takeoff angle are estimated in several cases of measured data. The results of the estimation show the effectiveness of the present solutions., 一般社団法人日本機械学会
日本機械学会論文集Ｃ編, 2008年07月, ［査読有り］

集団到着を考慮した個別輸送システム合流部の解析と制御戦略　　　　　　　　　　　　　　　
星野 貴弘,井上 康介,坪井 一洋,浜松 芳夫, This paper deals with a merging section of Personal Rapid Transit system with bulk arrival. Vehicles
arrive at a merging section from a main line and a sub line. Basically, vehicles on the main line have priority.
Therefore a queue is formed only on the sub line. In this situation, the queuing delay on the sub line may
become extremely long. We propose a control strategy whereby vehicles on the main line are usually not
stopped; vehicles on the main line are only stopped when the number of waiting vehicles on the sub line
exceeds some specific number. This specific number is a sort of threshold. Thus, queues are formed on the
main line and the sub line, respectively. The both of queuing delays are controlled by the threshold. The
optimal value of threshold is obtained from the result of analyzing the Markov chain model.
電気学会論文誌Ｄ, 2008年, ［査読有り］

Lateral Excitation between Dissimilar Orientation Columns for Ongoing Subthreshold Membrane Oscillations in Primary Visual Cortex
NAKAMURA Yuto; TSUBOI Kazuhiro and HOSHINO Osamu, Primary visual cortex (V1) is the first stage of visual information processing; detection of particularly oriented bars. These elemental visual features are sent to higher visual regions including V2; in which their combinations Such as "corners" and "junctions" are processed. A recent study, has demonstrated that corners and T-junctions could be processed even at the early visual stage (V1), raising the question of why the binding of bars by V1 is necessary. We stimulated a V1 neural network model, in which the so-called "orientation-columns" were connected through both excitatory and inhibitory synapses. The lateral excitatory connections contributed not only to binding paired bars constituting corners but also to making membrane oscillations near firing-threshold during ongoing (spontaneous) neuronal activity periods. This ongoing subthreshold neuronal state led to accelerating the reaction speed of neurons to paired bar-stimuli. The lateral inhibitory connections effectively enhanced the selective responsiveness of neurons to the stimuli. We suggest that coordinated lateral excitation and inhibition between orientation-columns in V1 could send angular information such as corners and junctions, presented to retina, rapidly to the next stage V2 for its full and precise analyses., SPRINGER-VERLAG BERLIN
Artificial Neural Networks - ICANN 2008, 2008年, ［査読有り］

An analytic solution of projectile motion with the quadratic resistance law using the homotopy analysis method
Kazuki Yabushita; Mariko Yamashita; Kazuhiro Tsuboi, We consider the problem of two-dimensional projectile motion in which the resistance acting on an object moving in air is proportional to the square of the velocity of the object ( quadratic resistance law). It is well known that the quadratic resistance law is valid in the range of the Reynolds number: 1 x 10(3) similar to 2 x 10(5) ( for instance, a sphere) for practical situations, such as throwing a ball. It has been considered that the equations of motion of this case are unsolvable for a general projectile angle, although some solutions have been obtained for a small projectile angle using perturbation techniques. To obtain a general analytic solution, we apply Liao's homotopy analysis method to this problem. The homotopy analysis method, which is different from a perturbation technique, can be applied to a problem which does not include small parameters. We apply the homotopy analysis method for not only governing differential equations, but also an algebraic equation of a velocity vector to extend the radius of convergence. Ultimately, we obtain the analytic solution to this problem and investigate the validation of the solution., IOP PUBLISHING LTD
JOURNAL OF PHYSICS A-MATHEMATICAL AND THEORETICAL, 2007年07月, ［査読有り］

２流体モデルによる粗氷形成の数値シミュレーション
坪井一洋,木村茂雄, 筆頭著者, Employing unsteady 2-fluid equation of multiphase flows as a mathematical model of flows including super-cooled droplets, we apply Eulerian approach to ice-accretion problem. A coupled system of equations for airstream and suspended water droplets is solved numerically with FDM and local impingement efficiency (mass flux of droplets on a body surface) is estimated. The temporal change of the shape of accreted ice is obtained from the local impingement efficiency and the renewal of grid distribution is performed. We also describe some difficulties on Eulerian approach in icing problems. These are mathematical property of the basic equations, the deterioration of representing the body shape due to the lack of grid point and stiffness due to large difference of characteristic time scales. In particular, the second difficulty is overcome by the method for grid interpolation. Hard-rime accretion on a NACA0015 air foil is simulated by means of the present procedure. In particular, the effect of two different definitions of local impingement efficiency and two different growth directions to the accreted ice is investigated. Comparison with an experimental observation shows that the simulated ice shape of local impingement efficiency with the uniform velocity direction and normal direction growth reproduces successfully the experimental results., 日本シミュレーション学会
シミュレーション, 2006年12月15日, ［査読有り］

風力タービンブレード上の着氷解析コードの開発　　　　　　　　　　　　　　　
木村茂雄,佐藤威,坪井一洋,森川浩司
日本雪氷学会誌, 2006年09月, ［査読有り］

２流体Saffman方程式の低Stokes数域における粒子相密度
坪井一洋, 筆頭著者, We investigate the density field of particle phase in low Stokes number region of 2-fluid Saffman equation which is a model of dilute particle flow on 2-fluid description with incompressible carrier flow. In the stagnation region of inviscid carrier phase, we find an exact solution of particle flow field and some features of the velocity and density fields are clarified. Furthermore, we reconsider Michael's criteria on the appearance of particle free layer. The order estimation of particle velocity field leads to the new criteria δ<κτ (δ,κ and τindicate the thickness of viscous boundary layer, the curvature and Stokes number) and computational results support the present criteria., 一般社団法人 日本応用数理学会
日本応用数理学会論文誌, 2006年03月, ［査読有り］

Computational Investigation on Flow Field of Particle Phase in 2-fluid Saffman Equation
TSUBOI Kazuhiro, 筆頭著者
AIAA paper, 2004年06月, ［査読有り］

Density Distribution in Stagnation Region of Saffman Eqaution for Dusty Gas
TSUBOI Kazuhiro, 筆頭著者, Proc. the 4th Int. Symp. on Computational Technologies for Fluid/Thermal/Chemical Systems with Industrial Applications, Volume 1
ASME PVP, 2002年08月, ［査読有り］

Aerodynamic characteristics of an iced cup-shaped body
KIMURA Shigeo; ABE ken'ichi; TSUBOI Kazuhiro; Tammelin Bengt; SUZUKI Kohji, The negative effects of ice accretion on wind speed measurements by a cup anemometer have been indicated by analyses of field test results amd already been studied numerically and experimentally by the authors of this paper. Implementation of calculations of the dynamic behavior of a cup anemometer with or without ice was carried out employing the fictitiously altered aerodynamic characteristics of an iced or clean cup based on NACA test results for a conical cup because no information covering the aerodynamic characteristics of an iced cup had been provided up to that point. In light of the absence of solid data, wind tunnel tests of the aerodynamic characteristics of iced cup-shaped bodies were performed. Imitated ice models were built and used for the test based on test results from an icing wind tunnel test conducted separately. We found that, depending on its amount, ice accretion in dry-growth conditions results in big penalties to the aerodynamics of cup-shaped bodies, while ice deposits of wet growth cause little effect, regardless of the thickness of ice. (C) 2001 Elsevier Science B.V. All rights reserved., ELSEVIER SCIENCE BV
Cold Regions Science and Technology, 2001年10月, ［査読有り］

Particle flow around an oscillating Circular Cylinder based on two-fluid model
TSUBOI Kazuhiro; KIMURA Shigeo, 筆頭著者, Proc. the 3rd Int. Symp. on Computational Technologies for Fluid/Thermal/Chemical Systems with Industrial Applications, Volume 1
ASME PVP, 2001年07月, ［査読有り］

Computational and Analytical Study on Particle Flow Field around a Circular Cylinder
TSUBOI Kazuhiro, 筆頭著者
AIAA paper, 2001年, ［査読有り］

伸長座標系における三重対角行列の条件数
坪井一洋, 筆頭著者, We investigate a condition number based on Enclidian norm of tridiagonal matrix which is obtained by discretizing Laplacian operator defined in one-dimensional general coordinate system. A particular stretching transformation of exponential function makes possible the analytical estimation of upper and lower bounds of the condition number. We also clarify the asymptotic behaviour of the condition number in large values of the transformation parameter, which shows the condition number comes close to the lower bound., 一般社団法人 日本応用数理学会
日本応用数理学会論文誌, 2000年, ［査読有り］

Computational Study on Local Impingement Efficiency of an Oscillating Circular Cylinder　　　　　　　　　　　　　　　
TSUBOI Kazuhiro; KIMURA Shigeo, 筆頭著者
AIAA paper, 2000年, ［査読有り］

振動物体への着氷速度に関する考察　　　　　　　　　　　　　　　
坪井一洋,木村茂雄, 筆頭著者
日本流体力学会誌, 1999年

Numerical Simulation of Ice Accretion on a Body with Droplet Flow Model
TSUBOI Kazuhiro; KIMURA Shigeo, 筆頭著者
AIAA paper, 1999年, ［査読有り］

Numerical Study of the Effect of Droplet Distribution in Incompressible Droplet Flows
Tsuboi Kazuhiro; Kimura Shigeo, 筆頭著者, We investigate numerically the effect of droplet distribution according to the diameter in incompressible droplet flows. The basic equations of droplet flows with the distribution effect is. presented and it is shown that the effect is expressed by a coefficient C. This coefficient gives non-dimensional form of the averaged relaxation time of droplets motion based on the characteristic time of main flow in employing the Stokes law. Computations with a finite difference method are performed to flows around a circular cylinder and a NACA0015 airfoil with the angle of attack 0°, which makes clear the variation of the flow field according to the value of C. In particular, the collision efficiency of droplet is obtained as a function of C from the computational results, which shows the drastic change in the range of 0.01<
C<
10. Our conclusion, therefore, is that the distribution of droplets becomes one of the important ingredients when the characteristic time of main flow is the same order as the relaxation time of droplets., American Institute of Aeronautics and Astronautics Inc.
AIAA paper, 1998年06月

装置周辺の空気の流れの可制御性とその制御法の一提案
白石昌武,出山浩行,柴田智和,坪井一洋, Air fluctuation around the fine equipment system seems to be a kind of extremely low velocity air flow and is one of the major causes in deteriorating system working accuracies. Unfortunately there is no active means to control air fluctuation at present. This paper describes the controllability of air flow around the equipment assuming that the air flow dynamics is approximated by a linear equation under the low velocity viscous flow condition. With the verification of its controllability, several experiments are performed to confirm the possibility of air flow control by using smoke particles. Specifically, smoke particles are stimulated by the optical radiation pressure of Ar laser beam. Their behaviors are regulated by the optical system such as a convex lens and an ND filter with the number of stimulated smoke particles. The experimental results basically show the possibility of air flow control using an optical radiation pressure., 公益社団法人精密工学会
精密工学会誌, 1998年05月05日, ［査読有り］

Formulation of Basic equations of Droplet flows with the Diameter Distribution　　　　　　　　　　　　　　　
筆頭著者
Proceediugs of The Third International Confenence on Fluid Mechanics, Beijing, 1998年

ゆっくり変動するモード間の共振的結合を伴う非線型偏微分方程式の数値シミュレーション
坪井一洋,ミハエル・トリベルスキー, 筆頭著者, A new approach to numerical integration of nonlinear PDE with two resonantly-coupling bands of slowly evolving modes, originated in the problem's symmetry, is developed. The bands are centered around a certain finite wavenumber and wavenumber equal to zero, respectively, that gives rise to a complicated spatiotemporal behaviour of solutions. The approach is based on reduction of the problem to a truncated set of coupled equations for Fourier-transforms of the order parameter supplemented, by exclusion of the slaved mode. Though the basic equation is stiff, an explicit time marching by Runge-Kutta method with some techniques is successfully performed. Application of the approach to the particular equation governing seismic waves yields solutions that display a new type of turbulence., 一般社団法人 日本応用数理学会
日本応用数理学会論文誌, 1997年

走り幅跳びの踏み切りにおける最適角について
坪井一洋, 筆頭著者, An approximate solution of orbit for long jump from takeoff to landing is obtained by considering the orbit as a projectile of a particle. The effect of air-drag is included as the linear form. Comparison with the exact solution of numerical integration shows the solution has good accuracy. Based on the solution, we can determine the optimum angle of takeoff successfully and esrimate the value from measured data., 一般社団法人 日本応用数理学会
日本応用数理学会論文誌, 1996年, ［査読有り］

New Scenario for Transition to Turbulence ?
Michael Tribelsky; Kazuhiro Tsuboi, Numerical study of the one-dimensional nonlinear partial differential equation, equivalent to that proposed [Recent Advances in Engineering Science (Springer-Verlag, Berlin, 1989)] to describe longitudinal seismic waves, is presented. The equation has a threshold of short-wave instability and symmetry, providing slow long-wave dynamics. It is shown that the threshold of the short-wave instability corresponds to a point of ''continuous'' (second order) transition from a spatially uniform state to a chaotic regime. Thus, contrary to the conventional scenarios, turbulence arises from the spatially uniform state as a result of one and the only one supercritical bifurcation., AMERICAN PHYSICAL SOC
Physical Review Letters, 1996年

天文台ドーム内の3次元流れに関する数値実験
坪井一洋,宮腰佳代,長澤幹夫,海部宣男,桑原邦郎, 筆頭著者
日本流体力学会誌「ながれ」, 1993年, ［査読有り］

Computational Study of unsteady viscous flow around a transversely and longitudinally oscillating circular cylinder in a uniform flow at high Reynolds numbers　　　　　　　　　　　　　　　
RAO P. M; KUWAHARA Kunio; TSUBOI Kazuhiro
Computational Mechanics, 1992年06月, ［査読有り］

Simulation of unsteady viscous flow around a lon-gitudinally oscillating cir-cular cylinder in a uniform flow.
P. M. Rao; Kunio KUWAHARA; Kazuhiro TSUBOI, A finite difference simulation method for the time-dependent viscous incompressible flow around a longitudinally oscillating circular cylinder at Reynolds numbers of Re = 1000 and 4000 is presented. The Navier-Stokes equations in finite difference form are solved on a moving grid system, based on a time-dependent coordinate transformation. Solution of the vortex street development behind the cylinder is obtained when the cylinder remains stationary and also when oscillating. Evolution with time of the flow field, by means of streamlines, pressure contours, and vorticity contours, is studied. The computer results predict the lock-in phenomenon, which occurs when the oscillation frequency is around double the natural vortex-shedding frequency. The time-dependent lift and drag coefficients are obtained by the integration of the pressure and shear forces around the body. The drag, lift, and the displacement relations are also discussed., BUTTERWORTH-HEINEMANN
Appl. Mech. Modelling, 1992年01月, ［査読有り］

An approach to incompressible flow problems with a massively parallel computer
TSUBOI Kazuhiro; MIYAKOSHI Kayo; KUWAHARA Kunio, 筆頭著者, Edited by K. G. Reinsch, W. Schmidt, A. Ecer, J.H user, J.Periaux, ELSEVIER SCIENCE PUBL B V
Proceedings of Parallel Computational Fluid Dynamics '91, Elsevier, 1992年, ［査読有り］

Incompressible flow simulations with solid models in rectangular grid system
TSUBOI Kazuhiro; MIYAKOSHI Kayo; KUWAHARA Kunio, 筆頭著者, Proceedings of Numerical Grid Generation in Computational Fluid Dynamics and Related Fields, ed. A.S.-Arcilla, J.H user, P.R.Eiseman, J.F.Thompson, ELSEVIER SCIENCE PUBL B V
Numerical Grid Generation in Computational Fluid Dynamics and Related Fields, North-Holand, 1991年, ［査読有り］

Incompressible Flow Simulation of Complicated Boundary Problems with Rectangular Grid System
TSUBOI Kazuhiro; MIYAKOSHI Kayo; KUWAHARA Kunio, 筆頭著者, UNIV TOKYO PRESS
Theoretical and Applied Mechanics, 1991年

INCOMPRESSIBLE-FLOW SIMULATION OF COMPLICATED BOUNDARY-PROBLEMS WITH RECTANGULAR GRID SYSTEM
K TSUBOI; K MIYAKOSHI; K KUWAHARA, Proceedings of Numerical Grid Generation in Computational Fluid Dynamics and Related Fields, ed. A.S.-Arcilla, J.H user, P.R.Eiseman, J.F.Thompson, UNIV TOKYO PRESS
THEORETICAL AND APPLIED MECHANICS, VOL 40, 1991年, ［査読有り］

Numerical Simulation of Unsteady Flow Patterns around a Vibrating Cylinder
T. Tamura; K.Tsuboi; K. Kuwahara, 26th Aerospace Sciences Meeting, ,AIAA-88-0128, pp. 1-11.
AIAA paper, 1988年01月

Computational Study of the Effect of Base Slant
TSUBOI Kazuhiro; SHIRAYAMA Susumu; OANA M; KUWAHARA Kunio, 筆頭著者, Proc. 2-nd Int. Conf. on Supercomputing Application and Automotive Industry, Cray Research, Inc., Cray Research Inc.
Supercomputer Applications in Automotive Research and Engineering Development, ed. C. Marino, North-Holland, 1988年

Numerical Study of Two-dimensional Vortex Street　　　　　　　　　　　　　　　
TSUBOI Kazuhiro; OHSHIMA Yuko, 筆頭著者
Numerical Methods in Fluid Mechanics, 1986年10月, ［査読有り］, ［招待有り］

Computational Study of the Effect of Base Slant　　　　　　　　　　　　　　　
K. Tsuboi; R. Himeno; S. Shirayama; K. Kuwahara, 筆頭著者, ASME/AIAA 4th Fluid Mechanics, Plasma Dy-namics and Lasers Conference
AIAA paper, 1986年05月, ［査読有り］

Merging of Two Dimensional Vortices by the Discrete Vortex Method
TSUBOI Kazuhiro; OHSHIMA Yuko, 筆頭著者, PHYSICAL SOCIETY JAPAN
Journal of the Physical Society of Japan, 1985年

書評「2点境界値問題の数理」
坪井 一洋
応用数理, 2007年03月
筆頭著者

着氷のシミュレーション
木村茂雄,坪井一洋
シミュレーション, 2001年03月, ［査読有り］

流体力学と並列計算
坪井一洋,桑原邦朗
コンピュートロール(コロナ社), 1992年
筆頭著者

Numerical Study of Two-dimensional Vortex Street
TSUBOI Kazuhiro; OSHIMA Koichi
The Institute of Space and Astronautical Science Report SP, 1985年03月, ［査読有り］
筆頭著者

〔主要な業績〕システムと微分方程式 第2版　　　　　　　　　　　　　　　
坪井一洋, 単著
三恵社, 2019年10月01日
9784866931159

システムと微分方程式　　　　　　　　　　　　　　　
坪井一洋, 単著
三恵社, 2011年04月01日
9784883618248

軌道データによる空力解析へのベイズ線形回帰の応用　　　　　　　　　　　　　　　
福岡 慎也; 杉原 皓大; 坪井 一洋
日本機械学会 スポーツ工学・ヒューマンダイナミクス部門講演会2024, 2024年11月16日
20241115, 20241117

軌道解析による空力評価への機械学習の応用　　　　　　　　　　　　　　　
関野 翔太; 坪井 一洋
日本機械学会 スポーツ工学・ヒューマンダイナミクス部門講演会2024, 2024年11月16日
20241115, 20241117

機械学習を用いたITRデータに基づく仕様分類　　　　　　　　　　　　　　　
関野 翔太; 坪井 一洋
第３２回 茨城講演会, 2024年08月30日
20240830, 20240830

ゴルフボールの三次元飛翔データに基づく非定常空力特性の評価　　　　　　　　　　　　　　　
岡﨑 健悟; 村上 航輝; 坪井 一洋
第３２回 茨城講演会, 2024年08月30日
20240830, 20240830

軌道解析に基づく飛翔物体の非定常空力特性の推定法　　　　　　　　　　　　　　　
村上 航輝,坪井 一洋
日本機械学会 2023 年度茨城講演会, 2023年08月18日, 日本機械学会関東支部，精密工学会
20230818, 20230818

機械学習による空力評価における特徴量の選定　　　　　　　　　　　　　　　
関野 翔太,坪井 一洋
日本機械学会 2023 年度茨城講演会, 2023年08月18日, 日本機械学会関東支部，精密工学会
20230818, 20230818

近似ベイズ計算による初期値問題のパラメータ推定　　　　　　　　　　　　　　　
福岡 慎也,坪井 一洋
日本機械学会 2023 年度茨城講演会, 2023年08月18日, 日本機械学会関東支部，精密工学会
20230818, 20230818

飛翔実験による空力特性の推定について
坪井 一洋
日本流体力学会年会2022, 2022年09月27日, 日本流体力学会
20220927, 20220929

飛翔軌道の幾何学量に基づくボールの非定常空力特性の推定　　　　　　　　　　　　　　　
森村 和馬; 坪井 一洋
日本機械学会シンポジウム：スポーツ工学・ヒューマンダイナミクス2021, 2021年11月13日, 日本機械学会
20211112, 20211114

野球ボール型加速度センサの出力シミュレーションによる飛翔時の加速度の検証　　　　　　　　　　　　　　　
桧山 莉央,坪井 一洋,野田 茂穂,姫野 龍太郎
日本機械学会シンポジウム：スポーツ工学・ヒューマンダイナミクス2021, 2021年11月13日, 日本機械学会
20211112, 20211114

近似ベイズ計算法による飛翔体の空力係数推定　　　　　　　　　　　　　　　
會澤 和希,坪井 一洋,竹田 晃人
日本機械学会 2021年度茨城講演会, 2021年08月21日, 日本機械学会
20210819, 20210819

軌道データに基づいた空力係数推定への近似ベイズ計算法の応用　　　　　　　　　　　　　　　
石塚 雅人,坪井 一洋,竹田 晃人,宮嵜 武
日本機械学会 シンポジウム： スポーツ工学・ヒューマンダイナミクス2020, 2020年11月14日, 日本機械学会
20201113, 20201115

ボール型加速度センサの出力シミュレーション　　　　　　　　　　　　　　　
桧山 莉央,長岡 大志; 坪井 一洋,野田 茂穂,姫野 龍太郎
日本機械学会 2020年度茨城講演会, 2020年08月21日, 日本機械学会
20200821, 20200821

軌道形状から推定した放物体の空力特性　　　　　　　　　　　　　　　
森村 和馬,坪井 一洋
日本機械学会 2020年度茨城講演会, 2020年08月21日, 日本機械学会
20200821, 20200821

加速度センサを内蔵したボールの飛翔状態での回転特性　　　　　　　　　　　　　　　
長岡大志,桧山莉央; 坪井一洋,野田 茂穂,姫野龍太郎
日本機械学会 シンポジウム： スポーツ工学・ヒューマンダイナミクス2019, 2019年10月27日, 日本機械学会

飛翔時のボール軌道の3次元計測と空力係数の推定　　　　　　　　　　　　　　　
三浦拓海,森村和馬,坪井一洋,梅津信幸
日本機械学会 シンポジウム： スポーツ工学・ヒューマンダイナミクス2019, 2019年10月27日, 日本機械学会

速度合成による質点投射モデルとその投射角　　　　　　　　　　　　　　　
坪井一洋
日本物理学会 2018年秋季大会, 2018年09月09日, 日本物理学会

加速度センサを用いた飛翔中の野球ボールの回転特性計測　　　　　　　　　　　　　　　
長岡 大志; 荻野 尚哉,坪井 一洋,野田 茂穂,姫野 龍太郎
Dynamics and Design Conference 2018, 2018年08月30日, 日本機械学会

飛翔実験によるボールの空力評価の誤差解析　　　　　　　　　　　　　　　
坪井一洋,渡辺卓馬,宮嵜武
日本機械学会 シンポジウム： スポーツ工学・ヒューマンダイナミクス2017, 2017年11月10日, 日本機械学会

軌道情報を用いた飛翔中のボールの空力係数推定　　　　　　　　　　　　　　　
渡辺卓馬,荒木田祐希,坪井一洋
日本機械学会 シンポジウム： スポーツ工学・ヒューマンダイナミクス2016, 2016年11月10日, 日本機械学会

拡張された質点投射モデルにおける最適角　　　　　　　　　　　　　　　
坪井;一洋
日本物理学会2016年秋季大会, 2016年09月15日, 日本物理学会

放物運動の最適角 －初等力学と最適化問題－　　　　　　　　　　　　　　　
坪井一洋
日本機械学会2015年度年次大会, 2015年09月15日, 日本機械学会

個別輸送システムにおける待ち時間を制限した合流制御戦略　　　　　　　　　　　　　　　
小室孝文,星野貴弘,坪井一洋,浜松芳夫
第22回電気学会東京支部茨城支所研究発表会, 2014年11月22日, 電気学会東京支部茨城支所

角度に依存する初速度をもつ放物運動の最大投射角　　　　　　　　　　　　　　　
坪井一洋,星野貴弘
日本応用数理学会2014年度年会, 2014年09月05日, 日本応用数理学会

軌道形状に基づく投射体の空力係数の推定法　　　　　　　　　　　　　　　
安田海人,坪井 一洋,田中晃平,宮嵜武
シンポジウム：スポーツアンドヒューマンダイナミクス2013, 2013年11月03日, 日本機械学会

サッカー競技のオブジェクト指向シミュレーション　　　　　　　　　　　　　　　
中村 翔吾,坪井 一洋
日本機械学会 2012 茨城講演会, 2012年08月24日, 日本機械学会

サッカーボールの空力係数の推定法　　　　　　　　　　　　　　　
安田 海人,坪井 一洋
日本機械学会 2012 茨城講演会, 2012年08月24日, 日本機械学会

２流体モデルによる着氷シミュレーションとその課題　　　　　　　　　　　　　　　
坪井一洋
第3回着氷研究会, 2012年07月18日

抵抗と揚力の働くボールの最大投射角　　　　　　　　　　　　　　　
坪井 一洋,鈴木 悠太
シンポジウム：スポーツ・アンド・ヒューマン・ダイナミクス2011, 2011年11月02日, 日本機械学会

助走を伴う投射モデルによる砲丸投げ投法の特徴付け　　　　　　　　　　　　　　　
坪井 一洋,小長谷 康明
シンポジウム：スポーツ・アンド・ヒューマン・ダイナミクス2010, 2010年11月05日, 日本機械学会

走り幅跳びにおける3自由度踏切モデルの検証　　　　　　　　　　　　　　　
小長谷 康明,坪井 一洋
シンポジウム：スポーツ・アンド・ヒューマン・ダイナミクス2010, 2010年11月04日, 日本機会学会

複数入力のある待ち行列モデルの解析　　　　　　　　　　　　　　　
星野貴弘; 坪井一洋; 浜松芳夫
日本応用数理学会2010年度年会, 2010年09月09日, 日本応用数理学会

A mathematical solution of the optimum takeoff angle in long jump　　　　　　　　　　　　　　　
TSUBOI Kazuhiro
8th Conference of the International Sports Engineering Association, 2010年07月15日, The International Sports Engineering Association

有限な保存量をもつ非粘性Burgers方程式の弱解　　　　　　　　　　　　　　　
李 菁,齋藤 将弘,坪井 一洋
日本機械学会関東支部第16期総会講演会, 2010年03月11日, 日本機械学会関東支部

Application of Optimization Technique to Solution of Incompressible Navier-Stokes Equation　　　　　　　　　　　　　　　
Kazuhiro Tsuboi; Kenji Sasajima
5th AIAA Theoretical Fluid Mechanics Conference, 2008年06月24日, 米国航空宇宙学会

最適化手法に基づく非圧縮性流れの計算法の平面Poiseuille問題への応用　　　　　　　　　　　　　　　
笹島健司,坪井一洋
第21回数値流体力学シンポジウム, 2007年12月21日, 日本流体力学会

着地時の重心低下を考慮した走り幅跳びの最適踏切角　　　　　　　　　　　　　　　
坪井一洋,高橋俊祐
ジョイント・シンポジウム2007 スポーツ工学シンポジウム／シンポジウム：ヒューマン・ダイナミクス, 2007年11月15日, （社）日本機械学会

Some Remarks on Basic Equations of Eulerian Approach in Icing Problems　　　　　　　　　　　　　　　
TSUBOI Kazuhiro; KIMURA Shigeo
第12回構造物の大気着氷雪に関する国際ワークショップ, 2007年10月10日

ある２流体方程式の分散相密度について　　　　　　　　　　　　　　　
坪井一洋
日本応用数理学会2007年度年会, 2007年09月17日, 日本応用数理学会

現代制御理論に基づく非圧縮性流れの新しい計算法　　　　　　　　　　　　　　　
笹島健司,坪井一洋
第26回日本シミュレーション学会大会, 2007年06月22日, 日本シミュレーション学会

最適化手法に基づく非圧縮性Navier-Stokes方程式の計算法　　　　　　　　　　　　　　　
笹島 健司; 坪井 一洋
第20回数値流体力学シンポジウム講演論文集, 2006年12月18日

空気抵抗が速度の２乗に比例する場合の放物運動の解析解とその性質　　　　　　　　　　　　　　　
山下真理子,藪下和樹,坪井一洋,岡畑豪
日本応用数理学会2006年度年会講演予稿集, 2006年09月17日

最適化手法に基づく定常非圧縮性流れの計算法　　　　　　　　　　　　　　　
坪井一洋,笹島健司
日本応用数理学会2006年度年会講演予稿集, 2006年09月16日

最適化手法を用いた非圧縮性流れの計算法　　　　　　　　　　　　　　　
笹島健司,佐藤陽介,坪井一洋
第25回日本シミュレーション学会大会発表論文集, 2006年06月23日

２流体モデルによる粗氷形成のシミュレーション　　　　　　　　　　　　　　　
坪井一洋,木村茂雄
第43回日本伝熱シンポジウム講演論文集, 2006年06月01日

ホモトピー解析法による放物運動の解析解 －空気抵抗が速度の２乗に比例する場合－　　　　　　　　　　　　　　　
山下真理子,藪下和樹,坪井一洋
第５５回理論応用力学講演会講演論文集, 2006年01月

回転型ＮＳ方程式による非圧縮性流れの計算　　　　　　　　　　　　　　　
坪井一洋
日本流体力学会年会2005講演論文集, 2005年09月07日

領域オブジェクトの分散化と非同期並列計算　　　　　　　　　　　　　　　
小林伊織,坪井一洋
第24回日本シミュレーション学会大会発表論文集, 2005年07月15日

ＴＶＤ差分法を用いた交通方程式の計算　　　　　　　　　　　　　　　
大塩晃之,坪井一洋
第24回日本シミュレーション学会大会発表論文集, 2005年07月14日

回転型の基礎方程式による非圧縮性流れの計算法　　　　　　　　　　　　　　　
坪井一洋
第54回理論応用力学講演会講演論文集, 2005年01月26日

非圧縮性流を主流とする２流体方程式の数値解　　　　　　　　　　　　　　　
坪井一洋,木村茂雄
日本流体力学会２００４講演論文集, 2004年08月11日

交通信号を伴う交通流の差分解　　　　　　　　　　　　　　　
冨沢哲朗,坪井一洋,仲田徹也
第23回日本シミュレーション学会大会発表論文集, 2004年06月16日

流体モデルによる交通流シミュレーションのためのオブジェクト指向モデリング　　　　　　　　　　　　　　　
坪井一洋,鈴木尊人
第23回日本シミュレーション学会大会発表論文集, 2004年06月16日

２流体Saffman方程式の数値計算　　　　　　　　　　　　　　　
坪井一洋,木村茂雄
第17回数値流体力学シンポジウム講演論文集, 2003年12月18日

２流体モデルによる物体への着氷シミュレーション　　　　　　　　　　　　　　　
坪井一洋,木村茂雄
日本流体力学会２００３講演論文集, 2003年07月28日

２流体方程式による物体への着氷シミュレーション　　　　　　　　　　　　　　　
坪井一洋,木村茂雄
第22回日本シミュレーション学会大会発表論文集, 2003年06月18日

流体モデルによる交通流シミュレーションへの領域オブジェクトの応用　　　　　　　　　　　　　　　
鈴木尊人.坪井一洋
第22回日本シミュレーション学会大会発表論文集, 2003年06月18日

分散させた領域オブジェクトによる並列計算　　　　　　　　　　　　　　　
三石智道,坪井一洋
第16回数値流体力学シンポジウム, 2002年12月16日, 日本数値流体力学会

２流体方程式による円柱への着氷シミュレーション　　　　　　　　　　　　　　　
坪井一洋,木村茂雄
第16回数値流体力学シンポジウム, 2002年12月16日, 日本数値流体力学会

2流体Saffman方程式による物体まわりの流れ計算　　　　　　　　　　　　　　　
坪井 一洋
日本流体力学会年会2002, 2002年07月23日, 日本流体力学会

物理学／電気と磁気　　　　　　　　　　　　　　　
2005年10月 - 2017年03月
茨城大学教養科目（自然／分野別基礎科目）

プログラミング演習　　　　　　　　　　　　　　　
茨城大学工学部システム工学科

コンピュータネットワーク　　　　　　　　　　　　　　　
茨城大学工学部システム工学科

シミュレーション　　　　　　　　　　　　　　　
茨城大学工学部システム工学科

情報処理概論　　　　　　　　　　　　　　　
茨城大学教養科目（情報関連科目）

1996年 - 2021年, 日本応用数理学会

1989年 - 2020年05月, 米国航空宇宙学会

日本機械学会

日本流体力学会

日本物理学会

センサー内蔵ボールに関する共同研究　　　　　　　　　　　　　　　
2020年04月 - 2022年03月

センサー内蔵ボールに関する共同研究　　　　　　　　　　　　　　　
2017年04月 - 2020年03月

野球ボール打球軌道シミュレータに関する研究　　　　　　　　　　　　　　　
2017年01月 - 2017年06月

着氷シミュレーションとその実用化に関する研究　　　　　　　　　　　　　　　
2013年09月

極地における新エネルギ導入によるエネルギ対策に関する研究　　　　　　　　　　　　　　　
2007年04月 - 2009年03月

ホモトピー解析法による非線形方程式の解法　　　　　　　　　　　　　　　
2005年06月 - 2008年03月

物体上の着氷雪に関する流体力学的研究　　　　　　　　　　　　　　　
2004年04月 - 2007年03月

物体への着氷現象の流体力学的研究　　　　　　　　　　　　　　　
2001年04月 - 2004年03月

振動物体への着氷予測に関する数値シミュレータの開発　　　　　　　　　　　　　　　
基盤研究(C)一般
2000年04月 - 2002年03月

物体への着氷現象の流体力学的研究　　　　　　　　　　　　　　　
1998年04月 - 2001年03月

振動翼上の着氷とその空力特性への影響に関する研究　　　　　　　　　　　　　　　
1997年04月 - 1998年03月

着氷したロータブレードの空気力学的特性に関する研究　　　　　　　　　　　　　　　
1996年04月 - 1998年03月

プロ野球「飛ばないボール問題」の議論が過熱 「縫い目が0.1ミリ高いだけで大きく影響」と専門家
株式会社朝日新聞出版, AERAdot., 2022年06月25日, インターネットメディア

博士論文審査委員会委員　　　　　　　　　　　　　　　
審査・評価
日本大学大学院理工学研究科, 2016年12月 - 2017年02月

日本機械学会関東支部第15期総会講演会　　　　　　　　　　　　　　　
企画立案・運営等
日本機械学会, 茨城大学水戸キャンパス, 2009年03月06日

第18回計算力学講演会　　　　　　　　　　　　　　　
企画立案・運営等
日本機械学会, 筑波大学春日キャンパス, 2005年11月19日

External examiner of Ph. D. thesis　　　　　　　　　　　　　　　
審査・評価
University of Manitoba, Canada, 2005年03月 - 2005年11月