GLOBAL EXISTENCE FOR PHASE TRANSITION PROBLEMS VIA A VARIATIONAL SCHEME

2004 ◽  
Vol 01 (04) ◽  
pp. 747-768
Author(s):  
CHRISTIAN ROHDE ◽  
MAI DUC THANH
Keyword(s):  
Phase Transition ◽  
Surface Tension ◽  
Initial Data ◽  
Time Discretization ◽  
Approximate Solutions ◽  
Implicit Time ◽  
Time Step ◽  
Dynamical Phase Transition ◽  
Transition Dynamics ◽  
Variational Scheme

We construct approximate solutions of the initial value problem for dynamical phase transition problems via a variational scheme in one space dimension. First, we deal with a local model of phase transition dynamics which contains second and third order spatial derivatives modeling the effects of viscosity and surface tension. Assuming that the initial data are periodic, we prove the convergence of approximate solutions to a weak solution which satisfies the natural dissipation inequality. We note that this result still holds for non-periodic initial data. Second, we consider a model of phase transition dynamics with only Lipschitz continuous stress–strain function which contains a non-local convolution term to take account of surface tension. We also establish the existence of weak solutions. In both cases the proof relies on implicit time discretization and the analysis of a minimization problem at each time step.

A Fast Semi-Implicit Level Set Method for Curvature Dependent Flows with an Application to Limit Cycles Extraction in Dynamical Systems

2015 ◽  
Vol 18 (1) ◽  
pp. 203-229 ◽  
Author(s):  
Guoqiao You ◽  
Shingyu Leung
Keyword(s):  
Dynamical Systems ◽  
Level Set Method ◽  
Level Set ◽  
Limit Cycles ◽  
Time Discretization ◽  
Numerical Approach ◽  
Initial Guess ◽  
Implicit Time ◽  
Time Step ◽  
Image Regularization

AbstractWe propose a new semi-implicit level set approach to a class of curvature dependent flows. The method generalizes a recent algorithm proposed for the motion by mean curvature where the interface is updated by solving the Rudin-Osher-Fatemi (ROF) model for image regularization. Our proposal is general enough so that one can easily extend and apply the method to other curvature dependent motions. Since the derivation is based on a semi-implicit time discretization, this suggests that the numerical scheme is stable even using a time-step significantly larger than that of the corresponding explicit method. As an interesting application of the numerical approach, we propose a new variational approach for extracting limit cycles in dynamical systems. The resulting algorithm can automatically detect multiple limit cycles staying inside the initial guess with no condition imposed on the number nor the location of the limit cycles. Further, we also propose in this work an Eulerian approach based on the level set method to test if the limit cycles are stable or unstable.

A Conservative Local Discontinuous Galerkin Method for the Schrödinger-KdV System

2014 ◽  
Vol 15 (4) ◽  
pp. 1091-1107 ◽  
Author(s):  
Yinhua Xia ◽  
Yan Xu
Keyword(s):  
Discontinuous Galerkin ◽  
Time Discretization ◽  
Local Discontinuous Galerkin Method ◽  
Implicit Time ◽  
Discrete Version ◽  
Time Step ◽  
Local Discontinuous Galerkin ◽  
Ldg Method ◽  
Spatial Derivatives ◽  
Number Of Particles

AbstractIn this paper we develop a conservative local discontinuous Galerkin (LDG) method for the Schrödinger-Korteweg-de Vries (Sch-KdV) system, which arises in various physical contexts as a model for the interaction of long and short nonlinear waves. Conservative quantities in the discrete version of the number of plasmons, energy of the oscillations and the number of particles are proved for the LDG scheme of the Sch-KdV system. Semi-implicit time discretization is adopted to relax the time step constraint from the high order spatial derivatives. Numerical results for accuracy tests of stationary traveling soliton, and the collision of solitons are shown. Numerical experiments illustrate the accuracy and capability of the method.

A multi-phase Mullins–Sekerka system: matched asymptotic expansions and an implicit time discretisation for the geometric evolution problem

1998 ◽  
Vol 128 (3) ◽  
pp. 481-506 ◽  
Author(s):  
Lia Bronsard ◽  
Harald Garcke ◽  
Barbara Stoth
Keyword(s):  
Asymptotic Expansions ◽  
Approximate Solutions ◽  
Implicit Time ◽  
Triple Junctions ◽  
Sharp Interface Model ◽  
Equilibrium Equations ◽  
Time Step ◽  
Time Discretisation ◽  
Geometric Evolution ◽  
Multi Phase

We propose a generalisation of the Mullins–Sekerka problem to model phase separation in multi-component systems. The model includes equilibrium equations in bulk, the Gibbs–Thomson relation on the interfaces, Young's law at triple junctions, together with a dynamic law of Stefan type. Using formal asymptotic expansions, we establish the relationship to a transition layer model known as the Cahn-Hilliard system. We introduce a notion of weak solutions for this sharp interface model based on integration by parts on manifolds, together with measure theoretical tools. Through an implicit time discretisation, we construct approximate solutions by stepwise minimisation. Under the assumption that there is no loss of area as the time step tends to zero, we show the existence of a weak solution.

scholarly journals All Speed Scheme for the Low Mach Number Limit of the Isentropic Euler Equations

2011 ◽  
Vol 10 (1) ◽  
pp. 1-31 ◽  
Author(s):  
Pierre Degond ◽  
Min Tang
Keyword(s):  
Mach Number ◽  
Euler Equations ◽  
Computational Cost ◽  
Time Discretization ◽  
Implicit Time ◽  
Elliptic Type ◽  
Implicit Methods ◽  
Time Step ◽  
Low Mach Number ◽  
Isentropic Euler Equations

AbstractAn all speed scheme for the Isentropic Euler equations is presented in this paper. When the Mach number tends to zero, the compressible Euler equations converge to their incompressible counterpart, in which the density becomes a constant. Increasing approximation errors and severe stability constraints are the main difficulty in the low Mach regime. The key idea of our all speed scheme is the special semi-implicit time discretization, in which the low Mach number stiff term is divided into two parts, one being treated explicitly and the other one implicitly. Moreover, the flux of the density equation is also treated implicitly and an elliptic type equation is derived to obtain the density. In this way, the correct limit can be captured without requesting the mesh size and time step to be smaller than the Mach number. Compared with previous semi-implicit methods, firstly, nonphysical oscillations can be suppressed by choosing proper parameter, besides, only a linear elliptic equation needs to be solved implicitly which reduces much computational cost. We develop this semi-implicit time discretization in the framework of a first order Local Lax-Friedrichs (or Rusanov) scheme and numerical tests are displayed to demonstrate its performances.

scholarly journals Comparison of Surface Tension Models for the Volume of Fluid Method

Processes ◽  
2019 ◽  
Vol 7 (8) ◽  
pp. 542 ◽  
Author(s):  
Kurian J. Vachaparambil ◽  
Kristian Etienne Einarsrud
Keyword(s):  
Surface Tension ◽  
Multiphase Flow ◽  
Capillary Rise ◽  
Surface Tension Force ◽  
Surface Force ◽  
Parallel Plates ◽  
Time Step ◽  
Mesh Resolution ◽  
Active Research ◽  
Spurious Currents

With the increasing use of Computational Fluid Dynamics to investigate multiphase flow scenarios, modelling surface tension effects has been a topic of active research. A well known associated problem is the generation of spurious velocities (or currents), arising due to inaccuracies in calculations of the surface tension force. These spurious currents cause nonphysical flows which can adversely affect the predictive capability of these simulations. In this paper, we implement the Continuum Surface Force (CSF), Smoothed CSF and Sharp Surface Force (SSF) models in OpenFOAM. The models were validated for various multiphase flow scenarios for Capillary numbers of 10 − 3 –10. All the surface tension models provide reasonable agreement with benchmarking data for rising bubble simulations. Both CSF and SSF models successfully predicted the capillary rise between two parallel plates, but Smoothed CSF could not provide reliable results. The evolution of spurious current were studied for millimetre-sized stationary bubbles. The results shows that SSF and CSF models generate the least and most spurious currents, respectively. We also show that maximum time step, mesh resolution and the under-relaxation factor used in the simulations affect the magnitude of spurious currents.

Complex network structure of flocks in the Vicsek Model with Vectorial Noise

2014 ◽  
Vol 25 (03) ◽  
pp. 1350095 ◽  
Author(s):  
Gabriel Baglietto ◽  
Ezequiel V. Albano ◽  
Julián Candia
Keyword(s):  
Phase Transition ◽  
Complex Network ◽  
Network Structure ◽  
Time Step ◽  
Vicsek Model ◽  
Disorder Phase Transition ◽  
First Order ◽  
Direction Of Movement ◽  
Long Range Ordering ◽  
Influence Of Noise

In the Vicsek Model (VM), self-driven individuals try to adopt the direction of movement of their neighbors under the influence of noise, thus leading to a noise-driven order–disorder phase transition. By implementing the so-called Vectorial Noise (VN) variant of the VM (i.e. the VM-VN model), this phase transition has been shown to be discontinuous (first-order). In this paper, we perform an extensive complex network study of VM-VN flocks and show that their topology can be described as highly clustered, assortative, and nonhierarchical. We also study the behavior of the VM-VN model in the case of "frozen flocks" in which, after the flocks are formed using the full dynamics, particle displacements are suppressed (i.e. only rotations are allowed). Under this kind of restricted dynamics, we show that VM-VN flocks are unable to support the ordered phase. Therefore, we conclude that the particle displacements at every time-step in the VM-VN dynamics are a key element needed to sustain long-range ordering throughout.

Sign in / Sign up

Export Citation Format

Share Document