June  2019, 6(1): 131-145. doi: 10.3934/jcd.2019006

Numerical efficacy study of data assimilation for the 2D magnetohydrodynamic equations

1. 

University of Maryland, Baltimore County, 1000 Hilltop Circle, Baltimore, MD 21055, USA

2. 

Johns Hopkins University Applied Physics Laboratory, 11100 Johns Hopkins Road, Laurel, MD 20723-6099, USA

3. 

Department of Mathematics, Indiana University, Bloomington, IN 47405, USA

* Corresponding author: Michael Jolly

Published  August 2019

We study the computational efficiency of several nudging data assimilation algorithms for the 2D magnetohydrodynamic equations, using varying amounts and types of data. We find that the algorithms work with much less resolution in the data than required by the rigorous estimates in [7]. We also test other abridged nudging algorithms to which the analytic techniques in [7] do not seem to apply. These latter tests indicate, in particular, that velocity data alone is sufficient for synchronization with a chaotic reference solution, while magnetic data alone is not. We demonstrate that a new nonlinear nudging algorithm, which is adaptive in both time and space, synchronizes at a super exponential rate.

Citation: Joshua Hudson, Michael Jolly. Numerical efficacy study of data assimilation for the 2D magnetohydrodynamic equations. Journal of Computational Dynamics, 2019, 6 (1) : 131-145. doi: 10.3934/jcd.2019006
References:
[1]

D. A. AlbanezH. J. Nussenzveig Lopes and E. S. Titi, Continuous data assimilation for the three-dimensional Navier–Stokes-$\alpha$ model, Asymptotic Anal., 97 (2016), 139-164.  doi: 10.3233/ASY-151351.  Google Scholar

[2]

M. U. AltafE. S. TitiT. GebraelO. M. KnioL. ZhaoM. F. McCabe and I. Hoteit, Downscaling the 2D Bénard convection equations using continuous data assimilation, Comput. Geosci., 21 (2017), 393-410.  doi: 10.1007/s10596-017-9619-2.  Google Scholar

[3]

A. AzouaniE. Olson and E. S. Titi, Continuous data assimilation using general interpolant observables, J. Nonlinear Sci., 24 (2014), 277-304.  doi: 10.1007/s00332-013-9189-y.  Google Scholar

[4]

A. Azouani and E. S. Titi, Feedback control of nonlinear dissipative systems by finite determining parameters—a reaction-diffusion paradigm, Evol. Equ. Control Theory, 3 (2014), 579-594.  doi: 10.3934/eect.2014.3.579.  Google Scholar

[5]

H. BessaihE. Olson and E. S. Titi, Continuous data assimilation with stochastically noisy data, Nonlinearity, 28 (2015), 729-753.  doi: 10.1088/0951-7715/28/3/729.  Google Scholar

[6]

A. BiswasC. FoiasC. F. Mondaini and E. S. Titi, Downscaling data assimilation algorithm with applications to statistical solutions of the Navier-Stokes equations, Ann. Inst. H. Poincaré Anal. Non Linéaire, 36 (2019), 295-326.  doi: 10.1016/j.anihpc.2018.05.004.  Google Scholar

[7]

A. BiswasJ. HudsonA. Larios and Y. Pei, Continuous data assimilation for the 2D magnetohydrodynamic equations using one component of the velocity and magnetic fields, Asymptot. Anal., 108 (2018), 1-43.   Google Scholar

[8]

K. J. Burns, G. M. Vasil, J. S. Oishi, D. Lecoanet and B. Brown, Dedalus: Flexible framework for spectrally solving differential equations, Astrophysics Source Code Library, 2016. Google Scholar

[9]

E. CelikE. Olson and E. S. Titi, Spectral Filtering of Interpolant Observables for a Discrete-in-Time Downscaling Data Assimilation Algorithm, SIAM J. Appl. Dyn. Syst., 18 (2019), 1118-1142.  doi: 10.1137/18M1218480.  Google Scholar

[10]

J. CharneyM. Halem and R. Jastrow, Use of incomplete historical data to infer the present state of the atmosphere, Journal of the Atmospheric Sciences, 26 (1969), 1160-1163.  doi: 10.1175/1520-0469(1969)026<1160:UOIHDT>2.0.CO;2.  Google Scholar

[11]

R. DascaliucC. Foias and M. Jolly, Relations between energy and enstrophy on the global attractor of the 2-d Navier-Stokes equations, J. Dynamics Differential Equations, 17 (2005), 643-736.  doi: 10.1007/s10884-005-8269-6.  Google Scholar

[12]

D. Drake, C. Kletzing, F. Skiff, G. Howes and S. Vincena, Design and use of an Elsässer probe for analysis of Alfén wave field according to wave direction, The Review of Scientific Instruments, 82 (2011), 103505. Google Scholar

[13]

G. Evensen, Data Assimilation, 2nd edition, Springer-Verlag, Berlin, 2009, The ensemble Kalman filter. doi: 10.1007/978-3-642-03711-5.  Google Scholar

[14]

A. FarhatH. JohnstonM. Jolly and E. S. Titi, Assimilation of nearly turbulent Rayleigh-Bénard flow through vorticity or local circulation measurements: a computational study, J. Sci. Comput., 77 (2018), 1519-1533.  doi: 10.1007/s10915-018-0686-x.  Google Scholar

[15]

A. FarhatM. S. Jolly and E. S. Titi, Continuous data assimilation for the 2D Bénard convection through velocity measurements alone, Phys. D, 303 (2015), 59-66.  doi: 10.1016/j.physd.2015.03.011.  Google Scholar

[16]

A. FarhatE. Lunasin and E. S. Titi, Abridged continuous data assimilation for the 2D Navier–Stokes equations utilizing measurements of only one component of the velocity field, J. Math. Fluid Mech., 18 (2016), 1-23.  doi: 10.1007/s00021-015-0225-6.  Google Scholar

[17]

A. FarhatE. Lunasin and E. S. Titi, Data assimilation algorithm for 3D Bénard convection in porous media employing only temperature measurements, Journal of Mathematical Analysis and Applications, 438 (2016), 492-506.  doi: 10.1016/j.jmaa.2016.01.072.  Google Scholar

[18]

A. Farhat, E. Lunasin and E. S. Titi, On the Charney conjecture of data assimilation employing temperature measurements alone: The paradigm of 3d planetary geostrophic model, arXiv: 1608.04770. Google Scholar

[19]

A. FarhatE. Lunasin and E. S. Titi, Continuous data assimilation for a 2D Bénard convection system through horizontal velocity measurements alone, J. Nonlinear Sci., 27 (2017), 1065-1087.  doi: 10.1007/s00332-017-9360-y.  Google Scholar

[20]

M. GeshoE. Olson and E. S. Titi, A computational study of a data assimilation algorithm for the two-dimensional Navier-Stokes equations, Commun. Comput. Phys., 19 (2016), 1094-1110.  doi: 10.4208/cicp.060515.161115a.  Google Scholar

[21]

M. GhilM. Halem and R. Atlas, Time-continuous assimilation of remote-sounding data and its effect on weather forecasting, Mon. Wea. Rev., 107 (1978), 140-171.   Google Scholar

[22]

K. HaydenE. Olson and E. S. Titi, Discrete data assimilation in the Lorenz and 2D Navier-Stokes equations, Phys. D, 240 (2011), 1416-1425.  doi: 10.1016/j.physd.2011.04.021.  Google Scholar

[23]

A. JacksonA. R. T. Jonkers and M. R. Walker, Four centuries of geomagnetic secular variation from historical records, Phil. Trans. R. Soc. Lond., 358 (2000), 957-990.  doi: 10.1098/rsta.2000.0569.  Google Scholar

[24]

M. S. JollyV. R. Martinez and E. S. Titi, A data assimilation algorithm for the subcritical surface quasi-geostrophic equation, Adv. Nonlinear Stud., 17 (2017), 167-192.   Google Scholar

[25]

M. S. JollyT. Sadigov and E. S. Titi, Determining form and data assimilation algorithm for weakly damped and driven Korteweg–de Vries equation—Fourier modes case, Nonlinear Anal. Real World Appl., 36 (2017), 287-317.  doi: 10.1016/j.nonrwa.2017.01.010.  Google Scholar

[26]

A. R. T. Jonkers, A. Jackson and A. Murray, Four centuries of geomagnetic data from historical records, Rev. Geophys., 41. Google Scholar

[27]

A. Larios and Y. Pei, Nonlinear continuous data assimilation, arXiv: 1703.03546. Google Scholar

[28]

D. Leunberger, An introduction to observers, IEEE Trans. Automat. Control, 16 (1971), 596-602.  doi: 10.1109/TAC.1971.1099826.  Google Scholar

[29]

P. A. MarkowichE. S. Titi and S. Trabelsi, Continuous data assimilation for the three-dimensional Brinkman–Forchheimer-extended Darcy model, Nonlinearity, 29 (2016), 1292-1328.  doi: 10.1088/0951-7715/29/4/1292.  Google Scholar

[30]

E. Olson and E. S. Titi, Determining modes for continuous data assimilation in 2D turbulence, J. Statist. Phys., 113 (2003), 799–840, Progress in statistical hydrodynamics (Santa Fe, NM, 2002). doi: 10.1023/A:1027312703252.  Google Scholar

[31] S. Reich and C. Cotter, Probabilistic Forecasting and Bayesian Data Assimilation, Cambridge University Press, New York, 2015.  doi: 10.1017/CBO9781107706804.  Google Scholar
[32]

F. E. Thau, Observing the state of non-linear dynamic systems, Int. J. Control, 17 (1973), 471-479.  doi: 10.1080/00207177308932395.  Google Scholar

show all references

References:
[1]

D. A. AlbanezH. J. Nussenzveig Lopes and E. S. Titi, Continuous data assimilation for the three-dimensional Navier–Stokes-$\alpha$ model, Asymptotic Anal., 97 (2016), 139-164.  doi: 10.3233/ASY-151351.  Google Scholar

[2]

M. U. AltafE. S. TitiT. GebraelO. M. KnioL. ZhaoM. F. McCabe and I. Hoteit, Downscaling the 2D Bénard convection equations using continuous data assimilation, Comput. Geosci., 21 (2017), 393-410.  doi: 10.1007/s10596-017-9619-2.  Google Scholar

[3]

A. AzouaniE. Olson and E. S. Titi, Continuous data assimilation using general interpolant observables, J. Nonlinear Sci., 24 (2014), 277-304.  doi: 10.1007/s00332-013-9189-y.  Google Scholar

[4]

A. Azouani and E. S. Titi, Feedback control of nonlinear dissipative systems by finite determining parameters—a reaction-diffusion paradigm, Evol. Equ. Control Theory, 3 (2014), 579-594.  doi: 10.3934/eect.2014.3.579.  Google Scholar

[5]

H. BessaihE. Olson and E. S. Titi, Continuous data assimilation with stochastically noisy data, Nonlinearity, 28 (2015), 729-753.  doi: 10.1088/0951-7715/28/3/729.  Google Scholar

[6]

A. BiswasC. FoiasC. F. Mondaini and E. S. Titi, Downscaling data assimilation algorithm with applications to statistical solutions of the Navier-Stokes equations, Ann. Inst. H. Poincaré Anal. Non Linéaire, 36 (2019), 295-326.  doi: 10.1016/j.anihpc.2018.05.004.  Google Scholar

[7]

A. BiswasJ. HudsonA. Larios and Y. Pei, Continuous data assimilation for the 2D magnetohydrodynamic equations using one component of the velocity and magnetic fields, Asymptot. Anal., 108 (2018), 1-43.   Google Scholar

[8]

K. J. Burns, G. M. Vasil, J. S. Oishi, D. Lecoanet and B. Brown, Dedalus: Flexible framework for spectrally solving differential equations, Astrophysics Source Code Library, 2016. Google Scholar

[9]

E. CelikE. Olson and E. S. Titi, Spectral Filtering of Interpolant Observables for a Discrete-in-Time Downscaling Data Assimilation Algorithm, SIAM J. Appl. Dyn. Syst., 18 (2019), 1118-1142.  doi: 10.1137/18M1218480.  Google Scholar

[10]

J. CharneyM. Halem and R. Jastrow, Use of incomplete historical data to infer the present state of the atmosphere, Journal of the Atmospheric Sciences, 26 (1969), 1160-1163.  doi: 10.1175/1520-0469(1969)026<1160:UOIHDT>2.0.CO;2.  Google Scholar

[11]

R. DascaliucC. Foias and M. Jolly, Relations between energy and enstrophy on the global attractor of the 2-d Navier-Stokes equations, J. Dynamics Differential Equations, 17 (2005), 643-736.  doi: 10.1007/s10884-005-8269-6.  Google Scholar

[12]

D. Drake, C. Kletzing, F. Skiff, G. Howes and S. Vincena, Design and use of an Elsässer probe for analysis of Alfén wave field according to wave direction, The Review of Scientific Instruments, 82 (2011), 103505. Google Scholar

[13]

G. Evensen, Data Assimilation, 2nd edition, Springer-Verlag, Berlin, 2009, The ensemble Kalman filter. doi: 10.1007/978-3-642-03711-5.  Google Scholar

[14]

A. FarhatH. JohnstonM. Jolly and E. S. Titi, Assimilation of nearly turbulent Rayleigh-Bénard flow through vorticity or local circulation measurements: a computational study, J. Sci. Comput., 77 (2018), 1519-1533.  doi: 10.1007/s10915-018-0686-x.  Google Scholar

[15]

A. FarhatM. S. Jolly and E. S. Titi, Continuous data assimilation for the 2D Bénard convection through velocity measurements alone, Phys. D, 303 (2015), 59-66.  doi: 10.1016/j.physd.2015.03.011.  Google Scholar

[16]

A. FarhatE. Lunasin and E. S. Titi, Abridged continuous data assimilation for the 2D Navier–Stokes equations utilizing measurements of only one component of the velocity field, J. Math. Fluid Mech., 18 (2016), 1-23.  doi: 10.1007/s00021-015-0225-6.  Google Scholar

[17]

A. FarhatE. Lunasin and E. S. Titi, Data assimilation algorithm for 3D Bénard convection in porous media employing only temperature measurements, Journal of Mathematical Analysis and Applications, 438 (2016), 492-506.  doi: 10.1016/j.jmaa.2016.01.072.  Google Scholar

[18]

A. Farhat, E. Lunasin and E. S. Titi, On the Charney conjecture of data assimilation employing temperature measurements alone: The paradigm of 3d planetary geostrophic model, arXiv: 1608.04770. Google Scholar

[19]

A. FarhatE. Lunasin and E. S. Titi, Continuous data assimilation for a 2D Bénard convection system through horizontal velocity measurements alone, J. Nonlinear Sci., 27 (2017), 1065-1087.  doi: 10.1007/s00332-017-9360-y.  Google Scholar

[20]

M. GeshoE. Olson and E. S. Titi, A computational study of a data assimilation algorithm for the two-dimensional Navier-Stokes equations, Commun. Comput. Phys., 19 (2016), 1094-1110.  doi: 10.4208/cicp.060515.161115a.  Google Scholar

[21]

M. GhilM. Halem and R. Atlas, Time-continuous assimilation of remote-sounding data and its effect on weather forecasting, Mon. Wea. Rev., 107 (1978), 140-171.   Google Scholar

[22]

K. HaydenE. Olson and E. S. Titi, Discrete data assimilation in the Lorenz and 2D Navier-Stokes equations, Phys. D, 240 (2011), 1416-1425.  doi: 10.1016/j.physd.2011.04.021.  Google Scholar

[23]

A. JacksonA. R. T. Jonkers and M. R. Walker, Four centuries of geomagnetic secular variation from historical records, Phil. Trans. R. Soc. Lond., 358 (2000), 957-990.  doi: 10.1098/rsta.2000.0569.  Google Scholar

[24]

M. S. JollyV. R. Martinez and E. S. Titi, A data assimilation algorithm for the subcritical surface quasi-geostrophic equation, Adv. Nonlinear Stud., 17 (2017), 167-192.   Google Scholar

[25]

M. S. JollyT. Sadigov and E. S. Titi, Determining form and data assimilation algorithm for weakly damped and driven Korteweg–de Vries equation—Fourier modes case, Nonlinear Anal. Real World Appl., 36 (2017), 287-317.  doi: 10.1016/j.nonrwa.2017.01.010.  Google Scholar

[26]

A. R. T. Jonkers, A. Jackson and A. Murray, Four centuries of geomagnetic data from historical records, Rev. Geophys., 41. Google Scholar

[27]

A. Larios and Y. Pei, Nonlinear continuous data assimilation, arXiv: 1703.03546. Google Scholar

[28]

D. Leunberger, An introduction to observers, IEEE Trans. Automat. Control, 16 (1971), 596-602.  doi: 10.1109/TAC.1971.1099826.  Google Scholar

[29]

P. A. MarkowichE. S. Titi and S. Trabelsi, Continuous data assimilation for the three-dimensional Brinkman–Forchheimer-extended Darcy model, Nonlinearity, 29 (2016), 1292-1328.  doi: 10.1088/0951-7715/29/4/1292.  Google Scholar

[30]

E. Olson and E. S. Titi, Determining modes for continuous data assimilation in 2D turbulence, J. Statist. Phys., 113 (2003), 799–840, Progress in statistical hydrodynamics (Santa Fe, NM, 2002). doi: 10.1023/A:1027312703252.  Google Scholar

[31] S. Reich and C. Cotter, Probabilistic Forecasting and Bayesian Data Assimilation, Cambridge University Press, New York, 2015.  doi: 10.1017/CBO9781107706804.  Google Scholar
[32]

F. E. Thau, Observing the state of non-linear dynamic systems, Int. J. Control, 17 (1973), 471-479.  doi: 10.1080/00207177308932395.  Google Scholar

Figure 1.  Properties of the reference solution. In (a) and (b) the evolution of the $ L^2 $ norm of the reference solution computed with $ 256^2 $ resolution is shown for $ t\in[0,100] $ and $ t\in[10,792.92] $. A comparison with $ 512^2 $ resolution is in (a). Plots of the $ H^1 $ semi-norm vs the $ L^2 $ norm of the solutions for $ t\in[10,90] $ with $ 256^2 $ and $ 512^2 $ resolutions are in (c) and (d) respectively
Figure 2.  Contour lines of the curl of the computed reference solution at time $ t = 729.92 $
Figure 3.  Dependence of the error (14) on $ \mu $ and $ N $. The solutions were computed over the time interval $ [729.9,734.9] $, and the error is at time $ t = 734.9 $
Figure 4.  Convergence results for Algorithms 2.2, 2.3, and 2.4, with damping $ \mu = 20 $
Figure 5.  Convergence results for Algorithms 4.1-4.4 with $ \mu = 20 $
Figure 6.  The evolution of the $ L^2 $ error is shown for solutions of nonlinear modifications of Algorithm 2.2. Each simulation was performed with $ N = 32 $
Figure 7.  The evolution of the $ L^2 $ error is shown for solutions of Algorithm 2.2 when subject to simulated measurement noise
[1]

Yuan Pei. Continuous data assimilation for the 3D primitive equations of the ocean. Communications on Pure & Applied Analysis, 2019, 18 (2) : 643-661. doi: 10.3934/cpaa.2019032

[2]

Alexandre J. Chorin, Fei Lu, Robert N. Miller, Matthias Morzfeld, Xuemin Tu. Sampling, feasibility, and priors in data assimilation. Discrete & Continuous Dynamical Systems - A, 2016, 36 (8) : 4227-4246. doi: 10.3934/dcds.2016.36.4227

[3]

Débora A. F. Albanez, Maicon J. Benvenutti. Continuous data assimilation algorithm for simplified Bardina model. Evolution Equations & Control Theory, 2018, 7 (1) : 33-52. doi: 10.3934/eect.2018002

[4]

Issam S. Strub, Julie Percelay, Olli-Pekka Tossavainen, Alexandre M. Bayen. Comparison of two data assimilation algorithms for shallow water flows. Networks & Heterogeneous Media, 2009, 4 (2) : 409-430. doi: 10.3934/nhm.2009.4.409

[5]

Juan Carlos De los Reyes, Estefanía Loayza-Romero. Total generalized variation regularization in data assimilation for Burgers' equation. Inverse Problems & Imaging, 2019, 13 (4) : 755-786. doi: 10.3934/ipi.2019035

[6]

Jishan Fan, Shuxiang Huang, Fucai Li. Global strong solutions to the planar compressible magnetohydrodynamic equations with large initial data and vacuum. Kinetic & Related Models, 2017, 10 (4) : 1035-1053. doi: 10.3934/krm.2017041

[7]

Yaobin Ou, Pan Shi. Global classical solutions to the free boundary problem of planar full magnetohydrodynamic equations with large initial data. Discrete & Continuous Dynamical Systems - B, 2017, 22 (2) : 537-567. doi: 10.3934/dcdsb.2017026

[8]

Yanmin Mu. Convergence of the compressible isentropic magnetohydrodynamic equations to the incompressible magnetohydrodynamic equations in critical spaces. Kinetic & Related Models, 2014, 7 (4) : 739-753. doi: 10.3934/krm.2014.7.739

[9]

Xiaoli Li, Dehua Wang. Global solutions to the incompressible magnetohydrodynamic equations. Communications on Pure & Applied Analysis, 2012, 11 (2) : 763-783. doi: 10.3934/cpaa.2012.11.763

[10]

Jishan Fan, Tohru Ozawa. An approximation model for the density-dependent magnetohydrodynamic equations. Conference Publications, 2013, 2013 (special) : 207-216. doi: 10.3934/proc.2013.2013.207

[11]

Elder Jesús Villamizar-Roa, Henry Lamos-Díaz, Gilberto Arenas-Díaz. Very weak solutions for the magnetohydrodynamic type equations. Discrete & Continuous Dynamical Systems - B, 2008, 10 (4) : 957-972. doi: 10.3934/dcdsb.2008.10.957

[12]

Tong Tang, Hongjun Gao. Local strong solutions to the compressible viscous magnetohydrodynamic equations. Discrete & Continuous Dynamical Systems - B, 2016, 21 (5) : 1617-1633. doi: 10.3934/dcdsb.2016014

[13]

Zhong Tan, Leilei Tong. Asymptotic stability of stationary solutions for magnetohydrodynamic equations. Discrete & Continuous Dynamical Systems - A, 2017, 37 (6) : 3435-3465. doi: 10.3934/dcds.2017146

[14]

Qing Chen, Zhong Tan. Global existence in critical spaces for the compressible magnetohydrodynamic equations. Kinetic & Related Models, 2012, 5 (4) : 743-767. doi: 10.3934/krm.2012.5.743

[15]

Chuangye Liu, Zhi-Qiang Wang. Synchronization of positive solutions for coupled Schrödinger equations. Discrete & Continuous Dynamical Systems - A, 2018, 38 (6) : 2795-2808. doi: 10.3934/dcds.2018118

[16]

Sho Nanao, Hiroyuki Masuyama, Shoji Kasahara, Yutaka Takahashi. Queueing analysis of data block synchronization mechanism in peer-to-peer based video streaming system. Journal of Industrial & Management Optimization, 2011, 7 (3) : 699-716. doi: 10.3934/jimo.2011.7.699

[17]

Dongfen Bian, Boling Guo. Global existence and large time behavior of solutions to the electric-magnetohydrodynamic equations. Kinetic & Related Models, 2013, 6 (3) : 481-503. doi: 10.3934/krm.2013.6.481

[18]

Youcef Amirat, Laurent Chupin, Rachid Touzani. Weak solutions to the equations of stationary magnetohydrodynamic flows in porous media. Communications on Pure & Applied Analysis, 2014, 13 (6) : 2445-2464. doi: 10.3934/cpaa.2014.13.2445

[19]

Zhong Tan, Qiuju Xu, Huaqiao Wang. Global existence and convergence rates for the compressible magnetohydrodynamic equations without heat conductivity. Discrete & Continuous Dynamical Systems - A, 2015, 35 (10) : 5083-5105. doi: 10.3934/dcds.2015.35.5083

[20]

Yong Zhou, Jishan Fan. Local well-posedness for the ideal incompressible density dependent magnetohydrodynamic equations. Communications on Pure & Applied Analysis, 2010, 9 (3) : 813-818. doi: 10.3934/cpaa.2010.9.813

 Impact Factor: 

Article outline

Figures and Tables

[Back to Top]