Solution of a heterogeneous modeling of a horizontal-flow anaerobic immobilized biomass (HAIB) reactor by the sequencing method

A modeling study was completed to develop a methodology that combines the sequencing and finite difference methods for the simulation of a heterogeneous model of a tubular reactor applied in the treatment of wastewater. The system included a liquid phase (convection diffusion transport) and a solid phase (diffusion reaction) that was obtained by completing a mass balance in the reactor and in the particle, respectively. The model was solved using a pilot-scale horizontal-flow anaerobic immobilized biomass (HAIB) reactor to treat domestic sewage, with the concentration results compared with the experimental data. A comparison of the behavior of the liquid phase concentration profile and the experimental results indicated that both the numerical methods offer a good description of the behavior of the concentration along the reactor. The advantage of the sequencing method over the finite difference method is that it is easier to apply and requires less computational time to model the dynamic simulation of outlet response of HAIB.

Influence of the Solvent Evaporation Rate on the Crystalline Phases of Solution-Cast Poly(Vinylidene Fluoride) Films

The influence of the solvent-evaporation rate on the formation of of. and P crystalline phases in solution-cast poly(vinylidene fluoride) (PVDF) films was systematically investigated. Films were crystallized from PVDF/N,N-dimethylformamide solutions with concentrations of 2.5, 5.0, 10, and 20 wt % at different temperatures. During crystallization, the solvent evaporation rate was monitored in situ by means of a semianalytic balance. With this system, it was possible to determine the evaporation rate for different concentrations and temperatures of the solution under specific ambient conditions (pressure, temperature, and humidity). Fourier-Transform InfraRed spectroscopy with Attenuated Total Reflectance revealed the P-phase content in the PVDF films and its dependence on previous evaporation rates. Based on the relation between the evaporation rate and the PVDF phase composition, a consistent explanation for the different amounts of P phase observed at the upper and lower sample surfaces is achieved. Furthermore, the role of the sample thickness has also been studied. The experimental results show that not only the temperature but also the evaporation rate have to be controlled to obtain the desired crystalline phases in solution-cast PVDF films. (C) 2009 Wiley Periodicals, Inc. J Appl Polym Sci 116: 785-791, 2010

A hybrid-mixed finite element formulation for the geometrically exact analysis of three-dimensional framed structures

This paper addresses the development of a hybrid-mixed finite element formulation for the quasi-static geometrically exact analysis of three-dimensional framed structures with linear elastic behavior. The formulation is based on a modified principle of stationary total complementary energy, involving, as independent variables, the generalized vectors of stress-resultants and displacements and, in addition, a set of Lagrange multipliers defined on the element boundaries. The finite element discretization scheme adopted within the framework of the proposed formulation leads to numerical solutions that strongly satisfy the equilibrium differential equations in the elements, as well as the equilibrium boundary conditions. This formulation consists, therefore, in a true equilibrium formulation for large displacements and rotations in space. Furthermore, this formulation is objective, as it ensures invariance of the strain measures under superposed rigid body rotations, and is not affected by the so-called shear-locking phenomenon. Also, the proposed formulation produces numerical solutions which are independent of the path of deformation. To validate and assess the accuracy of the proposed formulation, some benchmark problems are analyzed and their solutions compared with those obtained using the standard two-node displacement/ rotation-based formulation.

Shell curvature as an initial deformation: A geometrically exact finite element approach

An alternative approach for the analysis of arbitrarily curved shells is developed in this paper based on the idea of initial deformations. By `alternative` we mean that neither differential geometry nor the concept of degeneration is invoked here to describe the shell surface. We begin with a flat reference configuration for the shell mid-surface, after which the initial (curved) geometry is mapped as a stress-free deformation from the plane position. The actual motion of the shell takes place only after this initial mapping. In contrast to classical works in the literature, this strategy enables the use of only orthogonal frames within the theory and therefore objects such as Christoffel symbols, the second fundamental form or three-dimensional degenerated solids do not enter the formulation. Furthermore, the issue of physical components of tensors does not appear. Another important aspect (but not exclusive of our scheme) is the possibility to describe exactly the initial geometry. The model is kinematically exact, encompasses finite strains in a totally consistent manner and is here discretized under the light of the finite element method (although implementation via mesh-free techniques is also possible). Assessment is made by means of several numerical simulations. Copyright (C) 2009 John Wiley & Sons, Ltd.

An exact conserving algorithm for nonlinear dynamics with rotational DOFs and general hyperelasticity. Part 2: shells

Following the approach developed for rods in Part 1 of this paper (Pimenta et al. in Comput. Mech. 42:715-732, 2008), this work presents a fully conserving algorithm for the integration of the equations of motion in nonlinear shell dynamics. We begin with a re-parameterization of the rotation field in terms of the so-called Rodrigues rotation vector, allowing for an extremely simple update of the rotational variables within the scheme. The weak form is constructed via non-orthogonal projection, the time-collocation of which ensures exact conservation of momentum and total energy in the absence of external forces. Appealing is the fact that general hyperelastic materials (and not only materials with quadratic potentials) are permitted in a totally consistent way. Spatial discretization is performed using the finite element method and the robust performance of the scheme is demonstrated by means of numerical examples.

An exact conserving algorithm for nonlinear dynamics with rotational DOFs and general hyperelasticity. Part 1: Rods

A fully conserving algorithm is developed in this paper for the integration of the equations of motion in nonlinear rod dynamics. The starting point is a re-parameterization of the rotation field in terms of the so-called Rodrigues rotation vector, which results in an extremely simple update of the rotational variables. The weak form is constructed with a non-orthogonal projection corresponding to the application of the virtual power theorem. Together with an appropriate time-collocation, it ensures exact conservation of momentum and total energy in the absence of external forces. Appealing is the fact that nonlinear hyperelastic materials (and not only materials with quadratic potentials) are permitted without any prejudice on the conservation properties. Spatial discretization is performed via the finite element method and the performance of the scheme is assessed by means of several numerical simulations.

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.

A comparison of dynamic and quasi-static results for the flow towards a well

The present analysis takes into account the acceleration term in the differential equation of motion to obtain exact dynamic solutions concerning the groundwater flow towards a well in a confined aquifer. The results show that the error contained in the traditional quasi-static solution is very small in typical situations.

Hybrid and multi-field variational principles for geometrically exact three-dimensional beams

This paper addresses the development of several alternative novel hybrid/multi-field variational formulations of the geometrically exact three-dimensional elastostatic beam boundary-value problem. In the framework of the complementary energy-based formulations, a Legendre transformation is used to introduce the complementary energy density in the variational statements as a function of stresses only. The corresponding variational principles are shown to feature stationarity within the framework of the boundary-value problem. Both weak and linearized weak forms of the principles are presented. The main features of the principles are highlighted, giving special emphasis to their relationships from both theoretical and computational standpoints. (C) 2010 Elsevier Ltd. All rights reserved.

Numerical solution of the two-phase expansion of a metastable flashing liquid jet using the dispersion-controlled dissipative scheme

This paper presents a study of the stationary phenomenon of superheated or metastable liquid jets, flashing into a two-dimensional axisymmetric domain, while in the two-phase region. In general, the phenomenon starts off when a high-pressure, high-temperature liquid jet emerges from a small nozzle or orifice expanding into a low-pressure chamber, below its saturation pressure taken at the injection temperature. As the process evolves, crossing the saturation curve, one observes that the fluid remains in the liquid phase reaching a superheated condition. Then, the liquid undergoes an abrupt phase change by means of an oblique evaporation wave. Across this phase change the superheated liquid becomes a two-phase high-speed mixture in various directions, expanding to supersonic velocities. In order to reach the downstream pressure, the supersonic fluid continues to expand, crossing a complex bow shock wave. The balance equations that govern the phenomenon are mass conservation, momentum conservation, and energy conservation, plus an equation-of-state for the substance. A false-transient model is implemented using the shock capturing scheme: dispersion-controlled dissipative (DCD), which was used to calculate the flow conditions as the steady-state condition is reached. Numerical results with computational code DCD-2D vI have been analyzed. Copyright (C) 2009 John Wiley & Sons, Ltd.

Surface tension of flotation solution and its influence on the selectivity of the separation between apatite and gangue minerals

This work presents and discusses the influence of the surface tension (gamma(LV)) of methanol/water mixtures on the flotation response of apatite versus gangue minerals conditioned with flotation reagents (75 g/t cornstarch and 100 g/t Berol 867) at pH 10.6. Berol 867 is a collector composed of sodium alkyl sarcosinate plus nonionic surfactant. The highest Schulz efficiency of separation (recovery of apatite minus recovery of gangue) was achieved at approximate to 51.0 mN/m. The critical surface tension of wettability (gamma(C)) of apatite was found to occur at 34.7 mN/m when determined by means of gamma flotation experiments, , and it was 33.9 mN/m when determined by Zisman`s approach.

Post-failure behaviour of impulsively loaded circular plates

The dynamic plastic response of a simply supported circular plate is analysed. Emphasis is given to the plate behaviour after it has broken free from the supports due to a local material failure. The theoretical rigid plastic analysis predicts various features of the response such as the time to failure, residual kinetic energy and the critical velocity at failure. The residual kinetic energy of the plate could be significant enough to cause secondary impact damage. It is shown that the shape of the plate changes after breaking free from the supports, which is important for forensic investigations. The solution for various cases were proven to be exact in the context of the upper and lower bounds theorems of the theory of plasticity. (C) 2009 Elsevier Ltd. All rights reserved.

Evidences of the evolution from solid solution to surface segregation in Ni-doped SnO(2) nanoparticles using Raman spectroscopy

Ni-doped SnO(2) nanoparticles, promising for gas-sensing applications, have been synthesized by a polymer precursor method. X-ray diffraction (XRD) and transmission electron microscopy (TEM) data analyses indicate the exclusive formation of nanosized particles with rutile-type phase (tetragonal SnO(2)) for Ni contents below 10 mol%. The mean crystallite size shows a progressive reduction with the Ni content. Room-temperature Raman spectra of Ni-doped SnO(2) nanoparticles show the presence of Raman active modes and modes activated by size effects. From the evolution of the A(1g) mode with the Ni content, a solubility limit at similar to 2 mol% was estimated. Below that content, Raman results are consistent with the occurrence of solid solution (ss) and surface segregation (seg.) of Ni ions. Above similar to 2 mol% Ni, the redshift of A(1g) mode suggests that the surface segregation of Ni ions takes place. Disorder-activated bands were determined and their integrated intensity evolution with the Ni content suggest that the solid-solution regime favors the increase of disorder; meanwhile, that disorder becomes weaker as the Ni content is increased. Copyright (C) 2010 John Wiley & Sons, Ltd.

Selective corrosion of 550 degrees C aged Cu10Ni-3Al-1.3Fe alloy in NaCl aqueous solution

This work studied the electrochemical behavior of a solution treated or 550 degrees C aged Cu10Ni-3Al-1.3Fe alloy, in 0.01 M NaCl aqueous solution, through potentiodynamic polarization in both stagnant condition or under erosion process. Results showed the occurrence of a passivity break potential (E(pb)), related to the beginning of the denickelification process, which occurred as a localized attack under stagnant electrolyte. Under erosion conditions localized denickelification was not observed, despite of the presence of E(pb). This could indicate that selective corrosion of Ni, which caused the observed E(pb), occurred as a dissolution-redeposition process, with removal of the Cu deposits during erosion process. (C) 2011 Elsevier Ltd. All rights reserved.

Degradation of Poly(ethylene glycol) in Aqueous Solution by Photo-Fenton and H(2)O(2)/UV Processes

This article reports experimental results obtained in a laboratory-scale photochemical reactor on the photodegradation of poly(ethylene glycol) (PEG) in aqueous solutions by means of the photo-Fenton and H(2)O(2)/UV processes. Dilute water solutions of PEG were fed to a batch reactor, mixed with pertinent reactants, and allowed to react under different conditions. Reaction progress was evaluated by sampling and analyzing the concentration of the total organic carbon (TOC) in solution as a function of the reaction time. Organic acids formed during oxidation were determined by HPLC analyses. The main acids detected in both processes were acetic and formic. Glycolic acid was detected only in the photo-Fenton process, and malonic acid was detected only in the H(2)O(2)/UV treatment, indicating that different reaction paths occur in these processes. The characteristics of both processes are discussed, based on the evolution of the TOC-time curves and the concentration profiles of the monitored organic acids. The experimental results constitute a contribution to the design of industrial processes for the treatment of wastewaters containing soluble polymers with similar properties.