A moving-mesh finite element method and its application to the numerical solution of phase-change problems

A distributed Lagrangian moving-mesh finite element method is applied to problems involving changes of phase. The algorithm uses a distributed conservation principle to determine nodal mesh velocities, which are then used to move the nodes. The nodal values are obtained from an ALE (Arbitrary Lagrangian-Eulerian) equation, which represents a generalization of the original algorithm presented in Applied Numerical Mathematics, 54:450--469 (2005). Having described the details of the generalized algorithm it is validated on two test cases from the original paper and is then applied to one-phase and, for the first time, two-phase Stefan problems in one and two space dimensions, paying particular attention to the implementation of the interface boundary conditions. Results are presented to demonstrate the accuracy and the effectiveness of the method, including comparisons against analytical solutions where available.

Numerical modeling of inhaled charged aerosol deposition in human airways

A new numerical modeling of inhaled charge aerosol has been developed based on a modified Weibel's model. Both the velocity profiles (slug and parabolic flows) and the particle distributions (uniform and parabolic distributions) have been considered. Inhaled particles are modeled as a dilute dispersed phase flow in which the particle motion is controlled by fluid force and external forces acting on particles. This numerical study extends the previous numerical studies by considering both space- and image-charge forces. Because of the complex computation of interacting forces due to space-charge effect, the particle-mesh (PM) method is selected to calculate these forces. In the PM technique, the charges of all particles are assigned to the space-charge field mesh, for calculating charge density. The Poisson's equation of the electrostatic potential is then solved, and the electrostatic force acting on individual particle is interpolated. It is assumed that there is no effect of humidity on charged particles. The results show that many significant factors also affect the deposition, such as the volume of particle cloud, the velocity profile and the particle distribution. This study allows a better understanding of electrostatic mechanism of aerosol transport and deposition in human airways.

The distribution of solar wind speeds during solar minimum: calibration for numerical solar wind modeling constraints on the source of the slow solar wind

It took the solar polar passage of Ulysses in the early 1990s to establish the global structure of the solar wind speed during solar minimum. However, it remains unclear if the solar wind is composed of two distinct populations of solar wind from different sources (e.g., closed loops which open up to produce the slow solar wind) or if the fast and slow solar wind rely on the superradial expansion of the magnetic field to account for the observed solar wind speed variation. We investigate the solar wind in the inner corona using the Wang-Sheeley-Arge (WSA) coronal model incorporating a new empirical magnetic topology–velocity relationship calibrated for use at 0.1 AU. In this study the empirical solar wind speed relationship was determined by using Helios perihelion observations, along with results from Riley et al. (2003) and Schwadron et al. (2005) as constraints. The new relationship was tested by using it to drive the ENLIL 3-D MHD solar wind model and obtain solar wind parameters at Earth (1.0 AU) and Ulysses (1.4 AU). The improvements in speed, its variability, and the occurrence of high-speed enhancements provide confidence that the new velocity relationship better determines the solar wind speed in the outer corona (0.1 AU). An analysis of this improved velocity field within the WSA model suggests the existence of two distinct mechanisms of the solar wind generation, one for fast and one for slow solar wind, implying that a combination of present theories may be necessary to explain solar wind observations.

A method for merging flow-dependent forecast error statistics from an ensemble with static statistics for use in high resolution variational data assimilation

The background error covariance matrix, B, is often used in variational data assimilation for numerical weather prediction as a static and hence poor approximation to the fully dynamic forecast error covariance matrix, Pf. In this paper the concept of an Ensemble Reduced Rank Kalman Filter (EnRRKF) is outlined. In the EnRRKF the forecast error statistics in a subspace defined by an ensemble of states forecast by the dynamic model are found. These statistics are merged in a formal way with the static statistics, which apply in the remainder of the space. The combined statistics may then be used in a variational data assimilation setting. It is hoped that the nonlinear error growth of small-scale weather systems will be accurately captured by the EnRRKF, to produce accurate analyses and ultimately improved forecasts of extreme events.

A comparison of extra-tropical cyclones in recent re-analyses; ERA-Interim, NASA MERRA, NCEP CFSR and JRA-25

Extra-tropical cyclones are identified and compared using data from four recent re-analyses for the winter periods in both hemispheres. Results show the largest differences occur between the older lower resolution JRA25 re-analysis when compared with the newer high resolution re-analyses, in particular in the Southern Hemisphere (SH). Spatial differences between the newest re-analyses are small in both hemispheres and generally not significant except some common regions associated with cyclogenesis close to orography. Intensities are generally related to spatial resolution except NASA-MERRA which has larger intensities for several different measures. Matching storms between re-analyses shows the number matched between ERA-Interim and the other re-analyses are similar in the Northern Hemisphere (NH). In the SH the number matched between JRA25 and ERA-Interim is lower than in the NH, but for NASA-MERRA and NCEP-CFSR the number matched is similar to the NH. The mean separation of the identically same cyclones is typically less than 20 geodesic in both hemispheres for the latest re-analyses, whereas JRA25 compared with the other re-analyses has a broader distribution in the SH indicating greater uncertainty. The instantaneous intensity differences for matched storms shows narrow distributions for pressure while for winds and vorticity the distributions are much broader indicating larger uncertainty typical of smaller scale fields. Composite cyclone diagnostics show that cyclones are very similar between the re-analyses, with differences being related to the intensities, consistent with the intensity results. Overall, results show NH cyclones correspond well between re-analyses, with a significant improvement in the SH for the latest re-analyses, indicating a convergence between re-analyses for cyclone properties.

Tenant mix variety in regional shopping centres: some UK empirical analyses

The variety and quality of the tenant mix within a shopping centre is a key concern in shopping centre management. Tenant mix determines the extent of externalities between outlets in the centre, helps establish the image of the centre and, as a result, determines the attractiveness of the centre for consumers. This then translates into sales and rents. However, the management of tenant mix has largely been based on perceived “optimum” arrangements and industry rules of thumb. This paper attempts to model the impact of tenant mix on the rent paid by retailers in larger UK shopping centres and, hence, the returns made by shopping centre landlords. It extends work on shopping centre rent determination (see Working Paper 10/03) utilising a database of 148 regional shopping centres in the UK, with detailed data for over 1900 tenants. Econometric models test the relationship between rental levels and the levels of retail concentration and diversity, while controlling for a range of continuous and qualitative characteristics of each tenant, each retail product, and each shopping centre. Factor analysis is then used to extract the core retail and service categories from the tenant lists of the 148 shopping centres. The factor scores from these core retailer factors are then tested against rent payable. The results from the empirical analysis allow us to generate some clear analytical and empirical implications for optimal retail management.

Sensory and instrumental analyses of volatiles generated during the extrusion cooking of oat flours

Three batches of oats were extruded under four combinations of process temperature (150 or 180 °C) and process moisture (14.5 and 18%). Two of the extrudates were evaluated by a sensory panel, and three were analyzed by GC-MS. Maillard reaction products, such as pyrazines, pyrroles, furans, and sulfur-containing compounds, were found in the most severely processed extrudates (high-temperature, low-moisture). These extrudates were also described by the assessors as having toasted cereal attributes. Lipid degradation products, such as alkanals, 2-alkenals, and 2,4-alkadienals, were found at much higher levels in the extrudates of the oat flour that had been debranned. It contained lower protein and fiber levels than the others and showed increased lipase activity. Extrudates from these samples also had significantly lower levels of Maillard reaction products that correlated, in the sensory analysis, with terms such as stale oil and oatmeal. Linoleic acid was added to a fourth oat flour to simulate the result of increased lipase activity, and GC-MS analysis showed both an increase in lipid degradation products and a decrease in Maillard reaction products.

A study of the arrival over the United Kingdom in April 2010 of the Eyjafjallajökull ash cloud using ground-based lidar and numerical simulations

We make a qualitative and quantitative comparison of numericalsimulations of the ashcloud generated by the eruption of Eyjafjallajökull in April2010 with ground-basedlidar measurements at Exeter and Cardington in southern England. The numericalsimulations are performed using the Met Office’s dispersion model, NAME (Numerical Atmospheric-dispersion Modelling Environment). The results show that NAME captures many of the features of the observed ashcloud. The comparison enables us to estimate the fraction of material which survives the near-source fallout processes and enters into the distal plume. A number of simulations are performed which show that both the structure of the ashcloudover southern England and the concentration of ash within it are particularly sensitive to the height of the eruption column (and the consequent estimated mass emission rate), to the shape of the vertical source profile and the level of prescribed ‘turbulent diffusion’ (representing the mixing by the unresolved eddies) in the free troposphere with less sensitivity to the timing of the start of the eruption and the sedimentation of particulates in the distal plume.

Numerical convergence of the Block-Maxima approach to the generalized extreme value distribution

In this paper we perform an analytical and numerical study of Extreme Value distributions in discrete dynamical systems. In this setting, recent works have shown how to get a statistics of extremes in agreement with the classical Extreme Value Theory. We pursue these investigations by giving analytical expressions of Extreme Value distribution parameters for maps that have an absolutely continuous invariant measure. We compare these analytical results with numerical experiments in which we study the convergence to limiting distributions using the so called block-maxima approach, pointing out in which cases we obtain robust estimation of parameters. In regular maps for which mixing properties do not hold, we show that the fitting procedure to the classical Extreme Value Distribution fails, as expected. However, we obtain an empirical distribution that can be explained starting from a different observable function for which Nicolis et al. (Phys. Rev. Lett. 97(21): 210602, 2006) have found analytical results.

Estimation of systematic error in an equatorial ocean model using data assimilation

Assimilation of temperature observations into an ocean model near the equator often results in a dynamically unbalanced state with unrealistic overturning circulations. The way in which these circulations arise from systematic errors in the model or its forcing is discussed. A scheme is proposed, based on the theory of state augmentation, which uses the departures of the model state from the observations to update slowly evolving bias fields. Results are summarized from an experiment applying this bias correction scheme to an ocean general circulation model. They show that the method produces more balanced analyses and a better fit to the temperature observations.