Shock-induced two-phase flows in an aligned baffle system filled with suspended particles

This paper describes the shock propagation through a dilute gas-particle suspension in an aligned baffle system. Numerical solution to two-phase flows induced by a planar shock wave is given based on the two-continuum model with interphase coupling. The governing equations are numerically solved by using high-resolution schemes. The computational results show the shock reflection and diffraction patterns, and the shock-induced flow fields in the 4-baffle system filled with the dusty gas.

Self-consistent kinetic modeling of low-pressure inductively coupled radio-frequency discharges

An efficient method for solving the spatially inhomogeneous Boltzmann equation in a two-term approximation for low-pressure inductively coupled plasmas has been developed. The electron distribution function (EDF), a function of total electron energy and two spatial coordinates, is found self-consistently with the static space-charge potential which is computed from a 2D fluid model, and the rf electric field profile which is calculated from the Maxwell equations. The EDF and the spatial distributions of the electron density, potential, temperature, ionization rate, and the inductive electric field are calculated and discussed. (C) 1996 American Institute of Physics.

Self-consistent kinetic modeling of electron distribution function in a low-pressure inductively coupled plasma

Based upon the spatially inhomogeneous Boltzmann equation in two-term approximation coupled with electromagnetic and fluid model analysis for the recently developed inductively coupled plasma sources, a self-consistent electron kinetic model is developed. The electron distribution function, spatial distributions of the electron density and ionization rate are calculated and discussed.

利用Ghost Fluid方法模拟激波与柱形界面相互作用

利用Ghost Fluid方法(后面简称Ghost方法)和γ-model方法，在同样的时空离散精度条件下，对激波与柱形界面相互作用的二维可压缩流场进行了直接模拟，并与实验结果相比较．从模拟结果看，在短时间内， Ghost方法和γ-model方法模拟的结果与实验结果基本相同，两种方法均正确地模拟出界面的位置、激波的强度和速度．但随着时间的发展，具有较大数值耗散的γ-model方法的计算结果与实验差别越来越大；而数值耗散较小的Ghost方法能较为正确地模拟界面的运动．

Instabilities of a liquid film flowing down an inclined porous plane

The problem of a film flowing down an inclined porous layer is considered. The fully developed basic flow is driven by gravitation. A careful linear instability analysis is carried out. We use Darcy's law to describe the porous layer and solve the coupling equations of the fluid and the porous medium rather than the decoupled equations of the one-sided model used in previous works. The eigenvalue problem is solved by means of a Chebyshev collocation method. We compare the instability of the two-sided model with the results of the one-sided model. The result reveals a porous mode instability which is completely neglected in previous works. For a falling film on an inclined porous plane there are three instability modes, i.e., the surface mode, the shear mode, and the porous mode. We also study the influences of the depth ratio d, the Darcy number delta, and the Beavers-Joseph coefficient alpha(BJ) on the instability of the system.

Short-pulse laser absorption via J x B heating in ultrahigh intensity laser plasma interaction

An analytical fluid model for JxB heating during the normal incidence by a short ultraintense linearly polarized laser on a solid-density plasma is proposed. The steepening of an originally smooth electron density profile as the electrons are pushed inward by the laser is included self-consistently. It is shown that the JxB heating includes two distinct coupling processes depending on the initial laser and plasma conditions: for a moderate intensity (a <= 1), the ponderomotive force of the laser light can drive a large plasma wave at the point n(e)=4 gamma(0)n(c) resonantly. When this plasma wave is damped, the energy is transferred to the plasma. At higher intensity, the electron density is steepened to a high level by the time-independent ponderomotive force, n(e)> 4 gamma(0)n(c), so that no 2 omega resonance will occur, but the longitudinal component of the oscillating ponderomotive field can lead to an absorption mechanism similar to "vacuum heating." (c) 2006 American Institute of Physics.

Standing ultrashort relativistic solitons in overdense plasmas

By using a one-dimensional self-consistent relativistic fluid model, an investigation is made numerically on relativistic electromagnetic solitons with a high intensity in cold overdense plasmas with an electrons' initial velocity opposite to the laser propagating direction. Two types of standing solitons with zero group velocity are found at the given electrons' initial velocities. One is single-humped with a weakly relativistic intensity; the another is multi-humped with a strong relativistic amplitude. The properties of these two types of solitons are presented in detail.

Battery-Powered RF Pre-Ionization System for the Caltech Magnetohydrodynamically-Driven Jet Experiment: RF Discharge Properties and MHD-Driven Jet Dynamics

This thesis describes investigations of two classes of laboratory plasmas with rather different properties: partially ionized low pressure radiofrequency (RF) discharges, and fully ionized high density magnetohydrodynamically (MHD)-driven jets. An RF pre-ionization system was developed to enable neutral gas breakdown at lower pressures and create hotter, faster jets in the Caltech MHD-Driven Jet Experiment. The RF plasma source used a custom pulsed 3 kW 13.56 MHz RF power amplifier that was powered by AA batteries, allowing it to safely float at 4-6 kV with the cathode of the jet experiment. The argon RF discharge equilibrium and transport properties were analyzed, and novel jet dynamics were observed.

Although the RF plasma source was conceived as a wave-heated helicon source, scaling measurements and numerical modeling showed that inductive coupling was the dominant energy input mechanism. A one-dimensional time-dependent fluid model was developed to quantitatively explain the expansion of the pre-ionized plasma into the jet experiment chamber. The plasma transitioned from an ionizing phase with depressed neutral emission to a recombining phase with enhanced emission during the course of the experiment, causing fast camera images to be a poor indicator of the density distribution. Under certain conditions, the total visible and infrared brightness and the downstream ion density both increased after the RF power was turned off. The time-dependent emission patterns were used for an indirect measurement of the neutral gas pressure.

The low-mass jets formed with the aid of the pre-ionization system were extremely narrow and collimated near the electrodes, with peak density exceeding that of jets created without pre-ionization. The initial neutral gas distribution prior to plasma breakdown was found to be critical in determining the ultimate jet structure. The visible radius of the dense central jet column was several times narrower than the axial current channel radius, suggesting that the outer portion of the jet must have been force free, with the current parallel to the magnetic field. The studies of non-equilibrium flows and plasma self-organization being carried out at Caltech are relevant to astrophysical jets and fusion energy research.

A molecular dynamics study of the microscopic properties of simple dense fluids

The microscopic properties of a two-dimensional model dense fluid of Lennard-Jones disks have been studied using the so-called "molecular dynamics" method. Analyses of the computer-generated simulation data in terms of "conventional" thermodynamic and distribution functions verify the physical validity of the model and the simulation technique.

The radial distribution functions g(r) computed from the simulation data exhibit several subsidiary features rather similar to those appearing in some of the g(r) functions obtained by X-ray and thermal neutron diffraction measurements on real simple liquids. In the case of the model fluid, these "anomalous" features are thought to reflect the existence of two or more alternative configurations for local ordering.

Graphical display techniques have been used extensively to provide some intuitive insight into the various microscopic phenomena occurring in the model. For example, "snapshots" of the instantaneous system configurations for different times show that the "excess" area allotted to the fluid is collected into relatively large, irregular, and surprisingly persistent "holes". Plots of the particle trajectories over intervals of 2.0 to 6.0 x 10^-12 sec indicate that the mechanism for diffusion in the dense model fluid is "cooperative" in nature, and that extensive diffusive migration is generally restricted to groups of particles in the vicinity of a hole.

A quantitative analysis of diffusion in the model fluid shows that the cooperative mechanism is not inconsistent with the statistical predictions of existing theories of singlet, or self-diffusion in liquids. The relative diffusion of proximate particles is, however, found to be retarded by short-range dynamic correlations associated with the cooperative mechanism--a result of some importance from the standpoint of bimolecular reaction kinetics in solution.

A new, semi-empirical treatment for relative diffusion in liquids is developed, and is shown to reproduce the relative diffusion phenomena observed in the model fluid quite accurately. When incorporated into the standard Smoluchowski theory of diffusion-controlled reaction kinetics, the more exact treatment of relative diffusion is found to lower the predicted rate of reaction appreciably.

Finally, an entirely new approach to an understanding of the liquid state is suggested. Our experience in dealing with the simulation data--and especially, graphical displays of the simulation data--has led us to conclude that many of the more frustrating scientific problems involving the liquid state would be simplified considerably, were it possible to describe the microscopic structures characteristic of liquids in a concise and precise manner. To this end, we propose that the development of a formal language of partially-ordered structures be investigated.

Distributed Control and Optimization for Communication and Power Systems

We are at the cusp of a historic transformation of both communication system and electricity system. This creates challenges as well as opportunities for the study of networked systems. Problems of these systems typically involve a huge number of end points that require intelligent coordination in a distributed manner. In this thesis, we develop models, theories, and scalable distributed optimization and control algorithms to overcome these challenges.

This thesis focuses on two specific areas: multi-path TCP (Transmission Control Protocol) and electricity distribution system operation and control. Multi-path TCP (MP-TCP) is a TCP extension that allows a single data stream to be split across multiple paths. MP-TCP has the potential to greatly improve reliability as well as efficiency of communication devices. We propose a fluid model for a large class of MP-TCP algorithms and identify design criteria that guarantee the existence, uniqueness, and stability of system equilibrium. We clarify how algorithm parameters impact TCP-friendliness, responsiveness, and window oscillation and demonstrate an inevitable tradeoff among these properties. We discuss the implications of these properties on the behavior of existing algorithms and motivate a new algorithm Balia (balanced linked adaptation) which generalizes existing algorithms and strikes a good balance among TCP-friendliness, responsiveness, and window oscillation. We have implemented Balia in the Linux kernel. We use our prototype to compare the new proposed algorithm Balia with existing MP-TCP algorithms.

Our second focus is on designing computationally efficient algorithms for electricity distribution system operation and control. First, we develop efficient algorithms for feeder reconfiguration in distribution networks. The feeder reconfiguration problem chooses the on/off status of the switches in a distribution network in order to minimize a certain cost such as power loss. It is a mixed integer nonlinear program and hence hard to solve. We propose a heuristic algorithm that is based on the recently developed convex relaxation of the optimal power flow problem. The algorithm is efficient and can successfully computes an optimal configuration on all networks that we have tested. Moreover we prove that the algorithm solves the feeder reconfiguration problem optimally under certain conditions. We also propose a more efficient algorithm and it incurs a loss in optimality of less than 3% on the test networks.

Second, we develop efficient distributed algorithms that solve the optimal power flow (OPF) problem on distribution networks. The OPF problem determines a network operating point that minimizes a certain objective such as generation cost or power loss. Traditionally OPF is solved in a centralized manner. With increasing penetration of volatile renewable energy resources in distribution systems, we need faster and distributed solutions for real-time feedback control. This is difficult because power flow equations are nonlinear and kirchhoff's law is global. We propose solutions for both balanced and unbalanced radial distribution networks. They exploit recent results that suggest solving for a globally optimal solution of OPF over a radial network through a second-order cone program (SOCP) or semi-definite program (SDP) relaxation. Our distributed algorithms are based on the alternating direction method of multiplier (ADMM), but unlike standard ADMM-based distributed OPF algorithms that require solving optimization subproblems using iterative methods, the proposed solutions exploit the problem structure that greatly reduce the computation time. Specifically, for balanced networks, our decomposition allows us to derive closed form solutions for these subproblems and it speeds up the convergence by 1000x times in simulations. For unbalanced networks, the subproblems reduce to either closed form solutions or eigenvalue problems whose size remains constant as the network scales up and computation time is reduced by 100x compared with iterative methods.

Gas properties as a limit to gas turbine performance

Although increasing the turbine inlet temperature has traditionally proved the surest way to increase cycle efficiency, recent work suggests that the performance of future gas turbines may be limited by increased cooling flows and losses. Another limiting scenario concerns the effect on cycle performance of real gas properties at high temperatures. Cycle calculations of uncooled gas turbines show that when gas properties are modelled accurately, the variation of cycle efficiency with turbine inlet temperature at constant pressure ratio exhibits a maximum at temperatures well below the stoichiometric limit. Furthermore, the temperature at the maximum decreases with increasing compressor and turbine polytropic efficiency. This behaviour is examined in the context of a two-component model of the working fluid. The dominant influences come from the change of composition of the combustion products with varying air/fuel ratio (particularly the contribution from the water vapour) together with the temperature variation of the specific heat capacity of air. There are implications for future industrial development programmes, particularly in the context of advanced mixed gas-steam cycles.

Deposition of small particles from turbulent flows

This paper deals with particle deposition onto solid walls from turbulent flows. The aim of the study is to model particle deposition in industrial flows, such as the one in gas turbines. The numerical study has been carried out with a two fluid approach. The possible contribution to the deposition from Brownian diffusion, turbulent diffusion and shear-induced lift force are considered in the study. Three types of turbulent two-phase flows have been studied: turbulent channel flow, turbulent flow in a bent duct and turbulent flow in a turbine blade cascade. In the turbulent channel flow case, the numerical results from a two-dimensional code show good agreement with numerical and experimental results from other resources. Deposition problem in a bent duct flow is introduced to study the effect of curvature. Finally, the deposition of small particles on a cascade of turbine blades is simulated. The results show that the current two fluid models are capable of predicting particle deposition rates in complex industrial flows.

Computational fluid-structure interaction of DGB parachutes in compressible fluid flow

This paper describes large-scale simulations of compressible flows over a supersonic disk-gap-band parachute system. An adaptive mesh refinement method is used to resolve the coupled fluid-structure model. The fluid model employs large-eddy simulation to describe the turbulent wakes appearing upstream and downstream of the parachute canopy and the structural model employed a thin-shell finite element solver that allows large canopy deformations by using subdivision finite elements. The fluid-structure interaction is described by a variant of the Ghost-Fluid method. The simulation was carried out at Mach number 1.96 where strong nonlinear coupling between the system of bow shocks, turbulent wake and canopy is observed. It was found that the canopy oscillations were characterized by a breathing type motion due to the strong interaction of the turbulent wake and bow shock upstream of the flexible canopy. Copyright © 2010 by ASME.

Magnetically-actuated artificial cilia for microfluidic propulsion.

In this paper we quantitatively analyse the performance of magnetically-driven artificial cilia for lab-on-a-chip applications. The artificial cilia are fabricated using thin polymer films with embedded magnetic nano-particles and their deformation is studied under different external magnetic fields and flows. A coupled magneto-mechanical solid-fluid model that accurately captures the interaction between the magnetic field, cilia and fluid is used to simulate the cilia motion. The elastic and magnetic properties of the cilia are obtained by fitting the results of the computational model to the experimental data. The performance of the artificial cilia with a non-uniform cross-section is characterised using the numerical model for two channel configurations that are of practical importance: an open-loop and a closed-loop channel. We predict that the flow and pressure head generated by the artificial cilia can be as high as 18 microlitres per minute and 3 mm of water, respectively. We also study the effect of metachronal waves on the flow generated and show that the fluid propelled increases drastically compared to synchronously beating cilia, and is unidirectional. This increase is significant even when the phase difference between adjacent cilia is small. The obtained results provide guidelines for the optimal design of magnetically-driven artificial cilia for microfluidic propulsion.

Case study on a deep excavation in soft ground by centrifuge model tests

This paper provides a case study on the deepest excavation carried out so far in the construction of the metro network in Shanghai, which typically features soft ground. The excavation is 38 m deep with retaining walls 65 m deep braced by 9 levels of concrete props. To obtain a quick and rough prediction, two centrifuge model tests were conducted, in which one is for the 'standard' section with green field surrounding and the other with an adjacent piled building. The tests were carried out in a run-stop-excavation-run style, in which excavation was conducted manually. By analyzing the lateral wall displacement, ground deformation, bending moment and earth pressure, the test results are shown to be reasonably convincing and the design and construction were validated. Such industry orientated centrifuge modeling was shown to be useful in understanding the performance of geotechnical processes, especially when engineers lack relevant field experience. © 2010 Taylor & Francis Group, London.