Computing streamfunction and velocity potential in a limited domain of arbitrary shape. Part I: Theory and integral formulae
Tóm tắt
The non-uniqueness of solution and compatibility between the coupled boundary conditions in computing velocity potential and streamfunction from horizontal velocity in a limited domain of arbitrary shape are revisited theoretically with rigorous mathematic treatments. Classic integral formulas and their variants are used to formulate solutions for the coupled problems. In the absence of data holes, the total solution is the sum of two integral solutions. One is the internally induced solution produced purely and uniquely by the domain internal divergence and vorticity, and its two components (velocity potential and streamfunction) can be constructed by applying Green’s function for Poisson equation in unbounded domain to the divergence and vorticity inside the domain. The other is the externally induced solution produced purely but non-uniquely by the domain external divergence and vorticity, and the non-uniqueness is caused by the harmonic nature of the solution and the unknown divergence and vorticity distributions outside the domain. By setting either the velocity potential (or streamfunction) component to zero, the other component of the externally induced solution can be expressed by the imaginary (or real) part of the Cauchy integral constructed using the coupled boundary conditions and solvability conditions that exclude the internally induced solution. The streamfunction (or velocity potential) for the externally induced solution can also be expressed by the boundary integral of a double-layer (or single-layer) density function. In the presence of data holes, the total solution includes a data-hole-induced solution in addition to the above internally and externally induced solutions.
Tài liệu tham khảo
Batchelor, G. K., 1967: An Introduction to Fluid Dynamics. Cambridge University Press, 615pp.
Bijlsma, S. J., L. M. Hafkensheid, and P. Lynch, 1986: Computation of the streamfunction and velocity potential and reconstruction of the wind field. Mon. Wea. Rev., 114, 1547–1551.
Cao, J., and Q. Xu, 2011: Computing streamfunction and velocity potential in a limited domain of arbitrary shape. Part II: Numerical methods and test experiments. Adv. Atmos. Sci., doi: 10.1007/s00376-011-0186-5. (in press)
Chen, Q. S., and Y. H Kuo, 1992a: A harmonic-sine series expansion and its application to the partitioning and reconstruction problem in a limited area. Mon. Wea. Rev., 120, 91–112.
Chen, Q. S., and Y. H. Kuo, 1992b: A consistency condition for the wind field reconstruction in a limited area and a harmonic-cosine series expansion. Mon. Wea. Rev., 120, 2653–2670.
Daley, R., 1991: Atmospheric Data Analysis. Cambridge University Press, 457pp.
Daley, R. and E. Barker, 2001: NAVDAS: Formulation and Diagnostics. Mon. Wea. Rev., 129, 869–883.
Davies-Jones, R., 1988: On the formulation of surface geostrophic streamfunction. Mon. Wea. Rev., 116, 1824–1826.
Haltiner, G. J., and R. T. Williams, 1980: Numerical Prediction and Dynamic Meteorology. 2nd ed., John Wiley & Sons, 477pp.
Hawkins, H. F., and S. l. Rosenthal, 1965: On the computation of stream functions from the wind field. Mon. Wea. Rev., 93, 245–252.
Hollingsworth, A., and P. Lönnberg, 1986: The statistical structure of short-range forecast errors as determined from radiosonde data. Part I: The wind field. Tellus, 38A, 111–136.
Li, Z., Y. Chao, and J. C. McWilliams, 2006: Computation of the streamfunction and velocity potential for limited and irregular domains. Mon. Wea. Rev., 134, 3384–3394.
Lynch, P., 1988: Deducing the wind from vorticity and divergence. Mon. Wea. Rev., 116, 86–93.
Lynch, P., 1989: Partitioning the wind in a limited domain. Mon. Wea. Rev., 117, 1492–1500.
Parrish, D. F., and J. C. Derber, 1992: The National Meteorological Center’s spectral statistical-interpolation analysis system. Mon. Wea. Rev., 20, 1747–1763.
Sangster, W. E., 1960: A method of representing the horizontal pressure force without reduction of station pressure to sea level. J. Meteor., 17, 166–176.
Shukla, J., and K. R. Saha, 1974: Computation of nondivergent streamfunction and irrotational velocity potential velocity from the observed winds. Mon. Wea. Rev., 102, 419–425.
Smirnov, V. I., 1964a: A Course of Higher Mathematics. Vol. 3, D. E. Brown and J. N. Sneddon, Trans., Pergamon Press, Oxford, 702pp.
Smirnov, V. I., 1964b: A Course of Higher Mathematics. Vol. 5, D. E. Brown and J. N. Sneddon, Trans., Pergamon Press, Oxford, 814pp.
Wrobel, L. C., and M. H. Aliabadi, 2002: The Boundary Element Method. Wiley, New Jersey, 468pp.
Xu, Q., 2005: Representations of inverse covariances by differential operators. Adv. Atmos. Sci., 22, 181–198.
Xu, Q., and L. Wei, 2001: Estimation of three-dimensional error covariances. Part II: Analysis of wind innovation vectors. Mon. Wea. Rev., 129, 2939–2954.
Xu, Q., and J. Gong, 2003: Background error covariance functions for Doppler radial-wind analysis. Quart. J. Roy. Meteor. Soc., 129, 1703–1720.
Xu, Q., S. Liu, and M. Xue, 2006: Background error covariance functions for vector wind analyses using Doppler radar radial-velocity observations. Quart. J. Roy. Meteor. Soc., 132, 2887–2904.
Xu, Q., K. Nai, and L. Wei, 2007: An innovation method for estimating radar radial-velocity observation error and background wind error covariances. Quart. J. Roy. Meteor. Soc., 133, 407–415.
Zhao, Q., J. Cook, Q. Xu, and P. Harasti, 2006: Using radar wind observations to improve mesoscale numerical weather prediction. Wea. Forecasting, 21, 502–522.