Implementation of an interior point source in the ultra weak variational formulation through source extraction

The Ultra Weak Variational Formulation (UWVF) is a powerful numerical method for the approximation of acoustic, elastic and electromagnetic waves in the time-harmonic regime. The use of Trefftz-type basis functions incorporates the known wave-like behaviour of the solution in the discrete space, allowing large reductions in the required number of degrees of freedom for a given accuracy, when compared to standard finite element methods. However, the UWVF is not well disposed to the accurate approximation of singular sources in the interior of the computational domain. We propose an adjustment to the UWVF for seismic imaging applications, which we call the Source Extraction UWVF. Differing fields are solved for in subdomains around the source, and matched on the inter-domain boundaries. Numerical results are presented for a domain of constant wavenumber and for a domain of varying sound speed in a model used for seismic imaging.

A Lagrangian analysis of ice-supersaturated air over the North Atlantic

Understanding the nature of air parcels that exhibit ice-supersaturation is important because they are the regions of potential formation of both cirrus and aircraft contrails, which affect the radiation balance. Ice-supersaturated air parcels in the upper troposphere and lower stratosphere over the North Atlantic are investigated using Lagrangian trajectories. The trajectory calculations use ERA-Interim data for three winter and three summer seasons, resulting in approximately 200,000 trajectories with ice-supersaturation for each season. For both summer and winter, the median duration of ice-supersaturation along a trajectory is less than 6 hours. 5% of air which becomes ice-supersaturated in the troposphere, and 23% of air which becomes ice-supersaturated in the stratosphere will remain ice-supersaturated for at least 24 hours. Weighting the ice-supersaturation duration with the observed frequency indicates the likely overall importance of the longer duration ice-supersaturated trajectories. Ice-supersaturated air parcels typically experience a decrease in moisture content while ice-supersaturated, suggesting that cirrus clouds eventually form in the majority of such air. A comparison is made between short-lived (less than 24 h) and long-lived (greater than 24 h) ice-supersaturated air flows. For both air flows, ice-supersaturation occurs around the northernmost part of the trajectory. Short-lived ice-supersaturated air flows show no significant differences in speed or direction of movement to subsaturated air parcels. However, long-lived ice-supersaturated air occurs in slower moving air flows, which implies that they are not associated with the fastest moving air through a jet stream.

A laminated polymer film formulation for enteric delivery of live vaccine and probiotic bacteria

Live bacterial cells (LBC) are administered orally as attenuated vaccines, to deliver biopharmaceutical agents, and as probiotics to improve gastrointestinal health. However, LBC present unique formulation challenges and must survive gastrointestinal antimicrobial defenses including gastric acid after administration. We present a simple new formulation concept, termed Polymer Film Laminate (PFL). LBC are ambient dried onto cast acid-resistant enteric polymer films that are then laminated together to produce a solid oral dosage form. LBC of a model live bacterial vaccine and a probiotic were dried directly onto a cast film of enteric polymer. The effectiveness at protecting dried cells in a simulated gastric fluid (pH 2.0) depended on the composition of enteric polymer film used, with a blend of ethylcellulose plus Eudragit L100 55 providing greater protection from acid than Eudragit alone. However, although PFL made from blended polymers films completely released low molecular weight dye into intestinal conditions (pH 7.0), they failed to release LBC. In contrast, PFL made from Eudragit alone successfully protected dried probiotic or vaccine LBC from simulated gastric fluid for 2h, and subsequently released all viable cells within 60min of transfer into simulated intestinal fluid. Release kinetics could be controlled by modifying the lamination method.

Formulation and implementation of a residual-mean ocean circulation model

A parameterization of mesoscale eddies in coarse-resolution ocean general circulation models (GCM) is formulated and implemented using a residual-mean formalism. In that framework, mean buoyancy is advected by the residual velocity (the sum of the Eulerian and eddy-induced velocities) and modified by a residual flux which accounts for the diabatic effects of mesoscale eddies. The residual velocity is obtained by stepping forward a residual-mean momentum equation in which eddy stresses appear as forcing terms. Study of the spatial distribution of eddy stresses, derived by using them as control parameters to ‘‘fit’’ the residual-mean model to observations, supports the idea that eddy stresses can be likened to a vertical down-gradient flux of momentum with a coefficient which is constant in the vertical. The residual eddy flux is set to zero in the ocean interior, where mesoscale eddies are assumed to be quasi-adiabatic, but is parameterized by a horizontal down-gradient diffusivity near the surface where eddies develop a diabatic component as they stir properties horizontally across steep isopycnals. The residual-mean model is implemented and tested in the MIT general circulation model. It is shown that the resulting model (1) has a climatology that is superior to that obtained using the Gent and McWilliams parameterization scheme with a spatially uniform diffusivity and (2) allows one to significantly reduce the (spurious) horizontal viscosity used in coarse resolution GCMs.

Application of complex extreme learning machine to multiclass classification problems with high dimensionality: a THz spectra classification problem

We extend extreme learning machine (ELM) classifiers to complex Reproducing Kernel Hilbert Spaces (RKHS) where the input/output variables as well as the optimization variables are complex-valued. A new family of classifiers, called complex-valued ELM (CELM) suitable for complex-valued multiple-input–multiple-output processing is introduced. In the proposed method, the associated Lagrangian is computed using induced RKHS kernels, adopting a Wirtinger calculus approach formulated as a constrained optimization problem similarly to the conventional ELM classifier formulation. When training the CELM, the Karush–Khun–Tuker (KKT) theorem is used to solve the dual optimization problem that consists of satisfying simultaneously smallest training error as well as smallest norm of output weights criteria. The proposed formulation also addresses aspects of quaternary classification within a Clifford algebra context. For 2D complex-valued inputs, user-defined complex-coupled hyper-planes divide the classifier input space into four partitions. For 3D complex-valued inputs, the formulation generates three pairs of complex-coupled hyper-planes through orthogonal projections. The six hyper-planes then divide the 3D space into eight partitions. It is shown that the CELM problem formulation is equivalent to solving six real-valued ELM tasks, which are induced by projecting the chosen complex kernel across the different user-defined coordinate planes. A classification example of powdered samples on the basis of their terahertz spectral signatures is used to demonstrate the advantages of the CELM classifiers compared to their SVM counterparts. The proposed classifiers retain the advantages of their ELM counterparts, in that they can perform multiclass classification with lower computational complexity than SVM classifiers. Furthermore, because of their ability to perform classification tasks fast, the proposed formulations are of interest to real-time applications.

A parallel formulation for the simulation of a generic branch predictor

A parallel formulation for the simulation of a branch prediction algorithm is presented. This parallel formulation identifies independent tasks in the algorithm which can be executed concurrently. The parallel implementation is based on the multithreading model and two parallel programming platforms: pthreads and Cilk++. Improvement in execution performance by up to 7 times is observed for a generic 2-bit predictor in a 12-core multiprocessor system.

Effects of orange juice formulation on prebiotic functionality using an in vitro colonic model sytem

A three-stage continuous fermentative colonic model system was used to monitor in vitro the effect of different orange juice formulations on prebiotic activity. Three different juices with and without Bimuno, a GOS mixture containing galactooligosaccharides (B-GOS) were assessed in terms of their ability to induce a bifidogenic microbiota. The recipe development was based on incorporating 2.75g B-GOS into a 250 ml serving of juice (65°Brix of concentrate juice). Alongside the production of B-GOS juice, a control juice - orange juice without any additional Bimuno and a positive control juice, containing all the components of Bimuno (glucose, galactose and lactose) in the same relative proportions with the exception of B-GOS were developed. Ion Exchange Chromotography analysis was used to test the maintenance of bimuno components after the production process. Data showed that sterilisation had no significant effect on concentration of B-GOS and simple sugars. The three juice formulations were digested under conditions resembling the gastric and small intestinal environments. Main bacterial groups of the faecal microbiota were evaluated throughout the colonic model study using 16S rRNA-based fluorescence in situ hybridization (FISH). Potential effects of supplementation of the juices on microbial metabolism were studied measuring short chain fatty acids (SCFAs) using gas chromatography. Furthermore, B-GOS juices showed positive modulations of the microbiota composition and metabolic activity. In particular, numbers of faecal bifidobacteria and lactobacilli were significantly higher when B-GOS juice was fermented compared to controls. Furthermore, fermentation of B-GOS juice resulted in an increase in Roseburia subcluster and concomitantly increased butyrate production, which is of potential benefit to the host. In conclusion, this study has shown B-GOS within orange juice can have a beneficial effect on the fecal microbiota.

Optimization of a sponge cake formulation with inulin as fat replacer: structure, physicochemical, and sensory properties

The effects of several fat replacement levels (0%, 35%, 50%, 70%, and 100%) by inulin in sponge cake microstructure and physicochemical properties were studied. Oil substitution for inulin decreased significantly (P < 0.05) batter viscosity, giving heterogeneous bubbles size distributions as it was observed by light microscopy. Using confocal laser scanning microscopy the fat was observed to be located at the bubbles’ interface, enabling an optimum crumb cake structure development during baking. Cryo-SEM micrographs of cake crumbs showed a continuous matrix with embedded starch granules and coated with oil; when fat replacement levels increased, starch granules appeared as detached structures. Cakes with fat replacement up to 70% had a high crumb air cell values; they were softer and rated as acceptable by an untrained sensory panel (n = 51). So, the reformulation of a standard sponge cake recipe to obtain a new product with additional health benefits and accepted by consumers is achieved.

Optimisation of oat milk formulation to obtain fermented derivatives by using probiotic Lactobacillus reuteri microorganisms

Functional advantages of probiotics combined with interesting composition of oat were considered as an alternative to dairy products. In this study, fermentation of oat milk with Lactobacillus reuteri and Streptococcus thermophilus was analysed to develop a new probiotic product. Central composite design with response surface methodology was used to analyse the effect of different factors (glucose, fructose, inulin and starters) on the probiotic population in the product. Optimised formulation was characterised throughout storage time at 4 ℃ in terms of pH, acidity, β-glucan and oligosaccharides contents, colour and rheological behaviour. All formulations studied were adequate to produce fermented foods and minimum dose of each factor was considered as optimum. The selected formulation allowed starters survival above 107/cfu ml to be considered as a functional food and was maintained during the 28 days controlled. β-glucans remained in the final product with a positive effect on viscosity. Therefore, a new probiotic non-dairy milk was successfully developed in which high probiotic survivals were assured throughout the typical yoghurt-like shelf life.

Mechanisms underlying temperature extremes in Iberia: a Lagrangian perspective

The mechanisms underlying the occurrence of temperature extremes in Iberia are analysed considering a Lagrangian perspective of the atmospheric flow, using 6-hourly ERA-Interim reanalysis data for the years 1979–2012. Daily 2-m minimum temperatures below the 1st percentile and 2-m maximum temperatures above the 99th percentile at each grid point over Iberia are selected separately for winter and summer. Four categories of extremes are analysed using 10-d backward trajectories initialized at the extreme temperature grid points close to the surface: winter cold (WCE) and warm extremes (WWE), and summer cold (SCE) and warm extremes (SWE). Air masses leading to temperature extremes are first transported from the North Atlantic towards Europe for all categories. While there is a clear relation to large-scale circulation patterns in winter, the Iberian thermal low is important in summer. Along the trajectories, air mass characteristics are significantly modified through adiabatic warming (air parcel descent), upper-air radiative cooling and near-surface warming (surface heat fluxes and radiation). High residence times over continental areas, such as over northern-central Europe for WCE and, to a lesser extent, over Iberia for SWE, significantly enhance these air mass modifications. Near-surface diabatic warming is particularly striking for SWE. WCE and SWE are responsible for the most extreme conditions in a given year. For WWE and SCE, strong temperature advection associated with important meridional air mass transports are the main driving mechanisms, accompanied by comparatively minor changes in the air mass properties. These results permit a better understanding of mechanisms leading to temperature extremes in Iberia.

Data Assimilation and Inverse Methods in Terms of a Probabilistic Formulation

The weak-constraint inverse for nonlinear dynamical models is discussed and derived in terms of a probabilistic formulation. The well-known result that for Gaussian error statistics the minimum of the weak-constraint inverse is equal to the maximum-likelihood estimate is rederived. Then several methods based on ensemble statistics that can be used to find the smoother (as opposed to the filter) solution are introduced and compared to traditional methods. A strong point of the new methods is that they avoid the integration of adjoint equations, which is a complex task for real oceanographic or atmospheric applications. they also avoid iterative searches in a Hilbert space, and error estimates can be obtained without much additional computational effort. the feasibility of the new methods is illustrated in a two-layer quasigeostrophic model.