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.

Interface management for automating finite element analysis workflows

This paper outlines the importance of robust interface management for facilitating finite element analysis workflows. Topological equivalences between analysis model representations are identified and maintained in an editable and accessible manner. The model and its interfaces are automatically represented using an analysis-specific cellular decomposition of the design space. Rework of boundary conditions following changes to the design geometry or the analysis idealization can be minimized by tracking interface dependencies. Utilizing this information with the Simulation Intent specified by an analyst, automated decisions can be made to process the interface information required to rebuild analysis models. Through this work automated boundary condition application is realized within multi-component, multi-resolution and multi-fidelity analysis workflows.

Predicting low velocity impact damage and Compression-After-Impact (CAI) behaviour of composite laminates

Low-velocity impact damage can drastically reduce the residual strength of a composite structure even when the damage is barely visible. The ability to computationally predict the extent of damage and compression-after-impact (CAI) strength of a composite structure can potentially lead to the exploration of a larger design space without incurring significant time and cost penalties. A high-fidelity three-dimensional composite damage model, to predict both low-velocity impact damage and CAI strength of composite laminates, has been developed and implemented as a user material subroutine in the commercial finite element package, ABAQUS/Explicit. The intralaminar damage model component accounts for physically-based tensile and compressive failure mechanisms, of the fibres and matrix, when subjected to a three-dimensional stress state. Cohesive behaviour was employed to model the interlaminar failure between plies with a bi-linear traction–separation law for capturing damage onset and subsequent damage evolution. The virtual tests, set up in ABAQUS/Explicit, were executed in three steps, one to capture the impact damage, the second to stabilize the specimen by imposing new boundary conditions required for compression testing, and the third to predict the CAI strength. The observed intralaminar damage features, delamination damage area as well as residual strength are discussed. It is shown that the predicted results for impact damage and CAI strength correlated well with experimental testing without the need of model calibration which is often required with other damage models.

Impact of finite element idealisation on the prediction of welded fuselage stiffened panel buckling

Lap joints are widely used in the manufacture of stiffened panels and influence local panel sub-component stability, defining buckling unit dimensions and boundary conditions. Using the Finite Element method it is possible to model joints in great detail and predict panel buckling behaviour with accuracy. However, when modelling large panel structures such detailed analysis becomes computationally expensive. Moreover, the impact of local behaviour on global panel performance may reduce as the scale of the modelled structure increases. Thus this study presents coupled computational and experimental analysis, aimed at developing relationships between modelling fidelity and the size of the modelled structure, when the global static load to cause initial buckling is the required analysis output. Small, medium and large specimens representing welded lap-joined fuselage panel structure are examined. Two element types, shell and solid-shell, are employed to model each specimen, highlighting the impact of idealisation on the prediction of welded stiffened panel initial skin buckling.

Finite Element Analysis of the Enhanced Strength of Laterally Restrained RC Slabs

This paper presents the numerical simulation of the ultimate behaviour of 85 one-way and two-way spanning laterally restrained concrete slabs of variable thickness, span, reinforcement ratio, strength and boundary conditions reported in literature by different authors. The developed numerical model was described and all the assumptions were illustrated. ABAQUS, a Finite Element Analysis suite of software, was employed. Non-linear implicit static general analysis method offered by ABAQUS was used. Other analysis methods were also discussed in general in terms of application such as Explicit Dynamic Analysis and Riks method. The aim is to demonstrate the ability and efficacy of FEA to simulate the ultimate load behaviour of slabs considering different material properties and boundary conditions. The authors intended to present a numerical model that provides consistent predictions of the ultimate behaviour of laterally restrained slabs that could be used as an alternative for expensive real life testing as well as for the design and assessment of new and existing structures respectively. The enhanced strength of laterally-restrained slabs compared with conventional design methods predictions is believed to be due to compressive membrane action (CMA). CMA is an inherent phenomenon of laterally restrained concrete beams/slabs. The numerical predictions obtained from the developed model were in good correlation with the experimental results and with those obtained from the CMA method developed at the Queen’s University Belfast, UK.

The Reliability of the Arc-Length Method in the Analysis of Mode-Jumping Problems

The Arc-Length Method is a solution procedure that enables a generic non-linear problem to pass limit points. Some examples are provided of mode-jumping problems solutions using a commercial nite element package, and other investigations are carried out on a simple structure of which the numerical solution can be compared with an analytical one. It is shown that Arc-Length Method is not reliable when bifurcations are present in the primary equilibrium path; also the presence of very sharp snap-backs or special boundary conditions may cause convergence diÆculty at limit points. An improvement to the predictor used in the incremental procedure is suggested, together with a reliable criteria for selecting either solution of the quadratic arc-length constraint. The gap that is sometimes observed between the experimantal load level of mode-jumping and its arc-length prediction is explained through an example.

Finite-element investigation of cold-formed steel portal frames in fire

This paper presents the results of a full-scale site fire test performed on a cold-formed steel portal frame building with semi-rigid joints. The purpose of the study is to establish a performance-based approach for the design of such structures in fire boundary conditions. In the full-scale site fire test, the building collapsed asymmetrically at a temperature of 714°C. A non-linear elasto-plastic finite-element shell model is described and is validated against the results of the full-scale test. A parametric study is presented that highlights the importance of in-plane restraint from the side rails in preventing an outwards sway failure for both a single portal and full building geometry model. The study also demonstrates that the semi-rigidity of the joints should be taken into account in the design. The single portal and full building geometry models display a close match to site test results with failure at 682°C and 704°C, respectively. A design case is described in accordance with Steel Construction Institute design recommendations. The validated single portal model is tested with pinned bases, columns protected, realistic loading and rafters subject to symmetric uniform heating in accordance with the ISO 834 standard fire curve; failure occurs at 703°C.

Magnetospheric Structure and Atmospheric Joule Heating of Habitable Planets Orbiting M-dwarf Stars

We study the magnetospheric structure and the ionospheric Joule Heating of planets orbiting M-dwarf stars in the habitable zone using a set of magnetohydrodynamic models. The stellar wind solution is used to drive a model for the planetary magnetosphere, which is coupled with a model for the planetary ionosphere. Our simulations reveal that the space environment around close-in habitable planets is extreme, and the stellar wind plasma conditions change from sub- to super-Alfvénic along the planetary orbit. As a result, the magnetospheric structure changes dramatically with a bow shock forming in the super-Alfvénic sectors, while no bow shock forms in the sub-Alfvénic sectors. The planets reside most of the time in the sub-Alfvénic sectors with poor atmospheric protection. A significant amount of Joule Heating is provided at the top of the atmosphere as a result of the intense stellar wind. For the steady-state solution, the heating is about 0.1%-3% of the total incoming stellar irradiation, and it is enhanced by 50% for the time-dependent case. The significant Joule Heating obtained here should be considered in models for the atmospheres of habitable planets in terms of the thickness of the atmosphere, the top-side temperature and density, the boundary conditions for the atmospheric pressure, and particle radiation and transport. Here we assume constant ionospheric Pedersen conductance similar to that of the Earth. The conductance could be greater due to the intense EUV radiation leading to smaller heating rates. We plan to quantify the ionospheric conductance in future study.

Electron-impact ionization of diatomic molecules using a configuration-average distorted-wave method

Electron-impact ionization cross sections for diatomic molecules are calculated in a configuration-average distorted-wave method. Core bound orbitals for the molecular ion are calculated using a single-configuration self-consistent-field method based on a linear combination of Slater-type orbitals. The core bound orbitals are then transformed onto a two-dimensional (r,θ) numerical lattice from which a Hartree potential with local exchange is constructed. The single-particle Schrödinger equation is then solved for the valence bound orbital and continuum distorted-wave orbitals with S-matrix boundary conditions. Total cross section results for H2 and N2 are compared with those from semiempirical calculations and experimental measurements.

Existência de minimizantes relaxados e fenómeno de Lavrentiev

Neste trabalho prova-se a existência de minimizantes relaxados em problemas de controlo óptimo não convexos usando técnicas de compactificação. Faz-se a extensão do exemplo de Manià a dimensão dois, obtendo-se uma classe de problemas variacionais em 2D que apresentam Fenómeno de Lavrentiev. Prova-se que o fenómeno persiste a certas perturbações, obtendo- -se assim uma classe de funcionais cujos Lagrangianos são coercivos e convexos em relação ao gradiente. Adicionalmente, apresentam-se exemplos de problemas do cálculo das variações com diferentes condições de fronteira, e em diferentes tipos de domínios (incluindo domínios com fronteira fractal), que exibem Fenómeno de Lavrentiev.

Cálculo das variações fraccional generalizado

Nesta tese de doutoramento apresentamos um cálculo das variações fraccional generalizado. Consideramos problemas variacionais com derivadas e integrais fraccionais generalizados e estudamo-los usando métodos directos e indirectos. Em particular, obtemos condições necessárias de optimalidade de Euler-Lagrange para o problema fundamental e isoperimétrico, condições de transversalidade e teoremas de Noether. Demonstramos a existência de soluções, num espaço de funções apropriado, sob condições do tipo de Tonelli. Terminamos mostrando a existência de valores próprios, e correspondentes funções próprias ortogonais, para problemas de Sturm- Liouville.

Wind energy resource modelling in Portugal and its future large-scale alteration due to anthropogenic induced climate changes

The high dependence of Portugal from foreign energy sources (mainly fossil fuels), together with the international commitments assumed by Portugal and the national strategy in terms of energy policy, as well as resources sustainability and climate change issues, inevitably force Portugal to invest in its energetic self-sufficiency. The 20/20/20 Strategy defined by the European Union defines that in 2020 60% of the total electricity consumption must come from renewable energy sources. Wind energy is currently a major source of electricity generation in Portugal, producing about 23% of the national total electricity consumption in 2013. The National Energy Strategy 2020 (ENE2020), which aims to ensure the national compliance of the European Strategy 20/20/20, states that about half of this 60% target will be provided by wind energy. This work aims to implement and optimise a numerical weather prediction model in the simulation and modelling of the wind energy resource in Portugal, both in offshore and onshore areas. The numerical model optimisation consisted in the determination of which initial and boundary conditions and planetary boundary layer physical parameterizations options provide wind power flux (or energy density), wind speed and direction simulations closest to in situ measured wind data. Specifically for offshore areas, it is also intended to evaluate if the numerical model, once optimised, is able to produce power flux, wind speed and direction simulations more consistent with in situ measured data than wind measurements collected by satellites. This work also aims to study and analyse possible impacts that anthropogenic climate changes may have on the future wind energetic resource in Europe. The results show that the ECMWF reanalysis ERA-Interim are those that, among all the forcing databases currently available to drive numerical weather prediction models, allow wind power flux, wind speed and direction simulations more consistent with in situ wind measurements. It was also found that the Pleim-Xiu and ACM2 planetary boundary layer parameterizations are the ones that showed the best performance in terms of wind power flux, wind speed and direction simulations. This model optimisation allowed a significant reduction of the wind power flux, wind speed and direction simulations errors and, specifically for offshore areas, wind power flux, wind speed and direction simulations more consistent with in situ wind measurements than data obtained from satellites, which is a very valuable and interesting achievement. This work also revealed that future anthropogenic climate changes can negatively impact future European wind energy resource, due to tendencies towards a reduction in future wind speeds especially by the end of the current century and under stronger radiative forcing conditions.

Quasi steady state and dynamic hygrothermal performance of fibrous Hemp and Stone Wool insulations: Two innovative laboratory based investigations

Moisture and heat management properties of Hemp and Stone Wool insulations were studied by mounting them between a hot and a cold climate chamber. Both insulations were exposed to identical hygrothermal boundary conditions. Quasi steady state and dynamic tests were carried out at a range of relative humidity exposures. The likelihood of interstitial condensation was assessed and equivalent thermal conductivity values of the insulations were determined. The adsorption-desorption isotherms of the insulations were also determined in a dynamic vapour sorption (DVS) instrument. It was observed that the likelihood of condensation was higher in Stone Wool insulation than in Hemp insulation. Hemp insulation performed better in managing moisture due to its high hygric inertia and water absorption capacity. It was observed that the equivalent thermal conductivity of Stone Wool insulation was dependent on enthalpy flow and phase change of moisture. The equivalent thermal conductivity of Hemp insulation was close to its declared thermal conductivity in dynamic conditions when high relative humidity exposures were transient. In quasi steady state boundary conditions, when the insulation was allowed to reach the equilibrium moisture content at ranges of relative humidity, there was a moisture dependent increase of thermal conductivity in Hemp insulation.

Hydrodynamic model intercomparison for the patos lagoon (Brazil)

An intercomparison study is carried out between two models with different formulations and spatial discretizations in order to overcome the limitations posed by the standard calibration and validation procedures and improve confidence in the hydrodynamic results for the Patos Lagoon. Numerical simulations were carried out applying the TELEMAC and MOHID models, based on the same boundary conditions and identical calibration coefficients so differences in calculated flow conditions result from the formulations and parameterizations of each model. Results from both models are compared with measurements from three stations inside the lagoon. Preliminary results indicate that both models compare well with the measurements and with each other. These results increase the confidence on hydrodynamic results for the Patos Lagoon and provide the first step towards water quality studies for the area.

Towards a hydrodynamic operational model of the Algarve coast and St. Vincent Cape

This study attempts to implement a hydrodynamic operational model which can ultimately be used for projecting oil spill dispersal patterns and also sewage, pollution and can also be used in wave forecasting. A two layer nested model was created using MOHID Water, which is powerful ocean modelling software. The first layer (father) is used to impose the boundary conditions for the second layer (son). This was repeated for two different wind dominant regimes, Easterly and Westerly winds respectively. A qualitative comparison was done between measured tidal data and the tidal output. Sea surface temperature was also qualitatively compared with the model’s results. The results from both simulations were analysed and compared to historical literature. The comparison was done at the surface layer, 100 metre depth and at 800m depth. In the surface layer the first simulation generated an upwelling event near Cape St. Vincent and within the Algarve. The second simulation generated a non-upwelling event within which the surface was flow reversed and the warm water mass was along the Algarve coastline and evening turning clockwise around Cape St. Vincent. At the 100 metre depth for both simulations, velocity vortexes were observed near Cape St. Vincent travelling northerly and southerly at various instances. At 800metre depth a strong oceanic flow was observed moving north westerly along the continental shelf.