2017 ASME CFD Tutorial – Pressure Vessels and Piping Conference (PVP)

On Monday July 17, 2017, Sean McGuffie of Porter McGuffie, Inc. presented the tutorial “The Use of CFD in Design” at the ASME Pressure Vessels and Piping (PVP) Conference in Hawaii.  While this tutorial had the same title as the tutorial presented at the PVP Conference in Toronto in 2012, significant changes were incorporated.

These changes included giving the attendees the opportunity to review the 2012 tutorial to better understand the mathematical basis for the Navier-Stokes equations, online and outside of the presentation.  Removing the tedious mathematics from the presentation enabled Mr. McGuffie to include significantly more information on:

  • Turbulence modeling
  • The actual run process and judging convergence
  • Debugging models
  • The use of CFD in PVP-related applications, including coupling CFD with FEA

PowerPoint for session 1 is available here.

PowerPoint for session 2 is available here.

Energy & Process CFD Symposium 2017

PMI presented at the Energy & Process CFD Symposium in Houston on May 18, 2017. The meeting featured presentations from industry experts discussing the use of simulation for process engineering, energy, and oil and gas applications.

At the symposium, Anthony DeFilippo delivered a presentation, “The Use of Complex Physics Models within the Refining Industry.” Three case studies were presented, highlighting ways that Porter McGuffie combines complex physics models when solving problems. Specifically, the cases look at multi-phase interactions modeled using Eulerian Two-Phase (E2P), thermo-acoustics modeled with Detached Eddy Simulation (DES), and structural response modeled with finite element analysis (FEA.) The focus of these examples was not on how they were carried out, but rather, how CFD-derived information was used to support engineering decisions.

Download Featured PowerPoint

Video Card Options and Testing

Viewing and manipulating large models requires significant compute resources. For post-processing FE and CFD results three areas control the amount of wall time required to generate each image:  number of CPU’s, available system RAM and Video Cards. This post focuses on video cards with information provided from benchmarks performed by PMI.   It is shown that the latest generation of video cards provide significant improvements in performance over video cards that are a few years old.


Based on initial tests that we performed for benchmarking our machines we found large differences in machine performance depending on the hardware that was available. Additionally, with the release of Star CCM 10.02, CD-Adapco has started supporting GPUs. With measurable performance differences between machines with different hardware and the new support of GPUs we decided to perform additional benchmarks to determine a good machine configuration, primarily for pre and post-processing of CFD models.


Description of Cards Used for Benchmarks:
  1. Nvidia FX3800 (workstation card) (http://www.nvidia.com/object/product_quadro_fx_3800_us.html). This is the primary video card used by PMI for the last couple years. The cards have provided acceptable performance, but we were starting to get persistent OpenGL errors, primarily with SolidWorks, Hypermesh and Star-CCM+. Based on internet research, we believe this is because the card only supported OpenGL 3.1. 4gb Of memory.
  2. Nvidia GeForce GTX 660 (desktop card) (http://www.geforce.com/hardware/desktop-gpus/geforce-gtx-660). This is a great midrange video card with 2gb of memory.
  3. PNY Nvidia Quadro K4200 (workstation card) (https://www.pny.com/nvidia-quadro-k4200). This is a beast of a card with 4gb of memory.


The test model was about 34 million cells and had a file size of 15 Gb.   The test started by opening a Geometry scene followed by a velocity scene. A mesh scene was then opened and the model was rotated through the same increments using standard view controls in CCM+.

Below is a table of the recorded times with each video card.

                                                                  Version 10.04
FX3800 660 Windows Driver 660 Nvidia Driver K4200
Velocity Scene 3m16s 4m58s 2m45s 3:01
(Vel already open) Metal Temp Scene 2m55s 1m48s 1m48s(crashed) 1:50
MT, rotate, +x+y+z to -x+y-z 1m43s 1m28s (crashed) 0:02
MT, rotate, -x+y-z to +x+y 1m54s 1m25s (crashed) 0:02


CCM allows you to specify Nvidia Graphics or Windows.

The GTX660 with Nvidia driver would crash when trying to load the second scene. We could never figure out exactly why it wouldn’t use some of the pc’s resources for the left over memory needed.



You can see the K4200 is an amazing card. The initial load times are about the same but when it comes to rotation there’s nothing even close and yes that is 2 seconds. By the time I started the timer and looked up the rotation was complete. This is huge when exporting images to create a video. The code and sample video can be found HERE. Workstation cards are designed for specific types of work. In this case the workstation cards are a must.


When installing multiple cards make sure you are removing all of the old drivers. When installing the new drivers there is an option for clean install. In the test computer I used a software called Display Driver Uninstaller. Use this program at your own risk.

By, James Greeson

Is a Square or Triangular Pitch Better in 2-Phase Flow?

PMI is frequently asked which tube arrangement is better for shell-and-tube heat exchangers with two-phase flow on the shell-side: square or triangular pitch?  The short answer is: that’s an interesting question.  To shed some light on the shell-side flow patterns within the tube bundle, let’s take a look at two models.

Geometrically, the models have the same percentage of flow space blocked by the tubes ~30%.

  • The test sections are designed to model the flow characteristics, with no heat transfer, near the center of a typical bundle.
  • The models are analyzed using two-phase flow physics with gravity.
    • The two species considered were: water (r ~ 900 kg/m3, m ~ 0.2 cP) and steam (r ~ 3 kg/m3, m ~ 0.015 cP).
    • The approach (superficial) velocity is 0.5 m/s for both phases, with a bulk light phase volume fraction of 37.5%. This bulk volume fraction is typical near the center of the bundle.

The pictures below show the heavy phase velocities from each model.


Heavy Velocities TriangleHeavy Velocities Square

By the defined flow paths, channeling is apparent in both bundle types.  A higher level of channeling occurs in the square bundle, as there is a direct path for the flow to travel from bottom to top.

The pictures below show the light volume fractions from each model.

Light VF TriangleLight VF Square

As can be seen in the images, the volume fraction distribution in the triangular bundle is very regular, with blanketing of the light fluid on the downstream sides of the tubes indicated by the high volume fraction in the vortex downstream from the tube.  In the square bundle the distributions are not as regular, but there seems to exist blanketing of the light fluid on the top and bottom of the tubes.  The flow for this case cannot be characterized as steady, due to the irregularity of the volume fraction distributions.  Click here to see a short-animation of the quasi-steady convection of fluid through the square bundle.

While in this case it looks like the triangular bundle would exhibit less blanketing and the problems associated with that flow condition, several additional factors must be considered:

  • No heat transfer was modeled
  • A single flow condition was analyzed, and
  • The model was a simplified test section

Based on the small set of results presented here, it becomes obvious that a correct answer would need to address the following questions:

  • What are the volume fraction and velocity distributions within the bundle, because:
    • The flow is 3-dimensional
      • Flow between sections along the exchanger must be considered, and
      • The make-up water inflow into the bundle is highly dependent on bundle interaction with the shell
    • Phase-change must be considered
      • Where do 2-phase flow transitions (i.e., plug, slug, wispy, dispersed, etc.) occur, and
      • How does the shell-side flow regime affect the local heat transfer coefficients

As can be seen above, another answer to the initial question could be – it’s complicated.  The best answer to the question can be achieved through conducting CFD analyses on the whole exchanger.  PMI has conducted more than 20 such analyses, primarily on heat exchanger equipment used in sulfur recover processes.  For further information on these types of analyses please click the “Contact” button above.

By Sean McGuffie

Using Java Script To Create Animations Within Star CCM+

Adding motion to your simulation results can produce better understanding for all interested parties. Through exploiting the Java macro functionality offered in Star-CCM+ the camera’s position can be incrementally adjusted through multiple steps, while an image is saved at each step.  Using this functionality the model and results can be rotated, panned and zoomed automatically.  The saved images can then be stitched together using a wide variety of animation programs (Windows Movie Maker, ImageMagick, etc…).

PMI has developed a Java script that links to a local coordinate system within the Star-CCM+ model.  Using this coordinate system a wide variety of scene transformations can be implemented to create animations.  PMI has provided a copy of this script with this week’s update[1].  This script include a graphical user interface (GUI) to allow for most of the model transformations to be coded without the user needing to modify the underlying Java code.

The basics for using the script are as follows:

  • Copy the script into a text file and save it. Rename the file as ‘rotate_zoom_translate_GUI’ and change the extension from ‘.txt’ to ‘.java’.
  • Open your simulation file in Star CCM+. Run the ‘rotate_zoom_translate_GUI.java’ macro from ‘Play Macro…’ option.
  • Interactive dialogue boxes will pop-up to communicate information and will allow you to input variables for view transformations.

[1] NOTE:  This script was developed for PMI’s internal use and is provided without warranty.  PMI will not offer technical support for script modifications.

Java Code — Rotate_Zoom_Translate_GUI


By Lilunnahar Deju