You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('') and can be up to 35 characters long.
131 lines
6.8 KiB
131 lines
6.8 KiB


%% CMBE TEMPLATE FOR ABSTRACT SUBMISSION 

%% USE CLASS cmbe17.cls ** DO NOT MODIFY CLASS FILE ** 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% The extended abstract should consist of a short abstract (SUMMARY) 

% and 24 keywords, followed by a short article. The length of the 

% extended abstract must not be longer than four pages incl. references. 

% 

% Guidelines to complete each section are provided below. 

% References are managed using thebibliography with cite_key and 

% the \bibitem{cite_key} command. Make sure they are arranged by 

% order of appearance in the text. ** No page numbering. 

% 

% 

% Consult the conference website for submission information. 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 



\documentclass{cmbe21} 

% additional packages 

\usepackage{amsmath,amsfonts,amssymb} 

\usepackage{graphicx,wrapfig} 

% comment if not needed 

\usepackage[hyphens]{url} 

\usepackage{hyperref} 

\usepackage{color} 





% TITLE: replace text with your abstract title WITHOUT full stop 

\title{On monolithic and ChorinTemam schemes for incompressible flows in moving domains} 

% AUTHOR/AFFILIATION: handled by authblk. 

% Use only one of the two following methods for author listing. Delete or comment out the other. 

% Add/remove authors/affiliations as necessary, complete following the template without adding additional superscript/footnotes 



% 1 Authors have the same affiliation: 

% ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 

%\author{First A. Author} 

%\author{Second B. Author} 

%\author{Third C. Author} 

%\affil{Affiliation, Postal Address, \texttt{\{First,Second,Third\}@affil}} 



% 2 Multiple authors with multiple affiliations. Complete as follows: author[i] <> affiliation[i] 

% ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 

%\author[1]{Jerem\'ias Garay} 

%\author[2]{Second B. Author} 

%\author[3]{Third C. Author} 

%\author[2]{Fourth D. Author} 







% \author[1]{Jerem\'ias Garay} 

\author[1]{Reidmen Ar\'ostica} 

% \author[1]{David Nolte} 

\author[1]{Crist\'obal Bertoglio} 



\affil[1]{{Bernoulli Institute}, {University of Groningen}, 

{Groningen}, The Netherlands} 

%\affil[cmm]{{Center for Mathematical Modeling}, {Universidad de Chile}, {Santiago}, Chile} 

%\affil[tub]{{Department of Fluid Dynamics}, {Technische Universit\"at Berlin}, {Berlin}, Germany} 



%\affil[2]{{Bernoulli Institute}, {University of Groningen}, 

%{Groningen}, The Netherlands} 



%\affil[2]{{Biomedical Imaging Center}, {Pontificia Universidad Cat\'olica de Chile}, 

%{Santiago}, Chile} 





%\affil[3]{{School of Biomedical Engineering}, {Universidad de Valparaiso}, 

%{Valparaiso}, Chile} 

% 

%\affil[4]{{Department of Mathematical Engineering}, {Universidad de Concepci\'on}, 

%{Concepci\'on}, Chile} 

% 

% 

%\affil[5]{{Department of Mechanical Engineering}, {Universidad T\'ecnica Federico Santa Mar\'ia}, 

%{Santiago}, Chile} 

% 

%\affil[6]{{Joint last authors}, {in alphabetical order}} 









%\affil[1]{Affiliation 1, Postal Address, \texttt{First@affil1}} 

%\affil[2]{Affiliation 2, Postal Address, \texttt{\{Second,Fourth\}@affil2}} 

%\affil[3]{Affiliation 3, Postal Address, \texttt{Third@affil3}} 











% SUMMARY: replace text with a short summary 

\summary{Several time discretized domain for the incompressible NavierStokes equations (iNSE) in moving domains have been proposed in literature. Here, we introduce a unified formulation that combines different approaches found in literature, allowing a common well posedness and time stability analysis. It can be therefore shown that only a particular choice of numerical schemes ensure such properties under some restrictions. The analysis will be shown for ChorinTemam schemes using the insight found in the monolithic case. Results are supported from numerical simulations and its usage in fluidsolid interaction problems in cardiac geometries will be presented.} 



% KEYWORDS: replace text with 24 keywords, not capitalised, separated by comma, and without a full stop at the end. 

\keywords{numerical schemes, stability analysis, incompressible flows, fluidstructure interaction} 





\begin{document} 



Several works have been reported dealing with numerical solutions of the iNSE in moving domains within the Arbitrary Lagrangian Eulerian formulation (ALE), primarily in the context of fluidsolid coupling, e.g. \cite{astorinochoulyfernandez09, bertoglio2013sisc}. 



Different choices of time discretization schemes have been reported e.g. \cite{Basting2017, Hessenthaler2017}, nevertheless to the best of the authors knowledge, only a few monolithic schemes have been throughly analyzed, e.g. \cite{Lozovskiy2018, smaldone2014, letallecmouro01, Burtschell2017} while no analysis has been reported for ChorinTemam (CT) schemes, being a feasible alternative when requirements such a low time computations are needed. 



The goal of this talk is to present the finding of wellposedness and unconditional energy balance of the iNSEALE for several reported monolithic discretization schemes within a single formulation, recently published in \cite{arostica2021monolithic}. The main result is that under appropriate conditions on the rate of domain deformation, only some of first order time discretization schemes are unconditionally stable. 



We will show the extension to a CT scheme. Namely, in that case the following inequality can be shown: 

\begin{equation} 

\int_{\Omega^0} \frac{\rho J^{n+1}}{2\tau} \vert \tilde{\mathbf{u}}^{n+1} \vert^2 \, \text{d}\mathbf{X} \int_{\Omega^0} \frac{\rho J^{n}}{2\tau} \vert \tilde{\mathbf{u}}^{n} \vert^2 \, \text{d}\mathbf{X} \leq  \int_{\Omega^0} J^{\star} 2 \mu \vert \epsilon^{\star} (\tilde{\mathbf{u}}) \vert^2 \, \text{d}\mathbf{X}  \int_{\Omega^0} \frac{\tau J^n}{2 \rho} \vert Grad(p^n) H^n \vert^2 \, \text{d} \mathbf{X} 

\end{equation} 

for $(\mathbf{u}^n, p^n)$ the velocity/pressure pair solution at time $t^n$, in the reference domain $\Omega^0$, for operators to be specified in the talk. 



Our finding will be supplemented with an application to fluidsolid interaction in an idealized cardiac geometry, exploiting the splitting nature of the CT scheme with a wellknown coupling approach \cite{bertoglio2013sisc, fernandezgerbeaugrandmont06}, see example in Figure \ref{fig:comparison_figure}. 

%In such a case, we will exploit the fluid pressure projection step, coupling it with the solid problem in an efficient fashion. Simulations of such a case will be provided, as well as the current research done in more realistic geometries. 

\begin{figure}[!hbtp] 

\centering 

\includegraphics[width=0.8\textwidth]{figs/comparison_two_ways_to_one_way.png} 

\caption{Twoways FSI (left) and oneway FSI (right), on an ellipsoid. In arrows the fluid velocity magnitude and direction, driven by an hyperelastic active solid (in grey, with decreased opacity).} 

\label{fig:comparison_figure} 

\end{figure} 



\newpage 





\bibliography{biblio_merged.bib} 

\bibliographystyle{unsrt} 



\end{document} 









