Calculating Wave Speeds in Flownex

Categories:

Today, I will be discussing how to use a quick script component to calculate wave speed for a pipe in Flownex. We will be developing a script to determine the wave speed in two elastic tubes in different scenarios. Both networks are dependent on the 1D tube geometry and mechanical properties. Longitudinal pressure waves may be modeled in Flownex within a wide variety of materials or timescales.

First, we will calculate the pulse wave velocity (PWV). This is a non-invasive method of determining arterial stiffness using the blood pressure measured at to different arteries, used as an indicator of cardiovascular health. Next, we will apply the script to a pipe with a sudden valve closure to calculate the wave speed and the reflection time of the pressure wave. For this example, we will be using Flownex version 8.14.1.4845.

Choosing Components and Boundary Conditions

We will be simulating the wave speed of due to an influx of 60cc of blood (I will be creating a custom fluid for this) over a 0.1s interval (throughout our simulation) to produce a pressure wave. This arbitrary pressure wave will allow us to visualize the wave speed of the artery. We can simulate the artery by dragging and dropping a pipe component onto the canvas and placing boundary conditions on the inlet and outlet.

Next, we will specify a zero-mass flow at the inlet, and set the temperature to 35°C.

image 3

We will use a constant pressure for our outlet boundary condition.

image 4

Lastly, we will assign our inlet node to contain a volume of 60cm3.

image 5
image 14

Next, we can define the vessel geometry. We will check “Specify advanced inputs” in the pipe properties, allowing us to access the “Transient and Water Hammer Related Inputs” subset of mechanical properties. We will model a simplified aorta with a diameter of 25mm and a length of 30.48cm. Additional mechanical properties of the pipe are shown below:

image 9

Flownex will use the diameter to thickness ratio to determine if a thin wall model is necessary for the simulation. A diameter to thickness ratio above 40:1 will use a thin-walled pipe model.

Transient Actions Setup

Now that we have defined our boundary conditions, we can setup an action to perform this task by navigating to the configuration tab in the ribbon and clicking on Action Setup. Additional information on setting up transient actions form a table can be found here.

We will create a New Scenario in the actions setup and set the target property to the volume of the inlet node. We can do this by selecting the node and by dragging the volume from the inlet node properties (left) window into the actions setup (right).

image 8

Setting the value type to table will allow us to create a constant volume flow rate at the inlet after the no flow initial condition, using the input below.

Calculating Wave Speed

We can also easily calculate the wave speed of a component, using a quick script. Quick scripts are easily implemented, and we will specify unit types to allow Flownex to evaluate unit conversions for calculating wave speeds. I won’t cover quick script/data transfer link setup in this post, but you can find a previous example here. After inserting a quick script component, we can double click on it to modify the inputs, outputs, and calculation.

image 10

Rather than recomputing the speed of sound of our fluid, we will take it directly from the properties of the pipe. Additionally, we can calculate the time when the wave will reach our boundary. The PWV along multiple pipe components (such as carotid-femoral PWV or examining a pipe main in a complex network) can be also calculated with direct referencing to avoid using data transfer links to pipe components in a quick script.

Simulation Setup

We have defined the time constant (the time which the wave will reach the end of the pipe). We will use this to determine the duration of our simulation. We can specify our end time as our time constant (0.055435s).

image 11

Next, if we right click on the inlet node, we can hover over “Plot a Flow Path” and select Total Pressure. This can be plotted for any pathway specified through a network by explicitly specifying a flow path, covered in a previous post. Plotting the flow path will allow us to visualize the pressure profile across the length of the artery. We can solve steady state for our initial conditions and run the solution our network and run our solution.

Simulation Results

We can see a plot of the longitudinal pressure wave in the vessel below

FNX vessel wave

The pressure wave propagation will depend on the impedance mismatch of the downstream vasculature due to changes in diameter or arterial bifurcations.

Lastly, we can us the same procedure to model the wave reflection in the piping system used for Demo Network 4. This network simulated a sudden valve closure on the righthand side of a steel pipe. Using our script, we can determine that the wave will take 0.07525 seconds to reach the other side of the 101m pipe.

image 13

We can visualize the wave reflection in our steel pipe back to the valve by doubling our time constant and running a transient simulation.

FNX pipe wave

A flow path graph may also be used to visualize the placement of an accumulator along the pipe to dampen wave reflection to the valve.

We have seen that Flownex can easily model a wide variety of systems from pulse wave velocity to water hammer effects in piping systems. The pressure wave in the steel pipe(1329m/s) close to the speed of sound in open water, whereas the artery has a pulse wave velocity of 5.4983m/s, far below the speed of sound. Now, we can easily write a quick scripts to determine the speed which pressure waves will pass through a network pathway.

Categories

Get Your Ansys Products & Support from the Engineers who Contribute to this Blog.

Technical Expertise to Enable your Additive Manufacturing Success.

PADT’s Pulse Newsletter

Keep up to date on what is going on at PADT by subscribing to our newsletter.


By submitting this form, you are consenting to receive marketing emails from: . You can revoke your consent to receive emails at any time by using the SafeUnsubscribe® link, found at the bottom of every email. Emails are serviced by Constant Contact

Share this post:

Upcoming Events

09/27/2023

2023 AZ Bio Awards

09/26/2023

Experience Stratasys Truck Tour - Houston

09/22/2023

AIAA Rocky Mountain Section Technical Symposium 2023

09/22/2023

Experience Stratasys Truck Tour - Dallas, TX

09/21/2023

Accelerating the Energy Transition through Simulation

09/20/2023

3D Printing vs. CNC Machining - Webinar

09/13/2023

Maxwell Updates in Ansys 2023 R2 - Webinar

09/12/2023

Sandia Science & Technology Park 25th Anniversary

09/12/2023

Experience Stratasys Truck Tour - Tempe, AZ

09/08/2023

26th Annual New Mexico Flying 40 Awards

09/08/2023

New Mexico Tech Summit

09/07/2023

New Mexico Tech Summit

08/30/2023

Structures Updates in Ansys 2023 R2 (1) - Mechanical, Post & Graphics

08/23/2023

Improved Injection Molding with Additive - Webinar

08/22/2023

SPIE Optics & Photonics Exhibition 2023

08/16/2023

Fluids Updates in Ansys 2023 R2 - Webinar

08/04/2023

Experience Stratasys Truck Tour - Salt Lake City, Utah

08/01/2023

Experience Stratasys Truck Tour - Denver Colorado

07/26/2023

Solving Supply Chain Issues with Additive - Webinar

07/25/2023

Arizona Tech Leadership Golf Tournament

07/24/2023

Arizona Tech CEO Leadership Retreat

07/19/2023

System Automation & Optimization Updates in Ansys 2023 R1 - Webinar

07/13/2023

2023 AEROSPACE, AVIATION, DEFENSE AND MANUFACTURING CONFERENCE

07/12/2023

Materials Updates in Ansys Granta 2023 R1 - Webinar

06/30/2023

Turbo Expo 2023

06/29/2023

Turbo Expo 2023

06/28/2023

Turbo Expo 2023

06/28/2023

Revolutionize Packaging Design with Additive - Webinar

06/27/2023

Turbo Expo 2023

06/27/2023

2023 E-MOBILITY AND CLEAN ENERGY SUMMIT

06/26/2023

Turbo Expo 2023

06/21/2023

Optics Updates in Ansys 2023 R1 - Webinar

06/07/2023

LS-DYNA Updates in Ansys 2023 R1 - Webinar

05/31/2023

Driving Automotive Innovation with Additive - Webinar

05/24/2023

Hill Air Force Base Tech Expo

05/24/2023

Structural Updates in Ansys 2023 R1 (3) – Structural Optimization & Ex

05/23/2023

CROSSTALK 2023: Emerging Opportunities for Advanced Manufacturing Smal

05/10/2023

Signal & Power Integrity Updates in Ansys 2023 R1 - Webinar

04/26/2023

Additive Manufacturing Updates in Ansys 2023 R1 - Webinar

04/20/2023

38th Space Symposium Arizona Space Industry

More Info

04/19/2023

38th Space Symposium
Arizona Space Industry

04/19/2023

Additive Aids for Manufacturing - Webinar

04/18/2023

38th Space Symposium
Arizona Space Industry

04/17/2023

38th Space Symposium

04/13/2023

Venture Madness 2023

04/12/2023

Fluid Meshing & GPU-Solver Updates in Ansys 2023 R1 - Webinar

03/29/2023

8th Thermal and Fluids Engineering Conference

03/29/2023

Structural Updates in Ansys 2023 R1 - Composites, Fracture & MAPDL

03/28/2023

8th Thermal and Fluids Engineering Conference

03/27/2023

8th Thermal and Fluids Engineering Conference

03/26/2023

8TH Thermal and Fluids Engineering Conference

03/24/2023

Arizona BioPreneur Conference | Spring 2023

03/22/2023

2023 Arizona MedTech Conference

03/22/2023

Optimize Jigs & Fixtures with Additive - Webinar

03/15/2023

3D Design Updates in Ansys 2023 R1 - Webinar

03/08/2023

Competitive Advantages of 1D/3D Coupled Simulation - Webinar

03/01/2023

High Frequency Updates in Ansys 2023 R1 - Webinar

02/22/2023

Additive Advantages in Aerospace - Webinar

02/15/2023

Structural Updates in Ansys 2023 R1 (1) - Webinar

02/09/2023

IME 2023: MD&M | WestPack | ATX | D&M | Plastek

02/08/2023

IME 2023 MD&M | WestPack | ATX | D&M | Plastek

02/07/2023

IME 2023 MD&M | WestPack | ATX | D&M | Plastek

01/27/2023

Arizona Photonics Days, 2023

01/26/2023

Arizona Photonics Days, 2023

01/26/2023

TIPE 3D Printing | 2023

01/26/2023

Venture Cafe Phoenix Talent Night - Job Fari

01/26/2023

VFS 2023 Autonomous/Electric VTOL Symposium

01/25/2023

Arizona Photonics Days, 2023

01/25/2023

Building A.M.- Utah: Kickoff!

01/25/2023

TIPE 3D Printing | 2023

01/25/2023

VFS 2023 Autonomous/Electric VTOL Symposium

01/24/2023

VFS 2023 Autonomous/Electric VTOL Symposium

01/24/2023

TIPE 3D Printing | 2023

01/18/2023

2023 AZ Tech Council Golf Tournament

12/21/2022

Simulation Best Practices for 5G Technology - Webinar

12/14/2022

Digital Twins Updates in Ansys 2022 R2 - Webinar

12/08/2022

Tech the Halls - AZ Tech Council Holiday Mixer

12/07/2022

Electric Vehicle and Other Infrastructure Update Panel

11/30/2022

SPEOS Updates in Ansys 2022 R2 - Webinar

11/23/2022

Simulation Best Practices for Electronics Reliability - Webinar

11/16/2022

Discovery Updates in Ansys 2022 R2

11/10/2022

VentureCafe Phoenix Panel: Venture Capital in AZ

11/08/2022

2022 GOVERNOR’S CELEBRATION OF INNOVATION AWARDS + TECH SHOWCASE

11/03/2022

VentureCafe Phoenix Panel: Angel Investment in AZ

11/02/2022

High & Low Frequency Electromagnetics Updates in Ansys 2022 R2

10/26/2022

Simulation Best Practices For Chip-Package-System Design & Development

10/20/2022

Nerdtoberfest 2022

10/19/2022

2022 Southern Arizona Tech + Business Expo

10/19/2022

LS-DYNA Updates in Ansys 2022 R2 - Webinar

10/17/2022

Experience Stratasys Truck Tour - Clearfield Utah

10/14/2022

ASU School of Manufacturing Systems and Networks - Formal Opening Cele

10/14/2022

Experience Stratasys Truck Tour - Midvale Utah

10/12/2022

Experience Stratasys Truck Tour - Littleton Colorado

10/06/2022

Fluids Updates in Ansys 2022 R2 - Webinar

10/05/2022

Experience Stratasys Truck Tour - Colorado Springs

09/29/2022

White Hat Life Science Investor Conference - 2022

09/28/2022

2022 AZBio Awards

09/28/2022

Simulation Best Practices for Rotating Machinery Design & Development

09/21/2022

ExperienceIT NM 2022

Search in PADT site

Contact Us

Most of our customers receive their support over the phone or via email. Customers who are close by can also set up a face-to-face appointment with one of our engineers.

For most locations, simply contact us: