Curriculum Vitaes

Makoto Nagaoka

  (永岡 真)

Profile Information

Affiliation
Professor, Faculty of Engineering Department of Mechanical Engineering for Transportation, Osaka Sangyo University

Researcher number
90394600
J-GLOBAL ID
202201010673897653
researchmap Member ID
R000035619

Papers

 50
  • 永岡 真, 業天 祐治, 齋藤 崇志, 薮下 広高
    自動車技術会論文集, 53(5) 904-909, Sep, 2022  Peer-reviewedLead authorCorresponding author
  • 業天 祐治, 薮下 広高, 永岡 真, 齋藤 崇志, 吉岡 雅也
    自動車技術会論文集, 53(1) 114-118, Jan, 2022  Peer-reviewed
  • Hirotaka Yabushita, Makoto Nagaoka, Yuji Gyoten, Masaya Yoshioka, Yuichi Mori
    SAE International Journal of Advances and Current Practices in Mobility, 4(2) 583-591, Sep 21, 2021  Peer-reviewed
  • Yoshihiro Nomura, Seiji Yamamoto, Makoto Nagaoka, Stephan Diel, Kenta Kurihara, Rio Shimizu, Eiji Murase
    INTERNATIONAL JOURNAL OF ENGINE RESEARCH, 22(8) 2560-2574, Jul, 2021  Peer-reviewedCorresponding author
    A new predictive combustion model for a one-dimensional computational fluid dynamics tool in the multibody dynamics processes of gasoline engines was developed and validated. The model consists of (1) a turbulent burning velocity model featuring a flame radius-based transitional function, steady burning velocity that considers local quenching using the Karlovitz number and laminarization by turbulent Reynolds number, as well as turbulent flame thickness and its quenching model near the liner wall, and (2) a knock model featuring auto-ignition by the Livengood-Wu integration and ignition delay time obtained using a full-kinetic model. The proposed model and previous models were verified under a wide range of operating conditions using engines with widely different specifications. Good agreement was only obtained for combustion characteristics by the proposed model without requiring individual calibration of model constants. The model was also evaluated for utilization after prototyping. Improved accuracy, especially of ignition timing, was obtained after further calibration using a small amount of engine data. It was confirmed that the proposed model is highly accurate at the early stage of the engine development process, and is also applicable for engine calibration models that require higher accuracy.
  • YABUSHITA Hirotaka, NAGAOKA Makoto, YOSHIOKA Masaya, MORI Yuichi
    Transactions of the JSME (in Japanese), 87(895) 20-00358-20-00358, Mar, 2021  Peer-reviewed
    A novel total hydrocarbon (THC) emission concentration estimation model is proposed for reduction of engine development cost as well as simplification of measurement system. The model is based on machine learning algorithm including the least absolute shrinkage and selection operator (LASSO) regression and bagging techniques. Major features of the proposal model are taking the absorbance spectrum of Fourier transform infrared (FTIR) spectrometer as input and incorporating not only spectra of the engine exhaust gas but also those of individual hydrocarbon and inorganic gas components as training data set. This method was validated on the exhaust gas before the catalyst of a gasoline engine. The results show an error of less than 5% in both steady and transient operating conditions, outperforming the 20 % error of conventional regression model using only the reference hydrocarbon concentrations. We also evaluate the contribution to performance improvements in THC estimation of employing FTIR spectrum and incorporating spectrum information of gas components, respectively.
  • MORIYASU Ryuta, UEDA Matsuei, IKEDA Taro, NAGAOKA Makoto, JIMBO Tomohiko, MATSUNAGA Akio, NAKAMURA Toshihiro
    Transactions of the Society of Instrument and Control Engineers, 55(3) 172-180, Mar, 2019  Peer-reviewedCorresponding author
    This paper deals with the machine learning based controller design for realizing nonlinear model predictive control (MPC) with low computational cost, and the application for a diesel engine air path system is shown. Since the solution of an optimal control problem considered in MPC depends on several variables at each time, the relationship between the variables and the solution, that is, the control law could be approximated by a neural network. In the case of high dimensional inputs, however, some efficient sampling methods for the approximation are needed to avoid a combinatorial explosion. To reduce the sampling dimension, we newly applied an efficient sampling method which is combined with a design of experiment. Using the method, the neural network structured optimal controller that operates the valves in the air path system was designed, and its tracking capability to the reference value was demonstrated in the simulation. The computational time of the controller was approximately 0.022ms at each cycle, indicating that fast computation of MPC was achieved.
  • Moriyasu Ryuta, Ueda Matsuei, Nagaoka Makoto, Ikeda Taro, Nishikawa Kazuaki, Nojiri Sayaka, Jimbo Tomohiko, Matsunaga Akio, Nakamura Toshihiro
    Transactions of Society of Automotive Engineers of Japan, 49(6) 1162-1166, Nov, 2018  Peer-reviewedCorresponding author
    This paper considers machine learning based virtual design process of engine control system, and the demonstration in a diesel engine air path control is shown. This process contains two steps of machine learning. In the first step, a control-oriented forward model that predicts the transient behavior of the engine is learned from detailed engine model by using recurrent neural network (RNN). In the second step, an inverse model that determines the optimal control inputs to follow the references is learned from the numerical computation results of the offline model predictive control (MPC). The forward and inverse models could be used as a state observer and a controller, respectively, in a control system. An experiment of a diesel air path control system designed by the process was conducted using rapid control prototyping (RCP), and its following capability to the reference was demonstrated.
  • R. Moriyasu, M. Ueda, T. Ikeda, M. Nagaoka, T. Jimbo, A. Matsunaga, T. Nakamura
    IFAC-PapersOnLine, 51(31) 542-548, Sep, 2018  Peer-reviewedCorresponding author
    This study investigated a control design method based on machine learning to achieve non-linear model predictive control (MPC) at a low computational load. In addition, we describe examples of the application of this method to a diesel engine air path system. The solution to the optimal control problem determined at each point in time by MPC depends on several parameters at that time. Thus, if the relationship between the solution and the parameters could be approximated in advance using machine learning, solving this problem online would become unnecessary, and the control computation time could be reduced. We designed a controller that operates the valves of the air path system using this method and used a simulation to verify that this resulted in a favorable tracking performance of the target values. The computation time of the approximated MPC controller was 0.022 ms.
  • Masuda Ryo, Sayama Shogo, Fuyuto Takayuki, Nagaoka Makoto, Noguchi Yasushi, Sugiura Akimitsu
    Transactions of Society of Automotive Engineers of Japan, 49(4) 732-738, Jul, 2018  Peer-reviewedCorresponding author
    The submodels of spark channel short circuit and blow-out, which were described in the 1st report, were implemented into a spark ignition model. The short circuit, whose major factor is electric potential differences between two arbitrary locations, occurs in the early phase of the discharge period and the blow-out, whose major factor is discharge current, occurs in the latter phase of the discharge period respectively. The behavior of spark channel, current and voltage of the secondary circuit and the ignition limit by increase of EGR rate agreed with a measurement data of a spark ignition process in a chamber.
  • 微粒化 = Atomization : journal of the ILASS-Japan, 27(91) 43-49, Jul, 2018  Peer-reviewedCorresponding author
  • MORIYASU Ryuta, MATSUMORI Tadayoshi, NAGAOKA Makoto
    Transactions of the JSME (in Japanese), 83(854) 17-00144-17-00144, Sep, 2017  Peer-reviewedCorresponding author
    A topology optimization method is proposed for the design of shallow-flow channels based on quasi-three-dimensional flow models of laminar and turbulent flows. The models for laminar flow and turbulent flow are derived from the Navier-Stokes equations and the Reynolds-Averaged Navier-Stokes (RANS) equations, respectively, by integrating along the direction of channel thickness. The thickness is employed as the design variable in the topology optimization. The design variables are updated using a time-dependent diffusion equation with a design sensitivity which is calculated by a discrete adjoint approach. Numerical examples for minimizing dissipation energy or variance of flow velocity magnitude using the topology optimization demonstrates that the proposed method is capable of finding optimal solutions that satisfy the KKT conditions. In the former example, the design domain was clearly divided into domains where the thickness was either near the upper limit or near the lower limit. However, in the latter example, the thickness was at an intermediate level in almost the whole the design domain. The distribution of the thickness varied depending on the Reynolds number in both examples.
  • NAGATA Mitsuhiro, NAGAOKA Makoto
    Transactions of the JSME (in Japanese), 82(838) 16-00073-16-00073, Jun, 2016  Peer-reviewedCorresponding author
    A direct numerical simulation of compressible turbulent channel flow at low Mach number(Ma = 0.3) subject to strong temperature gradient is conducted. The wall temperature ratios which are defined by a temperature on the upper wall divided by that on the lower one, are set to 2 and 3. It is shown that the flow of high temperature wall side is laminarized as increasing the wall temperature ratio. Moreover, compressible and dilatational motions are induced on the low and high temperature wall sides, respectively. The identity of friction coefficient originally derived by Fukagata et al.(2002) (so-called ‘FIK identity’) is extended to those of friction coefficient and Nusselt number on the compressible channel flow. It is revealed that the viscous variation component attains to 16 % of the friction coefficient on the case of high temperature ratio. Furthermore, the pressure work component occupies about 20 % of the overall Nusselt number for the all cases. The pressure work has a redistributive effect between mean kinetic energy and internal energy. The energy transfer from mean kinetic energy to internal energy appears on the low temperature side whereas that toward the opposite direction on the high temperature side at the present low Mach number flow.
  • Naoki Baba, Hiroaki Yoshida, Makoto Nagaoka, Chikaaki Okuda, Shigehiro Kawauchi
    Journal of Power Sources, 252 214-228, Apr, 2014  Peer-reviewedCorresponding author
    To understand the thermal behavior of lithium-ion secondary batteries, distributed information related to local heat generation across the entire electrode plane, which is caused by the electrochemical reaction that results from lithium-ion intercalation or deintercalation, is required. To accomplish this, we first developed an enhanced single particle (ESP) model for lithium-ion batteries that provides a cost effective, timely, and accurate method for estimating the local heat generation rates without excessive computation costs. This model accounts for all the physical processes, including the solution phase limitation. Next, a two-way electrochemical-thermal coupled simulation method was established. In this method, the three dimensional (3D) thermal solver is coupled with the quasi-3D porous electrode solver that is applied to the unrolled plane of spirally wound electrodes, which allows both thermal and electrochemical behaviors to be reproduced simultaneously at every computational time-step. The quasi-3D porous electrode solver implements the ESP model. This two-way coupled simulation method was applied to a thermal behavior analysis of 18650-type lithium-ion cells where it was found that temperature estimates of the electrode interior and on the cell can wall obtained via the ESP model were in good agreement with actual experimental measurements. © 2013 Elsevier B.V. All rights reserved.
  • Hiroaki Yoshida, Makoto Nagaoka
    JOURNAL OF COMPUTATIONAL PHYSICS, 257 884-900, Jan, 2014  Peer-reviewedCorresponding author
    A lattice Boltzmann method for simulating convection and diffusion using a curvilinear grid system is presented. The proposed method does not require an interpolation or coarse-graining procedure, and thus maintains the algorithmic simplicity of the original lattice Boltzmann scheme. The lattice Boltzmann scheme is based on uniformly distributed lattice points in a transformed coordinate system, and the apparent anisotropy of diffusion that arises due to the coordinate transformation is properly handled using the multiple-relaxation-time collision operator. An asymptotic analysis of the lattice Boltzmann equation shows that the proposed method appropriately reproduces the transformed convection-diffusion equation. Several specific problems are numerically analyzed in order to validate the proposed method, including an axially symmetric (two-dimensional) problem in which the diffusion flux at an oblate hemispheroid is simulated using a body-fitted orthogonal curvilinear grid system. (C) 2013 Elsevier Inc. All rights reserved.
  • MASUDA Ryo, KAWAMURA Kiyomi, NAGAOKA Makoto, MIZOBUCHI Takeshi, SUZUKI Hisao
    微粒化 = Atomization : journal of the ILASS-Japan, 22(76) 130-136, Nov, 2013  Peer-reviewedCorresponding author
  • KAWAMURA Kiyomi, MASUDA Ryo, IDOTA Yoshinori, NAGAOKA Makoto
    TRANSACTIONS OF THE JAPAN SOCIETY OF MECHANICAL ENGINEERS Series B, 78(793) 1598-1608, Sep, 2012  Peer-reviewedCorresponding author
    A poppet type nozzle which promotes atomization by utilizing cavitation phenomena was designed to clarify the atomization mechanism. The configuration of the nozzle which has two orifices and an intermediate cavity between the orifices was decided by numerical analysis on the cavitating flow in the nozzle. In this report, atomization mechanism and spray characteristics of the nozzle are described. The following results are obtained: (1) The net-like fuel liquid sheet is observed closed to orifice. The net-like structure with cavitation bubbles is collapsed to droplets via ligaments. (2) The present nozzle has the features of generating smaller drop size and lower spray tip penetration than the conventional poppet type nozzle.
  • MASUDA Ryo, KAWAMURA Kiyomi, NAGAOKA Makoto
    TRANSACTIONS OF THE JAPAN SOCIETY OF MECHANICAL ENGINEERS Series B, 78(793) 1584-1597, Sep, 2012  Peer-reviewedCorresponding author
    A liquid fuel injection nozzle with a novel fuel atomization strategy which actively utilizes the effect of collapse of cavitation bubbles in a liquid jet is proposed. The nozzle comprises two orifices and an intermediate cavity between the orifices. The role of the upstream first orifice is to generate the cavitation bubbles at both inlet and exit of the orifice. The intermediate cavity is placed to hold the cavitation bubble within it. The role of the downstream second orifice is to mix the cavitation bubbles and liquid fuel, and to inject the mixture into atmosphere. Each size and relative positions of the first orifice, the intermediate cavity and the second orifice were studied using multi-phase computational fluid dynamics. Finally, a poppet type nozzle whose void fraction at the nozzle exit is 0.5 was designed.
  • H. Kosaka, Y. Nomura, M. Nagaoka, M. Inagaki, M. Kubota
    International Journal of Engine Research, 12(4) 393-401, Jun, 2011  Peer-reviewedCorresponding author
    A novel combustion model for large-eddy simulation (LES) for gasoline engines has been developed. Unlike conventional models based on Reynolds-averaged Navier-Stokes (RANS) models, the new model takes a unique approach; it is described by the fractal characteristics of flame front and a universal expression for the subgrid scale (SGS) flame speed. The present fractal combustion model was applied to calculations of a spark ignition engine. Both the 0-10 per cent and 10-90 per cent combustion periods agree well with the experimental data. Because the modelling of the SGS turbulent speed is based on fractal analysis with experimental observations, the SGS combustion model is able to apply a wide range of engine operating conditions. The present model was applied to a multi-cycle simulation of a singlecylinder engine. The fluctuations at the instant when the heat release rate peaked were compared with data that was obtained experimentally. The calculated magnitude of the fluctuations was found to be close to the experimental values. It is thought that the flow variation generated during the intake stroke significantly influences the cyclic variations. Copyright © 2011 by institution of Mechanical Engineers.
  • 山本征治, 永岡真, 植田玲子, 脇坂佳史, 野田進
    日本燃焼学会誌, 53(165), 2011  Peer-reviewedCorresponding author
  • Hiroaki Yoshida, Makoto Nagaoka
    JOURNAL OF COMPUTATIONAL PHYSICS, 229(20) 7774-7795, Oct, 2010  Peer-reviewedCorresponding author
    A lattice Boltzmann model with a multiple-relaxation-time (MRT) collision operator for the convection-diffusion equation is presented. The model uses seven discrete velocities in three dimensions (D3Q7 model). The off-diagonal components of the relaxation-time matrix, which originate from the rotation of the principal axes, enable us to take into account full anisotropy of diffusion. An asymptotic analysis of the model equation with boundary rules for the Dirichlet and Neumann-type (specified flux) conditions is carried out to show that the model is first- and second-order accurate in time and space, respectively. The results of the analysis are verified by several numerical examples. It is also shown numerically that the error of the MRT model is less sensitive to the variation of the relaxation-time coefficients than that of the classical BGK model. In addition, an alternative treatment for the Neumann-type boundary condition that improves the accuracy on a curved boundary is presented along with a numerical example of a spherical boundary. (C) 2010 Elsevier Inc. All rights reserved.
  • 41(4) 827-832, Sep, 2010  Peer-reviewedCorresponding author
  • M. Inagaki, M. Nagaoka, N. Horinouchi, K. Suga
    International Journal of Engine Research, 11(3) 229-241, Jun 1, 2010  Peer-reviewed
    Large eddy simulation (LES) using a mixed-time-scale (MTS) subgrid-scale (SGS) model is applied to the intake flows in simplified internal combustion engine geometry. A modified colocated grid system is employed to obtain results as precise as possible and to perform calculations in a stable way with a central difference scheme for convective terms. The results are compared with corresponding experimental data and the Reynolds averaged Navier-Stokes (RANS) equation model results obtained using the low-Reynolds-number linear k- In addition, it is made clear that when the QUICK scheme is used in LES for the convective terms instead of the central difference scheme, the result obtained deteriorates owing to the numerical viscosity. The importance of the discretization method in practical LES is also confirmed. © IMechE 2010.
  • S. Yamamoto, M. Nagaoka, R. Ueda, Y. Wakisaka, S. Noda
    International Journal of Engine Research, 11(1) 17-27, Feb, 2010  Peer-reviewedCorresponding author
    A series calculation method from injector nozzle internal flow to in-cylinder combustion in diesel engines has been developed. In this study, the calculation results of in-cylinder pressure, heat release rate, and NO concentration are compared with experimental results for conventional combustion conditions as well as an advanced combustion condition, which combines a high exhaust gas recirculation (EGR) rate, high boost pressure, and high pressure injection. The calculated pressure and heat release rate histories are in reasonable agreement with those of experiments, but the NO concentration is underestimated especially on the higher EGR rate condition. Therefore, an investigation to improve the quantitative accuracy of NO concentration is performed. The probability density function (PDF) method is one way to estimate the accurate mean reaction rate containing the reaction rate fluctuation by turbulence. In this study, the PDF method only applied for NO calculation is developed to balance the calculation accuracy and the computational cost. The result shows that the reaction rate fluctuation, on NO formation, has a small effect on conventional combustion conditions, but increases with decreasing combustion temperature. © IMechE 2010.
  • INAGAKI Masahide, NAGAOKA Makoto, HORINOUCHI Nariaki, SUGA Kazuhiko
    Transactions of the Japan Society of Mechanical Engineers. B, 73(725) 58-66, Jan 25, 2007  Peer-reviewed
    Large eddy simulation (LES) using a mixed-time-scale (MTS) SGS model is applied to the intake flows in simplified internal combustion engine geometry. A modified colocated grid system is employed so as to perform calculations stably with a central difference scheme for convection terms. The results are compared with corresponding experimental data and the computational results obtained by using the low-Reynolds-number linear κ-ε and cubic nonlinear κ-ε-A_2 turbulence models. The LES results show the best agreement with experimental data in three computational cases, not only in the mean velocity profiles but also in the profiles of turbulent energy. These results suggest that LES using the MTS SGS model is an effective method for accurately predicting the performance of a combustion engine involving turbulent diffusion of spray and combustion flame. In addition, the results are compared with those obtained by a quasi-direct simulation employing the QUICK scheme for the convection term and no turbulence model. It is clarified that the numerical viscosity from the QUICK scheme has a great influence on the computational results and decreases the prediction accuracy.
  • K. Suga, M. Nagaoka, N. Horinouchi
    Journal of Heat Transfer, 125(1) 200-203, Feb, 2003  Peer-reviewed
    A higher order version of the generalized gradient diffusion hypothesis (HOGGDH) for turbulent heat flux is applied to predict heat transfer in a square-sectioned U-bend duct. The flow field turbulence models coupled with are a cubic nonlinear eddy viscosity model and a full second moment closure. Both of them are low Reynolds number turbulence models. The benefits of using the HOGGDH heat flux model are presented through the comparison with the standard GGDH.
  • NAGAOKA Makoto, SUGA Kazuhiko, HORINOUCHI Nariaki
    TRANSACTIONS OF THE JAPAN SOCIETY OF MECHANICAL ENGINEERS Series B, 68(669) 1586-1592, May, 2002  Peer-reviewedLead authorCorresponding author
    This paper presents discussions on predicting turbulent flow in a simplified internal combustion engine geometry consisting of an intake port, valve and cylinder. The low-Reynolds-number linear k-ε, cubic nonlinear k-ε and k-ε-A2 turbulence models are applied to the flow field computation and evaluated against experimental data. The results suggest that the cubic nonlinear eddy viscosity modelling of turbulence is very successful for predicting mass flow rates and in-cylinder velocity distributions.
  • SUGA Kazuhiko, NAGAOKA Makoto, HORINOUCHI Nariaki
    Transactions of the Japan Society of Mechanical Engineers. B, 68(666) 495-503, Feb 25, 2002  Peer-reviewed
    This paper presents discussions on predicting turbulence and heat transfer in square sectioned 180° U-bend ducts whose curvature ratios are Rc/D=0.65 and 3.357. The low-Reynolds-number cubic nonlinear κ-ε and κ-ε-A_2 turbulence models are applied for the flow field computation. In the thermal field computations, the higher order gradient diffusion heat flux model is applied in comparison with the standard models. The results suggest that the cubic nonlinear eddy viscosity and the higher order gradient diffusion heat flux models are very successful for predicting such three dimensional turbulence and heat transfer fields while some more modification in the models may be needed to capture heat transfer in the strongly curved case.
  • Makoto Nagaoka, Kiyomi Kawamura
    SAE Transactions,JOURNAL OF ENGINES, 110(3) 1326-1335, Jan, 2002  Peer-reviewedLead authorCorresponding author
  • K. Suga, M. Nagaoka, N. Horinouchi, K. Abe, Y. Kondo
    International Journal of Heat and Fluid Flow, 22(3) 259-271, Jun, 2001  Peer-reviewed
    The three-equation cubic k-ε-A2 model proposed by Craft et al. (Int. J. Heat Fluid Flow 18 (1997) 15-28) is evaluated in three-dimensional (3-D) turbulent flows pertinent to engineering applications, especially in the automobile industry. For the computations of complex industrial flows, a numerical scheme has been developed using the cell vertex unstructured grid method. This scheme treats a mixture of tetrahedral, pyramidal, prismatic and hexahedral computational cells with high accuracy. The industrial flows chosen are internal combustion (IC) engine port-cylinder flows and flows around aerodynamic bluff bodies. The model performance in U-bend duct flows and a flow around a surface-mounted cubical obstacle is also examined. These fundamental flows include essential features of the industrial flows presently focused on. The model performs generally satisfactorily. However, the performance in a 3-D separating wake flow behind a bluff body suggests that the model needs further improvements. © 2001 Elsevier Science Inc. All rights reserved.
  • NAGAOKA Makoto
    66(647) 1885-1891, Jul 25, 2000  Peer-reviewedLead authorCorresponding author
    A theoretical model for a deforming single droplet with aerodynamic force is presented. The droplet shape is assumed to be spheroidal. A nonlinear oscillating equation is derived from energy conservation in a droplet. The present (OSD) model is compared with other droplet deformation models, Taylor Analogy Breakup (TAB), Droplet Deformation Breakup (DDB) and Synthesized Spheroidal Particle (SSP) models. It is found that the linearized OSD model is reduced to the TAB model. Numerical studies show that only present model and TAB model are compatible with Lamb's fundamental oscillating mode for a liquid droplet. On a droplet deformation at the initial stage of aerodynamic breakup, the calculation by the OSD model and TAB model agree well with the measurement.
  • Hiroshi Miyagawa, Makoto Nagaoka, Kazuhiro Akihama, Taketoshi Fujikawa
    SAE Transactions,JOURNAL OF FUELS AND LUBRICANTS, 108(4) 2062-2071, Jan, 2000  Peer-reviewed
  • Hiroshi Miyagawa, Makoto Nagaoka, Katsuyuki Ohsawa, Toshio Yamada
    JSAE Review, 19(4) 299-304, Oct 1, 1998  Peer-reviewed
    A multi-component droplet vaporization model for multi-dimensional calculation was studied in order to simulate the behavior of fuel droplets and mixture in a port injection gasoline engine. Calculation results for a single droplet show that the evaporation process of gasoline, which consists of more than 100 components, can be simulated using a model fuel composed of at least three representative species. Subsequently, three-dimensional calculations of flow and fuel spray for an actual engine intake port configuration were performed, using the developed multi-component vaporization model. It was shown that about 70% of the high-volatility component in injected gasoline flows into the cylinder, while about 70% of the low-volatility component stays in the intake port in the injected cycle. © 1998 Society of Automotive Engineers of Japan, Inc. and Elsevier Science B.V. All rights reserved.
  • Makoto Nagaoka, Katsuyuki Ohsawa, Brent Crary, Toshio Yamada, Shigeki Sugiura, Nobuo Imatake
    SAE Transactions,JOURNAL OF ENGINES, 106(3) 1369-1376, Jan, 1998  Peer-reviewedLead author
  • 永岡真, CRARY B, 大沢克幸, 山田敏生, 杉浦繁貴, 今竹信夫
    日本機械学会論文集 B編, 63(611), 1997  Peer-reviewedLead author
  • NAGAOKA Makoto, LOHNER Rainald
    TRANSACTIONS OF THE JAPAN SOCIETY OF MECHANICAL ENGINEERS Series B, 63(605) 153-160, Jan, 1997  Peer-reviewedLead authorCorresponding author
    Preconditioned iterative methods in an implicit unstructured grid solver of compressible flows are compared. The three-dimensional Euler and Navier-Stokes equations are discretized by the edge-based finite volume method. The artificial dissipation scheme for the inviscid numerical fluxes is developed for the unstructured grid from Jameson's scheme for the structured grid. The following methods are compared to the supersonic and transonic inviscid flows and subsonic viscous laminar and turbulent flows : the central difference methods with the scalar and matrix dissipation scheme and Roe's upwind method with the MUSCL sheme for the inviscid fluxes ; the scalar and matrix dissipation method for the left-hand side operator in the implicit scheme ; the ILU, DILU and SGS factorizations for the preconditioner ; and GMRES and Bi-CGSTAB methods for the linear system solver. The results show that a combination of the matrix dissipation in the inviscid fluxes and the scalar dissipation in the left-hand side works better than the other shemes. In the preconditioned iterative methods, a combination of DILU or SGS and the Bi-CGSTAB is recommended for the aspect of required memory.
  • NAGAOKA Makoto, NOMURA Naomi
    61(587) 2744-2750, Jul 25, 1995  Peer-reviewedLead author
    A calculation method for steady compressible flows using a solution adaptive unstructured mesh is proposed and it's applicability is demonstrated for two-dimensional and axisymmetric engine intake flows. The mesh is generated by the rule-based Delaunay triangulation method. The governing equations are discretized by the finite volume method. The second order accurate Roe's upwind scheme is used for the inviscid flux calculation. The viscous fluxes are calculated by the transformation to the local general coordinates. The utility of the adaptive unstructured mesh method is shown in the 2D engine intake flow calculations. The present method using a low Reynolds number k-c model is valid for the velocity profile in the turbulent boundary layer. The calculated discharge coefficients with valve lifts in an axisymmetric engine agree well with the measurements.
  • Makoto Nagaoka, Hiromitsu Kawazoe, Naomi Nomura
    SAE Transactions,JOURNAL OF ENGINES, 103(3) 878-896, Jan, 1995  Peer-reviewedLead author
  • Makoto Nagaoka, Nariaki Horinouchi
    日本数値流体力学会, CFD Journal, Vol.2(No.2) 169-180, Sep, 1993  Peer-reviewedLead author
    圧縮性流れの解法において、非構造格子上で有限体積法による離散化と2次精度の空間差分スキームを用いる際、時間解法には前処理付きBi-CGSTAB法が当時一般的なRunge-Kutta法やガウスザイデル法に比べて定常解への収束が早くかつ安定であることを示した。
  • Osamu Takata, Koukichi Nakanishi, Nariaki Horinouchi, Makoto Nagaoka
    (92号) 41-48, Nov, 1992  Peer-reviewed
  • KAWAZOE Hiromitsu, NAGAOKA Makoto, MIZUTA Junichi, OHSAWA Katsuyuki
    58(547) p965-971, Mar, 1992  Peer-reviewed
    Numerical prediction of a steady intake airflow was performed with an unstructured grid. Changing surface rougheness and configuration of an intakeport, airflow rate and velocity vectors were calculated. Comparing the calculated results with the experimental ones by LDV measurements, the accuracy was evaluated. It was found that airflow rate can be predicted with the accuracy of 5 percent, and the predicted velocities, especially in the intakeport, agreed with the measured ones. The predicted airjets at the valve annular passage were slower than the experimental ones. The application of the standard k-ε turbulence model to the airjet at the valve passage was found to be an essential factor for the discrepancy.
  • Y. Tsujikawa, M. Nagaoka
    Journal of Engineering for Gas Turbines and Power, 113(1) 100-105, Jan, 1992  Peer-reviewed
    This paper is devoted to the analyses and optimization of simple and sophisticated cycles, particularly for various gas turbine engines and aero-engines (including the scramjet engine) to achieve maximum performance. The optimization of such criteria as thermal efficiency, specific output, and total performance for gas turbine engines, and overall efficiency, nondimensional thrust, and specific impulse for aero-engines has been performed by the optimization procedure with the multiplier method. Comparison of results with analytical solutions establishes the validity of the optimization procedure.
  • NAGAOKA Makoto, OHSAWA Katsuyuki
    TRANSACTIONS OF THE JAPAN SOCIETY OF MECHANICAL ENGINEERS Series B, 57(540) 2678-2683, Aug, 1991  Peer-reviewedLead author
    The unstructured upwind method is proposed to solve compressible flow with complicated boundaries or an arbitrarily shaped mesh. The flux difference splitting (FDS) scheme and flux vector splitting (FVS) scheme, which are extended to the unstructured grid system, are compared. The Euler equation is discretized by the control volume method. All conserved variables are stored at the cell center. The MUSCL approach for the cell-centered unstructured grid is proposed to obtain higher-order accuracy. A three-stage Runge-Kutta method is used for the time stepping scheme. The present calculation methods are applied to two-dimensional front-facing step flows. The results demonstrate that the present MUSCL approach is effective to obtain high resolution, and the triangular mesh solver is comparable to the uniform orthogonal mesh solver. FDS gives slightly better resolution than FVS.
  • Kiyomi Nakakita, Makoto Nagaoka, Taketoshi Fujikawa, Katsuyuki Ohsawa, Shigeki Yamaguchi
    SAE Transactions,JOURNAL OF ENGINES,Part 2, 99(3) 2132-2144, Jan, 1991  Peer-reviewed
  • Ichiro Sakata, Kazuyoshi Ishisaka, Hiromichi Yanagihara, Hiroshi Sami, Makoto Nagaoka, Katsuyuki Ohsawa
    SAE Transactions,JOURNAL OF ENGINES, 99(3) 1488-1494, Jan, 1991  Peer-reviewed
  • KAWAZOE Hiromitsu, NAGAOKA Makoto, SUZUKI Sigeo, OHSAWA Katsuyuki
    57(533) p385-390, Jan, 1991  Peer-reviewed
    Fuel mixture formation processes were predicted using gasoline droplet data obtained by LDA measurements. In the intake period, air and fuel inlet boundary conditions were given at the valve annular passage. The last calculated result of the intake was used as the initial condition of the compression period. The effects of combustion chamber configuration on the mixture homogeneity were investigated in a lean-burn engine. Four types of piston cavities were investigated: pancake, Heron, off-cylinder-axis Heron type, and complicated configuration chamber. The complicated one improves homogeneity since the flow produces two inclined vortices in the cylinder, which play a great role in promoting convective mixing.
  • 永岡真, 大沢克幸, 阪田一郎
    自動車技術会論文集, No.44 3-8, Mar, 1990  Peer-reviewedLead author
    新開発の小型直噴ディーゼルエンジンにおける燃焼室内の流れと燃料噴霧衝突挙動を3次元シミュレーションで解析し、その燃焼室形状が乱れを強め、混合気形成を促進していることがわかった。
  • NAKAKITA Kiyomi, NAGAOKA Makoto, OHSAWA Katsuyuki
    56(521) p221-226, Jan, 1990  Peer-reviewed
    Diesel soot formation processes were simulated, using several models previously proposed for soot formation and oxidation. The models were assessed and improved by comparing the results with the experiments. Consequently, the model which represented qualitatively the total mass history, the formation position and the spacial distribution of soot was established. The outline and characteristics of the model are as follows. (1) The soot formation model consists of both the Tesner and Farmer models. The soot oxidation model consists of both the Nagle and Magnussen models. The model giving a smaller rate was selected for both models (formation and oxidation). (2) Some of the coefficients in the Tesner model were corrected depending on the in-cylindr volume change. (3) The Farmer model was used only to cover the deficiency of the Tesner model which must be improved due to tempareture dependence.
  • TSUJIKAWA Yoshiharu, SAWADA Teruo, NAGAOKA Makoto, TSUKAMOTO Yujiro
    53(493) p2915-2920, Sep, 1987  Peer-reviewed
    Many works on gas turbine engines contain much information concerning the performance of various engine schemes both in simpler and in more complicated forms. Such information enables the design performance to be estimated, but there appears to have been very little works on the optimum performance under non-design conditions. In this paper, we attempt to include a comparison of part load performance characteristics and to optimize them. The multiplier method is adopted as an optimization procedure. As a result of this work, it is possible to indicate the optimization technique by which the comparable part load performance of different gas turbine schemes may be assessed.
  • TSUJIKAWA Yoshiharu, SAWADA Teruo, NAGAOKA Makoto, TSUKAMOTO Yujiro
    53(491) p2219-2225, Jul, 1987  Peer-reviewed
    This paper is devoted to the analyses and optimization of simple and sophisticated cycles, particularly for various gas turbine engines and aero-engines, to achieve maximum efficiency or maximum specific output. The optimization of such criteria as thermal efficiency, specific output and total performance for gas turbine engines, and overall efficiency and non-dimensional thrust for aero-engines have been performed by the multiplier method. The comparisons of results with analytical solutions establishes the validity of the optimization technique.

Misc.

 20
  • Hirotaka Iseki, Makoto Nagaoka, Shuntaro Yokoi, Naoto Horibe, Hiroshi Kawanabe
    SAE Technical Papers 2021-01-0603, (2021), Apr, 2021  Peer-reviewed
    For the measurements of flow rate, pressure and/or temperature in an engine exhaust pipe, probes are often inserted into the exhaust pipe depending on the application. These measurement probes differ a lot in terms of their size and shape. The flow around the probes become further complicated due to the pulsation of engine exhaust flow. In this study, computational fluid dynamics (CFD) simulations were carried out and a zero-dimensional (0D) model was constructed to analyze the flow field around the probe and flow rate of a pulsating flow. The simulations and the measurements of the flow rate and pressure were performed on flows around a hexagonal prism inserted in a circular pipe which is intended to be a differential pressure flow meter. The velocity field was also measured using the particle image velocimetry (PIV) technique. The CFD simulation results were validated with the experiments for both steady and pulsating flows. In the 0D model for pulsating flow, the flow acceleration as well as pipe friction and prism drag losses were taken into account. The flow rates calculated using the model agreed well with the CFD simulation results. The relationship between the flow rate and the pressure was analyzed using the CFD and the 0D model. In the low flow rate and low pressure difference period, the relationship between the flow rate and the square root of pressure difference deviated from linear and exhibited hysteresis due to the flow acceleration. The cycle-averaged flow rates calculated using the 0D model were closer to those by the CFD simulations than those of a conventional steady flow correlation.
  • Ryo Masuda, Shogo Sayama, Takayuki Fuyuto, Makoto Nagaoka, Akimitsu Sugiura, Yasushi Noguchi
    SAE Technical Papers 2018-01-1727, Sep, 2018  Peer-reviewed
    This report describes the implementation of the spark channel short circuit and blow-out submodels, which were described in the previous report, into a spark ignition model. The spark channel which is modeled by a particle series is elongated by moving individual spark particles along local gas flows. The equation of the spark channel resistance developed by Kim et al. is modified in order to describe the behavior of the current and the voltage in high flow velocity conditions and implemented into the electrical circuit model of the electrical inductive system of the spark plug. Input parameters of the circuit model are the following: initial discharge energy, inductance, internal resistance and capacitance of the spark plug, and the spark channel length obtained by the spark channel model. The instantaneous discharge current and the voltage are obtained as outputs of the circuit model. When two arbitrary spark particles of the spark channel get close, the short circuit occurs if the electric potential differences between the two locations exceed a certain threshold voltage, which is raised with increasing distance between the two particles and decreasing discharge current. When the current falls below a lower limit current for maintenance of discharge, the spark blow-out occurs. A new spark channel is formed if the secondary circuit has the remaining energy which can break the electrical insulation between electrodes. Each line element of the spark channel particles heats and ignites the surrounding mixture gas. The turbulent flame speed and extinction are considered in the flame kernel behavior. The behavior of the spark channel, the current and voltage of the secondary circuit, and the ignition limit due to in-creases in the EGR rate were consistent with data measured from the spark ignition process in a combustion chamber.
  • Mitsuhiro Nagata, Makoto Nagaoka
    Proc. 10th International Symposium on Turbulence and Shear Flow Phenomena (TSFP10), Oct, 2017  Peer-reviewed
  • Ryo Masuda, Kiyomi Kawamura, Makoto Nagaoka
    R&D Review of Toyota CRDL, Vol.45(No.3) 73-75, Sep, 2014  
  • Makoto Nagaoka, Reiko Ueda, Ryo Masuda, Eberhard von Berg, Reinhard Tatschl
    R&D Review of Toyota CRDL, Vol.42(No.2) 73-84, Jun, 2011  

Presentations

 1

Professional Memberships

 3

Industrial Property Rights

 13

研究テーマ

 2
  • 研究テーマ(英語)
    自動車の内燃機関・パワートレイン要素のサロゲートモデルの研究
    研究期間(開始)(英語)
    2015
  • 研究テーマ(英語)
    Analysis and modeling of fuel behavior for the utilization of carbon-neutral fuel and renewable energy
    研究期間(開始)(英語)
    2022