Verification of a mixed high-order accurate DNS code for laminar turbulent transition by the method of manufactured solutions

This paper presents results on a verification test of a Direct Numerical Simulation code of mixed high-order of accuracy using the method of manufactured solutions (MMS). This test is based on the formulation of an analytical solution for the Navier-Stokes equations modified by the addition of a source term. The present numerical code was aimed at simulating the temporal evolution of instability waves in a plane Poiseuille flow. The governing equations were solved in a vorticity-velocity formulation for a two-dimensional incompressible flow. The code employed two different numerical schemes. One used mixed high-order compact and non-compact finite-differences from fourth-order to sixth-order of accuracy. The other scheme used spectral methods instead of finite-difference methods for the streamwise direction, which was periodic. In the present test, particular attention was paid to the boundary conditions of the physical problem of interest. Indeed, the verification procedure using MMS can be more demanding than the often used comparison with Linear Stability Theory. That is particularly because in the latter test no attention is paid to the nonlinear terms. For the present verification test, it was possible to manufacture an analytical solution that reproduced some aspects of an instability wave in a nonlinear stage. Although the results of the verification by MMS for this mixed-order numerical scheme had to be interpreted with care, the test was very useful as it gave confidence that the code was free of programming errors. Copyright (C) 2009 John Wiley & Sons, Ltd.

Recursive linear estimation for general discrete-time descriptor systems

This paper considers the optimal linear estimates recursion problem for discrete-time linear systems in its more general formulation. The system is allowed to be in descriptor form, rectangular, time-variant, and with the dynamical and measurement noises correlated. We propose a new expression for the filter recursive equations which presents an interesting simple and symmetric structure. Convergence of the associated Riccati recursion and stability properties of the steady-state filter are provided. (C) 2010 Elsevier Ltd. All rights reserved.

Long-term stability of hydrogen and organic acids production in an anaerobic fluidized-bed reactor using heat treated anaerobic sludge inoculum

This study evaluates the stability of hydrogen and organic acids production in an anaerobic fluidized-bed reactor (AFBR) that contains expanded clay (2.8-3.35 mm in diameter) as a support medium and is operated on a long-term basis. The reactor was inoculated with thermally pre-treated anaerobic sludge and operated with decreasing hydraulic retention time (HRT), from 8 h to 1 h, at a controlled temperature of 30 degrees C and a pH of about 3.8. Glucose (2000 mg L(-1)) was used as the substrate, generating conversion rates of 92-98%. Decreasing the HRT from 8 h to 1 h led to an increase in average hydrogen-production rates, with a maximum value of 1.28 L h(-1) L(-1) for an HRT of 1 h. In general, hydrogen yield production increased as HRT decreased, reaching 2.29 mol of H(2)/mol glucose at an HRT of 2 h and yielding a maximum hydrogen content of 37% in the biogas. No methane was detected in the biogas throughout the period of operation. The main soluble metabolites (SMP) were acetic acid (46.94-53.84% of SMP) and butyric acid (34.51-42.16% of SMP), with less than 15.49% ethanol. The steady performance of the AFBR may be attributed to adequate thermal treatment of the inoculum, the selection of a suitable support medium for microbial adhesion, and the choice of satisfactory environmental conditions imposed on the system. The results show that stable hydrogen production and organic acids production were maintained in the AFBR over a period of 178 days. (C) 2009 International Association for Hydrogen Energy. Published by Elsevier Ltd. All rights reserved.

Unsaturated Infinite Slope Stability Considering Surface Flux Conditions

A slope stability model is derived for an infinite slope subjected to unsaturated infiltration flow above a phreatic surface. Closed form steady state solutions are derived for the matric suction and degree of saturation profiles. Soil unit weight, consistent with the degree of saturation profile, is also directly calculated and introduced into the analyzes, resulting in closed-form solutions for typical soil parameters and an infinite series solution for arbitrary soil parameters. The solutions are coupled with the infinite slope stability equations to establish a fully realized safety factor function. In general, consideration of soil suction results in higher factor of safety. The increase in shear strength due to the inclusion of soil suction is analogous to making an addition to the cohesion, which, of course, increases the factor of safety against sliding. However, for cohesive soils, the results show lower safety factors for slip surfaces approaching the phreatic surface compared to those produced by common safety factor calculations. The lower factor of safety is due to the increased soil unit weight considered in the matric suction model but not usually accounted for in practice wherein the soil is treated as dry above the phreatic surface. The developed model is verified with a published case study, correctly predicting stability under dry conditions and correctly predicting failure for a particular storm.

An EFG method for the nonlinear analysis of plates undergoing arbitrarily large deformations

The applicability of a meshfree approximation method, namely the EFG method, on fully geometrically exact analysis of plates is investigated. Based on a unified nonlinear theory of plates, which allows for arbitrarily large rotations and displacements, a Galerkin approximation via MLS functions is settled. A hybrid method of analysis is proposed, where the solution is obtained by the independent approximation of the generalized internal displacement fields and the generalized boundary tractions. A consistent linearization procedure is performed, resulting in a semi-definite generalized tangent stiffness matrix which, for hyperelastic materials and conservative loadings, is always symmetric (even for configurations far from the generalized equilibrium trajectory). Besides the total Lagrangian formulation, an updated version is also presented, which enables the treatment of rotations beyond the parameterization limit. An extension of the arc-length method that includes the generalized domain displacement fields, the generalized boundary tractions and the load parameter in the constraint equation of the hyper-ellipsis is proposed to solve the resulting nonlinear problem. Extending the hybrid-displacement formulation, a multi-region decomposition is proposed to handle complex geometries. A criterium for the classification of the equilibrium`s stability, based on the Bordered-Hessian matrix analysis, is suggested. Several numerical examples are presented, illustrating the effectiveness of the method. Differently from the standard finite element methods (FEM), the resulting solutions are (arbitrary) smooth generalized displacement and stress fields. (c) 2007 Elsevier Ltd. All rights reserved.

Floquet stability analysis of the flow around an oscillating cylinder

This investigative work is concerned with the flow around a circular cylinder submitted to forced transverse oscillations. The goal is to investigate how the transition to turbulence is initiated in the wake for cases with different Reynolds numbers (Re) and displacement amplitudes (A). For each Re the motion frequency is kept constant, close to the Strouhal number of the flow around a fixed cylinder at the same Re. Stability analysis of two-dimensional periodic flows around a forced-oscillating cylinder is carried out with respect to three-dimensional infinitesimal perturbations. The procedure consists of performing a Floquet type analysis of time-periodic base flows, computed using the spectral/hp element method. With the results of the Floquet calculations, considerations regarding the stability of the system are drawn, and the form of the instability at its onset is obtained. The critical Reynolds number is observed to change with the amplitude of oscillation. With respect to instabilities, unstable modes with the same symmetry as mode A of a fixed cylinder are observed; however, they present different wavelengths. Also, the instabilities observed for the oscillating cylinder are distinctively stronger in the braid shear layers. Other unstable modes similar to mode B are found. Quasi-periodic modes are observed in the 2S wake, and subharmonic mode occurrences are reported in P + S wakes. (C) 2009 Elsevier Ltd. All rights reserved.

Metodologia para análise da confiabilidade estrutural de escavações em rocha.

A methodology for rock-excavation structural-reliability analysis that uses Distinct Element Method numerical models is presented. The methodology solves the problem of the conventional numerical models that supply only punctual results and use fixed input parameters, without considering its statistical errors. The analysis of rock-excavation stability must consider uncertainties from geological variability, from uncertainty in the choice of mechanical behaviour hypothesis, and from uncertainties in parameters adopted in numerical model construction. These uncertainties can be analyzed in simple deterministic models, but a new methodology was developed for numerical models with results of several natures. The methodology is based on Monte Carlo simulations and uses principles of Paraconsistent Logic. It will be presented in the analysis of a final slope of a large-dimensioned surface mine.

Interface Energy Measurement of MgO and ZnO: Understanding the Thermodynamic Stability of Nanoparticles

Nanomaterials have triggered excitement in both fundamental science and technological applications in several fields However, the same characteristic high interface area that is responsible for their unique properties causes unconventional instability, often leading to local collapsing during application Thermodynamically, this can be attributed to an increased contribution of the interface to the free energy, activating phenomena such as sintering and grain growth The lack of reliable interface energy data has restricted the development of conceptual models to allow the control of nanoparticle stability on a thermodynamic basis. Here we introduce a novel and accessible methodology to measure interface energy of nanoparticles exploiting the heat released during sintering to establish a quantitative relation between the solid solid and solid vapor interface energies. We exploited this method in MgO and ZnO nanoparticles and determined that the ratio between the solid solid and solid vapor interface energy is 11 for MgO and 0.7 for ZnO. We then discuss that this ratio is responsible for a thermodynamic metastable state that may prevent collapsing of nanoparticles and, therefore, may be used as a tool to design long-term stable nanoparticles.

Electron theoretical investigation of the stability of the B2-TiFe compound

The metastable phase diagram of the BCC-based ordering equilibria in the Ti-Fe system has been calculated using a truncated cluster expansion, through the combination of FP-LAPW and cluster variation method (CVM) in the irregular tetrahedron cluster approximation. The results are compared with phenomenological CVM assessments of the system and suggest that the value for the experimental formation enthalpy of the B2-TiFe compound should be significantly more negative than the currently assessed value. (C) 2008 Elsevier B.V. All rights reserved.

Interface excess and polymorphic stability of nanosized zirconia-magnesia

Controlling the phase stability of ZrO2 nanoparticles is of major importance in the development of new ZrO2-based nanotechnologies. Because of the fact that in nanoparticles the surface accounts for a larger fraction of the total atoms, the relative phase stability can be controlled throughout the surface composition, which can be toned by surface excess of one of the components of the system., The objective of this work is to delineate a relationship between surface excess (or solid solution) of MgO relative to ZrO2 and the polymorphic stability of (ZrO2)(1-x) - (MgO), nanopowders, where 0.0 <= x <= 0.6. The nanopowders were prepared by a liquid precursor method at 500 degrees C and characterized by N-2 adsorption (BET), X-ray diffraction (XRD), X-Ray photoelectron spectroscopy (XPS), and Raman spectroscopy. For pure ZrO2 samples, both tetragonal and monoclinic polymorphs were detected, as expected considering the literature. For MgO molar fractions varying from 0.05 to 0.10, extensive solid solution could not be detected, and a ZrO2 surface energy reduction, caused by Mg surface excess detected by XPS, promoted tetragonal polymorph thermodynamic stabilization with relation to monoclinic. For MgO molar fractions higher than 0.10 and up to 0.40, Mg solid solution could be detected and induced cubic phase stabilization. MgO periclase was observed only at x = 0.6. A discussion based on the relationship between the surface excess, surface energy, and polymorph stability is presented.

Evaluation of the effect of ammonium carbamate on the stability of proteins

BACKGROUND: The use of the volatile salt ammonium carbamate in protein downstream processing has recently been proposed. The main advantage of using volatile salts is that they can be removed from precipitates and liquid effluents through pressure reduction or temperature increase. Although previous studies showed that ammonium carbamate is efficient as a precipitant agent, there was evidence of denaturation in some enzymes. In this work, the effect of ammonium carbamate on the stability of five enzymes was evaluated. RESULTS: Activity assays showed that alpha-amylase (1,4-alpha-D-glucan glucanohydrolase, EC 3.2.1.1), lysozyme (1,4-beta-N-acetylmuramoylhydrolase, EC 3.2.1.17) and lipase (triacyl glycerol acyl hydrolase, EC 3.1.1.3) did not undergo activity loss in ammonium carbamate solutions with concentrations from 1.0 to 5.0 mol kg(-1), whereas cellulase complex (1,4-(1,3 : 14)-beta-D-glucan 4-glucano-hydrolase, EC 3.2.1.4) and peroxidase (hydrogen peroxide oxidoreductase, EC 1.11.1.7) showed an average activity loss of 55% and 44%, respectively. Precipitation assays did not show enzyme denaturation or phase separation for alpha-amylase and lipase, while celullase and peroxidase precipitated with some activity reduction. Analysis of similar experiments with ammonium and sodium sulfate did not affect the activity of enzymes. CONCLUSION: Celullase and peroxidase were denatured by ammonium carbamate. While more systematic studies are not available, care must be taken in designing a protein precipitation with this salt. The results suggest that the generally accepted idea that salts that denature proteins tend to solubilize them does not hold for ammonium carbamate. (C) 2010 Society of Chemical Industry

Stable Model Predictive Control for Integrating Systems with Optimizing Targets

This paper concern the development of a stable model predictive controller (MPC) to be integrated with real time optimization (RTO) in the control structure of a process system with stable and integrating outputs. The real time process optimizer produces Optimal targets for the system inputs and for Outputs that Should be dynamically implemented by the MPC controller. This paper is based oil a previous work (Comput. Chem. Eng. 2005, 29, 1089) where a nominally stable MPC was proposed for systems with the conventional control approach where only the outputs have set points. This work is also based oil the work of Gonzalez et at. (J. Process Control 2009, 19, 110) where the zone control of stable systems is studied. The new control for is obtained by defining ail extended control objective that includes input targets and zone controller the outputs. Additional decision variables are also defined to increase the set of feasible solutions to the control problem. The hard constraints resulting from the cancellation of the integrating modes Lit the end of the control horizon are softened,, and the resulting control problem is made feasible to a large class of unknown disturbances and changes of the optimizing targets. The methods are illustrated with the simulated application of the proposed,approaches to a distillation column of the oil refining industry.

Enlarging the domain of attraction of stable MPC controllers, maintaining the output performance

This work presents an alternative way to formulate the stable Model Predictive Control (MPC) optimization problem that allows the enlargement of the domain of attraction, while preserving the controller performance. Based on the dual MPC that uses the null local controller, it proposed the inclusion of an appropriate set of slacked terminal constraints into the control problem. As a result, the domain of attraction is unlimited for the stable modes of the system, and the largest possible for the non-stable modes. Although this controller does not achieve local optimality, simulations show that the input and output performances may be comparable to the ones obtained with the dual MPC that uses the LQR as a local controller. (C) 2009 Elsevier Ltd. All rights reserved.

Infinite horizon MPC with non-minimal state space feedback

In the MPC literature, stability is usually assured under the assumption that the state is measured. Since the closed-loop system may be nonlinear because of the constraints, it is not possible to apply the separation principle to prove global stability for the Output feedback case. It is well known that, a nonlinear closed-loop system with the state estimated via an exponentially converging observer combined with a state feedback controller can be unstable even when the controller is stable. One alternative to overcome the state estimation problem is to adopt a non-minimal state space model, in which the states are represented by measured past inputs and outputs [P.C. Young, M.A. Behzadi, C.L. Wang, A. Chotai, Direct digital and adaptative control by input-output, state variable feedback pole assignment, International journal of Control 46 (1987) 1867-1881; C. Wang, P.C. Young, Direct digital control by input-output, state variable feedback: theoretical background, International journal of Control 47 (1988) 97-109]. In this case, no observer is needed since the state variables can be directly measured. However, an important disadvantage of this approach is that the realigned model is not of minimal order, which makes the infinite horizon approach to obtain nominal stability difficult to apply. Here, we propose a method to properly formulate an infinite horizon MPC based on the output-realigned model, which avoids the use of an observer and guarantees the closed loop stability. The simulation results show that, besides providing closed-loop stability for systems with integrating and stable modes, the proposed controller may have a better performance than those MPC controllers that make use of an observer to estimate the current states. (C) 2008 Elsevier Ltd. All rights reserved.

Robust model predictive control with zone control

Model predictive control (MPC) is usually implemented as a control strategy where the system outputs are controlled within specified zones, instead of fixed set points. One strategy to implement the zone control is by means of the selection of different weights for the output error in the control cost function. A disadvantage of this approach is that closed-loop stability cannot be guaranteed, as a different linear controller may be activated at each time step. A way to implement a stable zone control is by means of the use of an infinite horizon cost in which the set point is an additional variable of the control problem. In this case, the set point is restricted to remain inside the output zone and an appropriate output slack variable is included in the optimisation problem to assure the recursive feasibility of the control optimisation problem. Following this approach, a robust MPC is developed for the case of multi-model uncertainty of open-loop stable systems. The controller is devoted to maintain the outputs within their corresponding feasible zone, while reaching the desired optimal input target. Simulation of a process of the oil re. ning industry illustrates the performance of the proposed strategy.