
Williams-Beuren syndrome (WBS) is a rare genetic disorder characterized by congenital heart defects, dysmorphic features, intellectual delay, and a distinctive social behavioral profile. This highly recurrent and homogeneous phenotype has been curiously reported to be associated with autism spectrum disorders (ASD). Both genetic and environmental origins have been implicated. This study aimed to describe Tunisian patients associated with WBS and ASD and explore the underlying etiologies.
Thirty-one clinically suspected WBS were referred for genetic exploration. A comprehensive evaluation using karyotyping, fluorescence in situ hybridization (FISH), and array comparative genomic hybridization (array-CGH) was performed.
All patients were clinically diagnosed and confirmed to have WBS through karyotyping and FISH analysis. Notably, six patients with complex or atypical clinical presentations underwent array-CGH. Two of these patients presented with ASD. Array CGH showed microdeletions ranging from 1.4 to 1.7 Mb in the 7q11.2 region. Further analysis of the extended region deletion identified a gene closely located in the deleted region, the HIP1 gene, involved in the central nervous system trafficking protein.
The recurrent deletion in WBS, as well as the mirror duplication, may contribute to ASD development in some cases, suggesting a potential involvement of the ASD genes pathway in this region. However, recessive genetic origins should also be considered, particularly in consanguineous families. Furthermore, our findings highlight the potential role of genetic factors and regulatory elements within the deleted region in modulating gene expression, notably the HIP1 gene. This underscores the implications of gene dosage and environmental factors in the broader WBS region, notably with language and social development.
The presence of ASD in WBS patients emphasizes the need to investigate all WBS patients for autistic traits to establish a better genotype–phenotype correlation. We underline the utility of array-CGH as a valuable genetic diagnostic tool for characterizing WBS cases, and we shed light on the complex interplay of behavioral disorders in the 7q11.2 region rearrangements.
Citation: Rim Khelifi, Afef Jelloul, Houda Ajmi, Wafa Slimani, Sarra Dimassi, Khouloud Rjiba, Manel Dardour, Moez Gribaa, Ali Saad, Soumaya Mougou-Zerelli. Toward autism spectrum disorders and Williams-Beuren syndrome co-occurrence condition in Tunisian patients: Genetic insights[J]. AIMS Molecular Science, 2024, 11(4): 379-394. doi: 10.3934/molsci.2024023
[1] | Meiting Qu, Shaohui Liu, Lei Li . Urban public health spatial planning using big data technology and visual communication in IoT. Mathematical Biosciences and Engineering, 2023, 20(5): 8583-8600. doi: 10.3934/mbe.2023377 |
[2] | Gayathri Vivekanandhan, Mahtab Mehrabbeik, Karthikeyan Rajagopal, Sajad Jafari, Stephen G. Lomber, Yaser Merrikhi . Applying machine learning techniques to detect the deployment of spatial working memory from the spiking activity of MT neurons. Mathematical Biosciences and Engineering, 2023, 20(2): 3216-3236. doi: 10.3934/mbe.2023151 |
[3] | Song Liu, Gusong Luo, Yonglong Cai, Wenjie Wu, Weitao Liu, Rong Zou, Wenxuan Tan . Determinants of consumer intention to adopt a self-service technology strategy for last-mile delivery in Guangzhou, China. Mathematical Biosciences and Engineering, 2024, 21(2): 3262-3280. doi: 10.3934/mbe.2024144 |
[4] | Jun Zhai, Bilin Xu . Research on meme transmission based on individual heterogeneity. Mathematical Biosciences and Engineering, 2021, 18(5): 5176-5193. doi: 10.3934/mbe.2021263 |
[5] | Xiang Wang, Yongcheng Wang, Limin He . An intelligent data analysis-based medical management method for lower limb health of football athletes. Mathematical Biosciences and Engineering, 2023, 20(8): 14005-14022. doi: 10.3934/mbe.2023624 |
[6] | Massimo Fioranelli, O. Eze Aru, Maria Grazia Roccia, Aroonkumar Beesham, Dana Flavin . A model for analyzing evolutions of neurons by using EEG waves. Mathematical Biosciences and Engineering, 2022, 19(12): 12936-12949. doi: 10.3934/mbe.2022604 |
[7] | Chrysovalantis Malesios, Nikoleta Jones, Alfie Begley, James McGinlay . Methodological approaches to exploring the spatial variation in social impacts of protected areas: An intercomparison of Bayesian regression modeling approaches and potential implications. Mathematical Biosciences and Engineering, 2024, 21(3): 3816-3837. doi: 10.3934/mbe.2024170 |
[8] | B. Spagnolo, D. Valenti, A. Fiasconaro . Noise in ecosystems: A short review. Mathematical Biosciences and Engineering, 2004, 1(1): 185-211. doi: 10.3934/mbe.2004.1.185 |
[9] | Zhipeng Ding, Hongxia Yun, Enze Li . A multimedia knowledge discovery-based optimal scheduling approach considering visual behavior in smart education. Mathematical Biosciences and Engineering, 2023, 20(3): 5901-5916. doi: 10.3934/mbe.2023254 |
[10] | Juncai Yao, Jing Shen, Congying Yao . Image quality assessment based on the perceived structural similarity index of an image. Mathematical Biosciences and Engineering, 2023, 20(5): 9385-9409. doi: 10.3934/mbe.2023412 |
Williams-Beuren syndrome (WBS) is a rare genetic disorder characterized by congenital heart defects, dysmorphic features, intellectual delay, and a distinctive social behavioral profile. This highly recurrent and homogeneous phenotype has been curiously reported to be associated with autism spectrum disorders (ASD). Both genetic and environmental origins have been implicated. This study aimed to describe Tunisian patients associated with WBS and ASD and explore the underlying etiologies.
Thirty-one clinically suspected WBS were referred for genetic exploration. A comprehensive evaluation using karyotyping, fluorescence in situ hybridization (FISH), and array comparative genomic hybridization (array-CGH) was performed.
All patients were clinically diagnosed and confirmed to have WBS through karyotyping and FISH analysis. Notably, six patients with complex or atypical clinical presentations underwent array-CGH. Two of these patients presented with ASD. Array CGH showed microdeletions ranging from 1.4 to 1.7 Mb in the 7q11.2 region. Further analysis of the extended region deletion identified a gene closely located in the deleted region, the HIP1 gene, involved in the central nervous system trafficking protein.
The recurrent deletion in WBS, as well as the mirror duplication, may contribute to ASD development in some cases, suggesting a potential involvement of the ASD genes pathway in this region. However, recessive genetic origins should also be considered, particularly in consanguineous families. Furthermore, our findings highlight the potential role of genetic factors and regulatory elements within the deleted region in modulating gene expression, notably the HIP1 gene. This underscores the implications of gene dosage and environmental factors in the broader WBS region, notably with language and social development.
The presence of ASD in WBS patients emphasizes the need to investigate all WBS patients for autistic traits to establish a better genotype–phenotype correlation. We underline the utility of array-CGH as a valuable genetic diagnostic tool for characterizing WBS cases, and we shed light on the complex interplay of behavioral disorders in the 7q11.2 region rearrangements.
The Cahn-Hilliard equation is solved in this study using the adaptive finite element approach.
{ut=−εΔ2u+1εΔf(u)in Ω×(0,T],∂nu=∂n(−εΔu+1εf(u))=0,on ∂Ω×(0,T],u(x,0)=u0in Ω. | (1.1) |
In the equation above, we use to Ω signify a limited domain in Rd(d=2 or 3) with Lipschitz boundary ∂Ω, and n stands for the unit outward normal to Ω, u also stands for the phase-field variable that will be solved. f(u) is a nonlinear function. ε is a constant. T is the ultimate time.
The Cahn-Hilliard (CH) equation, often known as the common phase field model, was first published in 1958 in a key paper by Cahn and Hillard's [2], which examined the thermodynamic phenomenon of mutual diffusion of two substances (such as alloys and polymers). In general, this equation can be used to express the intricate phase separation and coarsening phenomena in solids, especially in materials science and fluid dynamics, for examples, see [3,4,5,6] and the references therein. The Cahn-Hilliard equation, which holds multiple time scales and spatial scales, is a rigid and nonlinear fourth-order PDE. So it is challenging to identify the precise solution.
The majority of studies in recent years have focused on finite difference techniques or the Fourier-spectral approach (periodic boundary), and many authors have devoted their time to examining different Cahn-Hilliard equation versions, including the viscous Cahn-Hilliard equation, the surface Cahn-Hilliard equation, the CHNS equation, etc. and research issues from space and time that are related. Barrett and Elliott presented to study the Cahn-Hilliard system by using conforming and nonconforming finite element methods [7,8,9]. Shen and Xu described in [10] the treatment of the Cahn-Hilliard equation by the scalar auxiliary variable (SAV) method. Zhao and Xiao combined the surface finite element method (SFEM) with stabilized semi-implicit model to solve the surface Cahn-Hilliard model [11]. There are many parallels between the Allen-Cahn equation and the Cahn-Hilliard equation since they are phase field flows of the same energy in different regions, such as energy decreasing. In [12], Shen and Yang analyzed the discrete versions of the Allen-Cahn and Cahn-Hilliard equations. Numerous research has been conducted on the Allen-Cahn equation [13,14,15], and Chen and Huang applied the Superconvergent Cluster Recovery (SCR) approach for the Allen-Cahn equation [16]. The method was also used in [17] to solve the CHNS equation, but the Cahn-Hilliard equation has not yet seen this treatment. In [17], the Cahn-Hilliard component is treated as a first-order format combing with a time-space adaptive method. In our paper, we analyze the Cahn-Hilliard equation using a second-order time discrete method and convex splitting for the nonlinear factor f, combined with the SCR method.
There are two reasons why we use adaptive methods in this study. First, due to the presence of the tiny parameter ε, it produces the phenomenon of the interfacial layer. When meshing on this interfacial layer, different from other parts, we need fine encryption, and the adaptive strategy can effectively solve this problem. Second, it takes a very long time for the numerical simulation of the Cahn-Hilliard model to reach the steady state, therefore using adaptive mesh creation makes sense in terms of time and money savings. At present, numerous writers have examined the Cahn-Hilliard equation's adaptive technique, including [11], where the discrete chemical potential is employed as an error estimator to evaluate the fluctuation in numerical energy. Based on the a posteriori error estimations, the standard mesh refinement procedures were used in [18] for automatic mesh refining. In [17], the SCR based posteriori error estimators were constructed for space discretization of phase field variable and velocity function. The adaptive time-stepping method was applied to change the time step in [19,20].
In this work, we first provide an error estimate for the second-order fully discrete technique suggested in [1]. The fully discrete scheme is then introduced with an effective time-space adaptive technique called SCR. We need understand that the SCR approach has the same or higher accuracy while being more effective and affordable than other adaptive techniques. The SCR strategy was originally proposed by Huang and Yi [21], and they applied it to the Allen-Cahn system [16]. Moreover, for the definition of error estimator, the time discretization error estimator is given by the difference of numerical approximation between adjacent time steps, and the spatial discretization error estimator is defined by gradient, that is, the difference between the reconstructed gradient and the numerical gradient.
The remainder of the essay is structured as follows. For the Cahn-Hilliard equation, we offer a second-order fully discrete scheme in Section 2, after which we examine the error estimation and unconditional energy stability. In Section 3, we introduce the characteristics of the superconvergent cluster recovery operator, and define the error estimation operator in time and space. Several numerical tests illustrate the effectiveness of our adaptive procedure in Section 4.
First of all, we introduce a new variable μ=:−εΔu+1εf(u) and rewrite Eq (1.1) as follows.
{ut=Δμin Ω×(0,T],μ=−εΔu+1εf(u)in Ω×(0,T],∂nu=∂nμ=0on ∂Ω×(0,T],u(x,0)=u0in Ω, | (2.1) |
where ε>0 is a small parameter, which represents the thickness of the transition interface between materials, and u means the phase field variable, u0 ∈ H1(Ω) as the initial value at t=0. The solution is driven into the two pure states u=±1 by the nonlinear part f(u)=F′(u), and F(u) is a double-well potential function, which is given by F(u)=14(u2−1)2.
Actually, the Cahn-Hilliard equation could be viewed as the H−1−gradient flow of the free energy E(u):
E(u)=∫Ω(1εF(u)+ε2|∇u|2)dx, | (2.2) |
Furthermore, by a simple mathematical derivation, for the Cahn-Hilliard equation, the following energy law can be derived:
ddtE(u)=−∫Ω|∇(−εΔu+1εf(u))|2dx≤0, | (2.3) |
that the energy-decreasing property is satisfied by the free energy E(u), over time,
E(u)(tn+1)≤E(u)(tn),∀n∈N. | (2.4) |
Besides, the Cahn-Hilliard equation satisfies the law of conservation of mass [2].
ddt∫Ωu(x,t)dx=0,0≤t≤T. | (2.5) |
Next, we will present a few notations that will be used all across the paper before we get started. The standard Sobolev space will be indicated by the notation Hm=Wm,2(Ω), and its associated norm of Hm can be represented by ‖⋅‖m. As well as ‖⋅‖ and (⋅,⋅) respectively stand in for norm and the inner product of L2. Besides, we use the symbol ‖⋅‖∞ to signify the L∞ norm.
Let size h be the largest element diameter for every eh and Th={eh} be the set of regularly shaped triangles. Let Vh be the corresponding space in finite dimensions for piecewise linear continuous functions:
Vh={v∈H1(Ω):v|eh∈P1(eh),∀eh∈Th}, |
and its basis functions are the standard Lagrange basis functions ϕz(z∈Nh, where the Nh represents the set of the triangular vertices). Naturally, we denote V0h=Vh⋂H10(Ω).
In addition, for the real number p≥0 and satisfying the condition v(t), we define the following norm:
‖v‖Lp(0,T;X)=(∫T0‖v(t)‖pXdt)12. |
Then, by using the above symbol, we first introduce the weak form of Eq (2.1), which is as follow: find u∈L∞(0,T;H1(Ω)), w∈L2(0,T;H1(Ω))
{(ut,w)+(∇μ,∇w)=0∀w∈H1(Ω),(μ,v)=ε(∇u,∇v)+1ε(f(u),v)∀v∈H1(Ω),(u(0)−u0,q)=0∀q∈H1(Ω). | (2.6) |
To solve Eq (2.6), the semi-discrete system needs to find uh∈C1(0,T;Vh)
{(uh,t,wh)+(∇μh,∇wh)=0∀wh∈Vh,(μh,vh)=ε(∇uh,∇vh)+1ε(f(uh),vh)∀vh∈Vh,(uh(0)−u0h0,qh)=0∀qh∈Vh, | (2.7) |
where u0h is an approximation of u0 in Vh.
Divide [0,T] into Nst subintervals In:=(tn−1,tn], uniformly, 0=t0<t1<...<tN=T. Let △t represent the time step size so T=N△t, n=1,...,N. The fully discrete scheme for the nonlinear approximation, which is developed from a modified Crank-Nicolson approach, is as follows: find unh∈Vnh,n=1,2,...,N
{(unh−Πnun−1h△t,wh)+(∇μn−12h,∇wh)=0∀wh∈Vnh,(μn−12h,vh)=ε(∇^unh,∇vh)+1ε(f(unh,un−1h),vh)∀vh∈Vnh,(uh(0)−u0h0,qh)=0∀qh∈V0h, | (2.8) |
where ^unh, and f(unh,un−1h) are defined as
^unh=unh+Πnun−1h2,f(unh,un−1h)=(unh)3+(unh)2Πnun−1h+unh(Πnun−1h)2+(Πnun−1h)34−unh+Πnun−1h2, | (2.9) |
as well as μn−12 is the approximation at the midpoint tn−12=(tn−1+tn)/2 (directly computed), and Πn represents the interpolation into the finite element space Vh. Next, we consider the mass conservative, taking wh=1 in (2.8), the following results can be easily obtained
∫Ωun=∫Ωun−1=∫Ωu0. |
Obviously, the fully discrete scheme is mass conservative.
In addition, We use the convex partitioning method to divide f into two parts, i.e., the linear convex part and the nonlinear concave part. The linear convex part is then handled by the modified Crank-Nicolson method, while the linear convex part is handled by the Midpoint approximation (MP) method. So far, the modified Crank-Nicolson method has been applied several times to the Cahn-Hilliard equation to deal with the nonlinear concave terms, such as [22,23,24,25], and the linear convex terms in these works of literature are treated by BDF2.
Lemma 2.1. (Discrete Gronwall lemma) [26]. Let C and △t be non-negative constants, and ak,bk,ck,dk be non-negative sequences satisfying
ak+△tn∑k=0bk≤△tn−1∑k=0dkak+△tn−1∑k=0ck+C0,∀n≥1, |
then
an+△tn∑k=0bk≤exp(△tn−1∑k=0dk)(△tn−1∑k=0ck+C0),∀n≥1. |
Theorem 2.1. The solution of the fully discrete scheme Eqs (2.8) and (2.9), which is unconditionally energy stable, satisfies
E(unh)≤E(un−1h),∀n≥1. |
Proof. Taking wh=△tμn−12h and vh=unh−un−1h in Eq (2.8) we obtain
{(unh−un−1h,μn−12h)+△t‖∇μn−12h‖2=0,(μn−12h,unh−un−1h)=ε2(‖∇unh‖2−‖∇un−1h‖2)+1ε(f(unh,un−1h),unh−un−1h), |
where
1ε(f(unh,un−1h),unh−un−1h)=1ε((unh)3+(unh)2un−1h+unh(un−1h)2+(un−1h)34−unh+un−1h2,unh−un−1h)=1ε∫14(unh)4−14(un−1h)4−12(unh)2+12(un−1h)2dx, |
by integrating the upper formula, we obtain
0=△t‖∇μn−12h‖2+ε2(‖∇unh‖2−‖∇un−1h‖2)+1ε∫14(unh)4−14(un−1h)4−12(unh)2+12(un−1h)2dx. |
Also for energy we have
E(unh)−E(un−1h)=∫1εF(unh)+ε2|∇unh|2−1εF(un−1h)+ε2|∇un−1h|2dx=1ε∫14(unh)4−14(un−1h)4−12(unh)2+12(un−1h)2dx+ε2(‖∇unh‖2−‖un−1h‖2), |
the equation is finally obtained
△t‖∇μn−12h‖2+E(unh)−E(un−1h)=0, |
which means that the energy-decreasing feature is maintained by the provided fully discrete scheme Eqs (2.8) and (2.9). We omit Πn since it implies interpolation into the finite element space Vnh and has no bearing on the theoretical derivation procedure.
To derive the error estimates, we define the elliptic projection operator as Ph:H1→Vh, which satisfies
(∇(Phv−v),∇αh)=0,∀αh∈Vh. | (2.10) |
From the literature [27], we know the following two inequalities hold,
‖v−Phv‖≤Ch2‖v‖2,‖(v−Phv)t‖≤Ch2‖vt‖2. | (2.11) |
Then, we denote the error functions as,
˜en=Phu(tn)−unh,ˆen=u(tn)−Phu(tn),ˉen−12=Phμ(tn−12)−μn−12h,ˇen−12=μ(tn−12)−Phμ(tn−12). |
We may quickly derive by using the Taylor expansion with integral residuals and Young's inequality [12].
‖Rn1‖2s≤△t3∫tntn−1‖uttt(t)‖2sdt,s=−1,0, | (2.12) |
‖Rn2‖2s≤△t3∫tntn−1‖utt(t)‖2s+2dt,s=−1,0, | (2.13) |
where
Rn1=u(tn)−u(tn−1)△t−ut(tn−12), |
Rn2=−△(u(tn)+u(tn−1)2−u(tn−12)). |
Theorem 2.2. Let unh and u(tn) represent, respectively, the solutions of Eqs (2.8), (2.9) and (2.1). After that, for u∈C(0,T;H2(Ω)), ut∈L2(0,T;H2(Ω))⋂L2(0,T;L4(Ω)), utt∈L2(0,T;H1(Ω)), and uttt∈L2(0,T;H−1(Ω)) we get
‖ukh−u(tk)‖+(△tεk∑n=0‖μn−12h−μ(tn−12)‖2)12⩽C(ε,T)(K1(ε,u)△t2+K2(ε,u)h2), |
where
C(ε,T)∼exp(T∖ε),K1(ε,u)=√ε(‖uttt‖L2(0,T;H−1)+‖utt‖L2(0,T;H1))+1√ε(‖utt‖L2(0,T;L2)+‖ut‖L2(0,T;L4)),K2(ε,u)=‖u0‖2+√ε‖ut‖L2(0,T;H2)+1√ε‖μ‖C(0,T;H2)+‖u‖C(0,T;H2). |
Proof. Subtracting Eq (2.8) from the weak form of Eq (2.6) on Ω at tn−12, we obtain
1△t(˜en−˜en−1,wh)+(∇ˉen−12,∇wh)=(Rn1−1△t(I−Ph)(u(tn)−u(tn−1)),wh),(ˉen−12+ˇen−12,vh)=ε2(∇˜en+∇˜en−1,vh)+ε(△Rn2,vh)+1ε(f(u(tn−12))−f(unh,un−1h),vh), |
then taking wh=12ε△t(˜en+˜en−1) and vh=△tˉen−12 respectively, and summing up the two identities above, we gain
‖˜en‖2−‖˜en−1‖2+2△tε‖ˉen−12‖2=△t(Rn1,˜en+˜en−1)−(I−Ph)(u(tn)−u(tn−1),˜en+˜en−1)−2△t(Rn2,ˉen−12)+2△tε2(f(u(tn−12))−f(unh,unn−1),ˉen−12)−2△tε(ˇen−12,ˉen−12):=I+II+III+IV+V, |
where
I:=△t(Rn1,˜en+˜en−1),II:=−(I−Ph)(u(tn)−u(tn−1),˜en+˜en−1),III:=−2△t(Rn2,ˉen−12),IV:=2△tε2(f(u(tn−12))−f(unh,un−1h),ˉen−12),V:=−2△tε(ˇen−12,ˉen−12). |
the Cauchy inequality with ε and Young's inequality were used to estimate I, II,III, and V as following:
I≤△t‖Rn1‖‖˜en+˜en−1‖≤△t2(ε‖Rn1‖2+1ε‖˜en+˜en−1‖2)≤ε△t42∫tntn−1‖uttt(t)‖2−1dt+△t2ε‖˜en‖2+△t2ε‖˜en−1‖2,II≤‖(I−Ph)(u(tn)−u(tn−1))‖‖˜en+˜en−1‖≤ε2∫tntn−1‖(I−Ph)ut(t)‖2dt+12ε‖˜en‖2+12ε‖˜en−1‖2,III≤3△tε‖Rn2‖2+△t3ε‖ˉen−12‖2≤3ε△t4∫tntn−1‖utt(t)‖21dt+△t3ε‖ˉen−12‖2,V≤2△tε(√3‖ˇen−12‖)(1√3‖ˉen−12‖)≤3△tε‖ˇen−12‖2+△t3ε‖ˉen−12‖2. |
For simplicity, we denote u(tn)=un and u(tn−12)=un−12. We consider the fourth term Ⅳ as follows:
IV=2△tε2(f(un−12)−f(unh−un−1h),ˉen−12)=2△tε2(unh−un−1h2−un−12+(un−12)3−(unh)3+(unh)2un−1hh+unh(un−1h)2+(un−1h)34)=2△tε2(unh−un−1h2−un−12+(un−12)3−(un−un−12)3+(un−un−12)3−gn+gn−(unh)3+(unh)2un−1h+unh(un−1h)2+(un−1h)34,ˉen−12):=2△tε2(IV1+IV2+IV3+IV4,ˉen−12), |
where
gn=(un)3+(un)2un−1+un(un−1)2+(un−1)34,IV1:=unh+un−1h2,IV2:=(un−12)3−(un+un−12)3,IV3:=(un+un−12)3−gn,IV4:=gn−(unh)3+(unh)2un−1h+unh(un−1h)2+(un−1h)34. |
We continue our analysis of these four terms using the Taylor expansion with integral remainders, Cauchy-Schwarz inequality, and Eqs (2.12) and (2.13), we arrive at:
‖IV1‖=‖−˜en−ˆen−˜en−1−ˆen−1+un+un−12−un−12‖≤‖˜en+˜en−12‖+‖ˆen+ˆen−12‖+‖un+un−12−un−12‖≤‖˜en+˜en−12‖+‖ˆen+ˆen−12‖+△t3∫tntn−1‖utt(t)‖2dt, |
‖IV2‖=‖(un−12)3−(un+un−12)3‖≤‖3ξ21(un+un−12−un−12)‖≤3ξ21△t3∫tntn−1‖utt(t)‖2dt, |
‖IV3‖=‖(un)3−(un)2un−1−un(un−1)2+(un−1)38‖=‖((un)2−(un−1)2)(un−un−1)8‖≤‖ξ24(un−un−1)2‖≤ξ24△t3∫tntn−1‖u2t(t)‖2dt, |
where ξ1 lies between (un+un−1)/2 and un−12, ξ2 lies between un and un−1, ζ1 and ζ2 are located between un−1 and un.
‖IV4‖=14‖23(un)3+13(un+un−1)3+23(un−1)3−23(unh)3+13(unh+un−1h)3+23(un−1h)3‖=‖16((un)3−(unh)3)+16((un−1)3−(un−1h)3)+112((un+un−1)3−(unh+un−1h)3)‖≤12‖ξ23(un−unh)‖+12‖ξ24(un−1−un−1h)‖+14‖ξ25(un−unh+un−1−un−1h)‖≤C(‖un−1−un−1h‖+‖un−unh‖)≤C(‖ˆen−1+˜en−1‖+‖ˆen+˜en‖)≤C(‖ˆen−1‖+‖˜en−1‖+‖ˆen‖+‖˜en‖), |
where ξ3 lies between un and un−1, ξ4 lies between un−1 and un−1h, and ξ5 lies between un+un−1 and unh+un−1h. To guarantee the existence of ξ3,ξ4,ξ5, we need the condition [28] (Suppose C, a positive constant that is unaffected by △t or h, exists).
‖unh‖∞≤C,∀1≤n≤N. |
Further, combining the above inequalities IVi,i=1,2,3,4 into IV, we arrive at
IV=2△tε2(IV1+IV2+IV3+IV4,ˉen−12)≤2△tε2‖IV1+IV2+IV3+IV4‖‖ˉen−12‖≤2△tε2(‖IV1‖+‖IV2‖+‖IV3‖+‖IV4‖)‖ˉen−12‖≤2△tε2(‖˜en+˜en−12‖+‖ˆen+ˆen−12‖+△t3∫tntn−1‖utt(t)‖2dt+3ξ21△t3∫tntn−1‖utt(t)‖2dt+ξ24△t3∫tntn−1‖u2t(t)‖2dt+C(‖ˆen−1‖+‖˜en−1‖+‖ˆen‖+‖˜en‖))‖ˉen−12‖≤3△tε(‖˜en+˜en−12‖+‖ˆen+ˆen−12‖+△t3∫tntn−1‖utt(t)‖2dt+3ξ21△t3∫tntn−1‖utt(t)‖2dt+ξ24△t3∫tntn−1‖u2t(t)‖2dt+C(‖ˆen−1‖+‖˜en−1‖+‖ˆen‖+‖˜en‖))2+△t3ε‖ˉen−12‖2≤C△t4ε∫tntn−1‖utt(t)‖2dt+C△t4ε∫tntn−1‖u2t(t)‖2dt+C(‖ˆen−1‖2+‖˜en−1‖2+‖ˆen‖2+‖˜en‖2)+△t3ε‖ˉen−12‖2. |
Moreover, when we combine the terms I,II,III,IV and V, we arrive at the following estimates,
‖˜en‖2−‖˜en−1‖2+2△tε‖ˉen−12‖2≤ε△t42∫tntn−1‖uttt(t)‖2−1dt+ε2∫tntn−1‖(I−Ph)utt(t)‖2dt+3△tε‖ˇen−12‖2+△t3ε‖ˉen−12‖2+C△t4ε∫tntn−1‖utt(t)‖2dt+C△t4ε∫tntn−1‖u2t(t)‖2dt+3ε△t4∫tntn−1‖utt(t)‖21dt+△t3ε‖ˉen−12‖2+C(‖ˆen−1‖2+‖˜en−1‖2+‖ˆen‖2+‖˜en‖2)+△t3ε‖ˉen−12‖2. |
Using Eq (2.11), we may add up the aforementioned inequality for n=1,...,k(k≤T/△t), and we find,
‖˜ek‖2−‖˜e0‖2+△tεk∑n=1‖ˉen−12‖2≤12△t4ε‖uttt(t)‖2L2(0,T;H−1)+ε2Ch4‖ut(t)‖2L2(0,T;H2)+3△tεk∑n=1‖ˇen−12‖2+C△t4ε‖utt(t)‖2L2(0,T;L2)+C△t4ε‖ut(t)‖L2(0,T;L4)+3ε△t4‖utt(t)‖2L2(0,T;H1)+Ck∑n=1(‖ˆen−1‖2+‖˜en−1‖2+‖ˆen‖2+‖˜en‖2). |
Lastly, by using the triangular inequality and the discrete Gronwall lemma to the above inequality, we may conclude the proof.
The SCR technique takes derivatives to determine the recovered gradient at recovered locations after fitting a linear polynomial to solution values in a set of suitable sampling points around the vertex. In comparison to other methods, the SCR recovery technique implements adaptive algorithms more simply while saving on computing costs. The crucial component of the SCR technique is the introduction of a posterior error estimate operator. The process for getting the recovered gradient, which is broken down into three parts, is described in the sections that follow.
Step1: We plan to restore the gradient of an interior vertex a=a0∈Nh in Ω (Nh represents the mesh nodes), try to choose some points symmetrically distributed around z, and then form a new point set A={ai=(xi,yi)}, 1≤i≤n (n≥4). Specially, if select mesh nodes, we will get better results.
Step2: The recovery gradient operator Gh:Vh→Vh×Vh is defined as follows
(Ghuh)(a)=∇pa(x,y). |
Let Kz represent the sample points of convex polygons, the linear polynomial pa(x,y) can be found
pa(x,y):=argminp1∈P1n∑i=0(p1(ai)−uh(ai))2, |
In order to overcome the instability caused by small parameter h, we denote F by
F:(x,y)→(ψ,φ)=(x,y)−(x0,y0)h, |
where h:=max{|xi−x0|,|yi−y0|}, i=1,2,...,n. Then the fitting polynomial can be written as
pa(x,y)=PTm=ˆPTˆm, |
with
PT=(1,x,y),ˆPT=(2,ψ,φ),mT=(m1,m2,m3),ˆmT=(ˆm1,ˆm2,ˆm3). |
This ˆm is the coefficient vector, it satisfies the linear systems
ATAˆm=ATu, |
where
A=(1ψ0φ01ψ1φ1⋮⋮⋮1ψnφn)andu=(u(a0)u(a1)⋮u(an)). |
Finally, the recovered gradient can be derived
Ghu=∇pa=(m2m3)=1h(ˆm2ˆm3) |
Step3: the recovered gradient Ghu on Ω is obtained by the interpolation
Ghu=∑a∈NhGhu(a)ϕa. |
Next, the following are the properties of the operator Gh was proved in [21].
a. For eh∈Th, there is a constant C which is independent of the value of h, satisfied
‖Ghv‖0,eh≤C|v|1,K∀v∈Vh, |
where K=⋃3i=1Ki with Ki indicates the components containing the sample points coming from the ith vertex of eh.
b. The recovering point a=a0=(x0,y0) is the center of the circle, and the sampling points ai=(xi,yi),i=1,2,...,n(n≥4) equally distributed around it. we derive
|∇v(a)−Ghv(a)|≤Ch2,∀v∈W3∞(eh). |
We use ηnt and ηh,eh as discretization error indicators in time and space respectively, and then adjust the time step and mesh adaption in the algorithm.
ηnt:=‖unh−Πnun−1h‖,ηh,eh:=‖Ghuh−∇uh‖0,eh,η2h=∑eh∈Thη2h,eh. |
During SCR discretization, the time step changes with the error estimate, and the grid size changes according to the time step, so each tn constructs a new grid, and continues to construct a new space Vh on that grid and Πn represents the interpolation into the finite element space Vh.
Remark 3.1. If the gradient is on the boundary ∂Ω, we can treat it in some methods such as a higher order extension, extrapolation, and take the average value.
Remark 3.2. The sample points we selected are placed with maximum symmetry around a. On the basis of this, we talk about the sample points in the next two situations. One is that the approximation order can be raised if they are also nodes. The other is that it will only improve recovery accuracy if they are not nodes.
Remark 3.3. The SCR strategy can be used not only for the Crank-Nicolson approximation of time-discrete the Cahn-Hilliard equations, but also for other time-discrete formats of the Cahn-Hilliard equations, for example, we can apply to the BDF approximation mentioned in the literature [29][30], and its stability and unique solvability have been demonstrated. Unlike the format in this paper, the format has three levels of time, so we need unh and un−1h to approximate un+1h.
{(3un+1h−4Πnunh+Πnun−1h2△t,vh)+(∇wn+1h,∇vh)=0,∀vh∈Vh,(wn+1h,φh)=ε2(∇un+1h,∇φh)+((un+1h)3−2Πnunh+Πnun−1h,φh)+A△t(∇(un+1h−Πnunh),∇φh),∀φh∈Vh, |
To demonstrate the viability of our strategy, some numerical findings are offered in this section. We discuss a few characteristics of the Cahn-Hilliard Eqs (2.8) and (2.9), which are fully discrete and have various beginning value requirements.
Example 4.1. In the first test, we use the initial condition listed below to investigate Eq (2.1).
u0(x,y)=tanh(((x−0.3)2+y2−0.22)/ε)tanh(((x+0.3)2+y2−0.22)/ε)×tanh((x2+(y−0.3)2−0.22)/ε)tanh((x2+(y+0.3)2−0.22)/ε), | (4.1) |
where the parameter Ω=[−1,1]×[−1,1],andε=0.01.
Sequences of the mesh adaption results and the numerical solutions for six distinct time steps are shown in Figure 1. We can observe that when the node count decreases, the mesh refinement moves in accordance with the level-set zeros. Figure 2 shows the energy level dropping over time. Furthermore, it is also possible to perceive intuitively how the number of nodes and time steps change over time. We see that during the stage of energy's rapid evolution, a small time-step is chosen to capture the change in the numerical solution, and as the system settles into a stable state, the time steps increase.
Example 4.2. In terms of the Cahn-Hilliard Eq (2.1). Allowing for the following initial condition, let Ω=[−1,1]×[−1,1].
u0(x,y)=tanh((x2+y2−0.152)/ε)×tanh(((x−0.31)2+y2−0.152)/ε)×tanh(((x+0.31)2+y2−0.152)/ε)×tanh((x2+(y−0.31)2−0.152)/ε)×tanh((x2+(y+0.31)2−0.152)/ε)×tanh(((x−0.31)2+(y−0.31)2−0.152)/ε)×tanh(((x−0.31)2+(y+0.31)2−0.152)/ε)×tanh(((x+0.31)2+(y−0.31)2−0.152)/ε)×tanh(((x+0.31)2+(y+0.31)2−0.152)/ε), | (4.2) |
with the parameter ε=0.01.
In this illustration, Figure 3 shows a succession of mesh adaptation results and their numerical solutions. It is therefore evident via six different time steps that the mesh transformation adapts the rule of the zeros level-set and the number of nodes is decreasing. We infer from Figure 4 that the energy that decays over time, the number of nodes, and the number of time steps all change with time. This numerical result therefore shows that the time steps will grow greater than before when the system stabilizes.
Example 4.3. In the third example, we apply the following initial value to the Cahn-Hilliard Eq (2.1)
u0(x,y)=tanh(((x−0.3)2+y2−0.252)/ε)tanh(((x+0.3)2+y2−0.32)/ε), | (4.3) |
where the parameter Ω=[−1,1]×[−1,1],andε=0.01.
Figure 5 shows the various adaptive meshes and the associated numerical solutions. It is evident that the mesh transformation applies the level-set rule for zeros. The discrete-time history in Figure 6 also shows that the energy degrades over time. Additionally, Figure 6 shows the connection between time and the quantity of nodes (time-step). In the early stages of the simulation, small steps are used to record changes in the numerical solution, and when the system stabilizes in the latter stages, bigger time steps are used.
Example 4.4. In the last example, Eq (2.1) with the following initial value was analyzed
u0(x,y)=0.01rand(x,y), | (4.4) |
where the parameter Ω=[−1,1]×[−1,1],andε=0.01.
In Figure 7, we plot the different adaptive meshes and their contour plots of the numerical solutions. Figure 8, it displays energy is getting more and more small over time and the distribution of the time steps as well as nodes. We come to the conclusion quite similar to the previous examples.
Remark 4.1. When ε takes a value less than 0.01, the change of node number throughout the process still follows the law of zero level set, and the energy is still declining. But the time step will become smaller and the computational cost of the program will become larger, so finally, in the paper we choose ε=0.01 as the experimental result.
Table 1 displays the whole discrete format's spatial and temporal precision Eqs (2.8) and (2.9). We take the initial value u0(x,y)=0.01rand(x,y), the parameter ε=0.01, and T=1. Since the precise solution to the initial equation is unknown, the numerical error ehu:=uh(x,T)−uh2(x,T) will be determined by the difference between the coarse and fine grids, and the convergence order log2(‖ehu‖/‖eh2u‖) will be determined by the ratio of the errors. The table shows that the accuracy converges to 2, which is consistent with our theoretical findings for both spatial and temporal accuracy.
h | ‖umh−umh2‖0 | Rate | △t | Rate |
116 | 0.052470400 | 116 | ||
132 | 0.018185700 | 1.528700 | 132 | 2.18632 |
164 | 0.003449260 | 2.398450 | 164 | 2.52708 |
1128 | 0.000945365 | 1.867340 | 1128 | 1.78946 |
The fundamental idea of this paper is to obtain the spatial discretization operator by SCR method, and then use it as the index of spatial discretization of the Cahn-Hilliard equation. Additionally, the index of temporal discretization is determined by the variation in approximate solutions between adjacent time steps. First of all, in order to solve the Cahn-Hilliard equation, we first develop a second-order scheme that is unconditionally energy stable. Then, the fully discrete system has second-order accuracy both time and space, as determined by error estimation. Finally, these numerical findings are used to demonstrate the efficacy of this approach.
Tian's research was supported by Shanxi Scholarship Council of China (No. 2021029) and the 2021 Shanxi Science and Technology Cooperation and Exchange Special Program (No. 202104041101019). Chen's reserach was partially supported by NSFC Project (12201010).
No potential conflict of interest was reported by the authors.
[1] |
Ramírez-Velazco A, Aguayo-Orozco TA, Figuera L, et al. (2019) Williams–Beuren syndrome in Mexican patients confirmed by FISH and assessed by aCGH. J Genet 98: 34. https://doi.org/10.1007/s12041-019-1080-7 ![]() |
[2] |
Yuan SM (2017) Congenital heart defects in Williams syndrome. Turk J Pediatr 59: 225-231. https://doi.org/10.24953/turkjped.2017.03.001 ![]() |
[3] |
Debelle L, Tamburro AM (1999) Elastin: Molecular description and function. Int J Biochem Cell Biol 31: 261-272. https://doi.org/10.1016/S1357-2725(98)00098-3 ![]() |
[4] |
Ewart AK, Morris CA, Atkinson D, et al. (1993) Hemizygosity at the elastin locus in a developmental disorder, Williams syndrome. Nat Genet 5: 11-16. https://doi.org/10.1038/ng0993-11 ![]() |
[5] |
van der Bom T, Zomer AC, Zwinderman AH, et al. (2011) The changing epidemiology of congenital heart disease. Nat Rev Cardiol 8: 50-60. https://doi.org/10.1038/nrcardio.2010.166 ![]() |
[6] |
Schopler E, Reichler RJ, DeVellis RF, et al. (1980) Toward objective classification of childhood autism: Childhood Autism Rating Scale (CARS). J Autism Dev Disord 10: 91-103. https://doi.org/10.1007/BF02408436 ![]() |
[7] |
Edelmann L, Prosnitz A, Pardo S, et al. (2006) An atypical deletion of the Williams-Beuren syndrome interval implicates genes associated with defective visuospatial processing and autism. J Med Genet 44: 136-143. https://doi.org/10.1136/jmg.2006.044537 ![]() |
[8] |
Depienne C, Heron D, Betancur C, et al. (2007) Autism, language delay and mental retardation in a patient with 7q11 duplication. J Med Genet 44: 452-458. https://doi.org/10.1136/jmg.2006.047092 ![]() |
[9] |
Berg JS, Brunetti-Pierri N, Peters SU, et al. (2007) Speech delay and autism spectrum behaviors are frequently associated with duplication of the 7q11.23 Williams-Beuren syndrome region. Genet Med 9: 427-441. https://doi.org/10.1097/GIM.0b013e3180986192 ![]() |
[10] |
Fusco C, Micale L, Augello B, et al. (2014) Smaller and larger deletions of the Williams Beuren syndrome region implicate genes involved in mild facial phenotype, epilepsy and autistic traits. Eur J Hum Genet 22: 64-70. https://doi.org/10.1038/ejhg.2013.101 ![]() |
[11] |
Alesi V, Loddo S, Orlando V, et al. (2021) Atypical 7q11.23 deletions excluding ELN gene result in Williams–Beuren syndrome craniofacial features and neurocognitive profile. Am J Med Genet A 185: 242-249. https://doi.org/10.1002/ajmg.a.61937 ![]() |
[12] |
Kruszka P, Porras AR, de Souza DH, et al. (2018) Williams–Beuren syndrome in diverse populations. Am J Med Genet A 176: 1128-1136. https://doi.org/10.1002/ajmg.a.38672 ![]() |
[13] |
Lee CL, Lin SM, Chen MR, et al. (2022) Long-term cardiovascular findings in Williams syndrome: A single medical center experience in Taiwan. J Pers Med 12: 817. https://doi.org/10.3390/jpm12050817 ![]() |
[14] |
Niego A, Benítez-Burraco A (2022) Autism and Williams syndrome: truly mirror conditions in the socio-cognitive domain?. Int J Dev Disabil 68: 399-415. https://doi.org/10.1080/20473869.2020.1817717 ![]() |
[15] |
Richards C, Jones C, Groves L, et al. (2015) Prevalence of autism spectrum disorder phenomenology in genetic disorders: a systematic review and meta-analysis. Lancet Psychiat 2: 909-916. https://doi.org/10.1016/S2215-0366(15)00376-4 ![]() |
[16] |
Masson J, Demily C, Chatron N, et al. (2019) Molecular investigation, using chromosomal microarray and whole exome sequencing, of six patients affected by Williams Beuren syndrome and Autism Spectrum Disorder. Orphanet J Rare Dis 14: 121. https://doi.org/10.1186/s13023-019-1094-5 ![]() |
[17] |
Ramocki MB, Bartnik M, Szafranski P, et al. (2010) Recurrent distal 7q11.23 deletion including HIP1 and YWHAG identified in patients with intellectual disabilities, epilepsy, and neurobehavioral problems. Am J Hum Genet 87: 857-865. https://doi.org/10.1016/j.ajhg.2010.10.019 ![]() |
[18] |
Sanders SJ, Ercan-Sencicek AG, Hus V, et al. (2011) Multiple recurrent de novo CNVs, including duplications of the 7q11.23 Williams syndrome region, are strongly associated with autism. Neuron 70: 863-885. https://doi.org/10.1016/j.neuron.2011.05.002 ![]() |
h | ‖umh−umh2‖0 | Rate | △t | Rate |
116 | 0.052470400 | 116 | ||
132 | 0.018185700 | 1.528700 | 132 | 2.18632 |
164 | 0.003449260 | 2.398450 | 164 | 2.52708 |
1128 | 0.000945365 | 1.867340 | 1128 | 1.78946 |