isogeometric analysis of the cahn hilliard equation a convergence study herein we present a numerical convergence study of the cahn hilliard phase field model within an isogeometric finite element analysis framework using a manufactured solution a mixed formulation of the cahn hilliard equation and the direct discretisation of the weak form which requires a c 1 continuous approximation are compared in terms of convergence rates for approximations that are higher than second order in space the direct discretisation is found to be superior suboptimal convergence rates occur when splines of order p 2 are used this is validated with a priori error estimates for linear problems the convergence analysis is completed with an investigation of the temporal discretisation second order accuracy is found for the generalised method this ensures the functionality of an adaptive time stepping scheme which is required for the efficient numerical solution of the cahn hilliard equation the isogeometric finite element framework is eventually validated by two numerical examples of spinodal decomposition

conservative nonlinear difference scheme for the cahn hilliard equation ii abstract numerical solutions for the cahn hilliard equation is considered using the crank nicolson type finite difference method existence of the solution for the difference scheme has been shown by brouwer fixed point theorem stability convergence and error analysis of the scheme are shown we also show that the scheme preserves the discrete mass even though the linearized scheme in 1 is conditionally stable and does not preserve the mass

a quantitative comparison between c0 and c1 elements for solving the cahn hilliard equation the cahn hilliard ch equation is a time dependent fourth order partial differential equation pde when solving the ch equation via the finite element method fem the domain is discretized by c 1 continuous basis functions or the equation is split into a pair of second order pdes and discretized via c 0 continuous basis functions in the current work a quantitative comparison between c 1 hermite and c 0 lagrange elements is carried out using a continuous galerkin fem formulation the different discretizations are evaluated using the method of manufactured solutions solved with newton s method and jacobian free newton krylov it is found that the use of linear lagrange elements provides the fastest computation time for a given number of elements while the use of cubic hermite elements provides the lowest error the results offer a set of benchmarks to consider when choosing basis functions to solve the ch equation in addition an example of microstructure evolution demonstrates the different types of elements for a traditional phase field model

isogeometric analysis of the advective cahn hilliard equation spinodal decomposition under shear flow we present a numerical study of the spinodal decomposition of a binary fluid undergoing shear flow using the advective cahn hilliard equation a stiff nonlinear parabolic equation characterized by the presence of fourth order spatial derivatives our numerical solution procedure is based on isogeometric analysis an approximation technique for which basis functions of high order continuity are employed these basis functions allow us to directly discretize the advective cahn hilliard equation without resorting to a mixed formulation we present steady state solutions for rectangular domains in two dimensions and for the first time in three dimensions we also present steady state solutions for the two dimensional taylor couette cell to enforce periodic boundary conditions in this curved domain we derive and utilize a new periodic bezier extraction operator we present an extensive numerical study showing the effects of shear rate surface tension and the geometry of the domain on the phase evolution of the binary fluid theoretical and experimental results are compared with our simulations

finite element approximation of the cahn hilliard equation with degenerate mobility we consider a fully practical finite element approximation of the cahn hilliard equation with degenerate mobility textstyle frac partial u partial t del b u del gamma lap u psi u where b cdot geq 0 is a diffusional mobility and psi cdot is a homogeneous free energy in addition to showing well posedness and stability bounds for our approximation we prove convergence in one space dimension furthermore an iterative scheme for solving the resulting nonlinear discrete system is analyzed we also discuss how our approximation has to be modified in order to be applicable to a logarithmic homogeneous free energy finally some numerical experiments are presented

thermodynamically consistent phase field models of fracture variational principles and multi field fe implementations the computational modeling of failure mechanisms in solids due to fracture based on sharp crack discontinuities suffers in situations with complex crack topologies this can be overcome by a diffusive crack modeling based on the introduction of a crack phase field in this paper we outline a thermodynamically consistent framework for phase field models of crack propagation in elastic solids develop incremental variational principles and consider their numerical implementations by multi field finite element methods we start our investigation with an intuitive and descriptive derivation of a regularized crack surface functional that converges for vanishing length scale parameter to a sharp crack topology functional this functional provides the basis for the definition of suitable convex dissipation functions that govern the evolution of the crack phase field here we propose alternative rate independent and viscous over force models that ensure the local growth of the phase field next we define an energy storage function whose positive tensile part degrades with increasing phase field with these constitutive functionals at hand we derive the coupled balances of quasi static stress equilibrium and gradient type phase field evolution in the solid from the argument of virtual power here we consider a canonical two field setting for rate independent response and a time regularized three field formulation with viscous over force response it is then shown that these balances follow as the euler equations of incremental variational principles that govern the multi field problems these principles make the proposed formulation extremely compact and provide a perfect base for the finite element implementation including features such as the symmetry of the monolithic tangent matrices we demonstrate the performance of the proposed phase field formulations of fracture by means of representative numerical examples copyright 2010 john wiley sons ltd

a second order splitting method for the cahn hilliard equation a semi discrete finite element method requiring only continuous element is presented for the approximation of the solution of the evolutionary fourth order in space cahn hilliard equation optimal order error bounds are derived in various norms for an implementation which uses mass lumping the continuous problem has an energy based lyapunov functional it is proved that this property holds for the discrete problem

using angular dispersion of gradient direction for detecting edge ribbons edges in a scene generally project to smooth continuous curves nearly everywhere in the image which results in low angular deviation of the intensity gradients in small neighborhoods straddling edges angular deviation is shown to be a measure of snr a theoretical analysis of angular deviation arising due to independent and identically distributed n 0 sigma sup 2 random noise is presented angular deviation thresholds for neighborhood sizes from 3 3 to 11 11 pixels are determined both from this analysis and numerical examples the proposed gradient angular dispersion detection algorithm detects edge elements by comparing the measured angular deviation with values computed for the minimum acceptable snr low values of deviation violate the no edge hypothesis the algorithm is shown to make good use of the limited dynamic range of the imaging system the sensitivity and selectivity of the strategy are both shown to be high

bounded and periodic solutions of differential equations in banach space abstract we study the equation u t au t f t t 0 1 in a banach space x with a the generator of an analytic or a strongly continuous semigroup s and prove that if solutions of 1 are bounded and ultimate bounded with f t periodic and s t compact then 1 has a t periodic solution we also show that the existence of a proper liapunov function implies the boundedness and ultimate boundedness of solutions of 1 these results extend earlier results in finite dimensional spaces we then apply the results to a parabolic partial differential equation u t t x 2m c x d u t x f t x 2

comparison of high order finite element and discontinuous galerkin methods for phase field equations application to structural damage abstract phase field equations are used to model a wide range of multiphase problems such as separation of fluids solidification viscous fingering fracture and fatigue a wide variety of methods to numerically solve phase field equations can be found in the literature in particular high order methods are an effective option when accuracy improvement is desired in the first part of this work we analyze the accuracy and computational efficiency of the high order finite element method fem and discontinuous galerkin dg method applied to the second order allen cahn ac and fourth order cahn hilliard ch equations several schemes for time integration are used for these equations the explicit schemes are the forward euler classical fourth order runge kutta rk4 and the strong stability preserving ten stages fourth order runge kutta rkssp 10 4 described in gottlieb et al 2011 the backward euler and trapezoidal implicit methods are adopted in the full and semi implicit schemes as proposed in eyre unpublished manufactured solutions for one dimensional problems are used in order to evaluate the errors and to compare the different numerical methods by choosing an adequate discretization for ac equations resulting from the previous analysis of the first part of the work in the second part we propose a numerical semi implicit scheme to solve the damage and fracture model described in boldrini et al 2016 this procedure employs the fem for spatial discretization the newmark method for time integration of the kinematics equation and the backward euler for the damage phase field evolution finally results for 2d benchmark tests are presented for the fracture phase field model and the convergence to a sharp crack for a small width of the damage phase field layer is verified

boundedness and periodicity in impulsive ordinary and functional differential equations abstract in this paper we study the boundedness and periodicity of solutions for impulsive ordinary differential equations odes and functional differential equations fdes by using horn s fixed point theorem hale yoshizawa type criteria for the existence of t periodic solutions are established to use these criteria we also give two new boundedness theorems and establish new existence results for t periodic solutions which show that the impulsive perturbations do contribute to yielding periodic solutions even when the underlying systems do not enjoy periodic solutions

a higher order phase field model for brittle fracture formulation and analysis within the isogeometric analysis framework phase field models based on the variational formulation for brittle fracture have recently been gaining popularity these models have proven capable of accurately and robustly predicting complex crack behavior in both two and three dimensions in this work we propose a fourth order model for the phase field approximation of the variational formulation for brittle fracture we derive the thermodynamically consistent governing equations for the fourth order phase field model by way of a variational principle based on energy balance assumptions the resulting model leads to higher regularity in the exact phase field solution which can be exploited by the smooth spline function spaces utilized in isogeometric analysis this increased regularity improves the convergence rate of the numerical solution and opens the door to higher order convergence rates for fracture problems we present an analysis of our proposed theory and numerical examples that support this claim we also demonstrate the robustness of the model in capturing complex three dimensional crack behavior

precision edge contrast and orientation estimation the contrast and orientation estimation accuracy of several edge operators that have been proposed in the literature is examined both for the noiseless case and in the presence of additive gaussian noise the test image is an ideal step edge that has been sampled with a square aperture grid the effects of subpixel translations and rotations of the edge on the performance of the operators are studied it is shown that the effect of subpixel translations of an edge can generate more error than moderate noise levels methods with improved results are presented for sobel angle estimates and the nevatia babu operator and theoretical noise performance evaluations are also provided an edge operator based on two dimensional spatial moments is presented all methods are compared according to worst case and rms error in an ideal noiseless situation and rms error under various noise levels

isogeometric analysis of the cahn hilliard phase field model abstract the cahn hilliard equation involves fourth order spatial derivatives finite element solutions are not common because primal variational formulations of fourth order operators are only well defined and integrable if the finite element basis functions are piecewise smooth and globally c 1 continuous there are a very limited number of two dimensional finite elements possessing c 1 continuity applicable to complex geometries but none in three dimensions we propose isogeometric analysis as a technology that possesses a unique combination of attributes for complex problems involving higher order differential operators namely higher order accuracy robustness two and three dimensional geometric flexibility compact support and most importantly the possibility of c 1 and higher order continuity a nurbs based variational formulation for the cahn hilliard equation was tested on two and three dimensional problems we present steady state solutions in two dimensions and for the first time in three dimensions to achieve these results an adaptive time stepping method is introduced we also present a technique for desensitizing calculations to dependence on mesh refinement this enables the calculation of topologically correct solutions on coarse meshes opening the way to practical engineering applications of phase field methodology

robust and efficient image alignment based on relative gradient matching in this paper we present a robust image alignment algorithm based on matching of relative gradient maps this algorithm consists of two stages namely a learning based approximate pattern search and an iterative energy minimization procedure for matching relative image gradient the first stage finds some candidate poses of the pattern from the image through a fast nearest neighbor search of the best match of the relative gradient features computed from training database of feature vectors which are obtained from the synthesis of the geometrically transformed template image with the transformation parameters uniformly sampled from a given transformation parameter space subsequently the candidate poses are further verified and refined by matching the relative gradient images through an iterative energy minimization procedure this approach based on the matching of relative gradients is robust against nonuniform illumination variations experimental results on both simulated and real images are shown to demonstrate superior efficiency and robustness of the proposed algorithm over the conventional normalized correlation method

robust fft based scale invariant image registration with image gradients we present a robust fft based approach to scale invariant image registration our method relies on fft based correlation twice once in the log polar fourier domain to estimate the scaling and rotation and once in the spatial domain to recover the residual translation previous methods based on the same principles are not robust to equip our scheme with robustness and accuracy we introduce modifications which tailor the method to the nature of images first we derive efficient log polar fourier representations by replacing image functions with complex gray level edge maps we show that this representation both captures the structure of salient image features and circumvents problems related to the low pass nature of images interpolation errors border effects and aliasing second to recover the unknown parameters we introduce the normalized gradient correlation we show that using image gradients to perform correlation the errors induced by outliers are mapped to a uniform distribution for which our normalized gradient correlation features robust performance exhaustive experimentation with real images showed that unlike any other fourier based correlation techniques the proposed method was able to estimate translations arbitrary rotations and scale factors up to 6

computationally efficient solution to the cahn hilliard equation adaptive implicit time schemes mesh sensitivity analysis and the 3d isoperimetric problem we present an efficient numerical framework for analyzing spinodal decomposition described by the cahn hilliard equation we focus on the analysis of various implicit time schemes for two and three dimensional problems we demonstrate that significant computational gains can be obtained by applying embedded higher order runge kutta methods in a time adaptive setting this allows accessing time scales that vary by five orders of magnitude in addition we also formulate a set of test problems that isolate each of the sub processes involved in spinodal decomposition interface creation and bulky phase coarsening we analyze the error fluctuations using these test problems on the split form of the cahn hilliard equation solved using the finite element method with basis functions of different orders any scheme that ensures at least four elements per interface satisfactorily captures both sub processes our findings show that linear basis functions have superior error to cost properties this strategy coupled with a domain decomposition based parallel implementation let us notably augment the efficiency of a numerical cahn hillard solver and open new venues for its practical applications especially when three dimensional problems are considered we use this framework to address the isoperimetric problem of identifying local solutions in the periodic cube in three dimensions the framework is able to generate all five hypothesized candidates for the local solution of periodic isoperimetric problem in 3d sphere cylinder lamella doubly periodic surface with genus two lawson surface and triply periodic minimal surface p schwarz surface

a phase field description of dynamic brittle fracture abstract in contrast to discrete descriptions of fracture phase field descriptions do not require numerical tracking of discontinuities in the displacement field this greatly reduces implementation complexity in this work we extend a phase field model for quasi static brittle fracture to the dynamic case we introduce a phase field approximation to the lagrangian for discrete fracture problems and derive the coupled system of equations that govern the motion of the body and evolution of the phase field we study the behavior of the model in one dimension and show how it influences material properties for the temporal discretization of the equations of motion we present both a monolithic and staggered time integration scheme we study the behavior of the dynamic model by performing a number of two and three dimensional numerical experiments we also introduce a local adaptive refinement strategy and study its performance in the context of locally refined t splines we show that the combination of the phase field model and local adaptive refinement provides an effective method for simulating fracture in three dimensions

a phase field model for rate independent crack propagation robust algorithmic implementation based on operator splits the computational modeling of failure mechanisms in solids due to fracture based on sharp crack discontinuities suffers in situations with complex crack topologies this can be overcome by a diffusive crack modeling based on the introduction of a crack phase field following our recent work c miehe f welschinger m hofacker thermodynamically consistent phase field models of fracture variational principles and multi field fe implementations international journal for numerical methods in engineering doi 10 1002 nme 2861 on phase field type fracture we propose in this paper a new variational framework for rate independent diffusive fracture that bases on the introduction of a local history field it contains a maximum reference energy obtained in the deformation history which may be considered as a measure for the maximum tensile strain obtained in history it is shown that this local variable drives the evolution of the crack phase field the introduction of the history field provides a very transparent representation of the balance equation that governs the diffusive crack topology in particular it allows for the construction of a new algorithmic treatment of diffusive fracture here we propose an extremely robust operator split scheme that successively updates in a typical time step the history field the crack phase field and finally the displacement field a regularization based on a viscous crack resistance that even enhances the robustness of the algorithm may easily be added the proposed algorithm is considered to be the canonically simple scheme for the treatment of diffusive fracture in elastic solids we demonstrate the performance of the phase field formulation of fracture by means of representative numerical examples

the periodic solutions for general periodic impulsive population systems of functional differential equations and its applications in this paper the general periodic impulsive population systems of functional differential equations are investigated by using the method of poincare map and horn s fixed point theorem we prove that the ultimate boundedness of all solutions implies the existence of periodic solutions as applications of this result the existence of positive periodic solutions for the general periodic impulsive kolmogorov type population dynamical systems are discussed we further prove that as long as the system is permanent there must exist at least one positive periodic solution in addition the permanence and existence of positive periodic solutions are discussed for the periodic impulsive single species logistic models and the periodic impulsive n species lotka volterra competitive models with delays

new razumikhin type theorems for impulsive functional differential equations in this paper we investigate stability and boundedness for a class of impulsive functional differential equations by using liapunov functions and razumikhin technique some new razumikhin type theorems on stability and boundedness are obtained the corresponding results obtained recently in j shen j yan nonlinear anal 33 1998 519 531 x liu j shen to appear are improved

accurate efficient and iso geometrically flexible collocation methods for phase field models we propose new collocation methods for phase field models our algorithms are based on isogeometric analysis a new technology that makes use of functions from computational geometry such as for example non uniform rational b splines nurbs nurbs exhibit excellent approximability and controllable global smoothness and can represent exactly most geometries encapsulated in computer aided design cad models these attributes permitted us to derive accurate efficient and geometrically flexible collocation methods for phase field models the performance of our method is demonstrated by several numerical examples of phase separation modeled by the cahn hilliard equation we feel that our method successfully combines the geometrical flexibility of finite elements with the accuracy and simplicity of pseudo spectral collocation methods and is a viable alternative to classical collocation methods

a time adaptive finite volume method for the cahn hilliard and kuramoto sivashinsky equations this paper presents a complete finite volume method for the cahn hilliard and kuramoto sivashinsky type of equations the spatial discretization is high order accurate and suitable for general unstructured grids the time integration is addressed by means of implicit an implicit explicit fourth order runge kutta schemes with error control and adaptive time stepping the outcome is a practical accurate and efficient simulation tool which has been successfully applied to accuracy tests and representative simulations the use of adaptive time stepping is of paramount importance in problems governed by the cahn hilliard model an adaptive method may be several orders of magnitude more efficient than schemes using constant or heuristic time steps in addition to driving the simulations efficiently the time adaptive procedure provides a quantitative not just qualitative characterization of the rich temporal scales present in phase separation processes governed by the cahn hilliard phase field model

inpainting of binary images using the cahn hilliard equation image inpainting is the filling in of missing or damaged regions of images using information from surrounding areas we outline here the use of a model for binary inpainting based on the cahn hilliard equation which allows for fast efficient inpainting of degraded text as well as super resolution of high contrast images

continuous discontinuous finite element approximations of fourth order elliptic problems in structural and continuum mechanics with applications to thin beams and plates and strain gradient elasticity abstract a new finite element method for fourth order elliptic partial differential equations is presented and applied to thin bending theory problems in structural mechanics and to a strain gradient theory problem the method combines concepts from the continuous galerkin cg method the discontinuous galerkin dg method and stabilization techniques a brief review of the cg method the dg method and stabilization techniques highlights the advantages and disadvantages of these methods and suggests a new approach for the solution of fourth order elliptic problems a continuous discontinuous galerkin c dg method is proposed which uses c 0 continuous interpolation functions and is formulated in the primary variable only the advantage of this formulation over a more traditional mixed approach is that the introduction of additional unknowns and related difficulties can be avoided in the context of thin bending theory the c dg method leads to a formulation where displacements are the only degrees of freedom and no rotational degrees of freedom need to be considered the main feature of the c dg method is the weak enforcement of continuity of first and higher order derivatives through stabilizing terms on interior boundaries consistency stability and convergence of the method are shown analytically numerical experiments verify the theoretical results and applications are presented for bernoulli euler beam bending poisson kirchhoff plate bending and a shear layer problem using toupin mindlin strain gradient theory

a new approach to practical stability of impulsive functional differential equations in terms of two measures in this work we consider a new approach to the practical stability theory of impulsive functional differential equations with lyapunov functionals and razumikhin technique we use a new technique in the division of lyapunov functions given by shunian zhang and obtain conditions sufficient for the uniform practical asymptotical stability of impulsive delay differential equations an example is also discussed to illustrate the advantage of the proposed results

gradient orientation pattern matching with the hamming distance abstract this paper presents a novel pattern matching technique that is robust to illumination changes and the occlusion problem the technique is based on the matching of gradient orientations in place of traditional image features such as intensities or gradients gradient orientations depend on the texture in an image they are known to be insensitive to changes of image intensities that are often caused by time varying illuminations or the auto gain control agc function of the camera moreover the proposed method employs a voting strategy in the process of matching gradient orientations the method works remarkably well even when a large part of the pattern is occluded with a foreign object consequently the proposed method is robust to both irregular lighting conditions and the occlusion problem

impulsive differential equations this paper deals with the periodic solutions problem for impulsive differential equations by using lyapunov s second method and the contraction mapping principle some conditions ensuring the existence and global attractiveness of unique periodic solutions are derived which are given from impulsive control and impulsive perturbation points of view as an application the existence and global attractiveness of unique periodic solutions for hopfield neural networks are discussed finally two numerical examples are provided to demonstrate the effectiveness of the proposed results

hybrid image matching combining hausdorff distance with normalized gradient matching image matching has been a central problem in computer vision and image processing for decades most of the previous approaches to image matching can be categorized into the intensity based and edge based comparison hausdorff distance has been widely used for comparing point sets or edge maps since it does not require point correspondences in this paper we propose a new image similarity measure combining the hausdorff distance with a normalized gradient consistency score for image matching the normalized gradient consistency score is designed to compare the normalized image gradient fields between two images to alleviate the illumination variation problem in image matching by combining the edge based and intensity based information for image matching we are able to achieve robust image matching under different lighting conditions we show the superior robustness property of the proposed image matching technique through experiments on face recognition under different lighting conditions

a discontinuous galerkin method for the cahn hilliard equation a discontinuous galerkin finite element method has been developed to treat the high order spatial derivatives appearing in the cahn hilliard equation the cahn hilliard equation is a fourth order nonlinear parabolic partial differential equation originally proposed to model phase segregation of binary alloys the developed discontinuous galerkin approach avoids the need for mixed finite element methods coupled equations or interpolation functions with a high degree of continuity that have been employed in the literature to treat the fourth order spatial derivatives the variational formulation of the discontinuous galerkin method its implementation and numerical examples are presented in this communication it is also shown under what conditions the method is stable and an error estimate in an energy type norm is presented the method is evaluated by comparison with a standard finite element treatment in which the cahn hilliard equation is decomposed into two coupled partial differential equations

local discontinuous galerkin methods for the cahn hilliard type equations in this paper we develop local discontinuous galerkin ldg methods for the fourth order nonlinear cahn hilliard equation and system the energy stability of the ldg methods is proved for the general nonlinear case numerical examples for the cahn hilliard equation and the cahn hilliard system in one and two dimensions are presented and the numerical results illustrate the accuracy and capability of the methods