Home >

## CFD Blog Feeds

### Another Fine Meshtop

► This Week in CFD
27 Mar, 2020
You get the bare bones version of CFD news today because WordPress’ post editor lost my virtually completed first draft and, honestly, it used up all my authoring for today. Here’s an eye-catching image from a Simulia article about tire … Continue reading
► Mesh Adaptation With Python, Pointwise and Caelus
25 Mar, 2020
Mesh adaptation is one of the NASA 2030 goals, and has been one of the areas of research and development by Pointwise in past case studies and webinars. Applied CCM is the distributor for Pointwise in Australia and New Zealand, … Continue reading
► Meshing From Home
24 Mar, 2020
The city of Fort Worth finally caught up with its neighbor to the east and issued stay-home and business closure orders that go into effect tomorrow in order to slow the spread of the virus that causes COVID-19. So like … Continue reading
► This Week in CFD
20 Mar, 2020
This week’s CFD news includes high-order, GPU-accelerated simulation of turbomachinery, a sports car that seems to use a crossflow jet in place of a windshield, an opportunity to win $20,000, and more. Shown here is a screen capture of a … Continue reading ► This Week in CFD 13 Mar, 2020 It seems this week that CFD stands for Covid-19 Fueled Disarray. But enough of that. Even though this week’s CFD news is a tad light, there are several really cool grid pictures. Shown here is a model of the lungs … Continue reading ► Survey Results: Mesh Types, Part 2 12 Mar, 2020 The Mesh Types survey included several questions with short-form answers and we shared that data last month. Here we have answers to the long-form, open-ended question about what types of meshes you said you wanted to generate. So let’s start … Continue reading ### F*** Yeah Fluid Dynamicstop ► Dunes Avoid Collisions 27 Mar, 2020 The speed at which a dune migrates depends on its size; smaller dunes move faster than larger ones. That speed differential implies that small dunes should frequently collide into and merge with larger dunes, eventually forming one giant dune rather than a field of smaller separate ones. But that’s not what we observe in nature. To figure out why dunes aren’t colliding that often, researchers built a dune field of their own in the form of a rotating water tank. Inside the tank, their two artificial dunes can chase one another indefinitely while the researchers observe their interactions. What they found is that the dunes “communicate” with one another through the flow. As flow moves over the upstream dune, it generates turbulence in its wake, which the downstream dune then encounters. All that extra turbulence affects how sediment is picked up and transported for the downstream dune, ultimately changing its migration speed. For two dunes of initially equal size and close spacing, these interactions push the downstream dune further away until the separation between the dunes is large enough that they both migrate at the same speed. Even between dunes of unequal sizes, the researchers found that these repulsive interactions force the dunes away from collision and into migration at the same speed. (Image credit: dune field – G. Montani, others – K. Bacik et al.; research credit: K. Bacik et al.; via Cosmos; submitted by Kam-Yung Soh) ► Sunlight Is Older Than You Think 26 Mar, 2020 Joe Hanson over at “It’s Okay to Be Smart” has a great video on the random walk photons have to make to escape the core of the sun and other stars. Because the high-energy photons born in the star’s core have to bounce their way out rather than flying in a straight line, those photons can spend thousands of years escaping the sun. After that, the eight-and-a-half minute trip to Earth is nothing. But there’s a key element missing in this explanation: convection! That radiative random walk photons do doesn’t last all the way from the core of the sun to its surface. From a depth of about 200,000 km onward, the dominant mode of transport in the sun is convection, actual fluid motion that carries heat and light much faster than simple molecular diffusion, or Brownian motion, does. That’s why the surface of the sun shines with convection cells similar to the ones you’ll see in your skillet when heating a layer of oil. Fluid motion beyond molecular diffusion is also a big part of the other flows Joe describes in the video. If you had to wait on Brownian motion in order to smell your morning coffee, it would be cold long before you knew it was there! (Video and image credit: It’s Okay to Be Smart; sun surface image credit: Big Bear Solar Observatory/NJIT) ► Kicking Droplets 25 Mar, 2020 Moving the surface a droplet sits on creates some interesting dynamics, especially if the surface is hydrophobic. That’s what we see here with these droplets launched off an impulsively-moved plate. On the left, the drop has some limited contact with the plate and it takes time for the droplet to completely detach. When accelerated, the droplet first flattens into a pancake, the rim of which quickly leaves the plate. The center of the droplet is slower to detach, stretching the drop into a vase-like shape. When the drop does finally lose contact, it creates a fast-moving jet that shoots upward at several meters per second! In contrast the image on the left shows a levitating Leidenfrost droplet. Since this drop has no physical contact with the plate, the kick makes it leave the surface all at once, launching a pancake-like drop that quickly forms unstable lobes. (Image and research credit: M. Coux et al.) ► Listening to a Bubble’s Pop 24 Mar, 2020 Sound is an important aspect of many flows, from the scream of a rocket engine to the hum of electrical wires vibrating in the wind. Critically, those sounds carry important information about the flow. A new study extends these acoustic diagnostics to the popping of soap bubbles. When a hole opens in a soap bubble, it throws the surface-tension-driven capillary forces of the bubble into disarray. The rim around the hole retracts, pushing fluid away from the expanding hole. At the same time, air is pushed out of the collapsing bubble. Using microphone arrays, the researchers found they could measure and distinguish sound from both sources — the escaping air and the expanding hole. From the sound, they developed a model that predicts the rupture location, bubble thickness profile, and other properties of the bubble. They confirmed the model’s results by comparing with high-speed photography. The authors hope their new acoustic technique will shed light on bubble bursting events that are hard to observe visually, like the bubbling of magma. (Image and research credit: A. Bussonnière et al.; via Science News; submitted by Kam-Yung Soh) ► The Cricket’s Chirp 23 Mar, 2020 Growing up, my summer nights often featured a chorus of crickets and bull frogs. Even now, the sound of those chirps reminds me of home. So how do crickets make their calls? As this video shows, it’s a matter of scraping the hard edge of one wing along a tiny series of spines, similar to the teeth of a comb, that sit on the other wing. How fast the cricket’s wings move affects how frequently the chirps are heard. Being cold-blooded, the insects’ speed is affected by the external temperature, which is why you can count cricket chirps to estimate the temperature. Essentially, the chemical reactions necessary to regulate wing movement are temperature-dependent, so colder crickets produce slower chirps. (Video and image credit: Deep Look) ► Frozen Wavelets 20 Mar, 2020 Photographer Eric Gross captured these surreal alpine landscapes in Colorado’s Rocky Mountains. Although the lake’s surface appears to have frozen waves, the prevailing theory is that these mounds and divots occur when snowdrifts form atop the lake, melt and refreeze. Over multiple melting and freezing cycles, the lake builds up with what appear to be wind-driven waves frozen in time. (Image credit: E. Gross; via Colossal) ### Symscapetop ► CFD Simulates Distant Past 25 Jun, 2019 There is an interesting new trend in using Computational Fluid Dynamics (CFD). Until recently CFD simulation was focused on existing and future things, think flying cars. Now we see CFD being applied to simulate fluid flow in the distant past, think fossils. CFD shows Ediacaran dinner party featured plenty to eat and adequate sanitation read more ► Background on the Caedium v6.0 Release 31 May, 2019 Let's first address the elephant in the room - it's been a while since the last Caedium release. The multi-substance infrastructure for the Conjugate Heat Transfer (CHT) capability was a much larger effort than I anticipated and consumed a lot of resources. This lead to the relative quiet you may have noticed on our website. However, with the new foundation laid and solid we can look forward to a bright future. Conjugate Heat Transfer Through a Water-Air Radiator Simulation shows separate air and water streamline paths colored by temperature read more ► Long-Necked Dinosaurs Succumb To CFD 14 Jul, 2017 It turns out that Computational Fluid Dynamics (CFD) has a key role to play in determining the behavior of long extinct creatures. In a previous, post we described a CFD study of parvancorina, and now Pernille Troelsen at Liverpool John Moore University is using CFD for insights into how long-necked plesiosaurs might have swum and hunted. CFD Water Flow Simulation over an Idealized Plesiosaur: Streamline VectorsIllustration only, not part of the study read more ► CFD Provides Insight Into Mystery Fossils 23 Jun, 2017 Fossilized imprints of Parvancorina from over 500 million years ago have puzzled paleontologists for decades. What makes it difficult to infer their behavior is that Parvancorina have none of the familiar features we might expect of animals, e.g., limbs, mouth. In an attempt to shed some light on how Parvancorina might have interacted with their environment researchers have enlisted the help of Computational Fluid Dynamics (CFD). CFD Water Flow Simulation over a Parvancorina: Forward directionIllustration only, not part of the study read more ► Wind Turbine Design According to Insects 14 Jun, 2017 One of nature's smallest aerodynamic specialists - insects - have provided a clue to more efficient and robust wind turbine design. Dragonfly: Yellow-winged DarterLicense: CC BY-SA 2.5, André Karwath read more ► Runners Discover Drafting 1 Jun, 2017 The recent attempt to break the 2 hour marathon came very close at 2:00:24, with various aids that would be deemed illegal under current IAAF rules. The bold and obvious aerodynamic aid appeared to be a Tesla fitted with an oversized digital clock leading the runners by a few meters. 2 Hour Marathon Attempt read more ### CFD Onlinetop ► Physical Aspects of Lagrangian Particle Tracking (DPM) 1 Mar, 2020 This blog is first of the two part blog related to DPM Modeling. It is not to be considered as a professional document. Nor is it meant for those who consider themselves to be expert in DPM. It is meant as a guide for the beginners who want to understand various aspects related to DPM Modeling. Though terminology is meant to be general but it is rather closer to the one used in Fluent. Of course, experts are also welcome to go through it and share their comments. Assume that there is a need to track the evolution of contaminant in a canal. There are at least two ways to do it. Place a few sensors at a few locations. Or attach a sensor to multiple boats and let them traverse the canal. The latter is called Lagrangian. The boats may or may not have any drives to control their motion. Now these boats may be very small and light, just like the paper boats from the childhood. These will go with the flow without disturbing the flow of the water in any noticeable manner. This is known as One-Way coupling in DPM or LPT terminology. Or, these could be boats with rotors or big in size or heavy in mass and, consequently, the flow of the canal will be affected due to their presence. Even the distribution of the contaminant could be affected. This is called Two-Way coupling. Now, imagine having these big or heavy boats instead of light ones. Or it could be mix as well. But, if there are too many boats in the canal, a wave generated by one will affect the motion of all the nearby boats. That is three-way coupling. And if their number increases so much that it feels like traffic on the road, boats hitting each other, literally, to the extent that its not the flow rather the hits that decide the motion, you have what is called four-way coupling. This is where models like DEM come into play. Some users think of DEM as a separate model than DPM. The reality is that particles are always tracked using DPM but when there are so many of these that the forces generated due to their hits with each other affect their motion, some model is required to predict this force. DEM predicts that force. Lets consider one of these heavy boats. What are the bare minimum forces acting on it. 1. Its mass or weight, if you like the latter term, pulling it towards the center of the earth 2. The resistance from the air as well as water, only water if it is a submarine. This is called drag force and depends on various aspects, such as, contact area with the fluid, surface roughness, shape, and presence of other boats or submarines around it. The method to address three-way coupling is actually addition of extra forces or modification of drag due to this presence of extra particles around the particle being tracked. And there is lift along with the drag; always normal to the drag. 3. A large wave could also push the boat in a particular direction. Essentially, a wave brings a pressure difference because of its height and boat is pushed from higher pressure to lower pressure. This is called the pressure gradient force. 4. Other types of forces could be generated by rolling of the boat, its rotor at the back, or if it being towed. There is buoyancy that is pushing it up. For particles, that are not as big as boats, even smaller than paper boats, other forces also need to be considered. But before we go into that, lets take a look at the shape of the particles. As far as LPT is concerned, particles are assumed to have spherical shape. But particles are also termed as point particles. So, how is spherical shape important. Or how is shape important at all. To rephrase it, how is the shape taken into account. Its via the exchange coefficients, the models that define how a particle and the continuous fluid around it exchange various fields, such as, momentum, mass, thermal energy, etc. If drag model for a sphere is used, then it is obvious that the user is assuming the particles to be sphere. But lets say that the user wants to assume the particles to be oblong spheres or cylinders. Can he use drag coefficient for the cylinder or for oblong sphere? Ideally, yes. Practically, there are a few issues. Since LPT tracks only the positions of the particles, their orientations are unknown. Drag for a cylinder will certainly depend on its orientation. As far as sphere is concerned, it is invariant under rotations, i.e., it does not matter about which axis you rotate it, as long as the axis is passing through its center, its orientation will not vary. So, a single drag coefficient can be applied to all particles. For oblong spheres, their orientation is required but is unavailable. One plausible approach would be to use some probability distribution function for the orientation of the particles and apply drag or other forces based on that. But usually, it is not justified because the particles are much smaller than the characteristic length scale of the flow. Mass and heat transfer are also affected by the orientation. Though the area for drag force, or reaction, or heat convection remains the same, Nusselt and Sherwood numbers do not. They are affected by the orientation just like drag and lift coefficient. Therefore, particles are usually assumed or expected to be spherical or slightly non-spherical for applicability of DPM. Another important aspect is the volume fraction. Remember the waves from one boat affecting the flow of nearby boats. Well, standard DPM is valid until that happens. In other words, DPM is valid for use if the volume occupied by the particles is so less that they do not affect their neighbors in any significant manner. Usually, it is assumed that if the particles occupy less than 10% volume, then the LPT is valid. Do note that this condition needs to be satisfied locally as well as globally, i.e., if one cell has one particle then the cell volume should be at least 10 times the particle volume. There is no such limit on mass loading or mass fraction though. Now, these particles have to come from somewhere. That is what is called Injection. The injections require the important parameters to be specified, such as, initial position and velocity vectors, diameter for the calculation of mass and area (remember it is a sphere), total mass flow rate, temperature (if temperature is important for the simulation and thermal energy equation is being solved), and time duration for the injection. Now, there is no number specified here. CFD tools determine this number by taking a ratio of mass flow rate and mass of each particle, which is based on material property and diameter. Usually, this turns out to be in millions. Since each particle requires one equation, and this will be captured in the part two, millions of equations need to be solved. Now, that's expensive. Solution is to use representative particles, or otherwise known as parcels. Each parcel can depict any number of particles, varying from 1 to hundreds of thousands. This makes the situation workable. But what are basis for this? There are particles that have similar momentum or similar mass or similar position, or some other similarity that can be exploited to track them together. And that's how parcels are tracked. This number can be controlled by the user. So, once the tool knows the initial conditions, provided via injections, and it knows the forces acting on it, all that remains is to solve a rather simple ODE, called Newton's second law. As you observe, the equation has only one independent variable, time. So, the particle tracking is always transient. Solution procedure for this and the setup for DPM in Fluent will make the second part of this blog. Do not expect a tutorial, rather aspects that need to be looked after while setting it up. Until next time... ► Quick notes on testing optimization flags with OpenFOAM et al 21 Jan, 2020 Greetings to all! Tobi sent me an email earlier today related to this and I might as well leave a public note as well, to share with everyone the suggestions I had... so this is a quick copy-paste-adapt for future reference, until I or anyone else bothers with writing this at openfoamwiki.net I have no idea yet for the current generation of Ryzen CPUs (Ryzen 3000 series), but I do know of this report for EPYC: http://www.prace-ri.eu/best-practice-guide-amd-epyc If you look for Table 5, you will see the options they suggest for GCC/G++. However, the "znver1" architecture is possibly not the best for this generation of Ryzen/Threadripper... there is an alternative, which is to use: Code: -march=native -mtune=native It will only work properly with a recent GCC version for the more recent CPUs. Beyond this, it might take some trial and error. Some guidelines are given here: https://wiki.gentoo.org/wiki/GCC_optimization You can use the following strategy to test various builds with different optimization flags: 1. Go into the folder "wmake/rules/linux64Gcc" 2. Copy the files "cOpt" and "c++Opt" to another name, for example: "cOptNative" and "c++OptNative" 3. "cOPT" and "c++OPT" are the lines that need to be updated. 4. Create a new alias in your ".bashrc" file for this, for example: Code: alias ofdevNative='source$HOME/OpenFOAM/OpenFOAM-dev/etc/bashrc WM_COMPILE_OPTION=OptNative'
5. Start a new terminal and activate this alias ofdevNative.
6. Then run ./Allwmake inside "OpenFOAM-dev".
7. Repeat the same strategy for other names and therefore you can do several builds with small changes to the optimization flags.

Warning: Last time I checked, AVX and AVX2 are not used by OpenFOAM, so don't bother with them.

Best regards,
Bruno
► Mixing of Ammoniak and Exhaust
16 Aug, 2019
Dear Foamers,

in my thesis I worked with static mixers.
If you like to see my case you can see it here.
https://www.dropbox.com/sh/5rndjj0qs...Wci0dlNqa?dl=0
► Determination of mixing quality/ uniformity index
16 Aug, 2019
Dear guys,

for a long time I had problems determining the mixing quality of a mixing line. Now I've come across a usable formula. I would like to share this with you.
It is the degree of uniformity also called uniformity index.
The calculation is cell-based.
U = 1 - (SUM^{N}_{i=1}(Ci-Cm))/(2*N*Cm)
with N cells
and concentration of a cell Ci
and the arythmetic agent Cm
Cm = (SUM^{N}_{i=1}(Ci))/N
The easiest way is to export the cells with concentration of the considered region (outlet) and create an Excel file.
An example is shown in my public dropbox:
Greetings Philipp
► Connecting Fortran with VTK - the MPI way
24 May, 2019
I wrote a little couple of programs, respectively in Fortran and C++, as a proof of concept for connecting a Fortran program to a sort of visualization server based on VTK. The nice thing is that it uses MPI for the connection, so on the Fortran side nothing new and scary.

The code (you can find it at https://github.com/plampite/vtkForMPI) and the idea strongly predate a similar example in Using Advanced MPI by W. Gropp et al., but makes it more concrete by adding actual visualization based on VTK.

Of course, this is just a proof of concept, and nothing really interesting is really visualized (just a cylinder with parameters passed from Fortran side), but it is intended as an example to adapt to particular use cases (the VTK itself is taken from https://lorensen.github.io/VTKExamples/site/, where a lot of additional examples are present).
► Direct Numerical Simulation on a wing profile
14 May, 2019
Direct Numerical Simulation on a wing profile

1 billion points DNS (Direct Numerical Simulation) on a NACA4412 profile at 5 degrees angle of attack. Reynolds number is 350000 per airfoil chord and Mach number is 0.117. Both upper and lower turbulent boundary layers are tripped respectively at 15% and 50% by roughness elements evenly spaced in the boundary layer created by a zonal immersed boundary condition (Journal of Computational Physics, Volume 363, 15 June 2018, Pages 231-255, https://www.sciencedirect.com/science...). The spanwise extent is 0.3*chord. The computation has been performed on a structured multiblock mesh with the FastS compressible flow solver developed by ONERA on 1064 MPI cores. The video shows the early stages of the calculation (equivalent to 40000 time steps) highlighting the spatial development of fine-scale turbulence in both attached boundary layer and free wake. Post-processing and flow images have been made with Cassiopée (http://elsa.onera.fr/Cassiopee).

### curiosityFluidstop

► Creating curves in blockMesh (An Example)
29 Apr, 2019

In this post, I’ll give a simple example of how to create curves in blockMesh. For this example, we’ll look at the following basic setup:

As you can see, we’ll be simulating the flow over a bump defined by the curve:

$y=H*\sin\left(\pi x \right)$

First, let’s look at the basic blockMeshDict for this blocking layout WITHOUT any curves defined:

/*--------------------------------*- C++ -*----------------------------------*\
=========                 |
\\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
\\    /   O peration     | Website:  https://openfoam.org
\\  /    A nd           | Version:  6
\\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
version     2.0;
format      ascii;
class       dictionary;
object      blockMeshDict;
}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 1;

vertices
(
(-1 0 0)    // 0
(0 0 0)     // 1
(1 0 0)     // 2
(2 0 0)     // 3
(-1 2 0)    // 4
(0 2 0)     // 5
(1 2 0)     // 6
(2 2 0)     // 7

(-1 0 1)    // 8
(0 0 1)     // 9
(1 0 1)     // 10
(2 0 1)     // 11
(-1 2 1)    // 12
(0 2 1)     // 13
(1 2 1)     // 14
(2 2 1)     // 15
);

blocks
(
hex (0 1 5 4 8 9 13 12) (20 100 1) simpleGrading (0.1 10 1)
hex (1 2 6 5 9 10 14 13) (80 100 1) simpleGrading (1 10 1)
hex (2 3 7 6 10 11 15 14) (20 100 1) simpleGrading (10 10 1)
);

edges
(
);

boundary
(
inlet
{
type patch;
faces
(
(0 8 12 4)
);
}
outlet
{
type patch;
faces
(
(3 7 15 11)
);
}
lowerWall
{
type wall;
faces
(
(0 1 9 8)
(1 2 10 9)
(2 3 11 10)
);
}
upperWall
{
type patch;
faces
(
(4 12 13 5)
(5 13 14 6)
(6 14 15 7)
);
}
frontAndBack
{
type empty;
faces
(
(8 9 13 12)
(9 10 14 13)
(10 11 15 14)
(1 0 4 5)
(2 1 5 6)
(3 2 6 7)
);
}
);

// ************************************************************************* //

This blockMeshDict produces the following grid:

It is best practice in my opinion to first make your blockMesh without any edges. This lets you see if there are any major errors resulting from the block topology itself. From the results above, we can see we’re ready to move on!

So now we need to define the curve. In blockMesh, curves are added using the edges sub-dictionary. This is a simple sub dictionary that is just a list of interpolation points:

edges
(
polyLine 1 2
(
(0	0       0)
(0.1	0.0309016994    0)
(0.2	0.0587785252    0)
(0.3	0.0809016994    0)
(0.4	0.0951056516    0)
(0.5	0.1     0)
(0.6	0.0951056516    0)
(0.7	0.0809016994    0)
(0.8	0.0587785252    0)
(0.9	0.0309016994    0)
(1	0       0)
)

polyLine 9 10
(
(0	0       1)
(0.1	0.0309016994    1)
(0.2	0.0587785252    1)
(0.3	0.0809016994    1)
(0.4	0.0951056516    1)
(0.5	0.1     1)
(0.6	0.0951056516    1)
(0.7	0.0809016994    1)
(0.8	0.0587785252    1)
(0.9	0.0309016994    1)
(1	0       1)
)
);

The sub-dictionary above is just a list of points on the curve $y=H\sin(\pi x)$. The interpolation method is polyLine (straight lines between interpolation points). An alternative interpolation method could be spline.

The following mesh is produced:

Hopefully this simple example will help some people looking to incorporate curved edges into their blockMeshing!

Cheers.

This offering is not approved or endorsed by OpenCFD Limited, producer and distributor of the OpenFOAM software via http://www.openfoam.com, and owner of theOPENFOAM®  andOpenCFD®  trademarks.

► Creating synthetic Schlieren and Shadowgraph images in Paraview
28 Apr, 2019

Experimentally visualizing high-speed flow was a serious challenge for decades. Before the advent of modern laser diagnostics and velocimetry, the only real techniques for visualizing high speed flow fields were the optical techniques of Schlieren and Shadowgraph.

Today, Schlieren and Shadowgraph remain an extremely popular means to visualize high-speed flows. In particular, Schlieren and Shadowgraph allow us to visualize complex flow phenomena such as shockwaves, expansion waves, slip lines, and shear layers very effectively.

In CFD there are many reasons to recreate these types of images. First, they look awesome. Second, if you are doing a study comparing to experiments, occasionally the only full-field data you have could be experimental images in the form of Schlieren and Shadowgraph.

Without going into detail about Schlieren and Shadowgraph themselves, primarily you just need to understand that Schlieren and Shadowgraph represent visualizations of the first and second derivatives of the flow field refractive index (which is directly related to density).

In Schlieren, a knife-edge is used to selectively cut off light that has been refracted. As a result you get a visualization of the first derivative of the refractive index in the direction normal to the knife edge. So for example, if an experiment used a horizontal knife edge, you would see the vertical derivative of the refractive index, and hence the density.

For Shadowgraph, no knife edge is used, and the images are a visualization of the second derivative of the refractive index. Unlike the Schlieren images, shadowgraph has no direction and shows you the laplacian of the refractive index field (or density field).

In this post, I’ll use a simple case I did previously (https://curiosityfluids.com/2016/03/28/mach-1-5-flow-over-23-degree-wedge-rhocentralfoam/) as an example and produce some synthetic Schlieren and Shadowgraph images using the data.

#### So how do we create these images in paraview?

Well as you might expect, from the introduction, we simply do this by visualizing the gradients of the density field.

In ParaView the necessary tool for this is:

Once you’ve selected this, we then need to set the properties so that we are going to operate on the density field:

To do this, simply set the “Scalar Array” to the density field (rho), and change the name of the result Array name to SyntheticSchlieren. Now you should see something like this:

There are a few problems with the above image (1) Schlieren images are directional and this is a magnitude (2) Schlieren and Shadowgraph images are black and white. So if you really want your Schlieren images to look like the real thing, you should change to black and white. ALTHOUGH, Cold and Hot, Black-Body radiation, and Rainbow Desatured all look pretty amazing.

To fix these, you should only visualize one component of the Synthetic Schlieren array at a time, and you should visualize using the X-ray color preset:

The results look pretty realistic:

#### Vertical Knife Edge

The process of computing the shadowgraph field is very similar. However, recall that shadowgraph visualizes the Laplacian of the density field. BUT THERE IS NO LAPLACIAN CALCULATOR IN PARAVIEW!?! Haha no big deal. Just remember the basic vector calculus identity:

$\nabla^2\left[\right] = \nabla \cdot \nabla \left[\right]$

Therefore, in order for us to get the Shadowgraph image, we just need to take the Divergence of the Synthetic Schlieren vector field!

To do this, we just have to use the Gradient of Unstructured DataSet tool again:

This time, Deselect “Compute Gradient” and the select “Compute Divergence” and change the Divergence array name to Shadowgraph.

Visualized in black and white, we get a very realistic looking synthetic Shadowgraph image:

## So what do the values mean?

Now this is an important question, but a simple one to answer. And the answer is…. not much. Physically, we know exactly what these mean, these are: Schlieren is the gradient of the density field in one direction and Shadowgraph is the laplacian of the density field. But what you need to remember is that both Schlieren and Shadowgraph are qualitative images. The position of the knife edge, brightness of the light etc. all affect how a real experimental Schlieren or Shadowgraph image will look.

This means, very often, in order to get the synthetic Schlieren to closely match an experiment, you will likely have to change the scale of your synthetic images. In the end though, you can end up with extremely realistic and accurate synthetic Schlieren images.

Hopefully this post will be helpful to some of you out there. Cheers!

► Solving for your own Sutherland Coefficients using Python
24 Apr, 2019

Sutherland’s equation is a useful model for the temperature dependence of the viscosity of gases. I give a few details about it in this post: https://curiosityfluids.com/2019/02/15/sutherlands-law/

The law given by:

$\mu=\mu_o\frac{T_o + C}{T+C}\left(\frac{T}{T_o}\right)^{3/2}$

It is also often simplified (as it is in OpenFOAM) to:

$\mu=\frac{C_1 T^{3/2}}{T+C}=\frac{A_s T^{3/2}}{T+T_s}$

In order to use these equations, obviously, you need to know the coefficients. Here, I’m going to show you how you can simply create your own Sutherland coefficients using least-squares fitting Python 3.

So why would you do this? Basically, there are two main reasons for this. First, if you are not using air, the Sutherland coefficients can be hard to find. If you happen to find them, they can be hard to reference, and you may not know how accurate they are. So creating your own Sutherland coefficients makes a ton of sense from an academic point of view. In your thesis or paper, you can say that you created them yourself, and not only that you can give an exact number for the error in the temperature range you are investigating.

So let’s say we are looking for a viscosity model of Nitrogen N2 – and we can’t find the coefficients anywhere – or for the second reason above, you’ve decided its best to create your own.

By far the simplest way to achieve this is using Python and the Scipy.optimize package.

Step 1: Get Data

The first step is to find some well known, and easily cited, source for viscosity data. I usually use the NIST webbook (
https://webbook.nist.gov/), but occasionally the temperatures there aren’t high enough. So you could also pull the data out of a publication somewhere. Here I’ll use the following data from NIST:

 Temparature (K) Viscosity (Pa.s) 200 0.000012924 400 0.000022217 600 0.000029602 800 0.000035932 1000 0.000041597 1200 0.000046812 1400 0.000051704 1600 0.000056357 1800 0.000060829 2000 0.000065162

This data is the dynamics viscosity of nitrogen N2 pulled from the NIST database for 0.101 MPa. (Note that in these ranges viscosity should be only temperature dependent).

Step 2: Use python to fit the data

If you are unfamiliar with Python, this may seem a little foreign to you, but python is extremely simple.

First, we need to load the necessary packages (here, we’ll load numpy, scipy.optimize, and matplotlib):

import numpy as np
import matplotlib.pyplot as plt
from scipy.optimize import curve_fit

Now we define the sutherland function:

def sutherland(T, As, Ts):
return As*T**(3/2)/(Ts+T)

Next we input the data:

T=[200,
400,
600,
800,
1000,
1200,
1400,
1600,
1800,
2000]

mu=[0.000012924,
0.000022217,
0.000029602,
0.000035932,
0.000041597,
0.000046812,
0.000051704,
0.000056357,
0.000060829,
0.000065162]

Then we fit the data using the curve_fit function from scipy.optimize. This function uses a least squares minimization to solve for the unknown coefficients. The output variable popt is an array that contains our desired variables As and Ts.

popt = curve_fit(sutherland, T, mu)
As=popt[0]
Ts=popt[1]

Now we can just output our data to the screen and plot the results if we so wish:

print('As = '+str(popt[0])+'\n')
print('Ts = '+str(popt[1])+'\n')

xplot=np.linspace(200,2000,100)
yplot=sutherland(xplot,As,Ts)

plt.plot(T,mu,'ok',xplot,yplot,'-r')
plt.xlabel('Temperature (K)')
plt.ylabel('Dynamic Viscosity (Pa.s)')
plt.legend(['NIST Data', 'Sutherland'])
plt.show()

Overall the entire code looks like this:

import numpy as np
import matplotlib.pyplot as plt
from scipy.optimize import curve_fit

def sutherland(T, As, Ts):
return As*T**(3/2)/(Ts+T)

T=[200, 400, 600,
800,
1000,
1200,
1400,
1600,
1800,
2000]

mu=[0.000012924,
0.000022217,
0.000029602,
0.000035932,
0.000041597,
0.000046812,
0.000051704,
0.000056357,
0.000060829,
0.000065162]

popt, pcov = curve_fit(sutherland, T, mu)
As=popt[0]
Ts=popt[1]
print('As = '+str(popt[0])+'\n')
print('Ts = '+str(popt[1])+'\n')

xplot=np.linspace(200,2000,100)
yplot=sutherland(xplot,As,Ts)

plt.plot(T,mu,'ok',xplot,yplot,'-r')
plt.xlabel('Temperature (K)')
plt.ylabel('Dynamic Viscosity (Pa.s)')
plt.legend(['NIST Data', 'Sutherland'])
plt.show()



And the results for nitrogen gas in this range are As=1.55902E-6, and Ts=168.766 K. Now we have our own coefficients that we can quantify the error on and use in our academic research! Wahoo!

## Summary

In this post, we looked at how we can simply use a database of viscosity-temperature data and use the python package scipy to solve for our unknown Sutherland viscosity coefficients. This NIST database was used to grab some data, and the data was then loaded into Python and curve-fit using scipy.optimize curve_fit function.

This task could also easily be accomplished using the Matlab curve-fitting toolbox, or perhaps in excel. However, I have not had good success using the excel solver to solve for unknown coefficients.

► Tips for tackling the OpenFOAM learning curve
23 Apr, 2019

The most common complaint I hear, and the most common problem I observe with OpenFOAM is its supposed “steep learning curve”. I would argue however, that for those who want to practice CFD effectively, the learning curve is equally as steep as any other software.

There is a distinction that should be made between “user friendliness” and the learning curve required to do good CFD.

While I concede that other commercial programs have better basic user friendliness (a nice graphical interface, drop down menus, point and click options etc), it is equally as likely (if not more likely) that you will get bad results in those programs as with OpenFOAM. In fact, to some extent, the high user friendliness of commercial software can encourage a level of ignorance that can be dangerous. Additionally, once you are comfortable operating in the OpenFOAM world, the possibilities become endless and things like code modification, and bash and python scripting can make OpenFOAM worklows EXTREMELY efficient and powerful.

Anyway, here are a few tips to more easily tackle the OpenFOAM learning curve:

(1) Understand CFD

This may seem obvious… but its not to some. Troubleshooting bad simulation results or unstable simulations that crash is impossible if you don’t have at least a basic understanding of what is happening under the hood. My favorite books on CFD are:

(a) The Finite Volume Method in Computational Fluid Dynamics: An Advanced Introduction with OpenFOAM® and Matlab by
F. Moukalled, L. Mangani, and M. Darwish

(b) An introduction to computational fluid dynamics – the finite volume method – by H K Versteeg and W Malalasekera

(c) Computational fluid dynamics – the basics with applications – By John D. Anderson

(2) Understand fluid dynamics

Again, this may seem obvious and not very insightful. But if you are going to assess the quality of your results, and understand and appreciate the limitations of the various assumptions you are making – you need to understand fluid dynamics. In particular, you should familiarize yourself with the fundamentals of turbulence, and turbulence modeling.

(3) Avoid building cases from scratch

Whenever I start a new case, I find the tutorial case that most closely matches what I am trying to accomplish. This greatly speeds things up. It will take you a super long time to set up any case from scratch – and you’ll probably make a bunch of mistakes, forget key variable entries etc. The OpenFOAM developers have done a lot of work setting up the tutorial cases for you, so use them!

As you continue to work in OpenFOAM on different projects, you should be compiling a library of your own templates based on previous work.

(4) Using Ubuntu makes things much easier

This is strictly my opinion. But I have found this to be true. Yes its true that Ubuntu has its own learning curve, but I have found that OpenFOAM works seamlessly in the Ubuntu or any Ubuntu-like linux environment. OpenFOAM now has Windows flavors using docker and the like- but I can’t really speak to how well they work – mostly because I’ve never bothered. Once you unlock the power of Linux – the only reason to use Windows is for Microsoft Office (I guess unless you’re a gamer – and even then more and more games are now on Linux). Not only that- but the VAST majority of forums and troubleshooting associated with OpenFOAM you’ll find on the internet are from Ubuntu users.

I much prefer to use Ubuntu with a virtual Windows environment inside it. My current office setup is my primary desktop running Ubuntu – plus a windows VirtualBox, plus a laptop running windows that I use for traditional windows type stuff. Dual booting is another option, but seamlessly moving between the environments is easier.

(5) If you’re struggling, simplify

Unless you know exactly what you are doing, you probably shouldn’t dive into the most complicated version of whatever you are trying to solve/study. It is best to start simple, and layer the complexity on top. This way, when something goes wrong, it is much easier to figure out where the problem is coming from.

(6) Familiarize yourself with the cfd-online forum

If you are having trouble, the cfd-online forum is super helpful. Most likely, someone else is has had the same problem you have. If not, the people there are extremely helpful and overall the forum is an extremely positive environment for working out the kinks with your simulations.

(7) The results from checkMesh matter

If you run checkMesh and your mesh fails – fix your mesh. This is important. Especially if you are not planning on familiarizing yourself with the available numerical schemes in OpenFOAM, you should at least have a beautiful mesh. In particular, if your mesh is highly non-orthogonal, you will have serious problems. If you insist on using a bad mesh, you will probably need to manipulate the numerical schemes. A great source for how schemes should be manipulated based on mesh non-orthogonality is:

http://www.wolfdynamics.com/wiki/OFtipsandtricks.pdf

(8) CFL Number Matters

If you are running a transient case, the Courant-Freidrechs-Lewis (CFL) number matters… a lot. Not just for accuracy (if you are trying to capture a transient event) but for stability. If your time-step is too large you are going to have problems. There is a solid mathematical basis for this stability criteria for advection-diffusion problems. Additionally the Navier-Stokes equations are very non-linear and the complexity of the problem and the quality of your grid etc can make the simulation even less stable. When I have a transient simulation crash, if I know my mesh is OK, I decrease the timestep by a factor of 2. More often than not, this solves the problem.

For large time stepping, you can add outer loops to solvers based on the pimple algorithm, but you may end up losing important transient information. Excellent explanation of how to do this is given in the book by T. Holzmann:

https://holzmann-cfd.de/publications/mathematics-numerics-derivations-and-openfoam

For the record, this points falls into point (1) of Understanding CFD.

(9) Work through the OpenFOAM Wiki “3 Week” Series

If you are starting OpenFOAM for the first time, it is worth it to work through an organized program of learning. One such example (and there are others) is the “3 Weeks Series” on the OpenFOAM wiki:

https://wiki.openfoam.com/%223_weeks%22_series

If you are a graduate student, and have no job to do other than learn OpenFOAM, it will not take 3 weeks. This touches on all the necessary points you need to get started.

(10) OpenFOAM is not a second-tier software – it is top tier

I know some people who have started out with the attitude from the get-go that they should be using a different software. They think somehow Open-Source means that it is not good. This is a pretty silly attitude. Many top researchers around the world are now using OpenFOAM or some other open source package. The number of OpenFOAM citations has grown every year consistently (

In my opinion, the only place where mainstream commercial CFD packages will persist is in industry labs where cost is no concern, and changing software is more trouble than its worth. OpenFOAM has been widely benchmarked, and widely validated from fundamental flows to hypersonics (see any of my 17 publications using it for this). If your results aren’t good, you are probably doing something wrong. If you have the attitude that you would rather be using something else, and are bitter that your supervisor wants you to use OpenFOAM, when something goes wrong you will immediately think there is something wrong with the program… which is silly – and you may quit.

(11) Meshing… Ugh Meshing

For the record, meshing is an art in any software. But meshing is the only area where I will concede any limitation in OpenFOAM. HOWEVER, as I have outlined in my previous post (https://curiosityfluids.com/2019/02/14/high-level-overview-of-meshing-for-openfoam/) most things can be accomplished in OpenFOAM, and there are enough third party meshing programs out there that you should have no problem.

## Summary

Basically, if you are starting out in CFD or OpenFOAM, you need to put in time. If you are expecting to be able to just sit down and produce magnificent results, you will be disappointed. You might quit. And frankly, thats a pretty stupid attitude. However, if you accept that CFD and fluid dynamics in general are massive fields under constant development, and are willing to get up to speed, there are few limits to what you can accomplish.

Please take the time! If you want to do CFD, learning OpenFOAM is worth it. Seriously worth it.

This offering is notapproved or endorsed by OpenCFD Limited, producer and distributorof the OpenFOAM software via http://www.openfoam.com, and owner of theOPENFOAM®  andOpenCFD®  trade marks.

► Automatic Airfoil C-Grid Generation for OpenFOAM – Rev 1
22 Apr, 2019

Here I will present something I’ve been experimenting with regarding a simplified workflow for meshing airfoils in OpenFOAM. If you’re like me, (who knows if you are) I simulate a lot of airfoils. Partly because of my involvement in various UAV projects, partly through consulting projects, and also for testing and benchmarking OpenFOAM.

Because there is so much data out there on airfoils, they are a good way to test your setups and benchmark solver accuracy. But going from an airfoil .dat coordinate file to a mesh can be a bit of pain. Especially if you are starting from scratch.

The two main ways that I have meshed airfoils to date has been:

(a) Mesh it in a C or O grid in blockMesh (I have a few templates kicking around for this
(b) Generate a “ribbon” geometry and mesh it with cfMesh
(c) Or back in the day when I was a PhD student I could use Pointwise – oh how I miss it.

But getting the mesh to look good was always sort of tedious. So I attempted to come up with a python script that takes the airfoil data file, minimal inputs and outputs a blockMeshDict file that you just have to run.

The goals were as follows:
(a) Create a C-Grid domain
(b) be able to specify boundary layer growth rate
(c) be able to set the first layer wall thickness
(e) be mostly automatic (few user inputs)
(f) have good mesh quality – pass all checkMesh tests
(g) Quality is consistent – meaning when I make the mesh finer, the quality stays the same or gets better
(h) be able to do both closed and open trailing edges
(i) be able to handle most airfoils (up to high cambers)
(j) automatically handle hinge and flap deflections

In Rev 1 of this script, I believe I have accomplished (a) thru (g). Presently, it can only hand airfoils with closed trailing edge. Hinge and flap deflections are not possible, and highly cambered airfoils do not give very satisfactory results.

There are existing tools and scripts for automatically meshing airfoils, but I found personally that I wasn’t happy with the results. I also thought this would be a good opportunity to illustrate one of the ways python can be used to interface with OpenFOAM. So please view this as both a potentially useful script, but also something you can dissect to learn how to use python with OpenFOAM. This first version of the script leaves a lot open for improvement, so some may take it and be able to tailor it to their needs!

Hopefully, this is useful to some of you out there!

https://github.com/curiosityFluids/curiosityFluidsAirfoilMesher

Here you will also find a template based on the airfoil2D OpenFOAM tutorial.

## Instructions

(1) Copy curiosityFluidsAirfoilMesher.py to the root directory of your simulation case.
(2) Copy your airfoil coordinates in Selig .dat format into the same folder location.
(3) Modify curiosityFluidsAirfoilMesher.py to your desired values. Specifically, make sure that the string variable airfoilFile is referring to the right .dat file
(4) In the terminal run: python3 curiosityFluidsAirfoilMesher.py
(5) If no errors – run blockMesh

PS
You need to run this with python 3, and you need to have numpy installed

## Inputs

The inputs for the script are very simple:

ChordLength: This is simply the airfoil chord length if not equal to 1. The airfoil dat file should have a chordlength of 1. This variable allows you to scale the domain to a different size.

airfoilfile: This is a string with the name of the airfoil dat file. It should be in the same folder as the python script, and both should be in the root folder of your simulation directory. The script writes a blockMeshDict to the system folder.

DomainHeight: This is the height of the domain in multiples of chords.

WakeLength: Length of the wake domain in multiples of chords

firstLayerHeight: This is the height of the first layer. To estimate the requirement for this size, you can use the curiosityFluids y+ calculator

growthRate: Boundary layer growth rate

MaxCellSize: This is the max cell size along the centerline from the leading edge of the airfoil. Some cells will be larger than this depending on the gradings used.

The following inputs are used to improve the quality of the mesh. I have had pretty good results messing around with these to get checkMesh compliant grids.

BLHeight: This is the height of the boundary layer block off of the surfaces of the airfoil

inletGradingFactor: This is a grading factor that modifies the the grading along the inlet as a multiple of the leading edge grading and can help improve mesh uniformity

trailingBlockAngle: This is an angle in degrees that expresses the angles of the trailing edge blocks. This can reduce the aspect ratio of the boundary cells at the top and bottom of the domain, but can make other mesh parameters worse.

## Examples

### 12% Joukowski Airfoil

Inputs:

With the above inputs, the grid looks like this:

Mesh Quality:

These are some pretty good mesh statistics. We can also view them in paraView:

### Clark-y Airfoil

The clark-y has some camber, so I thought it would be a logical next test to the previous symmetric one. The inputs I used are basically the same as the previous airfoil:

With these inputs, the result looks like this:

Mesh Quality:

Visualizing the mesh quality:

### MH60 – Flying Wing Airfoil

Here is an example of a flying with airfoil (tested since the trailing edge is tilted upwards).

Inputs:

Again, these are basically the same as the others. I have found that with these settings, I get pretty consistently good results. When you change the MaxCellSize, firstLayerHeight, and Grading some modification may be required. However, if you just half the maxCell, and half the firstLayerHeight, you “should” get a similar grid quality just much finer.

Grid Quality:

Visualizing the grid quality

## Summary

Hopefully some of you find this tool useful! I plan to release a Rev 2 soon that will have the ability to handle highly cambered airfoils, and open trailing edges, as well as control surface hinges etc.

The long term goal will be an automatic mesher with an H-grid in the spanwise direction so that the readers of my blog can easily create semi-span wing models extremely quickly!

DISCLAIMER: This script is intended as an educational and productivity tool and starting point. You may use and modify how you wish. But I make no guarantee of its accuracy, reliability, or suitability for any use. This offering is not approved or endorsed by OpenCFD Limited, producer and distributor of the OpenFOAM software via http://www.openfoam.com, and owner of the OPENFOAM®  and OpenCFD®  trademarks.﻿

► Normal Shock Calculator
20 Feb, 2019

Here is a useful little tool for calculating the properties across a normal shock.

If you found this useful, and have the need for more, visit www.stfsol.com. One of STF Solutions specialties is providing our clients with custom software developed for their needs. Ranging from custom CFD codes to simpler targeted codes, scripts, macros and GUIs for a wide range of specific engineering purposes such as pipe sizing, pressure loss calculations, heat transfer calculations, 1D flow transients, optimization and more. Visit STF Solutions at www.stfsol.com for more information!

Disclaimer: This calculator is for educational purposes and is free to use. STF Solutions and curiosityFluids makes no guarantee of the accuracy of the results, or suitability, or outcome for any given purpose.

### Hanley Innovationstop

► Accurate Aircraft Performance Predictions using Stallion 3D
26 Feb, 2020

Stallion 3D uses your CAD design to simulate the performance of your aircraft.  This enables you to verify your design and compute quantities such as cruise speed, power required and range at a given cruise altitude. Stallion 3D is used to optimize the design before moving forward with building and testing prototypes.

The table below shows the results of Stallion 3D around the cruise angles of attack of the Cessna 402c aircraft.  The CAD design can be obtained from the OpenVSP hangar.

The results were obtained by simulating 5 angles of attack in Stallion 3D on an ordinary laptop computer running MS Windows 10 .  Given the aircraft geometry and flight conditions, Stallion 3D computed the CL, CD, L/D and other aerodynamic quantities.  With this accurate aerodynamics results, the preliminary performance data such as cruise speed, power, range and endurance can be obtained.

Lift Coefficient versus Angle of Attack computed with Stallion 3D

Lift to Drag Ratio versus True Airspeed at 10,000 feet

Power Required versus True Airspeed at 10,000 feet

The Stallion 3D results shows good agreement with the published data for the Cessna 402.  For example, the cruse speed of the aircraft at 10,000 feet is around 140 knots. This coincides with the speed at the maximum L/D (best range) shown in the graph and table above.

http://www.hanleyinnovations.com/stallion3d.html

Hanley Innovations is a pioneer in developing user friendly and accurate software that is accessible to engineers, designers and students.  For more information, please visit > http://www.hanleyinnovations.com

► 5 Tips For Excellent Aerodynamic Analysis and Design
8 Feb, 2020
Stallion 3D analysis of Uber Elevate eCRM-100 model

Being the best aerodynamics engineer requires meticulous planning and execution.  Here are 5 steps you can following to start your journey to being one of the best aerodynamicist.

1.  Airfoils analysis (VisualFoil) - the wing will not be better than the airfoil. Start with the best airfoil for the design.

2.  Wing analysis (3Dfoil) - know the benefits/limits of taper, geometric & aerodynamic twist, dihedral angles, sweep, induced drag and aspect ratio.

3. Stability analysis (3Dfoil) - longitudinal & lateral static & dynamic stability analysis.  If the airplane is not stable, it might not fly (well).

4. High Lift (MultiElement Airfoils) - airfoil arrangements can do wonders for takeoff, climb, cruise and landing.

5. Analyze the whole arrangement (Stallion 3D) - this is the best information you will get until you flight test the design.

Hanley Innovations is a pioneer in developing user friendly and accurate software the is accessible to engineers, designs and students.  For more information, please visit > http://www.hanleyinnovations.com

► Accurate Aerodynamics with Stallion 3D
17 Aug, 2019

Stallion 3D is an extremely versatile tool for 3D aerodynamics simulations.  The software solves the 3D compressible Navier-Stokes equations using novel algorithms for grid generation, flow solutions and turbulence modeling.

The proprietary grid generation and immersed boundary methods find objects arbitrarily placed in the flow field and then automatically place an accurate grid around them without user intervention.

Stallion 3D algorithms are fine tuned to analyze invisid flow with minimal losses. The above figure shows the surface pressure of the BD-5 aircraft (obtained OpenVSP hangar) using the compressible Euler algorithm.

Stallion 3D solves the Reynolds Averaged Navier-Stokes (RANS) equations using a proprietary implementation of the k-epsilon turbulence model in conjunction with an accurate wall function approach.

Stallion 3D can be used to solve problems in aerodynamics about complex geometries in subsonic, transonic and supersonic flows.  The software computes and displays the lift, drag and moments for complex geometries in the STL file format.  Actuator disc (up to 100) can be added to simulate prop wash for propeller and VTOL/eVTOL aircraft analysis.

Stallion 3D is a versatile and easy-to-use software package for aerodynamic analysis.  It can be used for computing performance and stability (both static and dynamic) of aerial vehicles including drones, eVTOLs aircraft, light airplane and dragons (above graphics via Thingiverse).

► Hanley Innovations Upgrades Stallion 3D to Version 5.0
18 Jul, 2017
The CAD for the King Air was obtained from Thingiverse

Stallion 3D is a 3D aerodynamics analysis software package developed by Dr. Patrick Hanley of Hanley Innovations in Ocala, FL. Starting with only the STL file, Stallion 3D is an all-in-one digital tool that rapidly validate conceptual and preliminary aerodynamic designs of aircraft, UAVs, hydrofoil and road vehicles.

Version 5.0 has the following features:
• Built-in automatic grid generation
• Built-in 3D compressible Euler Solver for fast aerodynamics analysis.
• Built-in 3D laminar Navier-Stokes solver
• Built-in 3D Reynolds Averaged Navier-Stokes (RANS) solver
• Multi-core flow solver processing on your Windows laptop or desktop using OpenMP
• Inputs STL files for processing
• Built-in wing/hydrofoil geometry creation tool
• Enables stability derivative computation using quasi-steady rigid body rotation
• Up to 100 actuator disc (RANS solver only) for simulating jets and prop wash
• Reports the lift, drag and moment coefficients
• Reports the lift, drag and moment magnitudes
• Plots surface pressure, velocity, Mach number and temperatures
• Produces 2-d plots of Cp and other quantities along constant coordinates line along the structure
The introductory price of Stallion 3D 5.0 is $3,495 for the yearly subscription or$8,000.  The software is also available in Lab and Class Packages.

► Airfoil Digitizer
18 Jun, 2017

Airfoil Digitizer is a software package for extracting airfoil data files from images. The software accepts images in the jpg, gif, bmp, png and tiff formats. Airfoil data can be exported as AutoCAD DXF files (line entities), UIUC airfoil database format and Hanley Innovations VisualFoil Format.

The following tutorial show how to use Airfoil Digitizer to obtain hard to find airfoil ordinates from pictures.

http:/www.hanleyinnovations.com/airfoildigitizerhelp.html

15 Feb, 2017

Have you ever wish for the power to solve your 3D aerodynamics analysis problems within your company just at the push of a button?  Stallion 3D gives you this very power using your MS Windows laptop or desktop computers. The software provides accurate CL, CD, & CM numbers directly from CAD geometries without the need for user-grid-generation and costly cloud computing.

Stallion 3D v 4 is the only MS windows software that enables you to solve turbulent compressible flows on your PC.  It utilizes the power that is hidden in your personal computer (64 bit & multi-cores technologies). The software simultaneously solves seven unsteady non-linear partial differential equations on your PC. Five of these equations (the Reynolds averaged Navier-Stokes, RANs) ensure conservation of mass, momentum and energy for a compressible fluid. Two additional equations captures the dynamics of a turbulent flow field.

Unlike other CFD software that require you to purchase a grid generation software (and spend days generating a grid), grid generation is automatic and is included within Stallion 3D.  Results are often obtained within a few hours after opening the software.

Do you need to analyze upwind and down wind sails?  Do you need data for wings and ship stabilizers at 10,  40, 80, 120 degrees angles and beyond? Do you need accurate lift, drag & temperature predictions at subsonic, transonic and supersonic flows? Stallion 3D can handle all flow speeds for any geometry all on your ordinary PC.

http://www.hanleyinnovations.com/stallion3d.html

Dr. Patrick Hanley is the owner of Hanley Innovations. He received his Ph.D. degree in fluid dynamics for Massachusetts Institute of Technology (MIT) department of Aeronautics and Astronautics (Course XVI). Dr. Hanley is the author of Stallion 3D, MultiSurface Aerodynamics, MultiElement Airfoils, VisualFoil and the booklet Aerodynamics in Plain English.

### CFD and others...top

► What Happens When You Run a LES on a RANS Mesh?
27 Dec, 2019

Surely, you will get garbage because there is no way your LES will have any chance of resolving the turbulent boundary layer. As a result, your skin friction will be way off. Therefore, your drag and lift will be a total disaster.

To actually demonstrate this point of view, we recently embarked upon a numerical experiment to run an implicit large eddy simulation (ILES) of the NASA CRM high-lift configuration from the 3rd AIAA High-Lift Prediction Workshop. The flow conditions are: Mach = 0.2, Reynolds number = 3.26 million based on the mean aerodynamic chord, and the angle of attack = 16 degrees.

A quadratic (Q2) mesh was generated by Dr. Steve Karman of Pointwise, and is shown in Figure 1.

Figure 1. Quadratic mesh for the NASA CRM high-lift configuration (generated by Pointwise)

The mesh has roughly 2.2 million mixed elements, and is highly clustered near the wall with an average equivalent y+ value smaller than one. A p-refinement study was conducted to assess the mesh sensitivity using our high-order LES tool based on the FR/CPR method, hpMusic. Simulations were performed with solution polynomial degrees of p = 1, 2 and 3, corresponding to 2nd, 3rd and 4th orders in accuracy respectively. No wall-model was used. Needless to say, the higher order simulations captured finer turbulence scales, as shown in Figure 2, which displays the iso-surfaces of the Q-criteria colored by the Mach number.

 p = 1

 p = 2

 p = 3
Figure 2. Iso-surfaces of the Q-criteria colored by the Mach number

Clearly the flow is mostly laminar on the pressure side, and transitional/turbulent on the suction side of the main wing and the flap. Although the p = 1 simulation captured the least scales, it still correctly identified the laminar and turbulent regions.

The drag and lift coefficients from the present p-refinement study are compared with experimental data from NASA in Table I. Although the 2nd order results (p = 1) are quite different than those of higher orders, the 3rd and 4th order results are very close, demonstrating very good p-convergence in both the lift and drag coefficients. The lift agrees better with experimental data than the drag, bearing in mind that the experiment has wind tunnel wall effects, and other small instruments which are not present in the computational model.

Table I. Comparison of lift and drag coefficients with experimental data
 CL CD p = 1 2.020 0.293 p = 2 2.411 0.282 p = 3 2.413 0.283 Experiment 2.479 0.252

This exercise seems to contradict the common sense logic stated in the beginning of this blog. So what happened? The answer is that in this high-lift configuration, the dominant force is due to pressure, rather than friction. In fact, 98.65% of the drag and 99.98% of the lift are due to the pressure force. For such flow problems, running a LES on a RANS mesh (with sufficient accuracy) may produce reasonable predictions in drag and lift. More studies are needed to draw any definite conclusion. We would like to hear from you if you have done something similar.

This study will be presented in the forthcoming AIAA SciTech conference, to be held on January 6th to 10th, 2020 in Orlando, Florida.

► Not All Numerical Methods are Born Equal for LES
15 Dec, 2018
Large eddy simulations (LES) are notoriously expensive for high Reynolds number problems because of the disparate length and time scales in the turbulent flow. Recent high-order CFD workshops have demonstrated the accuracy/efficiency advantage of high-order methods for LES.

The ideal numerical method for implicit LES (with no sub-grid scale models) should have very low dissipation AND dispersion errors over the resolvable range of wave numbers, but dissipative for non-resolvable high wave numbers. In this way, the simulation will resolve a wide turbulent spectrum, while damping out the non-resolvable small eddies to prevent energy pile-up, which can drive the simulation divergent.

We want to emphasize the equal importance of both numerical dissipation and dispersion, which can be generated from both the space and time discretizations. It is well-known that standard central finite difference (FD) schemes and energy-preserving schemes have no numerical dissipation in space. However, numerical dissipation can still be introduced by time integration, e.g., explicit Runge-Kutta schemes.

We recently analysed and compared several 6th-order spatial schemes for LES: the standard central FD, the upwind-biased FD, the filtered compact difference (FCD), and the discontinuous Galerkin (DG) schemes, with the same time integration approach (an Runge-Kutta scheme) and the same time step.  The FCD schemes have an 8th order filter with two different filtering coefficients, 0.49 (weak) and 0.40 (strong). We first show the results for the linear wave equation with 36 degrees-of-freedom (DOFs) in Figure 1.  The initial condition is a Gaussian-profile and a periodic boundary condition was used. The profile traversed the domain 200 times to highlight the difference.

Figure 1. Comparison of the Gaussian profiles for the DG, FD, and CD schemes

Note that the DG scheme gave the best performance, followed closely by the two FCD schemes, then the upwind-biased FD scheme, and finally the central FD scheme. The large dispersion error from the central FD scheme caused it to miss the peak, and also generate large errors elsewhere.

Finally simulation results with the viscous Burgers' equation are shown in Figure 2, which compares the energy spectrum computed with various schemes against that of the direct numerical simulation (DNS).

Figure 2. Comparison of the energy spectrum

Note again that the worst performance is delivered by the central FD scheme with a significant high-wave number energy pile-up. Although the FCD scheme with the weak filter resolved the widest spectrum, the pile-up at high-wave numbers may cause robustness issues. Therefore, the best performers are the DG scheme and the FCD scheme with the strong filter. It is obvious that the upwind-biased FD scheme out-performed the central FD scheme since it resolved the same range of wave numbers without the energy pile-up.

► Are High-Order CFD Solvers Ready for Industrial LES?
1 Jan, 2018
The potential of high-order methods (order > 2nd) is higher accuracy at lower cost than low order methods (1st or 2nd order). This potential has been conclusively demonstrated for benchmark scale-resolving simulations (such as large eddy simulation, or LES) by multiple international workshops on high-order CFD methods.

For industrial LES, in addition to accuracy and efficiency, there are several other important factors to consider:

• Ability to handle complex geometries, and ease of mesh generation
• Robustness for a wide variety of flow problems
• Scalability on supercomputers
For general-purpose industry applications, methods capable of handling unstructured meshes are preferred because of the ease in mesh generation, and load balancing on parallel architectures. DG and related methods such as SD and FR/CPR have received much attention because of their geometric flexibility and scalability. They have matured to become quite robust for a wide range of applications.

Our own research effort has led to the development of a high-order solver based on the FR/CPR method called hpMusic. We recently performed a benchmark LES comparison between hpMusic and a leading commercial solver, on the same family of hybrid meshes at a transonic condition with a Reynolds number more than 1M. The 3rd order hpMusic simulation has 9.6M degrees of freedom (DOFs), and costs about 1/3 the CPU time of the 2nd order simulation, which has 28.7M DOFs, using the commercial solver. Furthermore, the 3rd order simulation is much more accurate as shown in Figure 1. It is estimated that hpMusic would be an order magnitude faster to achieve a similar accuracy. This study will be presented at AIAA's SciTech 2018 conference next week.

(a) hpMusic 3rd Order, 9.6M DOFs
(b) Commercial Solver, 2nd Order, 28.7M DOFs
Figure 1. Comparison of Q-criterion and Schlieren

I certainly believe high-order solvers are ready for industrial LES. In fact, the commercial version of our high-order solver, hoMusic (pronounced hi-o-music), is announced by hoCFD LLC (disclaimer: I am the company founder). Give it a try for your problems, and you may be surprised. Academic and trial uses are completely free. Just visit hocfd.com to download the solver. A GUI has been developed to simplify problem setup. Your thoughts and comments are highly welcome.

Happy 2018!

► Sub-grid Scale (SGS) Stress Models in Large Eddy Simulation
17 Nov, 2017
The simulation of turbulent flow has been a considerable challenge for many decades. There are three main approaches to compute turbulence: 1) the Reynolds averaged Navier-Stokes (RANS) approach, in which all turbulence scales are modeled; 2) the Direct Numerical Simulations (DNS) approach, in which all scales are resolved; 3) the Large Eddy Simulation (LES) approach, in which large scales are computed, while the small scales are modeled. I really like the following picture comparing DNS, LES and RANS.

DNS (left), LES (middle) and RANS (right) predictions of a turbulent jet. - A. Maries, University of Pittsburgh

Although the RANS approach has achieved wide-spread success in engineering design, some applications call for LES, e.g., flow at high-angles of attack. The spatial filtering of a non-linear PDE results in a SGS term, which needs to be modeled based on the resolved field. The earliest SGS model was the Smagorinsky model, which relates the SGS stress with the rate-of-strain tensor. The purpose of the SGS model is to dissipate energy at a rate that is physically correct. Later an improved version called the dynamic Smagorinsky model was developed by Germano et al, and demonstrated much better results.

In CFD, physics and numerics are often intertwined very tightly, and one may draw erroneous conclusions if not careful. Personally, I believe the debate regarding SGS models can offer some valuable lessons regarding physics vs numerics.

It is well known that a central finite difference scheme does not contain numerical dissipation.  However, time integration can introduce dissipation. For example, a 2nd order central difference scheme is linearly stable with the SSP RK3 scheme (subject to a CFL condition), and does contain numerical dissipation. When this scheme is used to perform a LES, the simulation will blow up without a SGS model because of a lack of dissipation for eddies at high wave numbers. It is easy to conclude that the successful LES is because the SGS stress is properly modeled. A recent study with the Burger's equation strongly disputes this conclusion. It was shown that the SGS stress from the Smargorinsky model does not correlate well with the physical SGS stress. Therefore, the role of the SGS model, in the above scenario, was to stabilize the simulation by adding numerical dissipation.

For numerical methods which have natural dissipation at high-wave numbers, such as the DG, SD or FR/CPR methods, or methods with spatial filtering, the SGS model can damage the solution quality because this extra dissipation is not needed for stability. For such methods, there have been overwhelming evidence in the literature to support the use of implicit LES (ILES), where the SGS stress simply vanishes. In effect, the numerical dissipation in these methods serves as the SGS model. Personally, I would prefer to call such simulations coarse DNS, i.e., DNS on coarse meshes which do not resolve all scales.

I understand this topic may be controversial. Please do leave a comment if you agree or disagree. I want to emphasize that I support physics-based SGS models.
► 2016: What a Year!
3 Jan, 2017
2016 is undoubtedly the most extraordinary year for small-odds events. Take sports, for example:
• Leicester won the Premier League in England defying odds of 5000 to 1
• Cubs won World Series after 108 years waiting
In politics, I do not believe many people truly believed Britain would exit the EU, and Trump would become the next US president.

From a personal level, I also experienced an equally extraordinary event: the coup in Turkey.

The 9th International Conference on CFD (ICCFD9) took place on July 11-15, 2016 in the historic city of Istanbul. A terror attack on the Istanbul International airport occurred less than two weeks before ICCFD9 was to start. We were informed that ICCFD9 would still take place although many attendees cancelled their trips. We figured that two terror attacks at the same place within a month were quite unlikely, and decided to go to Istanbul to attend and support the conference.

Given the extraordinary circumstances, the conference organizers did a fine job in pulling the conference through. More than half of the attendees withdrew their papers. Backup papers were used to form two parallel sessions though three sessions were planned originally. We really enjoyed Istanbul with the beautiful natural attractions and friendly people.

Then on Friday evening, 12 hours before we were supposed to depart Istanbul, a military coup broke out. The government TV station was controlled by the rebels. However, the Turkish President managed to Facetime a private TV station, essentially turning around the event. Soon after, many people went to the bridge, the squares, and overpowered the rebels with bare fists.

A Tank outside my taxi

A beautiful night in Zurich

The trip back to the US was complicated by the fact that the FAA banned all direct flight from Turkey. I was lucky enough to find a new flight, with a stop in Zurich...

In 2016, I lost a very good friend, and CFD pioneer, Professor Jaw-Yen Yang. He suffered a horrific injury from tennis in early 2015. Many of his friends and colleagues gathered in Taipei on December 3-5 2016 to remember him.

This is a CFD blog after all, and so it is important to show at least one CFD picture. In a validation simulation [1] with our high-order solver, hpMusic, we achieved remarkable agreement with experimental heat transfer for a high-pressure turbine configuration. Here is a flow picture.

Computational Schlieren and iso-surfaces of Q-criterion

To close, I wish all of you a very happy 2017!

1. Laskowski GM, Kopriva J, Michelassi V, Shankaran S, Paliath U, Bhaskaran R, Wang Q, Talnikar C, Wang ZJ, Jia F. Future directions of high fidelity CFD for aerothermal turbomachinery research, analysis and design, AIAA-2016-3322.

20 Apr, 2016
The 64-bit version for the Linux operating system is now ready for you to download. Because of the complexities associated with various libraries, we experienced a delay of slightly more than a month. Here is the link again.

### ANSYS Blogtop

► How to Increase the Acceleration and Efficiency of Electric Cars for the Shell Eco Marathon
10 Oct, 2018

Illini EV Concept Team Photo at Shell Eco Marathon 2018

Weight is the enemy of all teams that design electric cars for the Shell Eco Marathon.

Reducing the weight of electric cars improves the vehicle’s acceleration and power efficiency. These performance improvements make all the difference come race day.

However, if the car’s weight is reduced too much, it could lead to safety concerns.

Illini EV Concept (Illini) is a Shell Eco Marathon team out of the University of Illinois. Team members use ANSYS academic research software to optimize the chassis of their electric car without compromising safety.

## Where to Start When Reducing the Weight of Electric Cars?

Front bump composite failure under a load of 2000N.

The first hurdle of the Shell Eco Marathon is an initial efficiency contest. Only the best teams from this efficiency assessment even make it into the race.

Therefore, Illini concentrates on reducing the most weight in the shortest amount of time to ensure it makes it to the starting line.

Illini notes that its focus is on reducing the weight of its electric car’s chassis.

“The chassis is by far the heaviest component of our car, so ANSYS was used extensively to help design our first carbon fiber monocoque chassis,” says Richard Mauge, body and chassis leader for Illini.

“Several loading conditions were tested to ensure the chassis was stiff enough and the carbon fiber did not fail using the composite failure tool,” he adds.

Competition regulations ensure the safety of all team members. These regulations state that each team must prove that their car is safe under various conditions. Simulation is a great tool to prove a design is within safety tolerances.

“One of these tests included ensuring the bulkhead could withstand a 700 N load in all directions, per competition regulations,” says Mauge. If the teams’ electric car designs can’t survive this simulation come race day, then their cars are not racing.

## Iterate and Optimize the Design of Electronic Cars with Simulation

Front bump deformation under a load of 2000N.

Simulations can do more than prove a design is safe. They can also help to optimize designs.

Illini uses what it learns from simulation to optimize the geometry of its electric car’s chassis.

The team found that its new designs have a torsional rigidity increase around 100 percent. This is after a 15 percent decrease in weight compared to last year’s model.

“Simulations ensure that the chassis is safe enough for our driver. It also proved that the chassis is lighter and stiffer than ever before. ANSYS composite analysis gave us the confidence to move forward with our radical chassis redesign,” notes Mauge.

The story optimization story continues from Illini. It plans to explore easier and more cost-effective ways to manufacture carbon fiber parts. For instance, the team wants to replace the core of its parts with foam and increase the number of bonded pieces.

If team members just go with their gut on these hunches, they could find themselves scratching their heads when something goes wrong. However, with simulations, the team makes better informed decisions about its redesigns and manufacturing process.

To get started with simulation, try our free student download. For student teams that need to solve in-depth problems, check out our software sponsorship program.

The post How to Increase the Acceleration and Efficiency of Electric Cars for the Shell Eco Marathon appeared first on ANSYS.

► Post-Processing Large Simulation Data Sets Quickly Over Multiple Servers
9 Oct, 2018

This engine intake simulation was post-processed using EnSight Enterprise. This allowed for the processing of a large data set to be shared among servers.

Simulation data sets have a funny habit of ballooning as engineers move through the development cycle. At some point, post-processing these data sets on a single machine becomes impractical.

Engineers can speed up post-processing by spatially or temporally decomposing large data sets so they can be post-processed across numerous servers.

The idea is to utilize the idle compute nodes you used to run the solver in parallel to now run the post-processing in parallel.

In ANSYS 19.2 Ensight Enterprise you can spatially or temporally decompose data sets. Ensignt Enterprise is an updated version of EnSight HPC.

## Post-Processing Using Spatial Decomposition

EnSight is a client/server architecture. The client program takes care of the graphical user interface (GUI) and rendering operations, while the server program loads the data, creates parts, extracts features and calculates results.

If your model is too large to post-process on a single machine, you can utilize the spatial decomposed parallel operation to assign each spatial partition to its own EnSight Server. A good server-to-model ratio is one server for every 50 million elements.

Each EnSight Server can be located on a separate compute node on any compute resource you’d like. This allows engineers to utilize the memory and processing power of heterogeneous high-performance computing (HPC) resources for data set post-processing.

The engineers effectively split the large data set up into pieces with each piece assigned to its own compute resource. This dramatically increases the data set sizes you can load and process.

Once you have loaded the model into EnSight Enterprise, there are no additional changes to your workflow, experience or operations.

## Post-Processing Using Temporal Decomposition

Keep in mind that this decomposition concept can also be applied to transient data sets. In this case, the dataset is split up temporally rather than spatially. In this scenario, each server receives its own set of time steps.

A turbulence simulation created using EnSight Enterprise post-processing

EnSight Enterprise offers performance gains when the server operations outweigh the communication and rendering time of each time step. Since it’s hard to predict network communication or rendering workloads, you can’t easily create a guiding principle for the server-to-model ratio.

However, you might want to use a few servers when your model has more than 10 million elements and over a hundred time steps. This will help keep the processing load of each server to a moderate level.

## How EnSight Speeds Up the Post-Processing of Large Simulation Data Sets

Another good tip to ensure you are post-processed optimally within EnSight Enterprise. Engineers achieve the best performance gains by pre-decomposing the data and locating it locally to the compute resources they anticipate using. Ideally, this data should be in EnSight Case format.

To learn more, check out Ensight or register for the webinar Analyze, Visualize and Communicate Your Simulation Data with ANSYS EnSight.

The post Post-Processing Large Simulation Data Sets Quickly Over Multiple Servers appeared first on ANSYS.

► Discovery AIM Offers Design Teams Rapid Results and Physics-Aware Meshing
8 Oct, 2018

Your design team will make informed decisions about the products they create when they bring detailed simulations up front in the development cycle.

The 19.2 release of ANSYS Discovery AIM facilitates the need of early simulations.

It does this by streamlining templates for physics-aware meshing and rapid results.

## High-Fidelity Simulation Through Physics-Aware Meshing

Discovery AIM user interface with a solution fidelity slide bar (top left), area of interest marking tool (left, middle), manual mesh controls (bottom, center) and a switch to turn the mesh display on and off (right, top).

Analysts have likely told your design team about the importance of a quality mesh to achieve accurate simulation results.

Creating high quality meshes takes time and specialized training. Your design team doesn’t likely have the time or patience to learn this art.

To account for this, Discovery AIM automatically incorporates physics-aware meshing behind the scenes. In fact, your design team doesn’t even need to see the mesh creation process to complete the simulation.

This workflow employs several meshing best practices analysts typically use. The tool even accounts for areas that require mesh refinements based on the physics being assessed.

For instance, areas with a sliding contact gain a finer mesh so the sliding behavior can be accurately simulated. Additionally, areas near the walls of fluid-solid interfaces are also refined to ensure this interaction is properly captured. Physics-aware meshing ensures small features and areas of interests won’t get lost in your design team’s simulation.

The simplified meshing workflow also lets your design team choose their desired solution fidelity. This input will help the software balance the time the solver takes to compute results with the accuracy of the results.

Though physics-aware meshing can create the mesh under the hood of the simulation process, it still has tools allowing user-control of the mesh. This way, if your design team chooses to dig into the meshing details — or an analyst decides to step in — they can finely tune the mesh.

Capabilities like this further empower designers as techniques and knowledge traditionally known only by analysts are automated in an easy-to-use fashion.

## Gain Rapid Results in Important Areas You Might Miss

The 19.2 release of Discovery AIM has seen improvements with its ability to enable your design team to explore simulation results.

Many analysts will know instinctively where to focus their post-processing, but without this experience, designers may miss areas of interest.

Discovery AIM enables the designer to interactively explore and identify these critical results. These initial results are rapidly displayed as contours, streamlines or field flow lines.

Field flow and streamlines for an electromagnetics simulation

Once your design team finds locations of interest within the results, they can create higher fidelity results to examine those area of interest in further detail. Designers can then save the results and revisit them when comparing design points or after changing simulation inputs.

To learn more about other changes to Discovery AIM — like the ability to directly access fluid results — watch the Discovery AIM 19.2 release recorded webinar or take it for a test drive.

The post Discovery AIM Offers Design Teams Rapid Results and Physics-Aware Meshing appeared first on ANSYS.

► Simulation Optimizes a Chemotherapy Implant to Treat Pancreatic Cancer
5 Oct, 2018

Traditional chemotherapy can often be blocked by a tumor’s stroma.

There are few illnesses as crafty as pancreatic cancer. It spreads like weeds and resists chemotherapy.

Pancreatic cancer is often asymptomatic, has a low survival rate and is often misdiagnosed as diabetes. And, this violent killer is almost always inoperable.

The pancreatic tumor’s resistance to chemotherapy comes from a shield of supporting connective tissue, or stroma, which it builds around itself.

Current treatments attempt to overcome this defense by increasing the dosage of intravenously administered chemotherapy. Sadly, this rarely works, and the high dosage is exceptionally hard on patients.

Nonetheless, doctors need a way to shrink these tumors so that they can surgically remove them without risking the numerous organs and vasculature around the pancreas.

“We say if you can’t get the drugs to the tumor from the blood, why not get it through the stroma directly?” asks William Daunch, CTO at Advanced Chemotherapy Technologies (ACT), an ANSYS Startup Program member. “We are developing a medical device that implants directly onto the pancreas. It passes drugs through the organ, across the stroma to the tumor using iontophoresis.”

By treating the tumor directly, doctors can theoretically shrink the tumor to an operable size with a smaller dose of chemotherapy. This should significantly reduce the effects of the drugs on the rest of the patient’s body.

## How to Treat Pancreatic Cancer with a Little Electrochemistry

Simplified diagram of the iontophoresis used by ACT’s chemotherapy medical device.

Most of the drugs used to treat pancreatic cancer are charged. This means they are affected by electromotive forces.

ACT has created a medical device that takes advantage of the medication’s charge to beat the stroma’s defenses using electrochemistry and iontophoresis.

The device contains a reservoir with an electrode. The reservoir connects to tubes that connect to an infusion pump. This setup ensures that the reservoir is continuously filled. If the reservoir is full, the dosage doesn’t change.

The tubes and wires are all connected into a port that is surgically implanted into the patient’s abdomen.

A diagram of ACT’s chemotherapy medical device.

The circuit is completed by a metal panel on the back of the patient.

“When the infusion pump runs, and electricity is applied, the electromotive forces push the medication into the stroma’s tissue without a needle. The medication can pass up to 10 to 15 mm into the stroma’s tissue in about an hour. This is enough to get through the stroma and into the tumor,” says Daunch.

“Lab tests show that the medical device was highly effective in treating human pancreatic cancer cells within mice,” added Daunch. “With conventional infusion therapy, the tumors grew 700 percent and with the device working on natural diffusion alone the tumors grew 200 percent. However, when running the device with iontophoresis, the tumor shrank 40 percent. This could turn an inoperable tumor into an operable one.” Subsequent testing of a scaled-up device in canines demonstrated depth of penetration and the low systemic toxicity required for a human device.

Daunch notes that the Food and Drug Administration (FDA) took notice of these results. ACT’s next steps are to develop a human clinical device and move onto to human safety trials.

## Simulation Optimized the Fluid Dynamics in the Pancreatic Cancer Chemotherapy Implant

Before these promising tests, ACT faced a few design challenges when coming up with their chemotherapy implant.

For example, “There was some electrolysis on the electrode in the reservoir. This created bubbles that would change the electrode’s impedance,” explains Daunch. “We needed a mechanism to sweep the bubbles from the surface.”

An added challenge is that ACT never knows exactly where doctors will place the device on the pancreas. As a result, the mechanism to sweep the bubbles needs to work from any orientation.

Simulations help ACT design their medical device so bubbles do not collect on the electrode.

“We used ANSYS Fluent and ANSYS Discovery Live to iterate a series of designs,” says Daunch. “Our design team modeled and validated our work very quickly. We also noticed that the bubbles didn’t need to leave the reservoir, just the electrode.”

“If we place the electrode on a protrusion in a bowl-shaped reservoir the bubbles move aside into a trough,” explains Daunch. “The fast fluid flow in the center of the electrode and the slower flow around it would push the bubbles off the electrode and keep them off until the bubbles floated to the top.”

As a result, the natural fluid flow within the redesigned reservoir was able to ensure the bubbles didn’t affect the electrode’s impedance.

To learn how your startup can use computational fluid dynamics (CFD) software to address your design challenges, please visit the ANSYS Startup Program.

The post Simulation Optimizes a Chemotherapy Implant to Treat Pancreatic Cancer appeared first on ANSYS.

► Making Wireless Multigigabit Data Transfer Reliable with Simulation
4 Oct, 2018

The demand for wireless communications with high data transfer rates is growing.

Consumers want wireless 4K video streams, virtual reality, cloud backups and docking. However, it’s a challenge to offer these data transfer hogs wirelessly.

Peraso aims to overcome this challenge with their W120 WiGig chipset. This device offers multigigabit data transfers, is as small as a thumb-drive and plugs into a USB 3.0 port.

The chipset uses the Wi-Fi Alliance’s new wireless networking standard, WiGig.

This standard adds a 60 GHz communication band to the 2.4 and 5 GHz bands used by traditional Wi-Fi. The result is higher data rates, lower latency and dynamic session transferring with multiband devices.

In theory, the W120 WiGig chipset could run some of the heaviest data transfer hogs on the market without a cord. Peraso’s challenge is to design a way for the chipset to dissipate all the heat it generates.

Peraso uses the multiphysics capabilities within the ANSYS Electronics portfolio to predict the Joule heating and the subsequent heat flow effects of the W120 WiGig chipset. This information helps them iterate their designs to better dissipate the heat.

## How to Design High Speed Wireless Chips That Don’t Overheat

Systems designers know that asking for high-power transmitters in a compact and cost-effective enclosure translates into a thermal challenge. The W120 WiGig chipset is no different.

A cross section temperature map of the W120 WiGig chipset’s PCB. The map shows hot spots where air flow is constrained by narrow gaps between the PCB and enclosure.

The chipset includes active/passive components and two main chips that are mounted on a printed circuit board (PCB). The system reaches considerably high temperatures due to the Joule heating effect.

To dissipate this heat, design engineers include a large heat sink that connects only to the chips and a smaller one that connects only to the PCB. The system is also enclosed in a casing with limited openings.

Simulation of the air flow around the W120 WiGig chipset without an enclosure. Simulation was made using ANSYS Icepak.

Traditionally, optimizing this set up takes a lot of trial and error as measuring the air flow within the enclosure would be challenging.

Instead, Peraso uses ANSYS SIwave to simulate the Joule heating effects of the system. This heat map is transferred to ANSYS Icepak, which then simulates the current heat flow, orthotropic thermal conductivity, heat sources and other thermal effects.

This multiphysics simulation enables Peraso to predict the heat distribution and the temperature at every point of the W120 WiGig chipset.

From there, Peraso engineers iterate their designs until they reached their coolest setup.

This simulation led design tactic helps Peraso optimize their system until they reached a heat transfer balance they need. To learn how Peraso performed this iteration, read Cutting the Cords.

The post Making Wireless Multigigabit Data Transfer Reliable with Simulation appeared first on ANSYS.

► Designing 5G Cellular Base Station Antennas Using Parametric Studies
3 Oct, 2018

There is only so much communication bandwidth available. This will make it difficult to handle the boost in cellular traffic expected from the 5G network using conventional cellular technologies.

In fact, cellular networks are already running out of bandwidth. This severely limits the number of users and data rates that can be accommodated by wireless systems.

One potential solution is to leverage beamforming antennas. These devices transmit different signals to different locations on the cellular network simultaneously over the same frequency.

Pivotal Commware is using ANSYS HFSS to design beamforming antennas for cellular base stations that are much more affordable than current technology.

## How 5G Networks Will Send More Signals on Existing Bandwidths

A 28 GHz antenna for a cellular base station.

Traditionally, cellular technologies — 3G and 4G LTE — crammed more signals on the existing bandwidth by dividing the frequencies into small segments and splitting the signal time into smaller pulses.

The problem is, there is only so much you can do to chop up the bandwidth into segments.

Alternatively, Pivotal’s holographic beamforming (HBF) antennas are highly directional. This means they can split up the physical space a signal moves through.

This way, two cells in two locations can use the same frequency at the same time without interfering with each other.

Additionally, these HBF antennas use varactor (variable capacitors) and electronic components that are simpler and more affordable than existing beamforming antennas.

## How to Design HBF Antennas for 5G Cellular Base Stations

A parametric study of Pivotal’s HBF designs allowed them to look at a large portion of their design space and optimize for C-SWaP and roll-off. This study looks at roll-off as a function of degrees from the centerline of the antenna.

Antenna design companies — like Pivotal — are always looking to design devices that optimize cost, size, weight and power (C-SWaP) and performance.

So, how was Pivotal able to account for C-SWaP and performance so thoroughly?

Traditionally, this was done by building prototypes, finding flaws, creating new designs and integrating manually.

Meeting a product launch with an optimized product using this manual method is grueling.

Pivotal instead uses ANSYS HFSS to simulate their 5G antennas digitally. This allows them to assess their HBF antennas and iterate their designs faster using parametric studies.

For instance, Pivotal wants to optimize their design for performance characteristics like roll-off. To do so they can plug in the parameter values, run simulations with these values and see how each parameter affects roll-off.

By setting up parametric studies, Pivotal assess which parameters affect performance and C-SWaP the most. From there they could weigh different trade-offs until they settled on an optimized design that accounted for all the factors they studied.

To see how Pivotal set up their parametric studies and optimize their antenna designs, read 5G Antenna Technology for Smart Products.

The post Designing 5G Cellular Base Station Antennas Using Parametric Studies appeared first on ANSYS.

### Convergent Science Blogtop

► The Search for Soot-free Diesel: Modeling Ducted Fuel Injection With CONVERGE
26 Mar, 2020

At the upcoming CONVERGE User Conference, which will be held online from March 31–April 1, Andrea Piano will present results from experimental and numerical studies of the effects of ducted fuel injection on fuel spray characteristics. Dr. Piano is a Research Assistant in the e3 group, coordinated by Prof. Federico Millo at Politecnico di Torino, and these are the first results to be reported from their ongoing collaboration with Prof. Lucio Postrioti at Università degli Studi di Perugia, Andrea Bianco at Powertech Engineering, and Francesco Pesce and Alberto Vassallo at General Motors Global Propulsion Systems. This work is a great example of how CONVERGE can be used in tandem with experimental methods to advance research at the cutting edge of engine technology. Keep reading for a preview of the results that Dr. Piano will discuss in greater detail in his online presentation.

The idea behind ducted fuel injection (DFI), originally conceived by Charles Mueller at Sandia National Laboratories, is to suppress soot formation in diesel engines by allowing the fuel to mix more thoroughly with air before it ignites1. Soot forms when a fuel doesn’t burn completely, which happens when the fuel-to-air ratio is too high. In DFI, a small tube, or duct, is placed near the nozzle of the fuel injector and directed along the axis of the fuel stream toward the autoignition zone. The fuel spray that travels through this duct is better mixed than it would be in a ductless configuration. Experiments at Sandia have shown that DFI can reduce soot formation by as much as 95%, demonstrating the enormous potential of this technology for curtailing harmful emissions from diesel engines.

While the Sandia researchers have focused on heavy-duty diesel applications, Dr. Piano and his collaborators are targeting smaller engines, such as those found in passenger cars and light-duty trucks. To understand how the fuel spray evolves in the presence of a duct, they first performed imaging and phase Doppler anemometry analyses of non-reacting sprays in a constant-volume test vessel. Figure 1 shows a sample of the experimental results. The video on the left corresponds to a free spray configuration with no duct, while the video on the right corresponds to a ducted configuration. Observe how the dark liquid breaks up and evaporates more quickly in the ducted configuration—this is the enhanced mixing that occurs in DFI.

Their next step was to develop a CFD model of the fuel spray that could be calibrated against the experimental results. Dr. Piano and his colleagues reproduced the geometry of the experimental setup in a CONVERGE environment, using physical models available in CONVERGE to simulate the processes of spray breakup, evaporation, and boiling, as well as the interactions between the spray and the duct. With fixed embedding and Adaptive Mesh Refinement, they were able to increase the grid resolution in the vicinity of the spray and the duct without a significant increase in computational cost. They simulated the spray penetration for both the free spray and the ducted configuration over a range of operating conditions and validated those results against the experimental data.

With a calibrated spray model in hand, the researchers were then able to run predictive simulations of DFI for reacting fuel sprays. They combined their spray model with the SAGE detailed chemical kinetics solver for combustion modeling, along with the Particulate Mimic model of soot formation. They ran simulations at different rail pressures and vessel temperatures to see how DFI would affect the amount of soot mass produced under engine-like operating conditions. Figures 2 and 3 show examples of the simulation results for a rail pressure of 1200 bar and a vessel temperature of 1000 K. Consistent with the findings of Mueller et al.1, these results show a dramatic reduction in the mass of soot produced during combustion in the ducted configuration as compared to the free spray configuration.

While these early results are promising, Dr. Piano and his collaborators are just getting started. They will continue using CONVERGE to investigate phenomena such as the duct thermal behavior and to explore the effects of different geometries and operating conditions, with the long-term goal of incorporating DFI into the design of a real engine. If you are interested in learning more about this work, be sure to sign up for the CONVERGE User Conference today!

References

[1] Mueller, C.J., Nilsen, C.W., Ruth, D.J., Gehmlich, R.K., Pickett, L.M., and Skeen, S.A., “Ducted fuel injection: A new approach for lowering soot emissions from direct-injection engines,” Applied Energy, 204, 206-220, 2017. DOI: 10.1016/j.apenergy.2017.07.001

► An Evening With the Experts: Scaling CFD With High-Performance Computing
25 Feb, 2020
Listen to the full audio of the panel discussion.

As computing technology continues to advance rapidly, running simulations on hundreds and even thousands of cores is becoming standard practice in the CFD industry. Likewise, CFD software is continually evolving to keep pace with the advances in hardware. For example, CONVERGE 3.0, the latest major release of our software, is specifically designed to scale well in parallel on modern high-performance computing (HPC) systems. It’s clear that HPC is the future of CFD, so how does this shift affect those of us running simulations and how can we make the most of the increased availability of computational resources? At the 2019 CONVERGE User Conference–North America, we assembled a panel of engineers from industry and government to share their expertise.

In the panel discussion, which you can listen to above, you’ll learn about the computing resources available on the cloud and at the U.S. national laboratories and how to take advantage of them. The panelists discuss the types of novel, one-of-a-kind studies that HPC enables and how to handle post-processing data from massive cases run across many cores. Additionally, you’ll get a look at where post-processing is headed in the future to manage the ever-increasing amounts of data generated form large-scale simulations. Listen to the full panel discussion above!

### Panelists

Alan Klug, Vice President of Customer Development, Tecplot

Sibendu Som, Manager of the Computational Multi-Physics Section, Argonne National Laboratory

Joris Poort, CEO and Founder, Rescale

Kelly Senecal, Co-Founder and Owner, Convergent Science

### Moderator

Tiffany Cook, Partner & Public Relations Manager, Convergent Science

19 Dec, 2019

2019 proved to be an exciting and eventful year for Convergent Science. We released the highly anticipated major rewrite of our software, CONVERGE 3.0. Our United States, European, and Indian offices all saw significant increases in employee count. We have also continued to forge ahead in new application areas, strengthening our presence in the pump, compressor, biomedical, aerospace, and aftertreatment markets, and breaking into the oil and gas industry. Of course, we remain dedicated to simulating internal combustion engines and developing new tools and resources for the automotive community. In particular, we are expanding our repertoire to encompass batteries and electric motors in addition to conventional engines. Our team at Convergent Science continues to be enthusiastic about advancing simulation capabilities and providing unmatched customer support to empower our users to tackle hard CFD problems.

### CONVERGE 3.0

As I mentioned above, this year we released a major new version of our software, CONVERGE 3.0. We have frequently discussed 3.0 in the past few months, including in my recent blog post, so I’ll keep this brief. We set out to make our code more flexible, enable massive parallel scaling, and expand CONVERGE’s capabilities. The results have been remarkable. CONVERGE 3.0 scales with near-ideal efficiencies on thousands of cores, and the addition of inlaid meshes, new physical models, and enhanced chemistry capabilities have opened the door to new applications. Our team invested a lot of effort into making 3.0 a reality, and we’re very proud of what we’ve accomplished. Of course, now that CONVERGE 3.0 has been released, we can all start eagerly anticipating our next major release, CONVERGE 3.1.

### Computational Chemistry Consortium

2019 was a big year for the Computational Chemistry Consortium (C3). In July, the first annual face-to-face meeting took place at the Convergent Science World Headquarters in Madison, Wisconsin. Members of industry and researchers from the National University of Ireland Galway, Lawrence Livermore National Laboratory, RWTH Aachen University, and Politecnico di Milano came together to discuss the work done during the first year of the consortium and establish future research paths. The consortium is working on the C3 mechanism, a gasoline and diesel surrogate mechanism that includes NOx and PAH chemistry to model emissions. The first version of the mechanism was released this fall for use by C3 members, and the mechanism will be refined over the coming years. Our goal is to create the most accurate and consistent reaction mechanism for automotive fuels. Stay tuned for future updates!

### Third Annual European User Conference

Barcelona played host to this year’s European CONVERGE User Conference. CONVERGE users from across Europe gathered to share their recent work in CFD on topics including turbulent jet ignition, machine learning for design optimization, urea thermolysis, ammonia combustion in SI engines, and gas turbines. The conference also featured some exciting networking events—we spent an evening at the beautiful and historic Poble Espanyol and organized a kart race that pitted attendees against each other in a friendly competition.

### Inaugural CONVERGE User Conference–India

This year we hosted our first-ever CONVERGE User Conference–India in Bangalore and Pune. The conference consisted of two events, each covering different application areas. The event in Bangalore focused on applications such as gas turbines, fluid-structure interaction, and rotating machinery. In Pune, the emphasis was on IC engines and aftertreatment modeling. We saw presentations from both companies and universities, including General Electric, Cummins, Caterpillar, and the Indian Institutes of Technology Bombay, Kanpur, and Madras. We had a great turnout for the conference, with more than 200 attendees across the two events.

### CONVERGE in the Big Easy

The sixth annual CONVERGE User Conference–North America took place in New Orleans, Louisiana. Attendees came from industry, academic institutions, and national laboratories in the U.S. and around the globe. The technical presentations covered a wide variety of topics, including flame spray pyrolysis, rotating detonation engines, machine learning, pre-chamber ignition, blood pumps, and aerodynamic characterization of unmanned aerial systems. This year, we hosted a panel of CFD and HPC experts to discuss scaling CFD across thousands of processors; how to take advantage of clusters, supercomputers, and the cloud to run large-scale simulations; and how to post-process large datasets. For networking events, we took a dinner cruise down the Mississippi River and encouraged our guests to explore the vibrant city of New Orleans.

### KAUST Workshop

In 2019, we hosted the First CONVERGE Training Workshop and User Meeting at the King Abdullah University of Science and Technology (KAUST) in Saudi Arabia. Attendees came from KAUST and other Saudi Arabian universities and companies for two days of keynote presentations, hands-on CONVERGE tutorials, and networking opportunities. The workshop focused on leveraging CONVERGE for a variety of engineering applications, and running CONVERGE on local workstations, clusters, and Shaheen II, a world-class supercomputer located at KAUST.

### Best Use of HPC in Automotive

We and our colleagues at Argonne National Laboratory and Aramco Research Center – Detroit received this year’s 2019 HPCwire Editors’ Choice Award in the category of Best Use of HPC in Automotive. We were incredibly honored to receive this award for our work using HPC and AI to quickly optimize the design of a clean, highly efficient gasoline compression ignition engine. Using CONVERGE, we tested thousands of engine design variations in parallel to improve fuel efficiency and reduce emissions. We ran the simulations in days, rather than months, on an IBM Blue Gene/Q supercomputer located at Argonne National Laboratory and employed machine learning to further reduce design time. After running the simulations, the best-performing engine design was built in the real world. The engine demonstrated a reduction in CO2 of up to 5%. Our work shows that pairing HPC and AI to rapidly optimize engine design has the potential to significantly advance clean technology for heavy-duty transportation.

### Convergent Science Around the Globe

2019 was a great year for CONVERGE and Convergent Science around the world. In the United States, we gained nearly 20 employees. We added a new Convergent Science office in Houston, Texas, to serve the oil and gas industry. In addition, we have continued to increase our market share in other areas, including automotive, gas turbine, and pumps and compressors.

In Europe, we had a record year for new license sales, up 70% from 2018. A number of new employees joined our European team, including new engineers, sales personnel, and office administrators. We attended and exhibited at tradeshows on a breadth of topics all over Europe, and we expanded our industry and university clientele.

Our Indian office celebrated its second anniversary in 2019. The employee count nearly doubled in size from 2018, with the addition of several new software developers and marketing and support engineers. The first Indian CONVERGE User Conference was a huge success–we had to increase the maximum number of registrants to accommodate everyone who wanted to attend. We have also grown our client base in the transportation sector, bringing new customers in the automotive industry on board.

In Asia, our partners at IDAJ continue to do a fantastic job supporting CONVERGE. CONVERGE sales significantly increased in 2019 compared to 2018. And at this year’s IDAJ CAE Solution Conference, speakers from major corporations presented CONVERGE results, including Toyota, Daihatsu, Mazda, and DENSO.

While we like to recognize the successes of the past year, we’re always looking toward the future. Computing technology is constantly evolving, and we are eager to keep advancing CONVERGE to make the most of the increased availability of computational resources. With the expanded functionality that CONVERGE 3.0 offers, we’re also looking forward to delving into untapped application areas and breaking into new markets. In the upcoming year, we are excited to form new collaborations and strengthen existing partnerships to promote innovation and keep CONVERGE on the cutting-edge of CFD software.

► CONVERGE 3.0: From Specialized Software to CFD Powerhouse
25 Nov, 2019

When Eric, Keith, and I first wrote CONVERGE back in 2001, we wrote it as a serial code. That probably sounds a little crazy, since practically all CFD simulations these days are run in parallel on multiple CPUs, but that’s how it started. We ended up taking our serial code and making it parallel, which is arguably not the best way to create a parallel code. As a side effect of writing the code this way, there were inherent parts of CONVERGE that did not scale well, both in terms of speed and memory. This wasn’t a real issue for our clients who were running engine simulations on relatively small numbers of cores. But as time wore on, our users started simulating many different applications beyond IC engines, and those simulating engines wanted to run finer meshes on more cores. At the same time, computing technology was evolving from systems with relatively few cores per node and relatively high memory per core to modern HPC clusters with more cores and nodes per system and relatively less memory per core. We knew at some point we would have to rewrite CONVERGE to take advantage of the advancements in computing technology.

We first conceived of CONVERGE 3.0 around five years ago. At that point, none of the limitations in the code were significantly affecting our clients, but we would get the occasional request that was simply not feasible in the current software. When we got those requests, we would categorize them as “3.0”—requests we deemed important, but would have to wait until we rewrote the code. After a few years, some of the constraints of the code started to become real limitations for our clients, so our developers got to work in earnest on CONVERGE 3.0. Much of the core framework and infrastructure was redesigned from the ground up in version 3.0, including a new mesh API, surface and grid manipulation tools, input and output file formats, and load balancing algorithms. The resulting code enables our users to run larger, faster, and more accurate simulations for a wider range of applications.

Scalability and Shared Memory

Two of our major goals in rewriting CONVERGE were to improve the scalability of the code and to reduce the memory requirements. Scaling in CONVERGE 2.x versions was limited in large part because of the parallelization method. In the 2.x versions, the simulation domain is partitioned using blocks coarser than the solution grid. This can cause a poor distribution of workload among processors if you have high levels of embedding or Adaptive Mesh Refinement (AMR). In 3.0, the solution grid is now partitioned directly, so you can achieve a good load balance even with very high levels of embedding and AMR. In addition, load balancing is now performed automatically instead of on a fixed schedule, so the case is well balanced throughout more of the run. With these changes, we’ve seen a dramatic improvement in scaling in 3.0, even on thousands of cores.

To reduce memory requirements, our developers moved to a shared memory strategy and removed redundancies that existed in previous versions of CONVERGE. For example, many data structures, like surface triangulation, that were stored once per core in the 2.x versions are now only stored once per compute node. Similarly, CONVERGE 3.0 no longer stores the entire grid connectivity on every core as was done in previous versions. The memory footprint in 3.0 is thus greatly reduced, and memory requirements also scale well into thousands of cores.

Inlaid Mesh

Apart from the codebase rewrite, another significant change we made was to incorporate inlaid meshes into CONVERGE. For years, users have been asking for the ability to add extrusion layers to boundaries, and it made sense to add this feature now. As many of you are probably aware, autonomous meshing is one of the hallmarks of our software. CONVERGE automatically generates an optimized Cartesian mesh at runtime and dynamically refines the mesh throughout the simulation using AMR. All of this remains the same in CONVERGE 3.0, and you can still use meshes exactly as they were in all previous versions of CONVERGE! However now we’ve added the option to create an inlaid mesh made up of cells of arbitrary shape, size, and orientation. The inlaid mesh can be extruded from a triangulated surface (e.g., a boundary layer) or it can be a shaped mesh away from a surface (e.g., a spray cone). For the remainder of the domain not covered by an inlaid mesh, CONVERGE uses our traditional Cartesian mesh technology.

Inlaid meshes are always optional, but in some cases they can provide accurate results with fewer cells compared to a traditional Cartesian mesh. In the example of a boundary layer, you can now refine the mesh in only the direction normal to the surface, instead of all three directions. You can also align an inlaid mesh with the direction of the flow, which wasn’t always possible when using a Cartesian mesh. This feature makes CONVERGE better suited for certain applications, like external aerodynamics, than it was previously.

Combustion and Chemistry

In CONVERGE 3.0, our developers have also enhanced and added to our combustion models and chemistry tools. For the SAGE detailed chemistry solver, we optimized the rate calculations, improved the procedure to assemble the sparse Jacobian matrix, and we introduced a new preconditioner. The result is significant speedup in the chemistry solver, especially for large reaction mechanisms (>150 species). If you thought our chemistry solver was fast before (and it was!), you will be amazed at the speed of the new version. In addition, 3.0 features two new combustion models. In most large eddy simulations (LES) of premixed flames, the cells are not fine enough to resolve the laminar flame thickness. The thickened flame model for LES allows you to increase the flame thickness without changing the laminar flamespeed. The second new model, the SAGE three-point PDF model, can be used to account for turbulence-chemistry interaction (more specifically, the commutation error) when modeling turbulent combusting flows with RANS.

On the chemistry tools side, we’ve added a number of new 0D chemical reactors, including variable volume with heat loss, well-stirred, plug flow, and 0D engine. The 1D laminar flamespeed solver has seen significant improvements in scalability and parallelization, and we have new table generation tools in CONVERGE Studio for tabulated kinetics of ignition (TKI), tabulated laminar flamespeed (TLF), and flamelet generated manifold (FGM).

To streamline our users’ workflow, we have implemented several updates in CONVERGE Studio, CONVERGE’s graphical user interface (GUI). We partnered with Spatial to allow users to directly import CAD files into CONVERGE Studio 3.0, and triangulate the geometry on the fly in a way that’s optimized for CONVERGE. Additionally, Tecplot for CONVERGE, CONVERGE’s post-processing and visualization software, can now read CONVERGE output files directly, for a smoother workflow from start to finish.

CONVERGE 3.0 was a long time in the making, and we’re very excited about the new capabilities and opportunities this version offers our users. 3.0 is a big step towards CONVERGE being a flexible toolbox for solving any CFD problem.

[1] The National Center for Supercomputing Applications (NCSA) at the University of Illinois at Urbana-Champaign provides supercomputing and advanced digital resources for the nation’s science enterprise. At NCSA, University of Illinois faculty, staff, students, and collaborators from around the globe use advanced digital resources to address research grand challenges for the benefit of science and society. The NCSA Industry Program is the largest Industrial HPC outreach in the world, and it has been advancing one third of the Fortune 50® for more than 30 years by bringing industry, researchers, and students together to solve grand computational problems at rapid speed and scale. The CONVERGE simulations were run on NCSA’s Blue Waters supercomputer, which is one of the fastest supercomputers on a university campus. Blue Waters is supported by the National Science Foundation through awards ACI-0725070 and ACI-1238993.

► Changing the CFD Conference Game with the CONVERGE UC
21 Aug, 2019

As the 2019 CONVERGE User Conference in New Orleans approaches, I’ve been thinking about the past five years of CONVERGE events. Let me take you back to the first CONVERGE User Conference. It was September 2014 in Madison, Wisconsin, and I was one of the first speakers. I talked about two-phase flows and the spray modeling we were doing at Argonne National Laboratory. Many of the people in the audience didn’t know you could do the kinds of calculations in CONVERGE that we were doing. Take needle wobble, for example. At the time, people didn’t know that you could not only move the needle up and down, but you could actually simulate it wobbling. After my talk, we had many interesting discussions with the other attendees. We made connections with international companies that we otherwise would not have had the chance to meet, and we formed collaborations with some of those companies that are still ongoing today.

At Argonne National Laboratory, I lead a team of more than 20 researchers, all of them focused on simulating either piston engines or gas turbines using high-performance computing. Our goal is to improve the predictive capability of piston engine and gas turbine simulations, and we do a lot of our work using CONVERGE. We develop physics-based models that we couple with CONVERGE to gain deeper insights from our simulations.

We routinely attend and present our work at conferences like SAE World Congress and ASME, and what really sets the CONVERGE User Conference apart is the focus of the event—it’s dedicated towards the people doing simulation work with piston engines, gas turbines, and other real-world applications. The user conference is the go-to place where we can meet all of the people doing 3D CFD simulations, so it’s a fantastic networking opportunity. We get to speak to people from academia and industry and learn about their research needs—understand what their pain points are, what their bottlenecks are, where the physics is not predictive enough. Then we take that information back to Argonne, and it helps us focus our research.

Apart from the networking, the CONVERGE User Conference is also a great venue for presenting. My team has presented at the CONVERGE conferences on a wide variety of topics, including lean blow-out in gas turbine combustors, advanced ignition systems, co-optimization of engines and fuels, predicting cycle-to-cycle variation, machine learning for design optimizations, and modeling turbulent combustion in compression ignition and spark ignition engines. The attendees are engaged and highly technical, so you get direct, focused feedback on your work that can help you find solutions to challenges you may be encountering or give ideas for future studies.

The presenters themselves take the conference seriously. The quality of the presentations and the work presented is excellent. If you’ve never attended a CONVERGE User Conference before, my advice to you is to try to be a sponge. Bring your notebooks, bring your laptops, and take as many notes as you can. The amount of useful information you will gain from this conference is enormous and more relevant than other conferences you may attend, since this event is tailored for a specific audience. The CONVERGE User Conference also draws speakers from all over the world, which provides a unique opportunity to hear about the challenges that automotive original equipment manufacturers (OEMs), for example, face in other countries, which are different challenges than those in the United States. Listening to their presentations and getting access to those speakers has been very helpful for us. And since there are plenty of opportunities for networking, you can interact with the speakers at the conference and connect with them later on if you have further questions.

Overall, the CONVERGE User Conference is a great opportunity for presenting, learning, and networking. This is a conference where you will gain a lot of useful knowledge, meet many interesting people, and have some fun at the evening networking events. If you haven’t yet come to a CONVERGE User Conference—I highly recommend making this year your first.

Interested in learning more about the CONVERGE User Conference? Check out our website for details and registration!

► Apollo 11 at 50: Balancing the Two-Legged Stool
15 Jul, 2019

On July 16th, I will look up at the night sky and celebrate the 50-year anniversary of the launch of Apollo 11. As I admire the full moon, the CFDer in me will think about the classic metaphor of the three-legged stool. Modern engineering efforts depend on theory, simulation, and experiment: Theory gives us basic understanding, simulation tells us how to apply this theoretical understanding to a practical problem, and experiment confirms that our applied understanding is in agreement with the physical world. One element does not seek to replace another; instead, each element reinforces the others. By modern standards, simulation did not exist in the 1960s⁠—NASA’s primary “computers” were the women we saw in Hidden Figures, and humans are limited to relatively simple calculations. When NASA sent people to the moon, it had to build a modern cathedral balanced atop a two-legged stool.

I like the cathedral metaphor for the Saturn V rocket because it expresses some unexpected similarities between the efforts. A medieval cathedral was a huge, societal construction effort. It required workers from all walks of life to contribute above and beyond, not just in scale but in care and diligence. Designers had to go past what they fully understood, overcoming unknown engineering physics through sheer persistence. The end product was a unique and breathtaking expression of craftsmanship on a colossal scale.

In aerospace, we are habituated to assembly lines, but each Saturn V was a one-off. The Apollo program as a whole employed some 400,000 people, and the Saturn family of launch vehicles was a major slice of the pie. Though their tools were certainly more advanced than a medieval artisan’s, these workers essentially built this 363-foot-tall rocket by hand. They had to, because the rocket had to be perfect. The rocket had to be perfect because there was so little margin for error, because engineers were reaching so far beyond the existing limits of understanding. Huge rockets are not routine today, but I want to highlight a few design challenges of the Saturn V as places where modern simulation tools would have had a program-altering effect.

The mighty F-1 remains the largest single-chambered liquid-fueled rocket engine ever fired. All aspects of the design process were challenging, but devising a practical combustion chamber was particularly torturous. Large rocket engines are prone to a complex interaction between combustion dynamics and aeroacoustics. Pressure waves within the chamber can locally enhance the combustion rate, which in turn alters the flow within the engine. If these physical processes occur at the wrong rates, the entire system can become self-exciting and unstable. From a design standpoint, engineers must control engine stability through chamber shaping, fuel and oxidizer injector design, and internal baffling.

Without any way to simulate the fuel injection, mixing, combustion, and outflow, engineers were left with few approaches other than scaling, experimentation, and doggedness. They started with engines they knew and understood, then tried to vary them and enlarge them. They built a special 2D transparent thrust chamber, then applied high-speed photography to measure the unsteadiness of the combustion region. They literally set off tiny bombs within an operating engine, at a variety of locations, monitoring the internal pressure to see whether the blast waves decayed or were amplified. Eventually they produced a workable design for the F-1, but, in the words of program manager Werner von Braun:

…lack of suitable design criteria has forced the industry to adopt almost a completely empirical approach to injector and combustor development… [which] does not add to our understanding because a solution suitable for one engine system is usually not applicable to another…

It was being performed by engineers, but in some senses, it wasn’t quite engineering. Persistence paid off in the end, but F-1 combustion instability almost derailed the whole Apollo program.

Imagine if Rocketdyne engineers had had access to modern simulation tools! A tool like CONVERGE can simulate liquid fuel spray impingement directly, allowing an engineer to parametrically vary the geometry and spray parameters. A tool like CONVERGE can calculate the local combustion enhancement of impinging pressure fluctuations, allowing an engineer to introduce different baffle shapes and structures to measure their moderating effect. And the engineer can, in von Braun’s words, add to his or her understanding of how to combat combustion instability.

Fuel slosh in the colossal lower-stage tanks presented another design challenge. The first-stage liquid oxygen tank was 33 feet in diameter and about 60 feet long. How do you study slosh in such an immense tank while subjecting it to what you think will be flight-representative vibration and acceleration? What about the behavior of leftover propellant in zero gravity? In the 1960s, the answer was you built the rocket and flew it! In fact, the early Saturn launches (uncrewed, of course) featured video cameras to monitor fuel flow within the tanks. Cameras of that era recorded to film, and these cameras were housed in ejectable capsules. After collecting their several minutes of footage, the capsules would deploy from the spent stage and parachute to safety. I bet those engineers would have been over the moon if you had presented them with modern volume of fluid simulation tools.

Readers who have watched Apollo 13 may recall that the center engine of the Saturn V second stage failed during the launch. This was due to pogo, another combustion instability problem. In a rocket experiencing pogo, a momentary increase in thrust causes the rocket structure to flex, which (at the wrong frequency) can cause the fuel flow to surge, causing another self-exciting momentary increase in thrust. In severe cases, this vibration can destroy the vehicle. Designers added various standpipes and accumulators to de-tune the system, but this was only performed iteratively, flying a rocket to measure the effects. Today, we can study the fluid-structure interaction before we build the structure! Modern simulation tools are dramatic aids to the design process.

Today’s aerospace engineering community is doing some amazing things. SpaceX and Blue Origin are landing rockets on their tails. The United Launch Alliance has compiled a perfect operational record with the Delta IV and Atlas V. Companies like Rocket Lab and Firefly Aerospace are demonstrating that you don’t need to have the resources of a multinational conglomerate to put payloads into orbit. But for me, nothing may ever surpass the incredible feat of engineers battling physical processes they didn’t fully understand, flying people to the moon on a two-legged stool.

Interested in reading more about the Saturn V launch vehicle? I recommend starting with Dr. Roger Bilstein’s Stages to Saturn.

### Numerical Simulations using FLOW-3Dtop

► FLOW-3D Workshops: Water Civil Infrastructure
18 Mar, 2020

## Online Workshop Schedule

Tuesday, March 31 — Wednesday, April 1

• 9:00am – 12:00pm EST

Wednesday, April 29 — Thursday, April 30

• 11:00am – 2:00pm EST

Tuesday, May 26 — Wednesday, May 27

• 9:00am – 12:00pm EST

Wednesday, June 24 — Thursday, June 25

• 11:00am – 2:00pm EST

Wednesday, July 22  Thursday, July 23

• 9:00am – 12:00pm EST

Our workshops are designed to deliver focused, hands-on instruction that will leave you with a thorough understanding of how FLOW-3D is used in key water infrastructure industries. You will explore the hydraulics of typical dam and weir cases, municipal conveyance and wastewater problems, and river and environmental applications. You will be introduced to more sophisticated physics models, such as air entrainment, sediment scour and transport, thermal plumes and density flows and particle dynamics. By the end of the workshop, you will have absorbed the user interface and steps that are common to three classes of hydraulic problems, and used the advanced post-processing tool FlowSight to analyze the results of your simulations. Our workshop materials are comprehensive yet accessible for engineers new to CFD methods.

## Workshop Details

• Registration is limited to 10 attendees
• Cost: $499 (private sector);$299 (government); $99 (academic) • Each workshop is broken into two 3-hour sessions • 30-day FLOW-3D license* Technical requirements • A laptop or desktop running Windows 7 or later • An external mouse (not a touchpad device) • Dual monitor setup strongly recommended • Webcam recommended ##### For more info on recommended hardware, see our Supported Platforms page. *Workshop licenses available to prospective or lapsed users in the US, Canada, and the UK. ### Who should attend? • Practicing engineers working in the water resources, environmental, energy and civil engineering industries • Regulators and decision makers looking to better understand what state-of-the-art tools are available to the modeling community • University students interested in using CFD in their research • All modelers working in the field of environmental hydraulics ### What will you learn? • How to import geometry and set up free surface hydraulic models, including meshing and initial and boundary conditions. • How to add complexity by including sediment transport and scour, particles, scalars and turbulence. • How to use sophisticated visualization tools such as FlowSight to effectively analyze and convey simulation results. ### You’ve completed the workshop, now what? We recognize that you may want to further explore the capabilities of FLOW-3D by setting up your own problem or comparing CFD results with prior measurements in the field or in the lab. After the workshop, your license will be extended for 30 days. During this time you will have the support of one of our CFD engineers who will help you work through your specifics. You will also have access to our web-based training videos covering introductory through advanced modeling topics. ### Workshop Cancellation and Licensing Policy Flow Science reserves the right to cancel a workshop at any time, due to reasons such as insufficient registrations or instructor unavailability. In such cases, a full refund will be given, or attendees may opt to transfer their registration to another workshop. Flow Science is not responsible for any costs incurred. Attendees who are unable to attend a workshop may cancel up to one week in advance to receive a full refund. Attendees must cancel their registration by 5:00 pm MST one week prior to the date of the workshop; after that date, no refunds will be given. Alternatively, an attendee can request to have their registration transferred to another workshop. Workshop licenses are available to prospective or lapsed users in the US, Canada, and the UK. Existing users are welcome to register for a workshop, and should contact sales@flow3d.com to discuss their licensing options. ## About the Instructors Brian Fox is a senior applications engineer with Flow Science who specializes in water and environmental modeling. Brian received an M.S. in Civil Engineering from Colorado State University with a focus on river hydraulics and sedimentation. He has over 10 years of combined experience working within private, public and academic sectors using 1D, 2D and 3D hydraulic models for projects including fish passage, river restoration, bridge scour analysis, sediment transport modeling and analysis of hydraulic structures. Karthik Ramaswamy is a senior CFD engineer with Flow Science, where he specializes in water and environmental applications such as stormwater conveyance, municipal infrastructure, and wastewater treatment, as well as coastal and marine applications. Karthik holds a M.S. in Aerospace Engineering from the University of Illinois at Urbana-Champaign. John Wendelbo, Director of Sales, focuses on computational fluid dynamics applied to water civil infrastructure. His range of interests includes dams, water resources, municipal conveyance and wastewater treatment, as well as environmental flows, aeration and scour problems. John graduated from Imperial College with an MEng in Aeronautics, and from Southampton University with an MSc in Maritime Engineering Science. ► FLOW-3D World Users Conference 2021 10 Mar, 2020 The FLOW-3D World Users Conference 2020 has been rescheduled for June 7-9, 2021. We invite our customers from around the world to join us at the FLOW-3D World Users Conference 2021 to celebrate 40 years of FLOW-3D. The conference will be held on June 7-9, 2021 at the Maritim Hotel in Munich, Germany. Join engineers, researchers and scientists from some of the world’s most renowned companies and institutions to hone your simulation skills, explore new modeling approaches and learn about the latest software developments. The conference will feature metal casting and water & environmental application tracks, advanced training sessions, in-depth technical presentations by our customers, and the latest product developments presented by Flow Science’s senior technical staff. The conference will be co-hosted by Flow Science Deutschland. We are extremely pleased to confirm that Hubert Lang of BMW will be the conference keynote speaker. ## Keynote Speaker Announced! ### 15 years of FLOW-3D at BMW Hubert Lang studied Mechanical Engineering with a focus on automotive engineering at Landshut University of Applied Sciences. In 1998, he started in BMW’s Light Metal Foundry in Landshut, working in their tool design department, where he oversaw the development of casting tools for six-cylinder engines. In 2005, Hubert moved to the foundry’s simulation department, where he was introduced to FLOW-3D’s metal casting capabilities. Since then, he has led considerable expansion in the use of FLOW-3D, both in the volume of simulations as well as the number of application areas. Today, BMW uses FLOW-3D for sand casting, permanent mold gravity casting, low pressure die casting, high pressure die casting, and lost foam casting. FLOW-3D has also been applied to several special projects at BMW, such as supporting the development of an inorganic binder system for sand cores through the development of a core drying model; calculation of the heat input during coating of cylinder liners; the development of the casting geometry for the injector casting procedure; and the layout and dimensioning of cooling systems for casting tools. ## BMW Museum Tour We are pleased to offer a tour of the BMW Museum as part of the conference offerings. The tour will take place at 17:30 after the technical proceedings on Tuesday, June 8. You can sign up for the tour when you register for the conference. ## Conference Information ### Important Dates • October 30: Abstracts Due • November 30: Abstracts Accepted • May 3: Presentations Due • June 7: Advanced Training Sessions • June 7: Opening Reception • June 8: Tour of the BMW Museum • June 8: Conference Dinner ### Registration Fees • Day 1 and 2 of the conference: 300 € • Day 1 of the conference: 200 € • Day 2 of the conference: 200 € • Guest Fee: 50 € • Opening Reception: included with registration • BMW Tour: included with registration • Conference Dinner: included with registration ## Advanced Training Topics Taught by senior technical staff and experts in their fields, advanced training topics include Version Up seminars for FLOW-3D CAST and FLOW-3D AM users, as well as sessions focused on Troubleshooting techniques and Municipal applications. The courses are designed so that everyone, regardless of their application, can participate in the Troubleshooting Session. You can sign up for these training sessions when you register online. ### Training Times and Fees • June 7 – 13:00 – 14:00 – Version Up: FLOW-3D CAST – 100 € • June 7 – 14:00 – 15:00 – Version Up: FLOW-3D AM – 100 € • June 7 – 13:00 – 15:00 – Municipal Applications – 200 € • June 7 – 15:00 – 17:00 – Troubleshooting – 200 € ## Call for Abstracts Share your experiences, present your success stories and obtain valuable feedback from the FLOW-3D user community and our senior technical staff. We welcome abstracts on all topics including those focused on the following applications: • Metal Casting • Additive Manufacturing • Civil & Municipal Hydraulics • Consumer Products • Micro/Nano/Bio Fluidics • Energy • Aerospace • Automotive • Coating • Coastal Engineering • Maritime • General Applications Abstracts should include a title, author(s) and a 200 word description. The new abstract deadline is October 30, 2020. Please email your abstract to info@flow3d.com. Registration and training fees will be waived for presenters. ## Presenter Information Each presenter will have a 30 minute speaking slot, including Q & A. All presentations will be distributed to the conference attendees and on our website after the conference. A full paper is not required for this conference. Please contact us if you have any questions about presenting at the conference. Flow Science Deutschland will sponsor Best Presentation Awards for each track. ## Conference Dinner This conference dinner will be held in the ever-popular Augustiner-Keller. All conference attendees and their guests are invited to join us on Tuesday, June 8 for a traditional German feast in a beautiful and famous beer garden. The conference dinner will take place following the BMW Tour. ## Travel ### Conference Hotel Maritim Hotel Munich +49 (0) 89 55235-0 info.mun@maritim.de ### Munich All about Munich Download the Munich City Map ► Training Sessions at the FLOW-3D World Users Conference 2021 9 Mar, 2020 # Advanced Training Sessions In conjunction with the FLOW-3D World Users Conference 2021, advanced training sessions will be held the afternoon of June 7 at the conference hotel. Taught by senior technical staff and experts in their fields, advanced training topics include version up seminars for FLOW-3D CAST and FLOW-3D AM users, as well as sessions focused on troubleshooting techniques and municipal applications using FLOW-3D. The courses are designed so that everyone, regardless of their application, can participate in the troubleshooting session. You can sign up for multiple training sessions when you register online. ## Version Up: FLOW-3D CAST Instructor: Dr.-Ing. Dipl.-Phys. Matthias Todte, Flow Science Germany This one-hour FLOW-3D CAST training course will begin with an introductory overview and a review of the new features and GUI design changes in FLOW-3D CAST v5.1. Through examples, new workspaces will be covered in detail, including Investment Casting, Continuous, Sand Core Making, Centrifugal as well as the new Exothermic Sleeve capabilities and database. We will also discuss the new solidification model available in FLOW-3D CAST v5.1. ### Training Details Date: Monday, June 7 Time: 13:00 – 14:00 Cost: 100 € ## Version Up: FLOW-3D AM Instructor: Raed Marwan, President, Flow Science Japan This one-hour course is open to FLOW-3D AM users as well as those interested in exploring the powerful capabilities of FLOW-3D AM for simulating additive manufacturing and laser welding processes. An overview of the additive manufacturing processes that can be simulated using FLOW-3D AM will be briefly introduced at the beginning of the training. The training will then focus on how to set up simulations for Selective Laser Melting (SLM) processes. The training will cover the powder laying process for single and multi-layer beds, powder spreading, and powder melting. ### Training Details Date: Monday, June 7 Time: 14:00 – 15:00 Cost: 100 € ## Municipal Applications Instructor: Brian Fox, MSc, Senior Water & Environmental Applications Engineer, Flow Science CFD is rapidly gaining use as an advanced tool for the design and analysis of municipal systems for stormwater conveyance and water/wastewater treatment. FLOW-3D’s well-known strengths in free surface simulation provide excellent capabilities for simulating the complex flows encountered in stormwater conveyance structures. Our multiphysics capabilities offer a powerful tool for linking the physical, chemical and biological processes that are critical for the design and analysis of water/wastewater treatment systems. In this two-hour training session we will explore FLOW-3D’s current capabilities along with recent and proposed developments for municipal applications. This class will be divided into four segments: • Review of air entrainment and two-fluid options for spiral, baffle and tangential dropshafts • Simulation of contact tanks with the reaction kinetics model • How to use the settling sludge model for clarifier technology applications • Activated sludge modeling: advanced chemistry in FLOW-3D Attendees will leave the training with in-depth knowledge of FLOW-3D’s modeling capabilities for municipal applications. For users interested in expanding their service offerings, this is an excellent opportunity to learn about capabilities for this exciting and fast growing market. ### Training Details Date: Monday, June 7 Time: 13:00 – 15:00 Cost: 200 € ## Troubleshooting Techniques This two hour training is intended for all users of FLOW-3D products, regardless of application. Instructor: Brian Fox, Senior Water & Environmental Applications Engineer, Flow Science Understanding how to identify and resolve simulation and setup issues is a critical skill for every serious CFD modeler. In this workshop we will discuss how to efficiently diagnose and address issues with FLOW-3D simulations to help keep projects moving forward on schedule. Beginning from troubleshooting techniques and the overall process, we will review the methods and practical tools available in FLOW-3D for identifying, investigating and diagnosing problems simulation errors. We will then proceed to discuss model setup options that can be used to address these issues. Throughout the class, we will apply the approach to interactively troubleshoot several real simulations to demonstrate the troubleshooting process and techniques that will help you work more efficiently on your own simulations. Finally, we will describe how to use the ideas from this training in a preventive manner in your workflow. ### Training Details Date: Monday, June 7 Time: 15:00 – 17:00 Cost: 200 € ► FLOW-3D Workshops: Brazil 25 Feb, 2020 ## Sobre nosso Workshops Curioso sobre o FLOW-3D? Deseja saber como modelamos as aplicações em hidráulicas com superfície livre mais desafiadoras: Nossos workshops são projetadas para fornecer instruções focadas, práticas e abrangentes, que deixarão você com uma compreensão completa de como o FLOW-3D é usado nas principais industrias de água & meio ambiente assim como fundição de metais. Através de exemplos práticos, você explorará o sistema hidráulico de casos típicos de barragens e represas, problemas de transporte de aguas pluviais e esgoto, rios e meio ambiente, projetos de moldes para fundição sob pressão, gravidade entre outros. Posteriormente, será apresentado os modelos físicos mais sofisticados, incluindo aprisionamento de ar, sedimentos e erosão, plumas térmicas, fluxos com densidade variável, dinâmica de partículas, perfil de movimento do pistão de injeção e enchimento do molde. No final do dia, você configurou seis modelos, entendeu a interface do usuário, as etapas comuns a três classes de problemas hidráulicos e usou a avançada ferramenta de pós-processamento FlowSight para analisar os resultados de suas simulações. ## 2020 FLOW-3D Calendário de Workshops ### FLOW-3D Workshops: Água & Meio Ambiente • 15 e 16 abril, 9:00am – 12:00pm (GMT-3) – Online • 29 e 30 abril, 9:00am – 12:00pm (GMT-3) – Online • 14 e 15 maio, 9:00am – 12:00pm (GMT-3) – Online • 14 julho – Belo Horizonte – MG • 11 agosto – São Paulo – SP • 13 outubro – Porto Alegre – RS ### FLOW-3D CAST Workshops: Fundição de Metais • 13 e 14 abril, 9:00am – 12:00pm (GMT-3) – Online • 12 e 13 maio, 9:00am – 12:00pm (GMT-3)- Online • 23 e 24 junho, 9:00am – 12:00pm (GMT-3) – Online • 28 julho – Caxias do Sul – RS ## Cadastro para o Workshop ### Política de cancelamento e licenciamento do workshop A Flow Science se reserva o direito de cancelar um workshop a qualquer momento, devido a motivos como mau tempo, registros insuficientes ou indisponibilidade do instrutor. Nesses casos, será concedido um reembolso total ou os participantes poderão optar por transferir sua inscrição para outro workshop. A Flow Science não é responsável por nenhum hotel, transporte ou outros custos incorridos. Os participantes que não puderem participar de um workshop podem cancelar com até uma duas semanas de antecedência para receber um reembolso total. Os participantes devem cancelar sua inscrição até 17:00 Brasília duas semana antes da data do workshop; após essa data, nenhum reembolso será concedido. Como alternativa, um participante pode solicitar que sua inscrição seja transferida para outro workshop. As licenças do workshop estão disponíveis apenas para usuários em potencial da área do workshop inscrito. Os usuários existentes podem se registrar em um workshop e devem entrar em contato com suporte@mettalforma.com.br para discutir suas opções de licenciamento. ### Detalhes do Workshop • O cadastro é limitado em 15 participantes • Custos:$175 (setor privado); $175 (setor público);$150 (acadêmico)
• Horário: 9:00 – 16:00
• 30 – Dias de licença para o FLOW-3D ou FLOW-3D CAST
• Traga seu laptop e um mouse e participe, ou só assista
• O almoço será providenciado por nós

### Quem deve comparecer?

• Profissionais de engenharia em recursos hídricos, meio ambiente, energia, indústrias de engenharia civil, projetistas de moldes para fundição e engenheiros de produto.
• Reguladores e tomadores de decisão que desejam entender melhor quais ferramentas de ponta estão disponíveis para a comunidade de modelos.
• Todos os modeladores que trabalham no campo de hidráulica ambiental.

### O que você irá aprender?

• Como importar geometrias e modelos de conjuntos hidráulicos de superfícies livres, incluindo malhas, condições iniciais e as condições de limites.
• Como adicionar complexibilidade, incluindo transporte de sedimentos, partículas, escalares e turbulência.
• Como usar ferramentas sofisticadas de visualização tal como FlowSight para uma análise eficiente e transmitir os resultados da simulação.
• Tópicos avançados, incluindo aprisionamento de ar e fenômenos de volume de ar, águas rasas, modelagem híbrida 3D / águas rasas e química.

### Você concluiu o workshop de um dia, e agora?

Reconhecemos que tudo não será absorvido em um dia, você pode querer explorar ainda mais os recursos do FLOW-3D configurando seu próprio problema ou comparando os resultados de CFD com medições anteriores no campo ou no laboratório. Após o workshop, sua licença será estendida por 30 dias. Durante esse período, você terá o apoio de um de nossos engenheiros de CFD que o ajudará a trabalhar com suas especificidades. Você também terá acesso aos nossos vídeos de treinamento baseados na Web, cobrindo tópicos de modelagem introdutórios a avançados.

► FLOW-3D World Users Conference 2020 Conference Announced
12 Dec, 2019

Santa Fe, NM, December 12, 2019 — In conjunction with its 40th anniversary, Flow Science, Inc. will hold the FLOW-3D World Users Conference 2020 on June 8-10, 2020 at the Maritim Hotel in Munich, Germany. Customers from around the world have been invited to the FLOW-3D World Users Conference 2020 to celebrate Flow Science’s milestone anniversary. Co-hosted by Flow Science Deutschland, this year’s conference features metal casting and water & environmental application tracks, advanced training sessions, in-depth technical customer presentations, and the latest product developments presented by Flow Science’s senior technical staff. Attendees will also enjoy a tour of the BMW Museum as part of the conference’s social events.

Flow Science has confirmed that Hubert Lang of BMW will be this year’s keynote speaker. Hubert Lang has worked in BMW’s Light Metal Foundry in Landshut, Germany since 1998. Introduced to FLOW-3D’s metal casting capabilities in 2005, Lang has led the expansion of BMW’s use of FLOW-3D. Today BMW uses FLOW-3D for a wide range of metal casting processes and special projects.

This year’s conference is particularly special. Not only are we celebrating our 40th anniversary with our customers around the world, but we are very pleased to welcome our keynote speaker Hubert Lang to honor our 15 years of partnership with BMW. Hubert will showcase some of BMW’s innovative designs for which FLOW-3D has played an indispensable role over the years. Since starting out as pioneers in computational fluid dynamics (CFD) 40 years ago, we continue to develop cutting edge software to enable customers like BMW to solve the toughest CFD problems around the world, said Dr. Amir Isfahani, CEO of Flow Science.

The call for abstracts is now open. Customers are encouraged to share their experiences, present their success stories, case studies and validations, and obtain valuable feedback from their peers and Flow Science staff. Topics include but are not limited to: metal casting, additive manufacturing, civil & municipal hydraulics, micro/nano/bio fluidics, aerospace and automotive applications. The deadline to submit an abstract is Friday, April 17.

Advanced training sessions for FLOW-3D’s family of products will be offered as part of the conference. Taught by senior technical staff and experts in their fields, advanced training topics include version up seminars for FLOW-3D CAST and FLOW-3D AM users, as well as sessions focused on troubleshooting techniques and municipal applications using FLOW-3D. Detailed information about these training sessions is available on the training page

Online registration for the conference is now available.

Flow Science, Inc. is a privately-held software company specializing in transient, free-surface CFD flow modeling software for industrial and scientific applications worldwide. Flow Science has distributors for FLOW-3D sales and support in nations throughout the Americas, Europe, and Asia. Flow Science is located in Santa Fe, New Mexico.

### Media Contact

Flow Science, Inc.
683 Harkle Rd.
Santa Fe, NM 87505
Attn: Amanda Ruggles
info@flow3d.com
+1 505-982-0088

► FLOW-3D World Users Conference Registration
25 Nov, 2019

Due to COVID-19, the FLOW-3D World Users Conference 2020 has been rescheduled for June 7-9, 2021.

# Register for the FLOW-3DWorld Users Conference 2021

### Registration Fees

• Day 1 and 2 of the conference: 300 €
• Day 1 of the conference: 200 €
• Day 2 of the conference: 200 €
• Guest Fee: 50 €
• Opening Reception: included with registration
• BMW Tour: included with registration
• Conference Dinner: included with registration

These courses are designed so that everyone, regardless of their application, can participate in the Troubleshooting Session.

• June 7 – 13:00 – 14:00 – Version Up: FLOW-3D CAST – 100 €
• June 7 – 14:00 – 15:00 – Version Up: FLOW-3D AM – 100 €
• June 7 – 13:00 – 15:00 – Municipal Applications – 200 €
• June 7 – 15:00 – 17:00 – Troubleshooting  – 200 €

Registration and training fees are waived for conference speakers (one per presentation).

### Mentor Blogtop

► Event: Aviation Electrification: Integrated design of electric motors
23 Mar, 2020

The Webinar is designed for anyone who wants to learn more about effectively approaching electric motor design with electromagnetic, thermal, system level simulation, and optimization software to reduce prototype costs and design cycle time.

► Product Demo: Using cut planes
5 Mar, 2020

Watch this short video to learn about cut planes in Simcenter Flotherm 2019.2 and how they can now work alongside with surface plots.

Have a question: Email: questions_mechanical@mentor.com

For more information on Simcenter Flotherm 2019.2, take a look at these:
https://blogs.sw.siemens.com/simcenter/introducing-simcenter-flotherm-2019-2/
https://youtu.be/ZXBFe-sr55Q

► Product Demo: Using annotations
5 Mar, 2020

Learn how to quickly and easily create the three types of annotations now available in Simcenter Flotherm 2019.2. These include plot, scalar field and object annotations.

Have a question: Email: questions_mechanical@mentor.com

For more information on Simcenter Flotherm 2019.2, take a look at these:
https://blogs.sw.siemens.com/simcenter/introducing-simcenter-flotherm-2019-2/
https://youtu.be/ZXBFe-sr55Q

► Technology Overview: How to assign static surfaces
5 Mar, 2020

How to assign static surfaces in a rotating region smartpart in Simcenter Flotherm XT. This feature allows you to define static parts within a rotating assembly.

Have a question: Email: questions_mechanical@mentor.com

► Product Demo: Using command center
5 Mar, 2020

Watch this video about using command center in Simcenter Flotherm 2019.2.

Have a question: Email: questions_mechanical@mentor.com

For more information on Simcenter Flotherm 2019.2, take a look at these:
https://blogs.sw.siemens.com/simcenter/introducing-simcenter-flotherm-2019-2/
https://youtu.be/ZXBFe-sr55Q

► Product Demo: Using automation – Part 1 of 2
5 Mar, 2020

Watch this video showcasing the new Floscript features which enable automation in Simcenter Flotherm 2019.2.

Have a question: Email: questions_mechanical@mentor.com

For more information on Simcenter Flotherm 2019.2, take a look at these:
https://blogs.sw.siemens.com/simcenter/introducing-simcenter-flotherm-2019-2/
https://youtu.be/ZXBFe-sr55Q

### Tecplot Blogtop

► How to Extract Data from a Surface, through Time
20 Mar, 2020

### Upcoming Webinar:

#### Comparing Simulation and Measured Data Through Extractions

Thursday, March 26, 2020
10:00 – 11:00 AM PDT

Register Now

## Extracting data from a surface, through time, using PyTecplot

In this example we want to extract Pressure values through time from the wall of a simulation at specific XYZ locations. After the extractions, we want to plot the results with Time on the X-Axis and Pressure on the Y-Axis for each XYZ location.

To do this we must use the PyTecplot function tecplot.data.probe_on_surface(). Another probe function, tecplot.data.probe_at_position(), is available but is likely to fail since it expects data to exist at the exact XYZ location (rather than near the location). Since a surface is infinitely thin it’s unlikely that the supplied XYZ location will be exactly on the surface, so we use the tecplot.data.probe_on_surface() function to find the surface value closest to the XYZ location.

### Extract Data Instructions:

Figure 1. View of the wing inside the boundary

1. Load data.szplt and switch to the 3D Cartesian plot type.
2. Activate the ‘wall’ zone in the Zone Style dialog.
3. Turn on the Contour zone layer and change the contouring variable to Pressure.
4. Fit the data (Ctrl+F) and rotate the view as seen in Figure 1.
5. Allow Python connections via Scripting>PyTecplot Connections…
6. From a command prompt, ‘cd’ to the directory that contains ‘probe_on_surface_through_time.py’ and run the script using the command:
> python –O probe_on_surface_through_time.py
7. The script extracts data from the XYZ locations specified in XYZ.txt and will create one zone for each XYZ location. These zones represent the Pressure values on the wall at each of those locations.
8. To plot the new data, create a new frame and switch to XY Line plot type. Create line mappings using the Create Mappings dialog (see Figure 2 below).
9. To show to timing marker in the XY plot, run the show_markergridline.mcr macro.
10. In the 3D plot, activate the three new zones and use the Scatter plot layer to display the XYZ locations in the 3D plot.
11. The final result is shown in Figure 4.

Figure 2. Create Mappings dialog

Figure 3. Pressure over time at each of the three XYZ locations on the wall

### Final Result

This Tecpot image shows the line data and the XYZ locations in context of the 3D plot.

Figure 4. Final result showing the line data and the XYZ locations in context of the 3D plot

Also see the related Tecplot 360 Video Tutorial, External Flow – Comparing a CFD Simulation with Experimental Data.

The post How to Extract Data from a Surface, through Time appeared first on Tecplot.

► Visualization of Higher Order Elements
12 Feb, 2020

## Part 1: A Higher-Order Element Primer

This blog was written by Dr. Scott Imlay, Chief Technical Officer, Tecplot, Inc.

“DON’T PANIC!” Douglas Adams, Hitchhikers Guide to the Galaxy

In the next few blogs I’ll discuss our recent research on the visualization of results from higher-order element computational fluid dynamics (HOE CFD). If you are not one of the practitioners of this specialized form of CFD you may be asking “what?” If so, it’s OK, please just heed the advice of Douglas Adams and “DON’T PANIC!”

### Partial Differential Equations

First, a quick primer on partial differential equations. Again, “DON’T PANIC!”

The behavior of a fluid flow is described by a set of partial differential equations which relate how the local state of the flow (density, momentum, energy; or alternatively pressure, temperature, velocity) changes over time and space. The simplest of these equations is the conservation of mass (continuity) equation

where ρ and u are the density and velocity of the fluid at a point in space and time (x and t). This equation is the mathematical version of the statement “mass cannot be created or destroyed”, or “if more mass flows into a region than flows out, the density must increase over time.”

If you’re not familiar with symbols in this equation, DON’T PANIC!

The symbol represents a partial derivative and is the rate at which the density, ρ, is changing in time. Likewise, is the rate at which momentum (density times velocity) increases with increasing spatial dimension x. Another way to think of the partial derivative is as the slope of the curve:

Figure 1. ∂ρu/∂x is the rate at which momentum (density times velocity) increases with increasing spatial dimension x, which is the slope of the curve.

In the above chart the red line is the variation of ρu with x and the slope (change in ρu over change in x) of the tangent green line is at the point of tangency.

### Approximating the Derivative

Of course, we don’t generally know a functional form for ρu in terms of x, so we need to approximate the derivative. This is where the order of accuracy comes into play. In solving these partial differential equations we use a grid, with Δx being the grid spacings between grid points. If our approximation to the derivative, call it , has an error term proportional to Δxn we call it an n-th order accurate scheme. For example, a second order scheme would approximate the derivative with an error proportional to Δx2

where the … represents terms containing higher exponents of ∆x like Δx3, Δx4, and so on. As the grid spacing Δx gets smaller the error terms with exponents higher than second order will vanish more quickly than the leading Δx2 term, so the approximation is called second order accurate. The majority of CFD codes are second-order accurate.

### Higher-order Schemes

Higher-order schemes have a leading error of greater than second order. For example, a third-order scheme will have a leading error term of Δx3. As the order of the scheme gets higher, the approximation error drops more rapidly as the grid spacing Δx gets smaller. A higher-order scheme gives a more accurate solution than a second-order scheme with a given grid spacing. Alternatively, you can say that the grid spacing doesn’t need to be as small to get a desired level of accuracy. It is the ability to use a coarser grid to get the same level of accuracy that makes the higher-order schemes attractive.

Figure 2. It is the ability to use a coarser grid to get the same level of accuracy that makes the higher-order element schemes attractive.

### Finite-element Methods

For finite-element methods, the variation of the solution over the element is generally approximated as a polynomial with unknown coefficients. For example, in two-dimensions, the element shape may be triangular with a linear variation of solution over the triangle. The goal of the finite-element method, therefore, is to solve for coefficients of the polynomials for all elements in such a way to give the best approximation to the solution of the partial-differential equations.

Enough about the method, let’s focus on the finite-elements themselves!

### Why Use Higher-order Elements?

The degree of the polynomial determines the order of accuracy of the scheme (at least the spatial derivatives). For example, linear elements are second order accurate, quadratic elements are third-order accurate, and cubic elements are fourth-order accurate. In general, if p is the degree of the polynomial approximation then p+1 is the order of accuracy.

Why use higher-order elements? In short, for a large class of problems, HOEs achieve a desired level of accuracy with less computation expense. The following chart (Figure 3) from the AIAA Higher-Order Workshop shows how this is true. The horizontal axis of this chart is “work units” (how much computational resources the code uses) and the vertical axis is the solution accuracy (not order-of-accuracy, but how close the solution is to the correct solution). As you can see, the “_P3_” results for each of the codes tend to perform better than the “_P1_” or “_P2_” results. Here “P3” indicates a cubic polynomial which would yield a fourth-order accurate scheme. Note that in the chart below, down is more accurate.

Figure 3. Higher-order elements achieve a desired level of accuracy with less computation expense. This chart, from the AIAA Higher-Order workshop, shows how this is true. The horizontal axis of this chart is “work units” (how much computational resources the code uses) and the vertical axis is the solution accuracy (not order-of-accuracy, but how close the solution is to the correct solution). Down is more accurate.

Legend for Figure 3.

Unsteady flows with propagating vortices are the class of flows that benefit the most from higher-order methods. This is the type of flow you often get with large-eddy simulations, which are becoming more popular as computers get more powerful. For this reason, the usage of higher-order methods will likely increase in the future. Visual Analysis codes like Tecplot 360 need to be ready.

In the next blog I’ll discuss the challenge of visualizing higher-order finite-element solutions.

### Subscribe to Tecplot 360

Get notified when the next blog is posted!

## Subscribe to Tecplot 360

The post Visualization of Higher Order Elements appeared first on Tecplot.

► Tecplot Macro Tutorials
22 Jan, 2020

This blog on Tecplot Macro Tutorials was written from a webinar entitled Ask the Expert About Tecplot 360 hosted by Scott Fowler, Tecplot 360 Product Manager. We received numerous questions about Tecplot macros and this blog pulls our macro resources together in one place.

What is a Tecplot macro? A Tecplot macro is a set of instructions, called macro commands, that perform actions in Tecplot 360. Macro commands can be used to accomplish virtually any task that can be done via the Tecplot 360 interface, offering an easy way to automate Tecplot 360 processes.

We recommend that you download the macros, datasets and layout files used in the descriptions below. If you do not already have Tecplot 360 running, you can download a Free Trial.

Macro, Data and Layout Files (17MB ZIP)

## Using a Macro to Produce Transient Plots

This macro reads and edits different time steps and produces transient plots. This macro first finds out how many time steps are in your dataset, then it loops over those time steps. At each time step the current time step is set and an image is exported. That’s It!

export_over_time.mcr
$!EXTENDEDCOMMAND COMMANDPROCESSORID='extend time mcr' COMMAND='QUERY.NUMTIMESTEPS NUMTIMESTEPS'$!LOOP |NUMTIMESTEPS|
$!EXTENDEDCOMMAND COMMANDPROCESSORID='extend time mcr' COMMAND='SET.CURTIMESTEP |LOOP|'$!EXPORTSETUP EXPORTFORMAT = PNG
$!EXPORTSETUP EXPORTFNAME = "FILE_|LOOP|.png"$!EXPORT

## Converting Binary Files to ASCII

A question asked during the webinar was “How do I convert TEC files to binary?” The .tec file extension has been adopted by the greater Tecplot community. These files are usually ASCII, but in this case, the TEC file was a binary file. So, the question becomes “How do I convert from binary to ASCII?” Well, that’s easy!

Let’s backtrack for a moment and review the canonical Tecplot file extensions:

• PLT (.plt) – Tecplot binary format.
• DAT (.dat) – Tecplot ASCII format.
• SZL (.szplt) – Tecplot Subzone Load-on-demand binary format. This format allows you to load large volume metric grids with very little RAM.

We don’t have a standalone utility to convert binary files to ASCII, but you can do it using the macro binary_to_ascii.mcr (also included in the ZIP file mentioned above). In the macro, you simply use the ReadDataSet command and then a WriteDataSet command.

binary_to_ascii.mcr
#!MC 1410
$!ReadDataSet '"VortexShedding.plt" ' ReadDataOption = New ResetStyle = No VarLoadMode = ByName AssignStrandIDs = Yes$!WriteDataSet  "VortexShedding.dat"
IncludeText = No
IncludeGeom = No
IncludeCustomLabels = No
Binary = No
UsePointFormat = No
Precision = 9
TecplotVersionToWrite = TecplotCurrent

From the command line you would call tec360.exe and pass this macro. Note that the file names in the macro need to be hard coded.

>tec360 -b -p binary_to_ascii.mcr

Here a plug for the power of Python!
Using PyTecplot to convert files, in this case binary to ASCII, is so much simpler as you can see in these few lines of code. This Python script (binary_to_ascii.py)is in the ZIP file mentioned above.

import tecplot as tp
tp.data.save_tecplot_ascii(‘myasciifile.dat’)​

## Converting ASCII Files to Binary

Another user wanted to convert ASCII files to binary.

You can use a utility called Preplot. Preplot is included with Tecplot 360 and it’s located in the bin directory in the Tecplot 360 installation. You just pass it an ASCII file (.dat) and tell it the output file name. Preplot will do the conversion for you.

> preplot myasciifile.dat mybinaryfile.plt

A note about ASCII data, whenever Tecplot 360 loads an ASCII file we actually convert it to binary as it’s being loaded. Because of this, it is in your best interest, if you need to load that file over and over, to do the conversion only once. You will save yourself a lot of time!

## Placing Streamlines and Streaklines Using Macros

During the webinar, one user was trying to create an arc of streamlines (which we call streamtraces in Tecplot 360). This is not possible using the onscreen “Rake” tool! The Rake tool only places straight lines of streamlines.

If you want to define an arc of streamlines, the best way to do it would be with PyTecplot (but that will be the topic of a future blog). I explain here how to do it with a Tecplot macro.

Here is the Tecplot macro for placing the streamtrace object associated with the frame. In the code below, I’m adding a streamtrace, defining it as a volume line, setting the direction to “both” (default is forward), and placing it at X= 0.5, Y=0.5, Z=0.5. You can create an arc by having multiple of these add commands at different locations.

NUMPTS = 1 ​
STREAMTYPE = VolumeLine​
STREAMDIRECTION = Both​
STARTPOS { X = 0.5 Y = 0.5 Z = 0.5 }​

How do I share data across frames without having to reload the data? In this example, I want to tile frames and then save the style from the upper left plot and load it into the lower right frame. You can do this, but not directly in the Tecplot 360 user interface. We will use a macro to do this (It can also be done with PyTecplot).

### Copy Frame Style Steps

1. First, create a new frame (click the Frame icon and draw a new frame). Then change the plot type to 2D Cartesian. (Now I could easily click on the upper left frame and save the frame style (Frame>Save Frame Style). Then, I could click in the new frame and load the frame style, but that will also copy the frame size and position which will overlay the original position. I want to retain the position. So, I don’t want to do it this way.)
2. Drag and drop the copy_frame_style.mcr macro (from the ZIP file above) onto the Tecplot 360 user interface. Two new entries will appear in the Quick Macro Panel: “Copy frame Style” and “Paste frame Style.” Be sure the stylesheet directory in the macro is valid (C:\TEMP\temp_style.sty”).
3. Click on the frame you want to copy, click the Quick Macro “Copy frame Style” and click Play (or simply double-click on the macro name).
4. Select the new frame and double click on the “Paste frame Style” macro.
5. Voila! You now have a new frame with the same style. This is one of the macros I use most often.

### Copy Frame Style Video Tutorial

The magic is in the macro’s last command “INCLUDEFRAMESIZEANDPOSITION=NO,” which loads the frame style but ignores the position.

## Making Macros Persistent in Tecplot 360

If you want these macro functions always available, you can put them in the tecplot.mcr file (located in Tecplot 360 installation directory).

If you are sharing the Tecplot 360 installation, for example in a Linux environment, you can put the “.tecplot.mcr” file in your Linux home directory so as not to impact other users on your network.

More information can be found in the Tecplot 360 User’s Manual (search for “tecplot.mcr”).

The post Tecplot Macro Tutorials appeared first on Tecplot.

10 Jan, 2020

In this Webinar, customers ask questions of Tecplot 360 Product Manager, Scott Fowler. Scott will help you develop your own internal expertise as he answers your questions interactively. You’ll learn new and improved techniques, current best practices and implementation considerations.

Topics include:

• Streamlines and Streaklines
• Macros, PyTecplot
• MATLAB
• Frames & Styles

Download the files used in the webinar containing the presentation, data files, macros and Pytecplot scripts.

The post Ask the Expert about Tecplot 360 – Learn New Techniques! appeared first on Tecplot.

► 11 Questions About Tecplot 360 Basics
3 Dec, 2019

These questions were asked during the Webinar, From Zero to Hero: Tecplot 360 Basics. Tecplot 360 Product Manager, Scott Fowler, provides the answers below. Most of the answers show Webinar timestamps so that you can follow along in the Webinar. Learn more about Tecplot 360.

### 1. How do you save all the steps to a macro?

The steps taken in the Webinar (timestamp 36:55) can be recorded to a macro file. Macros are the legacy scripting language of Tecplot 360, and the language on which many of our file formats are based.

• You can customize the Tecplot 360 user interface using the quick macro panel.
• You can customize defaults using what we call a configuration file, tecplot.cfg.
• You can save your current work by saving to a layout file. If you have created new variables, you will also be asked to save the data file. The saved layout file is in this macro language.

Macro Video Tutorials

Tecplot 360 also has a Python API, called PyTecplot, which is available to customers on TecPLUS maintenance. Learn more about PyTecplot.

### 2. Why is the cylinder surface not showing velocity magnitude?

The cylinder in the Webinar example (timestamp 38:07) is a no-slip wall. Because it is a boundary condition, there are no velocities on the wall. If I want to see the velocities near the wall, I could change J planes to, for example, J=2.

### 3. Can we plot stagnation energy?

Yes. That is a variable that can be computed. From the main menu (Webinar timestamp 38:46), select Analyze>Calculate Variables and click Select…, then choose Stagnation Energy. See Section 21-3.2 Identifying State Variables in the Tecplot 360 User’s Manual.

### 4. How do I export a figure without a borderline?

You can hide the border by editing the active frame. To edit a frame (Webinar timestamp 36:55), right-click on the edge of your frame, then choose Edit Active Frame. Uncheck Show border.

Some graphics cards will have a little drop shadow on the right and bottom of plots. This is graphics card dependent. If you see a drop shadow after turning off borders, please see this article on how to correct it in the Tecplot Knowledge Base.

### 5. When I save a layout and move the layout to another folder, it won’t load.

When you use the menu command File>Save Layout As, there is a toggle (under the Save as type) that says Use relative paths (Webinar timestamp 36:55).

• Checking the toggle saves the referenced data layout file as a relative path.
• Unchecking the toggle saves the layout file as an absolute path. You should be able to load your data when it is saved as an absolute path.

You can also edit the layout file because they are saved as human readable macro files. (timestamp 39:50). I’ll quickly open a macro file and you can see what the Tecplot macro language looks like.

Here is the layout that I saved earlier. You can see that I have a command to point to linedata.plt. To move this file to a different folder, I need to change the file path to an absolute path. I can update the path in the macro file.

Relative paths are often used by people doing optimization studies when they have a folder hierarchy with, for example, Mach alpha sweep and data within each sub folder. Using a relative path, they can copy a single layout into each subfolder and load that data by relative path.

### 6. Can you repeat extracting the line plot?

Yes. Use Tools>Probe To Create Time Series Plot, then single click on the plot, and that will update the time series plot (Webinar timestamp 41:20). You can see that as I single click, the line plot updates through time.

For more on extracting, see the Tecplot 360 User’s Manual

### 7. Is it possible to export a cropped image?

Images cannot be cropped in the traditional sense. However, selecting File>Export from the main menu gives you some options (Webinar timestamp 41:53). You can export three different regions of your plot:

• All frames – Exports all frames in your workspace.
• Current frame – Exports only the current frame.
• Work area – Exports entire “gray” region of your workspace.

To export only a portion of your image, you can do what we call a paper zoom. Zoom in on your plot by holding down the shift key, middle mouse button while moving the mouse up. A work area export will export the zoomed-in plot you see in the work area.

### 8. How can I get two plots in the same frame?

In XY line plots you can plot up to five X-axes and up to five Y-axes (Webinar timestamp 41:20).

In this case we’ll use multiple Y-axes. In the line plot, go into the Mapping Style dialog from the Plot sidebar, and you will see that there are multiple Y axes. For this example, select RHO and Velocity Magnitude to plot in the same frame, then press CNTL+F to fit the view. You can see that these have two very different scales. In the Mapping Style dialog, set RHO to Y2, and now row is on the Y2 axis.

This will give you two plots in the same frame. Select the adjuster tool to move the label.

### 9. Can I add LaTeX symbols to the legend?

You cannot use LaTeX symbols in a legend, but you can add LaTeX annotation to your plot and position it near the legend (Webinar timestamp 47:16). Here is how to do that:
• Select the text tool .
• On your plot, click where you want to position the text, which brings up the Text Details dialog.
• Press the LaTeX button (upper right side of Text Details dialog).
• Add the LaTeX annotation, set the Size in points, and press Accept.
• Click on the legend, which brings up the contour & Multi-Coloring Details dialog.
• Click on the Legend Box… and click No Box.
• Close the dialog.
• Click on your LaTeX annotation and move it above the legend.

That will mimic adding LaTeX symbols to the legend.

Watch the Video on LaTeX Fonts

### 10. What is the Shade used for?

For information on translucency and lighting zone effects refer to Chapter 13, and for information on shade, refer to section 12-1 in the Tecplot 360 User’s Manual.

### 11. How do I perform a Fourier Transform?

In XY line plots, navigate to the Data>Fourier Transform… menu option. Follow along in this video for a demonstration: Loading Excel Data and FFT in Tecplot 360.

The post 11 Questions About Tecplot 360 Basics appeared first on Tecplot.

► Multiple Domains, Multiple Scales, One Visualization Tool
13 Nov, 2019

Case study contributed by Michael Callaghan, PhD, P.Eng – Senior Applications Engineer, Aquanty

Aquanty is a leading-edge water resources science and technology firm specializing in predictive analytics, simulation and forecasting, research services, and IoT. Aquanty’s solutions and services are deployed globally across a broad range of industrial sectors including; agriculture, oil and gas, mining, watershed management, contaminant remediation, and nuclear storage and disposal. Aquanty’s flagship platform, HydroGeoSphere, is a class leader in fully integrated three-dimensional surface/subsurface modeling.

## What Is an Integrated Surface-Subsurface Hydrologic Cycle Model?

Figure 1. Hydrological Cycle.

HydroGeoSphere (HGS) is a three-dimensional control-volume finite element simulator which is designed to simulate the entire terrestrial portion of the hydrologic cycle. It uses a globally-implicit approach to simultaneously solve the 2D diffusive-wave equation and the 3D form of Richards’ equation.

The basis for HGS’ integrated computation is multiple 1D, 2D, 3D ‘domains’ that interact with each other, including:

• a 2D overland flow domain,
• a 3D subsurface flow domain that can include separate discrete fractures and dual permeability domains,
• as well as 1D surface flow channels,
• 1D subsurface tile drains,
• and 1D water wells.

Data intensive model output for the hydrologic cycle benefits from the flexibility of Tecplot 360 for visualization, which often requires plotting of multiple domains simultaneously and in 3D.

“We chose Tecplot 360 because of the quality of the plots. We need to present our results to clients, and for that plot quality there is no substitute to Tecplot 360.”

– Michael Callaghan, PhD, P.Eng, Senior Applications Engineer, Aquanty

## Visualization of Surface Water – Groundwater Interaction

Simulating groundwater-surface water interaction in complex topography such as hummocky terrain has traditionally been viewed as a significant challenge to the hydrologic modelling community, for example so-called fill and spill behavior.

However, with HGS, the complex processes by which water movement is influenced by the combination of surface topography and highly variable subsurface hydrostratigraphy or preferential flow pathways can be readily reproduced.

In the example shown in Figure 2, precipitation is falling on an upland area which then initiates overland flow, filling, and spilling of surface depressions. The Tecplot 360 animation illustrates depression-focused groundwater recharge occurring beneath the depressions, with both a perched water table and a fractured aquitard influencing sub-surface water movement.

Figure 2. Visualization of surface water-groundwater interaction.

## Flood Inundation Visualization

HGS models may range across a number of scales from centimeters to meters to kilometers to 100’s of kilometers. Use of an unstructured finite element mesh makes this possible.

Tecplot 360’s inherent flexibility with unstructured grids makes it a very useful visualization tool across many scales of problems.

In this application, HydroGeoSphere is being used to recreate the Southern Alberta, Canada flooding that occurred in June 2013. The simulation results presented in Figure 3 depict a flood pulse derived from basin scale hydrologic simulations being routed across a local scale model of the City of Medicine Hat, Alberta, Canada. LiDAR-derived topography was used as input for this model, and results show excellent agreement between simulated and observed high water marks.

Figure 3. Flood inundation visualization.

The above visualization cases are made possible with Tecplot 360. HGS is highly integrated with Tecplot ecosystem using a powerful post-processing tool to directly produce results in Tecplot file formats. Tecplot 360 is an essential tool for everyday 2D data plotting, high quality 3D model visualization, results inspection and evaluation.

The post Multiple Domains, Multiple Scales, One Visualization Tool appeared first on Tecplot.

### Schnitger Corporation, CAE Markettop

► How your tech skills can help in the pandemic
29 Mar, 2020

# How your tech skills can help in the pandemic

The world may be a sucky place right now, but people are (mostly) awesome. Protocol just published a list of projects that can use your/my/everyone’s tech wizardry. That inspired me to go through my inbox and pick out all of the other ones I’ve been emailed about so far. If you’re starting something and can use help, tell me and I’ll add it here, but also tell Protocol (their writers’ email address is at the bottom of the first link below), post it on GitHub, your alumni news site (if that’s appropriate), LinkedIn — anywhere you can to get the word out. We need to get a lot of people working on this NOW.

As of 29 March,

• Protocol’s list includes pandemic data visualizations at Johns Hopkins, COVID-19 test tracking at USC, a project to 3D print face shields at Columbia University (and they need drivers, so if you’re in NYC …)
• LinuxInsider has another great list, here, for epidemiological modeling, hospital capacity planning, apps to help people feel less isolated during lockdowns
• There’s a Help with Covid website, where projects are advertizing for help and some individuals are advertizing for projects
• NYU has established the OpenFacePPE project, which seeks to connect medical supply chain players who can 3D print an open-source face shield (PPE stands for personal protective equipment and includes all sorts of devices, not just face shields)
• Gerber Technologies founded the Gerber PPE Task Force to help manufacturers increase production or convert to producing PPEs
• Opensource.com has yet more software projects to get involved in
• GraphicSpeak (the Jon Peddie newsletter) published an article about an MIT project on open source respirators — perhaps there’s a way to get involved there
• Then there’s Folding@home, which uses your PC’s or gaming system’s unused GPUs to work to help scientists develop models and, we all hope, a vaccine. I haven’t done this, but the instructions are here
• Stratasys has put together a network of 3D print shops and universities to print face shields. If you can do at least 100, contact them!

I’ll be adding to the list as I learn about more projects.

And if you’re not a techie, there are lots of things you can do, too. This PBS article has some terrific ideas.

Stay safe, don’t hoard and wash those hands!

The post How your tech skills can help in the pandemic appeared first on Schnitger Corporation.

► AVEVA adds production accounting toolset
26 Mar, 2020

# AVEVA adds production accounting toolset

AVEVA just announced that it has acquired accounting software from MES Enter to, as the company puts it, “complete AVEVA’s value chain optimization solution”.

The new AVEVA Production Accounting (fka MES ENTER ErrorSolver) becomes part of AVEVA’s Value Chain Optimization offering, to help customers integrate business processes across their operations. I saw the Value Chain Optimization tools at AVEVA World Summit last year — the idea is that one can make better decisions by knowing all of the knowable bits of information beforehand, to mitigate the risk inherent in the unknowable. So, for example, say an integrated oil producer believes it needs x gallons/day of heating oil to satisfy demand in the Northeastern US in March 2020. Backing off from that forecast, what can it buy at what price on the spot market in December 2019? Alternatively, what can it produce in October 2019? where can it best refine that crude into heating oil, given utilization rates at its various facilities? Which is cheaper, including transportation costs? What should the risk multiple be in any of these scenarios? Answering these questions could lead to a giant Excel spreadsheet, with dozens of sensitivity analyses, formulas, data sources and the individual quirks of each planner — or you could use a commercial solution like AVEVAs.

ExxonMobil, as just one concrete example, spoke about using advanced planning tools across its operations (aka value chain) –planning, scheduling, operations (including IoT data to provide a real-time, as-is view)– using in-house developed tools as well as commercial tools from vendors like AVEVA. AVEVA says this acquisition can help customers improve the accuracy of planning models, and I can see how that could be true.

What’s interesting, to me, is that this specific addition to the value chain tools is about production accounting, not accounting accounting. Of course, it all comes down to money, but production accounting is about mass balances: we have this much raw material in our tanks and draw it down at a specific rate (as shown by gauges) and then we turn it into output products, as measured by these other gauges. What happens in between, adding catalysts or heat or pressure or other feedstocks, affects the mass balance of the equation. If a plant knows what its inputs and outputs are, and how quickly it is using up its stock of raw materials, it can better plan. And it can do all of the ancillary things, too, like look for pilfering, leaks, suppliers overstating what they deliver. Better visibility leads to all sorts of benefits. (I suppose this fascinated me because it’s physics meets finance. Hmm.)

Anyway, the deal appears to be done and no financial details were announced. It’s another sign that AVEVA is expanding outwards from the legacy AVEVA’s engineering domain to serve its customers with the solutions they need to efficiently run their businesses.

The first version of this was incorrectly headlined “project accounting”. Naveen Kumar over on LinkedIn pointed out the error.

The post AVEVA adds production accounting toolset appeared first on Schnitger Corporation.

► ESI on the impact of Covid-19: “It’s a 3-month event”
24 Mar, 2020

# ESI on the impact of Covid-19: “It’s a 3-month event”

ESI Group held an investor update today, in which they expanded on the 2019 results reported last month. (You can read my recap here.) I’ll write more about the company’s vision, strategy, etc. as discussed in today’s session, but I wanted to share one very important note ASAP, while I get my thoughts together on the rest.

In a press release yesterday, ESI CEO Cristel de Rouvray said,

In the short-term, the disastrous coronavirus pandemic is expected to somewhat impact our H1. However, the resilience of our business model largely anchored on renewable and mission-critical software licenses will help us manage full-year risks. When industry recovers from this exceptional crisis, digital ways of working will be accelerated globally, fully dependent on ESI’s solutions to virtually anticipate and manage asset performance in-service, much beyond the traditional PLM certification target of the brand-new product.”

Today, she and CFO Olfa Zorgati elaborated on that theme. Ms. de Rouvray said that ESI, as a global company, has already been through the “Covid cycle” in Asia, and saw delays (not canceling) of projects to later in the year. Ms. Zorgati said that she’s seen a temporary dip in software utilization as well as a delay in sales activity, but reiterated that nothing has been canceled. Ms. de Rouvray wrapped up her remarks on Covid by saying that (paraphrasing),

If Asia, which is further along, is a predictor for the rest of the world, this is a 3-month event. There is some business inertia, but our compelling value means we are a protected supplier.

I don’t think ESI is alone in this –compelling value, mission-critical, key to new digital ways of working once this is all behind us– and believe that other PLMish companies are seeing similar patterns.

This week is the end of the fiscal quarter for most of the companies in our little corner of the IT world. If Q1 was a disaster, we should know as early as next week, if companies pre-announce earnings shortfalls. But I honestly don’t see that happening, as China is a relatively small part of most PLMish companies’ revenue and business elsewhere was decent until just a few weeks ago. I think most will follow ESI’s lead: business pushed out into the second half of the year. We might even see an uptick as companies look back over the last few weeks and realize how unprepared they were; perhaps investing in new IT solutions (some PLMish) to help them with supply chain agility, remote working and other strategies to better cope when the next crisis hits.

In the meantime … I hope you’re well, taking care of family, friends, and neighbors — but from a decent distance. I am once again amazed at how our PLMish community comes together to support, well, everyone. Free training classes, extra cloud credits, expanded product try-out periods … I’ve heard from so many vendors who are trying to help that I can’t list them all. If you’re at all interested, now is a great time to try a product you’ve always wanted to, learn something new, or mess around with a design problem that’s been driving you nuts. Or you can watch puppy videos. Why does the one on the right strut as if he owns the world? Because he probably does. Alpha dog. Be well, people.

The post ESI on the impact of Covid-19: “It’s a 3-month event” appeared first on Schnitger Corporation.

► Totally OT: Capybaras balance oranges
24 Mar, 2020

I could (and probably should) write something PLMish but let’s save that for later. Instead, I’m going to share a moment of zen, courtesy of HIRO@sea:

If you speak Japanese, you can understand the naturalist’s explanations; if not, just watch the capybaras float serenely around their pool. The focus shifts to the swimmers at about 0:45, and this handsome gentleman (?) does his trick at 1:45. Such a ham!

That’s it for now – ESI holds its investor event (virtually) in a few hours, and I’ll share what I learn tomorrow.

The post Totally OT: Capybaras balance oranges appeared first on Schnitger Corporation.

► Funny & wise from Mel Brooks’ family as we all try to flatten the curve
17 Mar, 2020

Last one for today, I promise. Too good not to share:

The title image is of a much younger Mel Brooks, in the classic movie, Spaceballs. So goofy.

The post Funny & wise from Mel Brooks’ family as we all try to flatten the curve appeared first on Schnitger Corporation.

► PLMish ideas on whiling away social isolation
17 Mar, 2020

# PLMish ideas on whiling away social isolation

Right up front: I’m not qualified to give medical advice so you won’t see any of that here. Go to cdc.gov, who.int, or another reputable source for the protocols that will ensure your safety and that of elders and those with compromised immune systems.

That said, why am I in your inbox? Because I have a list of things to do and think about while you’re at home, flattening the curve:

• Get some training. No matter what you do, you can always learn more. Use this time to upskill. Work towards that next level SolidWorks certification, learn something about simulation, dive into an adjacent trade skill, explore 3D printing, get a demo license to something you’ve been wanting to try … PLM is a huge, complicated set of technologies and more is always coming at us. Use this time to take legit distance learning classes, watch some of the thousands of webinars available for replay or get current on that giant to-read pile.
• Teach others. Marvelous tech can connect us, virtually, to help spread what we know. Hook up with a local school that’s doing distance learning to show off your CAD wizardry. Use Skype to teach neighbor kids how to make the coolest Lego thingamabob ever (or an exploding volcano in the kitchen, but don’t mention my name). Run a coffee or lunch seminar for your colleagues, to stay connected and to share what you know.
• If you’re part of a supply chain, you’re likely struggling. Parts you need may be on their way from China, but coming to an idled plant in Europe or the US. Or they’re not being made, at all, and you’ll face huge disruption when you can get back to work. Use this time to ponder your supplier choices. Too many manufacturers single-source today, because the volume deals are attractive and because it’s hard to collaboratively design with too many partners. That’s great when things work to plan. When they don’t, it’s a disaster. Think about how to create a more resilient supply chain — and PLMish tools can help explore options. Plow through your bills of material — what alternate components have you used? Who supplied them? Can they be spun up again? Combine this with reshoring — who, more locally, can you work with, for the benefit of your community?
• If you’re part of a larger organization, even at the lowest rung on that ladder, think about how it all fits together. Most of us do, what we do, the way we do it, because of choices made decades ago. Marketing does concept design, engineering adds the details, manufacturing makes things. What if we changed those roles, maybe for upfront CAE? Would marketing now do CAE? (Mind blown. Yes.) What would need to change to make that happen? Think beyond that example; what benefits might you see?
• Also, think beyond your role. How can you use the data in your ERP system to make better choices about suppliers, customer requirements, which facility makes what … What info would make your job more satisfying and lead to a better result? Where does that live? How can you gain access to it, what would you do with it, how can you convince your boss to make it available to you? PLM is one bucket of data; there are other buckets in your company. How can you dip into multiple buckets?
• Imagine your next product. Look at this as the gift of time. At the very least, you’re gaining the time you would have been commuting. Don’t read more news; that’ll just add stress. Use the time to sketch out ideas, dream up new offerings, check out the competition’s website, muse on what could be when the world is back to some sort of new normal.
• Be personally strategic, too. What do you need to know to get that next job? It might be related to the training I suggested above, but now might also be a great time to think about what you need to learn/become to climb that corporate ladder, open your own business or follow other dreams.
• Contingency plan. I was at an investor event last week (foolish in hindsight, but we know so much more now than we did even a week ago) at which one CFO spoke about her company’s trial runs of business systems, in case of a prolonged shutdown. This particular crisis is here, but what can you do differently next time? Because there will be another blizzard, fire, hurricane or other crisis. My Congressman has this great plan on his website, with three threat levels and a simple set of actions for each, take a look. He’s not PLMish in the slightest, but what can you adapt to your situation? What can you ask people to do remotely? What tech do you need to make that happen?

Above all, breathe. Everything I’ve seen, from economists, consultants, and bankers, points to a rapid economic recovery once we’re all able to get out again. At first, it’ll be a bit panicky as we spin things back up*. Then it’ll be incredibly competitive, as everyone scrambles for consumer and business dollars. Spend this sort-of downtime getting ready.

*It occurs to me that not everyone can relate to that. Sigh. Back in the day, computers had disk drives –huge things that had flat platters on which data was stored– that had to spin at a specific number of revolutions per minute to read and write data. So turning on a computer meant, literally, waiting for those to “spin up”.

The title image is of a cherry blossom tree by Majaranda from Pixabay.

The post PLMish ideas on whiling away social isolation appeared first on Schnitger Corporation.

return

### Layout Settings:

 Entries per feed: 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40 Display dates: No Yes Width of titles: 300 350 400 450 500 550 600 650 700 750 800 850 900 950 1000 Width of content: 400 450 500 550 600 650 700 750 800 850 900 950 1000