Esquema CE/SE de Alta Ordem Para Simulação de Fluxo de Ruptura de Barragem

<!-- metadata commented in wiki content
==Esquema CE/SE de Alta Ordem Para Simulação de Fluxo de Ruptura de Barragem==
-->

==Abstract==

Este artigo apresenta um novo esquema explícito para a solução das equações de águas rasas em uma e duas  dimensões, desenvolvido a partir do método dos elementos de conservação e elementos de solução espaço-tempo,  aqui abreviado por método CE/SE. As funções de base utilizadas são expansões de Taylor de segunda  ordem no tempo e no espaço. Esse aumento na ordem das funções de aproximação produz o aumento  no número de variáveis de marcha no tempo, por isso, além das variáveis fluxo e suas inclinações, também são  incógnitas no presente esquema suas derivadas espaciais de segunda ordem. Um processo iterativo para o  cálculo das derivadas de primeira e segunda ordem é formulado para problemas com choques e descontinuidades.  Experimentos computacionais demonstram acurácia de terceira ordem. Os problemas de ruptura de  barragem unidimensional e bidimensional considerados validam a acurácia e robustez do esquema.

'''keywords'''

Leis de conservação,  Equações de águas rasas,  Volume de controle espaço-tempo,  Alta ordem de acurácia.

==High-Order CE/SE Scheme for Dam-Break Flow Simulation==

'''Abstract'''

This paper presents a new explicit scheme for the solution of shallow water equations in one and two  space dimensions, developed from the space-time conservation element and solution element (CE/SE) method. The  basis functions used are second-order Taylor expansions in time and space. This increase in the order of the  approximation functions produces an increase in the number of unknowns in the scheme, therefore, besides the  flow variables and their slopes, their second-order partial derivatives are also unknown in the present  scheme. An iterative process for the calculation of the first and second order derivatives is formulated for  problems with shocks and discontinuities. Computational experiments demonstrate third-order accuracy. The  one-dimensional and two-dimensional dam-break problems presented validate the accuracy and robustness of this  scheme.

''Keywords'': Conservation laws, Shallow water equations, Space-time  control volume, High-order accuracy.

==3.1 Introduco==

As equações de águas rasas constituem um dos modelos mais comumente usados na análise de fluxo de água em  rios ou áreas costais <span id='citeF-1'></span>[[#cite-1|[1]]], sobretudo no que tange à simulações de fluxo de ruptura de barragem  <span id='citeF-2'></span><span id='citeF-3'></span><span id='citeF-4'></span><span id='citeF-5'></span>[[#cite-2|[2,3,4,5]]], parte essencial do projeto e avaliação da  segurança de barragens, controle de inundação de rios, mitigação de desastres de bacias hidrográficas, etc   <span id='citeF-6'></span>[[#cite-6|[6]]]. Além da importância prática em hidráulica e engenharia costal, fornecem um ótimo modelo  matemático para equações diferenciais hiperbólicas não lineares que podem ter soluções como ondas de choque  <span id='citeF-7'></span>[[#cite-7|[7]]] e, por isso, também são frequentemente utilizadas como problema de referência para novos  esquemas de previsão numérica <span id='citeF-8'></span><span id='citeF-9'></span>[[#cite-8|[8,9]]].

Extensivos estudos numéricos foram feitos no sentido de analisar os fenômenos governados por essas equações e  diversos métodos foram desenvolvidos, como de volumes finitos <span id='citeF-10'></span><span id='citeF-11'></span><span id='citeF-12'></span>[[#cite-10|[10,11,12]]], elementos finitos <span id='citeF-13'></span><span id='citeF-14'></span><span id='citeF-15'></span><span id='citeF-16'></span>[[#cite-13|[13,14,15,16]]],  diferenças finitas <span id='citeF-17'></span><span id='citeF-18'></span>[[#cite-17|[17,18]]] e outros <span id='citeF-19'></span><span id='citeF-20'></span>[[#cite-19|[19,20]]].

Recentemente, um novo método baseado em volumes de controle, denominado método dos elementos de  conservação e elementos de solução espaço-tempo, do inglês, ''space-time conservation element  and solution element method (CE/SE)'', foi proposto por Chang e To (1991)  <span id='citeF-21'></span>[[#cite-21|[21]]] e aplicações do mesmo em problemas de ruptura de barragem logo foram feitas  <span id='citeF-22'></span><span id='citeF-23'></span><span id='citeF-24'></span>[[#cite-22|[22,23,24]]]. Este método foi desenvolvido com o objetivo de obter soluções  numéricas de leis de conservação e possui características importantes <span id='citeF-25'></span><span id='citeF-26'></span><span id='citeF-27'></span><span id='citeF-28'></span>[[#cite-25|[25,26,27,28]]]: sua construção combina informações de ambas as formas diferencial e integral das leis  de conservação; o conjunto de variáveis de marcha é formado pelas variáveis fluxo e suas derivadas espaciais;  não são utilizadas técnicas de interpolação ou extrapolação nos valores da malha; não são utilizadas técnicas  baseadas nas características, o que torna o método de fácil extensão à várias dimensões, bem como aplicável à  equações não hiperbólicas, como as equações de Navier-Stokes, por exemplo; possui uma simples notação  estêncil, o que facilita a programação. É caracterizado, também, por sua versatilidade, tendo em vista que  pode ser adaptado à malhas não uniformes e não estruturadas <span id='citeF-29'></span><span id='citeF-30'></span>[[#cite-29|[29,30]]], bem como pode ser  estendido à altas ordens <span id='citeF-31'></span><span id='citeF-32'></span>[[#cite-31|[31,32]]], razão pela qual possui aplicações nas mais diversas  áreas <span id='citeF-33'></span><span id='citeF-34'></span><span id='citeF-35'></span><span id='citeF-36'></span>[[#cite-33|[33,34,35,36]]].

O objetivo deste trabalho é apresentar a construção de um novo esquema CE/SE explícito de alta ordem para as  equações de águas rasas e verificar sua eficiência em problemas de ruptura de barragem. Para isso, este  artigo está organizado da seguinte maneira: na seção [[#3.2 Equações Governantes e Método CE/SE|3.2]] é apresentado o modelo matemático a ser  estudado (subseção [[#3.2.1 Equações de Águas Rasas|3.2.1]]), bem como é desenvolvido o esquema numérico discreto  (subseção [[#3.2.2 Esquema CE/SE|3.2.2]]); na seção [[#3.3 Exemplos Numéricos|3.3]], experimentos numéricos unidimensionais e bidimensionais  são considerados com a finalidade de validar o método; as conclusões são tecidas junto à seção [[#3.4 Conclusões|3.4]] e,  ao final, na seção [[#3.5 Referências|3.5]], dispõe-se a lista de referências.

==3.2 Equações Governantes e Método CE/SE==

===3.2.1 Equações de Águas Rasas===

As equações de águas rasas, também conhecidas como equações de Saint Venant, possuem a seguinte forma  vetorial conservativa

<span id="eq-3.1"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>\displaystyle \frac{\partial \textrm{'''Q'''}}{\partial t}+\frac{\partial  \textrm{'''F(Q)'''}}{\partial x}+\frac{\partial \textrm{'''G(Q)'''}}{\partial y} =  \textrm{'''S(Q)'''},  </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.1)
|}

em que

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>\textrm{'''Q'''} = \left[\begin{array}{c}q_1 \\ q_2 \\ q_3 \end{array}\right] = \left[\begin{array}{c}h \\ hu \\ hv \end{array}\right], \quad  \textrm{'''F'''} = \left[\begin{array}{c}f_1 \\ f_2 \\ f_3 \end{array}\right] =  \left[ \begin{array}{c}hu \\ hu^2+gh^2/2 \\ huv  \end{array} \right] =  \left[ \begin{array}{c}q_2 \\ q_2^2/q_1+gq_1^2/2 \\ q_2q_3/q_1  \end{array} \right], </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.2)
|}

<span id="eq-3.3"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>\textrm{'''G'''} = \left[\begin{array}{c}g_1 \\ g_2 \\ g_3 \end{array}\right] =  \left[ \begin{array}{c}hv \\ huv \\ hv^2+gh^2/2  \end{array} \right] =  \left[ \begin{array}{c}q_3 \\  q_2q_3/q_1\\ q_3^2/q_1+gq_1^2/2  \end{array} \right]\quad \textrm{e} \quad  \textrm{'''S'''} = \left[\begin{array}{c}S_1 \\ S_2 \\ S_3 \end{array}\right] =  \left[ \begin{array}{c}0 \\ gq_1(S_{0x}-S_{fx}) \\ gq_1(S_{0y}-S_{fy})  \end{array} \right],  </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.3)
|}

sendo <math display="inline">h</math> a profundidade do fluxo (m); <math display="inline">u</math> e <math display="inline">v</math> são as velocidades médias  do escoamento (m/s) nas direções <math display="inline">x</math> e <math display="inline">y</math>, respectivamente; <math display="inline">g</math> é a aceleração da gravidade (<math display="inline">9.812</math>m/s<math display="inline">^2</math>); <math display="inline">S_{0{x}} = \partial Z_b/\partial x</math>  e <math display="inline">S_{0{y}} = \partial Z_b/ \partial y</math> são as inclinações do fundo do  canal cuja topografia é dada por <math display="inline">Z_b</math>; <math display="inline">S_{f{x}}</math>  e <math display="inline">S_{f{y}}</math> são os termos de resistência ao escoamento,  nas direções <math display="inline">x</math> e <math display="inline">y</math>, respectivamente, definidos por

<span id="eq-3.4"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>S_{fx} = \displaystyle \frac{n^2u\sqrt{u^2+v^2}}{h^{4/3}} =  \frac{n^2q_2\sqrt{q_2^2+q_3^2}}{q_1^{10/3}}  \quad \textrm{e} \quad S_{fy} = \displaystyle \frac{n^2v\sqrt{u^2+v^2}}{h^{4/3}}  = \frac{n^2q_3\sqrt{q_2^2+q_3^2}}{q_1^{10/3}}, </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.4)
|}

em que <math display="inline">n</math> (s/m<math display="inline">^{1/3}</math>) é o coeficiente de fricção de Manning.

===3.2.2 Esquema CE/SE===

Considere, por simplicidade, a Eq. ([[#eq-3.1|3.1]]) como

<span id="eq-3.5"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>\displaystyle \frac{\partial q_m}{\partial{t}}+\frac{\partial f_m}{\partial x}+\frac{\partial g_m}{\partial  y} = S_m, \quad m = 1,2,3 </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.5)
|}

então, pelo teorema da divergência no espaço <math display="inline">\mathbb{R}^3</math>, tem-se que a Eq. ([[#eq-3.5|3.5]])  representa a forma diferencial da lei integral de conservação

<span id="eq-3.6"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>\oint _{S(V)}\textrm{'''H'''}_m\cdot d\textrm{'''s'''} = \int _V S_mdV, </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.6)
|}

em que '''H'''<math display="inline">_m = (f_m, g_m,q_m)</math>, <math display="inline">m = 1,2,3</math> e <math display="inline">S(V)</math> representa o contorno de uma região  espaço-tempo <math display="inline">V \subset \mathbb{R}^3</math>.

No método CE/SE, um elemento de conservação (CE) é uma região espaço-tempo em que a conservação do fluxo, Eq.  ([[#eq-3.6|3.6]]), é forçada, enquanto que um elemento de solução (SE) é uma região espaço-tempo, normalmente  distinta, em que as variáveis fluxo são supostamente suaves e a Eq. ([[#eq-3.5|3.5]]) é válida. Para definir estes  dois importante objetos, considere primeiramente uma malha no plano <math display="inline">x</math>-<math display="inline">y</math>, conforme Fig. ma. Existem  dois grupos de pontos, marcados por  círculos e cruzes, que representam nós da malha em dois níveis de tempo  diferentes. A cada ponto <math display="inline">(i,j,k)</math> da malha associa-se um CE e um SE. O CE é definido como o quadrilátero <math display="inline">EFGHE'F'G'H'</math> e o SE é a união do quadrilátero <math display="inline">P''Q''R''S''P'Q'R'S'</math> e o polígono <math display="inline">EFGH</math> (veja  Fig. [[#img-1a|1a]]).

<div id='img-1a'></div>
<div id='img-1'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:draft_García_699893715-CE_2D_2.png|240px|]]
|[[Image:draft_García_699893715-CE_2D.png|300px|a) Pontos da malha representativa no plano x-y e (b) as definições dos  elementos de solução SE(i,j,k) e elementos de conservação CE(i,j,k) no ponto (i,j,k) da malha  [2012_Zhang, 1999_Zhang].]]
|- style="text-align: center; font-size: 75%;"
| colspan="2" | '''Figure 1:''' a) Pontos da malha representativa no plano <math>x</math>-<math>y</math> e (b) as definições dos  elementos de solução SE<math>(i,j,k)</math> e elementos de conservação CE<math>(i,j,k)</math> no ponto <math>(i,j,k)</math> da malha  [2012_Zhang, 1999_Zhang].
|}

Para todo <math display="inline">(x,y,t) \in \mbox{SE}(i,j,k)</math>, aproxima-se <math display="inline">q_m</math> por um polinômio de Taylor de segunda ordem

<span id="eq-3.7"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>q_m^*(x,y,t;i,j,k) =(q_m)_{i,j}^k+[(q_m)_x]_{i,j}^k(x-x_i)+[(q_m)_y]_{i,j}^k(y-y_j)+[(q_m)_t]_{i,j}^k(t-t_k)+[(q_m)_{xt}]_{i,j}  ^k(x-x_i)(t-t_k) </math>
|-
| style="text-align: center;" | <math>   +[(q_m)_{yt}]_{i,j}^k(y-y_j)(t-t_k)+[(q_m)_{xy}]_{i,j}^k(x-x_i)(y-y_j)+\frac{1}{2}[(q_m)_{xx}]_{i,j} ^k(x-x_i)^2 </math>
|-
| style="text-align: center;" | <math>    +\frac{1}{2}[(q_m)_{yy}]_{i,j}^k(y-y_j)^2+\frac{1}{2}[(q_m)_{tt}]_{i,j}^k(t-t_k)^2, \quad m = 1,2,3. </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.7)
|}

Aproximações análogas são feitas sobre as funções <math display="inline">f_m</math> e <math display="inline">g_m</math>, <math display="inline">m = 1,2,3</math>. O termo fonte, por outro lado,  é aproximado por uma ordem a menos:

<span id="eq-3.8"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>S_m^*(x,y,t;i,j,k) =  (S_m)_{i,j}^k+[(S_m)_x]_{i,j}^k(x-x_i)+[(S_m)_y]_{i,j}^k(y-y_j)+[(S_m)_t]_{i,j}^k(t-t_k), \quad m = 1,2,3. </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.8)
|}

Dessa forma, a Eq. ([[#eq-3.5|3.5]]) pode ser  aproximada no SE<math display="inline">(i,j,k)</math> por div '''H'''<math display="inline">_m^* = S_m^*</math>, em  que '''H'''<math display="inline">_m^* = \left(f_m^*, g_m^*, q_m^*\right)</math>, <math display="inline">m = 1,2,3</math>, ao mesmo tempo em que a Eq.  ([[#eq-3.6|3.6]]) é aproximada no CE<math display="inline">(i,j,k)</math> por

<span id="eq-3.9"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>\oint _{S(\textrm{CE}(i,j,k))}\textrm{'''H'''}_m^*\cdot d\textrm{'''s'''} =  \int _{\textrm{CE}(i,j,k)}  S_m^* dV,  \quad  m = 1,2,3. </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.9)
|}

====3.2.2.0 Avaliação de (q_m)_{i,j}^k====

Esto es una prueba.

====3.2.2.1 Avaliação de <math>(q_m)_{i,j}^k</math>====

Substituindo as funções <math display="inline">q_m^*, f_m^*, g_m^*</math> e <math display="inline">S_m^*</math> em ([[#eq-3.9|3.9]]), após diversas simplificações,  obtém-se o esquema de avanço no tempo para a variável <math display="inline">(q_m)_{i,j}^k</math>:

<span id="eq-3.10"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>(q_m)_{i,j}^k + \displaystyle \frac{\Delta x^2}{48}[(q_m)_{xx}]_{i,j}^k+\frac{\Delta  y^2}{48}[(q_m)_{yy}]_{i,j}^k = [P_m^{y}(\Delta x, f_m,g_m)]_{i-1/2,j}^{k-1/2}+[P_m^{y}(-\Delta  x, f_m, g_m)]_{i+1/2,j}^{k-1/2}+ </math>
|-
| style="text-align: center;" | <math>   [P_m^{x}(\Delta y, g_m, f_m)]_{i,j-1/2}^{k-1/2}+[P_m^{x}(-\Delta y, g_m, f_m)]_{i,j+1/2}^{k-1/2}, \quad m  = 1,2,3 \quad \quad  </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.10)
|}

em que

<span id="eq-3.11"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>P_m^{\zeta }(\Delta \xi ,\phi ,\psi ) = \displaystyle  \frac{1}{4}\left\{q_m+\frac{\Delta  \xi }{4}(q_m)_{\xi }+\frac{\Delta t}{4}\left[\phi _{\xi }-\psi _{\zeta }\right]+\frac{2\Delta t}{\Delta  \xi }\left[\phi{+\displaystyle}\frac{\Delta t}{4}\phi _t \right]+\frac{\Delta t^2}{16}\left[\phi _{\xi  t}-\psi _{\zeta t} \right]\right.+ \frac{\Delta t}{48\Delta \xi }\left[\Delta \zeta ^2\phi _{\zeta  \zeta }+\Delta \xi ^2\phi _{\xi \xi }\right]</math>
|-
| style="text-align: center;" | <math>   +\frac{\Delta t}{12\Delta \xi }\left[\Delta t^2\phi _{tt}-\frac{\Delta  \xi ^2}{2}\psi _{\zeta \xi }\right]+ \frac{\Delta \zeta ^2}{192}(q_m)_{\zeta \zeta }+\frac{7\Delta  \xi ^2}{192}(q_m)_{\xi \xi }+\left.\frac{\Delta t}{2}\left[S_m+\frac{\Delta t}{4}(S_m)_t+\frac{\Delta  \xi }{4}(S_m)_{\xi }\right]\right\}.  </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.11)
|}

É importante notar que as Equações ([[#eq-3.10|3.10]]) e ([[#eq-3.11|3.11]]) dependem apenas das incógnitas <math display="inline">q_m,  (q_m)_x, (q_m)_y, (q_m)_{xx}, (q_m)_{yy}</math> e <math display="inline">(q_m)_{xy}</math> do tempo <math display="inline">t_{k-1/2}</math>, para <math display="inline">m = 1,2,3</math>  pois, pelas Equações ([[#eq-3.1|3.1]])-([[#eq-3.4|3.4]]), <math display="inline">f_m</math>, <math display="inline">g_m</math> e <math display="inline">S_m</math> também dependem destas últimas.  Por  outro lado, é preciso conhecer <math display="inline">[(q_m)_{xx}]_{i,j}^k</math> e <math display="inline">[(q_m)_{yy}]_{i,j}^k</math> previamente à obtenção de <math display="inline">(q_m)_{i,j}^k</math> no tempo <math display="inline">t_k</math>, conforme Eq. ([[#eq-3.10|3.10]]). A seção a seguir descreve a obtenção dessas  derivadas espaciais duplas.

====3.2.2.2 Avaliação de <math>[(q_m)_{xx}]_{i,j}^k</math>, <math>[(q_m)_{xy}]_{i,j}^k</math> e <math>[(q_m)_{yy}]_{i,j}^k</math>====

Para prosseguir com o cálculo de <math display="inline">[(q_m)_{xx}]_{i,j}^k</math>, constrói-se primeiramente uma equação auxiliar a  partir da Eq. ([[#eq-3.5|3.5]]), derivando-a duas vezes em relação a <math display="inline">x</math>, obtendo as seguintes formas diferencial  e integral

<span id="eq-3.12"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>\displaystyle \frac{\partial (q_m)_{xx}}{\partial{t}}+\frac{\partial }{\partial x}  \left[(f_m)_{xx}-(S_m)_x\right]+\frac{\partial (g_m)_{xx}}{\partial  y} = \textrm{div }\hat{\textrm{'''H'''}}_m = 0 \quad \textrm{e} \quad  \oint _{S(V)}\hat{\textrm{'''H'''}}_m \cdot d\textrm{'''s'''} = 0, \quad m = 1,2,3. </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.12)
|}

As Equações em ([[#eq-3.12|3.12]]) possuem, conforme Equações ([[#eq-3.7|3.7]]) e ([[#eq-3.8|3.8]]), para todo <math display="inline">(x,y,t) \in </math> SE<math display="inline">(i,j,k)</math>, os seguintes análogos numéricos:

<span id="eq-3.13"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>\textrm{div }\hat{\textrm{'''H'''}}_m^* = 0  \quad \textrm{e} \quad \oint _{S(\textrm{CE}(i,j,k))}\hat{\textrm{'''H'''}}_m^* \cdot  d\textrm{'''s'''} = 0, \quad \quad m = 1,2,3.   </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.13)
|}

O campo vetorial aproximado <math display="inline">\hat{\textrm{'''H'''}}_m^* = ((f_m^*)_{xx}-(S_m^*)_x,  (g_m^*)_{xx},(q_m^*)_{xx})</math> é constante, de modo que a avaliação da Eq.  ([[#eq-3.13|3.13]]) sobre o elemento de conservação descrito pela Fig. [[#img-1a|1a]] retorna, para <math display="inline">m = 1,2,3</math>:

<span id="eq-3.14"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>[(q_m)_{xx}]_{i,j}^k = \frac{1}{4}\left\{[(q_m)_{xx}]_{i-1/2,j}^{k-1/2}+[(q_m)_{xx}]_{i+1/2,j}^{k-1/2}+ [(q_m)_{xx}]_{i,j-1/2}^{k-1/2}+[(q_m)_{xx}]_{i,j+1/2}^{k-1/2}\right\}+\frac{\Delta t}{2 \Delta x}\left\{ [(f_m)_{xx}]_{i-1/2,j}^{k-1/2}\right.</math>
|-
| style="text-align: center;" | <math>  \left.-[(f_m)_{xx}]_{i+1/2,j}^{k-1/2}  +[(S_m)_x]_{i+1/2,j}^{k-1/2}-[(S_m)_x]_{i-1/2,j}^{k-1/2}\right\}+\frac{\Delta t}{2 \Delta y}\left\{ [(g_m)_{xx}]_{i,j-1/2}^{k-1/2}-[(g_m)_{xx}]_{i,j+1/2}^{k-1/2}\right\}.  </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.14)
|}

Para computar as variáveis de marcha <math display="inline">[(q_m)_{yy}]_{i,j}^k</math> e <math display="inline">[(q_m)_{xy}]_{i,j}^k</math> procede-se de modo  análogo, obtendo-se após todos os cálculos e simplificações:

<span id="eq-3.15"></span>
<span id="eq-3.16"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>{[}(q_m)_{yy}{]}_{i,j}^k = \frac{1}{4}\left\{[(q_m)_{yy}]_{i-1/2,j}^{k-1/2}+[(q_m)_{yy}]_{i+1/2,j}^{k-1/2}+ [(q_m)_{yy}]_{i,j-1/2}^{k-1/2}+[(q_m)_{yy}]_{i,j+1/2}^{k-1/2}\right\}+\frac{\Delta t}{2 \Delta y}\left\{ [(g_m)_{yy}]_{i,j-1/2}^{k-1/2}\right.</math>
|-
| style="text-align: center;" | <math>   \left.-[(g_m)_{yy}]_{i,j+1/2}^{k-1/2}+[(S_m)_y]_{i,j+1/2}^{k-1/2}-[(S_m)_y] _ {  i,j-1/2} ^{k-1/2}\right\} +\frac{\Delta t}{2 \Delta x}\left\{ [(f_m)_{yy}]_{i-1/2,j}^{k-1/2}-[(f_m)_{yy}]_{i+1/2,j}^{k-1/2}\right\}, </math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.15)
|-
| style="text-align: center;" | <math> {[}(q_m)_{xy}{]}_{i,j}^k= \frac{1}{4}\left\{ [(q_m)_{xy}]_{i-1/2,j}^{k-1/2}+[(q_m)_{xy}]_{i+1/2,j}^{k-1/2}+ [(q_m)_{xy}]_{i,j-1/2}^{k-1/2}+[(q_m)_{xy}]_{i,j+1/2}^{k-1/2}\right\}</math>
|-
| style="text-align: center;" | <math>  +\frac{\Delta t}{2 \Delta x}\left\{ [(f_m)_{xy}]_{i-1/2,j}^{k-1/2}-[(f_m)_{xy}]_{i+1/2,j}^{k-1/2}+\frac{1}{2}[(S_m)_y]_{i+1/2,j}^{k-1/2}-\frac{1}{ 2}[(S_m)_y] _ { i-1/2 , j } ^{k-1/2}\right\}</math>
|-
| style="text-align: center;" | <math>   +\frac{\Delta t}{2 \Delta y}\left\{[(g_m)_{xy}]_{i,j-1/2}^{k-1/2}-[(g_m)_{xy}]_{i,j+1/2}^{k-1/2}  +\frac{1}{2}[(S_m)_x]_{i,j+1/2}^{k-1/2}-\frac{1}{2}[(S_m)_x] _ {i,j-1/2} ^{k-1/2}\right\}, \quad  m = 1,2,3.  </math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.16)
|}
|}

Aborda-se na seção a seguir o cálculo das derivadas de primeira ordem das variáveis dinâmicas <math display="inline">q_m</math>, <math display="inline">m =  1,2,3</math>.

====3.2.2.3 Avaliação de <math>[(q_m)_{x}]_{i,j}^k</math> e <math>[(q_m)_{y}]_{i,j}^k</math>====

As variáveis <math display="inline">[(q_m)_x]_{i,j}^k</math> e <math display="inline">[(q_m)_y]_{i,j}^k</math> podem ser determinadas por meio de uma estratégia  análoga àquela apresentada na seção [[#3.2.2.2 Avaliação de <math>[(q_m)_{xx}]_{i,j}^k</math>, <math>[(q_m)_{xy}]_{i,j}^k</math> e <math>[(q_m)_{yy}]_{i,j}^k</math>|3.2.2.2]] anterior, bastando, para isso, construir leis diferenciais  de conservação (já nas formas aproximadas)

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>\displaystyle \frac{\partial (q_m^*)_{x}}{\partial{t}}+\frac{\partial (f_m^*)_{x}}{\partial x}  +\frac{\partial (g_m^*)_{x}}{\partial y} =\textrm{div }\check{\textrm{'''H'''}}_m^*= (S_m^*)_x, \quad  \frac{\partial (q_m^*)_{y}}{\partial{t}} +\frac{\partial (f_m^*)_{y}}{\partial x} +\frac{\partial  (g_m^*)_{y}}{\partial y} = \textrm{div }\tilde{\textrm{'''H'''}}_m^* =  (S_m^*)_y </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.17)
|}

e integrais

<span id="eq-3.18"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>\oint _{S(\textrm{CE}(i,j,k))}\check{\textrm{'''H'''}}_m^*\cdot d\textrm{'''s'''}=  \int _{\textrm{CE}(i,j,k)}(S_m^*)_xdV \quad  \textrm{e} \quad  \oint _{S(\textrm{CE}(i,j,k))}\tilde{\textrm{'''H'''}}_m^*\cdot  d\textrm{'''s'''}=  \int _{\textrm{CE}(i,j,k)}(S_m^*)_ydV, \quad m = 1,2,3. </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.18)
|}

As avaliações das equações em ([[#eq-3.18|3.18]]) fornecem as inclinações equacionadas por

<span id="eq-3.19"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>[(q_m)_{\xi }]_{i,j}^k = [Q_m^{(1)}(\Delta x,\xi )]_{i-1/2,j}^{k-1/2}+[Q_m^{(1)}(-\Delta  x,\xi )]_{i+1/2,j}^{k-1/2}+[Q_m^{(2)}(\Delta y,\xi )]_{i,j-1/2}^{k-1/2}+[Q_m^{(2)}(-\Delta  y,\xi )]_{i,j+1/2}^{k-1/2}, </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.19)
|}

para <math display="inline">\xi = x,y</math> e <math display="inline">m = 1,2,3</math>, em que

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>Q_m^{(1)}(\Delta x,\xi ) = \displaystyle  \frac{1}{4}\left\{(q_m)_{\xi }+\frac{\Delta  x}{4}(q_m)_{x\xi }+\frac{\Delta t}{4}\left[(f_m)_{x\xi }-(g_m)_{y\xi }\right]+\frac{2\Delta t}{\Delta  x}\left[(f_m)_{\xi }+\displaystyle \frac{\Delta t}{4}(f_m)_{t\xi }  \right]+\frac{\Delta t}{2} (S_m)_{\xi }\right\},</math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.20)
|-
| style="text-align: center;" | <math> Q_m^{(2)}(\Delta y,\xi ) = \displaystyle  \frac{1}{4}\left\{(q_m)_{\xi }+\frac{\Delta  y}{4}(q_m)_{y\xi }+\frac{\Delta t}{4}\left[(g_m)_{y\xi }-(f_m)_{x\xi }\right]+\frac{2\Delta t}{\Delta  y}\left[(g_m)_{\xi }+\displaystyle \frac{\Delta t}{4}(g_m)_{t\xi }  \right]+\frac{\Delta t}{2}(S_m)_{\xi }\right\}. </math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.21)
|}
|}

O esquema formado pelo conjunto de Equações ([[#eq-3.10|3.10]]),  ([[#eq-3.14|3.14]])-([[#eq-3.16|3.16]]) e  ([[#eq-3.19|3.19]]) é interessante apenas em problemas com soluções  de comportamento suave. Para problemas com descontinuidades ou formação de choque, é necessário modificar o  cálculo das variáveis em ([[#eq-3.14|3.14]])-([[#eq-3.16|3.16]]) e ([[#eq-3.19|3.19]]). A proposta apresentada  neste artigo é uma variação do procedimento adotado por <span id='citeF-37'></span><span id='citeF-38'></span>[[#cite-37|[37,38]]] e consiste num processo  iterativo.

Para prosseguir, considere as seguintes formas regressivas e progressivas:

<span id="eq-3.22"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>[(q_m^{\pm })_{\xi }]_{i,j}^k = \delta _{\xi }^{\pm }(q_m) \quad \textrm{e} \quad  [(q_m^{\pm })_{\zeta \xi }]_{i,j}^k = \delta _{\xi }^{\pm }[(q_m)_{\zeta }], \quad \xi = x, y, \quad \zeta = x,  y \quad \textrm{e} \quad m = 1,2,3 </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.22)
|}

tal que o operador <math display="inline">\delta _{\xi }^{\pm }</math> seja definido como

<span id="eq-3.23"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>\delta _{\xi }^{\pm }(\phi ) =  \displaystyle \left\{\displaystyle \begin{array}{cl}\displaystyle \pm \frac{2}{\Delta \xi }\left[(\phi ')_{i\pm{1}/2,j}^k-\phi _{i,j}^k\right],&  \textrm{se } \;\xi = x\\ \displaystyle \pm \frac{2}{\Delta \xi }\left[(\phi ')_{i,j\pm{1}/2}^k-\phi _{i,j}^k\right],& \textrm{se }\;  \xi = y \end{array}\right.\quad \textrm{e} \quad \phi' = \left\{\begin{array}{cl}\displaystyle \phi +\frac{\Delta t}{2}\phi _{t}+\frac{\Delta t^2}{8}\phi _{tt}, & \textrm{se }\; \phi =  (q_m)_x, (q_m)_y, \\ \displaystyle \phi +\frac{\Delta t}{2}\phi _{t} , & \textrm{se }\; \phi =  (q_m)_{xx}, (q_m)_{xy}, (q_m)_{yy}. \end{array}\right. </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.23)
|}

As aproximações para as derivadas serão definidas como

<span id="eq-3.24"></span>
<span id="eq-3.25"></span>
{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>[(q_m)_{\xi }]_{i,j}^k = W\left\{[(q_m^+)_{\xi }]_{i,j}^k, [(q_m^-)_{\xi }]_{i,j}^k,  \alpha \right\}, </math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.24)
|-
| style="text-align: center;" | <math> {[}(q_m)_{\zeta \xi }{]}_{i,j}^k = W\left\{[(q_m^+)_{\zeta \xi }]_{i,j}^k,  [(q_m^-)_{\zeta \xi }]_{i,j}^k, \alpha \right\}, \quad m = 1,2,3, \quad \xi = x, y, \quad \zeta = x,  y \quad \textrm{e}  \quad  \alpha \geq 0,  </math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.25)
|}
|}

e a função <math display="inline">W</math> <span id='citeF-26'></span><span id='citeF-24'></span>[[#cite-26|[26,24]]]:

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>W(x_{-}, x_{+}, \alpha ) =   \left\{  \begin{array}{cl}\displaystyle \frac{|x_{+}|^{\alpha }x_{-}+|x_{-}|^{\alpha }x_{+}}{|x_{+}|^{\alpha }+|x_{-}|^{\alpha }}, &  \textrm{se} \quad |x_{+}|^{\alpha }+|x_{-}|^{\alpha } \neq 0, \\ 0, & \textrm{caso contrario.} \end{array} \right. </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.26)
|}

Observe que neste caso, as Equações ([[#eq-3.24|3.24]]) e ([[#eq-3.25|3.25]]) representam médias ponderadas das diferenças  progressivas e regressivas descritas em ([[#eq-3.22|3.22]]) e ([[#eq-3.23|3.23]]). No caso de <math display="inline">\alpha = 0</math>,  ([[#eq-3.24|3.24]]) e ([[#eq-3.25|3.25]]) tornam-se diferenças finitas centrais <span id='citeF-23'></span>[[#cite-23|[23]]]. Adota-se neste trabalho <math display="inline">\alpha =1</math>.

Importante notar, no entanto, que as Equações ([[#eq-3.24|3.24]]) e ([[#eq-3.25|3.25]]) dependem da variável <math display="inline">(q_m)_{i,j}^k</math> que, por sua vez, depende de <math display="inline">[(q_m)_{xx}]_{i,j}^k</math> e <math display="inline">[(q_m)_{yy}]_{i,j}^k</math>, conforme Eq.  ([[#eq-3.10|3.10]]). Dito de outra forma, para obter <math display="inline">(q_m)_{i,j}^k</math> é necessário obter antes <math display="inline">[(q_m)_{xx}]_{i,j}^k</math> e <math display="inline">[(q_m)_{yy}]_{i,j}^k</math>, e vice-versa. Para sanar esta dificuldade, construímos o  processo iterativo descrito no Algoritmo [[#algorithm-1|1]].


{| style="margin: 1em auto;border: 1px solid darkgray;"
|-
|
:Defina uma tolerância <math>\tau </math>
:      Faça <math>m = 1</math>
: <math>m = 3</math>      Faça <math>l = 1</math>
:      Calcule <math>[(q_m)_{xx}]_{i,j}^k</math> e <math>[(q_m)_{yy}]_{i,j}^k</math> com as Equações ([[#eq-3.14|3.14]]) e  ([[#eq-3.15|3.15]]), respectivamente;
|-
| Calcule <math>(q_m^l)_{i,j}^k</math> com a  Eq. ([[#eq-3.10|3.10]]);
|-
| <math>||(q_m^{l+1})_{i,j}^k-(q_m^l)_{i,j}^k||< \tau </math>      Calcule <math>[(q_m)_x]_{i,j}^k</math> e <math>[(q_m)_y]_{i,j}^k</math> com a Eq. ([[#eq-3.24|3.24]]);
|-
| Calcule <math>[(q_m)_{xx}]_{i,j}^k</math>, <math>[(q_m)_{xy}]_{i,j}^k</math> e <math>[(q_m)_{yy}]_{i,j}^k</math> com  a Eq.([[#eq-3.25|3.25]]);
|-
| Calcule <math>(q_m^{l+1})_{i,j}^k</math> com a  Eq. ([[#eq-3.10|3.10]]);            


|-
| style="text-align: center; font-size: 75%;"|
<span id='algorithm-1'></span>'''Algorithm. 1''' Algoritmo CE/SE para o cálculo das variáveis de marcha sobre um determinado ponto <math>(i,j,k)</math> da  malha.
|}

A estabilidade do esquema CE/SE satisfaz a condição de Courant-Friedrichs-Lewy (CFL) que, para  as equações de Saint Venant unidimensional e bidimensional são as restrições <span id='citeF-23'></span>[[#cite-23|[23]]]

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>CFL = \displaystyle \frac{\Delta t}{\Delta x}\textrm{max}\left(|u|+\sqrt{gh}\right)\quad \textrm{e} \quad  CFL = \displaystyle \textrm{max}\left\{\frac{\Delta t}{\Delta x}\left(|u|+\sqrt{gh}\right),  \frac{\Delta t}{\Delta y} \left(|v|+\sqrt{gh}\right)\right\}, </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.27)
|}

respectivamente, onde o número de Courant satisfaz <math display="inline">0 < CFL \leq 1</math>. Maiores detalhes sobre  estabilidade do método CE/SE pode ser encontrado em <span id='citeF-39'></span>[[#cite-39|[39]]].

==3.3 Exemplos Numéricos==

Esta seção tem por objetivo avaliar e validar o esquema numérico desenvolvido. Para isso, alguns problemas  teste são avaliados, como testes de acurácia em uma e duas dimensões, bem como problemas clássicos de ruptura  de barragem unidimensionais (subseção [[#3.3.1 Exemplos Unidimensionais|3.3.1]]) e bidimensionais (subseção [[#3.3.2 Exemplos Bidimensionais|3.3.2]]).

===3.3.1 Exemplos Unidimensionais===

<span id='theorem-ex1'></span>Exemplo 1:  O objetivo deste exemplo teste é verificar experimentalmente a ordem de acurácia do esquema  CE/SE. Considere o sistema hiperbólico linear unidimensional

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>  \textrm{'''U'''}_t+\textrm{'''F'''}(\textrm{'''U'''})_x = \textrm{'''0'''}, \quad  \textrm{onde} \quad \textrm{'''U'''} = (\phi , u), \quad \textrm{'''F'''} = (cu, c\phi ) </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.28)
|}

e <math>c</math> é uma constante dada. Impondo as condições inicias <math>\phi (x,0) =  -c^{-1}\sin (2\pi x)</math> e <math>u(x,0) = 0</math>, é possível obter a seguinte solução exata <span id='citeF-40'></span><span id='citeF-41'></span>[[#cite-40|[40,41]]]:

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>  \phi (x,t) = -c^{-1}\cos (2\pi c t)\sin (2\pi x) \quad \textrm{e} \quad u(x,t) = c^{-1}\sin (2\pi ct)\cos (2\pi  x). </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.29)
|}

As soluções são computadas com as Equações ([[#eq-3.10|3.10]]), ([[#eq-3.14|3.14]])-([[#eq-3.16|3.16]]) e  ([[#eq-3.19|3.19]]). A  Tabela [[#table-1|1]] apresenta os erros calculados nas normas <math>L^1</math> e <math>L^2</math>, nas simulações  realizadas sobre o domínio computacional <math>-1 \leq x \leq 1</math>, <math>c = 1</math>, <math>CFL = 0.8</math>  até atingir o tempo <math>t = 1.2</math>s, sendo <math>n</math> o número de células no espaço. O experimento verifica a terceira  ordem de acurácia do esquema para ambas as variáveis <math>u</math> e <math>\phi </math>.      
{|  class="floating_tableSCP wikitable" style="text-align: center; margin: 1em auto;min-width:50%;"
|+ style="font-size: 75%;" |<span id='table-1'></span>Table. 1  <span style="text-align: center; font-size: 75%;">Ordem de acurácia experimental do método CE/SE para o sistema hiperbólico linear  unidimensional.</span>
|- style="border-top: 2px solid;"
| rowspan='2' style="border-bottom: 2px solid;" | <math>n</math>
| colspan='4' | Norma <math>L^1</math> 
|    
| colspan='4' |  Norma <math>L^2</math>
|-
|  Erro (<math>\phi </math>) 
| Ordem 
|  Erro (<math>u</math>)
|  Ordem  
| 
|  Erro (<math>\phi </math>)
| Ordem 
|  Erro  (<math>u</math>)
|   Ordem 
|- style="border-top: 2px solid;"
|  20 
| 8.322<math>\times 10^{-3}</math> 
| &#8211; 
| 1.662<math>\times 10^{-2}</math> 
| &#8211; 
|  
| 7.157<math>\times 10^{-3}</math> 
| &#8211; 
| 1.347<math>\times 10^{-2}</math> 
| &#8211; 
|-
|  40 
| 1.849<math>\times 10^{-3}</math> 
|  2.170 
| 2.524<math>\times 10^{-3}</math> 
| 2.719
|  
| 1.488<math>\times 10^{-3}</math>  
| 2.266 
| 2.042<math>\times 10^{-3}</math> 
| 2.722
|-
|  80 
| 3.067<math>\times 10^{-4}</math> 
|  2.592 
| 3.401<math>\times 10^{-4}</math> 
| 2.892
|  
| 2.443<math>\times 10^{-4}</math>  
| 2.607 
| 2.715<math>\times 10^{-4}</math> 
| 2.911
|-
|  160 
| 4.369<math>\times 10^{-5}</math> 
|  2.811 
| 4.386<math>\times 10^{-5}</math> 
| 2.955
|  
| 3.493<math>\times 10^{-5}</math>  
| 2.806 
| 3.465<math>\times 10^{-5}</math> 
| 2.970
|-
|  320 
| 5.813<math>\times 10^{-6}</math> 
|  2.910 
| 5.566<math>\times 10^{-6}</math> 
| 2.978
|  
| 4.665<math>\times 10^{-6}</math>  
| 2.904 
| 4.369<math>\times 10^{-6}</math> 
| 2.987
|- style="border-bottom: 2px solid;"
|  640 
| 7.490<math>\times 10^{-7}</math> 
|  2.956 
| 7.010<math>\times 10^{-7}</math> 
| 2.989
|  
| 6.030<math>\times 10^{-7}</math>  
| 2.952 
| 5.490<math>\times 10^{-7}</math> 
| 2.992

|}

<span id='theorem-fig_DBS01'></span>Exemplo 2:  Considera-se agora um problema ideal de ruptura de barragem sobre um domínio molhado, isto é, a quebra  de barragem é instantânea, o fundo é plano e não existe resistência ao escoamento. As condições iniciais para  esta configuração seguem o clássico problema de Riemann

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>  h(x,0) =   \left\{   \begin{array}{cl}    h_l, & \textrm{ para } 0 \textrm{m} \leq x \leq x_0 \textrm{m}\\    h_r , & \textrm{ para } x_0 \textrm{m} < x \leq L \textrm{m}   \end{array}  \right., \quad \textrm{ com } \quad h_l \geq h_r \quad \textrm{e} \quad u(x,0) = 0 \textrm{m/s}. </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.30)
|}

O domínio considerado é <math>0\mbox{m}</math> <math>\leq x \leq 2000\mbox{m}</math>, <math>x_0 = 1000\mbox{m}</math> e a solução analítica  pode ser encontrada em <span id='citeF-42'></span>[[#cite-42|[42]]] ou <span id='citeF-43'></span>[[#cite-43|[43]]]. Utiliza-se, para fins de análise  de comportamento, a mesma relação de equações do Ex. [[#theorem-ex1|1]] anterior, isto é, Equações ([[#eq-3.10|3.10]]),  ([[#eq-3.14|3.14]])-([[#eq-3.16|3.16]]) e ([[#eq-3.19|3.19]]). As Fig. [[#theorem-fig_DBS01|2]] e  Fig. [[#img-2a|2a]] demonstram o comportamento da solução numérica em relação à analítica, calculadas no  tempo <math>t = 52\mbox{s}</math>, com <math>n+1 = 201</math> pontos, incremento espacial <math>\Delta x = 10\mbox{m}</math>, <math>CFL = 0.8</math> e  profundidades iniciais a montante e a jusante <math>h_l = 10\mbox{m}</math> e <math>h_r = 5\mbox{m}</math>, respectivamente. Note  que a resposta numérica é coerente com a analítica, embora apresente suavidade que a distancie nas regiões  com mudanças abruptas.   <div id='img-2a'></div>
<div id='img-2'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:draft_García_699893715-fig_DB_S01.png|414px|]]
|[[Image:draft_García_699893715-fig_DB_S11.png|414px|Elevação da superfície da água h e velocidade u para o problema de ruptura de  barragem, computada no tempo t = 52\mboxs, utilizando uma malha uniforme com 201  pontos, incremento espacial ∆x = 10\mboxm e h<sub>r</sub>/hₗ= 0.5.]]
|- style="text-align: center; font-size: 75%;"
| colspan="2" | '''Figure 2:''' Elevação da superfície da água <math>h</math> e velocidade <math>u</math> para o problema de ruptura de  barragem, computada no tempo <math>t = 52\mbox{s}</math>, utilizando uma malha uniforme com <math>201</math>  pontos, incremento espacial <math>\Delta x = 10\mbox{m}</math> e <math>h_r/h_l = 0.5</math>.
|}
Utilizando o Algoritmo [[#algorithm-1|1]], sob as mesmas condições e com os mesmos parâmetros, obtém-se os  gráficos constantes na  Fig. [[#img-3|3]]. Observa-se que esta solução numérica é superior à anterior,  sobretudo no que tange às regiões de rápidas mudanças.  Conforme Zhang ''et al.'' (2012) <span id='citeF-23'></span>[[#cite-23|[23]]],  a razão <math>h_r/h_l</math> é um importante índice para  julgar a aplicabilidade e a acurácia de esquemas numéricos no modelo 1D de ruptura de barragem. Segundo os  mesmos autores, os regimes de escoamentos subcrítico e supercrítico existem simultaneamente num canal  sem fricção, horizontal e retangular, quando <math>h_r/h_l < 0.138</math>.  Altera-se, neste sentido, estes  parâmetros para uma razão <math>h_r/h_l = 0.001</math>. Os resultados simulados são mostrados nas Fig. fig_sv1d11  e Fig. [[#img-3a|3a]].   <div id='img-3a'></div>
<div id='img-3'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:draft_García_699893715-fig_DB_S0.png|414px|]]
|[[Image:draft_García_699893715-fig_DB_S1.png|414px|]]
|-
|[[Image:draft_García_699893715-fig_DB_S02.png|414px|]]
|[[Image:draft_García_699893715-fig_DB_S12.png|414px|Elevação da superfície da água h e velocidade u para o problema de ruptura de  barragem, no tempo t = 52\mboxs, numa malha uniforme com 201 pontos e respectivo incremento espacial ∆x = 10\mboxm, calculados com o Algoritmo [[#algorithm-1|1]], sendo: (a)-(b) hₗ= 10\mboxm e h<sub>r</sub>= 5\mboxm, (c)-(d) hₗ= 10\mboxm e h<sub>r</sub>= 0.01\mboxm.]]
|- style="text-align: center; font-size: 75%;"
| colspan="2" | '''Figure 3:''' Elevação da superfície da água <math>h</math> e velocidade <math>u</math> para o problema de ruptura de  barragem, no tempo <math>t = 52\mbox{s}</math>, numa malha uniforme com <math>201</math> pontos e respectivo incremento espacial <math>\Delta x = 10\mbox{m}</math>, calculados com o Algoritmo [[#algorithm-1|1]], sendo: (a)-(b) <math>h_l = 10\mbox{m}</math> e <math>h_r = 5\mbox{m}</math>, (c)-(d) <math>h_l = 10\mbox{m}</math> e <math>h_r = 0.01\mbox{m}</math>.
|}

===3.3.2 Exemplos Bidimensionais===

<span id='theorem-tabela02'></span>Exemplo 3:  Considere agora o sistema hiperbólico linear bidimensional

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>  \textrm{'''U'''}_t+\textrm{'''F'''}(\textrm{'''U'''})_x  +\textrm{'''G'''}(\textrm{'''U'''})_y= \textrm{'''0'''}, \quad  \textrm{onde} \quad \textrm{'''U'''} = (\phi , u,v), \quad \textrm{'''F'''} = (cu, c\phi ,0), \quad  \textrm{'''G'''} = (cv, 0,c\phi ) </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.31)
|}

e <math>c</math> é uma constante dada. Para os dados iniciais <math>\phi (x,y,0) =  -c^{-1}[\sin (2\pi x)+\sin (2\pi y)]</math> e <math>u(x,y,0) = v(x,y,0) = 0</math>, este problema admite a seguinte solução  exata <span id='citeF-40'></span><span id='citeF-41'></span>[[#cite-40|[40,41]]]:

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>  \phi = -c^{-1}\cos (2\pi c t)[\sin (2\pi x)+\sin (2\pi y)], \quad u = c^{-1}\sin (2\pi  ct)\cos (2\pi x) \quad \textrm{e} \quad  v = c^{-1}\sin (2\pi ct)\cos (2\pi y). </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.32)
|}


{|  class="floating_tableSCP wikitable" style="text-align: center; margin: 1em auto;min-width:50%;"
|+ style="font-size: 75%;" |Table. 2  <span style="text-align: center; font-size: 75%;">Ordem de acurácia experimental do método CE/SE para o sistema hiperbólico linear  bidimensional.</span>
|- style="border-top: 2px solid;"
| rowspan='2' style="border-bottom: 2px solid;" | <math>n_x \times n_y</math>
| colspan='4' | Norma <math>L^1</math> 
|    
| colspan='4' |  Norma <math>L^2</math>
|-
|  Erro (<math>\phi </math>) 
| Ordem 
|  Erro (<math>u, v</math>)
|  Ordem  
| 
|  Erro (<math>\phi </math>)
| Ordem 
|  Erro  (<math>u, v</math>)
|   Ordem 
|- style="border-top: 2px solid;"
| <math>20\times 20</math> 
|  2.595<math>\times 10^{-3}</math> 
| &#8211; 
|  1.283<math>\times 10^{-2}</math> 
| &#8211; 
|  
| 1.797<math>\times 10^{-3}</math> 
| &#8211; 
|  7.478<math>\times 10^{-3}</math> 
| &#8211;  
|-
| <math>40\times 40</math> 
|  4.231<math>\times 10^{-4}</math> 
|  2.617 
|  2.032<math>\times 10^{-3}</math> 
|  2.658 
|  
| 3.377<math>\times 10^{-4}</math> 
|   2.412 
|  1.169<math>\times 10^{-3}</math> 
|  2.678 
|-
| <math>80\times 80</math> 
|  7.519<math>\times 10^{-5}</math> 
|  2.492 
|  2.924<math>\times 10^{-4}</math> 
|  2.797 
|  
| 6.111<math>\times 10^{-5}</math> 
|   2.466 
|  1.653<math>\times 10^{-4}</math> 
|    2.822 
|-
| <math>160 \times 160</math> 
|  1.197<math>\times 10^{-5}</math> 
|  2.651 
|  3.938<math>\times 10^{-5}</math> 
|  2.892 
|  
| 9.602<math>\times  10^{-6}</math> 
|  2.670 
|  2.217<math>\times 10^{-5}</math> 
|  2.898 
|-
| <math>320\times 320</math> 
|  1.735<math>\times 10^{-6}</math> 
|  2.786 
|  5.128<math>\times 10^{-6}</math> 
|  2.941 
|  
| 1.397<math>\times 10^{-6}</math> 
|   2.781 
|  2.895<math>\times 10^{-6}</math> 
|  2.937 
|- style="border-bottom: 2px solid;"
| <math>640\times 640</math> 
|  2.364<math>\times 10^{-7}</math> 
|  2.876 
|  6.554<math>\times 10^{-7}</math> 
|  2.968 
|  
| 1.939<math>\times 10^{-7}</math> 
|   2.849 
|  3.723<math>\times 10^{-7}</math> 
|  2.959 

|}
A Tabela [[#theorem-tabela02|3]] apresenta os erros numéricos calculados nas normas <math>L^1</math> e <math>L^2</math>, no tempo <math>t =  0.2</math>s, com parâmetros especificados em <math>c = 1</math>, <math>CFL = 0.4</math> e domínio computacional definido em <math>-1 \leq x, y  \leq 1</math>, sendo <math>n_x \times n_y</math> o número de células utilizadas na discretização. O experimento confirma  novamente a acurácia de terceira ordem do esquema CE/SE.

<span id='theorem-fig_RB2D'></span>Exemplo 4: Este problema hipotético bidimensional é um exemplo  utilizado por <span id='citeF-23'></span><span id='citeF-44'></span><span id='citeF-45'></span>[[#cite-23|[23,44,45]]]. Neste problema as velocidades iniciais são todas  nulas, a profundidade a montante é de <math display="inline">10</math>m, enquanto que a profundidade a jusante é  assumida ser <math display="inline">5</math>m ou <math display="inline">0.1</math>m. O domínio computacional consiste de uma região de <math display="inline">200\mbox{m}\times  200\mbox{m}</math>, com uma parede que se estende paralelamente ao eixo <math display="inline">y</math>, tendo <math display="inline">10</math>m de largura e está centrada  em <math display="inline">x = 100</math>m. A falha é suposta ser instantânea, possui <math display="inline">75</math>m de extensão a partir de <math display="inline">y = 95</math>m. O canal é  horizontal e desconsidera-se resistência ao escoamento. Espera-se a formação de uma frente  de choque após o rompimento. Os resultados em <math display="inline">t = 7.2</math>s são computados a  partir de uma malha uniforme composta por <math display="inline">101 \times 101</math> células e <math display="inline">CFL = 0.4</math>. Constam nas Figuras  [[#theorem-fig_RB2D|4]] e fig_CVRB2D gráficos da profundidade, vetores velocidade e curvas de nível para o  problema com profundidade inicial a jusante de <math display="inline">5</math>m, enquanto que as Figuras fig_RB2D1 e  [[#img-4a|4a]] referem-se ao problema com profundidade inicial a jusante de <math display="inline">0.1</math>m. Os resultados são  consistentes com aqueles presentes na literatura <span id='citeF-23'></span><span id='citeF-44'></span><span id='citeF-45'></span>[[#cite-23|[23,44,45]]].

<div id='img-4a'></div>
<div id='img-4'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:draft_García_699893715-fig_RB2D.png|570px|]]
|[[Image:draft_García_699893715-fig_CV_RB2D.png|414px|]]
|-
|[[Image:draft_García_699893715-fig_RB2D1.png|570px|]]
|[[Image:draft_García_699893715-fig_CV1_RB2D.png|414px|Elevação da superfície da água h, (a) e (c), vetores velocidade e curvas de nível, (b) e  (d), para a solução calculada no tempo t = 7.2\mboxs com CFL = 0.4, do problema de ruptura de  barragem anti-simétrica em um domínio horizontal e sem fricção, com profundidades a montante de 10m e a  jusante de 5m (a)-(b) e 0.1m (c)-(d).]]
|- style="text-align: center; font-size: 75%;"
| colspan="2" | '''Figure 4:''' Elevação da superfície da água <math>h</math>, (a) e (c), vetores velocidade e curvas de nível, (b) e  (d), para a solução calculada no tempo <math>t = 7.2\mbox{s}</math> com <math>CFL = 0.4</math>, do problema de ruptura de  barragem anti-simétrica em um domínio horizontal e sem fricção, com profundidades a montante de <math>10</math>m e a  jusante de <math>5</math>m (a)-(b) e <math>0.1</math>m (c)-(d).
|}

<span id='theorem-fig_RBC2D'></span>Exemplo 5: Este problema teste visa avaliar a  habilidade do esquema em preservar simetria. Considera-se um domínio de <math display="inline">50\mbox{m}\times 50\mbox{m}</math> com  condições iniciais:

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>h(x,y,0) = \left\{   \begin{array}{cl}10\mbox{m}, & \mbox{ se } (x-25)^2+(y-25)^2\leq 10.5^2\\     1\mbox{ m}, & \mbox{ caso contrario}   \end{array} \right.\quad \textrm{e} \quad \quad u(x,y,0) = v(x,y,0) = 0. </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3.33)
|}

No instante da falha da barragem, supõe-se que a parede circular seja removida completamente e,  subsequentemente, formam-se ondas que se espalham radialmente. A solução numérica é computada numa malha  retangular uniforme com <math display="inline">101\times 101</math> células e o passo de tempo é tal que <math display="inline">CFL = 0.6</math>. A Fig.  [[#theorem-fig_RBC2D|5]] apresenta o perfil da superfície da água <math display="inline">0.72</math>s após a hipotética falha na barragem  circular. Vetores velocidade e curvas nível relativas à superfície <math display="inline">h</math> estão dispostos na Fig.  [[#img-5a|5a]]. A simetria da solução numérica é bem preservada e está de acordo com a  literatura <span id='citeF-7'></span><span id='citeF-46'></span>[[#cite-7|[7,46]]].

<div id='img-5a'></div>
<div id='img-5'></div>
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
|-
|[[Image:draft_García_699893715-fig_RBC2D.png|570px|]]
|[[Image:draft_García_699893715-fig_CV1_RBC2D.png|414px|Perfil da superfície da água h (a), vetores velocidade e curvas de nível  (b), para a  solução calculada no tempo t = 0.72\mboxs com CFL = 0.6, do problema de ruptura de  barragem circular em um domínio horizontal e sem fricção.]]
|- style="text-align: center; font-size: 75%;"
| colspan="2" | '''Figure 5:''' Perfil da superfície da água <math>h</math> (a), vetores velocidade e curvas de nível  (b), para a  solução calculada no tempo <math>t = 0.72\mbox{s}</math> com <math>CFL = 0.6</math>, do problema de ruptura de  barragem circular em um domínio horizontal e sem fricção.
|}

==3.4 Conclusões==

Este artigo apresentou o desenvolvimento de um novo esquema CE/SE explícito para a solução das equações de  águas rasas em uma e duas dimensões. Na formulação construída, as variáveis dinâmicas e, consequentemente, as  leis diferencial e integral de conservação, foram aproximadas localmente por expansões de Taylor de segunda  ordem. O conjunto de variáveis de marcha tornou-se constituído pelas variáveis fluxo <math display="inline">q_m</math> e suas  derivadas espaciais de primeira <math display="inline">(q_m)_x</math>, <math display="inline">(q_m)_y</math> e segunda ordem <math display="inline">(q_m)_{xx}, (q_m)_{xy}</math> e <math display="inline">(q_m)_{yy}</math>.  Os experimentos computacionais realizados indicaram acurácia de terceira ordem sobre os problemas  hiperbólicos lineares testados. As soluções numéricas dos problemas teste de ruptura de barragem  unidimensional e bidimensional, caracterizados pela formação de choque, descontinuidade ou simetria,  apresentaram concordância com as soluções analíticas e/ou com a literatura. Conclui-se, com isso, que o  esquema proposto possui considerável habilidade em capturar choques e descontinuidades fazendo com que seja  uma boa ferramenta para análise de fluxo de ruptura de barragem.

==Agradecimentos==

Os autores agradecem à CAPES pelo apoio financeiro, ao Programa de Pós-Graduação em Métodos Numéricos em Engenharia (PPGMNE - UFPR) e ao Instituto Federal de Educação, Ciência e Tecnologia Catarinense, Campus Araquari, pelo apoio à pesquisa.

==3.5 Referências==

===BIBLIOGRAPHY===

<div id="cite-1"></div>
'''[[#citeF-1|[1]]]''' Yang, Suo and Kurganov, Alexander and Liu, Yingjie. (2015) "Well-Balanced Central Schemes on Overlapping Cells with Constant Subtraction Techniques for the Saint-Venant Shallow Water System", Volume 63. Journal of Scientific Computing 3 678 - 698

<div id="cite-2"></div>
'''[[#citeF-2|[2]]]''' Jaswant Singh and Mustafa S. Altinakar and Yan Ding. (2011) "Two-dimensional numerical modeling of dam-break flows over natural terrain using a central explicit scheme", Volume 34. Advances in Water Resources 10 1366 - 1375

<div id="cite-3"></div>
'''[[#citeF-3|[3]]]''' Soler, Joan and Bladé, E and Sánchez-Juny, M. (2012) "Ensayo comparativo entre modelos unidimensionales y bidimensionales en la modelización de la rotura de una balsa de materiales sueltos erosionables", Volume 28. Revista Internacional de Métodos Numéricos para Cálculo y Diseño en Ingeniería 2 103 - 111

<div id="cite-4"></div>
'''[[#citeF-4|[4]]]''' Hatice Ozmen-Cagatay and Selahattin Kocaman and Hasan Guzel. (2014) "Investigation of dam-break flood waves in a dry channel with a hump", Volume 8. Journal of Hydro-environment Research 3 304 - 315

<div id="cite-5"></div>
'''[[#citeF-5|[5]]]''' Selahattin Kocaman and Hatice Ozmen-Cagatay. (2015) "Investigation of dam-break induced shock waves impact on a vertical wall", Volume 525. Journal of Hydrology  1 - 12

<div id="cite-6"></div>
'''[[#citeF-6|[6]]]''' Tsang-Jung Chang and Hong-Ming Kao and Kao-Hua Chang and Ming-Hsi Hsu. (2011) "Numerical simulation of shallow-water dam break flows in open channels using smoothed particle hydrodynamics", Volume 408. Journal of Hydrology 1-2 78 - 90

<div id="cite-7"></div>
'''[[#citeF-7|[7]]]''' Akoh, R and Li, S and Xiao, F. (2008) "A CIP/multi-moment finite volume method for shallow water equations with source terms", Volume 56. Int. J. Numer. Meth. Fluids  2245 - 2270

<div id="cite-8"></div>
'''[[#citeF-8|[8]]]''' Kunihiko Toda and Youichi Ogata and Takashi Yabe. (2009) "Multi-dimensional conservative semi-Lagrangian method of characteristics CIP for the shallow water equations ", Volume 228. Journal of Computational Physics 13 4917 - 4944

<div id="cite-9"></div>
'''[[#citeF-9|[9]]]''' Raul Borsche. (2015) "A well-balanced solver for the Saint Venant equations with variable cross-section", Volume 23. Journal of Numerical Mathematics 2 99 - 115

<div id="cite-10"></div>
'''[[#citeF-10|[10]]]''' Ji-Wen Wang and Ru-Xun Liu. (2005) "Combined finite volume-finite element method for shallow water equations", Volume 34. Computers & Fluids  10 1199 - 1222

<div id="cite-11"></div>
'''[[#citeF-11|[11]]]''' M.J. Castro Díaz and J.A. López-García and Carlos Parés. (2013) "High order exactly well-balanced numerical methods for shallow water systems", Volume 246. Journal of Computational Physics  242 - 264

<div id="cite-12"></div>
'''[[#citeF-12|[12]]]''' Fayssal Benkhaldoun and Saida Sari and Mohammed Seaid. (2015) "Projection finite volume method for shallow water flows", Volume 118. Mathematics and Computers in Simulation  87 - 101

<div id="cite-13"></div>
'''[[#citeF-13|[13]]]''' Liang, Shin-Jye and Hsu, Tai-Wen. (2009) "Least-squares finite-element method for shallow-water equations with source terms", Volume 25. Acta Mechanica Sinica 5 597&#8211;610

<div id="cite-14"></div>
'''[[#citeF-14|[14]]]''' Clint Dawson and Juha H. Videman. (2013) "A streamline diffusion finite element method for the viscous shallow water equations", Volume 251. Journal of Computational and Applied Mathematics  1 - 7

<div id="cite-15"></div>
'''[[#citeF-15|[15]]]''' Georges Kesserwani and Daniel Caviedes-Voullieme and Nils Gerhard and Siegfried Müller. (2015) "Multiwavelet discontinuous Galerkin h-adaptive shallow water model", Volume 294. Computer Methods in Applied Mechanics and Engineering  56 - 71

<div id="cite-16"></div>
'''[[#citeF-16|[16]]]''' D. Wirasaet and S.R. Brus and C.E. Michoski and E.J. Kubatko and J.J. Westerink and C. Dawson. (2015) "Artificial boundary layers in discontinuous Galerkin solutions to shallow water equations in channels ", Volume 299. Journal of Computational Physics  597 - 612

<div id="cite-17"></div>
'''[[#citeF-17|[17]]]''' Mahir Rasulov and Zafer Aslan and Ozkan Pakdil. (2005) "Finite differences method for shallow water equations in a class of discontinuous functions", Volume 160. Applied Mathematics and Computation 2 343 - 353

<div id="cite-18"></div>
'''[[#citeF-18|[18]]]''' Zhengfu, Xu and Chi-Wang, Shu. (2006) "ANTI-DIFFUSIVE FINITE DIFFERENCE WENO METHODS FOR SHALLOW WATER WITH TRANSPORT OF POLLUTANT", Volume 24. Journal of Computational Mathematics 3 239 - 251

<div id="cite-19"></div>
'''[[#citeF-19|[19]]]''' Yao-Hsin Hwang. (2013) "A characteristic particle method for the Saint Venant equations ", Volume 76. Computers & Fluids   58 - 72

<div id="cite-20"></div>
'''[[#citeF-20|[20]]]''' Haifei Liu and Hongda Wang and Shu Liu and Changwei Hu and Yu Ding and Jie Zhang. (2015) "Lattice Boltzmann method for the Saint-Venant equations ", Volume 524. Journal of Hydrology  411 - 416

<div id="cite-21"></div>
'''[[#citeF-21|[21]]]''' Chang, Sin-Chung and To, Wai-Wing. (1991) "A New Numerical Framework for Solving Conservation Laws - The Method of Space-Time Conservation Element and Solution Element". NASA TM 104495

<div id="cite-22"></div>
'''[[#citeF-22|[22]]]''' Thomas Molls and Frank Molls. (1998) "Space-Time Conservation Method Applied to Saint Venant Equations", Volume 124. Journal of Hydraulic Engineering 5 501-508

<div id="cite-23"></div>
'''[[#citeF-23|[23]]]''' Zhang, Yongxiang and Zeng, Zhong and Chen, Jingqiu "The improved space-time conservation element and solution element scheme for two-dimensional dam-break flow simulation", Volume 68. International Journal for Numerical Methods in Fluids 5

<div id="cite-24"></div>
'''[[#citeF-24|[24]]]''' Shamsul Qamar and Saqib Zia and Waqas Ashraf. (2014) "The space-time CE/SE method for solving single and two-phase shallow flow models", Volume 96. Computers & Fluids  136 - 151

<div id="cite-25"></div>
'''[[#citeF-25|[25]]]''' Chang, Sin-Chung. (1993) "New Developments in the Method of Space-Time Conservation Element and Solution Element - Aplications to the Euler and Navier-Stokes Equations". NASA TM 106226

<div id="cite-26"></div>
'''[[#citeF-26|[26]]]''' Chang, Sin-Chung. (1995) "The Method of Space-Time Conservation Element and Solution Elemen - A New Approach for Solving the Navier-Stokes and Euler Equations", Volume 119. Journal of Computational Physics 2 295&#8211;324

<div id="cite-27"></div>
'''[[#citeF-27|[27]]]''' Chang, Sin-Chung and Wang, Xiao-Yen and Chow, Chuen-Yen. (1999) "The Space-Time Conservation Element and Solution Element Method: A New High-Resolution and Genuinely Multidimensional Paradigm for Solving Conservation Laws", Volume 156. Journal of Computational Physics 1 89&#8211;136

<div id="cite-28"></div>
'''[[#citeF-28|[28]]]''' Chang, Sin-Chung and Wang, Xiao-Yen and To, Wai-Ming. (2000) "Application of the Space-Time Conservation Element and Solution Element Method to One-Dimensional Convection-Diffusion Problems", Volume 165. Journal of Computational Physics 1 189&#8211;215

<div id="cite-29"></div>
'''[[#citeF-29|[29]]]''' Zhang, Zeng-Chan and John Yu, S. T. and Wang, Xiao-Yen and Chang, Sin-Chung and Himansu, Ananda and Jorgenson, Philip C. E. (2000) "The CE/SE Method for Navier-Stokes Equations Using Unstructured Meshes for Flows at All Speeds", Volume . AIAA 2000-0393  7

<div id="cite-30"></div>
'''[[#citeF-30|[30]]]''' Zhang, Zeng-Chan and Yu, S. T. John and Chang, Sin-Chung. (2002) "A Space-time Conservation Element and Solution Element Method for Solving the Two- and Three-dimensional Unsteady Euler Equations Using Quadrilateral and Hexahedral Meshes", Volume 175. Academic Press Professional, Inc. J. Comput. Phys. 1 168&#8211;199

<div id="cite-31"></div>
'''[[#citeF-31|[31]]]''' Chang, Sin-Chung. (2010) "A New Approach for Constructing Highly Stable High Order CESE Schemes". NASA TM 2010-216766

<div id="cite-32"></div>
'''[[#citeF-32|[32]]]''' Bilyeu, D. L. and Yu, S.-T. J. and Chen, Y. -Y. and Cambier, J. -L. (2014) "A two-dimensional fourth-order unstrutured-meshed Euler solver based on the CESE method", Volume 257. J. Comput. Phys. 1 981&#8211;999

<div id="cite-33"></div>
'''[[#citeF-33|[33]]]''' S. Jerez and J.V. Romero and M.D. Roselló and F.J. Arnau. (2004) "A semi-implicit space-time CE-SE method to improve mass conservation through tapered ducts in internal combustion engines ", Volume 40. Mathematical and Computer Modelling 9-10 941 - 951

<div id="cite-34"></div>
'''[[#citeF-34|[34]]]''' Shamsul Qamar and Gerald Warnecke. (2006) "Application of space-time CE/SE method to shallow water magnetohydrodynamic equations ", Volume 196. Journal of Computational and Applied Mathematics 1 132 - 149

<div id="cite-35"></div>
'''[[#citeF-35|[35]]]''' Xisheng Luo and Meili Wang and Jiming Yang and Ge Wang. (2007) "The space-time CESE method applied to phase transition of water vapor in compressible flows", Volume 36. Computers & Fluids  7 1247 - 1258

<div id="cite-36"></div>
'''[[#citeF-36|[36]]]''' Yin Chou and Ruey-Jen Yang. (2008) "Application of CESE method to simulate non-Fourier heat conduction in finite medium with pulse surface heating", Volume 51. International Journal of Heat and Mass Transfer 13-14 3525 - 3534

<div id="cite-37"></div>
'''[[#citeF-37|[37]]]''' Young-Il Lim and Sin-Chung Chang and Sten Bay Jrgensen. (2004) "A novel partial differential algebraic equation (PDAE) solver: iterative space-time conservation element/solution element (CE/SE) method ", Volume 28. Computers & Chemical Engineering 8 1309 - 1324

<div id="cite-38"></div>
'''[[#citeF-38|[38]]]''' Sheng-Tao John Yu and Lixiang Yang and Robert L. Lowe and Stephen E. Bechtel. (2010) "Numerical simulation of linear and nonlinear waves in hypoelastic solids by the CESE method ", Volume 47. Wave Motion 3 168 - 182

<div id="cite-39"></div>
'''[[#citeF-39|[39]]]''' Sin-Chung Chang and Xiao-Yen Wang and Chuen-Yen Chow. (1998) "The space-time conservation element and solution element method - a new high-resolution and genuinely multidimensional paradigm for solving conservaton laws". NASA TM 1998-208843

<div id="cite-40"></div>
'''[[#citeF-40|[40]]]''' Lukácová, M. and Morton, K. W. and Warnecke, G. (2000) "Evolution Galerkin Methods for Hyperbolic Systems in Two Space Dimensions", Volume 69. Mathematics of Computation 232 1355-1384

<div id="cite-41"></div>
'''[[#citeF-41|[41]]]''' Qurrat-Ul-Ain and Qamar, Shamsul and Warnecke, Gerald. (2007) "A High-Resolution Space-Time Conservative Method for Non-Linear Hyperbolic Conservation Laws", Volume 4. International Journal of Computational Methods 2 223-247

<div id="cite-42"></div>
'''[[#citeF-42|[42]]]''' Stoker, J. J. (1957) "Water Waves: The Mathematical Theory with Applications", Volume . Interscience Publishers,  Edition

<div id="cite-43"></div>
'''[[#citeF-43|[43]]]''' Delestre, Olivier and Lucas, Carine and Ksinant, Pierre-Antoine and Darboux, Frédéric and Laguerre, Christian and Vo, T.-N.-Tuoi and James, Francois and Cordier, Stéphane. (2013) "SWASHES: a compilation of shallow water analytic solutions for hydraulic and environmental studies", Volume 72. International Journal for Numerical Methods in Fluids 3 269&#8211;300

<div id="cite-44"></div>
'''[[#citeF-44|[44]]]''' Robert J. Fennema and M. Hanif Chaudhry. (1989) "Implicit methods for two-dimensional unsteady free-surface flows", Volume 27. Journal of Hydraulic Research 3 321&#8211;332

<div id="cite-45"></div>
'''[[#citeF-45|[45]]]''' C. Zoppou and S. Roberts. (2000) "Numerical solution of the two-dimensional unsteady dam break", Volume 24. Applied Mathematical Modelling 7 457 - 475

<div id="cite-46"></div>
'''[[#citeF-46|[46]]]''' Francisco Alcrudo and Pilar Garcia-Navarro. (1993) "A high-resolution Godunov-type scheme in finite volumes for 2D shallow-water equations", Volume 16. International Journal for Numerical Methods in Fluids  489&#8211;505