Computational modes and grid imprinting on five quasi-uniform spherical C-grids

Currently, most operational forecasting models use latitude-longitude grids, whose convergence of meridians towards the poles limits parallel scaling. Quasi-uniform grids might avoid this limitation. Thuburn et al, JCP, 2009 and Ringler et al, JCP, 2010 have developed a method for arbitrarily-structured, orthogonal C-grids (TRiSK), which has many of the desirable properties of the C-grid on latitude-longitude grids but which works on a variety of quasi-uniform grids. Here, five quasi-uniform, orthogonal grids of the sphere are investigated using TRiSK to solve the shallow-water equations. We demonstrate some of the advantages and disadvantages of the hexagonal and triangular icosahedra, a Voronoi-ised cubed sphere, a Voronoi-ised skipped latitude-longitude grid and a grid of kites in comparison to a full latitude-longitude grid. We will show that the hexagonal-icosahedron gives the most accurate results (for least computational cost). All of the grids suffer from spurious computational modes; this is especially true of the kite grid, despite it having exactly twice as many velocity degrees of freedom as height degrees of freedom. However, the computational modes are easiest to control on the hexagonal icosahedron since they consist of vorticity oscillations on the dual grid which can be controlled using a diffusive advection scheme for potential vorticity.

A variational method to retrieve the extinction profile in liquid clouds using multiple field-of-view lidar

Liquid clouds play a profound role in the global radiation budget but it is difficult to remotely retrieve their vertical profile. Ordinary narrow field-of-view (FOV) lidars receive a strong return from such clouds but the information is limited to the first few optical depths. Wideangle multiple-FOV lidars can isolate radiation scattered multiple times before returning to the instrument, often penetrating much deeper into the cloud than the singly-scattered signal. These returns potentially contain information on the vertical profile of extinction coefficient, but are challenging to interpret due to the lack of a fast radiative transfer model for simulating them. This paper describes a variational algorithm that incorporates a fast forward model based on the time-dependent two-stream approximation, and its adjoint. Application of the algorithm to simulated data from a hypothetical airborne three-FOV lidar with a maximum footprint width of 600m suggests that this approach should be able to retrieve the extinction structure down to an optical depth of around 6, and total opticaldepth up to at least 35, depending on the maximum lidar FOV. The convergence behavior of Gauss-Newton and quasi-Newton optimization schemes are compared. We then present results from an application of the algorithm to observations of stratocumulus by the 8-FOV airborne “THOR” lidar. It is demonstrated how the averaging kernel can be used to diagnose the effective vertical resolution of the retrieved profile, and therefore the depth to which information on the vertical structure can be recovered. This work enables exploitation of returns from spaceborne lidar and radar subject to multiple scattering more rigorously than previously possible.

Some numerical experiments on the effect of the variation of static stability in two-layer quasi-geostrophic models

A method to solve a quasi-geostrophic two-layer model including the variation of static stability is presented. The divergent part of the wind is incorporated by means of an iterative procedure. The procedure is rather fast and the time of computation is only 60–70% longer than for the usual two-layer model. The method of solution is justified by the conservation of the difference between the gross static stability and the kinetic energy. To eliminate the side-boundary conditions the experiments have been performed on a zonal channel model. The investigation falls mainly into three parts: The first part (section 5) contains a discussion of the significance of some physically inconsistent approximations. It is shown that physical inconsistencies are rather serious and for these inconsistent models which were studied the total kinetic energy increased faster than the gross static stability. In the next part (section 6) we are studying the effect of a Jacobian difference operator which conserves the total kinetic energy. The use of this operator in two-layer models will give a slight improvement but probably does not have any practical use in short periodic forecasts. It is also shown that the energy-conservative operator will change the wave-speed in an erroneous way if the wave-number or the grid-length is large in the meridional direction. In the final part (section 7) we investigate the behaviour of baroclinic waves for some different initial states and for two energy-consistent models, one with constant and one with variable static stability. According to the linear theory the waves adjust rather rapidly in such a way that the temperature wave will lag behind the pressure wave independent of the initial configuration. Thus, both models give rise to a baroclinic development even if the initial state is quasi-barotropic. The effect of the variation of static stability is very small, qualitative differences in the development are only observed during the first 12 hours. For an amplifying wave we will get a stabilization over the troughs and an instabilization over the ridges.

Nonlinear stability of multilayer quasi-geostrophic flow

New nonlinear stability theorems are derived for disturbances to steady basic flows in the context of the multilayer quasi-geostrophic equations. These theorems are analogues of Arnol’d's second stability theorem, the latter applying to the two-dimensional Euler equations. Explicit upper bounds are obtained on both the disturbance energy and disturbance potential enstrophy in terms of the initial disturbance fields. An important feature of the present analysis is that the disturbances are allowed to have non-zero circulation. While Arnol’d's stability method relies on the energy–Casimir invariant being sign-definite, the new criteria can be applied to cases where it is sign-indefinite because of the disturbance circulations. A version of Andrews’ theorem is established for this problem, and uniform potential vorticity flow is shown to be nonlinearly stable. The special case of two-layer flow is treated in detail, with particular attention paid to the Phillips model of baroclinic instability. It is found that the short-wave portion of the marginal stability curve found in linear theory is precisely captured by the new nonlinear stability criteria.

On Arnol'd's second nonlinear stability theorem for two-dimensional quasi-geostrophic flow

Arnol'd's second hydrodynamical stability theorem, proven originally for the two-dimensional Euler equations, can establish nonlinear stability of steady flows that are maxima of a suitably chosen energy-Casimir invariant. The usual derivations of this theorem require an assumption of zero disturbance circulation. In the present work an analogue of Arnol'd's second theorem is developed in the more general case of two-dimensional quasi-geostrophic flow, with the important feature that the disturbances are allowed to have non-zero circulation. New nonlinear stability criteria are derived, and explicit bounds are obtained on both the disturbance energy and potential enstrophy which are expressed in terms of the initial disturbance fields. While Arnol'd's stability method relies on the second variation of the energy-Casimir invariant being sign-definite, the new criteria can be applied to cases where the second variation is sign-indefinite because of the disturbance circulations. A version of Andrews' theorem is also established for this problem.

A finite element method for second order nonvariational elliptic problems

We propose a numerical method to approximate the solution of second order elliptic problems in nonvariational form. The method is of Galerkin type using conforming finite elements and applied directly to the nonvariational (nondivergence) form of a second order linear elliptic problem. The key tools are an appropriate concept of “finite element Hessian” and a Schur complement approach to solving the resulting linear algebra problem. The method is illustrated with computational experiments on three linear and one quasi-linear PDE, all in nonvariational form.

Proposal of a floating mechanism based on magnetic levitation method toward EEG-driven robotic system

We discussed a floating mechanism based on quasi-magnetic levitation method that can be attached at the endpoint of a robot arm in order to construct a novel redundant robot arm for producing compliant motions. The floating mechanism can be composed of magnets and a constraint mechanism such that the repelling force of the magnets floats the endpoint part of the mechanism stable for the guided motions. The analytical and experimental results show that the proposed floating mechanism can produce stable floating motions with small inertia and viscosity. The results also show that the proposed mechanism can detect small force applied to the endpoint part because the friction force of the mechanism is very small.

Indirect Study of the (2)H(d,p)(3)H and (2)H(d,n)(3)He Reactions at Astrophysical Energies via the Trojan Horse Method

The (2)H(d,p)(3)H and (2)H(d,n)(3)He reactions have been indirectly studied by means of the Trojan Horse Method applied to the quasi-free (2)H((3)He, p(3)H)(1)H (2)H((3)He, n(3)He)(1)H reaction at 18 MeV of beam energy. This is the first experiment where the spectator (here (1)H) has been detected in coincidence with the charged participant, avoiding the limitations of standard neutron detectors. The d - d relative energy has been measured from 1.5 MeV down to 2 keV, at center of mass angles from 40A degrees to 170A degrees. Indirect angular distributions are compared with the direct behaviour in the overlapping regions.

Low-energy d plus d fusion reactions via the Trojan Horse Method

The bare nucleus S(E) factors for the (2)H(d, p)(3)H and (2)H(d.n)(3)He reactions have been measured for the first time via the Trojan Horse Method off the proton in (3)He from 1.5 MeV down to 2 key. This range overlaps with the relevant region for Standard Big Bang Nucleosynthesis as well as with the thermal energies of future fusion reactors and deuterium burning in the Pre-Main-Sequence phase of stellar evolution. This is the first pioneering experiment in quasi free regime where the charged spectator is detected. Both the energy dependence and the absolute value of the S(E) factors deviate by more than 15% from available direct data with new S(0) values of 57.4 +/- 1.8 MeVb for (3)H + p and 60.1 +/- 1.9 MeV b for (3)He + n. None of the existing fitting curves is able to provide the correct slope of the new data in the full range, thus calling for a revision of the theoretical description. This has consequences in the calculation of the reaction rates with more than a 25% increase at the temperatures of future fusion reactors. (C) 2011 Elsevier By. All rights reserved.

From decoherence-free channels to decoherence-free and quasi-free subspaces within bosonic dissipative networks

We present a technique to build, within a dissipative bosonic network, decoherence-free channels (DFCs): a group of normal-mode oscillators with null effective damping rates. We verify that the states protected within the DFC define the well-known decoherence-free subspaces (DFSs) when mapped back into the natural network oscillators. Therefore, our technique to build protected normal-mode channels turns out to be an alternative way to build DFSs, which offers advantages over the conventional method. It enables the computation of all the network-protected states at once, as well as leading naturally to the concept of the decoherence quasi-free subspace (DQFS), inside which a superposition state is quasi-completely protected against decoherence. The concept of the DQFS, weaker than that of the DFS, may provide a more manageable mechanism to control decoherence. Finally, as an application of the DQFSs, we show how to build them for quasi-perfect state transfer in networks of coupled quantum dissipative oscillators.

Carbon supported electrocatalysts prepared by the sol-gel method and their utilization for the oxidation of methanol in acid media

One of the key objectives in fuel-cell technology is to improve the performance of the anode catalyst for the alcohol oxidation and reduce Pt loading. Here, we show the use of six different electrocatalysts synthesized by the sol -gel method on carbon powder to promote the oxidation of methanol in acid media. The catalysts Pt-PbO(x) and Pt-(RuO(2)-PbO(x)) with 10% of catalyst load exhibited significantly enhanced catalytic activity toward the methanol oxidation reaction as compared to Pt-(RuO(2))/C and Pt/C electrodes. Cyclic voltammetry studies showed that the electrocatalysts Pt-PbO(x)/C and Pt-(RuO(2)-PbO(x))/C started the oxidation process at extremely low potentials and that they represent a good novelty to oxidize methanol. Furthermore, quasi-stationary polarization experiments and cronoamperometry studies showed the good performance of the Pt-PbO(x), Pt-(RuO(2)-PbO(x))/C and Pt-(RuO(2)-IrO(2))/C catalysts during the oxidation process. Thus, the addition of metallic Pt and PbO(x) onto high-area carbon powder, by the sol -gel route, constitutes an interesting way to prepare anodes with high catalytic activity for further applications in direct methanol fuel cell systems.

Ondas de spin em quasi-cristais magnônicos

In this paper we investigate the spectra of band structures and transmittance in magnonic quasicrystals that exhibit the so-called deterministic disorders, specifically, magnetic multilayer systems, which are built obeying to the generalized Fibonacci (only golden mean (GM), silver mean (SM), bronze mean (BM), copper mean (CM) and nickel mean (NM) cases) and k-component Fibonacci substitutional sequences. The theoretical model is based on the Heisenberg Hamiltonian in the exchange regime, together with the powerful transfer matrix method, and taking into account the RPA approximation. The magnetic materials considered are simple cubic ferromagnets. Our main interest in this study is to investigate the effects of quasiperiodicity on the physical properties of the systems mentioned by analyzing the behavior of spin wave propagation through the dispersion and transmission spectra of these structures. Among of these results we detach: (i) the fragmentation of the bulk bands, which in the limit of high generations, become a Cantor set, and the presence of the mig-gap frequency in the spin waves transmission, for generalized Fibonacci sequence, and (ii) the strong dependence of the magnonic band gap with respect to the parameters k, which determines the amount of different magnetic materials are present in quasicrystal, and n, which is the generation number of the sequence k-component Fibonacci. In this last case, we have verified that the system presents a magnonic band gap, whose width and frequency region can be controlled by varying k and n. In the exchange regime, the spin waves propagate with frequency of the order of a few tens of terahertz (THz). Therefore, from a experimental and technological point of view, the magnonic quasicrystals can be used as carriers or processors of informations, and the magnon (the quantum spin wave) is responsible for this transport and processing

Sistemas Nanoestruturados: Heteroestruturas Quasi-Periódicas de Nitretos e Cálculos Ab initio em Polimorfos CaCO3

The physical properties and the excitations spectrum in oxides and semiconductors materials are presented in this work, whose the first part presents a study on the confinement of optical phonons in artificial systems based on III-V nitrides, grown in periodic and quasiperiodic forms. The second part of this work describes the Ab initio calculations which were carried out to obtain the optoeletronic properties of Calcium Oxide (CaO) and Calcium Carbonate (CaCO3) crystals. For periodic and quasi-periodic superlattices, we present some dynamical properties related to confined optical phonons (bulk and surface), obtained through simple theories, such as the dielectric continuous model, and using techniques such as the transfer-matrix method. The localization character of confined optical phonon modes, the magnitude of the bands in the spectrum and the power laws of these structures are presented as functions of the generation number of sequence. The ab initio calculations have been carried out using the CASTEP software (Cambridge Total Sequential Energy Package), and they were based on ultrasoft-like pseudopotentials and Density Functional Theory (DFT). Two di®erent geometry optimizations have been e®ectuated for CaO crystals and CaCO3 polymorphs, according to LDA (local density approximation) and GGA (generalized gradient approximation) approaches, determining several properties, e. g. lattice parameters, bond length, electrons density, energy band structures, electrons density of states, e®ective masses and optical properties, such as dielectric constant, absorption, re°ectivity, conductivity and refractive index. Those results were employed to investigate the confinement of excitons in spherical Si@CaCO3 and CaCO3@SiO2 quantum dots and in calcium carbonate nanoparticles, and were also employed in investigations of the photoluminescence spectra of CaCO3 crystal

Dissipation-managed soliton in a quasi-one-dimensional Bose-Einstein condensate

We use the time-dependent mean-field Cross-Pitaevskii equation to study the formation of a dynamically-stabilized dissipation managed bright soliton in a quasi-one dimensional Bose-Einstein condensate (BEC). Because of three-body recombination of bosonic atoms to molecules, atoms are lost (dissipated) from a BEC. Such dissipation leads to the decay of a BEC soliton. We demonstrate by a perturbation procedure that an alimentation of atoms from an external source to the BEC may compensate for the dissipation loss and lead to a dynamically-stabilized soliton. The result of the analytical perturbation method is in excellent agreement with mean-field numerics. It seems possible to obtain such a dynamically stabilized BEC soliton without dissipation in laboratory.

Photoluminescence in quasi-amorphous Pb0.8X0.2Zr0.53Ti0.47O3 (X = Ca, Sr and Ba) powders: An optical and structural study

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