Automatic extraction of Reduced order models from CFD simulations for building energy modeling

Thermal comfort is defined as “that condition of mind which expresses satisfaction with the thermal environment’ [1] [2]. Field studies have been completed in order to establish the governing conditions for thermal comfort [3]. These studies showed that the internal climate of a room was the strongest factor in establishing thermal comfort. Direct manipulation of the internal climate is necessary to retain an acceptable level of thermal comfort. In order for Building Energy Management Systems (BEMS) strategies to be efficiently utilised it is necessary to have the ability to predict the effect that activating a heating/cooling source (radiators, windows and doors) will have on the room. The numerical modelling of the domain can be challenging due to necessity to capture temperature stratification and/or different heat sources (radiators, computers and human beings). Computational Fluid Dynamic (CFD) models are usually utilised for this function because they provide the level of details required. Although they provide the necessary level of accuracy these models tend to be highly computationally expensive especially when transient behaviour needs to be analysed. Consequently they cannot be integrated in BEMS. This paper presents and describes validation of a CFD-ROM method for real-time simulations of building thermal performance. The CFD-ROM method involves the automatic extraction and solution of reduced order models (ROMs) from validated CFD simulations. The test case used in this work is a room of the Environmental Research Institute (ERI) Building at the University College Cork (UCC). ROMs have shown that they are sufficiently accurate with a total error of less than 1% and successfully retain a satisfactory representation of the phenomena modelled. The number of zones in a ROM defines the size and complexity of that ROM. It has been observed that ROMs with a higher number of zones produce more accurate results. As each ROM has a time to solution of less than 20 seconds they can be integrated into the BEMS of a building which opens the potential to real time physics based building energy modelling.

REDUCED ORDER MODELLING OF THE THERMAL BEHAVIOUR OF AN OFFICE SPACE

Reduced Order Models (ROMs) have proven to be a valid and efficient approach to model the thermal behaviour of building zones. The main issues associated with the use of zonal/lumped models are how to (1) divide the domain (lumps) and (2) evaluate the pa- rameters which characterise the lump-to-lump exchange of energy and momentum. The object of this research is to develop a methodology for the generation of ROMs from CFD models. The lumps of the ROM and their average property values are automatically ex- tracted from the CFD models through user defined constraints. This methodology has been applied to validated CFD models of a zone of the Environmental Research Insti- tute (ERI) Building in University College Cork (UCC). The ROM predicts temperature distribution in the domain with an average error lower than 2%. It is computationally efficient with an execution time of 3.45 seconds. Future steps in this research will be the development of the procedure to automatically extract the parameters which define lump-to-lump energy and momentum exchange. At the moment these parameters are evaluated through the minimisation of a cost function. The ROMs will also be utilised to predict the transient thermal behaviour of the building zone.

Automatic extraction of reduced-order models from CFD simulations for building energy modelling

Accurate modelling of the internal climate of buildings is essential if Building Energy Management Systems (BEMS) are to efficiently maintain adequate thermal comfort. Computational fluid dynamics (CFD) models are usually utilised to predict internal climate. Nevertheless CFD models, although providing the necessary level of accuracy, are highly computationally expensive, and cannot practically be integrated in BEMS. This paper presents and describes validation of a CFD-ROM method for real-time simulations of building thermal performance. The CFD-ROM method involves the automatic extraction and solution of reduced order models (ROMs) from validated CFD simulations. ROMs are shown to be adequately accurate with a total error below 5% and to retain satisfactory representation of the phenomena modelled. Each ROM has a time to solution under 20seconds, which opens the potential of their integration with BEMS, giving real-time physics-based building energy modelling. A parameter study was conducted to investigate the applicability of the extracted ROM to initial boundary conditions different from those from which it was extracted. The results show that the ROMs retained satisfactory total errors when the initial conditions in the room were varied by ±5°C. This allows the production of a finite number of ROMs with the ability to rapidly model many possible scenarios.

Transonic aeroelastic simulation for instability searches and uncertainty analysis

In this paper the use of eigenvalue stability analysis of very large dimension aeroelastic numerical models arising from the exploitation of computational fluid dynamics is reviewed. A formulation based on a block reduction of the system Jacobian proves powerful to allow various numerical algorithms to be exploited, including frequency domain solvers, reconstruction of a term describing the fluid–structure interaction from the sparse data which incurs the main computational cost, and sampling to place the expensive samples where they are most needed. The stability formulation also allows non-deterministic analysis to be carried out very efficiently through the use of an approximate Newton solver. Finally, the system eigenvectors are exploited to produce nonlinear and parameterised reduced order models for computing limit cycle responses. The performance of the methods is illustrated with results from a number of academic and large dimension aircraft test cases.

On how structural model variability influence transonic aeroelasticity stability.

This paper considers the ways in which structural model parameter variability can in?uence aeroelastic stability. Previous work on formulating the stability calculation (with the Euler equations providing the aerodynamic predictions) is exploited to use Monte Carlo, Interval and Perturbation calculations to allow this question to be investigated. Three routes are identi?ed. The ?rst involves variable normal mode frequencies only. The second involves normal mode frequencies and mode shapes. Finally, the third, in addition to normal mode frequencies and mode shapes, also includes their in?uence on the static equilibrium. Previous work has suggested only considering route 1, which allows signi?cant gains in computational e?ciency if reduced order models can be built for the aerodynamics. However, results in the current paper show that neglecting route 2 can give misleading results for the ?utter onset prediction.

How Structural Model Variability Influences Transonic Aeroelastic Stability

This paper considers the ways in which structural model parameter variability can influence aeroelastic stability. Previous work on formulating the stability calculation (with the Euler equations providing the aerodynamic predictions) is exploited to use Monte Carlo, interval, and perturbation calculations to allow this question to be investigated. Three routes are identified. The first involves variable normal-mode frequencies only. The second involves normal-mode frequencies and shapes. Finally, the third, in addition to normal-mode frequencies and shapes, also includes their influence on the static equilibrium. Previous work has suggested only considering the first route, which allows significant gains in computational efficiency if reduced-order models can be built for the aerodynamics. However, results in the current paper show that neglecting the mode-shape variation can give misleading results for the flutter-onset prediction, complicating the development of reduced aerodynamic models for variability analysis.

Development and Validation of a Compact Thermal Model for an Aircraft Compartment

The development of accurate structural/thermal numerical models of complex systems, such as aircraft fuselage barrels, is often limited and determined by the smallest scales that need to be modelled. The development of reduced order models of the smallest scales and consequently their integration with higher level models can be a way to minimise the bottle neck present, while still having efficient, robust and accurate numerical models. In this paper a methodology on how to develop compact thermal fluid models (CTFMs) for compartments where mixed convection regimes are present is demonstrated. Detailed numerical simulations (CFD) have been developed for an aircraft crown compartment and validated against experimental data obtained from a 1:1 scale compartment rig. The crown compartment is defined as the confined area between the upper fuselage and the passenger cabin in a single aisle commercial aircraft. CFD results were utilised to extract average quantities (temperature and heat fluxes) and characteristic parameters (heat transfer coefficients) to generate CTFMs. The CTFMs have then been compared with the results obtained from the detailed models showing average errors for temperature predictions lower than 5%. This error can be deemed acceptable when compared to the nominal experimental error associated with the thermocouple measurements.

The CTFMs methodology developed allows to generate accurate reduced order models where accuracy is restricted to the region of Boundary Conditions applied. This limitation arises from the sensitivity of the internal flow structures to the applied boundary condition set. CTFMs thus generated can be then integrated in complex numerical modelling of whole fuselage sections.

Further steps in the development of an exhaustive methodology would be the implementation of a logic ruled based approach to extract directly from the CFD simulations numbers and positions of the nodes for the CTFM.

Model reduction for nonlinear aerodynamics and aeroelasticity using a Discrete Empirical Interpolation Method

A novel surrogate model is proposed in lieu of Computational Fluid Dynamics (CFD) solvers, for fast nonlinear aerodynamic and aeroelastic modeling. A nonlinear function is identified on selected interpolation points by
a discrete empirical interpolation method (DEIM). The flow field is then reconstructed using a least square approximation of the flow modes extracted
by proper orthogonal decomposition (POD). The aeroelastic reduce order
model (ROM) is completed by introducing a nonlinear mapping function
between displacements and the DEIM points. The proposed model is investigated to predict the aerodynamic forces due to forced motions using
a N ACA 0012 airfoil undergoing a prescribed pitching oscillation. To investigate aeroelastic problems at transonic conditions, a pitch/plunge airfoil
and a cropped delta wing aeroelastic models are built using linear structural models. The presence of shock-waves triggers the appearance of limit
cycle oscillations (LCO), which the model is able to predict. For all cases
tested, the new ROM shows the ability to replicate the nonlinear aerodynamic forces, structural displacements and reconstruct the complete flow
field with sufficient accuracy at a fraction of the cost of full order CFD
model.

Ion-acoustic waves in a plasma consisting of adiabatic warm ions, non-isothermal electrons and a weakly relativistic electron beam: linear and higher-order nonlinear effects

The nonlinear propagation of finite amplitude ion acoustic solitary waves in a plasma consisting of adiabatic warm ions, nonisothermal electrons, and a weakly relativistic electron beam is studied via a two-fluid model. A multiple scales technique is employed to investigate the nonlinear regime. The existence of the electron beam gives rise to four linear ion acoustic modes, which propagate at different phase speeds. The numerical analysis shows that the propagation speed of two of these modes may become complex-valued (i.e., waves cannot occur) under conditions which depend on values of the beam-to-background-electron density ratio , the ion-to-free-electron temperature ratio , and the electron beam velocity v0; the remaining two modes remain real in all cases. The basic set of fluid equations are reduced to a Schamel-type equation and a linear inhomogeneous equation for the first and second-order potential perturbations, respectively. Stationary solutions of the coupled equations are derived using a renormalization method. Higher-order nonlinearity is thus shown to modify the solitary wave amplitude and may also deform its shape, even possibly transforming a simple pulse into a W-type curve for one of the modes. The dependence of the excitation amplitude and of the higher-order nonlinearity potential correction on the parameters , , and v0 is numerically investigated.

Higher order statistics for lognormal small-scale fading in mobile radio channels

This letter reports the statistical characterization and modeling of the indoor radio channel for a mobile wireless personal area network operating at 868 MHz. Line of sight (LOS) and non-LOS conditions were considered for three environments: anechoic chamber, open office area and hallway. Overall, the Nakagami-m cdf best described fading for bodyworn operation in 60% of all measured channels in anechoic chamber and open office area environments. The Nakagami distribution was also found to provide a good description of Rician distributed channels which predominated in the hallway. Multipath played an important role in channel statistics with the mean recorded m value being reduced from 7.8 in the anechoic chamber to 1.3 in both the open office area and hallway.

Even order harmonic Series-L/Parallel-Tuned Class-E frequency multiplier

In this brief, we propose a new Class-E frequency multiplier based on the recently introduced Series-L/Parallel-Tuned Class-E amplifier. The proposed circuit produces even-order output harmonics. Unlike previously reported solutions the proposed circuit can operate under 50% duty ratio which minimizes the conduction losses. The circuit also offers the possibility for increased maximum operating frequency, reduced peak switch voltage, higher load resistance and inherent bond wire absorption; all potentially useful in monolithic microwave integrated circuit implementations. In addition, the circuit topology suggested large transistors with high output capacitances can be deployed. Theoretical design equations are given and the predictions made using these are shown to agree with harmonic balance circuit simulation results.

Differential effects of reduced glycoprotein VI levels on activation of murine platelets by glycoprotein VI ligands

We have investigated the effects of decreased levels of the complex between glycoprotein VI (GPVI) and the Fc receptor gamma-chain (FcRgamma) on responses to collagen and GPVI-specific ligands in murine platelets. We show that levels of GPVI-FcRgamma of the order of 50 % and 20 % of wild-type levels caused 2- and 5-fold shifts to the right respectively in the dose-response curve for aggregation in response to collagen, the snake toxin convulxin and the monoclonal antibody JAQ1. In addition, there is a delay in the onset of aggregation in response to collagen. In contrast, the stimulation of protein tyrosine phosphorylation by collagen (as measured after 150 s) and adhesion to a collagen-coated surface under static conditions were unaffected in platelets with 50 % and 20 % of wild-type levels of GPVI. In contrast, responses to a collagen-related peptide (CRP), made up of repeat glycine-proline-hydroxyproline motifs, were markedly inhibited and abolished in platelets expressing 50 % and 20 % of wild-type levels of GPVI respectively. We suggest that the marked effect of a reduction in GPVI levels on the CRP-induced activation of platelets is due to the multivalent nature of CRP and the fact that GPVI is its sole receptor on platelets. Thus it appears that the interaction of CRP with GPVI is determined by a combination of affinity and avidity. The observation that collagen does not behave like CRP in platelets expressing reduced levels of GPVI, even in the combined presence of blocking antibodies against integrin alpha2beta1 and GPV, suggests that collagen has a greater affinity than CRP for GPVI, and/or that other receptors are involved in its binding to platelets. The clinical significance of these results is discussed.

Entanglement Spectrum, Critical Exponents and Order Parameters in Quantum Spin Chains

We investigate the entanglement spectrum near criticality in finite quantum spin chains. Using finite size scaling we show that when approaching a quantum phase transition, the Schmidt gap, i.e., the difference between the two largest eigenvalues of the reduced density matrix ?1, ?2, signals the critical point and scales with universal critical exponents related to the relevant operators of the corresponding perturbed conformal field theory describing the critical point. Such scaling behavior allows us to identify explicitly the Schmidt gap as a local order parameter.

Newly Acquired Fear of Falling Leads to Altered Eye Movement Patterns and Reduced Stepping Safety: A Case Study

This opportune case study describes visual and stepping behaviours of an 87 year old female (P8), both prior to, and following two falls. Before falling, when asked to walk along a path containing two stepping guides positioned before and after an obstacle, P8 generally visually fixated the first stepping guide until after foot contact inside it. However, after falling P8 consistently looked away from the stepping guide before completing the step into it in order to fixate the upcoming obstacle in her path. The timing of gaze redirection away from the target (in relation to foot contact inside it) correlated with absolute stepping error. No differences in eyesight, cognitive function, or balance were found between pre- and post-fall recordings. However, P8 did report large increases in fall-related anxiety and reduced balance confidence, supporting previously suggested links between anxiety/increased fear or falling and maladaptive visual/stepping behaviours. The results represent a novel insight into how psychological and related behavioural factors can change in older adults following a fall, and provide a possible partial rationalisation for why recent fallers are more likely to fall again in the following 12 months. These findings highlight novel possibilities for falls prevention and rehabilitation.

Loss of Synaptic Connectivity, Particularly in Second Order Neurons Is a Key Feature of Diabetic Retinal Neuropathy in the Ins2Akita Mouse

Retinal neurodegeneration is a key component of diabetic retinopathy (DR), although the detailed neuronal damage remains ill-defined. Recent evidence suggests that in addition to amacrine and ganglion cell, diabetes may also impact on other retinal neurons. In this study, we examined retinal degenerative changes in Ins2^Akita diabetic mice. In scotopic electroretinograms (ERG), b-wave and oscillatory potentials were severely impaired in 9-month old Ins2^Akita mice. Despite no obvious pathology in fundoscopic examination, optical coherence tomography (OCT) revealed a progressive thinning of the retina from 3 months onwards. Cone but not rod photoreceptor loss was observed in 3-month-old diabetic mice. Severe impairment of synaptic connectivity at the outer plexiform layer (OPL) was detected in 9-month old Ins2Akita mice. Specifically, photoreceptor presynaptic ribbons were reduced by 25% and postsynaptic boutons by 70%, although the density of horizontal, rod- and cone-bipolar cells remained similar to non-diabetic controls. Significant reductions in GABAergic and glycinergic amacrine cells and Brn3a+ retinal ganglion cells were also observed in 9-month old Ins2^Akita mice. In conclusion, the Ins2^Akita mouse develops cone photoreceptor degeneration and the impairment of synaptic connectivity at the OPL, predominately resulting from the loss of postsynaptic terminal boutons. Our findings suggest that the Ins2^Akita mouse is a good model to study diabetic retinal neuropathy.