|
|
Numerical study of the steady airflow in the human respiratory system during inhaling and exhaling
Lancmanová, Anna ; Bodnár, Tomáš
This paper presents some of the initial results of the numerical simulations of a steady turbulent flow in human upper airways during inhalation and exhalation. The mathematical model is based on the system of Reynolds-Averaged incompressible Navier-Stokes equations complemented by the SST k − ω turbulence model. The simulations were performed using finite-volume open source solver OpenFOAM on a realistic three-dimensional geometry. The main aim of this particular study is to verify the computational setup with special focus on appropriate choice and implementation of boundary conditions. The prescribed boundary conditions are chosen to mimic the physiological conditions during normal breathing cycle. This study aims to gain an insight into the airflow behavior during the inhalation and exhalation process by comparing the results of two distinct simulations corresponding to two different (opposite) flow rates . The obtained local flow rates and flow fields for both cases are presented and mutually compared. This initial work should serve as a foundation for future more complex simulations that will include the time-dependent and compressible effects.
|
|
|
Numerical evaluation of mass-diffusive compressible fluids flows models
Bodnár, Tomáš ; Fraunié, P.
This contribution presents first numerical tests of some recently published alternative models for solution of viscous compressible and nearly incompressible models. All models are solved by high resolution compact finite difference scheme with strong stability preserving RungeKutta time stepping. The two simple but challenging computational test cases are presented, based on the double-periodic shear layer and the Kelvin-Helmholtz instability. The obtained time-dependent flow fields are showing pronounced shear and vorticity layers being resolved by the standard as well as by the new mass-diffusive modified models. The preliminary results show that the new models are viable alternative to the well established classical models.
|
|
|
Interpolation with restrictions -- role of the boundary conditions and individual restrictions
Valášek, Jan ; Sváček, P.
The contribution deals with the remeshing procedure between two computational finite element meshes. The remeshing represented by the interpolation of an approximate solution onto a new mesh is needed in many applications like e.g. in aeroacoustics, here we are particularly interested in the numerical flow simulation of a gradual channel collapse connected with a~severe deterioration of the computational mesh quality. Since the classical Lagrangian projection from one mesh to another is a dissipative method not respecting conservation laws, a conservative interpolation method introducing constraints is described. The constraints have form of Lagrange multipliers enforcing conservation of desired flow quantities, like e.g. total fluid mass, flow kinetic energy or flow potential energy. Then the interpolation problem turns into an error minimization problem, such that the resulting quantities of proposed interpolation satisfy these physical properties while staying as close as possible to the results of Lagrangian interpolation in the L2 norm. The proposed interpolation scheme does not impose any restrictions on mesh generation process and it has a relatively low computational cost. The implementation details are discussed and test cases are shown.
|
|
|
On the development of a numerical model for the simulation of air flow in the human airways
Lancmanová, Anna ; Bodnár, Tomáš ; Sequeira, A.
This contribution reports on an ongoing study focusing on reduced order models for incompressible viscous fluid flow in two dimensional channels. A finite difference solver was developed using a simple implementation of the immersed boundary method to represent the channel geometry. The solver was validated for unsteady flow by comparing the obtained two-dimensional numerical solutions with analytical profiles computed from the Womersley solution. Finally the 2D model was coupled to a simple 1D extension simulating the flow in axisymmetric elastic vessel (tube). Some of the coupling principles and implementation issues are discussed in detail.
|
|
| |
guest ::
