The use of object oriented modelling and simulation tools in constraint satisfaction for embodiment design

The software package Dymola, which implements the new, vendor-independent standard modelling language Modelica, exemplifies the emerging generation of object-oriented modelling and simulation tools. This paper shows how, in addition to its simulation capabilities, it may be used as an embodiment design tool, to size automatically a design assembled from a library of generic parametric components. The example used is a miniature model aircraft diesel engine. To this end, the component classes contain extra algebraic equations calculating the overload factor (or its reciprocal, the safety factor) for all the different modes of failure, such as buckling or tensile yield. Thus the simulation results contain the maximum overload or minimum safety factor for each failure mode along with the critical instant and the device state at which it occurs. The Dymola "Initial Conditions Calculation" function, controlled by a simple software script, may then be used to perform automatic component sizing. Each component is minimised in mass, subject to a chosen safety factor against failure, over a given operating cycle. Whilst the example is in the realm of mechanical design, it must be emphasised that the approach is equally applicable to the electrical or mechatronic domains, indeed to any design problem requiring numerical constraint satisfaction.

Physical and numerical modeling of horizontal gas-hydrate wells

Methane hydrate, which is usually found under deep seabed or permafrost zones, is a potential energy resource for future years. Depressurization of horizontalwells bored in methane hydrate layer is considered as one possible method for hydrate dissociation and methane extraction from the hosting soil. Since hydrate is likely to behave as a bonding material to sandy soils, supported well construction is necessary to avoid wellcollapse due to the loss of the apparent cohesion during depressurization. This paper describes both physical and numerical modeling of such horizontal support wells. The experimental part involves depressurization of small well models in a large pressure cell, while the numerical part simulates the corresponding problem. While the experiment models simulate only gas saturated initial conditions, the numerical analysis simulates both gas-saturated and more realistic water-saturated conditions based on effective stress coupled flow-deformation formulation of these three phases. © 2006 Taylor & Francis Group, London.

2D PIV measurements in the near field of grid turbulence using stitched fields from multiple cameras

We present measurements of grid turbulence using 2D particle image velocimetry taken immediately downstream from the grid at a Reynolds number of Re M = 16500 where M is the rod spacing. A long field of view of 14M x 4M in the down- and cross-stream directions was achieved by stitching multiple cameras together. Two uniform biplanar grids were selected to have the same M and pressure drop but different rod diameter D and crosssection. A large data set (10 4 vector fields) was obtained to ensure good convergence of second-order statistics. Estimations of the dissipation rate ε of turbulent kinetic energy (TKE) were found to be sensitive to the number of meansquared velocity gradient terms included and not whether the turbulence was assumed to adhere to isotropy or axisymmetry. The resolution dependency of different turbulence statistics was assessed with a procedure that does not rely on the dissipation scale η. The streamwise evolution of the TKE components and ε was found to collapse across grids when the rod diameter was included in the normalisation. We argue that this should be the case between all regular grids when the other relevant dimensionless quantities are matched and the flow has become homogeneous across the stream. Two-point space correlation functions at x/M = 1 show evidence of complex wake interactions which exhibit a strong Reynolds number dependence. However, these changes in initial conditions disappear indicating rapid cross-stream homogenisation. On the other hand, isotropy was, as expected, not found to be established by x/M = 12 for any case studied. © Springer-Verlag 2012.

Self-similar breakup of near-inviscid liquids

The final stages of pinchoff and breakup of dripping droplets of near-inviscid Newtonian fluids are studied experimentally for pure water and ethanol. High-speed imaging and image analysis are used to determine the angle and the minimum neck size of the cone-shaped extrema of the ligaments attached to dripping droplets in the final microseconds before pinchoff. The angle is shown to steadily approach the value of 18.0 ±0.4, independently of the initial flow conditions or the type of breakup. The filament thins and necks following a τ2 /3 law in terms of the time remaining until pinchoff, regardless of the initial conditions. The observed behavior confirms theoretical predictions. © 2012 American Physical Society.

Self-similar breakup of near-inviscid liquids.

The final stages of pinchoff and breakup of dripping droplets of near-inviscid Newtonian fluids are studied experimentally for pure water and ethanol. High-speed imaging and image analysis are used to determine the angle and the minimum neck size of the cone-shaped extrema of the ligaments attached to dripping droplets in the final microseconds before pinchoff. The angle is shown to steadily approach the value of 18.0 ± 0.4°, independently of the initial flow conditions or the type of breakup. The filament thins and necks following a τ(2/3) law in terms of the time remaining until pinchoff, regardless of the initial conditions. The observed behavior confirms theoretical predictions.

Modelling of settlement induced building damage

This thesis focuses on the modelling of settlement induced damage to masonry buildings. In densely populated areas, the need for new space is nowadays producing a rapid increment of underground excavations. Due to the construction of new metro lines, tunnelling activity in urban areas is growing. One of the consequences is a greater attention to the risk of damage on existing structures. Thus, the assessment of potential damage of surface buildings has become an essential stage in the excavation projects in urban areas (Chapter 1). The current damage risk assessment procedure is based on strong simplifications, which not always lead to conservative results. Object of this thesis is the development of an improved damage classification system, which takes into account the parameters influencing the structural response to settlement, like the non-linear behaviour of masonry and the soil-structure interaction. The methodology used in this research is based on experimental and numerical modelling. The design and execution of an experimental benchmark test representative of the problem allows to identify the principal factors and mechanisms involved. The numerical simulations enable to generalize the results to a broader range of physical scenarios. The methodological choice is based on a critical review of the currently available procedures for the assessment of settlement-induced building damage (Chapter 2). A new experimental test on a 1/10th masonry façade with a rubber base interface is specifically designed to investigate the effect of soil-structure interaction on the tunnelling-induced damage (Chapter 3). The experimental results are used to validate a 2D semi-coupled finite element model for the simulation of the structural response (Chapter 4). The numerical approach, which includes a continuum cracking model for the masonry and a non-linear interface to simulate the soil-structure interaction, is then used to perform a sensitivity study on the effect of openings, material properties, initial damage, initial conditions, normal and shear behaviour of the base interface and applied settlement profile (Chapter 5). The results assess quantitatively the major role played by the normal stiffness of the soil-structure interaction and by the material parameters defining the quasi-brittle masonry behaviour. The limitation of the 2D modelling approach in simulating the progressive 3D displacement field induced by the excavation and the consequent torsional response of the building are overcome by the development of a 3D coupled model of building, foundation, soil and tunnel (Chapter 6). Following the same method applied to the 2D semi-coupled approach, the 3D model is validated through comparison with the monitoring data of a literature case study. The model is then used to carry out a series of parametric analyses on geometrical factors: the aspect ratio of horizontal building dimensions with respect to the tunnel axis direction, the presence of adjacent structures and the position and alignment of the building with respect to the excavation (Chapter 7). The results show the governing effect of the 3D building response, proving the relevance of 3D modelling. Finally, the results from the 2D and 3D parametric analyses are used to set the framework of an overall damage model which correlates the analysed structural features with the risk for the building of being damaged by a certain settlement (Chapter 8). This research therefore provides an increased experimental and numerical understanding of the building response to excavation-induced settlements, and sets the basis for an operational tool for the risk assessment of structural damage (Chapter 9).

Self-consistent analysis of double-delta-doped InAlAs/InGaAs/InP HEMTs

We have carried out a theoretical study of double-delta-doped InAlAs/InGaAs/InP high electron mobility transistor (HEMT) by means of the finite differential method. The electronic states in the quantum well of the HEMT are calculated self-consistently. Instead of boundary conditions, initial conditions are used to solve the Poisson equation. The concentration of two-dimensional electron gas (2DEG) and its distribution in the HEMT have been obtained. By changing the doping density of upper and lower impurity layers we find that the 2DEG concentration confined in the channel is greatly affected by these two doping layers. But the electrons depleted by the Schottky contact are hardly affected by the lower impurity layer. It is only related to the doping density of upper impurity layer. This means that we can deal with the doping concentrations of the two impurity layers and optimize them separately. Considering the sheet concentration and the mobility of the electrons in the channel, the optimized doping densities are found to be 5 x 10(12) and 3 x 10(12) cm(-2) for the upper and lower impurity layers, respectively, in the double-delta-doped InAlAs/InGaAs/InP HEMTs.

Dark incoherent spatial solitons in logarithmically saturable nonlinear media

This paper studies numerically the dark incoherent spatial solitons propagating in logarithmically saturable nonlinear media by using a coherent density approach and a split-step Fourier approach for the first time. Under odd and even initial conditions, a soliton triplet and a doublet are obtained respectively for given parameters. Simultaneously, coherence properties associated with the soliton triplet and doublet are discussed. In addition, if the values of the parameters are properly chosen, five and four splittings from the input dark incoherent spatial solitons can also form. Lastly, the grayness of the soliton triplet and that of the doublet are studied, in detail.

Kinetics and transport model for the chemical vapor epitaxy of GexSi1-x

A numerical model that combines mass transport and surface kinetics was applied, for the first time, to the chemical vapor epitaxy of GexSi1-x. The temperature, velocity and concentration fields were calculated from the conservation equations for energy, momentum and species coupled with the boundary conditions on the growth surface which were determined by surface kinetics. The deposition rates of Si and Ge were assumed to be limited, respectively, by surface kinetics and mass transport. A theoretical relation between the initial conditions and the Ge composition in the solid was established. The calculated growth rate as well as the Ge composition in the solid and its dependence on growth temperature agree well with experimental data.

Non-modal instabilities of two-dimensional disturbances in plane Couette flow of a power-law fluid

Instabilities of fluid flows have traditionally been investigated by normal mode analysis, i.e. by linearizing the equations of flow and testing for unstable eigenvalues of the linearized problem. However, the results of eigenvalue analysis agree poorly in many cases with experiments, especially for shear flows. In this paper we study the instabilities of two-dimensional Couette flow of a polymeric fluid in the framework of non-modal stability theory rather than normal mode analysis. A power-law model is used to describe the polymeric liquid. We focus on the response to external excitations and initial conditions by examining the pseudospectra structures and the transient energy growths. For both Newtonian and non-Newtonian flows, the results show that there can be a rather large transient growth even though the linear operator of Couette flow has no unstable eigenvalue. The effects of non-Newtonian viscosity on the transient behaviors are examined in this study. The results show that the "shear-thinning/shear-thickening" effect increases/decreases the amplitude of responses to external excitations and initial conditions. (C) 2010 Elsevier B.V. All rights reserved.

Instabilities and transient behaviors of a liquid film flowing down a porous inclined plane

The nonmodal linear stability of a falling film over a porous inclined plane has been investigated. The base flow is driven by gravity. We use Darcy's law to describe the flow in the porous medium. A simplified one-sided model is used to describe the fluid flow. In this model, the influence of the porous layer on the flow in the film can be identified by a parameter beta. The instabilities of a falling film have traditionally been investigated by linearizing the governing equations and testing for unstable eigenvalues of the linearized problem. However, the results of eigenvalue analysis agree poorly in many cases with experiments, especially for shear flows. In the present paper, we have studied the linear stability of three-dimensional disturbances using the nonmodal stability theory. Particular attentions are paid to the transient behavior rather than the long time behavior of eigenmodes predicted by traditional normal mode analysis. The transient behaviors of the response to external excitations and the response to initial conditions are studied by examining the pseudospectral structures and the energy growth function G(t) Before we study the nonmodal stability of the system, we extend the results of long-wave analysis in previous works by examining the linear stabilities for streamwise and spanwise disturbances. Results show that the critical conditions of both the surface mode and the shear mode instabilities are dependent on beta for streamwise disturbances. However, the spanwise disturbances have no unstable eigenvalue. 2010 American Institute of Physics. [doi:10.1063/1.3455503]

Properties and applications of cold supersonic gas jet

By analyzing the formation mechanism of a supersonic gas jet, a set of equations which describe the atomic beam properties were established. The influence of initial temperature, initial pressure, background gas pressure and pumping speed was discussed in detail. A simulation program was developed based on the equations, and the results under different initial conditions were obtained. The results are in good agreement with the experimental data, and suggest that, in order to get much smaller transverse momentum in collision experiments, it is necessary to lower the initial temperature and the initial pressure of the supersonic gas jet, together with increasing the pumping speed. These results are very instructive for construction of a new generation of cold supersonic gas jets.

激光微地貌扫描仪的开发研制及在坡面侵蚀研究应用初步

利用光线的反射与聚焦成象原理 ,将地表高程转换成不同物象点位置的电信号 ,再经计算机处理成数字高程模型 ,即形成用于获取地表糙度数字高程模型的激光微地貌扫描仪。实验测量得到 :在相同初始条件下 ,不同坡度土槽放水和降雨实验后地表形态的数字地形模型。通过对侵蚀地表形态的分析 ,得出坡度对坡面侵蚀及沟蚀影响的某些初步规律。

Dynamics of nonlinear error growth and season-dependent predictability of El Nino events in the Zebiak-Cane model

With the intermediate-complexity Zebiak-Cane model, we investigate the 'spring predictability barrier' (SPB) problem for El Nino events by tracing the evolution of conditional nonlinear optimal perturbation (CNOP), where CNOP is superimposed on the El Nino events and acts as the initial error with the biggest negative effect on the El Nino prediction. We show that the evolution of CNOP-type errors has obvious seasonal dependence and yields a significant SPB, with the most severe occurring in predictions made before the boreal spring in the growth phase of El Nino. The CNOP-type errors can be classified into two types: one possessing a sea-surface-temperature anomaly pattern with negative anomalies in the equatorial central-western Pacific, positive anomalies in the equatorial eastern Pacific, and a thermocline depth anomaly pattern with positive anomalies along the Equator, and another with patterns almost opposite to those of the former type. In predictions through the spring in the growth phase of El Nino, the initial error with the worst effect on the prediction tends to be the latter type of CNOP error, whereas in predictions through the spring in the decaying phase, the initial error with the biggest negative effect on the prediction is inclined to be the former type of CNOP error. Although the linear singular vector (LSV)-type errors also have patterns similar to the CNOP-type errors, they cover a more localized area than the CNOP-type errors and cause a much smaller prediction error, yielding a less significant SPB. Random errors in the initial conditions are also superimposed on El Nino events to investigate the SPB. We find that, whenever the predictions start, the random errors neither exhibit an obvious season-dependent evolution nor yield a large prediction error, and thus may not be responsible for the SPB phenomenon for El Nino events. These results suggest that the occurrence of the SPB is closely related to particular initial error patterns. The two kinds of CNOP-type error are most likely to cause a significant SPB. They have opposite signs and, consequently, opposite growth behaviours, a result which may demonstrate two dynamical mechanisms of error growth related to SPB: in one case, the errors grow in a manner similar to El Nino; in the other, the errors develop with a tendency opposite to El Nino. The two types of CNOP error may be most likely to provide the information regarding the 'sensitive area' of El Nino-Southern Oscillation (ENSO) predictions. If these types of initial error exist in realistic ENSO predictions and if a target method or a data assimilation approach can filter them, the ENSO forecast skill may be improved. Copyright (C) 2009 Royal Meteorological Society

平面自由运动半物理仿真系统的初始条件实现

提出了一种平面自由运动半物理仿真系统的设计,对系统仿真试验的初始条件形成过程进行了具体阐述。系统中应用了一种平面两自由度直角坐标运动装置,用于完成系统驱动及速度和位置控制,同时这种运动装置能够进行直角坐标跟随运动,并与其末端测量机构共同实现对平面自由运动物体的高精度大范围位姿测量。建立了系统的运动学模型,并设计相应的控制算法实现试验所需的运动过程,采用半闭环和闭环相结合的方法有效控制系统的末端累积误差;针对提出的复合测量方法,建立了测量原理的数学模型,并进行了精度分析和仿真计算。实验证明这种设计和相应的控制测量方法合理可行。