July 3, 2024

158 - Fire Fundamentals pt. 9 - Know you boundaries (in CFD)

158 - Fire Fundamentals pt. 9 - Know you boundaries (in CFD)

In today's Fire Science Show, we talk about how boundary conditions can make or break your fire simulation models. We'll explore boundary conditions' fundamental role in defining how simulations interact with their environments and how mastering these can lead to more accurate and reliable fire simulation models. I hope we break down some complex topics into manageable insights. Also, I hope we've turned something really boring into an interesting and fun episode. We discuss:

1. Boundaries interacting with flow:

  • pressure inlets/outlets [open];
  • velocity boundary conditions [vents];
  • mass flow inlets;
  • fans and HVAC models.

2. Boundaries containing the flow - walls:

  • heat transfer at walls;
  • shear at walls.

And other minor boundaries such as symmetry and interior. The episode is largely based on my expertise with ANSYS Fluent, but I've tried to make it relatable to FDS as well, ensuring that the content is practical and valuable for your work in fire science and simulation modeling.

----
The Fire Science Show is produced by the Fire Science Media in collaboration with OFR Consultants. Thank you to the podcast sponsor for their continuous support towards our mission.

Chapters

00:00 - Understanding CFD Boundary Conditions

10:52 - Understanding CFD Velocity Boundaries

23:26 - Modeling Flow and Solid Boundaries

38:27 - Solid and Fluid Boundary Conditions

Transcript
WEBVTT

00:00:00.221 --> 00:00:02.126
Hello everybody, welcome to the Fire Science Show.

00:00:02.126 --> 00:00:09.252
This is Fire Science Fundamentals, episode 9, and I know that you'll really enjoy this series.

00:00:09.252 --> 00:00:16.390
I know it because the first episode of Fire Fundamentals, with Rory Haddon on ignition has just reached 2000 downloads.

00:00:16.390 --> 00:00:18.353
That is incredible.

00:00:18.353 --> 00:00:22.109
You, the audience, were paramount in coming up with this series.

00:00:22.208 --> 00:00:30.653
I was asked over and over again to do some basic physics episodes in a kind of approachable way and it kind of took off and I really enjoyed doing it.

00:00:30.653 --> 00:00:33.588
So let's continue doing fire fundamentals.

00:00:33.588 --> 00:00:46.790
The last time I did it solo was on compartment fires and previously, before that, we talked about CFD, and the CFD episode was one that many people really enjoyed and yeah, unfortunately for me, you have asked for more.

00:00:46.790 --> 00:00:50.130
I was not really sure what's going to happen after that episode.

00:00:50.130 --> 00:01:00.088
I was pretty sure that I'm going to get a lot of rage emails that I am not covering fluid dynamics in a way that I should, not in a scientific way.

00:01:00.088 --> 00:01:01.771
But actually people have enjoyed that.

00:01:01.771 --> 00:01:03.134
There were no complaints.

00:01:03.134 --> 00:01:04.585
There were questions for more.

00:01:04.585 --> 00:01:05.701
So here we are.

00:01:05.701 --> 00:01:24.631
I'm going to do another episode on the fundamentals of CFD modeling and today we will be discussing boundary conditions and, as boring as it sounds, I'll try to make it exciting and interesting for you that even such a simple thing can be taught in a way that is fun and that really highlights the importance of the subject.

00:01:24.631 --> 00:01:32.200
For the goodness of your models, it's critical to know your boundary conditions, so let's spin the intro and learn our boundaries.

00:01:37.503 --> 00:01:39.144
Welcome to the Firesize Show.

00:01:39.144 --> 00:01:42.605
My name is Wojciech Wigrzyński and I will be your host.

00:01:42.605 --> 00:02:00.656
This podcast is brought to you in collaboration with OFR Consultants.

00:02:00.656 --> 00:02:03.597
Ofr is the UK's leading fire risk consultancy.

00:02:03.597 --> 00:02:14.448
Its globally established team has developed a reputation for preeminent fire engineering expertise, with colleagues working across the world to help protect people, property and environment.

00:02:14.448 --> 00:02:30.245
Established in the UK in 2016 as a startup business of two highly experienced fire engineering consultants, the business has grown phenomenally in just seven years, with offices across the country in seven locations, from Edinburgh to Bath, and now employing more than a hundred professionals.

00:02:30.245 --> 00:02:41.905
Colleagues are on a mission to continually explore the challenges that fire creates for clients and society, applying the best research experience and diligence for effective, tailored fire safety solutions.

00:02:41.905 --> 00:02:52.550
In 2024, ofr will grow its team once more and is always keen to hear from industry professionals who would like to collaborate on fire safety futures.

00:02:52.550 --> 00:02:55.865
This year, get in touch at ofrconsultantscom.

00:02:57.168 --> 00:03:06.520
I wonder if you have this feeling that you wake up and you realize that what you have planned is perhaps a really bad idea and you should not be pursuing that.

00:03:06.520 --> 00:03:15.467
I go through that every second podcast episode, but I'm very bad at listening to my voice of consciousness, or perhaps my imposter.

00:03:15.467 --> 00:03:19.760
And here I am talking on podcasts about boundary conditions.

00:03:19.760 --> 00:03:34.419
You could claim perhaps one of most boring topics one could pick, but I really see some funny and interesting things related to boundary conditions, on a very fundamental level of their use actually, which can make or break your simulation.

00:03:34.419 --> 00:03:48.275
And actually every now and then I see some people posting on LinkedIn their simulations and you immediately realize that their boundaries are not really the ones they should use or they're perhaps misused in a way.

00:03:48.560 --> 00:03:53.733
So let's try and talk about boundary conditions in fluid dynamics, in simulations.

00:03:53.733 --> 00:03:59.372
And note, before we start, you have to understand that I am ANSYS fluent person.

00:03:59.372 --> 00:04:00.765
I am an ANSYS fluent user.

00:04:00.765 --> 00:04:08.663
I'm trained and, let's say, skilled at FDS, but it's not my daily software that I would be using for everything.

00:04:08.663 --> 00:04:16.273
So if I confuse the name of a boundary condition in FDS, please feel free to correct me and I'll try to fix it.

00:04:16.273 --> 00:04:17.377
I'm learning every day.

00:04:17.377 --> 00:04:21.769
Hopefully I'll get my boundaries done in Ansys and correct If I mess that one.

00:04:21.769 --> 00:04:22.711
That's going to be hilarious.

00:04:23.199 --> 00:04:26.009
Anyway, if we talk about boundaries, what are boundary conditions?

00:04:26.009 --> 00:04:28.168
Why do we need them in our simulations?

00:04:28.168 --> 00:04:35.005
If you would like to attempt a fluid, dynamic simulation of an endless space, let's think planet Earth.

00:04:35.005 --> 00:04:41.466
You swim to the east, you eventually reach the point from where you've started An endless continuity.

00:04:41.466 --> 00:04:44.279
If you look down, well, there's a boundary condition.

00:04:44.279 --> 00:04:46.086
There's a solid surface or an ocean.

00:04:46.086 --> 00:04:51.002
If you look up, it ends in with a vacuum of space, also something you could call a boundary condition.

00:04:51.002 --> 00:04:55.982
So even in an endless space of Earth, you are still constrained by something.

00:04:55.982 --> 00:05:04.687
And obviously, when we simulate our fire cases, it's not that we're going to simulate an endless atmosphere around our buildings.

00:05:04.687 --> 00:05:08.274
I mean time is money and resources are everything.

00:05:08.274 --> 00:05:11.187
They dictate what we can and what we cannot do.

00:05:11.329 --> 00:05:34.089
In our simulations we have to narrow the field of what we're investigating into a part of the space and, as you remember from the previous CFD episode of the Far Fundamentals, cfd works basically on solving multiple equations momentum transfer, heat transfer, species conservation and so on, so on in a confined space.

00:05:34.089 --> 00:05:42.170
And for it to work, it must understand what's happening when it reaches the end of that space.

00:05:42.170 --> 00:05:44.206
Does it magically disappear?

00:05:44.206 --> 00:05:47.387
Does it bounce back from a solid obstacle?

00:05:47.387 --> 00:05:48.591
These are your boundaries.

00:05:48.591 --> 00:05:58.632
These are the elements that constrain your fluid, that constrain your model and through which you can actually act and define your model.

00:05:58.632 --> 00:06:06.761
These are critical parts of your model where you have the ability to interact with the fluid that you are modeling.

00:06:06.761 --> 00:06:18.072
Actually, if you think about fluid itself, you can define fluid as a substance that continuously deforms under the application of sheer stress and adapts its shape to the constraints.

00:06:18.072 --> 00:06:19.221
So here you are.

00:06:19.221 --> 00:06:22.550
That's your boundary conditions, things that act on your fluid.

00:06:22.750 --> 00:06:29.192
Now, if we want to narrow the discussion into particular types of boundary conditions I'm not going to list them, that's pointless.

00:06:29.192 --> 00:06:34.007
I would say I would define two major categories of the boundaries.

00:06:34.007 --> 00:06:43.012
One would be the boundaries that actually interact with your fluid, so I would call them the flow boundaries or the fluid boundaries.

00:06:43.012 --> 00:06:48.651
It's a very unprofessional way to call them, but it kind of delivers the message.

00:06:48.651 --> 00:06:55.952
It's the boundaries that interact with your fluid and act on it, and boundaries that are pretty much solid walls, like solid boundaries of your model.

00:06:55.952 --> 00:06:59.105
And there are some special boundaries which we'll talk in the end.

00:06:59.105 --> 00:07:04.605
So stay tuned for a surprise hopefully interesting For fluid boundaries.

00:07:05.307 --> 00:07:13.701
Into this big bucket I put all the boundaries that you use as a CFD user to interact with the fluid you have in your model.

00:07:13.701 --> 00:07:17.987
That most likely is air or some mixture of air and smoke.

00:07:17.987 --> 00:07:35.069
I often call the smoke spoiled air because in CFD modeling the smoke properties are actually quite same as the air has, and the only difference being that there are species in it, like soot and other products of combustion, and due to that it has different emissivity.

00:07:35.069 --> 00:07:41.408
However, for most cases the smoke and air are pretty much the same thing from the fluid dynamics perspective.

00:07:41.408 --> 00:07:47.108
Anyway, from my perspective, the first, most basic boundary condition I would always have in my model.

00:07:47.108 --> 00:07:51.584
It's hard for me to figure out a single model in which I would never use that boundary condition.

00:07:52.125 --> 00:07:58.461
These are conditions that simulate an open end to your numerical domain.

00:07:58.461 --> 00:08:02.805
In ANSYS we call them pressure boundaries, pressure inlets, pressure outlets.

00:08:02.805 --> 00:08:05.629
We call them pressure boundaries, pressure inlets, pressure outlets.

00:08:05.629 --> 00:08:12.196
In FTS I think we commonly refer to them open, and those are types of vents that are used in FTS.

00:08:12.196 --> 00:08:27.480
Regardless of the name, the point of this open or pressure boundary is to define space through which the fluid can flow freely and there's a defined pressure dynamic pressure on that boundary condition.

00:08:27.480 --> 00:08:28.983
Now, what does it mean?

00:08:28.983 --> 00:08:38.100
If the dynamic pressure on that boundary is positive, then the flow is coming from that boundary into your model, unless the pressure in your model is higher than that one.

00:08:38.100 --> 00:08:41.610
If the pressure is negative, there's going to be a suction effect.

00:08:41.610 --> 00:08:47.684
If it's zero, then the flow will be an outcome of the pressure value inside your domain.

00:08:47.684 --> 00:08:48.787
Pretty simple, isn't it?

00:08:48.787 --> 00:08:58.062
But there are some issues that you can easily incorporate into your model by simply applying these boundary conditions without thinking.

00:08:58.863 --> 00:09:06.541
If you think about pressure in your domain, well, it's not that the pressure in every point of your domain is exactly same.

00:09:06.541 --> 00:09:08.585
There's hydrostatic pressure.

00:09:08.585 --> 00:09:13.313
That is simply an outcome of your fluid having some weight.

00:09:13.313 --> 00:09:17.068
It just weights and you'll have a hydrostatic pressure buildup.

00:09:17.068 --> 00:09:20.625
In normal conditions that would be roughly 12 pascals per meter of height.

00:09:20.625 --> 00:09:22.610
So that's quite a significant pressure.

00:09:22.610 --> 00:09:32.562
Now imagine, on your other side of your boundary condition you always have zero and on the domain in which you are simulating your fluid flow, the pressure is growing.

00:09:32.562 --> 00:09:46.354
With every meter down, it simply increases, which means that you could end up with artificial pressure difference between your domain and exterior just because you've defined the pressure on the other side as zero.

00:09:46.354 --> 00:09:51.671
For the CFD software, along the height of your pressure boundary, the pressure is zero.

00:09:51.671 --> 00:09:53.062
It's not hydrostatic.

00:09:53.062 --> 00:09:59.923
Perhaps in other softwares than ANSYS it's predefined as hydrostatic, being included automatically.

00:09:59.923 --> 00:10:04.802
For us, there's actually a necessity to consider that and include that in our pressure boundaries.

00:10:04.802 --> 00:10:08.631
So here, even the most simple condition has a challenge for you.

00:10:09.139 --> 00:10:17.010
Now another thing that's interesting in pressure boundaries or open boundaries, sometimes we would be only interested in modeling the interior of our building.

00:10:17.010 --> 00:10:24.211
Let's say you're modeling a warehouse and you want to put some sort of inlets, outlets, to that warehouse.

00:10:24.211 --> 00:10:31.011
You're working with a natural smoke control, so you're obviously using things like natural smoke extraction ventilators.

00:10:31.011 --> 00:10:34.971
You're using doors as your makeup air sources.

00:10:34.971 --> 00:10:38.390
We've talked about that in previous Fire Fundamentals episodes.

00:10:38.390 --> 00:10:45.250
Now you would like to model them, and there is an urge that the simplest way to model them was a natural opening.

00:10:45.250 --> 00:10:49.054
That's like literally something that connects your interior to the exterior, right?

00:10:49.054 --> 00:10:52.504
So perhaps a pressure boundary is a perfect way to simulate that.

00:10:52.904 --> 00:11:08.282
Now the thing is that when you have natural ventilators, especially smaller or ones that do not open fully, there are different types of natural ventilators that are less efficient In those devices.

00:11:08.282 --> 00:11:13.113
It's not that the air flows freely through the entire cross-section of the device.

00:11:13.113 --> 00:11:15.368
It may be blocking the pathway for the air.

00:11:15.368 --> 00:11:18.669
It will definitely include some resistance for air to go.

00:11:18.669 --> 00:11:40.591
So if you have, let's say, a square meter of a natural ventilator on your roof and in that space you put a square meter of a natural ventilator on your roof and in that space you put a square meter of a pressure boundary, the effects of pressure boundary will be more pronounced, more significant than the effects of your natural ventilator, because the natural ventilator is less efficient than an idealized hole.

00:11:40.591 --> 00:11:43.402
Another thing outside you would have different temperatures.

00:11:43.402 --> 00:11:51.645
So you also have to define this temperature on the other side of your open boundary to really capture the chimney effects and all the buoyancy effects that will come in play.

00:11:51.966 --> 00:11:57.969
So just dropping your open boundaries, because they are not openings, they're pressure boundaries.

00:11:57.969 --> 00:12:03.028
You have to understand how you use them and be conscious in where you place them.

00:12:03.028 --> 00:12:06.662
If you want to simulate natural ventilators, they're not the way to go.

00:12:06.662 --> 00:12:15.898
You need to simulate some sort of a roof structure and some sort of openings in that roof that at least try to mimic the effectiveness of natural ventilators.

00:12:15.898 --> 00:12:23.524
Perhaps you can go away in the most simplistic case by using smaller, but I'm not going to give you a recommendation of how small you should make them.

00:12:23.524 --> 00:12:30.583
If you want to play like that, you need to do your own validation and figure out values that would work for yourself.

00:12:30.583 --> 00:12:34.775
Now we've briefly talked about the pressure boundaries.

00:12:35.378 --> 00:12:41.350
Another boundary that I would very, very commonly use in my models would be velocity boundaries.

00:12:41.350 --> 00:12:44.764
So here we are talking about mechanical vents.

00:12:44.764 --> 00:12:47.581
In ANSYS we would call them velocity inlets.

00:12:47.581 --> 00:13:00.994
In this case, instead of a pressure on my boundary, which again is something that fluid can flow through, I define velocity of the flow and I define a vector which way the flow goes.

00:13:00.994 --> 00:13:12.869
So I am capable of defining a fixed suction or inlet point in my model through which air will come or be extracted through with a specified velocity.

00:13:12.869 --> 00:13:15.701
I also define stuff like temperatures.

00:13:15.701 --> 00:13:22.082
I also define in Ansys stuff like turbulent intensity and other important parameters for my simulation.

00:13:22.504 --> 00:13:26.461
But the most important thing, what's the velocity at the boundary condition?

00:13:26.461 --> 00:13:29.327
Now, again, we're talking about boundary condition.

00:13:29.327 --> 00:13:34.649
You have a patch of your model that this boundary condition is applied to.

00:13:34.649 --> 00:13:57.472
When I say that the boundary condition of velocity is three meters per second, every single cell, that's a part of this patch that I called my velocity inlet now will have three meters of velocity introduced into my model, exactly this value, and the model will make sure that this value is the value introduced to my model.

00:13:57.980 --> 00:13:59.626
Now there's an issue with that.

00:13:59.626 --> 00:14:06.010
It's pretty robust and simple and in most cases it works without playing too much with it.

00:14:06.010 --> 00:14:19.725
But if you would like to be really, really precise in modeling your vents, it's not that the entire ventilator has the same velocity at its surface, whatever the surface is.

00:14:19.725 --> 00:14:22.509
It's actually a virtual concept the surface of a ventilator.

00:14:22.509 --> 00:14:41.532
Anyway, through the cross-section from which the air enters your domain, in a real case, in a real building, you will have some sort of velocity profile, whereas in your CFD model you will have idealized averaged flow over the entirety of the velocity patch that you put into your model.

00:14:42.261 --> 00:14:43.419
How big is the issue?

00:14:43.419 --> 00:14:59.850
So if you have a very narrow duct work, you have we call them pancake ducts very flat ducts in which one dimension is pretty long and the other is pretty narrow and you have a very big vent on such a pancake duct.

00:14:59.850 --> 00:15:05.130
You usually end up with extremely uneven flow through that duct.

00:15:05.130 --> 00:15:19.552
So you could have 2 m per second average on such an opening, which could sound acceptable in some sort of cases for smoke control, whereas in reality you would have no flow on most of the opening and 5 m at the bottom.

00:15:19.552 --> 00:15:32.654
Now those 5 m will cause you trouble with your smoke control system because that's a lot of kinetic energy introduced into your building and a good chance that this flow will mix something in a place that you don't want it to be mixed.

00:15:32.654 --> 00:15:56.620
And especially if you're modeling stuff like corridors or some narrow, tight spaces, this error can occur, lead to a system which works perfectly in simulations by applying this averaged velocity over your inlet, whereas in reality the velocity is skewed and you get much higher velocities and the system actually does not work.

00:15:56.620 --> 00:16:00.028
We've seen that happen and that's a pretty big challenge.

00:16:00.934 --> 00:16:19.043
The ways to overcome that we usually try to model the ductwork and so the the profile on the inlet is actually some sort of an outcome of the calculations of the model, and my velocity inlet is actually where my extraction or inlet fans are physically.

00:16:19.043 --> 00:16:26.917
So instead of simulating just grills that are in my compartment, I would would simulate the ductwork For the more challenging projects.

00:16:26.917 --> 00:16:30.376
I wouldn't say that we did it every single time, but when we need to pay attention.

00:16:30.376 --> 00:16:31.299
We do it like that.

00:16:31.299 --> 00:16:34.572
Another way you can define profiles.

00:16:34.572 --> 00:16:37.519
This is a very interesting way to work with those boundary conditions.

00:16:37.519 --> 00:16:48.441
So you can actually force CFD software to apply a specific profile to the velocity introduced so it's no longer even but actually conforms to however you define it.

00:16:48.441 --> 00:16:50.565
So yeah that, that that's the way.

00:16:51.226 --> 00:17:07.556
Another thing interesting with velocity boundaries is that when you think you introduce flow from a mechanical ventilator into your compartment, you would most likely think about your horizontal velocity component from that fan.

00:17:07.556 --> 00:17:11.328
Right, because it's just blowing air straight in front right.

00:17:11.328 --> 00:17:22.018
But in fact, because the fan is a rotating machine, it actually introduces quite a significant tangential component to velocity.

00:17:22.018 --> 00:17:31.403
It swirls and in some cases, in some aspects of modeling, this swirl may be actually quite important to the outcomes of your simulation.

00:17:31.403 --> 00:17:31.875
Here.

00:17:31.875 --> 00:17:34.765
I wouldn't say it's something we would use commonly.

00:17:34.765 --> 00:17:37.022
I would say it's very rarely used.

00:17:37.022 --> 00:17:48.105
But especially in tunneling projects when you're modeling big jet fans, you're not usually doing that with velocity inlets, but other flow conditions also can apply as well.

00:17:48.105 --> 00:17:56.343
This could actually make a big difference in the outcomes of your modeling of the flow profile of your jet fan.

00:17:56.343 --> 00:18:05.780
So yeah, sometimes you have to go in your understanding beyond just the basic characteristic that you're pushing air forward.

00:18:05.780 --> 00:18:08.135
I know some cases.

00:18:08.135 --> 00:18:12.965
There's a gentleman, fahd Itrada, who is doing amazing modeling of chat fans.

00:18:12.965 --> 00:18:27.598
They've developed a fan called Mojet and truly the modeling they do to model their fans is remarkable and the level of details they go into to actually capture the behavior of that chat fan is stunning.

00:18:27.598 --> 00:18:36.746
So, as with most of the things you do in fire safety science and fire safety engineering, when you go deep inside it becomes pretty complicated.

00:18:37.326 --> 00:18:42.503
Another aspect of velocity boundary conditions is when we try to apply wind.

00:18:42.503 --> 00:18:43.325
That's a big trap.

00:18:43.325 --> 00:18:44.748
A lot of people fall into that.

00:18:44.748 --> 00:18:50.324
So wind is not just a boundary condition that blows air at a constant speed.

00:18:50.324 --> 00:18:58.366
Wind is a complicated phenomenon, physical phenomenon, and you really have to introduce wind first with a profile.

00:18:58.366 --> 00:19:00.623
So it has to follow some profile.

00:19:00.623 --> 00:19:06.942
The wind at the ground level will be weaker than the wind at some height.

00:19:06.942 --> 00:19:11.228
There are logarithmic profiles that define the velocity of wind.

00:19:11.228 --> 00:19:15.500
With the height, wind will introduce turbulence, so you have to capture the turbulence.

00:19:15.500 --> 00:19:32.333
If you're modeling it with large eddy simulation, you may be even in need to get some sort of periodic boundary condition that increases, decreases velocity continuously so you capture the formation of large eddies or large vortices that come with winds.

00:19:32.333 --> 00:19:34.317
So you get the spikes in velocity correct.

00:19:34.317 --> 00:19:43.300
It's actually quite challenging, but just dropping a massive velocity inlet boundary condition in front of your model, it it's not really cutting.

00:19:43.300 --> 00:19:55.724
That that's not wind, that's a big chunk of velocity, not a physical phenomenon of wind, which is beautiful and complex, and I have many other podcast episodes that go deeper into that.

00:19:55.924 --> 00:19:57.127
There's another thing about velocity.

00:19:57.127 --> 00:19:58.660
I mean velocity inlets.

00:19:58.660 --> 00:20:06.981
This is very simple but it could be a very problematic boundary condition or it can have a very far-reaching consequences for your models.

00:20:06.981 --> 00:20:25.240
So if you define the velocity inlet boundary condition, you define the speed at which the air is put into your model or extracted from your model, which basically defines the exact value of volumetric flow that's flowing through that space, and this will be constrained by the model.

00:20:25.240 --> 00:20:28.460
Now imagine you have your model.

00:20:28.460 --> 00:20:30.727
It's a cube, it's a compartment.

00:20:30.727 --> 00:20:33.935
You have one velocity inlet that brings air into the room.

00:20:33.935 --> 00:20:37.304
You have velocity outlet that extract the air from the room.

00:20:37.304 --> 00:20:54.821
You've defined both velocities to match each other and now you introduce the fire into the room and what will happen is that you will observe a pressure increase into the room, perhaps even up to some unphysical levels of pressure and your CFD model will crash.

00:20:54.821 --> 00:20:57.242
It will give you very odd results.

00:20:57.242 --> 00:20:57.564
Why?

00:20:57.564 --> 00:21:04.923
Because the system is not able to blow more air through the velocity inlet, which has a defined velocity.

00:21:04.923 --> 00:21:10.484
It has to extract or input exact the number you've put into it.

00:21:10.484 --> 00:21:18.665
If there's a pressure increase around that boundary condition, it's not changing the outcomes, so the pressure increases, decreases.

00:21:18.705 --> 00:21:31.944
You can get into big trouble in your models if you only define very fixed values on those types of boundary conditions, and that's why I said I use pressure boundaries in almost every single model.

00:21:31.944 --> 00:21:42.684
You almost always need a pressure boundary somewhere to improve the convergence of your simulation, to have a way to manage the decrease, increase of pressures.

00:21:42.684 --> 00:21:46.598
Of course this has to represent some physical leakage in your building.

00:21:46.598 --> 00:21:57.669
If you're designing a, the most sealed building in the planet, you perhaps should do velocity boundaries, but then again the fans don't work like that and they would adjust.

00:21:57.669 --> 00:21:59.398
Oh yeah, let's go to fans.

00:21:59.901 --> 00:22:06.281
Actually, that's accidentally a very good segue to the next boundary condition I wanted to talk about which, which is fan boundary condition.

00:22:06.281 --> 00:22:13.240
And this is a boundary condition in which we have a slice of space through which the air can flow through.

00:22:13.240 --> 00:22:20.637
So it's a boundary condition that does not mark the outer regions of our simulation.

00:22:20.637 --> 00:22:24.942
It's placed inside the fluid and it accelerates the fluid.

00:22:24.942 --> 00:22:28.164
It acts as a fan, as a ventilator would.

00:22:28.164 --> 00:22:33.990
Now there are so many ways you can design this boundary condition.

00:22:33.990 --> 00:22:51.103
Actually, my friends from Silesian University of Technology, alexander and Gosia Kurl, have wrote a really nice paper on all the different ways you can use to model fan boundary condition to capture the behavior of a fan, especially jet fans, in your simulation better.

00:22:51.103 --> 00:22:55.840
And yeah, I really like to use this, this boundary condition in ansys.

00:22:55.840 --> 00:23:08.637
It gives me very good, uh, first, it gives me very good control over the flow of my fan and secondly, because it's based on decades of answers being used to model fans for fan manufacturers.

00:23:08.637 --> 00:23:15.057
They actually got very convenient and good and accurate models of fans built within the software.

00:23:15.057 --> 00:23:26.190
So for me it's just convenient and if you think about it, it's kind of funny because it's a very simple answer to a problem that have been quite complicated in the past.

00:23:26.692 --> 00:23:41.195
So if I remember correctly, the first times we tried to do jet fans we were mostly using velocity inlets and you know, having air sucked into the velocity inlet on one end of the jet fan, having it released on the other end.

00:23:41.195 --> 00:23:48.240
The issue was it was a different air because it kind of disappeared when it reached the other end of the boundary condition.

00:23:48.240 --> 00:23:49.684
It was not modeled anymore.

00:23:49.684 --> 00:24:03.243
Then we went into recirculation inlets, outlets, which were velocity inlets connected with, let's say, a magical duct that connects one end to another and that allowed us to transfer air from one side to another.

00:24:03.243 --> 00:24:13.537
But this did not allow us to model the behavior of the ventilator where it receives different pressures on the inlet or the losses change on the outlet side.

00:24:13.958 --> 00:24:15.942
The fan, of course, is a mechanical device.

00:24:15.942 --> 00:24:19.558
It has its power curve, it has its fan curve.

00:24:19.558 --> 00:24:24.999
Those are the things you can actually model with a good fan model and by introducing that.

00:24:24.999 --> 00:24:40.128
It takes a lot of work actually to introduce that correctly, but when you introduce that you are actually modeling a real ventilation device in your system and I think in some projects it's absolutely necessary, just like capturing the profiles on your velocity inlets.

00:24:40.128 --> 00:24:58.175
This is sometimes also critical to get your smoke control system modeled in a way that you can be fairly sure that what you've done is accurate and in the building, when they finally build the building and they finally do the smoke control system in it, it's going to work like you've simulated.

00:24:58.175 --> 00:25:07.317
It's a very stressful situation when reality validates or verifies your numerical modeling and you really want to be on the safe side.

00:25:07.317 --> 00:25:13.040
And you get on the safe side by knowing your boundaries, knowing your boundary conditions and applying them correctly.

00:25:13.842 --> 00:25:22.146
An interesting thing for modeling fans and, in general, smoke control systems or ventilation systems is the hvc model in fds.

00:25:22.146 --> 00:25:28.459
Massive shout out to jason floyd who was spearheading the efforts to write this model.

00:25:28.459 --> 00:25:44.876
I think it's a brilliant addition to FDS and it's simply a nodal model where you can basically build a network of nodes that interconnect to each other and simulate the ductworks, but instead of physically simulating the ducts, you know the shape and the flows inside.

00:25:44.876 --> 00:26:01.624
There are physical relations that describe the flow of air within the ducts, through ventilation devices into your rooms and so on, and this network is modeled along your CFD model and the boundary conditions are defined by that network.

00:26:01.624 --> 00:26:11.642
So it takes a lot of hard work from you onto the network and it actually calculates what are the flows, what are the pressure points on your fans and so on.

00:26:11.642 --> 00:26:16.128
A brilliant way to actually make the lives of fire engineers simpler.

00:26:16.128 --> 00:26:19.022
So highly appreciate that this was included in FDS.

00:26:19.022 --> 00:26:26.898
I actually wish we had something like that in Ansys, like if we wanted something like this, we would have to program it ourselves.

00:26:26.898 --> 00:26:32.555
And uh, yeah, sometimes we we do some attempts not as fancy as the solver in fds.

00:26:32.555 --> 00:26:35.041
I'm highly envious in this regard.

00:26:35.041 --> 00:26:38.090
Perhaps we need to recruit jason to write one for us.

00:26:38.792 --> 00:26:41.804
The hvc nodes are also something you can use to model your jet fan.

00:26:41.804 --> 00:26:52.444
So, as I said previously, we have the fan boundary condition in answers, which allows us to define the performance of the fan based on what comes in, what comes out, what's the losses and so on.

00:26:52.444 --> 00:26:58.281
You could actually do the same thing with HVAC nodes in FDS to efficiently model your jet fans.

00:26:58.281 --> 00:27:02.577
I'm not sure if you can add the tangential swill to your flow in FDS.

00:27:02.577 --> 00:27:06.567
Probably you can, but please don't quote me on that.

00:27:06.567 --> 00:27:11.244
You have to submit an issue on issue tracker and figure out if you can or not.

00:27:11.244 --> 00:27:12.816
If someone knows, let me know.

00:27:12.816 --> 00:27:16.730
I'm actually curious now because I have not fact-checked it before.

00:27:16.730 --> 00:27:19.840
I've never done it in this detail in FDS.

00:27:20.382 --> 00:27:24.714
Anyway, this would be a quick summary of the important flow boundaries.

00:27:24.714 --> 00:27:25.799
No wait, there's one more.

00:27:25.799 --> 00:27:29.240
There's one more very important flow boundary, that's mass flow inlets.

00:27:29.240 --> 00:27:35.488
So again, very similar to velocity inlets, very similar to your mechanical vents.

00:27:35.488 --> 00:27:43.522
But in this case, instead of prescribing a specific velocity or volumetric flow rate, you would specify the mass flow rate.

00:27:43.522 --> 00:27:49.703
Now you could think it's the same thing and yes, that's true if you have your flow in your ambient air.

00:27:49.703 --> 00:27:59.486
So if you're talking about delivering air to your model at the same temperature over the course of your simulation, yes, that actually is true that it's the same boundary condition.

00:27:59.968 --> 00:28:07.957
But if you define extraction in a form of mass flow, then you would get completely different outcomes as the fire grows.

00:28:07.957 --> 00:28:16.520
The reason for that is that as fire grows, the air heats around, and as the air is heated it changes its density.

00:28:16.520 --> 00:28:26.306
So if you say you're always removing five kilograms of air from your model, in a cold situation, that would be like four something ish cubic meters.

00:28:26.306 --> 00:28:30.201
But if the air is at 600 degrees, that would be twice the number.

00:28:30.201 --> 00:28:34.417
So you have to be careful with those boundary conditions.

00:28:34.478 --> 00:28:45.367
Actually, we in the past we found a very clever way to use mass flow boundary conditions to define a completely new way of extracting smoke from buildings, extremely efficient.

00:28:45.367 --> 00:28:47.481
We coined it as a smart smoke control.

00:28:47.481 --> 00:28:49.019
There are papers on that.

00:28:49.019 --> 00:28:53.142
If you wish to dig in deeper, perhaps one day I'll do an episode on that.

00:28:53.142 --> 00:28:56.065
That was my pet peeve in 2017.

00:28:56.065 --> 00:29:00.143
My biggest discovery you really deserve to learn that.

00:29:00.143 --> 00:29:02.041
Perhaps that's a topic for episode 200.

00:29:02.041 --> 00:29:02.423
We'll see.

00:29:03.636 --> 00:29:12.868
Anyway, there are clever uses for mass flow inlets and also one common use of mass flow inlet boundary conditions is related to your burners, to the fire source.

00:29:12.868 --> 00:29:18.888
So we would very often like to define the inlet of fuel to our model with a mass flow inlet.

00:29:18.888 --> 00:29:19.990
Why?

00:29:19.990 --> 00:29:25.363
Because it translates very easily to the heat release rate on paper.

00:29:25.363 --> 00:29:31.278
You just multiply by calorific value and efficiency of combustion and you're there, you got your heat release rate.

00:29:31.278 --> 00:29:35.201
So it's very easy to control it in that way.

00:29:35.201 --> 00:29:45.057
I think when you do burners in FTS and you define the amount of heat release, it actually calculates it down to the mass and then it's a mass flow inlet boundary't.

00:29:45.057 --> 00:29:45.656
Quote me on that.

00:29:45.656 --> 00:29:54.619
I'm not 100 sure, but the outcome is like that you you bringa specific value amount of mass into your model in every second of your simulation.

00:29:55.161 --> 00:29:56.894
Are there any challenges with that?

00:29:56.894 --> 00:29:57.678
Yeah, of course there are.

00:29:57.678 --> 00:29:59.606
With everything there's challenges.

00:29:59.606 --> 00:30:16.759
In fire safety engineering we've actually found one very interesting challenge If you're having very small burners and very small heat release rates, and that's the case when you have modeling of scaled down experiments in fire science.

00:30:16.759 --> 00:30:20.172
So we have this technique called the Froude number scaling.

00:30:20.172 --> 00:30:30.042
Previous episode of fire science show was very rough and here's another pebble to the pile saying that it's not the best technique to study fire science.

00:30:30.042 --> 00:30:58.430
So if you have a very small mass flow rates and the fuel is released from the surface with very low velocity, the efficiency of combustion actually changes, and my student, jakub, is working very hard on a paper that summarizes our findings in that regard and once this sees the daylight, I will be first to to tell you about it, because I'm very proud of this work and, yeah it, it highly relates to the mass flow inlet boundary condition.

00:30:58.430 --> 00:31:03.941
So, yeah, that was uh, my good summary of uh, boundary related to flow.

00:31:04.342 --> 00:31:14.039
But there's a second group boundary conditions that I would call the solid boundary conditions, or the physical constraints of your model, basically, the walls in your buildings, the walls in your model.

00:31:14.039 --> 00:31:34.310
So, even though you could consider it as one boundary condition, a wall, there are a lot of things that go into it and I think we can simplify them into things that are related to heat transfer and things that are related to sheer stress and forces and resistance.

00:31:34.310 --> 00:31:36.555
So let's try heat transfer first.

00:31:36.555 --> 00:31:40.343
How does wall interact with your fluid?

00:31:40.343 --> 00:31:43.596
First of all, the fluid is constrained to your walls.

00:31:43.596 --> 00:31:46.772
It takes the shape of whatever container is put in.

00:31:46.772 --> 00:31:54.859
So, yeah, walls define the shape of your fluid and, of course, the fluid cannot simply penetrate your walls unless they're porous.

00:31:54.859 --> 00:31:58.195
But we'll reach that as it touches the fluid.

00:31:58.195 --> 00:31:59.701
And we're talking fire science.

00:31:59.701 --> 00:32:06.433
So obviously the fluid is not ambient air, it's some outcome of the fire physics happening inside your model.

00:32:06.994 --> 00:32:18.809
There are heat transfer phenomena that are happening and actually there are multiple ways how a wall boundary condition can handle the heat transfer and can interact with the fluid.

00:32:18.809 --> 00:32:22.401
So there are these different types of boundary conditions.

00:32:22.401 --> 00:32:28.759
Formerly they're called the kinds the first, second, third kind, drichlet, neumann, robin, I think.

00:32:28.759 --> 00:32:36.771
But that's a very challenging way to explain it to people who may not be fluid mechanics For me.

00:32:36.771 --> 00:32:39.601
Think about the modes of heat transfer that we have.

00:32:39.601 --> 00:32:46.135
We have radiation, conduction, convection, and that's basically what we try to represent with the model of our boundaries.

00:32:46.738 --> 00:33:00.882
So in the simplest way you can just say that the boundary does not participate in heat transfer, it's adiabatic, it simply does not participate in the heat exchange and whatever heat is released into the fluid, it's not going anywhere.

00:33:00.882 --> 00:33:02.152
So that's the simplest way.

00:33:02.152 --> 00:33:21.441
When you omit the process and in some simulations I know this is practiced I don't see huge benefit in that, but I know in some models, some cases in which adiabatic walls were chosen to make the numerical case more simple, another way is to simply define the temperature of your wall.

00:33:21.441 --> 00:33:23.836
So yeah, that's a way.

00:33:23.836 --> 00:33:28.414
That's a way structural engineers would be defining boundary conditions for their structural models.

00:33:28.414 --> 00:33:31.980
They would tell what temperature the wall has and that's it.

00:33:31.980 --> 00:33:40.820
You can also define the heat flux on the wall, convective heat flux, so how the wall takes the heat from the air.

00:33:40.820 --> 00:33:57.851
Or you can actually make a combination of all of those boundary conditions and make a mixed boundary condition where the wall takes the heat flux from the fluid but its temperature is defined by a heat transfer phenomena inside the wall, and that's the correct way.

00:33:57.851 --> 00:34:00.259
That's the proper way to model them in fire science.

00:34:00.259 --> 00:34:00.903
In my opinion.

00:34:01.325 --> 00:34:03.432
There are some challenges into that.

00:34:03.432 --> 00:34:04.236
Of course there are.

00:34:04.236 --> 00:34:11.695
First is with the convective heat transfer.

00:34:11.695 --> 00:34:15.568
So to assess the convective heat transfer, you need to know the convective heat transfer coefficient, and that's not a simple number to know.

00:34:15.568 --> 00:34:18.855
Eurocode tells you 35 or 8, depending on where the wall is.

00:34:18.855 --> 00:34:23.182
That's watts per kelvin per square meter, if I'm not wrong.

00:34:23.182 --> 00:34:25.728
That's watts per Kelvin per square meter, if I'm not wrong.

00:34:25.728 --> 00:34:35.282
But there are also mathematical models that let you figure out the value of that coefficient, based on Nusselt number or some other simplifications of flow phenomena in boundary layer.

00:34:35.282 --> 00:34:39.239
Anyway, for simple fire modeling we have a value.

00:34:39.239 --> 00:34:43.016
That's a pretty good guess If you want to be a real scientist.

00:34:43.016 --> 00:34:52.016
I know people who are making PhDs of calculating that value and it's pretty fun to actually go deep into what goes into convective heat transfer coefficient.

00:34:52.016 --> 00:34:56.378
We had to go that way when we were modeling linear heat detectors.

00:34:56.378 --> 00:35:05.891
Another paper that's currently being written in my laboratory by my Jakub, who seems to write too many papers at the same time and doesn't have time to finish any of them.

00:35:05.891 --> 00:35:14.215
I will chase him, then you'll have a nice source to learn about convective heat transfer for thin elements.

00:35:14.775 --> 00:35:18.021
Anyway, there's also heat transfer inside into the solid.

00:35:18.021 --> 00:35:20.512
So the heat goes somewhere.

00:35:20.512 --> 00:35:22.195
It doesn't magically disappear.

00:35:22.195 --> 00:35:27.175
So you could imagine that in your fluid dynamics model you also have a solid model.

00:35:27.175 --> 00:35:35.378
So you're actually modeling the solid with discrete cells and you're solving the heat transfer phenomena that go into that solid wall.

00:35:35.378 --> 00:35:40.878
But that would be very computationally heavy and perhaps not necessary for fire science.

00:35:40.878 --> 00:35:47.835
In this case we have simplified models for heat transfer, very simple models based on Fourier's law.

00:35:47.835 --> 00:35:49.882
It's very well known.

00:35:49.882 --> 00:35:53.173
The heat transfer through solids is actually quite well known.

00:35:53.173 --> 00:35:59.956
We have one-dimensional models in which you just simulate the heat penetrating the depth of your material.

00:35:59.956 --> 00:36:09.981
We have a simplified three-dimensional models where the heat actually transfers through a virtual assembly, let's say, and you can track it.

00:36:10.322 --> 00:36:19.164
It's especially useful for structural engineering to have this type of an assessment and with this we are capable to say where the heat went.

00:36:19.164 --> 00:36:20.152
And it's critical.

00:36:20.152 --> 00:36:27.103
Solving heat transfer from your fluid to your structure is critical for two reasons.

00:36:27.103 --> 00:36:30.679
First, for your smoke to cool down.

00:36:30.679 --> 00:36:43.740
Like you have to capture how your smoke is cooling down through heat transfer because this influences the buoyancy of the smoke and buoyancy influences the smoke behavior, the flow velocities, the ceiling jet and everything.

00:36:43.740 --> 00:36:46.775
It's critical to capture that and it's critical.

00:36:46.775 --> 00:36:49.230
The second reason it's critical if you're studying structure.

00:36:49.230 --> 00:36:58.945
So if you want to solve a structural fire engineering case, you need to understand how much heat has went to your structure and then assess what that heat does to your structure.

00:36:59.610 --> 00:37:01.978
Sometimes you would get with extremely complicated models.

00:37:01.978 --> 00:37:06.742
Another proud paper that I participated in with Michał Malendowski.

00:37:06.742 --> 00:37:09.739
We wrote a paper that was Michał's model.

00:37:09.739 --> 00:37:12.057
I cannot take much credit for that.

00:37:12.057 --> 00:37:15.981
I delivered the experimental validation of his crazy ideas.

00:37:15.981 --> 00:37:34.092
The concept was that if you have a very complicated shape, let's say an I-beam, you cannot really model that within your fluid because you cannot afford putting one, two, three centimeter cells to solve the fluid around your very complicated structure of an I-beam.

00:37:34.092 --> 00:37:43.219
However, miho came up with the concept that you could build a virtual box around an I-beam, simply solve the boundaries of the virtual box.

00:37:43.219 --> 00:37:59.201
So now, instead of a very complicated shape of a beam, you have basically a rectangular shape of the virtual box and then with another set of equations you can calculate all the modes of heat transfer within that box to that exact shape of an I-beam.

00:37:59.201 --> 00:38:12.744
And that's like seven pages of very hardcore mathematics that give you a very quick and robust solution to heat transfer problem from fluid into a complicated thin wall structure.

00:38:12.744 --> 00:38:18.302
So yeah, I'll link the paper into show notes if you want to see it for yourself.

00:38:18.302 --> 00:38:19.775
It's quite a sight.

00:38:19.775 --> 00:38:23.460
I'm still waiting for someone to implement that model in a CFD tool.

00:38:23.460 --> 00:38:24.791
But when one does, I think it's going to be implement that model in a CFD tool.

00:38:24.791 --> 00:38:26.724
But when one does, I think it's going to be very, very useful.

00:38:27.425 --> 00:38:28.246
I said thin walls.

00:38:28.246 --> 00:38:34.157
So another aspect of wall boundaries related to heat transfer, how thin or thick they are.

00:38:34.157 --> 00:38:39.737
So if you have a very thick wall, it will behave differently than if you have a very thin wall.

00:38:39.737 --> 00:38:49.619
We sometimes even call them thermally thick or thermally thin, which doesn't really define their width or depth or whatever dimension is called.

00:38:49.619 --> 00:38:59.351
It defines how quickly the heat transferred through them and the point is whether you can assess that the entirety of the wall is at one temperature.

00:38:59.351 --> 00:39:15.286
So if you have a thin sheet of metal think trapezoidal steel sheet that would be a very thin boundary because you can say that the exterior surfaces would have a very similar temperature or the same temperature as your internal surface.

00:39:15.286 --> 00:39:17.117
There's no gradient inside that wall.

00:39:17.117 --> 00:39:25.565
Or you can have a thermally thick wall, let's say a thick concrete wall, where you would have a spectrum, a gradient of temperature going through that wall.

00:39:25.565 --> 00:39:36.976
And also because how long the fire phenomena are that we model, it's very unlikely that in your short CFD simulation, the heat would actually penetrate to the other side of the boundary.

00:39:37.530 --> 00:39:40.018
Now it's important in some cases.

00:39:40.018 --> 00:39:43.079
There were cases, for example, in the Iris project.

00:39:43.079 --> 00:39:53.405
I know they came up with some new observations regarding flashover in compartments with thermally thin walls.

00:39:53.405 --> 00:40:05.590
If I'm not wrong that was a paper authored by Mohamed Bashir and if I'm not wrong he's part of AFAR now, so good career choices there, not wrong.

00:40:05.590 --> 00:40:07.094
He's a part of our farm now, so good career choices there.

00:40:07.094 --> 00:40:12.376
Anyway, they've observed that if you had thermally thin walls, which kind of represents an informal settlement that's what the iris project was studying.

00:40:12.376 --> 00:40:23.233
In that setting, the the flashover cue is in a slightly different manner that you would predict with, let's say, mqh correlations or or some classical compartment fire science.

00:40:23.233 --> 00:40:24.719
Very, very interesting science.

00:40:24.909 --> 00:40:31.579
Perhaps another thing to bring to your attention in the podcast I'm noting down the ideas they come up so quickly.

00:40:31.579 --> 00:40:37.003
So I hope I've summarized the thermal aspects of the solid boundaries.

00:40:37.003 --> 00:40:42.097
But there's one more aspect to the solid boundary and that's their roughness.

00:40:42.097 --> 00:40:52.842
So the solid boundaries not only interact with the flow by exchanging heat with the fluid, but they also interact with the fluid by sheer stress.

00:40:52.842 --> 00:40:59.822
So the fluid flows against the wall and there's some momentum transfer between the wall and the fluid.

00:40:59.822 --> 00:41:03.237
It's because the walls are not perfectly slip.

00:41:03.237 --> 00:41:10.155
They cause resistance to the flow of air and this actually can be quite impactful.

00:41:10.235 --> 00:41:13.643
It's impactful in, I think, two situations for me.

00:41:13.643 --> 00:41:20.152
First situation is where I'm modeling ductwork, so it's critical to capture the pressure losses in my ductwork.

00:41:20.152 --> 00:41:32.271
I need to get this shear stress correct Because if I don't, my flows in the duct work will be completely off, my pressure point of my fan will be completely off, my CFD will be bullshit.

00:41:32.271 --> 00:41:34.516
So I don't want my CFD to be bullshit.

00:41:34.516 --> 00:41:37.170
I want my CFD to capture reality as close as possible.

00:41:37.170 --> 00:41:42.365
So that's one case where we really need to control that, that slip condition.

00:41:43.268 --> 00:41:45.342
And the other one is when we're doing wind engineering.

00:41:45.342 --> 00:41:55.500
So in wind engineering we would introduce a very, very large numerical domains, so the wind profile is not affected by the building that is inside of that domain.

00:41:55.500 --> 00:42:00.175
It's something called the blockade effect that we want to avoid.

00:42:00.175 --> 00:42:07.195
So our domain must be significantly larger than the building that's subject to the wind or neighborhood.

00:42:07.195 --> 00:42:20.659
If we're modeling a neighborhood, then it's kilometers in size, and when I have kilometers of space between my inlet point on which I meticulously put my wind profile with a profile function logarithmic wind profile.

00:42:20.659 --> 00:42:36.940
Usually I want that wind profile to reach my building, which I'm modeling, but when the flow is flowing through my domain it interacts with the ground because the ground has roughness in it and this roughness actually changes the wind profile.

00:42:36.940 --> 00:42:45.456
So if the wind flows for one kilometer until it reaches the building of my interest, I can have a completely different wind profile than the one that I've defined.

00:42:46.179 --> 00:42:47.403
There are ways to handle that.

00:42:47.403 --> 00:42:54.585
We have in-house modifications of ANSYS code which my postdoc, paulina, has published.

00:42:54.585 --> 00:43:06.530
We've took some clever models from the literature and applied them into our ANSYS fluent simulations and this allows us to precisely control what's happening between the air and the ground.

00:43:06.530 --> 00:43:12.061
How does the roughness of the ground, the surface, the affects the flow?

00:43:12.061 --> 00:43:23.251
And actually it's quite important because in wind engineering that's a well-known technique to uh to introduce the turbulence intensity and some shape of wind profile that you would expect in a very specific terrain case.

00:43:23.652 --> 00:43:33.262
Actually, if you go deeper in that, what's happening between the house, between the flow and the wall, is something we often refer to as a boundary layer problem.

00:43:33.262 --> 00:43:42.690
So a boundary layer problem is exactly that interface between the fluid and the solid and it's something that we're very ignorant in fire science about.

00:43:42.690 --> 00:44:02.181
So you don't really see good boundary layer control in cfd simulations for fires, perhaps because the velocities are very small, perhaps because it usually is not a part of interest, as I said, it's for me it's wind engineering or ducts, and this would not be the main parts of cfd analysis in fire now if I think about it.

00:44:02.181 --> 00:44:05.820
But yeah, boundary layer problems are are a hell of a problem.

00:44:05.820 --> 00:44:08.338
They are very, very challenging to solve.

00:44:08.338 --> 00:44:14.596
They require fine meshes, they require near-wall models, they require modifications to your turbulence modeling.

00:44:14.596 --> 00:44:22.815
They are quite a fun of a problem to solve with and if you ever have problems with them you can reach me out.

00:44:22.815 --> 00:44:27.896
Perhaps the way we've solved it with paulina for our wind studies will will help you out.

00:44:27.896 --> 00:44:32.128
Um, that would be it quickly for the wall boundaries.

00:44:32.309 --> 00:44:35.597
Okay, there's one more wall boundary I've kind of teased.

00:44:35.597 --> 00:44:50.353
That's a porous boundary condition, so boundary condition in which some of the flow can flow through your porous, medium porous wall, and we sometimes use this boundary condition to simulate inlet grills.

00:44:50.353 --> 00:44:57.617
So you can actually control the flows on your complex duct network by those boundary conditions, and they would be very efficient in doing that.

00:44:57.617 --> 00:45:14.217
However, it's quite difficult to get the data necessary to define those, and in my career I've also used the porous boundary condition to actually model some baffles on the roofs of buildings to break wind acting on the roof of the building.

00:45:14.217 --> 00:45:16.003
So that was an interesting case study.

00:45:16.525 --> 00:45:25.016
I know my friends also use porous boundary conditions to simulate resistance of a tunnel network, like if you're modeling a road tunnel or something.

00:45:25.016 --> 00:45:26.536
It's very similar to the duct.

00:45:26.536 --> 00:45:32.523
In that case the pressure loss along the length of the tunnel is important and you want to capture that.

00:45:32.523 --> 00:45:43.177
You can do it by simulating the entirety of the tunnel and the roughness, and you can also build some surrogate model in form of porous boundary condition and it actually kind of works.

00:45:43.177 --> 00:45:52.137
It solves the problem, allows you to input an artificial resistance to your tunnel that mimics the pressure losses that you would have in your tunneling system.

00:45:52.137 --> 00:45:54.429
So we've reached the end of the solid boundaries.

00:45:54.590 --> 00:45:59.577
I've promised you special boundary conditions, so we have a magical boundary condition called the symmetry.

00:45:59.577 --> 00:46:03.219
It basically is what it says it is.

00:46:03.219 --> 00:46:04.373
It's a symmetry plane.

00:46:04.373 --> 00:46:11.125
So if you're modeling a duct, you can model half of the duct and assume that the other half is symmetrical.

00:46:11.125 --> 00:46:13.632
It's interesting boundary conditions.

00:46:13.632 --> 00:46:15.657
A lot of people are using it in a funny way.

00:46:15.657 --> 00:46:25.514
So some people are using it as a top of their wind tunnel in their modeling, which is kind of ridiculous because that means that you have a mirror image of your wind tunnel in your simulation.

00:46:25.514 --> 00:46:33.041
That's not really great, capturing physics of your model, but it's a boundary condition that came to life because of savings.

00:46:33.041 --> 00:46:36.614
Necessary Time is money and CFD simulations cost a lot of money.

00:46:36.614 --> 00:46:53.456
When you have a symmetrical case, especially you're simulating something simple like a duct, then symmetric case is your best friend there in answers we have a very interesting boundary conditions called interiors, though that's basically the connections between every single cell in your model.

00:46:53.856 --> 00:47:00.500
But I can also artificially put that condition on any other uh, let's say physical obstacle in my model which makes them disappear.

00:47:00.500 --> 00:47:08.297
In fds you have the magical boundary condition called the hole, where you put a hole in the wall and that allows the flow to go through.

00:47:08.297 --> 00:47:15.199
In ANSYS we call that interior, so it changes the boundary condition for something that's completely permeable by the air.

00:47:15.199 --> 00:47:19.101
And there's also a variant of that boundary condition called the interface.

00:47:19.101 --> 00:47:27.041
So I can actually have two different meshes in my model and connect them with interface, so they don't match each other perfectly.

00:47:27.041 --> 00:47:39.483
But the interface boundary condition would work out how the flow from one domain flows into the other domain, very, very convenient, especially in wind engineering when we have to rotate our models.

00:47:39.483 --> 00:47:47.143
So I have my domain changing every time I change the wind angle because I change it through the rotation of the internal domain.

00:47:47.871 --> 00:47:54.681
Anyway, I was rushing at the end because I see the timer and I'm already talking about boundary conditions for like 15 minutes.

00:47:54.681 --> 00:47:55.813
That's crazy.

00:47:55.813 --> 00:47:58.195
But yeah, that's it, we've reached the end of my list.

00:47:58.195 --> 00:48:00.277
So we've went through the flow conditions.

00:48:00.277 --> 00:48:03.860
We've talked about the simple pressure inlets or open conditions.

00:48:03.860 --> 00:48:05.335
We've talked about velocity conditions.

00:48:05.335 --> 00:48:06.740
We've talked about the simple pressure inlets or open conditions.

00:48:06.740 --> 00:48:09.931
We've talked about velocity conditions, the velocity inlets, and we've touched a bit the mass flow inlets.

00:48:09.931 --> 00:48:13.340
We've talked about fans, fan boundary conditions.

00:48:13.340 --> 00:48:17.378
We talked about HVC capabilities of FDS, which are magnificent.

00:48:17.378 --> 00:48:23.202
Then we moved into the walls, the solid boundaries of your model, and we've discussed modeling the heat transfer phenomena model.

00:48:23.202 --> 00:48:28.277
And we've discussed modeling the heat transfer phenomena and we've discussed modeling shear stress on those boundaries.

00:48:28.277 --> 00:48:34.260
So I hope I took you on an interesting adventure through boundary conditions in numerical modeling.

00:48:34.320 --> 00:48:42.175
I've tried to give you as much real life examples as I could that are relevant to use the good, correct use of those boundary conditions.

00:48:42.175 --> 00:48:59.139
And what I hope to achieve with this episode is that the next time you're building your model, you perhaps will reflect on what boundary conditions are you using, what physics they represent, what role do they play in your model and if everything is perfectly fine with them.

00:48:59.139 --> 00:49:03.817
I think such a reflection is necessary for every CFD model.

00:49:03.817 --> 00:49:09.443
If there's an obvious error and you see it and you fix it, that's not a big problem for CFD.

00:49:09.443 --> 00:49:21.179
If there's an error that you don't see and don't recognize, it is there and you get the faulty results which you are not aware that they are faulty, that's an issue and to avoid that we need to know our boundaries.

00:49:21.789 --> 00:49:24.278
Thank you very much for being here with me.

00:49:24.278 --> 00:49:28.161
If you stayed till the end of this very difficult episode.

00:49:28.161 --> 00:49:41.097
I thank you twice Great job, thank you for staying here with me and I know you're a big fan of the Fire Science Show and I am sure I will see you here next Wednesday with the next episode of the Fire Science Show.

00:49:41.097 --> 00:49:46.675
And once again, congratulations, rory, for being the first episode ever in the fire science show to break.

00:49:46.675 --> 00:49:50.224
2 000 listens, insane 2 000 people.

00:49:50.224 --> 00:49:51.447
That's, that's crazy.

00:49:51.447 --> 00:49:52.590
Thank you very much.

00:49:52.590 --> 00:49:53.530
See you bye.

00:49:53.530 --> 00:50:20.840
Thank you.