ترغب بنشر مسار تعليمي؟ اضغط هنا

An energy stable $C^0$ finite element scheme for a quasi-incompressible phase-field model of moving contact line with variable density

66   0   0.0 ( 0 )
 نشر من قبل Shixin Xu
 تاريخ النشر 2019
  مجال البحث الهندسة المعلوماتية
والبحث باللغة English




اسأل ChatGPT حول البحث

In this paper, we focus on modeling and simulation of two-phase flow with moving contact lines and variable density. A thermodynamically consistent phase-field model with General Navier Boundary Condition is developed based on the concept of quasi-incompressibility and the energy variational method. Then a mass conserving and energy stable C0 finite element scheme is developed to solve the PDE system. Various numerical simulation results show that the proposed schemes are mass conservative, energy stable and the 2nd order for P1 element and 3rd order for P2 element convergence rate in the sense of L2 norm.



قيم البحث

اقرأ أيضاً

A thermodynamically consistent phase-field model is introduced for simulating motion and shape transformation of vesicles under flow conditions. In particular, a general slip boundary condition is used to describe the interaction between vesicles and the wall of the fluid domain. A second-order accurate in both space and time C0 finite element method is proposed to solve the model governing equations. Various numerical tests confirm the convergence, energy stability, and conservation of mass and surface area of cells of the proposed scheme. Vesicles with different mechanical properties are also used to explain the pathological risk for patients with sickle cell disease.
In this article, we present and analyze a finite element numerical scheme for a three-component macromolecular microsphere composite (MMC) hydrogel model, which takes the form of a ternary Cahn-Hilliard-type equation with Flory-Huggins-deGennes energ y potential. The numerical approach is based on a convex-concave decomposition of the energy functional in multi-phase space, in which the logarithmic and the nonlinear surface diffusion terms are treated implicitly, while the concave expansive linear terms are explicitly updated. A mass lumped finite element spatial approximation is applied, to ensure the positivity of the phase variables. In turn, a positivity-preserving property can be theoretically justified for the proposed fully discrete numerical scheme. In addition, unconditional energy stability is established as well, which comes from the convexity analysis. Several numerical simulations are carried out to verify the accuracy and positivity-preserving property of the proposed scheme.
121 - Yifei Li , Weizhu Bao 2020
We propose an energy-stable parametric finite element method (ES-PFEM) to discretize the motion of a closed curve under surface diffusion with an anisotropic surface energy $gamma(theta)$ -- anisotropic surface diffusion -- in two dimensions, while $ theta$ is the angle between the outward unit normal vector and the vertical axis. By introducing a positive definite surface energy (density) matrix $G(theta)$, we present a new and simple variational formulation for the anisotropic surface diffusion and prove that it satisfies area/mass conservation and energy dissipation. The variational problem is discretized in space by the parametric finite element method and area/mass conservation and energy dissipation are established for the semi-discretization. Then the problem is further discretized in time by a (semi-implicit) backward Euler method so that only a linear system is to be solved at each time step for the full-discretization and thus it is efficient. We establish well-posedness of the full-discretization and identify some simple conditions on $gamma(theta)$ such that the full-discretization keeps energy dissipation and thus it is unconditionally energy-stable. Finally the ES-PFEM is applied to simulate solid-state dewetting of thin films with anisotropic surface energies, i.e. the motion of an open curve under anisotropic surface diffusion with proper boundary conditions at the two triple points moving along the horizontal substrate. Numerical results are reported to demonstrate the efficiency and accuracy as well as energy dissipation of the proposed ES-PFEM.
97 - Wenjia Liu , Shuo Zhang 2021
In this paper, we propose a finite element pair for incompressible Stokes problem. The pair uses a slightly enriched piecewise linear polynomial space for velocity and piecewise constant space for pressure, and is illustrated to be a lowest-degree co nservative stable pair for the Stokes problem on general triangulations.
We introduce a new method for the numerical approximation of time-harmonic acoustic scattering problems stemming from material inhomogeneities. The method works for any frequency $omega$, but is especially efficient for high-frequency problems. It is based on a time-domain approach and consists of three steps: emph{i)} computation of a suitable incoming plane wavelet with compact support in the propagation direction; emph{ii)} solving a scattering problem in the time domain for the incoming plane wavelet; emph{iii)} reconstruction of the time-harmonic solution from the time-domain solution via a Fourier transform in time. An essential ingredient of the new method is a front-tracking mesh adaptation algorithm for solving the problem in emph{ii)}. By exploiting the limited support of the wave front, this allows us to make the number of the required degrees of freedom to reach a given accuracy significantly less dependent on the frequency $omega$, as shown in the numerical experiments. We also present a new algorithm for computing the Fourier transform in emph{iii)} that exploits the reduced number of degrees of freedom corresponding to the adapted meshes.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا