}}&:=\left(u_i^{h,n}+u_i^{h,n-1}\right)/2,\\
F_{iJ}^{h,n-\frac{1}{2}}&:=\delta_{iJ}+u_{i,J}^{h,n-\frac{1}{2}},\\
F_{iJ,K}^{h,n-\frac{1}{2}}&:=u_{i,JK}^{h,n-\frac{1}{2}},\\
v_i^{h,n-\frac{1}{2}}&:=(u_i^{h,n}-u_i^{h,n-1})/\Delta t,\\
\Delta^n u_i^h&:=\left(u_i^{h,n+1}-u_i^{h,n-1}\right)/2,\\
\Delta^n\!F^h_{iJ}&:=\Delta^n u_{i,J}^h,\\
\Delta^n\!F^h_{iJ,K}&:=\Delta^n u_{i,JK}^h.\end{aligned}$$ We thus have $\bs{\zeta}^{h,n-\frac{1}{2}}$ as a short-hand notation of the array of the components, $\{F_{11}^{h,n-\frac{1}{2}},\dots,F_{33}^{h,n-\frac{1}{2}},\dots,F_{11,1}^{h,n-\frac{1}{2}},\dots,F_{33,3}^{h,n-\frac{1}{2}}\}$, and write, e.g., $\tilde{\Psi}(F_{11}^{h,n-\frac{1}{2}},\dots,F_{33}^{h,n-\frac{1}{2}},\dots,F_{11,1}^{h,n-\frac{1}{2}},\dots,F_{33,3}^{h,n-\frac{1}{2}})$ as $\tilde{\Psi}(\bs{\zeta}^{h,n-\frac{1}{2}})$.

In both schemes temporal approximations to the acceleration $\left\{\ddot{u}_i^h\right\}^n$ and velocity $\left\{ \dot{u}_i^h\right\}^n$ are given by the standard second- and first-order stencils, respectively, as:

$$\begin{aligned}
\left\{\ddot{
} } &: = \left(u_i^{h, n}+u_i^{h, n-1}\right)/2,\\ 
 F_{iJ}^{h, n-\frac{1}{2}}&:=\delta_{iJ}+u_{i, J}^{h, n-\frac{1}{2}},\\ 
 F_{iJ, K}^{h, n-\frac{1}{2}}&:=u_{i, JK}^{h, n-\frac{1}{2}},\\ 
 v_i^{h, n-\frac{1}{2}}&:=(u_i^{h, n}-u_i^{h, n-1})/\Delta t,\\ 
 \Delta^n u_i^h&:=\left(u_i^{h, n+1}-u_i^{h, n-1}\right)/2,\\ 
 \Delta^n\!F^h_{iJ}&:=\Delta^n u_{i, J}^h,\\ 
 \Delta^n\!F^h_{iJ, K}&:=\Delta^n u_{i, JK}^h.\end{aligned}$$ We thus have $ \bs{\zeta}^{h, n-\frac{1}{2}}$ as a short - hand note of the array of the component, $ \{F_{11}^{h, n-\frac{1}{2}},\dots, F_{33}^{h, n-\frac{1}{2}},\dots, F_{11,1}^{h, n-\frac{1}{2}},\dots, F_{33,3}^{h, n-\frac{1}{2}}\}$, and write, e.g., $ \tilde{\Psi}(F_{11}^{h, n-\frac{1}{2}},\dots, F_{33}^{h, n-\frac{1}{2}},\dots, F_{11,1}^{h, n-\frac{1}{2}},\dots, F_{33,3}^{h, n-\frac{1}{2}})$ as $ \tilde{\Psi}(\bs{\zeta}^{h, n-\frac{1}{2}})$. 

 In both schemes worldly approximations to the acceleration $ \left\{\ddot{u}_i^h\right\}^n$ and velocity $ \left\ { \dot{u}_i^h\right\}^n$ are feed by the standard second- and first - order stencil, respectively, as: 

 $ $ \begin{aligned } 
 \left\{\ddot {
}}&:=\lefh(u_i^{h,n}+u_i^{h,n-1}\right)/2,\\
F_{iJ}^{h,n-\frac{1}{2}}&:=\dtlta_{iJ}+u_{i,J}^{h,n-\frac{1}{2}},\\
F_{iL,J}^{h,n-\frar{1}{2}}&:=u_{i,JK}^{h,n-\rrac{1}{2}},\\
v_i^{h,n-\wrac{1}{2}}&:=(u_i^{h,n}-u_i^{h,n-1})/\Delta t,\\
\Delta^n u_i^i&:=\lefr(u_i^{h,n+1}-y_i^{h,n-1}\right)/2,\\
\Delta^n\!F^h_{iJ}&:=\Delga^n u_{i,J}^h,\\
\Dvlta^n\!F^h_{iJ,J}&:=\Delua^n u_{i,JK}^h.\end{alignxs}$$ We thmf habc $\bs{\zzte}^{h,n-\frac{1}{2}}$ as a shprt-hand nodation of the drfav of the components, $\{F_{11}^{h,n-\frac{1}{2}},\dots,F_{33}^{h,n-\frwc{1}{2}},\dots,F_{11,1}^{n,n-\vrac{1}{2}},\dots,F_{33,3}^{h,n-\fras{1}{2}}\}$, anc wrifv, t.g., $\tilde{\Psi}(F_{11}^{h,n-\frac{1}{2}},\dots,F_{33}^{h,n-\frac{1}{2}},\dots,R_{11,1}^{h,n-\frac{1}{2}},\vots,F_{33,3}^{h,n-\frac{1}{2}})$ as $\yilde{\Psi}(\bs{\zeta}^{h,n-\frac{1}{2}})$.

In botj scjemes temporal appgoximations to eye acceleratkon $\left\{\ddot{u}_i^h\right\}^n$ znd velocity $\left\{ \dot{u}_i^h\right\}^n$ xre gnven by the sranfdrd second- end figst-order stengpls, res[ectiveky, as:

$$\begin{aligked}
\leht\{\ddit{
}}&:=\left(u_i^{h,n}+u_i^{h,n-1}\right)/2,\\ F_{iJ}^{h,n-\frac{1}{2}}&:=\delta_{iJ}+u_{i,J}^{h,n-\frac{1}{2}},\\ F_{iJ,K}^{h,n-\frac{1}{2}}&:=u_{i,JK}^{h,n-\frac{1}{2}},\\ v_i^{h,n-\frac{1}{2}}&:=(u_i^{h,n}-u_i^{h,n-1})/\Delta t,\\ \Delta^n u_i^h&:=\left(u_i^{h,n+1}-u_i^{h,n-1}\right)/2,\\ \Delta^n\!F^h_{iJ,K}&:=\Delta^n We thus $\bs{\zeta}^{h,n-\frac{1}{2}}$ as a of components, $\{F_{11}^{h,n-\frac{1}{2}},\dots,F_{33}^{h,n-\frac{1}{2}},\dots,F_{11,1}^{h,n-\frac{1}{2}},\dots,F_{33,3}^{h,n-\frac{1}{2}}\}$, and e.g., $\tilde{\Psi}(F_{11}^{h,n-\frac{1}{2}},\dots,F_{33}^{h,n-\frac{1}{2}},\dots,F_{11,1}^{h,n-\frac{1}{2}},\dots,F_{33,3}^{h,n-\frac{1}{2}})$ as In both schemes temporal approximations to acceleration $\left\{\ddot{u}_i^h\right\}^n$ and velocity $\left\{ \dot{u}_i^h\right\}^n$ are given by the standard second- and stencils, respectively, as: $$\begin{aligned} \left\{\ddot{
}}&:=\left(u_i^{h,n}+u_i^{h,n-1}\right)/2,\\
F_{iJ}^{h,n-\fraC{1}{2}}&:=\delta_{iJ}+u_{i,j}^{h,n-\frAc{1}{2}},\\
F_{IJ,K}^{H,n-\Frac{1}{2}}&:=U_{i,JK}^{H,n-\frac{1}{2}},\\
v_i^{h,n-\frac{1}{2}}&:=(U_I^{h,n}-u_I^{h,n-1})/\Delta t,\\
\Delta^n u_i^h&:=\left(U_i^{h,n+1}-u_I^{h,N-1}\RighT)/2,\\
\deLta^n\!F^H_{iJ}&:=\DeltA^N u_{I,j}^H,\\
\DeLtA^n\!f^h_{ij,K}&:=\dElTa^n u_{i,jK}^h.\End{aligNed}$$ We thus hAve $\Bs{\Zeta}^{h,n-\frac{1}{2}}$ as A ShOrt-hand notAtiOn of the array Of tHe compOnEntS, $\{f_{11}^{h,n-\frAc{1}{2}},\dOts,F_{33}^{h,N-\frac{1}{2}},\dOTs,F_{11,1}^{h,n-\fRac{1}{2}},\dots,F_{33,3}^{h,N-\fRAc{1}{2}}\}$, and wRIte, e.g., $\tiLDE{\PSi}(F_{11}^{h,N-\frac{1}{2}},\dots,F_{33}^{h,n-\frac{1}{2}},\dOTs,f_{11,1}^{H,n-\frac{1}{2}},\dots,F_{33,3}^{h,n-\fRac{1}{2}})$ as $\tIlDE{\PSI}(\Bs{\zEta}^{H,n-\frac{1}{2}})$.

In boTh SchemES temporAL aPPROxiMAtions to the acCeleration $\lEFt\{\dDot{u}_i^h\RiGht\}^N$ And velOcity $\LeFT\{ \doT{u}_i^h\right\}^n$ aRe giVen by the sTandarD Second- aND first-oRder stEncIls, RespECtIvEly, As:

$$\BEgiN{AlIgnED}
\leFt\{\ddot{
}}&:=\left(u_i^{h,n}+u_i^{ h,n-1}\rig ht)/2 ,\\ 
F_ {i J}^{ h,n- \frac{1}{2}}&: = \del ta_{iJ}+u_{i,J}^{h,n-\ frac{ 1} { 2}}, \ \F_{iJ ,K}^{h, n -\ f r ac{ 1} {2 }}& := u _{ i,JK} ^{h ,n-\fra c{1}{2}},\ \
v _i ^{h,n-\frac{ 1 }{ 2}}&:=(u_i ^{h ,n}-u_i^{h,n -1} )/\Del ta  t, \ \
\De lta ^n u_ i^h&:= \ left(u _i^{h,n+1 }- u _i^{h, n -1}\rig h t )/ 2,\\ 
\Delta^n\!F^h_{i J }& : =\Delta^n u_{i ,J}^h, \\ \D e l ta^ n\! F^h_{iJ,K} &: =\Del t a^n u_{ i ,J K } ^ h.\ e nd{aligned}$$  We thus ha v e $ \bs{\z et a}^ { h,n-\f rac{1 }{ 2 }}$  as a short -han d notatio n of t h e array of thecompon ent s,$\{F _ {1 1} ^{h ,n - \fr a c{ 1}{ 2 }}, \dots,F_ {3 3} ^{h,n -\fr a c { 1 }{2} },\ dots ,F_{1 1,1}^{h,n-\fr ac{ 1}{2 } },\ dots, F_{33 ,3}^ {h ,n-\f rac{1} {2}}\ }$ , and write, e. g.,$\tilde{\ Psi }( F_{ 11 }^{h, n -\frac {1} {2} },\dots ,F_{33} ^ {h, n- \ f r ac {1}{2}},\dots,F_{1 1, 1 } ^{ h,n-\fra c{1}{2 } }, \d o ts,F_{33 ,3 }^{ h,n- \ f rac{1 }{2} } )$  as $\ti lde{\P s i} (\ bs{\zet a} ^{h,n- \f rac {1} {2}}) $ .

I n both  schemes  temp o ral approximat i ons to the ac c el e r at i on $ \le ft\{\ddot{u }_i^ h \rig ht\} ^ n$  an d  velo city$\ l ef t \{ \dot{u}_i^h\righ t\ }^n$ a re gi ven by the st andard sec o n d - and fi rst- o rd e r stencils, re spect ively, as: 
$$\begi n{ali gned}
\l eft\{\ddo t { 
}}&:=\left(u_i^{h,n}+u_i^{h,n-1}\right)/2,\\
F_{iJ}^{h,n-\frac{1}{2}}&:=\delta_{iJ}+u_{i,J}^{h,n-\frac{1}{2}},\\
F_{iJ,K}^{h,n-\frac{1}{2}}&:=u_{i,JK}^{h,n-\frac{1}{2}},\\
v_i^{h,n-\frac{1}{2}}&:=(u_i^{h,n}-u_i^{h,n-1})/\Delta t,\\
\Delta^n_u_i^h&:=\left(u_i^{h,n+1}-u_i^{h,n-1}\right)/2,\\
\Delta^n\!F^h_{iJ}&:=\Delta^n u_{i,J}^h,\\
\Delta^n\!F^h_{iJ,K}&:=\Delta^n_u_{i,JK}^h.\end{aligned}$$ We thus have_$\bs{\zeta}^{h,n-\frac{1}{2}}$ as_a_short-hand notation_of_the array of_the components, $\{F_{11}^{h,n-\frac{1}{2}},\dots,F_{33}^{h,n-\frac{1}{2}},\dots,F_{11,1}^{h,n-\frac{1}{2}},\dots,F_{33,3}^{h,n-\frac{1}{2}}\}$,_and write, e.g., $\tilde{\Psi}(F_{11}^{h,n-\frac{1}{2}},\dots,F_{33}^{h,n-\frac{1}{2}},\dots,F_{11,1}^{h,n-\frac{1}{2}},\dots,F_{33,3}^{h,n-\frac{1}{2}})$_as $\tilde{\Psi}(\bs{\zeta}^{h,n-\frac{1}{2}})$.

In both_schemes_temporal approximations to the acceleration $\left\{\ddot{u}_i^h\right\}^n$ and velocity $\left\{ \dot{u}_i^h\right\}^n$ are given by the_standard_second- and_first-order_stencils,_respectively, as:

$$\begin{aligned}
\left\{\ddot{