Tutorials › Calculus III › Pressure-Vessel Mass Under Volume And Stress Limits

Applied Multivariable Optimization · Tutorial 493 of 999

Pressure-Vessel Mass Under Volume And Stress Limits

Derive the stress-limited mass model and determine what the capacity constraint implies for the vessel’s least possible mass.

Advanced 10 min read

What You'll Learn

  • Write the volume and thin-wall surface area of a cylinder with hemispherical ends
  • Translate allowable hoop stress into a minimum wall thickness
  • Reduce the mass formula using the pressure and stress constraints
  • Find the least mass among capacity-equality designs
  • Explain why excess-capacity designs do not invalidate the limiting result
  • Distinguish an unattained infimum from a minimum

Introduction

A pressure vessel must hold a specified amount of fluid without exceeding the stress that its material can safely withstand. Increasing the wall thickness reduces stress but adds mass; changing the radius and cylindrical length changes both the capacity and the required amount of material.

We will model a thin-walled cylinder closed by two hemispherical ends. The design variables are the cylinder radius, cylindrical length, and wall thickness. The aim is to minimize modeled mass subject to a required capacity and an allowable-stress limit. Under the idealized model used here, the analysis has a subtle outcome: there is a finite lower bound on mass, but no design attains it.

Geometry and Assumptions

Let \(r>0\) be the vessel’s internal radius, \(L\geq 0\) the length of its cylindrical section, and \(t>0\) the uniform wall thickness. The vessel’s two hemispherical ends form one sphere. Thus the internal capacity is the sum of the cylindrical and spherical volumes.

$$ V(r,L)=\pi r^2L+\frac{4}{3}\pi r^3. $$

The cylindrical wall has lateral area \(2\pi rL\), and the two hemispherical ends have the area of a sphere, \(4\pi r^2\). In a thin-wall model, the material volume is approximated by wall area times thickness. If the material density is \(\rho>0\), the modeled mass is therefore:

$$ m(r,L,t)=\rho t\left(2\pi rL+4\pi r^2\right). $$

Let \(V_0>0\) be the required capacity, \(p>0\) the internal pressure, and \(S>0\) the allowable tensile stress. We use the standard thin-wall pressure-vessel stress model: the cylinder’s hoop stress is \(pr/t\); its longitudinal stress is \(pr/(2t)\); and the spherical ends’ membrane stress is \(pr/(2t)\). The hoop stress is the largest of these, so satisfying its limit also satisfies the other two in this model.

Model constraints: A design is feasible when \(r>0\), \(L\geq 0\), \(t>0\), \(V(r,L)\geq V_0\), and \(pr/t\leq S\). The model assumes thin walls, uniform thickness and material, and the stated idealized stress formulas. It does not impose a maximum vessel length, a minimum radius, or manufacturing constraints.

Use the Stress Limit to Choose Thickness

For fixed \(r\) and \(L\), mass increases with thickness because the surface area and density are positive. The hoop-stress condition is equivalent to \(t\geq pr/S\). Consequently, the least-mass thickness for those dimensions is the smallest allowable one:

$$ t_{\min}(r)=\frac{pr}{S}. $$

Substituting this thickness into the mass model gives the least mass associated with each pair \((r,L)\):

$$ m(r,L)=\frac{\rho p}{S}\left(2\pi r^2L+4\pi r^3\right). $$

Here and below, \(m(r,L)\) means the mass with the stress constraint active. Any thicker wall has greater mass for the same dimensions and cannot improve the objective.

Capacity Equality and Its Limitation

For a fixed radius, a longer cylinder has more capacity and more wall area. Let \(L_0\) be the length that gives capacity exactly \(V_0\) at radius \(r\). When \(L_0\ge0\) and a feasible design has \(L>L_0\), shortening it to \(L_0\) preserves the radius and thickness, keeps the required capacity, and strictly decreases mass; if \(L=L_0\), the mass is unchanged. Solving \(V(r,L)=V_0\) for \(L\) gives:

$$ L_0=\frac{V_0}{\pi r^2}-\frac{4r}{3}. $$

This length is physically allowed only when it is nonnegative. Equivalently, the radius must satisfy \(r\leq (3V_0/(4\pi))^{1/3}\). If \(r\) is larger, the hemispherical ends alone already exceed the required capacity; setting \(L=0\) cannot produce a capacity-equality design. Thus it is not valid simply to assume every feasible design can be shortened to capacity equality.

Deriving the Limiting Mass

When a capacity-equality design is possible, substitute its length into the stress-limited mass formula. The substitution simplifies the expression as follows:

$$ \begin{aligned} m(r) &=\frac{\rho p}{S}\left[ 2\pi r^2\left(\frac{V_0}{\pi r^2}-\frac{4r}{3}\right)+4\pi r^3 \right]\\ &=\frac{\rho p}{S}\left[ 2V_0-\frac{8}{3}\pi r^3+4\pi r^3 \right]\\ &=\frac{\rho p}{S}\left(2V_0+\frac{4}{3}\pi r^3\right). \end{aligned} $$

For positive \(r\), this is greater than \(2\rho pV_0/S\), and it approaches that value as \(r\) approaches zero. But the capacity constraint then requires a very long cylinder. The limiting dimensions do not form a finite vessel.

We can establish the same lower bound for every feasible design, including ones whose hemispherical ends alone exceed the target capacity. First, the stress constraint gives \(t\geq pr/S\). Since area is positive, this yields a lower bound on mass. Then use the capacity constraint to rearrange that bound:

$$ \begin{aligned} m &=\rho t(2\pi rL+4\pi r^2)\\ &\geq \frac{\rho p}{S}(2\pi r^2L+4\pi r^3)\\ &=\frac{\rho p}{S}\left(2V(r,L)+\frac{4}{3}\pi r^3\right)\\ &\geq \frac{\rho p}{S}\left(2V_0+\frac{4}{3}\pi r^3\right)\\ &> \frac{2\rho pV_0}{S}. \end{aligned} $$

The strict final inequality follows from \(r>0\). It applies whether or not the vessel has excess capacity. To see that the lower bound can be approached, choose capacity-equality designs with positive radii tending to zero. For sufficiently small \(r\), the length \(V_0/(\pi r^2)-4r/3\) is positive, and the mass formula tends to \(2\rho pV_0/S\).

Conclusion: Under this model, the mass has infimum \(2\rho pV_0/S\), but no feasible design attains that value. Every design has strictly greater mass. A true minimum would require additional restrictions, such as a lower bound on radius or an upper bound on length.

Worked Example: A Capacity-Equality Vessel

Worked Example: A Capacity-Equality Vessel

Take \(V_0=100\text{ m}^3\), \(p=1\text{ MPa}\), \(S=100\text{ MPa}\), \(\rho=8000\text{ kg/m}^3\), and choose \(r=2\text{ m}\). Find the length, thickness, and modeled mass for a vessel whose capacity is exactly the target.

First calculate the length from the capacity-equality formula:

$$ \begin{aligned} L &=\frac{100\text{ m}^3}{\pi(2\text{ m})^2}-\frac{4(2\text{ m})}{3}\\ &=\left(\frac{25}{\pi}-\frac{8}{3}\right)\text{ m}\\ &\approx 5.2911\text{ m}. \end{aligned} $$

This length is positive, so the dimensions are feasible. The stress-limited thickness is:

$$ t=\frac{p r}{S} =\frac{1\text{ MPa}}{100\text{ MPa}}(2\text{ m}) =0.020\text{ m}. $$

Using the simplified capacity-equality mass formula:

$$ \begin{aligned} m &=\frac{(8000\text{ kg/m}^3)(1\text{ MPa})}{100\text{ MPa}} \left(2(100\text{ m}^3)+\frac{4}{3}\pi(2\text{ m})^3\right)\\ &=80\text{ kg/m}^3\left(200+\frac{32\pi}{3}\right)\text{ m}^3\\ &=16000+\frac{2560\pi}{3}\text{ kg}\\ &\approx 18680.83\text{ kg}. \end{aligned} $$

Check the capacity by substituting the exact length: \(\pi(2\text{ m})^2L=100\text{ m}^3-\frac{32}{3}\pi\text{ m}^3\), while the ends contribute \(\frac{32}{3}\pi\text{ m}^3\). Their sum is \(100\text{ m}^3\), as required. The thickness gives hoop stress \(pr/t=(1\text{ MPa})(2\text{ m})/(0.020\text{ m})=100\text{ MPa}=S\). Both checks agree with the design conditions.

Worked Example: Comparing Two Feasible Radii

Worked Example: Comparing Two Feasible Radii

Keep the data from the previous example, but choose \(r=1\text{ m}\). Calculate the capacity-equality length and mass, then compare with the \(r=2\text{ m}\) design.

The length is:

$$ L=\frac{100\text{ m}^3}{\pi(1\text{ m})^2}-\frac{4(1\text{ m})}{3} =\left(\frac{100}{\pi}-\frac{4}{3}\right)\text{ m} \approx30.4977\text{ m}. $$

The required thickness is \(t=(1/100)(1\text{ m})=0.010\text{ m}\). The mass is:

$$ \begin{aligned} m &=80\text{ kg/m}^3\left(200+\frac{4\pi}{3}\right)\text{ m}^3\\ &=16000+\frac{320\pi}{3}\text{ kg}\\ &\approx16335.10\text{ kg}. \end{aligned} $$

To check capacity, the cylinder contributes \(100\text{ m}^3-\frac{4}{3}\pi\text{ m}^3\), and the ends contribute \(\frac{4}{3}\pi\text{ m}^3\); the total is \(100\text{ m}^3\). The stress check is \((1\text{ MPa})(1\text{ m})/(0.010\text{ m})=100\text{ MPa}\). This design is lighter than the \(r=2\text{ m}\) design, consistent with the increasing term \(4\pi r^3/3\) in the capacity-equality mass formula.

Worked Example: Approaching the Infimum

Worked Example: Approaching the Infimum

Use the same data with \(r=0.5\text{ m}\). This illustrates how a smaller radius reduces mass while making the cylinder longer.

The capacity-equality length is:

$$ L=\frac{100\text{ m}^3}{\pi(0.5\text{ m})^2}-\frac{4(0.5\text{ m})}{3} =\left(\frac{400}{\pi}-\frac{2}{3}\right)\text{ m} \approx126.6573\text{ m}. $$

The stress-limited thickness is \(t=(1/100)(0.5\text{ m})=0.005\text{ m}\). The mass becomes:

$$ \begin{aligned} m &=80\text{ kg/m}^3\left(200+\frac{\pi}{6}\right)\text{ m}^3\\ &=16000+\frac{40\pi}{3}\text{ kg}\\ &\approx16041.89\text{ kg}. \end{aligned} $$

For the capacity check, the cylindrical contribution is \(100\text{ m}^3-\frac{\pi}{6}\text{ m}^3\), and the spherical ends contribute \(\frac{\pi}{6}\text{ m}^3\), giving \(100\text{ m}^3\) in total. The hoop stress is \((1\text{ MPa})(0.5\text{ m})/(0.005\text{ m})=100\text{ MPa}\). The mass is closer to the lower bound \(2\rho pV_0/S=16000\text{ kg}\), but it remains greater; the corresponding cylinder is much longer.

Worked Example: Excess Capacity Does Not Break the Bound

Worked Example: Excess Capacity Does Not Break the Bound

With the same data, take \(r=3\text{ m}\) and \(L=0\). The vessel has no cylindrical section, only the two hemispherical ends. Check its capacity and stress-limited mass.

The capacity is:

$$ V(3,0)=\frac{4}{3}\pi(3\text{ m})^3=36\pi\text{ m}^3\approx113.10\text{ m}^3. $$

This exceeds \(100\text{ m}^3\), so the design meets the capacity requirement. Its thickness at the allowable hoop stress is \(t=(1/100)(3\text{ m})=0.030\text{ m}\). Its mass is:

$$ \begin{aligned} m &=8000\text{ kg/m}^3(0.030\text{ m})\,4\pi(3\text{ m})^2\\ &=8640\pi\text{ kg}\\ &\approx27143.36\text{ kg}. \end{aligned} $$

The spherical-shell stress is \(pr/(2t)=(1\text{ MPa})(3\text{ m})/(2(0.030\text{ m}))=50\text{ MPa}\), below the allowable stress. The general lower bound also applies: \(m\geq80\text{ kg/m}^3(2(100\text{ m}^3)+\frac{4}{3}\pi(3\text{ m})^3)>16000\text{ kg}\). This excess-capacity design is therefore feasible but does not challenge the limiting lower bound.

A Common Mistake

A tempting argument says, “If the vessel holds too much, shorten the cylinder until the capacity is exactly \(V_0\).” That step is not always possible. For example, when \(L=0\) and \(r=3\text{ m}\) in the preceding example, the ends already have capacity \(36\pi\text{ m}^3>100\text{ m}^3\). Shortening the cylindrical section cannot reduce the vessel below this capacity.

The correction is to use two complementary arguments. For designs where \(V_0/(\pi r^2)-4r/3\geq0\), shortening to the capacity-equality length can reduce mass. For all feasible designs, whether or not that length is nonnegative, the stress-and-capacity bound proves \(m>2\rho pV_0/S\). The bound is approached by a family of small-radius, long-cylinder designs, so it characterizes the infimum without incorrectly claiming that every design can be reduced to exact capacity.

What the Result Means for Design

The limiting result is a mathematical consequence of the assumptions, not a recommendation to build an indefinitely long vessel. As \(r\) decreases, the required thickness \(pr/S\) decreases in direct proportion, while the cylindrical length needed for fixed capacity grows approximately like \(1/r^2\). The model allows that length to grow without penalty, and it omits practical limits such as fabrication, transport, stability, connections, and minimum gauge thickness.

Adding any appropriate design restriction can change the optimization problem and may produce an attainable minimum. The first task in an applied optimization problem is therefore to identify the model’s constraints precisely. Here, the distinction between a minimum and an infimum is essential: the mass can be made arbitrarily close to its lower bound, but no finite positive radius reaches the bound.

Check Your Understanding

Use the model and reasoning above to answer these questions.

  1. Why does the allowable hoop-stress condition lead to the smallest mass when \(t=pr/S\)?
  2. For which radii is a capacity-equality length \(L\geq0\) possible?
  3. Why does the capacity-equality mass decrease as \(r\) decreases?
  4. How does the lower-bound proof handle vessels whose capacity exceeds the target even when \(L=0\)?
  5. Why is \(2\rho pV_0/S\) an infimum rather than an attained minimum in this model?