Numerical solution of some nonlocal, nonlinear dispersive wave equations

We use a spectral method to solve numerically two nonlocal, nonlinear, dispersive, integrable wave equations, the Benjamin-Ono and the Intermediate Long Wave equations. The proposed numerical method is able to capture well the dynamics of the solutions; we use it to investigate the behaviour of solitary wave solutions of the equations with special attention to those, among the properties usually connected with integrability, for which there is at present no analytic proof. Thus we study in particular the resolution property of arbitrary initial profiles into sequences of solitary waves for both equations and clean interaction of Benjamin-Ono solitary waves. We also verify numerically that the behaviour of the solution of the Intermediate Long Wave equation as the model parameter tends to the infinite depth limit is the one predicted by the theory.

Error estimates for a fully discrete spectral scheme for a class of nonlinear, nonlocal dispersive wave equations

We analyze a fully discrete spectral method for the numerical solution of the initial- and periodic boundary-value problem for two nonlinear, nonlocal, dispersive wave equations, the Benjamin–Ono and the Intermediate Long Wave equations. The equations are discretized in space by the standard Fourier–Galerkin spectral method and in time by the explicit leap-frog scheme. For the resulting fully discrete, conditionally stable scheme we prove an L2-error bound of spectral accuracy in space and of second-order accuracy in time.

Asymptotic approach for the nonlinear equatorial long wave interactions

In the present work we use an asymptotic approach to obtain the long wave equations. The shallow water equation is put as a function of an external parameter that is a measure of both the spatial scales anisotropy and the fast to slow time ratio. The values given to the external parameters are consistent with those computed using typical values of the perturbations in tropical dynamics. Asymptotically, the model converge toward the long wave model. Thus, it is possible to go toward the long wave approximation through intermediate realizable states. With this approach, the resonant nonlinear wave interactions are studied. To simplify, the reduced dynamics of a single resonant triad is used for some selected equatorial trios. It was verified by both theoretical and numerical results that the nonlinear energy exchange period increases smoothly as we move toward the long wave approach. The magnitude of the energy exchanges is also modified, but in this case depends on the particular triad used and also on the initial energy partition among the triad components. Some implications of the results for the tropical dynamics are disccussed. In particular, we discuss the implications of the results for El Nĩo and the Madden-Julian in connection with other scales of time and spatial variability. © Published under licence by IOP Publishing Ltd.

Runge-Kutta IMEX schemes for the Horizontally Explicit/Vertically Implicit (HEVI) solution of wave equations

Many operational weather forecasting centres use semi-implicit time-stepping schemes because of their good efficiency. However, as computers become ever more parallel, horizontally explicit solutions of the equations of atmospheric motion might become an attractive alternative due to the additional inter-processor communication of implicit methods. Implicit and explicit (IMEX) time-stepping schemes have long been combined in models of the atmosphere using semi-implicit, split-explicit or HEVI splitting. However, most studies of the accuracy and stability of IMEX schemes have been limited to the parabolic case of advection–diffusion equations. We demonstrate how a number of Runge–Kutta IMEX schemes can be used to solve hyperbolic wave equations either semi-implicitly or HEVI. A new form of HEVI splitting is proposed, UfPreb, which dramatically improves accuracy and stability of simulations of gravity waves in stratified flow. As a consequence it is found that there are HEVI schemes that do not lose accuracy in comparison to semi-implicit ones. The stability limits of a number of variations of trapezoidal implicit and some Runge–Kutta IMEX schemes are found and the schemes are tested on two vertical slice cases using the compressible Boussinesq equations split into various combinations of implicit and explicit terms. Some of the Runge–Kutta schemes are found to be beneficial over trapezoidal, especially since they damp high frequencies without dropping to first-order accuracy. We test schemes that are not formally accurate for stiff systems but in stiff limits (nearly incompressible) and find that they can perform well. The scheme ARK2(2,3,2) performs the best in the tests.

Long-wave and short-wave asymptotics in nonlinear dispersive systems

In this paper we study the interplay between short- and long-space scales in the context of conservative dispersive systems. We consider model systems in (1 + 1) dimensions that admit both long- and short-wavelength solutions in the linear regime. A nonlinear analysis of these systems is constructed, making use of multiscale expansions. We show that the equations governing the lowest order involve only short-wave properties and that the long-wave effects to leading order are determined by a secularity elimination procedure. © 1999 The American Physical Society.

Analytical solutions of the multi-term modified power law wave equations in a finite domain

Fractional partial differential equations with more than one fractional derivative term in time, such as the Szabo wave equation, or the power law wave equation, describe important physical phenomena. However, studies of these multi-term time-space or time fractional wave equations are still under development. In this paper, multi-term modified power law wave equations in a finite domain are considered. The multi-term time fractional derivatives are defined in the Caputo sense, whose orders belong to the intervals (1, 2], [2, 3), [2, 4) or (0, n) (n > 2), respectively. Analytical solutions of the multi-term modified power law wave equations are derived. These new techniques are based on Luchko’s Theorem, a spectral representation of the Laplacian operator, a method of separating variables and fractional derivative techniques. Then these general methods are applied to the special cases of the Szabo wave equation and the power law wave equation. These methods and techniques can also be extended to other kinds of the multi term time-space fractional models including fractional Laplacian.

Effect of long-wave UV radiation on mouse melanoma: an in vitro and in vivo study

The skin cancer incidence has increased substantially over the past decades and the role of ultraviolet (UV) radiation in the etiology of skin cancer is well established. Ultraviolet B radiation (280-320 nm) is commonly considered as the more harmful part of the UV-spectrum due to its DNA-damaging potential and well-known carcinogenic effects. Ultraviolet A radiation (320-400 nm) is still regarded as a relatively low health hazard. However, UVA radiation is the predominant component in sunlight, constituting more than 90% of the environmentally relevant solar ultraviolet radiation. In the light of the recent scientific evidence, UVA has been shown to have genotoxic and immunologic effects, and it has been proposed that UVA plays a significant role in the development of skin cancer. Due to the popularity of skin tanning lamps, which emit high intensity UVA radiation and because of the prolonged sun tanning periods with the help of effective UVB blockers, the potential deleterious effects of UVA has emerged as a source of concern for public health. The possibility that UV radiation may affect melanoma metastasis has not been addressed before. UVA radiation can modulate various cellular processes, some of which might affect the metastatic potential of melanoma cells. The aim of the present study was to investigate the possible role of UVA irradiation on the metastatic capacity of mouse melanoma both in vitro and in vivo. The in vitro part of the study dealt with the enhancement of the intercellular interactions occurring either between tumor cells or between tumor cells and endothelial cells after UVA irradiation. The use of the mouse melanoma/endothelium in vitro model showed that a single-dose of UVA to melanoma cells causes an increase in melanoma cell adhesiveness to non-irradiated endothelium after 24-h irradiation. Multiple-dose irradiation of melanoma cells already increased adhesion at a 1-h time-point, which suggests the possible cumulative effect of multiple doses of UVA irradiation. This enhancement of adhesiveness might lead to an increase in binding tumor cells to the endothelial lining of vasculature in various internal organs if occurring also in vivo. A further novel observation is that UVA induced both decline in the expression of E-cadherin adhesion molecule and increase in the expression of the N-cadherin adhesion molecule. In addition, a significant decline in homotypic melanoma-melanoma adhesion (clustering) was observed, which might result in the reduction of E-cadherin expression. The aim of the in vivo animal study was to confirm the physiological significance of previously obtained in vitro results and to determine whether UVA radiation might increase melanoma metastasis in vivo. The use of C57BL/6 mice and syngeneic melanoma cell lines B16-F1 and B16-F10 showed that mice, which were i.v. injected with B16-F1 melanoma cells and thereafter exposed to UVA developed significantly more lung metastases when compared with the non-UVA-exposed group. To study the mechanism behind this phenomenon, the direct effect of UVA-induced lung colonization capacity was examined by the in vitro exposure of B16-F1 cells. Alternatively, the UVA-induced immunosuppression, which might be involved in increased melanoma metastasis, was measured by standard contact hypersensitivity assay (CHS). It appears that the UVA-induced increase of metastasis in vivo might be caused by a combination of UVA-induced systemic immunosuppression, and to the lesser extent, it might be caused by the increased adhesiveness of UVA irradiated melanoma cells. Finally, the UVA effect on gene expression in mouse melanoma was determined by a cDNA array, which revealed UVA-induced changes in the 9 differentially expressed genes that are involved in angiogenesis, cell cycle, stress-response, and cell motility. These results suggest that observed genes might be involved in cellular response to UVA and a physiologically relevant UVA dose have previously unknown cellular implications. The novel results presented in this thesis offer evidence that UVA exposure might increase the metastatic potential of the melanoma cells present in blood circulation. Considering the wellknown UVA-induced deleterious effects on cellular level, this study further supports the notion that UVA radiation might have more potential impact on health than previously suggested. The possibility of the pro-metastatic effects of UVA exposure might not be of very high significance for daily exposures. However, UVA effects might gain physiological significance following extensive sunbathing or solaria tanning periods. Whether similar UVA-induced pro-metastatic effects occur in people sunbathing or using solaria remains to be determined. In the light of the results presented in this thesis, the avoidance of solaria use could be well justified.

Retrieval of dust aerosols during night: improved assessment of long wave dust radiative forcing over Afro-Asian regions

Several investigators in the past have used the radiance depression (with respect to clear-sky infrared radiance), resulting from the presence of mineral dust aerosols in the atmosphere, as an index of dust aerosol load in the atmosphere during local noon. Here, we have used a modified approach to retrieve dust index during night since assessment of diurnal average infrared dust forcing essentially requires information on dust aerosols during night. For this purpose, we used infrared radiance (10.5-12.5 mu m), acquired from the METEOSAT-5 satellite (similar to 5 km resolution). We found that the `dust index' algorithm, valid for daytime, will no longer hold during the night because dust is then hotter than the theoretical dust-free reference. Hence we followed a `minimum reference' approach instead of a conventional `maximum reference' approach. A detailed analysis suggests that the maximum dust load occurs during the daytime. Over the desert regions of India and Africa, maximum change in dust load is as much as a factor of four between day and night and factor of two variations are commonly observed. By realizing the consequent impact on long wave dust forcing, sensitivity studies were carried out, which indicate that utilizing day time data for estimating the diurnally averaged long-wave dust radiative forcing results in significant errors (as much as 50 to 70%). Annually and regionally averaged long wave dust radiative forcing (which account for the diurnal variation of dust) at the top of the atmosphere over Afro-Asian region is 2.6 +/- 1.8 W m(-2), which is 30 to 50% lower than those reported earlier. Our studies indicate that neglecting diurnal variation of dust while assessing its radiative impact leads to an overestimation of dust radiative forcing, which in turn result in underestimation of the radiative impact of anthropogenic aerosols.

A characteristic based numerical method with tracking for nonlinear wave equations

Many physical problems can be modeled by scalar, first-order, nonlinear, hyperbolic, partial differential equations (PDEs). The solutions to these PDEs often contain shock and rarefaction waves, where the solution becomes discontinuous or has a discontinuous derivative. One can encounter difficulties using traditional finite difference methods to solve these equations. In this paper, we introduce a numerical method for solving first-order scalar wave equations. The method involves solving ordinary differential equations (ODEs) to advance the solution along the characteristics and to propagate the characteristics in time. Shocks are created when characteristics cross, and the shocks are then propagated by applying analytical jump conditions. New characteristics are inserted in spreading rarefaction fans. New characteristics are also inserted when values on adjacent characteristics lie on opposite sides of an inflection point of a nonconvex flux function, Solutions along characteristics are propagated using a standard fourth-order Runge-Kutta ODE solver. Shocks waves are kept perfectly sharp. In addition, shock locations and velocities are determined without analyzing smeared profiles or taking numerical derivatives. In order to test the numerical method, we study analytically a particular class of nonlinear hyperbolic PDEs, deriving closed form solutions for certain special initial data. We also find bounded, smooth, self-similar solutions using group theoretic methods. The numerical method is validated against these analytical results. In addition, we compare the errors in our method with those using the Lax-Wendroff method for both convex and nonconvex flux functions. Finally, we apply the method to solve a PDE with a convex flux function describing the development of a thin liquid film on a horizontally rotating disk and a PDE with a nonconvex flux function, arising in a problem concerning flow in an underground reservoir.

Frequency Domain Based Solution for Certain Class of Wave Equations: An exhaustive study of Numerical Solutions

The paper discusses the frequency domain based solution for a certain class of wave equations such as: a second order partial differential equation in one variable with constant and varying coefficients (Cantilever beam) and a coupled second order partial differential equation in two variables with constant and varying coefficients (Timoshenko beam). The exact solution of the Cantilever beam with uniform and varying cross-section and the Timoshenko beam with uniform cross-section is available. However, the exact solution for Timoshenko beam with varying cross-section is not available. Laplace spectral methods are used to solve these problems exactly in frequency domain. The numerical solution in frequency domain is done by discretisation in space by approximating the unknown function using spectral functions like Chebyshev polynomials, Legendre polynomials and also Normal polynomials. Different numerical methods such as Galerkin Method, Petrov- Galerkin method, Method of moments and Collocation method or the Pseudo-spectral method in frequency domain are studied and compared with the available exact solution. An approximate solution is also obtained for the Timoshenko beam with varying cross-section using Laplace Spectral Element Method (LSEM). The group speeds are computed exactly for the Cantilever beam and Timoshenko beam with uniform cross-section and is compared with the group speeds obtained numerically. The shear mode and the bending modes of the Timoshenko beam with uniform cross-section are separated numerically by applying a modulated pulse as the shear force and the corresponding group speeds for varying taper parameter in are obtained numerically by varying the frequency of the input pulse. An approximate expression for calculating group speeds corresponding to the shear mode and the bending mode, and also the cut-off frequency is obtained. Finally, we show that the cut-off frequency disappears for large in, for epsilon > 0 and increases for large in, for epsilon < 0.

Long-Wave Dynamics of Sea Level Variations during Indian Ocean Dipole Events

Long-wave dynamics of the interannual variations of the equatorial Indian Ocean circulation are studied using an ocean general circulation model forced by the assimilated surface winds and heat flux of the European Centre for Medium-Range Weather Forecasts. The simulation has reproduced the sea level anomalies of the Ocean Topography Experiment (TOPEX)/Poseidon altimeter observations well. The equatorial Kelvin and Rossby waves decomposed from the model simulation show that western boundary reflections provide important negative feedbacks to the evolution of the upwelling currents off the Java coast during Indian Ocean dipole (IOD) events. Two downwelling Kelvin wave pulses are generated at the western boundary during IOD events: the first is reflected from the equatorial Rossby waves and the second from the off-equatorial Rossby waves in the southern Indian Ocean. The upwelling in the eastern basin during the 1997-98 IOD event is weakened by the first Kelvin wave pulse and terminated by the second. In comparison, the upwelling during the 1994 IOD event is terminated by the first Kelvin wave pulse because the southeasterly winds off the Java coast are weak at the end of 1994. The atmospheric intraseasonal forcing, which plays an important role in inducing Java upwelling during the early stage of an IOD event, is found to play a minor role in terminating the upwelling off the Java coast because the intraseasonal winds are either weak or absent during the IOD mature phase. The equatorial wave analyses suggest that the upwelling off the Java coast during IOD events is terminated primarily by western boundary reflections.

An estimate of the global impact of multiple scattering by clouds on outgoing long-wave radiation

Although the potential importance of scattering of long-wave radiation by clouds has been recognised, most studies have concentrated on the impact of high clouds and few estimates of the global impact of scattering have been presented. This study shows that scattering in low clouds has a significant impact on outgoing long-wave radiation (OLR) in regions of marine stratocumulus (-3.5 W m(-2) for overcast conditions) where the column water vapour is relatively low. This corresponds to an enhancement of the greenhouse effect of such clouds by 10%. The near-global impact of scattering on OLR is estimated to be -3.0 W m(-2), with low clouds contributing -0.9 W m(-2), mid-level cloud -0.7 W m(-2) and high clouds -1.4 W m(-2). Although this effect appears small compared to the global mean OLR of 240 W m(-2), it indicates that neglect of scattering will lead to an error in cloud long-wave forcing of about 10% and an error in net cloud forcing of about 20%.

Examination of long-wave radiative bias in general circulation models over North Africa during May–July

Satellite data are used to quantify and examine the bias in the outgoing long-wave (LW) radiation over North Africa during May–July simulated by a range of climate models and the Met Office global numerical weather prediction (NWP) model. Simulations from an ensemble-mean of multiple climate models overestimate outgoing clear-sky long-wave radiation (LWc) by more than 20 W m−2 relative to observations from Clouds and the Earth's Radiant Energy System (CERES) for May–July 2000 over parts of the west Sahara, and by 9 W m−2 for the North Africa region (20°W–30°E, 10–40°N). Experiments with the atmosphere-only version of the High-resolution Hadley Centre Global Environment Model (HiGEM), suggest that including mineral dust radiative effects removes this bias. Furthermore, only by reducing surface temperature and emissivity by unrealistic amounts is it possible to explain the magnitude of the bias. Comparing simulations from the Met Office NWP model with satellite observations from Geostationary Earth Radiation Budget (GERB) instruments suggests that the model overestimates the LW by 20–40 W m−2 during North African summer. The bias declines over the period 2003–2008, although this is likely to relate to improvements in the model and inhomogeneity in the satellite time series. The bias in LWc coincides with high aerosol dust loading estimated from the Ozone Monitoring Instrument (OMI), including during the GERBILS field campaign (18–28 June 2007) where model overestimates in LWc greater than 20 W m−2 and OMI-estimated aerosol optical depth (AOD) greater than 0.8 are concurrent around 20°N, 0–20°W. A model-minus-GERB LW bias of around 30 W m−2 coincides with high AOD during the period 18–21 June 2007, although differences in cloud cover also impact the model–GERB differences. Copyright © Royal Meteorological Society and Crown Copyright, 2010