A probabilistic energy model for non-domestic building sectors applied to analysis of school buildings in greater London

The diversity of non-domestic buildings at urban scale poses a number of difficulties to develop models for large scale analysis of the stock. This research proposes a probabilistic, engineering-based, bottom-up model to address these issues. In a recent study we classified London's non-domestic buildings based on the service they provide, such as offices, retail premise, and schools, and proposed the creation of one probabilistic representational model per building type. This paper investigates techniques for the development of such models. The representational model is a statistical surrogate of a dynamic energy simulation (ES) model. We first identify the main parameters affecting energy consumption in a particular building sector/type by using sampling-based global sensitivity analysis methods, and then generate statistical surrogate models of the dynamic ES model within the dominant model parameters. Given a sample of actual energy consumption for that sector, we use the surrogate model to infer the distribution of model parameters by inverse analysis. The inferred distributions of input parameters are able to quantify the relative benefits of alternative energy saving measures on an entire building sector with requisite quantification of uncertainties. Secondary school buildings are used for illustrating the application of this probabilistic method. © 2012 Elsevier B.V. All rights reserved.

Material model of lung parenchyma based on living precision-cut lung slice testing.

We describe a novel constitutive model of lung parenchyma, which can be used for continuum mechanics based predictive simulations. To develop this model, we experimentally determined the nonlinear material behavior of rat lung parenchyma. This was achieved via uni-axial tension tests on living precision-cut rat lung slices. The resulting force-displacement curves were then used as inputs for an inverse analysis. The Levenberg-Marquardt algorithm was utilized to optimize the material parameters of combinations and recombinations of established strain-energy density functions (SEFs). Comparing the best-fits of the tested SEFs we found Wpar = 4.1 kPa(I1-3)2 + 20.7 kPa(I1 - 3)3 + 4.1 kPa(-2 ln J + J2 - 1) to be the optimal constitutive model. This SEF consists of three summands: the first can be interpreted as the contribution of the elastin fibers and the ground substance, the second as the contribution of the collagen fibers while the third controls the volumetric change. The presented approach will help to model the behavior of the pulmonary parenchyma and to quantify the strains and stresses during ventilation.

Energy use of buildings at urban scale: A case study of london school buildings

The diversity of non-domestic buildings at urban scale poses a number of difficulties to develop building stock models. This research proposes an engineering-based bottom-up stock model in a probabilistic manner to address these issues. School buildings are used for illustrating the application of this probabilistic method. Two sampling-based global sensitivity methods are used to identify key factors affecting building energy performance. The sensitivity analysis methods can also create statistical regression models for inverse analysis, which are used to estimate input information for building stock energy models. The effects of different energy saving measures are analysed by changing these building stock input distributions.

Wave propagation under vertically excited surface foundation

Recently developed equipment allows measurement of the shear modulus of soil in situ as a function of level of strain. In these field experiments, the excitation is applied on the ground surface using large scale shakers, and the response of the soil deposit is recorded through embedded receivers. The focus of this paper is on the simulation of signals which would be recorded at the receiver locations in idealized conditions to provide guidelines on the interpretation of field measurements. Discrete and finite element methods are employed to model one dimensional and three dimensional geometries, respectively, under various lateral boundary conditions. When the first times of arrival are detected by receivers under the vertical impulse, they coincide with the arrival of the P wave, related to the constrained modulus of the material, regardless of lateral boundary conditions. If one considers, on the other hand, phase differences between the motions at two receivers the picture is far more complicated and one would obtain propagation velocities, function of frequency and depth, which do not correspond to either the constrained modulus or Young's modulus. It is thus necessary to apply some care when interpreting the data from field tests based on vertical steady state vibrations. The use of inverse analysis can be considered as a way of extracting the shear modulus of soil from the field test measurements. © 2008 ASCE.

Application of the fixed point method for solution in time stepping finite element analysis using the inverse vector Jiles-Atherton model

An implementation of the inverse vector Jiles-Atherton model for the solution of non-linear hysteretic finite element problems is presented. The implementation applies the fixed point method with differential reluctivity values obtained from the Jiles-Atherton model. Differential reluctivities are usually computed using numerical differentiation, which is ill-posed and amplifies small perturbations causing large sudden increases or decreases of differential reluctivity values, which may cause numerical problems. A rule based algorithm for conditioning differential reluctivity values is presented. Unwanted perturbations on the computed differential reluctivity values are eliminated or reduced with the aim to guarantee convergence. Details of the algorithm are presented together with an evaluation of the algorithm by a numerical example. The algorithm is shown to guarantee convergence, although the rate of convergence depends on the choice of algorithm parameters. © 2011 IEEE.

Finite element analysis of compressor disc heat transfer using the transient conduction method with sparse boundary collocation points

Surface temperature measurements from two discs of a gas turbine compressor rig are used as boundary conditions for the transient conduction solution (inverse heat transfer analysis). The disc geometry is complex, and so the finite element method is used. There are often large radial temperature gradients on the discs, and the equations are therefore solved taking into account the dependence of thermal conductivity on temperature. The solution technique also makes use of a multigrid algorithm to reduce the solution time. This is particularly important since a large amount of data must be analyzed to obtain correlations of the heat transfer. The finite element grid is also used for a network analysis to calculate the radiant heat transfer in the cavity formed between the two compressor discs. The work discussed here proved particularly challenging as the disc temperatures were only measured at four different radial locations. Four methods of surface temperature interpolation are examined, together with their effect on the local heat fluxes. It is found that the choice of interpolation method depends on the available number of data points. Bessel interpolation gives the best results for four data points, whereas cubic splines are preferred when there are considerably more data points. The results from the analysis of the compressor rig data show that the heat transfer near the disc inner radius appears to be influenced by the central throughflow. However, for larger radii, the heat transfer from the discs and peripheral shroud is found to be consistent with that of a buoyancy-induced flow.