The composite-tendency RAW filter in semi-implicit integrations

The time discretization in weather and climate models introduces truncation errors that limit the accuracy of the simulations. Recent work has yielded a method for reducing the amplitude errors in leapfrog integrations from first-order to fifth-order. This improvement is achieved by replacing the Robert--Asselin filter with the RAW filter and using a linear combination of the unfiltered and filtered states to compute the tendency term. The purpose of the present paper is to apply the composite-tendency RAW-filtered leapfrog scheme to semi-implicit integrations. A theoretical analysis shows that the stability and accuracy are unaffected by the introduction of the implicitly treated mode. The scheme is tested in semi-implicit numerical integrations in both a simple nonlinear stiff system and a medium-complexity atmospheric general circulation model, and yields substantial improvements in both cases. We conclude that the composite-tendency RAW-filtered leapfrog scheme is suitable for use in semi-implicit integrations.

The composite-tendency Robert-Asselin-Williams (RAW) filter in semi-implicit integrations

Timediscretization in weatherandclimate modelsintroduces truncation errors that limit the accuracy of the simulations. Recent work has yielded a method for reducing the amplitude errors in leap-frog integrations from first-order to fifth-order.This improvement is achieved by replacing the Robert–Asselin filter with the Robert–Asselin–Williams (RAW) filter and using a linear combination of unfiltered and filtered states to compute the tendency term. The purpose of the present article is to apply the composite-tendency RAW-filtered leapfrog scheme to semi-implicit integrations. A theoretical analysis shows that the stability and accuracy are unaffected by the introduction of the implicitly treated mode. The scheme is tested in semi-implicit numerical integrations in both a simple nonlinear stiff system and a medium-complexity atmospheric general circulation model and yields substantial improvements in both cases. We conclude that the composite-tendency RAW-filtered leap-frog scheme is suitable for use in semi-implicit integrations.

Energy partitioning for "fuzzy" atoms

The total energy of molecule in terms of 'fuzzy atoms' presented as sum of one- and two-atomic energy components is described. The divisions of three-dimensional physical space into atomic regions exhibit continuous transition from one to another. The energy components are on chemical energy scale according to proper definitions. The Becke's integration scheme and weight function determines realization of method which permits effective numerical integrations

Front and domain growth in the presence of gravity

Front and domain growth of a binary mixture in the presence of a gravitational field is studied. The interplay of bulk- and surface-diffusion mechanisms is analyzed. An equation for the evolution of interfaces is derived from a time-dependent Ginzburg-Landau equation with a concentration-dependent diffusion coefficient. Scaling arguments on this equation give the exponents of a power-law growth. Numerical integrations of the Ginzburg-Landau equation corroborate the theoretical analysis.

Phase-field model of Hele-Shaw flows in the high-viscosity contrast regime

A one-sided phase-field model is proposed to study the dynamics of unstable interfaces of Hele-Shaw flows in the high viscosity contrast regime. The corresponding macroscopic equations are obtained by means of an asymptotic expansion from the phase-field model. Numerical integrations of the phase-field model in a rectangular Hele-Shaw cell reproduce finger competition with the final evolution to a steady-state finger.

Dynamical scaling of imbibition in columnar geometries

compatible with the usual nonlocal model governed by surface tension that results from a macroscopic description. To explore this discrepancy, we exhaustively analyze numerical integrations of a phase-field model with dichotomic columnar disorder. We find that two distinct behaviors are possible depending on the capillary contrast between the two values of disorder. In a high-contrast case, where interface evolution is mainly dominated by the disorder, an inherent anomalous scaling is always observed. Moreover, in agreement with experimental work, the interface motion has to be described through a local model. On the other hand, in a lower-contrast case, the interface is dominated by interfacial tension and can be well modeled by a nonlocal model. We have studied both spontaneous and forced-flow imbibition situations, giving a complete set of scaling exponents in each case, as well as a comparison to the experimental results.

One- and two-center energy components in the atoms in molecules theory

The energy decomposition scheme proposed in a recent paper has been realized by performing numerical integrations. The sample calculations carried out for some simple molecules show excellent agreement with the chemical picture of molecules, indicating that such an energy decomposition analysis can be useful from the point of view of connecting quantum mechanics with the genuine chemical concepts

Energy partitioning for "fuzzy" atoms

The total energy of molecule in terms of 'fuzzy atoms' presented as sum of one- and two-atomic energy components is described. The divisions of three-dimensional physical space into atomic regions exhibit continuous transition from one to another. The energy components are on chemical energy scale according to proper definitions. The Becke's integration scheme and weight function determines realization of method which permits effective numerical integrations

Chirikov diffusion in the asteroidal three-body resonance (5,-2,-2)

The theory of diffusion in many-dimensional Hamiltonian system is applied to asteroidal dynamics. The general formulation developed by Chirikov is applied to the NesvornA1/2-Morbidelli analytic model of three-body (three-orbit) mean-motion resonances (Jupiter-Saturn-asteroid). In particular, we investigate the diffusion along and across the separatrices of the (5, -2, -2) resonance of the (490) Veritas asteroidal family and their relationship to diffusion in semi-major axis and eccentricity. The estimations of diffusion were obtained using the Melnikov integral, a Hadjidemetriou-type sympletic map and numerical integrations for times up to 10(8) years.

FROZEN ORBITS AROUND EUROPA

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Instability zones for satellites of asteroids: The example of the (87) Sylvia system

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Resonance and chaos .1. First-order interior resonances

Analytical models for studying the dynamical behaviour of objects near interior, mean motion resonances are reviewed in the context of the planar, circular, restricted three-body problem. The predicted widths of the resonances are compared with the results of numerical integrations using Poincare surfaces of section with a mass ratio of 10(-3) (similar to the Jupiter-Sun case). It is shown that for very low eccentricities the phase space between the 2:1 and 3:2 resonances is predominantly regular, contrary to simple theoretical predictions based on overlapping resonance. A numerical study of the 'evolution' of the stable equilibrium point of the 3:2 resonance as a function of the Jacobi constant shows how apocentric libration at the 2:1 resonance arises; there is evidence of a similar mechanism being responsible for the centre of the 4:3 resonance evolving towards 3:2 apocentric libration. This effect is due to perturbations from other resonances and demonstrates that resonances cannot be considered in isolation. on theoretical grounds the maximum libration width of first-order resonances should increase as the orbit of the perturbing secondary is approached. However, in reality the width decreases due to the chaotic effect of nearby resonances.

Stability regions around the components of the triple system 2001 SN263

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Resonance and chaos: I. First-order interior resonances

Analytical models for studying the dynamical behaviour of objects near interior, mean motion resonances are reviewed in the context of the planar, circular, restricted threebody problem. The predicted widths of the resonances are compared with the results of numerical integrations using Poincaré surfaces of section with a mass ratio of 10-3 (similar to the Jupiter-Sun case). It is shown that for very low eccentricities the phase space between the 2:1 and 3:2 resonances is predominantly regular, contrary to simple theoretical predictions based on overlapping resonance. A numerical study of the 'evolution' of the stable equilibrium point of the 3:2 resonance as a function of the Jacobi constant shows how apocentric libration at the 2:1 resonance arises; there is evidence of a similar mechanism being responsible for the centre of the 4:3 resonance evolving towards 3:2 apocentric libration. This effect is due to perturbations from other resonances and demonstrates that resonances cannot be considered in isolation. On theoretical grounds the maximum libration width of first-order resonances should increase as the orbit of the perturbing secondary is approached. However, in reality the width decreases due to the chaotic effect of nearby resonances.

Dynamics of some fictitious satellites of Venus and Mars

The dynamics of some fictitious satellites of Venus and Mars are studied considering only solar perturbation and the oblateness of the planet, as disturbing forces. Several numerical integrations of the averaged system, taking different values of the obliquity of ecliptic (a), show the existence of strong chaotic motion, provided that the semi major axis is near a critical value. As a consequence, large increase of eccentricities occur and the satellites may collide with the planet or cross possible internal orbits. Even starting from almost circular and equatorial orbits, most satellites can easily reach prohibitive values. The extension of the chaotic zone depends clearly on the value ε, so that, previous regular regions may become chaotic, provided ε increases sufficiently. © 1999 Elsevier Science Ltd. All rights reserved.