Multiscale Schemes for the BGK--Vlasov--Poisson System in the Quasi-Neutral and Fluid Limits. Stability Analysis and First Order Schemes

This paper deals with the development and the analysis of asymptotically stable and consistent schemes in the joint quasi-neutral and fluid limits for the collisional Vlasov-Poisson system. In these limits, the classical explicit schemes suffer from time step restrictions due to the small plasma period and Knudsen number. To solve this problem, we propose a new scheme stable for choices of time steps independent from the small scales dynamics and with comparable computational cost with respect to standard explicit schemes. In addition, this scheme reduces automatically to consistent discretizations of the underlying asymptotic systems. In this first work on this subject, we propose a first order in time scheme and we perform a relative linear stability analysis to deal with such problems. The framework we propose permits to extend this approach to high order schemes in the next future. We finally show the capability of the method in dealing with small scales through numerical experiments.

Improving a global model from the boundary layer: total turbulent energy and the neutral limit Prandtl number

Model intercomparisons have identified important deficits in the representation of the stable boundary layer by turbulence parametrizations used in current weather and climate models. However, detrimental impacts of more realistic schemes on the large-scale flow have hindered progress in this area. Here we implement a total turbulent energy scheme into the climate model ECHAM6. The total turbulent energy scheme considers the effects of Earth’s rotation and static stability on the turbulence length scale. In contrast to the previously used turbulence scheme, the TTE scheme also implicitly represents entrainment flux in a dry convective boundary layer. Reducing the previously exaggerated surface drag in stable boundary layers indeed causes an increase in southern hemispheric zonal winds and large-scale pressure gradients beyond observed values. These biases can be largely removed by increasing the parametrized orographic drag. Reducing the neutral limit turbulent Prandtl number warms and moistens low-latitude boundary layers and acts to reduce longstanding radiation biases in the stratocumulus regions, the Southern Ocean and the equatorial cold tongue that are common to many climate models.

Nonexistence of looping trajectories in hydromagnetic waves of finite amplitude. Breaking of waves in a cold collision-free plasma in a magnetic field. On stabilities of periodic waves in a cold collision-free plasma in a magnetic field.

This dissertation consists of three parts. In Part I, it is shown that looping trajectories cannot exist in finite amplitude stationary hydromagnetic waves propagating across a magnetic field in a quasi-neutral cold collision-free plasma. In Part II, time-dependent solutions in series expansion are presented for the magnetic piston problem, which describes waves propagating into a quasi-neutral cold collision-free plasma, ensuing from magnetic disturbances on the boundary of the plasma. The expansion is equivalent to Picard's successive approximations. It is then shown that orbit crossings of plasma particles occur on the boundary for strong disturbances and inside the plasma for weak disturbances. In Part III, the existence of periodic waves propagating at an arbitrary angle to the magnetic field in a plasma is demonstrated by Stokes expansions in amplitude. Then stability analysis is made for such periodic waves with respect to side-band frequency disturbances. It is shown that waves of slow mode are unstable whereas waves of fast mode are stable if the frequency is below the cutoff frequency. The cutoff frequency depends on the propagation angle. For longitudinal propagation the cutoff frequency is equal to one-fourth of the electron's gyrofrequency. For transverse propagation the cutoff frequency is so high that waves of all frequencies are stable.

Effective medium model for refractive indices of thin films with oblique columnar structure

The refractive indices of thin films, containing dielectric and voids in an oblique columnar structure, are, modelled in the quasi-static limit. The dielectric function is shown to be strongly dependent on the angle of incidence and on the columnar orientation for p-polarized light. This model is applied to model ZnS thin films with oblique columnar structures and the computed results have been given.

Extended effective medium model for refractive indices of thin films with oblique columnar structure

The refractive indices of thin films, containing dielectric and voids in an oblique columnar structure, are modeled by extended effective medium in the quasi-static limit. The dielectric function is shown to be strongly dependent on the angle of incidence and on the columnar orientation for p-polarized light. This model is applied to model ZrO2 thin films with oblique columnar structures and the computed results, with the Maxwell Garnett, the Bragg-Pippard, and the Bruggeman formalisms, have been given. (c) 2004 Elsevier B.V. All rights reserved.

Simplification of structured population dynamics in complex ecological communities

Goldstone's idea of slow dynamics resulting from spontaneously broken symmetries is applied to Hubbell's neutral hypothesis of community dynamics, to efficiently simplify stage-structured multi-species models-introducing the quasi-neutral approximation (QNA). Rather than assuming population-dynamical neutrality in the QNA, deviations from ideal neutrality, thought to be small, drive dynamics. The QNA is systematically derived to first and second order in a two-scale singular perturbation expansion. The total reproductive value of species, as computed from the effective life-history parameters resulting from the non-linear interactions with the surrounding community, emerges as the new dynamic variables in this aggregated description. Using a simple stage-structured community-assembly model, the QNA is demonstrated to accurately reproduce population dynamics in large, complex communities. Further, the utility of the QNA in building intuition for management problems is illustrated by estimating the responses of a fish stock to harvesting and variations in fecundity.

A complete analytic theory for structure and dynamics of populations and communities spanning wide ranges in body size

The prediction and management of ecosystem responses to global environmental change would profit from a clearer understanding of the mechanisms determining the structure and dynamics of ecological communities. The analytic theory presented here develops a causally closed picture for the mechanisms controlling community and population size structure, in particular community size spectra, and their dynamic responses to perturbations, with emphasis on marine ecosystems. Important implications are summarised in non-technical form. These include the identification of three different responses of community size spectra to size-specific pressures (of which one is the classical trophic cascade), an explanation for the observed slow recovery of fish communities from exploitation, and clarification of the mechanism controlling predation mortality rates. The theory builds on a community model that describes trophic interactions among size-structured populations and explicitly represents the full life cycles of species. An approximate time-dependent analytic solution of the model is obtained by coarse graining over maturation body sizes to obtain a simple description of the model steady state, linearising near the steady state, and then eliminating intraspecific size structure by means of the quasi-neutral approximation. The result is a convolution equation for trophic interactions among species of different maturation body sizes, which is solved analytically using a novel technique based on a multiscale expansion.

A table-top laser-based source of short, collimated, ultra-relativistic positron beams

An ultra-relativistic electron beam passing through a thick, high-Z solid target triggers an electromagnetic cascade, whereby a large number of high energy photons and electron-positron pairs are produced. By exploiting this physical process, we present here the first experimental evidence of the generation of ultra-short, highly collimated and ultra-relativistic positron beams following the interaction of a laser-wake field accelerated electron beam with high-Z solid targets. Clear evidence has also been obtained of the generation of GeV electron-positron jets with variable composition depending on the solid target material and thickness. The percentage of positrons in the overall leptonic beam has been observed to vary from a few per cent up to almost fifty per cent, implying a quasi-neutral electron-positron beam. We anticipate that these beams will be of direct relevance to the laboratory study of astrophysical leptonic jets and their interaction with the interstellar medium.

On vertical electric fields at lunar magnetic anomalies

We study the interaction between a magnetic dipole mimicking the Gerasimovich magnetic anomaly on the lunar surface and the solar wind in a self-consistent 3-D quasi-neutral hybrid simulation where ions are modeled as particles and electrons as a charge-neutralizing fluid. Especially, we consider the origin of the recently observed electric potentials at lunar magnetic anomalies. An antimoonward Hall electric field forms in our simulation resulting in a potential difference of <300V on the lunar surface, in which the value is similar to observations. Since the hybrid model assumes charge neutrality, our results suggest that the electric potentials at lunar magnetic anomalies can be formed by decoupling of ion and electron motion even without charge separation.

The ion beam shepherd: A new concept for asteroid deflection

A novel slow push asteroid deflection strategy has been recently proposed in which an Earth threatening asteroid can be deflected by exploiting the momentum transmitted by a collimated beam of quasi-neutral plasma impinging against the asteroid surface. The beam can be generated with state-of-the art ion engines from a hovering spacecraft with no need for physical attachment or gravitational interaction with the celestial body. The spacecraft, placed at a distance of a few asteroid diameters, would need an ion thruster pointed at the asteroid surface as well as a second propulsion system to compensate for the ion engine reaction and keep the distance between the asteroid and the shepherd satellite constant throughout the deflection phase. A comparison in terms of required spacecraft mass per total imparted deflection impulse shows that the method outperforms the gravity tractor concept by more than one order of magnitude for asteroids up to about 200 m diameter. The two methods would yield comparable performance for asteroids larger than about 2 km

Bearing capacity of circular footings over rock mass by using axisymmetric quasi lower bound finite element limit analysis

The ultimate bearing capacity of a circular footing, placed over rock mass, is evaluated by using the lower bound theorem of the limit analysis in conjunction with finite elements and nonlinear optimization. The generalized Hoek-Brown (HB) failure criterion, but by keeping a constant value of the exponent, alpha = 0.5, was used. The failure criterion was smoothened both in the meridian and pi planes. The nonlinear optimization was carried out by employing an interior point method based on the logarithmic barrier function. The results for the obtained bearing capacity were presented in a non-dimensional form for different values of GSI, m(i), sigma(ci)/(gamma b) and q/sigma(ci). Failure patterns were also examined for a few cases. For validating the results, computations were also performed for a strip footing as well. The results obtained from the analysis compare well with the data reported in literature. Since the equilibrium conditions are precisely satisfied only at the centroids of the elements, not everywhere in the domain, the obtained lower bound solution will be approximate not true. (C) 2015 Elsevier Ltd. All rights reserved.

Reply to Discussion on ``Bearing capacity of circular footings over rock mass by using axisymmetric quasi lower bound finite element limit analysis'' Comput. Geotech. 70 (2015) 138-149]

A discussion has been provided for the comments raised by the discusser (Clausen, 2015)1] on the article recently published by the authors (Chakraborty and Kumar, 2015). The effect of exponent alpha for values of GSI approximately smaller than 30 becomes more critical. On the other hand, for greater values of GSI, the results obtained by the authors earlier remain primarily independent of alpha and can be easily used. (C) 2015 Elsevier Ltd. All rights reserved.

Energy levels, wavefunction compositions and electric dipole transitions in neutral Ca

A configuration-interaction approach, based on the use of B-spline basis sets combined with a model potential including monoelectronic and dielectronic core polarization effects, is employed to calculate term energies and wavefunctions for neutral Ca. Results are reported for singlet and triplet bound states, and some quasi-bound states above the lowest ionization limit, with angular momentum up to L = 4. Comparison with experiment and with other theoretical results shows that this method yields the most accurate energy values for neutral Ca obtained to date. Wavefunction compositions, necessary for labelling the levels, and the effects of semi-empirical polarization potentials on the wavefunctions are discussed, as are some recent identifications of doubly-excited states. It is shown that taking into account dielectronic core polarization changes the energies of the lowest terms in Ca significantly, in general by a few hundred cm(-1), the effect decreasing rapidly for the higher bound states. For Rydberg states with n approximate to 7 the accuracy of the results is often better than a few cm(-1). For series members (or perturbers) with a pronounced 3d character the error can reach 150 cm(-1). The wavefunctions are used to calculate oscillator strengths and lifetimes for a number of terms and these are compared with existing measurements. The agreement is good but points to a need for improved measurements.

Oscillation Criteria of Second-Order Quasi-Linear Neutral Delay Difference Equations

2000 Mathematics Subject Classification: 39A10.

Semi-empirical procedures for estimation of residual velocity and ballistic limit for impact on mild steel plates by projectiles

This paper deals with the development of simplified semi-empirical relations for the prediction of residual velocities of small calibre projectiles impacting on mild steel target plates, normally or at an angle, and the ballistic limits for such plates. It has been shown, for several impact cases for which test results on perforation of mild steel plates are available, that most of the existing semi-empirical relations which are applicable only to normal projectile impact do not yield satisfactory estimations of residual velocity. Furthermore, it is difficult to quantify some of the empirical parameters present in these relations for a given problem. With an eye towards simplicity and ease of use, two new regression-based relations employing standard material parameters have been discussed here for predicting residual velocity and ballistic limit for both normal and oblique impact. The latter expressions differ in terms of usage of quasi-static or strain rate-dependent average plate material strength. Residual velocities yielded by the present semi-empirical models compare well with the experimental results. Additionally, ballistic limits from these relations show close correlation with the corresponding finite element-based predictions.