Interface problems of mixed spatial order

Resumo Técnico: Problemas de Interface de Ordem Mista via o Método da Transformada Unificada

Enunciado do Problema
Este trabalho aborda o problema de valor inicial e de contorno (IBVP) para equações diferenciais parciais (EDPs) de evolução linear definidas em um domínio espacial dividido em duas semirretas, $(-\infty, 0)$ e $(0, \infty)$, separadas por uma interface em $x=0$. Diferentemente de estudos anteriores, onde as equações governantes em cada lado da interface eram idênticas ou diferiam apenas por coeficientes escalares, este artigo investiga problemas de interface de "ordem mista". Especificamente, considera a interação entre equações de evolução linear de ordens espaciais e comportamentos físicos fundamentalmente diferentes (dispersivos vs. dissipativos).

Os exemplos prototípicos analisados envolvem:

1. Semirreta Esquerda ($x < 0$): a equação de Schrödinger linear (lS), a equação do calor (h), ou a equação de Schrödinger biarmônica (biS).
2. Semirreta Direita ($x > 0$): a equação de Airy ou a equação de Korteweg-de Vries linearizada (lKdV).

Estes sistemas são acoplados via condições de interface em $x=0$. O artigo considera dois conjuntos de condições:

• Condição (IfC.1): continuidade da solução e de sua primeira derivada espacial ($C^1$ continuidade).
• Condição (IfC.2): continuidade da solução e de suas duas primeiras derivadas espaciais ($C^2$ continuidade), necessária para a equação de Schrödinger biarmônica de quarta ordem.

Os valores de interface (dados de contorno em $x=0$) não são prescritos a priori; são funções desconhecidas determinadas pelo acoplamento das duas EDPs distintas.

Metodologia
Os autores empregam o Método da Transformada Unificada (UTM), também conhecido como método de Fokas, que serve como uma generalização da transformada de Fourier para IBVPs em domínios com fronteiras. O procedimento de solução envolve as seguintes etapas:

1. Relações Globais: Para cada problema de semirreta, os autores derivam "relações globais" (identidades espectrais) que ligam os dados iniciais, a solução e os valores de contorno desconhecidos no plano espectral complexo.
2. Mudança de Variáveis: Para lidar com as ordens mistas, os autores introduzem funções específicas do tipo raiz ($\sigma, \rho, \tau, \gamma$) e um biholomorfismo $\nu$ (para o caso lKdV), que mapeiam o parâmetro espectral $\lambda$ para o plano complexo $k$ de modo a alinhar as relações de dispersão das diferentes equações (por exemplo, relacionando $k^2, k^3$ e $k^4$).
3. Construção do Sistema Linear: Ao aplicar estas mudanças de variáveis às relações globais, os autores constroem um sistema linear (tipicamente $3\times3$ ou $4\times4$) para as transformadas dos valores de interface desconhecidos.
4. Representação da Solução: Resolver este sistema linear permite a eliminação dos valores de interface desconhecidos. Os autores então substituem estas expressões de volta nas representações integrais das soluções para cada semirreta.
5. Deformação de Contorno e Assintótica: Uma contribuição técnica crítica é a justificativa rigorosa para a deformação dos contornos de integração (por exemplo, da reta real para contornos como $\tilde{R}$ no semiplano inferior) para contornar singularidades (especificamente em $\lambda=0$) e para eliminar integrais envolvendo os termos dependentes do tempo da solução ($u(x,t), v(x,t)$). Isto é alcançado através de argumentos específicos utilizando o comportamento assintótico das funções raiz e dos núcleos exponenciais, em vez do padrão do Lema de Jordan, devido ao descompasso nos termos de ordem principal.

Contribuições Principais e Resultados
O artigo fornece fórmulas de solução explícitas para quatro problemas de interface distintos, apresentados nas Proposições 1 através 4:

• Proposição 1: Resolve a interface entre a equação de Schrödinger linear e a equação de Airy sob condições de interface $C^1$.
• Proposição 2: Resolve a interface entre a equação de Schrödinger linear e a equação de KdV linearizada (incluindo um termo de ordem inferior $av_x$) sob condições $C^1$. Isto requer a introdução do biholomorfismo $\nu$ para lidar com a relação de dispersão não homogênea.
• Proposição 3: Resolve a interface entre a equação do calor (dissipativa) e a equação de Airy (dispersiva) sob condições $C^1$.
• Proposição 4: Resolve a interface entre a equação de Schrödinger biarmônica (dispersiva de quarta ordem) e a equação de Airy sob condições $C^2$. Este caso envolve um sistema linear $4\times4$ devido à ordem superior da EDP.

As fórmulas explícitas são expressas como integrais de contorno envolvendo as transformadas dos dados iniciais e funções de núcleo específicas derivadas dos mapeamentos do tipo raiz.

Avaliação Numérica
Para demonstrar a utilidade das fórmulas derivadas, os autores realizam avaliações numéricas usando quadratura de Filon. Este método é escolhido para lidar com a natureza altamente oscilatória dos integrandos. O artigo apresenta gráficos das soluções para várias configurações de dados iniciais:

• Schrödinger-Airy: Mostra a transição de ondas dispersivas do regime unidirecional de Airy para o regime bidirecional de Schrödinger, e vice-versa, destacando os padrões de "efeito roleta" e efeitos de arrasto causados pela interface.
• Biharmonic Schrödinger-Airy: Ilustra interações de ondas semelhantes, mas com a dinâmica governada por uma relação de dispersão quártica, levando à emergência mais rápida de padrões complexos.
• Heat-Airy: Demonstra a interação entre ondas difusivas e dispersivas, mostrando decaimento exponencial no lado do calor e dispersão no lado de Airy.

Significância
O artigo afirma sua significância ao estender o Método da Transformada Unificada para uma nova classe de problemas onde as equações governantes de cada lado de uma interface são de tipos fundamentalmente diferentes (diferentes ordens espaciais e naturezas físicas). Aplicações anteriores do UTM a problemas de interface foram amplamente restritas a casos onde as equações eram idênticas ou diferiam apenas por coeficientes escalares.

Ao resolver com sucesso estes problemas de ordem mista, os autores demonstram a robustez do UTM em lidar com estruturas espectrais complexas decorrentes de relações de dispersão díspares. Este trabalho fornece o primeiro quadro de solução analítica para tais interações, que são relevantes para modelar fenômenos físicos envolvendo meios heterogêneos, como o crescimento de tumores em tecido cerebral (conforme citado na introdução em relação aos gliomas) ou materiais compostos com propriedades físicas variadas. As fórmulas explícitas derivadas servem como base para análises futuras, incluindo potenciais extensões para problemas não lineares via técnicas de contração.