New Algorithm for Computing LQ Suboptimal Output Feedback Gains of Decentralized Control Systems

In this paper, the problem of computing the suboptimal output feedback gains of decentralized control systems is investigated. First, the problem is formulated. Then, the gradient matrices based on the index function are derived and a new algorithm is established based on some nice properties. This algorithm shows that a suboptimal gain can be computed by solving several ordinary differential equations (ODEs). In order to find an initial condition for the ODEs, an algorithm for finding a stabilizing output feedback gain is exploited, and the convergence of this algorithm is discussed. Finally, an example is given to illustrate the proposed algorithm.     

Long Time Behavior of the Cahn-Hilliard Equation with Irregular Potentials and Dynamic Boundary Conditions

The Cahn-Hilliard equation with irregular potentials and dynamic boundary conditions is considered.The existence of the global attractor is proved and the long time behavior of the trajectories,namely,the convergence to steady states,is studied.     

Fluid membranes in hydrodynamic flow fields: Formalism and an application to fluctuating quasispherical vesicles in shear flow

The dynamics of a single fluid bilayer membrane in an external hydrodynamic flow field is considered. The deterministic equation of motion for the configuration is derived taking into account both viscous dissipation in the surrounding liquid and local incompressibility of the membrane. For quasi-spherical vesicles in shear flow, thermal fluctuations can be incorporated in a Langevin-type equation of motion for the deformation amplitudes. The solution to this equation shows an overdamped oscillatory approach to a stationary tanktreading shape. Inclination angle and ellipticity of the contour are determined as a function of excess area and shear rate. Comparisons to numerical results and experiments are discussed. Received 20 August 1998     

Influence of the surface deformability and variable viscosity on buoyant - thermocapillary instability in a liquid layer

The primary stationary and oscillatory Bénard-Marangoni instability is investigated in a fluid layer of infinite horizontal extent, bounded below by a rigid plane and above by a deformable upper surface, subjected to a vertical temperature gradient. Since the viscosity is temperature-dependent the consequences of relaxing Oberbeck-Boussinesq approximation and free surface deformability are theoretically examined by means of small disturbance analysis. The problem has been solved numerically by the Taylor series expansion method. The results obtained confirm that when the free surface is undeformable, stationary convection develops in the form of polygonal cells, and oscillatory motion cannot be detected. When the surface deformability is considered, stationary convection sets in, either as a short-wavelength hexagonal instability or as a long-wavelengh mode or as both, and oscillatory convection is also possible. The stability threshold for the short-wavelength mode depends mainly on the viscosity variation while the long-wavelength mode is determined by the surface deformation. Numerically, it is found that the neutral oscillatory Marangoni numbers are only negative. When a variable-viscosity model is used the theoretical and experimental results are in better agreement. Received 15 May 1997     

Effect of parametric modulation on pattern formation in reaction diffusion systems

Reaction diffusion systems can exhibite both spatial and temporal patterns. We show that the effect of spatial variation of the removal rate can have significant effect on the stability boundaries. In particular there can be a case of parametric resonance. Received 11 March 1998     

Maximum flow problem on dynamic generative network flows with time-varying bounds

This paper considers a new class of network flows, called dynamic generative network flows in which, the flow commodity is dynamically generated at a source node and dynamically consumed at a sink node and the arc-flow bounds are time dependent. Then the maximum dynamic flow problem in such networks for a pre-specified time horizon T is defined and mathematically formulated in both arc flow and path flow presentations. By exploiting the special structure of the problem, an efficient algorithm is developed to solve the general form of the dynamic problem as a minimum cost static flow problem.     

A matrix-free high-order discontinuous Galerkin compressible Navier-Stokes solver: A performance comparison of compressible and incompressible formulations for turbulent incompressible flows

Both compressible and incompressible Navier-Stokes solvers can be used and are used to solve incompressible turbulent flow problems. In the compressible case, the Mach number is then considered as a solver parameter that is set to a small value, M ≈0.1, in order to mimic incompressible flows. This strategy is widely used for high-order discontinuous Galerkin (DG) discretizations of the compressible Navier-Stokes equations. The present work raises the question regarding the computational efficiency of compressible DG solvers as compared to an incompressible formulation. Our contributions to the state of the art are twofold: Firstly, we present a high-performance DG solver for the compressible Navier-Stokes equations based on a highly efficient matrix-free implementation that targets modern cache-based multicore architectures with Flop/Byte ratios significantly larger than 1. The performance results presented in this work focus on the node-level performance, and our results suggest that there is great potential for further performance improvements for current state-of-the-art DG implementations of the compressible Navier-Stokes equations. Secondly, this compressible Navier-Stokes solver is put into perspective by comparing it to an incompressible DG solver that uses the same matrix-free implementation. We discuss algorithmic differences between both solution strategies and present an in-depth numerical investigation of the performance. The considered benchmark test cases are the three-dimensional Taylor-Green vortex problem as a representative of transitional flows and the turbulent channel flow problem as a representative of wall-bounded turbulent flows. The results indicate a clear performance advantage of the incompressible formulation over the compressible one.     

Beyond direct simulation Monte Carlo (DSMC) modelling of collision environments

ABSTRACT

Direct simulation Monte Carlo (DSMC) models have been successfully adopted and adapted to describe gas flows in a wide range of environments since the method was first introduced by Bird in the 1960s. We propose a new approach to modelling collisions between gas-phase particles in this work – operating in a similar way to the DSMC model, but with one key difference. Particles move in a mean field, generated by all previously propagated particles, which removes the requirement that all particles be propagated simultaneously. This yields a significant reduction in computation effort and lends itself to applications for which DSMC becomes intractable, such as when a species of interest is only a minor component of a large gas mixture.     

Simulating bubble dynamics in a buoyant system

Multiphase flows are critical components of many physical systems; however, numerical models of multiphase flows with large parameter gradients can be challenging. Here, two different numerical methods, volume of fluid (VOF) and smoothed particle hydrodynamics (SPH), are used to model the buoyant rise of isolated gas bubbles through quiescent fluids for a range of Bond and Reynolds numbers. The VOF is an Eulerian grid–based method, whereas the SPH is Lagrangian and mesh free. Each method has unique strengths and weaknesses, and a comparison of the two approaches as applied to multiphase phenomena has not previously been performed. The VOF and SPH simulations are compared, verified, and validated. Results using two-dimensional VOF and SPH simulations are similar to each other and are able to reproduce numerical benchmarks and experimental results for sufficiently large Morton and Reynolds numbers. It is also shown that at low Reynolds numbers, the two methods, SPH and VOF, diverge in the transient regime of the bubble rise. Regimes that require simulations capable of representing three-dimensional drag are identified as well as regimes in which results from VOF and SPH diverge.     

Precise measurement and control of the pressure-driven flows for microfluidic systems

The pressure-driven device is designed and the flow rates of the microfluidic systems can be supplied by the pressure-driven flows, which can significantly reduce the flow-rate fluctuations coming from the pump source. For pressure-driven flows, the flow rates of the fluids can be predicted by measuring the pressure drop along a polytetrafluoroethylene (PTFE) tubing. Especially, by varying the geometrical parameters of the PTFE tubing, the predicted flow rates of the fluids are compared with the experimental measurements, and the testing precision of the pressure-driven flows can be obtained. Meanwhile, the dynamic characteristics of the open-loop and closed-loop control pressure-driven device are comparatively studied. Particularly, a proportional and integral (PI) controller is integrated with the closed-loop control pressure-driven device, and the effects of the parameters of the PI controller on the dynamic characteristics of the pressure-driven devices are mainly discussed. Most importantly, by improving the dynamic characteristics of the pressure-driven devices, precise measurement and control of the pressure-driven flows can be achieved for microfluidic systems.