# Ideal Gas Compression in an Actuating Piston

In this application, AcuSolve is used to simulate pressure and temperature inside an actuating piston using the ideal gas relationship and fully defined mesh motion. AcuSolve results are compared with analytical results as described in Moran and Shapiro (2000). The close agreement of AcuSolve results with analytical results validates the ability of AcuSolve to model cases with material properties defined by the ideal gas law subjected to significant mesh distortion.

## Problem Description

The problem consists of air compressed by a piston within a cylinder 10 m long with a diameter of π (3.14159) m, as shown in the following image, which is not drawn to scale. The piston is driven by a crankshaft 4 m long attached to a connecting rod 6 m long. The crankshaft rotates at 60o/s around a center point. At time (t) = 0.0 s the crank and piston are at top dead center (TDC) with an initial cylinder volume of 80 m3. At t = 3.0 s the piston is at bottom dead center (BDC), at which point the volume in the cylinder is reduced to 16 m3. The piston then moves back to TDC at t = 6.0 s. Air at atmospheric conditions, with material properties of an ideal gas and a dynamic viscosity of 1.781 X 105 kg/m-sec is the fluid modeled in the domain.

The motion of the piston is described by the following equation.

$x=r*\mathrm{cos}A+\sqrt{{l}^{2}-{r}^{2}{\mathrm{sin}}^{2}A}$
Where
• x is the distance of the piston from TDC
• r is the length of the connecting rod, 4 m
• A is the crank angle as a function of time (t)
• l is the connecting rod length, 6 m
During the simulation, the top surface of the cylinder, representing the piston, moves with an applied mesh motion that mimics the motion of the reciprocating crank, rod, and piston assembly. The ideal gas is compressed by contracting the mesh according to the motion of that assembly and by allowing the fluid properties to change as a function of time.

## AcuSolve Results

The simulation was run as transient for a total of 360 time steps with a duration of 6.0 seconds. The AcuSolve solution shows time dependent pressure and temperature in a cylinder volume that compresses and expands in a 6.0 second cycle. Pressure and temperature increase as the cylinder volume decreases and returns toward initial conditions as the cylinder volume expands. The contours of pressure and temperature within the cylinder at the initial stage (t = 0.0 s) and the half cycle (t = 3.0 s) are shown in the following images.
Comparison of time history plots of the AcuSolve and analytical solutions shows that the AcuSolve solution is nearly identical to the analytical solution. The AcuSolve computed internal pressure reaches a maximum of 8.62 X 105 Pa and the temperature reaches a maximum of 246.45 degrees C.

## Summary

The AcuSolve solution compares well with analytical results for ideal gas compression in an actuating piston. The AcuSolve transient solution compares nearly identically with the analytical solution, with 0.1 percent and 0.16 percent error for the maximum pressure and temperature, respectively. As a result, AcuSolve demonstrates the ability to predict ideal gas compression in an idealized actuating piston.

## Simulation Settings for Ideal Gas Compression in an Actuating Piston

HyperMesh CFD database file: <your working directory>\cylinder_piston_compression\cylinder_piston_compression.hm

Global

• Problem Description
• Analysis type - Transient
• Temperature equation - Advective Diffusive
• Turbulence equation - Laminar
• Mesh type - Fully Specified
• Auto Solution Strategy
• Max time steps - 360
• Initial time increment - 0.01667 sec
• Convergence tolerance - 1.0e-5 sec
• Max stagger iterations - 10
• Multiplier Function
• Piston
• Type - Piecewise Linear
• Curve fit variable - Time
• Curve fit values - (included in database as defined above)
• Material Model
• Air
• Density
• Type - Ideal Gas
• Gas constant - 286.9 J/kg-K
• Abs. pressure offset - 101325.0 N/m2
• Abs. temperature offset - 273.15 K
• Viscosity - 1.781e-5 kg/m-sec
• Mesh Motion
• Piston
• Type - Translation
• Translation velocity - 0.0, -1.0, 0.0
• Velocity variable - Multiplier Function
• Translation multiplier function - Piston

Model

• Volumes
• Fluid
• Element Set
• Material model - Air
• Element Output - on
• Surfaces
• BDC
• Simple Boundary Condition
• Type - Wall
• Wall velocity type - Match Mesh Velocity
• Mesh displacement BC type - Fixed
• SideWall
• Simple Boundary Condition
• Type - Wall
• Wall velocity type - Match Mesh Velocity
• Mesh displacement BC type - Fixed
• TDC
• Simple Boundary Condition
• Type - Wall
• Mesh displacement BC type - Fixed
• Mesh motion - Piston
• Nodes
• VolumeNodes
• Mesh X-Displacement
• Type - Zero
• Mesh Y-Displacement
• Type - Linear
• Mesh motion - Piston
• Curve fit variable - Y reference coordinate
• X1 - 10.0 m
• Y1 - 1.0 m
• Y2 - 0.0 m
• Mesh Z-Displacement
• Type - Zero

## References

M.J. Moran and H.N Shapiro. "Fundamentals of Engineering Thermodynamics". Wiley. 4th Ed. 2000.