Shock wave reflection from a wedge in a dusty gas

The present paper contains a detailed study of shock wave reflection from a wedge placed in various suspensions. In past works, the incident shock propagated initially in pure gas and the suspension started only at the leading edge of the deflecting wedge. However, in the present case the entire flow field is filled with a gas-dust suspension and the initial shock wave has steady-state structure relative to the shock front. In former studies the transmitted shock wave starts its propagation into the suspension and is reflected from the wedge at the same time. It is therefore obvious that the two unrelated processes of (2D) reflection and (1D) "transitional" relaxation occur simultaneously. In the present case the suspension behind the incident shock wave has reached steady state (i.e., it is a traveling wave) before the shock reaches the wedge leading edge. The reflection process from the deflecting wedge is studied for different dust mass loadings and different dust-particle diameter. It is shown that when the dust loading is low and the dust particle diameter is small the wave reflection pattern is similar to that observed in a similar pure gas case. In addition, an equilibrium state is reached, behind the evolved waves, very quickly. On the other hand, when the dust loading is relatively high and/or the dust particle diameter is relatively large, the observed reflection wave pattern is very different from that seen in a similar pure gas case. In such cases it takes much longer time to reach an equilibrium state behind the reflecting waves. It is also shown that the dust presence significantly affects the (gas) pressure on the wedge surface. The higher the dust loading is, the higher the pressure on the wedge surface. Suspensions composed of solid particle of different size, but having the same dust mass loading, will approach the same equilibrium pressure. However, it will take longer time to reach an equilibrium state for suspensions having large diameter particles. (C) 2004 Elsevier Ltd. All rights reserved.

Concentration distribution in crystallization from solution under microgravity

Concentration distribution in crystallization from solution under microgravity is numerically studied. A quasi-steady state growth and dissolution in a 2D rectangular enclosure filled with sodium chlorate (NaClO3) aqueous solution, in which one wall is the growth surface of the crystal and the opposite one is the dissolution surface, is considered. The solute transport process at the growth surface is described by the diffusion-reaction theory with finite interface kinetics coefficient. The results show that the concentration at the growth surface is supersaturated and the supersaturation distribution is of non-uniformity, i.e. the supersaturation in a region facing an incoming flow is high. On the other hand, the non-uniformity of supersaturation at the growth surface is closely related to the gravity level even under microgravity, it exponentially increases as the thermal Rayleigh number on behalf of the gravity level rises.

Jet:progress in performance and understanding

In 1990 JET operated with a number of technical improvements which led to advances in performance and permitted the carrying out of experiments specifically aimed at improving physics understanding of selected topics relevant to the "NEXT STEP". The new facilities include beryllium antenna screens, a prototype lower hybrid current drive system, and modification of the NI system to enable the injection of He-3 and He-4. Continued investigation of the hot-ion H-mode produced a value of n(D)(0)tau-E(T)(i)(0) = 9 x 10(20)m-3s keV, which is near conditions required for Q(DT) = 1, while a new peaked density profile H-mode was developed with only slightly lower performance. Progress towards steady state operation has been made by achieving ELMy H-modes under certain operating conditions, while maintaining good tau-E values. Experimental simulation of He ash transport indicates effective removal of alpha-particles from the plasma core for both L and H mode plasmas. Detailed analyses of particle and energy transport have helped establish a firmer link between particle and energy transport, and have suggested a connection between reduced energy transport and reversed shear. Numerical and analytic studies of divertor physics carried out for the pumped divertor phase of JET have helped clarify the key parameters governing impurity retention, and an intensive model validation effort has begun. Experimental simulation of alpha-particle effects with beta-fast up to 8% have shown that the slowing down processes are classical, and have given no evidence of deleterious collective effects.

Kinetic theory for particle concentration distribution in two phase flow

The results of experiments in open channels and closed pipelines show two kinds of patterns for the vertical distribution of particle concentration (i.e., pattern I and pattern II). The former shows a pattern of maximum concentration at some location above the bottom and the downward decay of the concentration below the location. The latter always shows an increase of the particle concentration downward over the whole vertical, with the maximum value at the bottom. Many investigations were made on the pattern II, but few were made on pattern I. In this paper, a particle velocity distribution function is first obtained in the equilibrium state or in dilute steady state for the particle in two-phase flows, then a theoretical model for the particle concentration distribution is derived from the kinetic theory. More attention is paid to the predictions of the concentration distribution of pattern I and comparisons of the present model are made with the data measured by means of laser doppler anemometry (LDA). Very good agreements are obtained between the measured and calculated results.

Fiscal Policy and Economic Stability: Does PIGS stand for Procyclicality In Government Spending?

The Financial Crisis has hit particularly hard countries like Ireland or Spain. Procyclical fiscal policy has contributed to a boom-bust cycle that undermined fiscal positions and deepened current account deficits during the boom. We set up an RBC model of a small open economy, following Mendoza (1991), and introduce the effect of fiscal policy decisions that change over the cycle. We calibrate the model on data for Ireland, and simulate the effect of different spending policies in response to supply shocks. Procyclical fiscal policy distorts intertemporal allocation decisions. Temporary spending boosts in booms spur investment, and hence the need for external finance, and so generates very volatile cycles in investment and the current account. This economic instability is also harmful for the steady state level of output. Our model is able to replicate the relation between the degree of cyclicality of fiscal policy, and the volatility of consumption, investment and the current account observed in OECD countries.

Reference Points and Optimal Management in Stochastic Age-Structured Fisheries Models

The purpose of this article is to characterize dynamic optimal harvesting trajectories that maximize discounted utility assuming an age-structured population model, in the same line as Tahvonen (2009). The main novelty of our study is that uses as an age-structured population model the standard stochastic cohort framework applied in Virtual Population Analysis for fish stock assessment. This allows us to compare optimal harvesting in a discounted economic context with standard reference points used by fisheries agencies for long term management plans (e.g. Fmsy). Our main findings are the following. First, optimal steady state is characterized and sufficient conditions that guarantees its existence and uniqueness for the general case of n cohorts are shown. It is also proved that the optimal steady state coincides with the traditional target Fmsy when the utility function to be maximized is the yield and the discount rate is zero. Second, an algorithm to calculate the optimal path that easily drives the resource to the steady state is developed. And third, the algorithm is applied to the Northern Stock of hake. Results show that management plans based exclusively on traditional reference targets as Fmsy may drive fishery economic results far from the optimal.

Growth in Overlapping Generation Economies with Non-Renewable Resources

Published as an article in: Journal of Environmental Economics and Management, 2005, vol. 50, issue 2, pages 387-407.

The Evaluation of Fisheries Management: A Dynamic Stochastic Approach

In this article, we analyze how to evaluate fishery resource management under “ecological uncertainty”. In this context, an efficient policy consists of applying a different exploitation rule depending on the state of the resource and we could say that the stock is always in transition, jumping from one steady state to another. First, we propose a method for calibrating the growth path of the resource such that observed dynamics of resource and captures are matched. Second, we apply the calibration procedure proposed in two different fishing grounds: the European Anchovy (Division VIII) and the Southern Stock of Hake. Our results show that the role played by uncertainty is essential for the conclusions. For European Anchovy fishery (Division VIII) we find, in contrast with Del Valle et al. (2001), that this is not an overexploited fishing ground. However, we show that the Southern Stock of Hake is in a dangerous situation. In both cases our results are in accordance with ICES advice.

在激波波面中氧分子的非平衡离解

对强激波作用下双原子分子振动与离解耦合的非平衡离解过程进行了理论计算.本工作的特点是将计算起点建立在分子基本参数上,采用主方程理论处理振动与离解的耦合,振动跃迁几率用SSH理论计算,在离解限附近考虑多量子数跃迁并计及原子复合的影响.对O2-Ar体系,计算给出了在正激波后O2分子振动能级分布、振动弛豫时间、离解孕育时间、离解产物浓度、离解速率系数等物理量随时间的演化.计算结果分别与Camac和Wray的实验相符.计算显示,在激波作用的后期,有准稳态的振动能级布居分布.计算结果显示,Park模型低估了非平衡离解速率系数,Hansen模型则高估了非平衡离解速率系数.

Reference Points Based on Dynamic Optimisation: A Versatil Algorithm for Mixed Fishery Management with Bio-economic Agestructured Models

Single-species management objectives may not be consistent within mixed fisheries. They may lead species to unsafe situations, promote discarding of over-quota and/or misreporting of catches. We provide an algorithm for characterising bio-economic reference points for a mixed fishery as the steady-state solution of a dynamic optimal management problem. The optimisation problem takes into account: i) that species are fishing simultaneously in unselective fishing operations and ii)intertemporal discounting and fleet costs to relate reference points to discounted economic profits along optimal trajectories. We illustrate how the algorithm can be implemented by applying it to the European Northern Stock of Hake (Merluccius merluccius), where fleets also capture Northern megrim (Lepidorhombus whiffiagonis) and Northern anglerfish (Lophius piscatorius and Lophius budegassa). We find that optimal mixed management leads to a target reference point that is quite similar to the 2/3 of the Fmsy single-species (hake) target. Mixed management is superior to singlespecies management because it leads the fishery to higher discounted profits with higher long-term SSB for all species. We calculate that the losses due to the use of the Fmsy single-species (hake) target in this mixed fishery account for 11.4% of total discounted profits.

The Propagation Characteristics of Microwave in a Non-Equilibrium Plasma Layer

The influence of non-equilibrium plasma layer pressure and thickness on the transmission of microwave is considered when the incidence of wave is at an arbitrary angle. The plasma is cold, weakly ionized, and steady-state. It is assumed that it is a layered media with a kind of distribution of electron number density and the microwave is a plane wave. The results show that the pressure of plasma affects the absorption of microwave deeply, and the thickness relatively weakly in a non-equilibrium plasma slab.

3D Finite Volume Scheme for Czochralski Single Crystal Growth

Czochralski (Cz) technique, which is used for growing single crystals, has dominated the production of single crystals for electronic applications. The Cz growth process involves multiple phases, moving interface and three-dimensional behavior. Much has been done to study these phenomena by means of numerical methods as well as experimental observations. A three-dimensional curvilinear finite volume based algorithm has been developed to model the Cz process. A body-fitted transformation based approach is adopted in conjunction with a multizone adaptive grid generation (MAGG) technique to accurately handle the three-dimensional problems of phase-change in irregular geometries with free and moving surfaces. The multizone adaptive model is used to perform a three-dimensional simulation of the Cz growth of silicon single crystals.Since the phase change interface are irregular in shape and they move in response to the solution, accurate treatment of these interfaces is important from numerical accuracy point of view. The multizone adaptive grid generation (MAGG) is the appropriate scheme for this purpose. Another challenge encountered is the moving and periodic boundary conditions, which is essential to the numerical solution of the governing equations. Special treatments are implemented to impose the periodic boundary condition in a particular direction and to determine the internal boundary position and shape varying with the combination of ambient physicochemical transport process and interfacial dynamics. As indicated above that the applications and processes characterized by multi-phase, moving interfaces and irregular shape render the associated physical phenomena three-dimensional and unsteady. Therefore a generalized 3D model rather than a 2D simulation, in which the governing equations are solved in a general non-orthogonal coordinate system, is constructed to describe and capture the features of the growth process. All this has been implemented and validated by using it to model the low pressure Cz growth of silicon. Accuracy of this scheme is demonstrated by agreement of simulation data with available experimental data. Using the quasi-steady state approximation, it is shown that the flow and temperature fields in the melt under certain operating conditions become asymmetric and unsteady even in the absence of extrinsic sources of asymmetry. Asymmetry in the flow and temperature fields, caused by high shear initiated phenomena, affects the interface shape in the azimuthal direction thus results in the thermal stress distribution in the vicinity, which has serious implications from crystal quality point of view.

A Synoptic View of a Coastal Plain Estuary.

During October, 1972 the Patuxent River Estuary was monitored intensively and synoptically over two tidal cycles to determine the spatial and temporal patterns of various hydrodynamic, chemical and biological features. Forty-one depths at eleven stations along nine transects were sampled simultaneously at hourly intervals for salinity, temperature, dissolved oxygen, chlorohyll a, particulate nitrogen, nitrate, nitrite, total kjeldahl nitrogen, ammonia, particulate carbohydrate, dissolved organic carbon, total hydrolizable phosphorous, dissolved inorganic phosphorous, suspended sediment, particle size distribution, and zooplankton. Tidal velocity was continuously monitored at each depth by recording current meters. Riverine input and meteorological conditions were relatively stable for two weeks preceeding the deployment. This communication describes the calculation of the intrinsic rates of change of the observed variables from their measured distributions in the Estuary. The steady-state, one-dimensional equation of species continuity is employed to separate the advection and tidal dispersion of a hydrodynamically passive substance frbm its intrinsic rate of change at point. A new spatial transform is introduced for the purpose of interpolation and extrapolation of data.The intrinsic rate of change profiles reveal a region of heavy bloom activity in the upper estuary and a secondary bloom near the point in the River that most of the suspended material settles out. The changes in ammonia and nitrates are highly correlated to the productivity patterns. Phosphorous rates are less closely correlated to productivity. The perturbations that the Chalk Point steam electric power plant have on the heat and oxygen balances are easily discernible.

Information Preservation Modeling of Rayleigh-Benard Transition from Thermal Conduction to Convection

Onset and evolution of the Rayleigh-Benard (R-B) convection are investigated using the Information Preservation (IP) method. The information velocity and temperature are updated using the Octant Flux Splitting (OFS) model developed by Masters & Ye based on the Maxwell transport equation suggested by Sun & Boyd. Statistical noise inherent in particle approaches such as the direct simulation Monte Carlo (DSMC) method is effectively reduced by the IP method, and therefore the evolutions from an initial quiescent fluid to a final steady state are shown clearly. An interesting phenomenon is observed: when the Rayleigh number (Ra) exceeds its critical value, there exists an obvious incubation stage. During the incubation stage, the vortex structure clearly appears and evolves, whereas the Nusselt number (Nu) of the lower plate is close to unity. After the incubation stage, the vortex velocity and Nu rapidly increase, and the flow field quickly reaches a steady, convective state. A relation of Nu to Ra given by IP agrees with those given by DSMC, the classical theory and experimental data.

Investigations of cohesive zone properties of the interface between metal film and ceramic substrate at nano- and micro-scales

Cohesive zone characterizations of the interface between metal film and ceramic substrate at micro- and nano-scales are performed in the present research. At the nano-scale, a special potential for special material interface (Ag/MgO) is adopted to investigate the interface separation mechanism by using MD simulation, and stress-separation relationship will be obtained. At the micro-scale, peeling experiment is performed for the Al film/Al2O3 substrate system with an adhesive layer at the interface. Adhesive is a mixture of epoxy and polyimide with mass ratio 1:1, by which a brittle cohesive property is obtained. The relationships between energy release rate, the film thickness and the adhesive layer thickness are measured during the steady-state peeling. The experimental result has a similar trend as modeling result for a weak adhesion interface case.