Thermoeconomic analysis method for optimization of combined heat and power systems - part II

In this paper, a thermoeconomic functional analysis method based on the Second Law of Thermodynamics and applied to analyze four cogeneration systems is presented. The objective of the developed technique is to minimize the operating costs of the cogeneration plant, namely exergetic production cost (EPC), assuming fixed rates of electricity production and process steam in exergy base. In this study a comparison is made between the same four configurations of part I. The cogeneration system consisting of a gas turbine with a heat recovery steam generator, without supplementary firing, has the lowest EPC. (C) 2004 Published by Elsevier Ltd.

Thermoeconomic analysis method for optimization of combined heat and power systems. Part I

This article presents a thermoeconomic analysis of cogeneration plants, applied as a rational technique to produce electric power and saturated steam. The aim of this new methodology is the minimum Exergetic Production Cost (EPC), based on the Second Law of Thermodynamics. The variables selected for the optimization are the pressure and the temperature of the steam leaving the boiler in the case of using steam turbine, and the pressure ratio, turbine exhaust temperature and mass flow in the case of using gas turbines. The equations for calculating the capital costs of the components and products are formulated as a function of these decision variables. An application of the method using real data of a multinational chemical industry located in São Paulo state is presented. The conditions which establish the minimum cost are presented as final output. (C) 2003 Elsevier Ltd. All rights reserved.

Generalized Coupled Algebraic Riccati Equations for Discrete-time Markov Jump with Multiplicative Noise Systems

In this paper we consider the existence of the maximal and mean square stabilizing solutions for a set of generalized coupled algebraic Riccati equations (GCARE for short) associated to the infinite-horizon stochastic optimal control problem of discrete-time Markov jump with multiplicative noise linear systems. The weighting matrices of the state and control for the quadratic part are allowed to be indefinite. We present a sufficient condition, based only on some positive semi-definite and kernel restrictions on some matrices, under which there exists the maximal solution and a necessary and sufficient condition under which there exists the mean square stabilizing solution fir the GCARE. We also present a solution for the discounted and long run average cost problems when the performance criterion is assumed be composed by a linear combination of an indefinite quadratic part and a linear part in the state and control variables. The paper is concluded with a numerical example for pension fund with regime switching.

Assessment of skinfold thickness equations in estimating body composition in children with inflammatory bowel

Introduction: Growth is a central process in paediatrics. Weight and height evaluation are therefore routine exams for every child but in some situation, particularly inflammatory bowel disease (IBD), a wider evaluation of nutritional status needs to be performed. The assessment of body composition is essential in order to maintain acceptable growth using the following techniques: Dual-energy X-ray absorptiometry (DEXA), bio-impedance-analysis (BIA) and anthropometric measurements (skinfold thickness skin), the latter being most easily available and most cost effective. Objectives: To assess the accuracy of skinfold equations in estimating percentage body fat (%BF) in children with inflammatory bowel disease (IBD), compared with assessment of body fat dual energy X-ray absorptiometry (DEXA). Methods: Twenty-one patients (11 females, 10 males; mean age: 14.3 years, range 12 - 16 years) with IBD (Crohn's disease n = 15, ulcerative colitis n = 6)). Estimated%BF was computed using 6 established equations based on the triceps, biceps, subscapular and suprailiac skinfolds (Deurenberg, Weststrate, Slaughter, Durnin & Rahaman, Johnston, Brook) and compared to DEXA. Concordance analysis was performed using Lin's concordance correlation and the Bland-Altman limits of agreement method. Results: Durnin & Rahaman's equation shows a higher Lin's concordance coefficient with a small difference amongst raw values for skinfolds and DEXA compared to the other equations. Correlation coefficient between mean and difference is close to zero with a non-significant Bradley-Blackwood test. Conclusion: Body composition in paediatric IBD patients using the Durnin & Rahaman skinfold-equation adequately reflects values obtained by DEXA.

Voronoi, Delaunay, and Block-Structured Mesh Refinement for Solution of the Shallow-Water Equations on the Sphere

Alternative meshes of the sphere and adaptive mesh refinement could be immensely beneficial for weather and climate forecasts, but it is not clear how mesh refinement should be achieved. A finite-volume model that solves the shallow-water equations on any mesh of the surface of the sphere is presented. The accuracy and cost effectiveness of four quasi-uniform meshes of the sphere are compared: a cubed sphere, reduced latitude–longitude, hexagonal–icosahedral, and triangular–icosahedral. On some standard shallow-water tests, the hexagonal–icosahedral mesh performs best and the reduced latitude–longitude mesh performs well only when the flow is aligned with the mesh. The inclusion of a refined mesh over a disc-shaped region is achieved using either gradual Delaunay, gradual Voronoi, or abrupt 2:1 block-structured refinement. These refined regions can actually degrade global accuracy, presumably because of changes in wave dispersion where the mesh is highly nonuniform. However, using gradual refinement to resolve a mountain in an otherwise coarse mesh can improve accuracy for the same cost. The model prognostic variables are height and momentum collocated at cell centers, and (to remove grid-scale oscillations of the A grid) the mass flux between cells is advanced from the old momentum using the momentum equation. Quadratic and upwind biased cubic differencing methods are used as explicit corrections to a fast implicit solution that uses linear differencing.

Comparison of the computational cost of a Monte Carlo and deterministic algorithm for computing bilinear forms of matrix powers

In this paper we consider bilinear forms of matrix polynomials and show that these polynomials can be used to construct solutions for the problems of solving systems of linear algebraic equations, matrix inversion and finding extremal eigenvalues. An almost Optimal Monte Carlo (MAO) algorithm for computing bilinear forms of matrix polynomials is presented. Results for the computational costs of a balanced algorithm for computing the bilinear form of a matrix power is presented, i.e., an algorithm for which probability and systematic errors are of the same order, and this is compared with the computational cost for a corresponding deterministic method.

Solving complex optimal control problems at no cost with PSOPT

This paper introduces PSOPT, an open source optimal control solver written in C++. PSOPT uses pseudospectral and local discretizations, sparse nonlinear programming, automatic differentiation, and it incorporates automatic scaling and mesh refinement facilities. The software is able to solve complex optimal control problems including multiple phases, delayed differential equations, nonlinear path constraints, interior point constraints, integral constraints, and free initial and/or final times. The software does not require any non-free platform to run, not even the operating system, as it is able to run under Linux. Additionally, the software generates plots as well as LATEX code so that its results can easily be included in publications. An illustrative example is provided.

A fast two-grid and finite section method for a class of integral equations on the real line with application to an acoustic scattering problem in the half-plane

We consider the numerical treatment of second kind integral equations on the real line of the form ∅(s) = ∫_(-∞)^(+∞)▒〖κ(s-t)z(t)ϕ(t)dt,s=R〗 (abbreviated ϕ= ψ+K_z ϕ) in which K ϵ L_1 (R), z ϵ L_∞ (R) and ψ ϵ BC(R), the space of bounded continuous functions on R, are assumed known and ϕ ϵ BC(R) is to be determined. We first derive sharp error estimates for the finite section approximation (reducing the range of integration to [-A, A]) via bounds on (1-K_z )^(-1)as an operator on spaces of weighted continuous functions. Numerical solution by a simple discrete collocation method on a uniform grid on R is then analysed: in the case when z is compactly supported this leads to a coefficient matrix which allows a rapid matrix-vector multiply via the FFT. To utilise this possibility we propose a modified two-grid iteration, a feature of which is that the coarse grid matrix is approximated by a banded matrix, and analyse convergence and computational cost. In cases where z is not compactly supported a combined finite section and two-grid algorithm can be applied and we extend the analysis to this case. As an application we consider acoustic scattering in the half-plane with a Robin or impedance boundary condition which we formulate as a boundary integral equation of the class studied. Our final result is that if z (related to the boundary impedance in the application) takes values in an appropriate compact subset Q of the complex plane, then the difference between ϕ(s)and its finite section approximation computed numerically using the iterative scheme proposed is ≤C_1 [kh log⁡〖(1⁄kh)+(1-Θ)^((-1)⁄2) (kA)^((-1)⁄2) 〗 ] in the interval [-ΘA,ΘA](Θ<1) for kh sufficiently small, where k is the wavenumber and h the grid spacing. Moreover this numerical approximation can be computed in ≤C_2 N log⁡N operations, where N = 2A/h is the number of degrees of freedom. The values of the constants C1 and C2 depend only on the set Q and not on the wavenumber k or the support of z.

A low-cost 3D reconstruction system using a single-shot projection of a pattern matrix

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Development of a methodology for cost determination of wastewater treatment based on functional diagram

This work describes a methodology developed for determination of costs associated to products generated in a small wastewater treatment station for sanitary wastewater from a university campus. This methodology begins with plant component units identification, relating their fluid and thermodynamics features for each point marked in its process diagram. Following, its functional diagram is developed and its formulation is elaborated, in exergetic base, describing all equations for these points, which are the constraints for exergetic production cost problem and are used in equations to determine the costs associated to products generated in SWTS. This methodology was applied to a hypothetical system based on SWTS former parts and presented consistent results when compared to expected values based on previous exergetic expertise. (C) 2008 Elsevier Ltd. All rights reserved.

Development of equations to estimate microbial contamination in ruminal incubation residues of forage produced under tropical conditions using 15N as a label1

The objective of this study was to use 15N to label microbial cells to allow development of equations for estimating the microbial contamination in ruminal in situ incubation residues of forage produced under tropical conditions. A total of 24 tropical forages were ruminal incubated in 3 steers at 3 separate times. To determine microbial contamination of the incubated residues, ruminal bacteria were labeled with 15N by continuous intraruminal infusion 60 h before the first incubation and continued until the last day of incubation. Ruminal digesta was collected for the isolation of bacteria before the first infusion of 15N on adaptation period and after the infusion of 15N on collection period. To determine the microbial contamination of CP fractions, restricted models were compared with the full model using the model identity test. A value of the corrected fraction A was estimated from the corresponding noncorrected fraction by this equation: Corrected A fraction (ACPC) = 1.99286 + 0.98256 × A fraction without correction (ACPWC). The corrected fraction B was estimated from the corresponding noncorrected fraction and from CP, NDF, neutral detergent insoluble protein (NDIP), and indigestible NDF (iNDF) using the equation corrected B fraction (BCPC) = -17.2181 - 0.0344 × fraction B without correction (BCPWC) + 0.65433 × CP + 1.03787 × NDF + 2.66010 × NDIP - 0.85979 × iNDF. The corrected degradation rate of B fraction (kd)was estimated using the equation corrected degradation rate of B fraction (kdCPC) = 0.04667 + 0.35139 × degradation rate of B fraction without correction (kdCPWC) + 0.0020 × CP - 0.00055839 × NDF - 0.00336 × NDIP + 0.00075089 × iNDF. This equation was obtained to estimate the contamination using CP of the feeds: %C = 79.21 × (1 - e-0.0555t) × e-0.0874CP. It was concluded that A and B fractions and kd of CP could be highly biased by microbial CP contamination, and therefore these corrected values could be obtained mathematically, replacing the use of microbial markers. The percentage of contamination and the corrected apparent degradability of CP could be obtained from values of CP and time of incubation for each feed, which could reduce cost and labor involved when using 15N. © 2013 American Society of Animal Science. All rights reserved.

Matrix-vector multiplication and triangular linear solver using GPGPU for symmetric positive definite matrices derived from elliptic equations

The modern GPUs are well suited for intensive computational tasks and massive parallel computation. Sparse matrix multiplication and linear triangular solver are the most important and heavily used kernels in scientific computation, and several challenges in developing a high performance kernel with the two modules is investigated. The main interest it to solve linear systems derived from the elliptic equations with triangular elements. The resulting linear system has a symmetric positive definite matrix. The sparse matrix is stored in the compressed sparse row (CSR) format. It is proposed a CUDA algorithm to execute the matrix vector multiplication using directly the CSR format. A dependence tree algorithm is used to determine which variables the linear triangular solver can determine in parallel. To increase the number of the parallel threads, a coloring graph algorithm is implemented to reorder the mesh numbering in a pre-processing phase. The proposed method is compared with parallel and serial available libraries. The results show that the proposed method improves the computation cost of the matrix vector multiplication. The pre-processing associated with the triangular solver needs to be executed just once in the proposed method. The conjugate gradient method was implemented and showed similar convergence rate for all the compared methods. The proposed method showed significant smaller execution time.

Analysis of optimal control problems for the incompressible MHD equations and implementation in a finite element multiphysics code

This thesis deals with the study of optimal control problems for the incompressible Magnetohydrodynamics (MHD) equations. Particular attention to these problems arises from several applications in science and engineering, such as fission nuclear reactors with liquid metal coolant and aluminum casting in metallurgy. In such applications it is of great interest to achieve the control on the fluid state variables through the action of the magnetic Lorentz force. In this thesis we investigate a class of boundary optimal control problems, in which the flow is controlled through the boundary conditions of the magnetic field. Due to their complexity, these problems present various challenges in the definition of an adequate solution approach, both from a theoretical and from a computational point of view. In this thesis we propose a new boundary control approach, based on lifting functions of the boundary conditions, which yields both theoretical and numerical advantages. With the introduction of lifting functions, boundary control problems can be formulated as extended distributed problems. We consider a systematic mathematical formulation of these problems in terms of the minimization of a cost functional constrained by the MHD equations. The existence of a solution to the flow equations and to the optimal control problem are shown. The Lagrange multiplier technique is used to derive an optimality system from which candidate solutions for the control problem can be obtained. In order to achieve the numerical solution of this system, a finite element approximation is considered for the discretization together with an appropriate gradient-type algorithm. A finite element object-oriented library has been developed to obtain a parallel and multigrid computational implementation of the optimality system based on a multiphysics approach. Numerical results of two- and three-dimensional computations show that a possible minimum for the control problem can be computed in a robust and accurate manner.

The determinants of improvements in health outcomes and of cost reduction in hospital inpatient care

This study aims to address two research questions. First, ‘Can we identify factors that are determinants both of improved health outcomes and of reduced costs for hospitalized patients with one of six common diagnoses?’ Second, ‘Can we identify other factors that are determinants of improved health outcomes for such hospitalized patients but which are not associated with costs?’ The Healthcare Cost and Utilization Project (HCUP) Nationwide Inpatient Sample (NIS) database from 2003 to 2006 was employed in this study. The total study sample consisted of hospitals which had at least 30 patients each year for the given diagnosis: 954 hospitals for acute myocardial infarction (AMI), 1552 hospitals for congestive heart failure (CHF), 1120 hospitals for stroke (STR), 1283 hospitals for gastrointestinal hemorrhage (GIH), 979 hospitals for hip fracture (HIP), and 1716 hospitals for pneumonia (PNE). This study used simultaneous equations models to investigate the determinants of improvement in health outcomes and of cost reduction in hospital inpatient care for these six common diagnoses. In addition, the study used instrumental variables and two-stage least squares random effect model for unbalanced panel data estimation. The study concluded that a few factors were determinants of high quality and low cost. Specifically, high specialty was the determinant of high quality and low costs for CHF patients; small hospital size was the determinant of high quality and low costs for AMI patients. Furthermore, CHF patients who were treated in Midwest, South, and West region hospitals had better health outcomes and lower hospital costs than patients who were treated in Northeast region hospitals. Gastrointestinal hemorrhage and pneumonia patients who were treated in South region hospitals also had better health outcomes and lower hospital costs than patients who were treated in Northeast region hospitals. This study found that six non-cost factors were related to health outcomes for a few diagnoses: hospital volume, percentage emergency room admissions for a given diagnosis, hospital competition, specialty, bed size, and hospital region.^