Railway Vehicle Aerodynamics
Theory and Physics
Overview
Professor, what makes the aerodynamic analysis of high-speed trains like the Shinkansen so difficult?
Railway vehicle aerodynamics has many unique challenges. These include pressure waves (micro-pressure waves) upon tunnel entry, crosswind stability, passing aerodynamic effects, and the reduction of running resistance.
In high-speed railways like the Shinkansen, aerodynamic drag accounts for over 80% of the total running resistance. For vehicles exceeding 300 km/h, aerodynamic design directly impacts power consumption.
Governing Equations
The aerodynamic drag of a train is proportional to the square of its speed.
Here, the first term is pressure drag (dependent on the nose/tail shape), and the second term is friction drag (on the vehicle sides, where $P$ is the perimeter and $L$ is the train length).
It seems like friction drag would contribute significantly for long trains.
For a 16-car formation like the N700S (total length approx. 400m), friction drag accounts for 60--70% of the total resistance. Therefore, it's not just the nose shape; reducing surface irregularities like gaps between cars, pantographs, and bogie fairings is also crucial.
Tunnel Micro-Pressure Wave
When a high-speed train enters a tunnel, a compression wave is generated and released like a shock wave at the opposite exit. This is the micro-pressure wave (tunnel boom).
The magnitude of the micro-pressure wave is proportional to the cube of the train speed. If the speed doubles, the micro-pressure wave becomes eight times larger, making nose shape optimization essential for higher speeds.
So that's why Shinkansen noses keep getting longer. The 500 Series is 15m, and the N700S is 10.7m, right?
Exactly. By making the rate of change of the nose cross-sectional area $dA/dx$ more gradual, the micro-pressure wave is reduced. Optimizing the cross-sectional area distribution using CFD is a modern design method.
Crosswind Stability
Overturning moment coefficient under crosswind:
Here, $V_{rel}$ is the resultant wind speed from the train speed and crosswind, and $M_y$ is the rolling moment about the rail surface.
The European standard (EN 14067-6) permits the use of CFD analysis to calculate the Characteristic Wind Curve (CWC) for crosswind conditions. Aerodynamic coefficients are obtained by calculating for yaw angles $\beta$ in the range of 0--90 degrees.
| Yaw Angle $\beta$ | Aerodynamic Characteristics |
|---|---|
| 0 degrees | Pure headwind. Only $C_D$ |
| 10--30 degrees | Side force and rolling moment increase sharply |
| 30--60 degrees | Maximum side force region. Maximum overturning risk |
| 90 degrees | Pure crosswind. The train behaves like a prism |
The Shinkansen's "Long Nose" is Designed for Tunnel Micro-Pressure Waves
When a Shinkansen enters a tunnel at high speed, a "boom!" sound (micro-pressure wave, also called a sonic boom) occurs near the exit. This problem became serious after the opening of the Sanyo Shinkansen in the 1970s, and CFD was used to identify the cause. The more abrupt the nose shape, the steeper the compression wave front, leading to larger pressure fluctuations at the tunnel exit. Therefore, by making the nose shape longer and with a gentler curve, the compression wave is "smoothed out," reducing noise. The 500 Series Shinkansen's 15m-long nose is the crystallization of this design philosophy.
Physical Meaning of Each Term
- Temporal Term $\partial(\rho\phi)/\partial t$: Imagine the moment you turn on a faucet. At first, water comes out in an unstable, splashing manner, but after a while, it becomes a steady flow, right? This "period of change" is described by the temporal term. The pulsation of blood flow from a heartbeat, or the flow fluctuation each time an engine valve opens and closes—all are unsteady phenomena. So what is steady-state analysis? It looks only at "after sufficient time has passed and the flow has settled down"—meaning setting this term to zero. This significantly reduces computational cost, so starting with a steady-state solution is a basic CFD strategy.
- Convection Term $\nabla \cdot (\rho \mathbf{u} \phi)$: What happens if you drop a leaf into a river? It gets carried downstream by the flow, right? This is "convection"—the effect where fluid motion transports things. Warm air from a heater reaching the far side of a room is also because the air, as a "carrier," transports heat via convection. Here's the interesting part—this term contains "velocity × velocity," making it nonlinear. That is, as the flow becomes faster, this term rapidly strengthens, making control difficult. This is the root cause of turbulence. A common misconception: "Convection and conduction are similar" → They are completely different! Convection is carried by flow, conduction is transmitted by molecules. There is an order of magnitude difference in efficiency.
- Diffusion Term $\nabla \cdot (\Gamma \nabla \phi)$: Have you ever put milk in coffee and left it? Even without stirring, after a while, it naturally mixes, right? That's molecular diffusion. Now a question—honey and water, which flows more easily? Obviously water, right? Honey has high viscosity ($\mu$), so it flows poorly. Higher viscosity strengthens the diffusion term, making the fluid move in a "thick" manner. In low Reynolds number flow (slow, viscous), diffusion is dominant. Conversely, in high Re number flow, convection overwhelms, and diffusion becomes a minor player.
- Pressure Term $-\nabla p$: When you push the plunger of a syringe, liquid shoots out forcefully from the needle tip, right? Why? Because the piston side is high pressure, and the needle tip is low pressure—this pressure difference provides the force that pushes the fluid. Dam discharge works on the same principle. On a weather map, where isobars are tightly packed? That's right, strong winds blow. "Flow occurs where there is a pressure difference"—this is the physical meaning of the pressure term in the Navier-Stokes equations. A point of confusion here: "Pressure" in CFD is often gauge pressure, not absolute pressure. If results become strange when switching to compressible analysis, it might be due to confusion between absolute and gauge pressure.
- Source Term $S_\phi$: Heated air rises—why? Because it becomes lighter (lower density) than its surroundings, so it's pushed up by buoyancy. This buoyancy is added to the equation as a source term. Other examples: chemical reaction heat from a gas stove flame, Lorentz force acting on molten metal in a factory's electromagnetic pump... These are all actions that "inject energy or force into the fluid from the outside," expressed by the source term. What happens if you forget the source term? In natural convection analysis, forgetting buoyancy means the fluid doesn't move at all—a physically impossible result where warm air doesn't rise in a room with the heater on in winter.
Assumptions and Applicability Limits
- Continuum Assumption: Valid for Knudsen number Kn < 0.01 (molecular mean free path ≪ characteristic length)
- Newtonian Fluid Assumption: Shear stress and strain rate have a linear relationship (non-Newtonian fluids require viscosity models)
- Incompressibility Assumption (for Ma < 0.3): Treat density as constant. For Mach numbers above 0.3, compressibility effects must be considered
- Boussinesq Approximation (Natural Convection): Density variation is considered only in the buoyancy term; constant density is used in other terms
- Non-applicable Cases: Rarefied gases (Kn > 0.1), supersonic/hypersonic flow (requires shock wave capturing), free surface flow (requires VOF/Level Set, etc.)
Dimensional Analysis and Unit Systems
| Variable | SI Unit | Notes / Conversion Memo |
|---|---|---|
| Velocity $u$ | m/s | When converting from volumetric flow rate for inlet conditions, pay attention to cross-sectional area units |
| Pressure $p$ | Pa | Distinguish between gauge and absolute pressure. Use absolute pressure for compressible analysis |
| Density $\rho$ | kg/m³ | Air: approx. 1.225 kg/m³ @20°C, Water: approx. 998 kg/m³ @20°C |
| Viscosity Coefficient $\mu$ | Pa·s | Be careful not to confuse with kinematic viscosity $\nu = \mu/\rho$ [m²/s] |
| Reynolds Number $Re$ | Dimensionless | $Re = \rho u L / \mu$. Criterion for laminar/turbulent transition |
| CFL Number | Dimensionless | $CFL = u \Delta t / \Delta x$. Directly related to time step stability |
Numerical Methods and Implementation
Mesh and Computational Domain
How large-scale is CFD for trains?
Trains are very long (over 400m total length), so the computational scale becomes large.
| Analysis Target | Typical Cell Count | Computational Domain Size |
|---|---|---|
| Lead car only | 20--50 million | 5 times car length |
| 3-car formation | 50--100 million | 3--5 times formation length |
| Full formation | 100--500 million | Formation length + wake region |
| Tunnel entry | 50--200 million | Full tunnel length + front/rear |
Calculating a full formation is often impractical, so an approach using a model of the lead + 1-2 intermediate cars + tail, and interpolating the friction drag of intermediate cars with empirical formulas, is common.
Turbulence Model
Which turbulence models are used in railway vehicle CFD?
Selected according to the application.
Application Recommended Model Reason Steady aerodynamic drag SST k-omega Prediction accuracy for separation and reattachment Crosswind stability SST k-omega / DDES Unsteadiness of large-scale separation Tunnel micro-pressure wave Compressible RANS Capturing pressure wave propagation Passing aerodynamic effects Unsteady RANS / LES Rapid pressure fluctuations In-cabin pressure fluctuations Unsteady RANS Passenger ear discomfort
Tunnel Entry Analysis Methods
How is CFD analysis for tunnel micro-pressure waves done?
Solve the process of a train entering a tunnel unsteadily using a compressible solver. There are two approaches.
1. Sliding Mesh Method
- The train physically moves. Most faithful reproduction
- High computational cost
- STAR-CCM+'s Overset Mesh is suitable
2. Moving Reference Frame Method
- In a train-fixed coordinate system, the tunnel approaches
- No mesh movement required, but inlet/outlet treatment is complex
Pressure wave propagation occurs at the speed of sound, so the time step is determined from the CFL condition:
$$ \Delta t < \frac{\Delta x}{c + V_{train}} $$
Here, $c \approx 340$ m/s is the speed of sound, and $V_{train}$ is the train speed.
A very small time step is needed.
That's right. For a cell size of 0.1m, $\Delta t < 0.0002$ seconds. Since tunnel passage takes several seconds, tens of thousands of time steps are needed.
Ground Effect and Flow Around Bogies
Since trains run close to the ground, ground effect is important.
- Moving Ground: Wall condition moving at the same speed as the train
- Ballast Track Bed: Set as a rough wall with roughness
- Bogie Fairing: Significant streamlining effect, can reduce resistance by 10--15%
- Inter-car Gap Cover: Reduces gap wind, cuts friction drag by 5--8%
The effect of bogie fairings is significant.
The N700S reduced aerodynamic drag by about 7% compared to the N700A by adopting full-skirt and bogie fairings. This is the result of shape optimization using CFD.
Coffee Break Trivia
Tunnel Micro-Pressure Wave—The "Shock Wave" Created by the Shinkansen
When a Shinkansen enters a tunnel, the air in front of the train is compressed, and a weak shock wave called a "micro-pressure wave" radiates from the exit. This creates the "boom!" explosion sound, which became a serious noise problem for residents along the lines in the 1990s. The countermeasure adopted was a significant lengthening of the nose shape—the 500 Series Shinkansen's 15m-long nose is for this reason. The method of numerically analyzing compression waves inside tunnels using CFD and optimizing the nose shape is one of the most important challenges in railway aerodynamic CFD. This analysis requires handling unsteady compressible fluids, and the computational cost is correspondingly high.
Upwind Scheme
First-order upwind: Large numerical diffusion but stable. Second-order upwind: Improved accuracy but risk of oscillations. Essential for high Reynolds number flows.
Central Differencing
Second-order accurate, but numerical oscillations occur for Pe number > 2. Suitable for low Reynolds number, diffusion-dominated flows.
TVD Scheme (MUSCL, QUICK, etc.)
Maintains high accuracy while suppressing numerical oscillations via limiter functions. Effective for capturing shock waves or steep gradients.
Finite Volume Method vs Finite Element Method
FVM: Naturally satisfies conservation laws. Mainstream in CFD. FEM: Advantageous for complex shapes and multi-physics. Mesh-free methods like SPH are also developing.
CFL Condition (Courant Number)
Explicit method: CFL ≤ 1 is the stability condition. Implicit method: Stable even for CFL > 1, but affects accuracy and iteration count. LES: CFL ≈ 1 is recommended. Physical meaning: Information should not travel more than one cell per time step.
Residual Monitoring
Convergence is judged when residuals for the continuity equation, momentum, and energy decrease by 3-4 orders of magnitude. The mass conservation residual is particularly important.
Relaxation Factor
Pressure: 0.2-0.3, Velocity: 0.5-0.7 are typical initial values. If diverging, lower the relaxation factor. After convergence, increase to accelerate.
Internal Iterations for Unsteady Calculations
Related Topics
この記事の評価ご回答ありがとうございます！参考に
なったもっと
詳しく誤りを
報告
Which turbulence models are used in railway vehicle CFD?
Selected according to the application.
| Application | Recommended Model | Reason |
|---|---|---|
| Steady aerodynamic drag | SST k-omega | Prediction accuracy for separation and reattachment |
| Crosswind stability | SST k-omega / DDES | Unsteadiness of large-scale separation |
| Tunnel micro-pressure wave | Compressible RANS | Capturing pressure wave propagation |
| Passing aerodynamic effects | Unsteady RANS / LES | Rapid pressure fluctuations |
| In-cabin pressure fluctuations | Unsteady RANS | Passenger ear discomfort |
How is CFD analysis for tunnel micro-pressure waves done?
Solve the process of a train entering a tunnel unsteadily using a compressible solver. There are two approaches.
1. Sliding Mesh Method
- The train physically moves. Most faithful reproduction
- High computational cost
- STAR-CCM+'s Overset Mesh is suitable
2. Moving Reference Frame Method
- In a train-fixed coordinate system, the tunnel approaches
- No mesh movement required, but inlet/outlet treatment is complex
Pressure wave propagation occurs at the speed of sound, so the time step is determined from the CFL condition:
Here, $c \approx 340$ m/s is the speed of sound, and $V_{train}$ is the train speed.
A very small time step is needed.
That's right. For a cell size of 0.1m, $\Delta t < 0.0002$ seconds. Since tunnel passage takes several seconds, tens of thousands of time steps are needed.
Since trains run close to the ground, ground effect is important.
- Moving Ground: Wall condition moving at the same speed as the train
- Ballast Track Bed: Set as a rough wall with roughness
- Bogie Fairing: Significant streamlining effect, can reduce resistance by 10--15%
- Inter-car Gap Cover: Reduces gap wind, cuts friction drag by 5--8%
The effect of bogie fairings is significant.
The N700S reduced aerodynamic drag by about 7% compared to the N700A by adopting full-skirt and bogie fairings. This is the result of shape optimization using CFD.
Tunnel Micro-Pressure Wave—The "Shock Wave" Created by the Shinkansen
When a Shinkansen enters a tunnel, the air in front of the train is compressed, and a weak shock wave called a "micro-pressure wave" radiates from the exit. This creates the "boom!" explosion sound, which became a serious noise problem for residents along the lines in the 1990s. The countermeasure adopted was a significant lengthening of the nose shape—the 500 Series Shinkansen's 15m-long nose is for this reason. The method of numerically analyzing compression waves inside tunnels using CFD and optimizing the nose shape is one of the most important challenges in railway aerodynamic CFD. This analysis requires handling unsteady compressible fluids, and the computational cost is correspondingly high.
Upwind Scheme
First-order upwind: Large numerical diffusion but stable. Second-order upwind: Improved accuracy but risk of oscillations. Essential for high Reynolds number flows.
Central Differencing
Second-order accurate, but numerical oscillations occur for Pe number > 2. Suitable for low Reynolds number, diffusion-dominated flows.
TVD Scheme (MUSCL, QUICK, etc.)
Maintains high accuracy while suppressing numerical oscillations via limiter functions. Effective for capturing shock waves or steep gradients.
Finite Volume Method vs Finite Element Method
FVM: Naturally satisfies conservation laws. Mainstream in CFD. FEM: Advantageous for complex shapes and multi-physics. Mesh-free methods like SPH are also developing.
CFL Condition (Courant Number)
Explicit method: CFL ≤ 1 is the stability condition. Implicit method: Stable even for CFL > 1, but affects accuracy and iteration count. LES: CFL ≈ 1 is recommended. Physical meaning: Information should not travel more than one cell per time step.
Residual Monitoring
Convergence is judged when residuals for the continuity equation, momentum, and energy decrease by 3-4 orders of magnitude. The mass conservation residual is particularly important.
Relaxation Factor
Pressure: 0.2-0.3, Velocity: 0.5-0.7 are typical initial values. If diverging, lower the relaxation factor. After convergence, increase to accelerate.
Internal Iterations for Unsteady Calculations
Related Topics
なった
詳しく
報告