Analysis of a finite volume method for a cross-diffusion model in population dynamics

The main goal of this paper is to propose a convergent finite volume method for a reactionâeuro"diffusion system with cross-diffusion. First, we sketch an existence proof for a class of cross-diffusion systems. Then the standard two-point finite volume fluxes are used in combination with a nonlinear positivity-preserving approximation of the cross-diffusion coefficients. Existence and uniqueness of the approximate solution are addressed, and it is also shown that the scheme converges to the corresponding weak solution for the studied model. Furthermore, we provide a stability analysis to study pattern-formation phenomena, and we perform two-dimensional numerical examples which exhibit formation of nonuniform spatial patterns. From the simulations it is also found that experimental rates of convergence are slightly below second order. The convergence proof uses two ingredients of interest for various applications, namely the discrete Sobolev embedding inequalities with general boundary conditions and a space-time $L^1$ compactness argument that mimics the compactness lemma due to Kruzhkov. The proofs of these results are given in the Appendix.

Cardiovascular effects of short term exposure to electric and magnetic fields of electricity power transmission

In the electrical industry the 50 Hz electric and magnetic fields are often higher than in the average working environment. The electric and magnetic fields can be studied by measuring or by calculatingthe fields in the environment. For example, the electric field under a 400 kV power line is 1 to 10 kV/m, and the magnetic flux density is 1 to 15 µT. Electricand magnetic fields of a power line induce a weak electric field and electric currents in the exposed body. The average current density in a human being standing under a 400 kV line is 1 to 2 mA/m2. The aim of this study is to find out thepossible effects of short term exposure to electric and magnetic fields of electricity power transmission on workers' health, in particular the cardiovascular effects. The study consists of two parts; Experiment I: influence on extrasystoles, and Experiment II: influence on heart rate. In Experiment I two groups, 26 voluntary men (Group 1) and 27 transmission-line workers (Group 2), were measured. Their electrocardiogram (ECG) was recorded with an ambulatory recorder both in and outside the field. In Group 1 the fields were 1.7 to 4.9 kV/m and 1.1 to 7.1 pT; in Group 2 they were 0.1 to 10.2 kV/m and 1.0 to 15.4 pT. In the ECG analysis the only significant observation was a decrease in the heart rate after field exposure (Group 1). The drop cannot be explained with the first measuring method. Therefore Experiment II was carried out. In Experiment II two groups were used; Group 1 (26 male volunteers) were measured in real field exposure, Group 2 (15 male volunteers) in "sham" fields. The subjects of Group 1 spent 1 h outside the field, then 1 h in the field under a 400 kV transmission line, and then again 1 h outside the field. Under the 400 kV linethe field strength varied from 3.5 to 4.3 kV/m, and from 1.4 to 6.6 pT. Group 2spent the entire test period (3 h) in a 33 kV outdoor testing station in a "sham" field. ECG, blood pressure, and electroencephalogram (EEG) were measured by ambulatory methods. Before and after the field exposure, the subjects performed some cardiovascular autonomic function tests. The analysis of the results (Experiments I and II) showed that extrasystoles or arrythmias were as frequent in the field (below 4 kV/m and 4 pT) as outside it. In Experiment II there was no decrease detected in the heart rate, and the systolic and diastolic blood pressure stayed nearly the same. No health effects were found in this study.

A novel approachfor assessing the fatigue strength of ultrasonic impact treated welded structures.

It is commonly observed that complex fabricated structures subject tofatigue loading fail at the welded joints. Some problems can be corrected by proper detail design but fatigue performance can also be improved using post-weld improvement methods. In general, improvement methods can be divided into two main groups: weld geometry modification methods and residual stress modification methods. The former remove weld toe defects and/or reduce the stress concentrationwhile the latter introduce compressive stress fields in the area where fatigue cracks are likely to initiate. Ultrasonic impact treatment (UIT) is a novel post-weld treatment method that influences both the residual stress distribution andimproves the local geometry of the weld. The structural fatigue strength of non-load carrying attachments in the as-welded condition has been experimentally compared to the structural fatigue strength of ultrasonic impact treated welds. Longitudinal attachment specimens made of two thicknesses of steel S355 J0 have been tested for determining the efficiency of ultrasonic impacttreatment. Treated welds were found to have about 50% greater structural fatigue strength, when the slope of the S-N-curve is three. High mean stress fatigue testing based on the Ohta-method decreased the degree of weld improvement only 19%. This indicated that the method could be also applied for large fabricated structures operating under high reactive residual stresses equilibrated within the volume of the structure. The thickness of specimens has no significant effect tothe structural fatigue strength. The fatigue class difference between 5 mm and 8 mm specimen was only 8%. It was hypothesized that the UIT method added a significant crack initiation period to the total fatigue life of the welded joints. Crack initiation life was estimated by a local strain approach. Material parameters were defined using a modified Uniform Material Law developed in Germany. Finite element analysis and X-ray diffraction were used to define, respectively, the stress concentration and mean stress. The theoretical fatigue life was found to have good accuracy comparing to experimental fatigue tests.The predictive behaviour of the local strain approach combined with the uniformmaterial law was excellent for the joint types and conditions studied in this work.

Theoretical and numerical study of thermomagnetic convection in magnetic fluids

In this thesis, the magnetic field control of convection instabilities and heat and mass transfer processesin magnetic fluids have been investigated by numerical simulations and theoretical considerations. Simulation models based on finite element and finite volume methods have been developed. In addition to standard conservation equations, themagnetic field inside the simulation domain is calculated from Maxwell equations and the necessary terms to take into account for the magnetic body force and magnetic dissipation have been added to the equations governing the fluid motion.Numerical simulations of magnetic fluid convection near the threshold supportedexperimental observations qualitatively. Near the onset of convection the competitive action of thermal and concentration density gradients leads to mostly spatiotemporally chaotic convection with oscillatory and travelling wave regimes, previously observed in binary mixtures and nematic liquid crystals. In many applications of magnetic fluids, the heat and mass transfer processes including the effects of external magnetic fields are of great importance. In addition to magnetic fluids, the concepts and the simulation models used in this study may be applied also to the studies of convective instabilities in ordinary fluids as well as in other binary mixtures and complex fluids.

Energy-momentum tensor for scalar fields coupled to the dilaton in two dimensions

We clarify some issues related to the evaluation of the mean value of the energy-momentum tensor for quantum scalar fields coupled to the dilaton field in two-dimensional gravity. Because of this coupling, the energy-momentum tensor for matter is not conserved and therefore it is not determined by the trace anomaly. We discuss different approximations for the calculation of the energy-momentum tensor and show how to obtain the correct amount of Hawking radiation. We also compute cosmological particle creation and quantum corrections to the Newtonian potential.

Technical Note: Dissolved organic matter fluorescence a finite mixture approach to deconvolve excitation-emission matrices

The analysis of the shape of excitation-emission matrices (EEMs) is a relevant tool for exploring the origin, transport and fate of dissolved organic matter (DOM) in aquatic ecosystems. Within this context, the decomposition of EEMs is acquiring a notable relevance. A simple mathematical algorithm that automatically deconvolves individual EEMs is described, creating new possibilities for the comparison of DOM fluorescence properties and EEMs that are very different from each other. A mixture model approach is adopted to decompose complex surfaces into sub-peaks. The laplacian operator and the Nelder-Mead optimisation algorithm are implemented to individuate and automatically locate potential peaks in the EEM landscape. The EEMs of a simple artificial mixture of fluorophores and DOM samples collected in a Mediterranean river are used to describe the model application and to illustrate a strategy that optimises the search for the optimal output.

Single acquisition electrical property mapping based on relative coil sensitivities: A proof-of-concept demonstration.

PURPOSE: All methods presented to date to map both conductivity and permittivity rely on multiple acquisitions to compute quantitatively the magnitude of radiofrequency transmit fields, B1+. In this work, we propose a method to compute both conductivity and permittivity based solely on relative receive coil sensitivities ( B1-) that can be obtained in one single measurement without the need to neither explicitly perform transmit/receive phase separation nor make assumptions regarding those phases. THEORY AND METHODS: To demonstrate the validity and the noise sensitivity of our method we used electromagnetic finite differences simulations of a 16-channel transceiver array. To experimentally validate our methodology at 7 Tesla, multi compartment phantom data was acquired using a standard 32-channel receive coil system and two-dimensional (2D) and 3D gradient echo acquisition. The reconstructed electric properties were correlated to those measured using dielectric probes. RESULTS: The method was demonstrated both in simulations and in phantom data with correlations to both the modeled and bench measurements being close to identity. The noise properties were modeled and understood. CONCLUSION: The proposed methodology allows to quantitatively determine the electrical properties of a sample using any MR contrast, with the only constraint being the need to have 4 or more receive coils and high SNR. Magn Reson Med, 2014. © 2014 Wiley Periodicals, Inc.

Thomas-Fermi theory for atomic nuclei revisited

The recently developed semiclassical variational Wigner-Kirkwood (VWK) approach is applied to finite nuclei using external potentials and self-consistent mean fields derived from Skyrme inter-actions and from relativistic mean field theory. VWK consist s of the Thomas-Fermi part plus a pure, perturbative h 2 correction. In external potentials, VWK passes through the average of the quantal values of the accumulated level density and total en energy as a function of the Fermi energy. However, there is a problem of overbinding when the energy per particle is displayed as a function of the particle number. The situation is analyzed comparing spherical and deformed harmonic oscillator potentials. In the self-consistent case, we show for Skyrme forces that VWK binding energies are very close to those obtained from extended Thomas-Fermi functionals of h 4 order, pointing to the rapid convergence of the VWK theory. This satisfying result, however, does not cure the overbinding problem, i.e., the semiclassical energies show more binding than they should. This feature is more pronounced in the case of Skyrme forces than with the relativistic mean field approach. However, even in the latter case the shell correction energy for e.g.208 Pb turns out to be only ∼ −6 MeV what is about a factor two or three off the generally accepted value. As an adhoc remedy, increasing the kinetic energy by 2.5%, leads to shell correction energies well acceptable throughout the periodic table. The general importance of the present studies for other finite Fermi systems, self-bound or in external potentials, is pointed out.

Behavior of propagating and evanescent components in azimuthally polarized non-paraxial fields

The contribution of the propagating and the evanescent waves associated with freely propagating non-paraxial light fields whose transverse component is azimuthally polarized at some plane is investigated. Analytic expressions are derived for describing both the spatial shape and the relative weight of the propagating and the evanescent components integrated over the transverse plane. The analysis is carried out within the framework of the plane-wave angular spectrum approach. These results are used to illustrate the behavior of a kind of donut-like beams with transverse azimuthal polarization at some plane.

Ecuaciones de recurrencia estocásticas en el cálculo de la prima de reaseguro Finite Risk

RESUMEN: El objetivo de este trabajo es calcular el importe de la prima pura periódica que debe cobrar el reasegurador a la cedente en un reaseguro finite risk en ambiente financiero estocástico. El problema de la convolución de las diferentes variables aleatorias que intervienen en el cálculo de la prima lo hemos solucionado simulando, por Monte-Carlo, trayectorias de siniestralidad para el reasegurador aplicando posteriormente, en cada trayectoria simulada, los criterios de decisión financieros, esperanza, varianza y desviación. En los criterios de la varianza y de la desviación proponemos utilizar una ecuación de recurrencia estocástica para evitar el problema de la dependencia que existe entre los factores de capitalización estocásticos, obteniendo la prima de reaseguro en función del nivel de aversión al riesgo del reasegurador y de la volatilidad del tipo de interés. Palabras clave: Finite risk, ambiente estocástico, ecuación de recurrencia, simulación de Monte-Carlo, prima pura periódica.

Finite Element Analysis of Electrically Excited Quartz Tuning Fork Devices

Quartz Tuning Fork (QTF)-based Scanning Probe Microscopy (SPM) is an important field of research. A suitable model for the QTF is important to obtain quantitative measurements with these devices. Analytical models have the limitation of being based on the double cantilever configuration. In this paper, we present an electromechanical finite element model of the QTF electrically excited with two free prongs. The model goes beyond the state-of-the-art of numerical simulations currently found in the literature for this QTF configuration. We present the first numerical analysis of both the electrical and mechanical behavior of QTF devices. Experimental measurements obtained with 10 units of the same model of QTF validate the finite element model with a good agreement.

Synthesis of highly focused fields with circular polarization at any transverse plane

We develop a method for generating focused vector beams with circular polarization at any transverse plane. Based on the Richards-Wolf vector model, we derive analytical expressions to describe the propagation of these set of beams near the focal area. Since the polarization and the amplitude of the input beam are not uniform, an interferometric system capable of generating spatially-variant polarized beams has to be used. In particular, this wavefront is manipulated by means of spatial light modulators displaying computer generated holograms and subsequently focused using a high numerical aperture objective lens. Experimental results using a NA=0.85 system are provided: irradiance and Stokes images of the focused field at different planes near the focal plane are presented and compared with those obtained by numerical simulation.

On the longitudinal component of paraxial fields

The analysis of paraxial Gaussian beams features in most undergraduate courses in laser physics, advanced optics and photonics. These beams provide a simple model of the field generated in the resonant cavities of lasers, thus constituting a basic element for understanding laser theory. Usually, uniformly polarized beams are considered in the analytical calculations, with the electric field vibrating at normal planes to the propagation direction. However, such paraxial fields do not verify the Maxwell equations. In this paper we discuss how to overcome this apparent contradiction and evaluate the longitudinal component that any paraxial Gaussian beam should exhibit. Despite the fact that the assumption of a purely transverse paraxial field is useful and accurate, the inclusion of the above issue in the program helps students to clarify the importance of the electromagnetic nature of light, thus providing a more complete understanding of the paraxial approach.

Transverse and longitudinal components of the propagating and evanescent waves associated to radially-polarized nonparaxial fields

A comparison is established between the contributions of transverse and longitudinal components of both the propagating and the evanescent waves associated to freely propagating radially polarized nonparaxial beams. Attention is focused on those fields that remain radially polarized upon propagation. In terms of the plane-wave angular spectrum of these fields, analytical expressions are given for determining both the spatial shape of the above components and their relative weight integrated over the whole transverse plane. The results are applied to two kinds of doughnut-like beams with radial polarization, and we compare the behavior of such fields at two transverse planes.