Numerical simulation of a bubble rising in shear-thinning fluids

The motion of a single bubble rising freely in quiescent non-Newtonian viscous fluids was investigated experimentally and computationally. The non-Newtonian effects in the flow of viscous inelastic fluids are modeled by the Carreau theological model. An improved level set approach for computing the incompressible two-phase flow with deformable free interface is used. The control volume formulation with the SIMPLEC algorithm incorporated is used to solve the governing equations on a staggered Eulerian grid. The simulation results demonstrate that the algorithm is robust for shear-thinning liquids with large density (rho(1)/rho(g) up to 10(3)) and high viscosity (eta(1)/eta(g) up to 10(4)). The comparison of the experimental measurements of terminal bubble shape and velocity with the computational results is satisfactory. It is shown that the local change in viscosity around a bubble greatly depends on the bubble shape and the zero-shear viscosity of non-Newtonian shear-thinning liquids. The shear-rate distribution and velocity fields are used to elucidate the formation of a region of large viscosity at the rear of a bubble as a result of the rather stagnant flow behind the bubble. The numerical results provide the basis for further investigations, such as the numerical simulation of viscoelastic fluids. (C) 2010 Elsevier B.V. All rights reserved.

Simulation of the Lanzhou Penning Trap LPT

The LPT (Lanzhou Penning Trap) is under construction and its task is to perforin direct mass measurement of fusion-evaporation residue, and if possible for heavy isotopes. Detailed simulations have been clone for a good understanding to the ion's movement and mechanics in the trap. The optimizal ion of the LPT is also performed based on the simulation. With a scale of 0.5 mm per grid used in the, simulation and many other limitations a highest mass resolution has been achieved to be 1.9 x 10(-5). An unexpected behavioin in the simulation related to magnetron motion has been found.

Numerical simulation of coil-helix transition processes of gelatin

Gelatin is widely used in food, pharmaceutical, and photographic industries due to the coil-helix transition, whereas the structural inhomogeneity considerably affects its essential properties closely connecting with the industrial applications. The spatially structural inhomogeneity of the gelatin caused by the uneven and unstable temperature field is analyzed by the finite element method during the cooling-induced coil-helix transition process. The helix conversion and the crosslinking density as functions of time and spatial grid are calculated by the incremental method. A length distribution density function is introduced to describe the continuous length distributions of two kinds of triple helices.

Numerical simulation of tides and tidal currents in Liaodong Bay with POM

A three-dimensional ocean circulation model, called Princeton Ocean Model (POM), is employed to simulate tides and tidal currents in Liaodong Bay. The nested grid technique is adopted to improve the computation precision. Computed harmonic constants of M-1, M-2 tides at five tidal gauge stations and surface elevations at two oil platforms are compared with those observed, and show good agreements with them. Based on the calculated results, the co-amplitude and co-phase tag lines of nil and M-2 tidal constituents, the residual current field of M-2 constituent, tidal form, tidal Current ellipse and the moving style of tidal current are given. It is found that diurnal tidal constituents have no amphidromic point whereas semi-diurnal constituents have one in the region of interest. Meanwhile, some meaningful results are concluded and presented, which are conducive to a thorough knowledge of the characteristics of tides and tidal currents in the Liaodong Bay.

Monotone scheme and boundary conditions for finite volume simulation of magnetohydrodynamic internal flows at high Hartmann number

A monotone scheme for finite volume simulation of magnetohydrodynamic internal flows at high Hartmann number is presented. The numerical stability is analysed with respect to the electromagnetic force. Standard central finite differences applied to finite volumes can only be numerically stable if the vector products involved in this force are computed with a scheme using a fully staggered grid. The electromagnetic quantities (electric currents and electric potential) must be shifted by half the grid size from the mechanical ones (velocity and pressure). An integral treatment of the boundary layers is used in conjunction with boundary conditions for electrically conducting walls. The simulations are performed with inhomogeneous electrical conductivities of the walls and reach high Hartmann numbers in three-dimensional simulations, even though a non-adaptive grid is used.

Feature-oriented grid topology design for aerospace configurations

The last few years have seen a substantial increase in the geometric complexity for 3D flow simulation. In this paper we describe the challenges in generating computation grids for 3D aerospace configuations and demonstrate the progress made to eventually achieve a push button technology for CAD to visualized flow. Special emphasis is given to the interfacing from the grid generator to the flow solver by semi-automatic generation of boundary conditions during the grid generation process. In this regard, once a grid has been generated, push button technology of most commercial flow solvers has been achieved. This will be demonstrated by the ad hoc simulation for the Hopper configuration.

Numerical simulation of flow-induced cavity noise in self-sustained oscillations

The generation and near-field radiation of aerodynamic sound from a low-speed unsteady flow over a two-dimensional automobile door cavity is simulated by using a source-extraction-based coupling method. In the coupling procedure, the unsteady cavity flow field is first computed solving the Reynolds averaged Navier–Stokes (RANS) equations. The radiated sound is then calculated by using a set of acoustic perturbation equations with acoustic source terms which are extracted from the time-dependent solutions of the unsteady flow. The aerodynamic and its resulting acoustic field are computed for the Reynolds number of 53,266 based on the base length of the cavity. The free stream flow velocity is taken to be 50.9m/s. As first stage of the numerical investigation of flow-induced cavity noise, laminar flow is assumed. The CFD solver is based on a cell-centered finite volume method. A dispersion-relation-preserving (DRP), optimized, fourth-order finite difference scheme with fully staggered-grid implementation is used in the acoustic solver

Coordinated control of DFIG’s rotor and grid side converters during network unbalance

This paper proposes a coordinated control of the rotor and grid side converters (RSC & GSC) of doubly-fed induction generator (DFIG) based wind generation systems under unbalanced voltage conditions. System behaviors and operations of the RSC and GSC under unbalanced voltage are illustrated. To provide enhanced operation, the RSC is controlled to eliminate the torque oscillations at double supply frequency under unbalanced stator supply. The oscillation of the stator output active power is then cancelled by the active power output from the GSC, to ensure constant active power output from the overall DFIG generation system. To provide the required positive and negative sequence currents control for the RSC and GSC, a current control strategy containing a main controller and an auxiliary controller is analyzed. The main controller is implemented in the positive (dq)+ frame without involving positive/negative sequence decomposition whereas the auxiliary controller is implemented in the negative sequence (dq)? frame with negative sequence current extracted. Simulation results using EMTDC/PSCAD are presented for a 2MW DFIG wind generation system to validate the proposed control scheme and to show the enhanced system operation during unbalanced voltage supply.

Improved Direct Power Control of Grid Connected DC/AC Converters

This paper proposes new direct power control (DPC) strategies for three-phase DC/AC converters with improved dynamic response and steady-state performance. As with an electrical machine, source and converter flux which equal the integration of the respective source and converter voltage are used to define active and reactive power flow. Optimization of the look-up-table used in conventional DPC is outlined first, to improve the power control and reduce the current distortion. Then constant switching frequency DPC is developed where the required converter voltage vector within a fixed half switching period is calculated directly from the active and reactive power errors. Detailed angle compensation due to the finite sampling frequency and the use of integral controller to further improve the power control accuracy, are described. Both simulation and experimental results are used to compare conventional DPC and vector control, and to demonstrate the effectiveness and robustness of the proposed control strategies during active and reactive power steps, and line inductance variations.

Coordinated Control of DFIG and FSIG-based Wind Farms under Unbalanced Grid Conditions

This paper investigates the control and operation of doubly-fed induction generator (DFIG) and fixed-speed induction generator (FSIG) based wind farms under unbalanced grid conditions. A DFIG system model suitable for analyzing unbalanced operation is developed, and used to assess the impact of an unbalanced supply on DFIG and FSIG operation. Unbalanced voltage at DFIG and FSIG terminals can cause unequal heating on the stator windings, extra mechanical stresses and output power fluctuations. These problems are particularly serious for the FSIG-based wind farm without a power electronic interface to the grid. To improve the stability of a wind energy system containing both DFIG and FSIG based wind farms during network unbalance, a control strategy of unbalanced voltage compensation by the DFIG systems is proposed. The DFIG system compensation ability and the impact of transmission network impedance are illustrated. The simulation results implemented in Matlab/Simulink show that the proposed DFIG control system improves not only its own performance, but also the stability of the FSIG system with the same grid connection point during network unbalance.

Wideband and Isotropic Room Acoustics Simulation Using 2-D Interpolated FDTD Schemes

In this paper, a complete method for finite-difference time-domain modeling of rooms in 2-D using compact explicit schemes is presented. A family of interpolated schemes using a rectilinear, nonstaggered grid is reviewed, and the most accurate and isotropic schemes are identified. Frequency-dependent boundaries are modeled using a digital impedance filter formulation that is consistent with locally reacting surface theory. A structurally stable and efficient boundary formulation is constructed by carefully combining the boundary condition with the interpolated scheme. An analytic prediction formula for the effective numerical reflectance is given, and a stability proof provided. The results indicate that the identified accurate and isotropic schemes are also very accurate in terms of numerical boundary reflectance, and outperform directly related methods such as Yee's scheme and the standard digital waveguide mesh. In addition, one particular scheme-referred to here as the interpolated wideband scheme-is suggested as the best scheme for most applications.

Room acoustics simulation using 3-D compact explicit FDTD schemes

This paper presents methods for simulating room acoustics using the finite-difference time-domain (FDTD) technique, focusing on boundary and medium modeling. A family of nonstaggered 3-D compact explicit FDTD schemes is analyzed in terms of stability, accuracy, and computational efficiency, and the most accurate and isotropic schemes based on a rectilinear grid are identified. A frequency-dependent boundary model that is consistent with locally reacting surface theory is also presented, in which the wall impedance is represented with a digital filter. For boundaries, accuracy in numerical reflection is analyzed and a stability proof is provided. The results indicate that the proposed 3-D interpolated wideband and isotropic schemes outperform directly related techniques based on Yee's staggered grid and standard digital waveguide mesh, and that the boundary formulations generally have properties that are similar to that of the basic scheme used.

A Physically-Constrained Source Model For FDTD Acoustic Simulation

The Finite Difference Time Domain (FDTD) method is becoming increasingly popular for room acoustics simulation. Yet, the literature on grid excitation methods is relatively sparse, and source functions are traditionally implemented in a hard or additive form
using arbitrarily-shaped functions which do not necessarily obey the physical laws of sound generation. In this paper we formulate
a source function based on a small pulsating sphere model. A physically plausible method to inject a source signal into the grid
is derived from first principles, resulting in a source with a near-flat spectrum that does not scatter incoming waves. In the final
discrete-time formulation, the source signal is the result of passing a Gaussian pulse through a digital filter simulating the dynamics of the pulsating sphere, hence facilitating a physically correct means to design source functions that generate a prescribed sound field.

Impact of Cyber-Security Issues on Smart Grid

Greater complexity and interconnectivity across systems embracing Smart Grid technologies has meant that cyber-security issues have attracted significant attention. This paper describes pertinent cyber-security requirements, in particular cyber attacks and countermeasures which are critical for reliable Smart Grid operation. Relevant published literature is presented for critical aspects of Smart Grid cyber-security, such as vulnerability, interdependency, simulation, and standards. Furthermore, a preliminary study case is given which demonstrates the impact of a cyber attack which violates the integrity of data on the load management of real power system. Finally, the paper proposes future work plan which focuses on applying intrusion detection and prevention technology to address cyber-security issues. This paper also provides an overview of Smart Grid cyber-security with reference to related cross-disciplinary research topics.