fluid–structure interaction

Fluid–structure interaction (FSI) is the interaction of some movable or deformable structure with an internal or surrounding fluid flow.{{cite book | editor1-last=Bungartz | editor1-first=Hans-Joachim | editor2-last=Schäfer | editor2-first=Michael | year=2006 | title=Fluid-structure Interaction: Modelling, Simulation, Optimization | publisher=Springer-Verlag | isbn=978-3-540-34595-4}} Fluid–structure interactions can be stable or oscillatory. In oscillatory interactions, the strain induced in the solid structure causes it to move such that the source of strain is reduced, and the structure returns to its former state only for the process to repeat. Image:Pressure Wave In Tube.png

Examples

Fluid–structure interactions are a crucial consideration in the design of many engineering systems, e.g. automobile, aircraft, spacecraft, engines and bridges. Failing to consider the effects of oscillatory interactions can be catastrophic, especially in structures comprising materials susceptible to fatigue. Tacoma Narrows Bridge (1940), the first Tacoma Narrows Bridge, is probably one of the most infamous examples of large-scale failure. Aircraft wings and turbine blades can break due to FSI oscillations. A reed actually produces sound because the system of equations governing its dynamics has oscillatory solutions. The dynamic of reed valves used in two strokes engines and compressors is governed by FSI. The act of "blowing a raspberry" is another such example. The interaction between tribological machine components, such as bearings and gears, and lubricant is also an example of FSI.{{Cite journal|last1=Singh|first1=Kushagra|last2=Sadeghi|first2=Farshid|last3=Russell|first3=Thomas|last4=Lorenz|first4=Steven J.|last5=Peterson|first5=Wyatt|last6=Villarreal|first6=Jaret|last7=Jinmon|first7=Takumi|date=2021-09-01|title=Fluid–Structure Interaction Modeling of Elastohydrodynamically Lubricated Line Contacts|url=https://asmedigitalcollection.asme.org/tribology/article/143/9/091602/1091936/Fluid-Structure-Interaction-Modeling-of|journal=Journal of Tribology|language=en|volume=143|issue=9|doi=10.1115/1.4049260|s2cid=230619508 |issn=0742-4787|url-access=subscription}} The lubricant flows between the contacting solid components and causes elastic deformation in them during this process. Fluid–structure interactions also occur in moving containers, where liquid oscillations due to the container motion impose substantial magnitudes of forces and moments to the container structure that affect the stability of the container transport system in a highly adverse manner.{{Cite journal|last1=Kolaei|first1=Amir|last2=Rakheja|first2=Subhash|last3=Richard|first3=Marc J.|date=2016-01-25|title=An efficient methodology for simulating roll dynamics of a tank vehicle coupled with transient fluid slosh|journal=Journal of Vibration and Control|volume=23|issue=19|pages=3216–3232|doi=10.1177/1077546315627565|s2cid=123621791 |issn=1077-5463}}{{Cite journal|last1=Kolaei|first1=Amir|last2=Rakheja|first2=Subhash|last3=Richard|first3=Marc J.|date=2015-09-01|title=Three-dimensional dynamic liquid slosh in partially-filled horizontal tanks subject to simultaneous longitudinal and lateral excitations|journal=European Journal of Mechanics B|volume=53|pages=251–263|doi=10.1016/j.euromechflu.2015.06.001|bibcode=2015EJMF...53..251K}}{{Cite journal|last1=Kolaei|first1=Amir|last2=Rakheja|first2=Subhash|last3=Richard|first3=Marc J.|date=2014-01-06|title=Range of applicability of the linear fluid slosh theory for predicting transient lateral slosh and roll stability of tank vehicles|journal=Journal of Sound and Vibration|volume=333|issue=1|pages=263–282|doi=10.1016/j.jsv.2013.09.002|bibcode=2014JSV...333..263K}}{{Cite journal|last1=Kolaei|first1=Amir|last2=Rakheja|first2=Subhash|last3=Richard|first3=Marc J.|date=2014-07-01|title=Effects of tank cross-section on dynamic fluid slosh loads and roll stability of a partly-filled tank truck|journal=European Journal of Mechanics B|volume=46|pages=46–58|doi=10.1016/j.euromechflu.2014.01.008|bibcode=2014EJMF...46...46K}} Another prominent example is the start up of a rocket engine, e.g. Space Shuttle main engine (SSME), where FSI can lead to considerable unsteady side loads on the nozzle structure.{{Cite thesis|last=Manuel|first=Frey|date=2001|title=Behandlung von Strömungsproblemen in Raketendüsen bei Überexpansion|url=http://elib.uni-stuttgart.de/handle/11682/3667|language=de|doi=10.18419/opus-3650}} In addition to pressure-driven effects, FSI can also have a large influence on surface temperatures on supersonic and hypersonic vehicles.{{Cite journal|last1=Daub|first1=Dennis|last2=Esser|first2=Burkard|last3=Gülhan|first3=Ali|date=April 2020|title=Experiments on High-Temperature Hypersonic Fluid–Structure Interaction with Plastic Deformation|journal=AIAA Journal|language=en|volume=58|issue=4|pages=1423–1431|doi=10.2514/1.J059150|bibcode=2020AIAAJ..58.1423D |issn=0001-1452|doi-access=free}}

Fluid–structure interactions also play a major role in appropriate modeling of blood flow. Blood vessels act as compliant tubes that change size dynamically when there are changes to blood pressure and velocity of flow.{{Cite journal|last=Taren|first=James A.|date=1965|title=Cerebral Aneurysm|journal=The American Journal of Nursing|volume=65|issue=4|pages=88–91|doi=10.2307/3453223|issn=0002-936X|jstor=3453223|pmid=14258014 |s2cid=31190911 }} Failure to take into account this property of blood vessels can lead to a significant overestimation of resulting wall shear stress (WSS). This effect is especially imperative to take into account when analyzing aneurysms. It has become common practice to use computational fluid dynamics to analyze patient specific models. The neck of an aneurysm is the most susceptible to changes in to WSS. If the aneurysmal wall becomes weak enough, it becomes at risk of rupturing when WSS becomes too high. FSI models contain an overall lower WSS compared to non-compliant models. This is significant because incorrect modeling of aneurysms could lead to doctors deciding to perform invasive surgery on patients who were not at a high risk of rupture. While FSI offers better analysis, it comes at a cost of highly increased computational time. Non-compliant models have a computational time of a few hours, while FSI models could take up to 7 days to finish running. This leads to FSI models to be most useful for preventative measures for aneurysms caught early, but unusable for emergency situations where the aneurysm may have already ruptured.{{Cite journal|last1=Sforza|first1=Daniel M.|last2=Putman|first2=Christopher M.|last3=Cebral|first3=Juan R.|date=June 2012|title=Computational fluid dynamics in brain aneurysms|journal=International Journal for Numerical Methods in Biomedical Engineering|volume=28|issue=6–7|pages=801–808|doi=10.1002/cnm.1481|issn=2040-7939|pmc=4221804|pmid=25364852}}{{Cite journal|last1=Khe|first1=A K|last2=Cherevko|first2=A A|last3=Chupakhin|first3=A P|last4=Bobkova|first4=M S|last5=Krivoshapkin|first5=A L|last6=Orlov|first6=K Yu|date=June 2016|title=Haemodynamics of giant cerebral aneurysm: A comparison between the rigid-wall, one-way and two-way FSI models|journal=Journal of Physics: Conference Series|volume=722|issue=1|pages=012042|doi=10.1088/1742-6596/722/1/012042|issn=1742-6588|bibcode=2016JPhCS.722a2042K|doi-access=free}}{{Cite journal|last1=Torii|first1=Ryo|last2=Oshima|first2=Marie|last3=Kobayashi|first3=Toshio|last4=Takagi|first4=Kiyoshi|last5=Tezduyar|first5=Tayfun E.|date=2009-09-15|title=Fluid–structure interaction modeling of blood flow and cerebral aneurysm: Significance of artery and aneurysm shapes|journal=Computer Methods in Applied Mechanics and Engineering|series=Models and Methods in Computational Vascular and Cardiovascular Mechanics|volume=198|issue=45|pages=3613–3621|doi=10.1016/j.cma.2008.08.020|issn=0045-7825|bibcode=2009CMAME.198.3613T}}{{Cite journal|last1=Razaghi|first1=Reza|last2=Biglari|first2=Hasan|last3=Karimi|first3=Alireza|date=2019-07-01|title=Risk of rupture of the cerebral aneurysm in relation to traumatic brain injury using a patient-specific fluid-structure interaction model|journal=Computer Methods and Programs in Biomedicine|volume=176|pages=9–16|doi=10.1016/j.cmpb.2019.04.015|pmid=31200915 |s2cid=155305862 |issn=0169-2607}}

Analysis

Fluid–structure interaction problems and multiphysics problems in general are often too complex to solve analytically and so they have to be analyzed by means of experiments or numerical simulation. Research in the fields of computational fluid dynamics and computational structural dynamics is still ongoing but the maturity of these fields enables numerical simulation of fluid-structure interaction.J.F. Sigrist (2015). Fluid-Structure Interaction: An Introduction to Finite Element Coupling. Wiley ({{ISBN|978-1-119-95227-5}}) Two main approaches exist for the simulation of fluid–structure interaction problems:

  • Monolithic approach: the equations governing the flow and the displacement of the structure are solved simultaneously, with a single solver
  • Partitioned approach: the equations governing the flow and the displacement of the structure are solved separately, with two distinct solvers

The monolithic approach requires a code developed for this particular combination of physical problems whereas the partitioned approach preserves software modularity because an existing flow solver and structural solver are coupled. Moreover, the partitioned approach facilitates solution of the flow equations and the structural equations with different, possibly more efficient techniques which have been developed specifically for either flow equations or structural equations. On the other hand, development of stable and accurate coupling algorithm is required in partitioned simulations. In conclusion, the partitioned approach allows reusing existing software which is an attractive advantage. However, stability of the coupling method needs to be taken into consideration. This is especially difficult, if the mass of the moving structure is small in comparison to the mass of fluid which is displaced by the structure movement.

In addition, the treatment of meshes introduces other classifications of FSI analysis. For example,one can classify them as the conforming mesh methods and the non-conforming mesh methods.{{Cite web |url=http://www.global-sci.com/openaccess/v12_337.pdf |title=Archived copy |access-date=2014-11-28 |archive-date=2014-10-31 |archive-url=https://web.archive.org/web/20141031190349/http://www.global-sci.com/openaccess/v12_337.pdf |url-status=dead }} Other classifications can be mesh-based methods and meshless methods.{{Cite journal|last1=Zhang|first1=Chi|last2=Rezavand|first2=Massoud|last3=Hu|first3=Xiangyu|date=2021-03-15|title=A multi-resolution SPH method for fluid-structure interactions|journal=Journal of Computational Physics|volume=429|pages=110028|doi=10.1016/j.jcp.2020.110028|arxiv=1911.13255 |bibcode=2021JCoPh.42910028Z |s2cid=208513116 }}

Numerical simulation

The Newton–Raphson method or a different fixed-point iteration can be used to solve FSI problems. Methods based on Newton–Raphson iteration are used in both the monolithic

{{cite journal |doi=10.1016/j.cma.2003.09.006 |author=M. Heil |title=An efficient solver for the fully coupled solution of large-displacement fluid-structure interaction problems |journal=Computer Methods in Applied Mechanics and Engineering |volume=193 |issue=1–2 |pages=1–23 |year=2004|bibcode = 2004CMAME.193....1H }}{{cite journal |doi=10.1002/nme.959 |author1=K.-J. Bathe |author-link=Klaus-Jürgen Bathe | author2=H. Zhang |title=Finite element developments for general fluid flows with structural interactions |journal=International Journal for Numerical Methods in Engineering |volume=60 |issue=1 |pages=213–232 |year=2004|bibcode = 2004IJNME..60..213B |citeseerx=10.1.1.163.1531 |s2cid=17143434 }}{{cite book |author=J. Hron, S. Turek |title=A monolithic FEM/multigrid solver for ALE formulation of fluid-structure interaction with application in biomechanics |editor=H.-J. Bungartz | editor2=M. Schäfer | volume=Fluid–Structure Interaction – Modelling, Simulation, Optimisation | series=Lecture Notes in Computational Science and Engineering | publisher=Springer-Verlag | year=2006 | pages=146–170 | isbn=978-3-540-34595-4}} and the partitioned {{cite journal |doi=10.1016/S0045-7949(02)00409-1 |author1=H. Matthies |author2=J. Steindorf |title=Partitioned strong coupling algorithms for fluid-structure interaction |journal=Computers and Structures |volume=81 |issue=8–11 |pages=805–812 |year=2003|citeseerx=10.1.1.487.5577 }}{{cite journal |doi=10.1016/j.cma.2004.11.032 |author1=H. Matthies |author2=R. Niekamp |author3=J. Steindorf |title=Algorithms for strong coupling procedures |journal=Computer Methods in Applied Mechanics and Engineering |volume=195 |issue=17–18 |pages=2028–2049 |year=2006|bibcode = 2006CMAME.195.2028M }} approach. These methods solve the nonlinear flow equations and the structural equations in the entire fluid and solid domain with the Newton–Raphson method. The system of linear equations within the Newton–Raphson iteration can be solved without knowledge of the Jacobian with a matrix-free iterative method, using a finite difference approximation of the Jacobian-vector product.

Whereas Newton–Raphson methods solve the flow and structural problem for the state in the entire fluid and solid domain, it is also possible to reformulate an FSI problem as a system with only the degrees of freedom in the interface’s position as unknowns. This domain decomposition condenses the error of the FSI problem into a subspace related to the interface.{{cite journal | doi=10.1016/j.cma.2005.01.018 |author1=C. Michler |author2=E. van Brummelen |author3=R. de Borst | title=Error-amplification analysis of subiteration-preconditioned GMRES for fluid-structure interaction | journal=Computer Methods in Applied Mechanics and Engineering | volume=195 |issue=17–18 | year=2006 | pages=2124–2148|bibcode = 2006CMAME.195.2124M }} The FSI problem can hence be written as either a root finding problem or a fixed point problem, with the interface’s position as unknowns.

Interface Newton–Raphson methods solve this root-finding problem with Newton–Raphson iterations, e.g. with an approximation of the Jacobian from a linear reduced-physics model.{{cite journal | doi=10.1051/m2an:2003049 |author1=J.-F. Gerbeau |author2=M. Vidrascu | title=A quasi-Newton algorithm based on a reduced model for fluid-structure interaction problems in blood flows | journal=ESAIM: Mathematical Modelling and Numerical Analysis | volume=37 | issue=4 | year=2003 | pages=631–648| url=https://hal.inria.fr/inria-00071895/file/RR-4691.pdf }}{{cite journal | doi=10.1016/j.compstruc.2004.03.083 |author1=J.-F. Gerbeau |author2=M. Vidrascu |author3=P. Frey | title=Fluid-structure interaction in blood flows on geometries based on medical imaging | journal=Computers and Structures | volume=83 | issue=2–3 | year=2005 | pages=155–165}} The interface quasi-Newton method with approximation for the inverse of the Jacobian from a least-squares model couples a black-box flow solver and structural solver {{cite journal |author1=J. Degroote |author2=K.-J. Bathe |author3=J. Vierendeels | title=Performance of a new partitioned procedure versus a monolithic procedure in fluid–structure interaction | journal=Computers and Structures | year=2009 | doi=10.1016/j.compstruc.2008.11.013 | volume=87 |issue=11–12 | pages=793–801|citeseerx=10.1.1.163.827 }} by means of the information that has been gathered during the coupling iterations. This technique is based on the interface block quasi-Newton technique with an approximation for the Jacobians from least-squares models which reformulates the FSI problem as a system of equations with both the interface’s position and the stress distribution on the interface as unknowns. This system is solved with block quasi-Newton iterations of the Gauss–Seidel type and the Jacobians of the flow solver and structural solver are approximated by means of least-squares models.{{cite journal | doi=10.1016/j.compstruc.2006.11.006 |author1=J. Vierendeels |author2=L. Lanoye |author3=J. Degroote |author4=P. Verdonck | title=Implicit coupling of partitioned fluid-structure interaction problems with reduced order models | journal=Computers and Structures | volume=85 | issue=11–14 | year=2007 | pages=970–976}}

The fixed-point problem can be solved with fixed-point iterations, also called (block) Gauss–Seidel iterations, which means that the flow problem and structural problem are solved successively until the change is smaller than the convergence criterion. However, the iterations converge slowly if at all, especially when the interaction between the fluid and the structure is strong due to a high fluid/structure density ratio or the incompressibility of the fluid.{{cite journal | doi=10.1016/j.cma.2004.12.005 |author1=P. Causin |author2=J.-F. Gerbeau |author3=F. Nobile | title=Added-mass effect in the design of partitioned algorithms for fluid-structure problems | journal=Computer Methods in Applied Mechanics and Engineering | volume=194 | issue=42–44 | year=2005 | pages=4506–4527|bibcode = 2005CMAME.194.4506C |s2cid=122528121 |url=https://hal.inria.fr/inria-00071499/file/RR-5084.pdf }} The convergence of the fixed point iterations can be stabilized and accelerated by Aitken relaxation and steepest descent relaxation, which adapt the relaxation factor in each iteration based on the previous iterations.{{cite journal |author1=U. Küttler |author2=W. Wall | title=Fixed-point fluid-structure interaction solvers with dynamic relaxation | journal=Computational Mechanics | volume=43 | issue=1 | year=2008 | pages=61–72 | doi=10.1007/s00466-008-0255-5|bibcode = 2008CompM..43...61K |s2cid=122209351 }}

If the interaction between the fluid and the structure is weak, only one fixed-point iteration is required within each time step. These so-called staggered or loosely coupled methods do not enforce the equilibrium on the fluid–structure interface within a time step but they are suitable for the simulation of aeroelasticity with a heavy and rather stiff structure.

Several studies have analyzed the stability of partitioned algorithms for the simulation of fluid-structure interaction{{cite journal | doi=10.1016/j.compstruc.2008.05.005 |author1=J. Degroote |author2=P. Bruggeman |author3=R. Haelterman |author4=J. Vierendeels | title=Stability of a coupling technique for partitioned solvers in FSI applications | journal=Computers and Structures | volume=86 | issue=23–24 | year=2008 | pages=2224–2234|hdl=1854/LU-533350 |url=https://biblio.ugent.be/publication/533350 | hdl-access=free }}

{{cite journal | doi=10.1016/j.jcp.2006.02.016 |author1=R. Jaiman |author2=X. Jiao |author3=P. Geubelle |author4=E. Loth | title=Conservative load transfer along curved fluid-solid interface with non-matching meshes | journal=Journal of Computational Physics | volume=218 | issue=1 | year=2006 | pages=372–397|bibcode = 2006JCoPh.218..372J |citeseerx=10.1.1.147.4391 }}

.{{cite journal | doi=10.2514/1.3660 |author1=J. Vierendeels |author2=K. Dumont |author3=E. Dick |author4=P. Verdonck | title=Analysis and stabilization of fluid-structure interaction algorithm for rigid-body motion | journal=AIAA Journal | volume=43 | issue=12 | year=2005 | pages=2549–2557|bibcode = 2005AIAAJ..43.2549V |url=https://zenodo.org/record/894358}}{{cite conference | conference=European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 |editor1=P. Wesseling |editor2=E. Oñate |editor3=J. Périaux | location=The Netherlands | year=2006 | title=The artificial added mass effect in sequential staggered fluid-structure interaction algorithms |author1=Christiane Förster |author2=Wolfgang A. Wall |author3=Ekkehard Ramm }}{{cite journal | doi=10.1016/j.cma.2006.09.002 | journal=Computer Methods in Applied Mechanics and Engineering | volume=196 | issue=7 | year=2007 | pages=1278–1293 | title=Artificial added mass instabilities in sequential staggered coupling of nonlinear structures and incompressible viscous flows |author1=Christiane Förster |author2=Wolfgang A. Wall |author3=Ekkehard Ramm |bibcode = 2007CMAME.196.1278F }}

See also

Open source codes

  • [https://bitbucket.org/philip_cardiff/solids4foam-release/src/master/solids4foam solids4Foam], a toolbox for OpenFOAM with capabilities for solid mechanics and fluid solid interactions
  • [http://www.oomph-lib.org/ oomph-lib]
  • [https://web.archive.org/web/20081203094902/http://www.csc.fi/english/pages/elmer/examples/Coupled_problems/hemodyn Elmer FSI page]
  • [https://launchpad.net/cbc.solve/ CBC.solve Biomedical Solvers]
  • [http://www.precice.org preCICE Coupling Library]
  • [https://www.sphinxsys.org/ SPHinXsys multi-physics library] It provides C++ APIs for physical accurate simulation and aims to model coupled industrial dynamic systems including fluid, solid, multi-body dynamics and beyond with SPH (smoothed particle hydrodynamics), a meshless computational method using particle discretization.

Academic Codes

  • [http://www.math.ucsb.edu/~atzberg/SIB_Codes/index.html Stochastic Immersed Boundary Methods in 3D, P. Atzberger, UCSB]
  • [http://www.math.nyu.edu/~griffith/IBAMR/ Immersed Boundary Method for Adaptive Meshes in 3D, B. Griffith, NYU.]
  • [http://www.math.utah.edu/IBIS/ Immersed Boundary Method for Uniform Meshes in 2D, A. Fogelson, Utah]
  • [https://www.tu-braunschweig.de/ifl/simulationswerkzeuge/ifls/ IFLS, IFL, TU Braunschweig]

Commercial Codes

  • [http://www.simulia.com/products/multiphysics_coupling.html Abaqus Multiphysics Coupling]
  • [http://www.altairhyperworks.com/Product,54,AcuSolve.aspx?AspxAutoDetectCookieSupport=1 AcuSolve FSI applications]
  • [http://www.adina.com/fluid-structure-interaction.shtml ADINA FSI homepage] {{Webarchive|url=https://web.archive.org/web/20210428202131/http://www.adina.com/fluid-structure-interaction.shtml |date=2021-04-28 }}
  • [https://www.ansys.com/products/platform/multiphysics-simulation/fluid-structure-interaction Ansys' FSI homepage]
  • [https://web.archive.org/web/20110128234413/http://epsilon-sa.com/ Altair RADIOSS]
  • [https://web.archive.org/web/20131223074755/http://www.autodesk.com/products/autodesk-simulation-family/features/simulation-cfd/all/gallery-view Autodesk Simulation CFD]
  • Simcenter STAR-CCM+ from Siemens Digital Industries Software
  • [https://web.archive.org/web/20110620191613/http://www.even-ag.ch/index.php/software/colyx CoLyX - FSI and mesh-morphing from EVEN - Evolutionary Engineering AG]
  • [http://www.fluidyn.com/fluidyn/multi-physics/coupling-fluids-structures Fluidyn-MP FSI Multiphysics Coupling] {{Webarchive|url=https://web.archive.org/web/20140517122229/http://www.fluidyn.com/fluidyn/multi-physics/coupling-fluids-structures |date=2014-05-17 }}
  • [http://www.comsol.com/showroom/gallery/361/ COMSOL FSI homepage]
  • [http://www.mpcci.de/ MpCCI homepage]
  • [http://www.mscsoftware.com/Products/CAE-Tools/MD-Nastran.aspx MSC Software MD Nastran]
  • [http://www.mscsoftware.com/Products/CAE-Tools/Dytran.aspx MSC Software Dytran]
  • [http://www.open-engineering.com/index.php/eng/Applications/FSI FINE/Oofelie FSI: Fully integrated and strongly coupled for better convergence] {{Webarchive|url=https://web.archive.org/web/20120401020239/http://www.open-engineering.com/index.php/eng/Applications/FSI |date=2012-04-01 }}
  • [http://www.lstc.com/ LS-DYNA Home Page]
  • [https://archive.today/20130123235636/http://www.fluidyn.com/fluidyn/coupling-fluid-and-structures Fluidyn-MP FSI: Fluid-Structure Interaction]
  • [http://www.compassis.com/compass/en/Productos/Tdyn CompassFEM Tdyn]
  • [http://www.compassis.com/compass/en/Productos/SeaFEM CompassFEM SeaFEM]
  • [http://www.cradle-cfd.com/products/sctetra/index.html Cradle SC/Tetra CFD Software]
  • [http://www.cimne.com/parachutes PARACHUTES FSI HomePage]

References

{{reflist}}

Further reading

Modarres-Sadeghi, Yahya: Introduction to Fluid-Structure Interactions, 2021, Springer Nature, 978-3-030-85882-7, http://dx.doi.org/10.1007/978-3-030-85884-1

Introduces the subject of Fluid-Structure Interactions (FSI) to students and professionals and discusses the major ideas in FSI with the goal of providing the fundamental understanding to the readers who possess limited or no understanding of the subject.

{{Physics-footer}}

{{DEFAULTSORT:Fluid-structure interaction}}

Category:Fluid mechanics

Category:Fluid dynamics