A three-dimensional dynamical model of Amundsenisen icefield, Svalbard.

Amundsenisen is an ice field, 80 km2 in area, located in Southern Spitsbergen, Svalbard. Radio-echo sounding measurements at 20 MHz show high intensity returns from a nearly flat basal reflector at four zones, all of them with ice thickness larger than 500m. These reflections suggest possible subglacial lakes. To determine whether basal liquid water is compatible with current pressure and temperature conditions, we aim at applying a thermo mechanical model with a free boundary at the bed defined as solution of a Stefan problem for the interface ice-subglaciallake. The complexity of the problem suggests the use of a bi-dimensional model, but this requires that well-defined flowlines across the zones with suspected subglacial lakes are available. We define these flow lines from the solution of a three-dimensional dynamical model, and this is the main goal of the present contribution. We apply a three-dimensional full-Stokes model of glacier dynamics to Amundsenisen icefield. We are mostly interested in the plateau zone of the icefield, so we introduce artificial vertical boundaries at the heads of the main outlet glaciers draining Amundsenisen. At these boundaries we set velocity boundary conditions. Velocities near the centres of the heads of the outlets are known from experimental measurements. The velocities at depth are calculated according to a SIA velocity-depth profile, and those at the rest of the transverse section are computed following Nye’s (1952) model. We select as southeastern boundary of the model domain an ice divide, where we set boundary conditions of zero horizontal velocities and zero vertical shear stresses. The upper boundary is a traction-free boundary. For the basal boundary conditions, on the zones of suspected subglacial lakes we set free-slip boundary conditions, while for the rest of the basal boundary we use a friction law linking the sliding velocity to the basal shear stress,in such a way that, contrary to the shallow ice approximation, the basal shear stress is not equal to the basal driving stress but rather part of the solution.

Analysis of Lead-Acid batteries with a third-order dynamical model

Electrical energy storage is a really important issue nowadays. As electricity is not easy to be directly stored, it can be stored in other forms and converted back to electricity when needed. As a consequence, storage technologies for electricity can be classified by the form of storage, and in particular we focus on electrochemical energy storage systems, better known as electrochemical batteries. Largely the more widespread batteries are the Lead-Acid ones, in the two main types known as flooded and valve-regulated. Batteries need to be present in many important applications such as in renewable energy systems and in motor vehicles. Consequently, in order to simulate these complex electrical systems, reliable battery models are needed. Although there exist some models developed by experts of chemistry, they are too complex and not expressed in terms of electrical networks. Thus, they are not convenient for a practical use by electrical engineers, who need to interface these models with other electrical systems models, usually described by means of electrical circuits. There are many techniques available in literature by which a battery can be modeled. Starting from the Thevenin based electrical model, it can be adapted to be more reliable for Lead-Acid battery type, with the addition of a parasitic reaction branch and a parallel network. The third-order formulation of this model can be chosen, being a trustworthy general-purpose model, characterized by a good ratio between accuracy and complexity. Considering the equivalent circuit network, all the useful equations describing the battery model are discussed, and then implemented one by one in Matlab/Simulink. The model has been finally validated, and then used to simulate the battery behaviour in different typical conditions.

A population dynamical model of Operophtera brumata, L. extended by climatic factors

Setting out from the database of Operophtera brumata, L. in between 1973 and 2000 due to the Light Trap Network in Hungary, we introduce a simple theta-logistic population dynamical model based on endogenous and exogenous factors, only. We create an indicator set from which we can choose some elements with which we can improve the fitting results the most effectively. Than we extend the basic simple model with additive climatic factors. The parameter optimization is based on the minimized root mean square error. The best model is chosen according to the Akaike Information Criterion. Finally we run the calibrated extended model with daily outputs of the regional climate model RegCM3.1, regarding 1961-1990 as reference period and 2021-2050 with 2071-2100 as future predictions. The results of the three time intervals are fitted with Beta distributions and compared statistically. The expected changes are discussed.

Avoiding healthy cells extinction in a cancer model

We consider a dynamical model of cancer growth including three interacting cell populations of tumor cells, healthy host cells and immune effector cells. For certain parameter choice, the dynamical system displays chaotic motion and by decreasing the response of the immune system to the tumor cells, a boundary crisis leading to transient chaotic dynamics is observed. This means that the system behaves chaotically for a finite amount of time until the unavoidable extinction of the healthy and immune cell populations occurs. Our main goal here is to apply a control method to avoid extinction. For that purpose, we apply the partial control method, which aims to control transient chaotic dynamics in the presence of external disturbances. As a result, we have succeeded to avoid the uncontrolled growth of tumor cells and the extinction of healthy tissue. The possibility of using this method compared to the frequently used therapies is discussed. (C) 2014 Elsevier Ltd. All rights reserved.

Cournot model with investments to change the market size

We present a new deterministic dynamical model on the market size of Cournot competitions, based on Nash equilibria of R&D investment strategies to increase the size of the market of the firms at every period of the game. We compute the unique Nash equilibrium for the second subgame and the profit functions for both firms. Adding uncertainty to the R&D investment strategies, we get a new stochastic dynamical model and we analyse the importance of the uncertainty to reverse the initial advantage of one firm with respect to the other.

Hammerstein model identification algorithm using Bezier-Bernstein approximation

A new identification algorithm is introduced for the Hammerstein model consisting of a nonlinear static function followed by a linear dynamical model. The nonlinear static function is characterised by using the Bezier-Bernstein approximation. The identification method is based on a hybrid scheme including the applications of the inverse of de Casteljau's algorithm, the least squares algorithm and the Gauss-Newton algorithm subject to constraints. The related work and the extension of the proposed algorithm to multi-input multi-output systems are discussed. Numerical examples including systems with some hard nonlinearities are used to illustrate the efficacy of the proposed approach through comparisons with other approaches.

Relaxing the Tropics to an ‘observed’ state: analysis using a simple baroclinic model

The technique of relaxation of the tropical atmosphere towards an analysis in a month-season forecast model has previously been successfully exploited in a number of contexts. Here it is shown that when tropical relaxation is used to investigate the possible origin of the observed anomalies in June–July 2007, a simple dynamical model is able to reproduce the observed component of the pattern of anomalies given by an ensemble of ECMWF forecast runs. Following this result, the simple model is used for a range of experiments on time-scales of relaxation, variables and regions relaxed based on a control model run with equatorial heating in a zonal flow. A theory based on scale analysis for the large-scale tropics is used to interpret the results. Typical relationships between scales are determined from the basic equations, and for a specified diabatic heating a chain of deductions for determining the dependent variables is derived. Different critical time-scales are found for tropical relaxation of different dependent variables to be effective. Vorticity has the longest critical time-scale, typically 1.2 days. For temperature and divergence, the time-scales are 10 hours and 3 hours, respectively. However not all the tropical fields, in particular the vertical motion, are reproduced correctly by the model unless divergence is heavily damped. To obtain the correct extra-tropical fields, it is crucial to have the correct rotational flow in the subtropics to initiate the Rossby wave propagation from there. It is sufficient to relax vorticity or temperature on a time-scale comparable or less than their critical time-scales to obtain this. However if the divergent advection of vorticity is important in the Rossby Wave Source then strong relaxation of divergence is required to accurately represent the tropical forcing of Rossby waves.

A Hamiltonian weak-wave model for shallow-water flow

A reduced dynamical model is derived which describes the interaction of weak inertia–gravity waves with nonlinear vortical motion in the context of rotating shallow–water flow. The formal scaling assumptions are (i) that there is a separation in timescales between the vortical motion and the inertia–gravity waves, and (ii) that the divergence is weak compared to the vorticity. The model is Hamiltonian, and possesses conservation laws analogous to those in the shallow–water equations. Unlike the shallow–water equations, the energy invariant is quadratic. Nonlinear stability theorems are derived for this system, and its linear eigenvalue properties are investigated in the context of some simple basic flows.

Electromagnetic structure and weak decay of meson K in a light-front QCD-inspired model

The kaon electromagnetic (e.m.) form factor is reviewed considering a light-front constituent quark model. In this approach, it is discussed the relevance of the quark-antiquark pair terms for the full covariance of the e.m. current. It is also verified, by considering a QCD dynamical model, that a good agreement with experimental data can be obtained for the kaon weak decay constant once a probability of about 80% of the valence component is taken into account.

Electromagnetic structure and weak decay of pseudoscalar mesons in a light-front QCD-inspired model

We study the scaling of the S-3(1)-S-1(0) meson mass splitting and the pseudoscalar weak-decay constants with the mass of the meson, as seen in the available experimental data. We use an effective light-front QCD-inspired dynamical model regulated at short distances to describe the valence component of the pseudoscalar mesons. The experimentally known values of the mass splitting, decay constants (from global lattice-QCD averages) and the pion charge form factor up to 4 [GeV/c](2) are reasonably described by the model.

Positronium-hydrogen-atom scattering in a five-state model

The scattering of orthopositronium (Ps) by hydrogen atoms has been investigated in a five-state coupled-channel model allowing for Ps(1s)H(2s,2p) and Ps(2s,2p)H(1s) excitations using a recently proposed electron-exchange model potential. The higher (n greater than or equal to 3) excitations and ionization of the Ps atom are calculated using the first Born approximation. Calculations are reported of scattering lengths, phase shifts. elastic, Ps and H excitation, and total cross sections. Remarkable correlations are observed between the S-wave Ps-H binding energy and the singlet scattering length, effective range, and resonance energy obtained in various model calculations. These correlations suggest that if a Ps-H dynamical model yields the correct result for one of these four observables, it is expected to lead to the correct result for the other three. The present model, which is constructed so as to reproduce the Ps-H resonance at 4.01 eV, automatically yields a Ps-H bound state at - 1.05 eV that compares well with the accurate value of - 1.067 eV. The model leads to a singlet scattering length of 3.72a(0) and effective range of 1.67a(0), whereas the correlations suggest the precise values of 3.50a(0) and 1.65a(0) for these observables, respectively. [S1050-2947(99)07703-3].

Gap solitons in a model of a superfluid fermion gas in optical lattices

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

NUCLEAR ALIGNMENT - CLASSICAL DYNAMIC-MODEL FOR THE U-238 U-238 SYSTEM

Dynamical properties of the U-238-U-238 system at the classical turning point, specifically the distance of closest approach, the relative orientations of the nuclei, and deformations have been studied at the sub-Coulomb energy of E(lab) = 6.07 MeV/nucleon using a classical dynamical model with a variable moment of inertia. Probability of favorable alignment for anomalous positron-electron pair emission through vacuum decay is calculated. The calculated small favorable alignment probability value of 0.116 is found to be enhanced by about 16% in comparison with the results of a similar study using a fixed moment of inertia as well as the results from a semiquantal calculation reported earlier.

Decomposition driven interface evolution for layers of binary mixtures: I. Model dirivation and base states

A dynamical model is proposed to describe the coupled decomposition and profile evolution of a free surfacefilm of a binary mixture. An example is a thin film of a polymer blend on a solid substrate undergoing simultaneous phase separation and dewetting. The model is based on model-H describing the coupled transport of the mass of one component (convective Cahn-Hilliard equation) and momentum (Navier-Stokes-Korteweg equations) supplemented by appropriate boundary conditions at the solid substrate and the free surface. General transport equations are derived using phenomenological nonequilibrium thermodynamics for a general nonisothermal setting taking into account Soret and Dufour effects and interfacial viscosity for the internal diffuse interface between the two components. Focusing on an isothermal setting the resulting model is compared to literature results and its base states corresponding to homogeneous or vertically stratified flat layers are analyzed.