GPU Acceleration of Partial Differential Equation Solvers

Differential equations are often directly solvable by analytical means only in their one dimensional version. Partial differential equations are generally not solvable by analytical means in two and three dimensions, with the exception of few special cases. In all other cases, numerical approximation methods need to be utilized. One of the most popular methods is the finite element method. The main areas of focus, here, are the Poisson heat equation and the plate bending equation. The purpose of this paper is to provide a quick walkthrough of the various approaches that the authors followed in pursuit of creating optimal solvers, accelerated with the use of graphical processing units, and comparing them in terms of accuracy and time efficiency with existing or self-made non-accelerated solvers.

Efficiency analysis of flat plate collectors for building facade integration

A theoretical approach has been used to evaluate the performance of facade integrated solar collectors based on the physical collector parameters such as absorber plate absorptance, transmittance of the glazed cover plate and insulation thickness. A 1D steady state model, based on the Hottel Whillier Bliss equation, was employed to determine the effect of changing parameters to meet façade integration criteria.

Difference equation approach to two-thermocouple sensor characterisation in constant velocity flow environments

Thermocouples are one of the most popular devices for temperature measurement due to their robustness, ease of manufacture and installation, and low cost. However, when used in certain harsh environments, for example, in combustion systems and engine exhausts, large wire diameters are required, and consequently the measurement bandwidth is reduced. This article discusses a software compensation technique to address the loss of high frequency fluctuations based on measurements from two thermocouples. In particular, a difference equation sDEd approach is proposed and compared with existing methods both in simulation and on experimental test rig data with constant flow velocity. It is found that the DE algorithm, combined with the use of generalized total least squares for parameter identification, provides better performance in terms of time constant estimation without any a priori assumption on the time constant ratios of the thermocouples.

Exploiting a priori time constant ratio information in difference equation two-thermocouple sensor characterisation

The characterization of thermocouple sensors for temperature measurement in varying-flow environments is a challenging problem. Recently, the authors introduced novel difference-equation-based algorithms that allow in situ characterization of temperature measurement probes consisting of two-thermocouple sensors with differing time constants. In particular, a linear least squares (LS) lambda formulation of the characterization problem, which yields unbiased estimates when identified using generalized total LS, was introduced. These algorithms assume that time constants do not change during operation and are, therefore, appropriate for temperature measurement in homogenous constant-velocity liquid or gas flows. This paper develops an alternative ß-formulation of the characterization problem that has the major advantage of allowing exploitation of a priori knowledge of the ratio of the sensor time constants, thereby facilitating the implementation of computationally efficient algorithms that are less sensitive to measurement noise. A number of variants of the ß-formulation are developed, and appropriate unbiased estimators are identified. Monte Carlo simulation results are used to support the analysis.

Improved Equivalent Frame Analysis Method for Flat Plate Structures in Vicinity of Edge Columns

This paper proposes a modification to the ACI 318-02 equivalent frame method of analysis of reinforced concrete flat plate exterior panels. Two existing code methods were examined: ACI 318 and BS 8110. The derivation of the torsional stiffness of the edge strip as proposed by ACI 318 is examined and a more accurate estimate of this value is proposed, based on both theoretical analysis and experimental results. A series of 1/3-scale models of flat plate exterior panels have been tested. Unique experimental results were obtained by measuring strains in reinforcing bars at approximately 200 selected locations in the plate panel throughout the entire loading history. The measured strains were used to calculate curvature and, hence, bending moments; these were used along with moments in the columns to assess the accuracy of the equivalent frame methods. The proposed method leads to a more accurate prediction of the moments in the plate at the column front face, at the panel midspan, and in the edge column. Registered Subscribers: View the full article. This document is available as a free download to qualified members. An electronic (PDF) version is available for purchase and download. Click on the Order Now button to continue with the download.

Solving the advection-dispersion equation using the discrete puff particle method

A flexible, mass-conservative numerical technique for solving the advection-dispersion equation for miscible contaminant transport is presented. The method combines features of puff transport models from air pollution studies with features from the random walk particle method used in water resources studies, providing a deterministic time-marching algorithm which is independent of the grid Peclet number and scales from one to higher dimensions simply. The concentration field is discretised into a number of particles, each of which is treated as a point release which advects and disperses over the time interval. The dispersed puff is itself discretised into a spatial distribution of particles whose masses can be pre-calculated. Concentration within the simulation domain is then calculated from the mass distribution as an average over some small volume. Comparison with analytical solutions for a one-dimensional fixed-duration concentration pulse and for two-dimensional transport in an axisymmetric flow field indicate that the algorithm performs well. For a given level of accuracy the new method has lower computation times than the random walk particle method.

A Van der Pol-Mathieu equation for the dynamics of dust grain charge in dusty plasmas

The chaotic profile of dust grain dynamics associated with dust-acoustic oscillations in a dusty plasma is considered. The collective behaviour of the dust plasma component is described via a multi-fluid model, comprising Boltzmann distributed electrons and ions, as well as an equation of continuity possessing a source term for the dust grains, the dust momentum and Poisson's equations. A Van der Pol–Mathieu-type nonlinear ordinary differential equation for the dust grain density dynamics is derived. The dynamical system is cast into an autonomous form by employing an averaging method. Critical stability boundaries for a particular trivial solution of the governing equation with varying parameters are specified. The equation is analysed to determine the resonance region, and finally numerically solved by using a fourth-order Runge–Kutta method. The presence of chaotic limit cycles is pointed out.

Difference equation sensor characterization algorithms for two-thermocouple probes

The characterization of thermocouple sensors for temperature measurement in variable flow environments is a challenging problem. In this paper, novel difference equation-based algorithms are presented that allow in situ characterization of temperature measurement probes consisting of two-thermocouple sensors with differing time constants. Linear and non-linear least squares formulations of the characterization problem are introduced and compared in terms of their computational complexity, robustness to noise and statistical properties. With the aid of this analysis, least squares optimization procedures that yield unbiased estimates are identified. The main contribution of the paper is the development of a linear two-parameter generalized total least squares formulation of the sensor characterization problem. Monte-Carlo simulation results are used to support the analysis.

Stiffened panel stability behaviour and performance gains with plate prismatic sub-stiffening

To increase the structural efficiency of integrally machined aluminium alloy stiffened panels, it is plausible to introduce plate sub-stiffening to increase the local stability and thus panel static strength performance. Reported herein is the experimental validation of prismatic sub-stiffening, and the computational verification of such concepts within larger recurring structure. The experimental work demonstrates the potential to 'control' plate buckling modes. For the tested sub-stiffening design, an initial plate buckling performance gain of +89% over an equivalent mass design was measured. The numerical simulations, modelling the tested sub-stiffening design, demonstrate equivalent behaviour and performance gains (+66%) within larger structures consisting of recurring panels. (C) 2009 Elsevier Ltd. All rights reserved.