Methods and Models for Accelerating Dynamic Simulation of Fluid Power Circuits

The objective of this dissertation is to improve the dynamic simulation of fluid power circuits. A fluid power circuit is a typical way to implement power transmission in mobile working machines, e.g. cranes, excavators etc. Dynamic simulation is an essential tool in developing controllability and energy-efficient solutions for mobile machines. Efficient dynamic simulation is the basic requirement for the real-time simulation. In the real-time simulation of fluid power circuits there exist numerical problems due to the software and methods used for modelling and integration. A simulation model of a fluid power circuit is typically created using differential and algebraic equations. Efficient numerical methods are required since differential equations must be solved in real time. Unfortunately, simulation software packages offer only a limited selection of numerical solvers. Numerical problems cause noise to the results, which in many cases leads the simulation run to fail. Mathematically the fluid power circuit models are stiff systems of ordinary differential equations. Numerical solution of the stiff systems can be improved by two alternative approaches. The first is to develop numerical solvers suitable for solving stiff systems. The second is to decrease the model stiffness itself by introducing models and algorithms that either decrease the highest eigenvalues or neglect them by introducing steady-state solutions of the stiff parts of the models. The thesis proposes novel methods using the latter approach. The study aims to develop practical methods usable in dynamic simulation of fluid power circuits using explicit fixed-step integration algorithms. In this thesis, twomechanisms whichmake the systemstiff are studied. These are the pressure drop approaching zero in the turbulent orifice model and the volume approaching zero in the equation of pressure build-up. These are the critical areas to which alternative methods for modelling and numerical simulation are proposed. Generally, in hydraulic power transmission systems the orifice flow is clearly in the turbulent area. The flow becomes laminar as the pressure drop over the orifice approaches zero only in rare situations. These are e.g. when a valve is closed, or an actuator is driven against an end stopper, or external force makes actuator to switch its direction during operation. This means that in terms of accuracy, the description of laminar flow is not necessary. But, unfortunately, when a purely turbulent description of the orifice is used, numerical problems occur when the pressure drop comes close to zero since the first derivative of flow with respect to the pressure drop approaches infinity when the pressure drop approaches zero. Furthermore, the second derivative becomes discontinuous, which causes numerical noise and an infinitely small integration step when a variable step integrator is used. A numerically efficient model for the orifice flow is proposed using a cubic spline function to describe the flow in the laminar and transition areas. Parameters for the cubic spline function are selected such that its first derivative is equal to the first derivative of the pure turbulent orifice flow model in the boundary condition. In the dynamic simulation of fluid power circuits, a tradeoff exists between accuracy and calculation speed. This investigation is made for the two-regime flow orifice model. Especially inside of many types of valves, as well as between them, there exist very small volumes. The integration of pressures in small fluid volumes causes numerical problems in fluid power circuit simulation. Particularly in realtime simulation, these numerical problems are a great weakness. The system stiffness approaches infinity as the fluid volume approaches zero. If fixed step explicit algorithms for solving ordinary differential equations (ODE) are used, the system stability would easily be lost when integrating pressures in small volumes. To solve the problem caused by small fluid volumes, a pseudo-dynamic solver is proposed. Instead of integration of the pressure in a small volume, the pressure is solved as a steady-state pressure created in a separate cascade loop by numerical integration. The hydraulic capacitance V/Be of the parts of the circuit whose pressures are solved by the pseudo-dynamic method should be orders of magnitude smaller than that of those partswhose pressures are integrated. The key advantage of this novel method is that the numerical problems caused by the small volumes are completely avoided. Also, the method is freely applicable regardless of the integration routine applied. The superiority of both above-mentioned methods is that they are suited for use together with the semi-empirical modelling method which necessarily does not require any geometrical data of the valves and actuators to be modelled. In this modelling method, most of the needed component information can be taken from the manufacturer’s nominal graphs. This thesis introduces the methods and shows several numerical examples to demonstrate how the proposed methods improve the dynamic simulation of various hydraulic circuits.

Numerical study on mixed convection in a horizontal flow past a square porous cavity using UNIFAES scheme

Mixed convection on the flow past a heated length and past a porous cavity located in a horizontal wall bounding a saturated porous medium is numerically simulated. The cavity is heated from below. The steady-state regime is studied for several intensities of the buoyancy effects due to temperature variations. The influences of Péclet and Rayleigh numbers on the flow pattern and the temperature distributions are examined. Local and global Nusselt numbers are reported for the heated surface. The convective-diffusive fluxes at the volume boundaries are represented using the UNIFAES, Unified Finite Approach Exponential-type Scheme, with the Power-Law approximation to reduce the computing time. The conditions established by Rivas for the quadratic order of accuracy of the central differencing to be maintained in irregular grids are shown to be extensible to other quadratic schemes, including UNIFAES, so that accuracy estimates could be obtained.

Novel technology for preparation of optically active chemicals

Asymmetric synthesis using modified heterogeneous catalysts has gained lots of interest in the production of optically pure chemicals, such as pharmaceuticals, nutraceuticals, fragrances and agrochemicals. Heterogeneous modified catalysts capable of inducing high enantioselectivities are preferred in industrial scale due to their superior separation and handling properties. The topic has been intensively investigated both in industry and academia. The enantioselective hydrogenation of ethyl benzoylformate (EBF) to (R)-ethyl mandelate over (-)-cinchonidine (CD)-modified Pt/Al2O3 catalyst in a laboratory-scale semi-batch reactor was studied as a function of modifier concentration, reaction temperature, stirring rate and catalyst particle size. The main product was always (R)-ethyl mandelate while small amounts of (S)-ethyl mandelate were obtained as by product. The kinetic results showed higher enantioselectivity and lower initial rates approaching asymptotically to a constant value as the amount of modifier was increased. Additionally, catalyst deactivation due to presence of impurities in the feed was prominent in some cases; therefore activated carbon was used as a cleaning agent of the raw material to remove impurities prior to catalyst addition. Detailed characterizations methods (SEM, EDX, TPR, BET, chemisorption, particle size distribution) of the catalysts were carried out. Solvent effects were also studied in the semi-batch reactor. Solvents with dielectric constant (e) between 2 and 25 were applied. The enantiomeric excess (ee) increased with an increase of the dielectric coefficient up to a maximum followed by a nonlinear decrease. A kinetic model was proposed for the enantioselectivity dependence on the dielectric constant based on the Kirkwood treatment. The non-linear dependence of ee on (e) successfully described the variation of ee in different solvents. Systematic kinetic experiments were carried out in the semi-batch reactor. Toluene was used as a solvent. Based on these results, a kinetic model based on the assumption of different number of sites was developed. Density functional theory calculations were applied to study the energetics of the EBF adsorption on pure Pt(1 1 1). The hydrogenation rate constants were determined along with the adsorption parameters by non-linear regression analysis. A comparison between the model and the experimental data revealed a very good correspondence. Transient experiments in a fixed-bed reactor were also carried out in this work. The results demonstrated that continuous enantioselective hydrogenation of EBF in hexane/2-propanol 90/10 (v/v) is possible and that continuous feeding of (-)-cinchonidine is needed to maintain a high steady-state enantioselectivity. The catalyst showed a good stability and high enantioselectivity was achieved in the fixed-bed reactor. Chromatographic separation of (R)- and (S)-ethyl mandelate originating from the continuous reactor was investigated. A commercial column filled with a chiral resin was chosen as a perspective preparative-scale adsorbent. Since the adsorption equilibrium isotherms were linear within the entire investigated range of concentrations, they were determined by pulse experiments for the isomers present in a post-reaction mixture. Breakthrough curves were measured and described successfully by the dispersive plug flow model with a linear driving force approximation. The focus of this research project was the development of a new integrated production concept of optically active chemicals by combining heterogeneous catalysis and chromatographic separation technology. The proposed work is fundamental research in advanced process technology aiming to improve efficiency and enable clean and environmentally benign production of enantiomeric pure chemicals.

Synechococcus sp. PCC 7002 CpcB lyase null mutations alter phycocyanin chromophore function and, consequently, affect the redistribution of excitation energy via the light state transition /

Phycobilisomes are the major light harvesting complexes for cyanobacteria and phycocyanin is the primary phycobiliprotein of the phycobilisome rod. The phycocyanobilin lyases responsible for chromophorylating the phycocyanin p subunit (CpcB) have been recently identified in the cyanobacterium Synechococcus sp. PCC 7002. Surprisingly, mutants missing the CpcB lyases were nevertheless capable of producing pigmented phycocyanin. 10K absorbance measurements revealed that the energy states of the p phycocyanin chromophores were only subtly shifted; however, 77K steady state fluorescence emission spectroscopy showed excitation energy transfer involving the targeted chromophores to be highly disrupted. Such evidence suggests that phycobilin orientation within the binding domain is specifically modified. We hypothesized that alternate, less specific lyases are able to act on the p binding sites. A phycocyanin linker-polypeptide deficient mutant was similarly characterized. The light state transition, a short term adaptation of the photosynthetic light harvesting apparatus resulting in the redistribution of excitation energy among the photo systems, was shown to be dominated by the reallocation of phycocyanin-absorbed excitation energy. Treatment with a high M phosphate buffer effectively prevented the redistribution of both chlorophyll a- and phycobilisome- absorbed excitation energy, suggesting that the two effects are not strictly independent. The mutant strains required a larger redistribution of excitation energy between light states, perhaps to compensate for their loss in phycobilisome antenna function.

FastAero – A Precorrected FFT – Fast Multipole Tree Steady and Unsteady Potential Flow Solver

In this paper a precorrected FFT-Fast Multipole Tree (pFFT-FMT) method for solving the potential flow around arbitrary three dimensional bodies is presented. The method takes advantage of the efficiency of the pFFT and FMT algorithms to facilitate more demanding computations such as automatic wake generation and hands-off steady and unsteady aerodynamic simulations. The velocity potential on the body surfaces and in the domain is determined using a pFFT Boundary Element Method (BEM) approach based on the Green’s Theorem Boundary Integral Equation. The vorticity trailing all lifting surfaces in the domain is represented using a Fast Multipole Tree, time advected, vortex participle method. Some simple steady state flow solutions are performed to demonstrate the basic capabilities of the solver. Although this paper focuses primarily on steady state solutions, it should be noted that this approach is designed to be a robust and efficient unsteady potential flow simulation tool, useful for rapid computational prototyping.

Solid-state fermentation: a continuous process for fungal tannase production

Truly continuous solid-state fermentations with operating times of 2-3 weeks were conducted in a prototype bioreactor for the production of fungal (Penicillium glabrum) tannase from a tannin-containing model substrate. Substantial quantities of the enzyme were synthesized throughout the operating periods and (imperfect) steady-state conditions seemed to be achieved soon after start-up of the fermentations. This demonstrated for the first time the possibility of conducting solid-state fermentations in the continuous mode and with a constant noninoculated feed. The operating variables and fermentation conditions in the bioreactor were sufficiently well predicted for the basic reinoculation concept to succeed. However, an incomplete understanding of the microbial mechanisms, the experimental system, and their interaction indicated the need for more research in this novel area of solid-state fermentation. (C) 2004 Wiley Periodicals, Inc.

An LMS style variable tap-length algorithm for structure adaptation

Searching for the optimum tap-length that best balances the complexity and steady-state performance of an adaptive filter has attracted attention recently. Among existing algorithms that can be found in the literature, two of which, namely the segmented filter (SF) and gradient descent (GD) algorithms, are of particular interest as they can search for the optimum tap-length quickly. In this paper, at first, we carefully compare the SF and GD algorithms and show that the two algorithms are equivalent in performance under some constraints, but each has advantages/disadvantages relative to the other. Then, we propose an improved variable tap-length algorithm using the concept of the pseudo fractional tap-length (FT). Updating the tap-length with instantaneous errors in a style similar to that used in the stochastic gradient [or least mean squares (LMS)] algorithm, the proposed FT algorithm not only retains the advantages from both the SF and the GD algorithms but also has significantly less complexity than existing algorithms. Both performance analysis and numerical simulations are given to verify the new proposed algorithm.

Modeling electrocortical activity through improved local approximations of integral neural field equations

Neural field models of firing rate activity typically take the form of integral equations with space-dependent axonal delays. Under natural assumptions on the synaptic connectivity we show how one can derive an equivalent partial differential equation (PDE) model that properly treats the axonal delay terms of the integral formulation. Our analysis avoids the so-called long-wavelength approximation that has previously been used to formulate PDE models for neural activity in two spatial dimensions. Direct numerical simulations of this PDE model show instabilities of the homogeneous steady state that are in full agreement with a Turing instability analysis of the original integral model. We discuss the benefits of such a local model and its usefulness in modeling electrocortical activity. In particular, we are able to treat “patchy” connections, whereby a homogeneous and isotropic system is modulated in a spatially periodic fashion. In this case the emergence of a “lattice-directed” traveling wave predicted by a linear instability analysis is confirmed by the numerical simulation of an appropriate set of coupled PDEs.

Coalescence of particles by differential sedimentation

We consider a three dimensional system consisting of a large number of small spherical particles, distributed in a range of sizes and heights (with uniform distribution in the horizontal direction). Particles move vertically at a size-dependent terminal velocity. They are either allowed to merge whenever they cross or there is a size ratio criterion enforced to account for collision efficiency. Such a system may be described, in mean field approximation, by the Smoluchowski kinetic equation with a differential sedimentation kernel. We obtain self-similar steady-state and time-dependent solutions to the kinetic equation, using methods borrowed from weak turbulence theory. Analytical results are compared with direct numerical simulations (DNS) of moving and merging particles, and a good agreement is found.

Excited-State Dynamics in fac-[Re(CO)(3)(Me(4)phen)(L)](+)

Excited-state dynamics in fac-[Re(CO)(3)(Me(4)phen)(cis-L)](+) (Me(4)phen = 3,4,7,8-tetramethyl-1,10-phenanthroline, L = 4-styrylpyridine (stpy) or 1,2-bis(4-pyridyl)ethylene (bpe)) were investigated by steady-state and time-resolved techniques. A complex equilibrium among three closely lying excited states, 3IL(cis-L), (3)MLCT(Re -> me4phen), and (3)IL(Me4phen), has been established. Under UV irradiation, cis-to-trans isomerization of coordinated cis-L is observed with a quantum yield of 0.15 in acetonitrile solutions. This photoreaction competes with radiative decay from (3)MLCT(Re -> Me4phen) and (3)IL(Me4phen) excited states, leading to a decrease in the emission quantum yield relative to the nonisomerizable complex fac-[Re(CO)(3)(Me(4)phen)(bpa)](+) (bpa = 1,2-bis(4-pyridyl)ethane). From temperature-dependent time-resolved emission measurements in solution and in poly(methyl methacrylate) (PMMA) films, energy barriers (Delta E(a)) for interconversion between (3)MLCT(Re -> me4Phen) and (3)IL(Me4phen) emitting states were determined. For L = cis-stpy, Delta E(a) = 11 (920 cm(-1)) and 15 kJ mol(-1) (1254 cm(-1)) in 5:4 propionitrile/butyronitrile and PMMA, respectively. For L = cis-bpe, Delta E(a) = 13 kJ mol(-1) (1087 cm(-1)) in 5:4 propionitrile/butyronitrile. These energy barriers are sufficient to decrease the rate constant for internal conversion from higher-lying (3)IL(me4phen) state to (3)MLCT(Re -> Me4phen), k(i) congruent to 10(6) s(-1). The decrease in rate allows for the observation of intraligand phosphorescence, even in fluid medium at room temperature. Our results provide additional insight into the role of energy gap and excited-state dynamics on the photochemical and photophysical properties of Re(I) polypyridyl complexes.

Dynamic whole system testing of combined renewable heating systems : the current state of the art

Objective: For the evaluation of the energetic performance of combined renewable heating systems that supply space heat and domestic hot water for single family houses, dynamic behaviour, component interactions, and control of the system play a crucial role and should be included in test methods. Methods: New dynamic whole system test methods were developed based on “hardware in the loop” concepts. Three similar approaches are described and their differences are discussed. The methods were applied for testing solar thermal systems in combination with fossil fuel boilers (heating oil and natural gas), biomass boilers, and/or heat pumps. Results: All three methods were able to show the performance of combined heating systems under transient operating conditions. The methods often detected unexpected behaviour of the tested system that cannot be detected based on steady state performance tests that are usually applied to single components. Conclusion: Further work will be needed to harmonize the different test methods in order to reach comparable results between the different laboratories. Practice implications: A harmonized approach for whole system tests may lead to new test standards and improve the accuracy of performance prediction as well as reduce the need for field tests.

Otimização topológica de estruturas termoelásticas tridimensionais

This work presents an optimization technique based on structural topology optimization methods, TOM, designed to solve problems of thermoelasticity 3D. The presented approach is based on the adjoint method of sensitivity analysis unified design and is intended to loosely coupled thermomechanical problems. The technique makes use of analytical expressions of sensitivities, enabling a reduction in the computational cost through the use of a coupled field adjoint equation, defined in terms the of temperature and displacement fields. The TOM used is based on the material aproach. Thus, to make the domain is composed of a continuous distribution of material, enabling the use of classical models in nonlinear programming optimization problem, the microstructure is considered as a porous medium and its constitutive equation is a function only of the homogenized relative density of the material. In this approach, the actual properties of materials with intermediate densities are penalized based on an artificial microstructure model based on the SIMP (Solid Isotropic Material with Penalty). To circumvent problems chessboard and reduce dependence on layout in relation to the final optimal initial mesh, caused by problems of numerical instability, restrictions on components of the gradient of relative densities were applied. The optimization problem is solved by applying the augmented Lagrangian method, the solution being obtained by applying the finite element method of Galerkin, the process of approximation using the finite element Tetra4. This element has the ability to interpolate both the relative density and the displacement components and temperature. As for the definition of the problem, the heat load is assumed in steady state, i.e., the effects of conduction and convection of heat does not vary with time. The mechanical load is assumed static and distributed

Kinetics of interface state-limited hole injection in alpha-naphthylphenylbiphenyl diamine (alpha-NPD) thin layers

Injection-limited operation is identified in thin-film, alpha-NPD-based diodes. A detailed model for the impedance of the injection process is provided which considers the kinetics of filling/releasing of interface states as the key factor behind the injection mechanism. The injection model is able to simultaneously account for the steady-state, current-voltage (J-V) characteristics and impedance response. and is based on the sequential injection of holes mediated by energetically distributed surface states at the metal-organic interface. The model takes into account the vacuum level offset caused by the interface dipole, along with the partial shift of the interface level distribution with bias voltage. This approach connects the low-frequency (similar to 1 Hz) capacitance spectra, which exhibits a transition between positive to negative values, to the change in the occupancy of interface states with voltage. Simulations based on the model allow to derive the density of interface states effectively intervening in the carrier injection (similar to 5 x 10(12) cm(-2)), which exhibit a Gaussian-like distribution. A kinetically determined hole barrier is calculated at levels located similar to 0.4 eV below the contact work function. (C) 2008 Elsevier B.V. All rights reserved.

UNUSUALLY STRONG ATTRACTION IN THE PRESENCE OF CONTINUUM BOUND-STATE

The result of few-particle ground-state calculation employing a two-particle nonlocal potential supporting a continuum bound state in addition to a negative-energy bound state has occasionally revealed unusually strong attraction in producing a very strongly bound ground state. In the presence of the continuum bound state the difference of phase shift between zero and infinite energies has an extra jump of pi as in the presence of an additional bound state. The wave function of the continuum bound state is identical with that of a strongly bound negative-energy state, which leads us to postulate a pseudo bound state in the two-particle system in order to explain the unexpected attraction. The role of the Pauli forbidden states is expected to be similar to these pseudo states.

Pharmacokinetics of dexmedetomidine administered intravenously in isofluraneanesthetized cats

Objective-To determine the pharmacokinetics of dexmedetomidine administered as a short-duration IV infusion in isoflurane-anesthetized cats. Animals-6 healthy adult domestic female cats. Procedures-Dexmedetomidine hydrochloride was injected IV (10 μg/kg over 5 minutes [rate, 2 μg/kg/min]) in isoflurane-anesthetized cats. Blood samples were obtained immediately prior to and at 1, 2, 5, 6, 7, 10, 15, 30, 60, 90, 120, 240, and 480 minutes following the start of the IV infusion. Collected blood samples were transferred to tubes containing EDTA, immediately placed on ice, and then centrifuged at 3,901 X g for 10 minutes at 4°C. The plasma was harvested and stored at -20°C until analyzed. Plasma dexmedetomidine concentrations were determined by means of liquid chromatography-mass spectrometry. Dexmedetomidine plasma concentration-time data were fitted to compartmental models. Results-A 2-compartment model with input in and elimination from the central compartment best described the disposition of dexmedetomidine administered via short-duration IV infusion in isoflurane-anesthetized cats. Weighted mean ± SEM apparent volume of distribution of the central compartment and apparent volume of distribution at steady-state were 402 ± 47 mL/kg and 1,701 ± 200 mL/kg, respectively; clearance and terminal half-life (harmonic mean ± jackknife pseudo-SD) were 6.3 ± 2.8 mL/min/kg and 198 ± 75 minutes, respectively. The area under the plasma concentration curve and maximal plasma concentration were 1,061 ± 292 min·ng/mL and 17.6 ± 1.8 ng/mL, respectively. Conclusions and Clinical Relevance-Disposition of dexmedetomidine administered via short-duration IV infusion in isoflurane-anesthetized cats was characterized by a moderate clearance and a long terminal half-life.