二维抛物化稳定性方程的特征和次特征

特征分析表明：对原始扰动量的抛物化稳定性方程组（PSE），它在亚、超音速区分别具有椭圆和抛物特性，给出PSE特征对马赫数的依赖关系，阐明PSE仅把信息对流-扩散传播特性抛物化，而保留了信息对流-扰动传播特性，因此PSE应称为扩散抛物化稳定性方程（DPSE）。

A Multiscale-Domain Decomposition Method for High Reynolds Number Flows

The high Reynolds number flow contains a wide range of length and time scales, and the flow domain can be divided into several sub-domains with different characteristic scales. In some sub-domains, the viscosity dissipation scale can only be considered in a certain direction; in some sub-domains, the viscosity dissipation scales need to be considered in all directions; in some sub-domains, the viscosity dissipation scales are unnecessary to be considered at all. For laminar boundary layer region, the characteristic length scales in the streamwise and normal directions are L and L Re-1/ 2 , respectively. The characteristic length scale and the velocity scale in the outer region of the boundary layer are L and U, respectively. In the neighborhood region of the separated point, the length scale l<the velocity scale u<the light of this reason there exists a great disparity of the grid Reynolds numbers RDxi between different cells in Navier-Stokes (NS) equations computations for high Reynolds number flows, an idea of solving the conservation equations for discrete cells was proposed and named the discrete fluid dynamics (DFD) algorithm. Analysis shows that the basic conservative equations for discrete cells are the Euler equations, NS- and diffusion parabolized (DP) NS equations. In this paper, a new multiscale-domain decomposition method is developed for the high Reynolds number flow. First, the whole domain is decomposed to different sub-domains with the different characteristic scales. Then the different dominant equation of all sub-domains is defined according to the diffusion parabolized (DP) theory of viscous flow. Finally these different equations are solved simultaneously in whole computational region. For numerical tests of high Reynolds numerical flows, two-dimensional supersonic flows over rearward and frontward steps as well as an interaction flow between shock wave and boundary layer were solved numerically. The pressure distributions and local coefficients of skin friction on the wall are given. The numerical results obtained by the multiscale-domain decomposition algorithm are well agreement with those by NS equations. Comparing with the usual method of solving the Navier-Stokes equations in the whole flow, under the same numerical accuracy, the present multiscale domain decomposition method decreases CPU consuming about 20% and reflects the physical mechanism of practical flow more accurately.

Effects of Thermocapillary Convection on the Growth Process in Industrial 8-inch Silicon Czochralski Growth

Czochralski (CZ) crystal growth process is a widely used technique in manufacturing of silicon crystals and other semiconductor materials. The ultimate goal of the IC industry is to have the highest quality substrates, which are free of point defect, impurities and micro defect clusters. The scale up of silicon wafer size from 200 mm to 300 mm requires large crucible size and more heat power. Transport phenomena in crystal growth processes are quite complex due to melt and gas flows that may be oscillatory and/or turbulent, coupled convection and radiation, impurities and dopant distributions, unsteady kinetics of the growth process, melt crystal interface dynamics, free surface and meniscus, stoichiometry in the case of compound materials. A global model has been developed to simulate the temperature distribution and melt flow in an 8-inch system. The present program features the fluid convection, magnetohydrodynamics, and radiation models. A multi-zone method is used to divide the Cz system into different zones, e.g., the melt, the crystal and the hot zone. For calculation of temperature distribution, the whole system inside the stainless chamber is considered. For the convective flow, only the melt is considered. The widely used zonal method divides the surface of the radiation enclosure into a number of zones, which has a uniform distribution of temperature, radiative properties and composition. The integro-differential equations for the radiative heat transfer are solved using the matrix inversion technique. The zonal method for radiative heat transfer is used in the growth chamber, which is confined by crystal surface, melt surface, heat shield, and pull chamber. Free surface and crystal/melt interface are tracked using adaptive grid generation. The competition between the thermocapillary convection induced by non-uniform temperature distributions on the free surface and the forced convection by the rotation of the crystal determines the interface shape, dopant distribution, and striation pattern. The temperature gradients on the free surface are influenced by the effects of the thermocapillary force on the free surface and the rotation of the crystal and the crucible.

Nonlinear Dynamic Responses of the Tensioned Tether under Parametric Excitations

The nonlinear dynamic responses of the tensioned tether subjected to combined surge and heave motions of floating platform are investigated using 2-D nonlinear beam model. It is shown that if the transverse-axial coupling of nonlinear beam model and the combined surge-heave motions of platform are considered, the governing equation is not Mathieu equation any more, it becomes nonlinear Hill equation. The Hill stability chart is obtained by using the Hill's infinite determinant and harmonic balance method. A parameter M, which is the function of tether length, the surge and heave amplitude of platform, is defined. The Hill stability chart is obviously different from Mathieu stability chart which is the specific case as M=0. Some case studies are performed by employing linear and nonlinear beam model respectively. It can be found that the results differences between nonlinear and linear model are apparent.

Basic equations for ferroelectrics

A set of new formula of energy functions for ferroelectrics was proposed, and then the new basic equations were derived in this paper. The finite element formulation based on the new basic equations was improved to avoid the equivalent nodal load produced by remnant polarization. With regard to the fundamentals of mathematics and physics, the new energy functions and basic equations are reasonable for the material element of ferroelectrics in finite element analysis.

Seasonal Stability Tests in gretl. An Application to International Tourism Data

The seasonal stability tests of Canova & Hansen (1995) (CH) provide a method complementary to that of Hylleberg et al. (1990) for testing for seasonal unit roots. But the distribution of the CH tests are unknown in small samples. We present a method to numerically compute critical values and P-values for the CH tests for any sample size and any seasonal periodicity. In fact this method is applicable to the types of seasonality which are commonly in use, but also to any other.

Seasonal, horizontal, and vertical distribution of phytoplankton chlorophyll a in the Northeast U.S. Continental Shelf Ecosystem

The broad scale features in the horizontal, vertical, and seasonal distribution of phytoplankton chlorophyll a on the northeast U.S. continental shelf are described based on 57,088 measurements made during 78 oceanographic surveys from 1977 through 1988. Highest mean water column chlorophyll concentration (Chlw,) is usually observed in nearshore areas adjacent to the mouths of the estuaries in the Middle Atlantic Bight (MAB), over the shallow water on Georges Bank, and a small area sampled along the southeast edge of Nantucket Shoals. Lowest Chlw «0.125 ug l-1) is usually restricted to the most seaward stations sampled along the shelf-break and the central deep waters in the Gulf of Maine. There is at least a twofold seasonal variation in phytoplankton biomass in all areas, with highest phytoplankton concentrations (m3) and highest integrated standing stocks (m2) occurring during the winter-spring (WS) bloom, and the lowest during summer, when vertical density stratification is maximal. In most regions, a secondary phytoplankton biomass pulse is evident during convective destratification in fall, usually in October. Fall bloom in some areas of Georges Bank approaches the magnitude of the WS-bloom, but Georges Bank and Middle Atlantic Bight fall blooms are clearly subordinate to WS-blooms. Measurements of chlorophyll in two size-fractions of the phytoplankton, netplankton (>20 um) and nanoplankton «20 um), revealed that the smaller nanoplankton are responsible for most of the phytoplankton biomass on the northeast U.S. shelf. Netplankton tend to be more abundant in nearshore areas of the MAB and shallow water on Georges Bank, where chlorophyll a is usually high; nanoplankton dominate deeper water at the shelf-break and deep water in the Gulf of Maine, where Chlw is usually low. As a general rule, the percent of phytoplankton in the netplankton size-fraction increases with increasing depth below surface and decreases proceeding offshore. There are distinct seasonal and regional patterns in the vertical distribution of chlorophyll a and percent netplankton, as revealed in composite vertical profiles of chlorophyll a constructed for 11 layers of the water column. Subsurface chlorophyll a maxima are ubiquitous during summer in stratified water. Chlorophyll a in the subsurface maximum layer is generally 2-8 times the concentration in the overlying and underlying water and approaches 50 to 75% of the levels observed in surface water during WS-bloom. The distribution of the ratio of the subsurface maximum chlorophyll a to surface chlorophyll a (SSR) during summer parallels the shelfwide pattern for stability, indexed as the difference in density (sigma-t) between 40 m and surface (stability 40. The weakest stability and lowest SSR's are found in shallow tidally-mixed water on Georges Bank; the greatest stability and highest SSR's (8-12:1) are along the mid and outer MAB shelf, over the winter residual water known as the "cold band." On Georges Bank, the distribution of SSR and the stability40 are roughly congruent with the pattern for maximum surface tidal current velocity, with values above 50 cms-1 defining SSR's less than 2:1 and the well-mixed area. Physical factors (bathymetry, vertical mixing by strong tidal currents, and seasonal and regional differences in the intensity and duration of vertical stratification) appear to explain much of the variability in phytoplankton chlorophyll a throughout this ecosystem. (PDF file contains 126 pages.)

Emergence of Cooperation in Heterogeneous Population: A Discrete-Time Replicator Dynamics Analysis

The emergence of cooperation is analyzed in heterogeneous populations where individuals can be classified in two groups according to their phenotypic appearance. Phenotype recognition is assumed for all individuals: individuals are able to identify the type of every other individual, but fail to recognize their own type, and thus behave under partial information conditions. The interactions between individuals are described by 2 × 2 symmetric games where individuals can either cooperate or defect. The evolution of such populations is studied in the framework of evolutionary game by means of the replicator dynamics. Overlapping generations are considered, so the replicator equations are formulated in discrete-time form. The well-posedness conditions of the system are derived. Depending on the parameters of the game, a restriction may exist for the generation length. The stability analysis of the dynamical system is carried out and a detailed description of the behavior of trajectories starting from the interior of the state-space is given. We find that, provided the conditions of well-posedness are verified, the linear stability of monomorphic states in the discrete-time replicator coincides with the one of the continuous case. Specific from the discrete-time case, a relaxed restriction for the generation length is derived, for which larger time-steps can be used without compromising the well-posedness of the replicator system.

Sliding Mode Control Strategy for Wind Turbine Power Maximization

The efficiency of the wind power conversions systems can be greatly improved using an appropriate control algorithm. In this work, a sliding mode control for variable speed wind turbine that incorporates a doubly fed induction generator is described. The electrical system incorporates a wound rotor induction machine with back-to-back three phase power converter bridges between its rotor and the grid. In the presented design the so-called vector control theory is applied, in order to simplify the electrical equations. The proposed control scheme uses stator flux-oriented vector control for the rotor side converter bridge control and grid voltage vector control for the grid side converter bridge control. The stability analysis of the proposed sliding mode controller under disturbances and parameter uncertainties is provided using the Lyapunov stability theory. Finally simulated results show, on the one hand, that the proposed controller provides high-performance dynamic characteristics, and on the other hand, that this scheme is robust with respect to the uncertainties that usually appear in the real systems.

Robust speed control for a variable speed wind turbine

Modern wind turbines are designed in order to work in variable speed opera-tions. To perform this task, these turbines are provided with adjustable speed generators, like the double feed induction generator (DFIG). One of the main advantages of adjustable speed generators is improving the system efficiency compared with _xed speed generators, because turbine speed can be adjusted as a function of wind speed in order to maximize the output power. However, this system requires a suitable speed controller in order to track the optimal reference speed of the wind turbine. In this work, a sliding mode control for variable speed wind turbines is proposed. The proposed design also uses the vector oriented control theory in order to simplify the DFIG dynamical equations. The stability analysis of the proposed controller has been carried out under wind variations and pa-rameter uncertainties using the Lyapunov stability theory. Finally, the simulated results show on the one hand that the proposed controller provides a high-performance dynamic behavior, and on the other hand that this scheme is robust with respect to parameter uncertainties and wind speed variations, which usually appear in real systems.

Sliding mode control law for a variable speedwind turbine

Modern wind turbines are designed in order to work in variable speed operations. To perform this task, wind turbines are provided with adjustable speed generators, like the double feed induction generator. One of the main advantage of adjustable speed generators is improving the system efficiency compared to fixed speed generators, because turbine speed can be adjusted as a function of wind speed in order to maximize the output power. However this system requires a suitable speed controller in order to track the optimal reference speed of the wind turbine. In this work, a sliding mode control for variable speed wind turbines is proposed. An integral sliding surface is used, because the integral term avoids the use of the acceleration signal, which reduces the high frequency components in the sliding variable. The proposed design also uses the vector oriented control theory in order to simplify the generator dynamical equations. The stability analysis of the proposed controller has been carried out under wind variations and parameter uncertainties by using the Lyapunov stability theory. Finally simulated results show, on the one hand that the proposed controller provides a high-performance dynamic behavior, and on the other hand that this scheme is robust with respect to parameter uncertainties and wind speed variations, that usually appear in real systems.

An adaptive sliding mode control law for the power maximization of the wind turbine system

POWERENG 2011

Energy Functions And Basic Equations For Ferroelectrics

Energy functions (or characteristic functions) and basic equations for ferroelectrics in use today are given by those for ordinary dielectrics in the physical and mechanical communications. Based on these basic equations and energy functions, the finite element computation of the nonlinear behavior of the ferroelectrics has been carried out by several research groups. However, it is difficult to process the finite element computation further after domain switching, and the computation results are remarkably deviating from the experimental results. For the crack problem, the iterative solution of the finite element calculation could not converge and the solutions for fields near the crack tip oscillate. In order to finish the calculation smoothly, the finite element formulation should be modified to neglect the equivalent nodal load produced by spontaneous polarization gradient. Meanwhile, certain energy functions for ferroelectrics in use today are not compatible with the constitutive equations of ferroelectrics and need to be modified. This paper proposes a set of new formulae of the energy functions for ferroelectrics. With regard to the new formulae of the energy functions, the new basic equations for ferroelectrics are derived and can reasonably explain the question in the current finite element analysis for ferroelectrics.

Numerical simulation of condensation of bubbles under microgravity conditions by moving mesh method in the double-staggered grid

A numerical 2D method for simulation of two-phase flows including phase change under microgravity conditions is presented in this paper, with a level set method being coupled with the moving mesh method in the double-staggered grid systems. When the grid lines bend very much in a curvilinear grid, great errors may be generated by using the collocated grid or the staggered grid. So the double-staggered grid was adopted in this paper. The level set method is used to track the liquid-vapor interface. The numerical analysis is fulfilled by solving the Navier-Stokes equations using the SIMPLER method, and the surface tension force is modeled by a continuum surface force approximation. A comparison of the numerical results obtained with different numerical strategies shows that the double-staggered grid moving-mesh method presented in this paper is more accurate than that used previously in the collocated grid system. Based on the method presented in this paper, the condensation of a single bubble in the cold water under different level of gravity is simulated. The results show that the condensation process under the normal gravity condition is different from the condensation process under microgravity conditions. The whole condensation time is much longer under the normal gravity than under the microgravity conditions.

A numerical investigation on the flow physics of koi carp (Cyprinus carpio koi) burst-and-coast swimming

Abstract—Burst-and-coast is the most common locomotion type in freely routine swimming of koi carps (Cyprinus carpio koi), which consists of a burst phase and a coast phase in each cycle and mostly leads to a straight-line trajectory. Combining with the tracking experiment, the flow physics of koi carp’s burst-andcoast swimming is investigated using a novel integrated CFD method solving the body-fluid interaction problem. The dynamical equations of a deforming body are formulated. Following that, the loose-coupled equations of the body dynamics and the fluid dynamics are numerically solved with the integrated method. The two burst modes, MT (Multiple Tail-beat) and HT (Half Tail-beat), which have been reported by the experiments, are investigated by numerical simulations in this paper. The body kinematics is predicted and the flow physics is visualized, which are in good agreement with the corresponding experiments. Furthermore, the optimization on the energy cost and several critical control mechanisms in burst-and-coast swimming of koi carps are explored, by varying the parameters in its selfpropelled swimming. In this paper, energetics is measured by the two mechanical quantities, total output power CP and Froude efficiency Fr. Results and discussion show that from the standpoint of mechanical energy, burst-and-coast swimming does not actually save energy comparing with steady swimming at the same average speed, in that frequently changing of speed leads to decrease of efficiency.