2025-11-18T01:22:13.423549

Existence and numerical approximation of a one-dimensional Boussinesq system with variable coefficients on a finite interval

Grajales, Pizo
In this paper, we investigate the well-posedness of a nonlinear dispersive model with variable coefficients that describes the evolution of surface waves propagating through a one-dimensional shallow water channel of finite length with irregular bottom topography. To complement the theoretical analysis, we utilize the numerical solver developed by the authors in \cite{PizoMunoz} to approximate solutions of the model on a finite spatial interval, considering various parameter values and forms of the variable coefficients in the Boussinesq system under study. Additionally, we present preliminary numerical experiments addressing an inverse problem: the reconstruction of the initial wave elevation and fluid velocity from measurements taken at a final time. This is achieved by formulating an optimization problem in which the initial conditions are estimated as minimizers of a functional that quantifies the discrepancy between the observed final state and the numerical solution evolved from a trial initial state.
academic

Esistenza e approssimazione numerica di un sistema di Boussinesq unidimensionale con coefficienti variabili su un intervallo finito

Informazioni Fondamentali

  • ID Articolo: 2510.10829
  • Titolo: Esistenza e approssimazione numerica di un sistema di Boussinesq unidimensionale con coefficienti variabili su un intervallo finito
  • Autori: Deissy Marcela Pizo, Juan Carlos Muñoz Grajales (Universidad del Valle, Colombia)
  • Classificazione: math.NA, cs.NA, math.AP
  • Data di Pubblicazione: 14 ottobre 2025 (preprint arXiv)
  • Link dell'Articolo: https://arxiv.org/abs/2510.10829

Riassunto

Il presente articolo affronta il problema della buona positura di un modello non lineare dispersivo con coefficienti variabili, che descrive l'evoluzione delle onde di superficie in un canale di acque basse unidimensionale di lunghezza finita con topografia del fondo irregolare. A complemento dell'analisi teorica, gli autori utilizzano un risolutore numerico sviluppato in lavori precedenti (metodo degli elementi finiti) per approssimare le soluzioni del modello su intervalli spaziali limitati, considerando vari valori parametrici e forme funzionali dei coefficienti variabili nel sistema di Boussinesq. Inoltre, l'articolo presenta esperimenti numerici preliminari per affrontare un problema inverso: la ricostruzione dell'altezza d'onda iniziale e della velocità del fluido da dati misurati al tempo finale. Ciò viene realizzato formulando un problema di ottimizzazione in cui le condizioni iniziali sono stimate come soluzioni che minimizzano un funzionale obiettivo, che quantifica la discrepanza tra lo stato finale osservato e la soluzione numerica evoluta da condizioni iniziali di prova.

Contesto di Ricerca e Motivazione

  1. Problemi da Affrontare:
    • Investigare la buona positura matematica (esistenza, unicità, dipendenza continua) di un sistema di Boussinesq unidimensionale con coefficienti variabili
    • Sviluppare metodi numerici efficienti per risolvere tale sistema con coefficienti variabili
    • Esplorare il problema inverso di ricostruzione delle condizioni iniziali da dati osservati al tempo finale
  2. Importanza del Problema:
    • Il sistema di Boussinesq è un modello matematico fondamentale per descrivere la propagazione di onde in acque basse, con importanti applicazioni nell'ingegneria marina e nella protezione costiera
    • I bacini idrici reali presentano spesso topografie del fondo non uniformi, pertanto i modelli con coefficienti variabili rispecchiano meglio la realtà fisica
    • I problemi inversi hanno significativo valore pratico nella previsione delle onde e nel monitoraggio oceanico
  3. Limitazioni dei Metodi Esistenti:
    • La maggior parte degli studi si concentra su sistemi di Boussinesq con coefficienti costanti
    • L'analisi teorica per il caso con coefficienti variabili è relativamente scarsa
    • La ricerca su problemi inversi per tali sistemi dispersivi è ancora in fase iniziale
  4. Motivazione della Ricerca:
    • Colmare il divario nell'analisi teorica dei sistemi di Boussinesq con coefficienti variabili
    • Fornire modelli matematici e strumenti numerici più accurati per applicazioni pratiche in ingegneria marina
    • Aprire nuove direzioni nella ricerca di problemi inversi per sistemi di onde dispersive

Contributi Principali

  1. Stabilimento della teoria di buona positura per il sistema di Boussinesq unidimensionale con coefficienti variabili: dimostrazione dell'esistenza, unicità e dipendenza continua della soluzione dai dati iniziali
  2. Sviluppo di metodi di analisi teorica basati sulla tecnica della funzione di Green: trasformazione del problema in un problema di punto fisso per operatori non lineari
  3. Fornitura di leggi di conservazione dell'energia: dimostrazione della conservazione del funzionale energetico del sistema nel tempo
  4. Implementazione di un risolutore numerico agli elementi finiti efficiente: sviluppo di uno schema numerico stabile combinando la libreria FEniCS
  5. Primo studio del problema inverso per tale classe di sistemi dispersivi: proposizione di un metodo di ricostruzione delle condizioni iniziali basato su ottimizzazione variazionale

Dettagli dei Metodi

Definizione del Problema

Si consideri il seguente sistema di Boussinesq con coefficienti variabili:

N_t + [(1 + \alpha c^2(\xi)N)V]_\xi - \frac{\beta}{6}N_{\xi\xi t} = 0 \\ V_t + (c(\xi)N)_\xi + \frac{1}{2}\alpha(c^2V^2)_\xi - \frac{\beta}{6}V_{\xi\xi t} = 0 \end{cases}$$ dove $(\xi,t) \in [0,L] \times [0,T]$, con condizioni iniziali: $$N(\xi,0) = N_0(\xi), \quad V(\xi,0) = V_0(\xi)$$ e condizioni al contorno di Dirichlet: $$N(0,t) = N(L,t) = 0, \quad V(0,t) = V(L,t) = 0$$ ### Metodi di Analisi Teorica **1. Tecnica della Funzione di Green** Il sistema viene riscritto come: $$\left(I - \frac{\beta}{6}\partial_\xi^2\right)N_t = -\partial_\xi[(1 + \alpha c^2(\xi)N)V]$$ $$\left(I - \frac{\beta}{6}\partial_\xi^2\right)V_t = -\partial_\xi\left[c(\xi)N + \frac{1}{2}\alpha c^2V^2\right]$$ Utilizzando la funzione di Green dell'operatore $P = I - \frac{\beta}{6}\partial_\xi^2$: $$G(\xi,s) = \frac{1}{2\sqrt{\beta/6}} \frac{\cosh\left(\frac{L-|\xi-s|}{\sqrt{\beta/6}}\right) - \cosh\left(\frac{L-(\xi+s)}{\sqrt{\beta/6}}\right)}{\sinh\left(\frac{L}{\sqrt{\beta/6}}\right)}$$ **2. Teorema del Punto Fisso** Il problema viene trasformato in un sistema di equazioni integrali, definendo un operatore non lineare $A$, e l'esistenza e l'unicità della soluzione vengono provate mediante il teorema del punto fisso di Banach. ### Metodi Numerici **Discretizzazione agli Elementi Finiti** Si utilizza il metodo di Galerkin, con forma debole: $$\partial_t\left(\langle N,v_1\rangle_{L^2} + \frac{\beta}{6}\langle N_\xi, v_{1\xi}\rangle_{L^2}\right) = \langle[(1+\alpha c^2N)V], v_{1\xi}\rangle_{L^2}$$ **Discretizzazione Temporale** Si impiega il metodo $\theta$ (con $\theta = 1/2$) per l'integrazione temporale, garantendo la stabilità numerica. ### Metodo per il Problema Inverso Si definisce il funzionale obiettivo: $$J(N,V) := \frac{1}{2}\int_0^L \left(|N(\xi,T;N_0,V_0) - N_T(\xi)|^2 + |V(\xi,T;N_0,V_0) - V_T(\xi)|^2\right)d\xi$$ Si utilizza l'algoritmo L-BFGS-B per minimizzare tale funzionale, combinato con la libreria Dolfin-Adjoint per l'implementazione della differenziazione automatica. ## Configurazione degli Esperimenti ### Parametri Numerici - **Dominio spaziale**: $[-20, 40]$ - **Numero di punti della griglia**: 3000 punti spaziali - **Risoluzione spaziale**: $\Delta\xi \approx 0.02$ - **Passo temporale**: $\Delta t = 8/3000 \approx 2.7 \times 10^{-3}$ - **Parametri del modello**: $\alpha = \beta = 0.1$ ### Condizioni Iniziali Si utilizza un profilo gaussiano localizzato: $$N(\xi,0) = V(\xi,0) = e^{-(\xi-18)^2}$$ ### Funzioni con Coefficienti Variabili **Caso 1**: Oscillazione liscia con perturbazione gaussiana $$c(\xi) = 1 + 0.3\sin\left(\frac{\pi}{5}\xi\right) + 0.6e^{-(\xi-8)^2}$$ **Caso 2**: Funzione costante a tratti (simulazione di mezzo stratificato) $$c(\xi) = \begin{cases} 0.8, & \xi < 20 \\ 1.5, & 20 < \xi \leq 22 \\ 2.0, & 22 < \xi \leq 27 \\ \vdots \end{cases}$$ ## Risultati degli Esperimenti ### Risultati Teorici **Teorema 3.1 (Buona Positura Locale)**: Esiste $T_0 > 0$ tale che il sistema ammette un'unica soluzione $(N,V)$ in $L^2_{T_0} \times L^2_{T_0}$, soddisfacendo la stima: $$\|(N,V)\|_{L^2_{T_0} \times L^2_{T_0}} \leq \|(N_0,V_0)\|_{H^1_0 \times H^1_0} T_0^{1/2} \exp\left(\frac{L^{1/2}}{(\beta/6)^{1/2}}DT_0\right)$$ **Teorema 3.2 (Conservazione dell'Energia)**: $$\|V\|_{L^2}^2 + \alpha\|c|N|^{1/2}V\|_{L^2}^2 + \||c|^{1/2}N\|_{L^2}^2 = \text{costante}$$ ### Risultati degli Esperimenti Numerici 1. **Risoluzione del Problema Diretto**: - Simulazione riuscita della propagazione delle onde in mezzo con coefficienti variabili - Osservazione di variazioni della velocità d'onda, modulazione dell'ampiezza ed effetti dispersivi dovuti alle variazioni della topografia del fondo - Il metodo numerico dimostra buona stabilità e accuratezza 2. **Ricostruzione del Problema Inverso**: - L'algoritmo L-BFGS-B migliora significativamente la ricostruzione delle condizioni iniziali in 4 iterazioni - L'altezza d'onda e la velocità del fluido ricostruite convergono progressivamente verso il profilo target - Viene provata l'efficacia della strategia di assimilazione variazionale ### Verifica Numerica - La legge di conservazione dell'energia è ben preservata nei calcoli numerici - La coerenza dei risultati con diverse risoluzioni di griglia verifica la convergenza del metodo numerico ## Lavori Correlati 1. **Teoria dei Sistemi di Boussinesq**: - Bona & Chen (1998, 2002, 2004): stabilimento della teoria fondamentale per il caso con coefficienti costanti - Il presente articolo estende il risultato al caso con coefficienti variabili, colmando un divario teorico 2. **Metodi Numerici**: - Applicazione del metodo degli elementi finiti alle equazioni di onde dispersive - Il presente articolo implementa un risolutore efficiente combinando la libreria FEniCS 3. **Problemi Inversi**: - Metodi tradizionali per problemi inversi di equazioni di onde - Il presente articolo applica per la prima volta l'ottimizzazione variazionale ai problemi inversi del sistema di Boussinesq ## Conclusioni e Discussione ### Conclusioni Principali 1. Stabilimento riuscito delle fondamenta matematiche della teoria per il sistema di Boussinesq unidimensionale con coefficienti variabili 2. Sviluppo di metodi numerici efficienti per la risoluzione 3. Apertura di nuove direzioni nella ricerca di problemi inversi per sistemi dispersivi 4. Fornitura di strumenti teorici e computazionali per applicazioni pratiche in ingegneria marina ### Limitazioni 1. **Aspetto Teorico**: - È provata solo la buona positura locale; l'esistenza globale rimane da investigare - Manca l'analisi teorica dell'unicità e della stabilità del problema inverso 2. **Aspetto Numerico**: - Gli esperimenti sul problema inverso considerano solo il caso ideale senza rumore - Manca l'analisi della sensibilità al rumore di osservazione 3. **Aspetto Applicativo**: - Limitato al caso unidimensionale; le applicazioni pratiche richiedono estensione al caso bidimensionale - La scelta delle funzioni con coefficienti variabili è relativamente semplice ### Direzioni Future 1. **Sviluppo Teorico**: - Teoria della buona positura globale - Analisi dell'unicità e della stabilità del problema inverso - Teoria della regolarizzazione 2. **Miglioramenti Numerici**: - Metodi robusti di ricostruzione con dati rumorosi - Estensione a dimensioni superiori - Integrazione di tecniche di apprendimento automatico 3. **Applicazioni Pratiche**: - Validazione con dati oceanici reali - Integrazione con sistemi di osservazione oceanica ## Valutazione Approfondita ### Punti di Forza 1. **Contributo Teorico Significativo**: primo studio sistematico della buona positura del sistema di Boussinesq con coefficienti variabili 2. **Innovazione Metodologica**: applicazione ingegnosa della tecnica della funzione di Green per affrontare problemi con coefficienti variabili 3. **Implementazione Numerica Completa**: fornitura di un quadro completo di risoluzione numerica 4. **Novità del Problema**: apertura di un nuovo campo nella ricerca di problemi inversi per sistemi dispersivi 5. **Chiarezza della Presentazione**: derivazioni matematiche rigorose e esposizione chiara ### Insufficienze 1. **Profondità Teorica Limitata**: stabilimento solo della teoria locale, mancanza di analisi globale 2. **Verifica Sperimentale Insufficiente**: esperimenti sul problema inverso relativamente semplici, mancanza di verifica in scenari complessi 3. **Praticità da Verificare**: necessità di ulteriore verifica con dati reali 4. **Mancanza di Analisi della Complessità Computazionale**: non è fornita l'analisi della complessità computazionale dell'algoritmo ### Impatto 1. **Valore Accademico**: posa le fondamenta per lo sviluppo della teoria dei sistemi dispersivi con coefficienti variabili 2. **Prospettive Applicative**: potenziale valore applicativo in ingegneria marina, previsione delle onde e altri campi 3. **Contributo Metodologico**: la combinazione del metodo della funzione di Green e dell'ottimizzazione variazionale fornisce nuove prospettive per problemi correlati ### Scenari di Applicabilità 1. Modellazione della propagazione di onde in acque basse 2. Analisi delle onde nell'ingegneria costiera 3. Problemi di inversione nel monitoraggio ambientale marino 4. Assimilazione dati nell'oceanografia numerica ## Bibliografia L'articolo cita 20 importanti riferimenti bibliografici, coprendo la teoria dei sistemi di Boussinesq, metodi numerici, analisi funzionale e altri campi correlati, fornendo una solida base teorica per la ricerca.