透平叶栅三维粘性气动反问题的控制理论方法

将基于控制理论的形状优化设计方法应用于粘性可压流动条件下的透平叶栅三维气动反设计,详细推导了三维N-S方程伴随系统的偏微分方程组及其各类边界条件.讨论了伴随系统的解的适定性条件,并由此给出应用N-S方程进行气动优化的目标函数的选取限制.研究了伴随方程的数值求解技术,给出敏感性导数的最终计算式,结合拟牛顿算法发展了三维透平叶栅粘性反问题的气动设计方法.     

The analytical transfer matrix method for quantum reflection

In this paper, the analytical transfer matrix method (ATMM) is applied to study the properties of quantum reflection in three systems: a sech$^{2}$ barrier, a ramp potential and an inverse harmonic oscillator. Our results agree with those obtained by Landau and Lifshitz [Landau L D and Lifshitz E M 1977 \wx{Quantum Mechanics (Non-relativistic Theory)}{} (New York: Pergamon)], which proves that ATMM is a simple and effective method for quantum reflection.     

Statistics of resonances and of delay times in quasiperiodic Schrodinger equations

We study the distributions of the resonance widths P(gamma) and of delay times P(tau) in one-dimensional quasiperiodic tight-binding systems at critical conditions with one open channel. Both quantities are found to decay algebraically as gamma(-alpha) and tau(-gamma) on small and large scales, respectively. The exponents alpha and gamma are related to the fractal dimension D(E)(0) of the spectrum of the closed system as alpha = 1+D(E)(0) and gamma = 2-D(E)(0). Our results are verified for the Harper model at the metal-insulator transition and for Fibonacci lattices.     

Data completion method for the characterization of sound sources

This paper presents a different approach to solve the inverse acoustic problem. This problem is an "ill-posed" problem since the solution is very sensitive to measurement precision. A classical way to solve this problem consists in inversing a propagation operator which relates structure quantities (acoustic pressures or gradients) to near-field quantities (acoustic pressures or gradients). This can be achieved by using near-field acoustical holography (NAH) in separable coordinate systems. In order to overcome this limitation, the inverse boundary element method (IBEM) can be implemented to recover all acoustic quantities in a three-dimensional space and on an arbitrary three-dimensional source surface. In this paper, the data completion method (DCM) is developed: the acoustic gradients and pressures are known on a surface surrounding the source, but are unknown on its structure. The solution is given by the resolution of the Helmholtz formulation applied on the empty domain between the two boundaries made by the measurements quantities and the structure of the source. The conventional method applies directly the integral formulation for the empty domain. Another way of solving this Helmholtz formulation can be achieved by splitting it in two well-posed subproblems in a Steklov-Poincare?'s formulation. The data completion method allows one to solve the problem with acoustic perturbations due to sources on the exterior domain, or due to a confined domain, without altering the results.     

Improved adaptive model control

The Adaptive Model Control methodology is based on obtaining the model of a plant adaptively by using a FIR (Finite Impulse Response) filter, and inverting this model algebraically for use as the controller. This will practically cancel the plant dynamics when connected in cascade with the plant. The Adaptive Inverse Control methodology has the same aim but a different means for obtaining the plant inverse; the inverse is obtained by a second adaptive process. Adaptive Inverse Control systems have been studied with the "disturbance cancelation" and the "reference model" features added, while the algebraically inverting Adaptive Model Control systems could also support the same advanced features with a simpler and faster adapting structure. In this work, a scheme for adding the above features to an adaptive Model Control system is proposed, and using the A-4D aircraft as a test platform in simulations, the performance of this system is evaluated and compared to that of the Adaptive Inverse Control under the same operating conditions.     

Stability of inverse bicontinuous cubic phases in lipid-water mixtures

We investigate the stability of seven inverse bicontinuous cubic phases [ G, D, P, C(P), S, I-WP, F-RD] in lipid-water mixtures based on a curvature model of membranes. Lipid monolayers are described by parallel surfaces to triply periodic minimal surfaces. The phase behavior is determined by the distribution of the Gaussian curvature on the minimal surface and the porosity of each structure. Only G, D, and P are found to be stable, and to coexist along a triple line. The calculated phase diagram agrees very well with experimental results for 2:1 lauric acid/DLPC.     

Randomly chosen chaotic maps can give rise to nearly ordered behavior

Parrondo’s paradox [J.M.R. Parrondo, G.P. Harmer, D. Abbott, New paradoxical games based on Brownian ratchets, Phys. Rev. Lett. 85 (2000), 5226–5229] (see also [O.E. Percus, J.K. Percus, Can two wrongs make a right? Coin-tossing games and Parrondo’s paradox, Math. Intelligencer 24 (3) (2002) 68–72]) states that two losing gambling games when combined one after the other (either deterministically or randomly) can result in a winning game: that is, a losing game followed by a losing game = a winning game. Inspired by this paradox, a recent study [J. Almeida, D. Peralta-Salas, M. Romera, Can two chaotic systems give rise to order? Physica D 200 (2005) 124–132] asked an analogous question in discrete time dynamical system: can two chaotic systems give rise to order, namely can they be combined into another dynamical system which does not behave chaotically? Numerical evidence is provided in [J. Almeida, D. Peralta-Salas, M. Romera, Can two chaotic systems give rise to order? Physica D 200 (2005) 124–132] that two chaotic quadratic maps, when composed with each other, create a new dynamical system which has a stable period orbit. The question of what happens in the case of random composition of maps is posed in [J. Almeida, D. Peralta-Salas, M. Romera, Can two chaotic systems give rise to order? Physica D 200 (2005) 124–132] but left unanswered. In this note we present an example of a dynamical system where, at each iteration, a map is chosen in a probabilistic manner from a collection of chaotic maps. The resulting random map is proved to have an infinite absolutely continuous invariant measure (acim) with spikes at two points. From this we show that the dynamics behaves in a nearly ordered manner. When the foregoing maps are applied one after the other, deterministically as in [O.E. Percus, J.K. Percus, Can two wrongs make a right? Coin-tossing games and Parrondo’s paradox, Math. Intelligencer 24 (3) (2002) 68–72], the resulting composed map has a periodic orbit which is stable.     

Reconstructing the density of states by history-dependent metadynamics

We present a novel method for the calculation of the energy density of states D(E) for systems described by classical statistical mechanics. The method builds on an extension of a recently proposed strategy that allows the free-energy profile of a canonical system to be recovered within a preassigned accuracy [Proc. Natl. Acad. Sci. U.S.A. 99, 12562 (2002)]]. The method allows a good control over the error on the recovered system entropy. This fact is exploited to obtain D(E) more efficiently by combining measurements at different temperatures. The accuracy and efficiency of the method are tested for the two-dimensional Ising model (up to size 50 x 50) by comparison with both exact results and previous studies. This method is a general one and should be applicable to more realistic model systems.     

Full three-dimensional imaging via ground penetrating radar: assessment in controlled conditions and on field for archaeological prospecting

This paper deals with an advanced microwave tomographic approach capable of providing full 3D images of buried targets from scattered field data gathered by means of Ground Penetrating Radar (GPR) systems. The approach is based on an approximated model of the scattering phenomenon and it is capable of accounting for the vectorial nature of the interactions occurring between electromagnetic waves and probed materials. Moreover, the Truncated Singular Value Decomposition inversion scheme is exploited to solve the involved linear inverse scattering problem in a stable and accurate way. The advantages offered by the full 3D inversion algorithm with respect to a commonly adopted strategy, which produces 3D images by interpolating 2D reconstructions, are assessed against GPR data gathered in laboratory controlled conditions. Moreover, to provide an example of the full 3D imaging capabilities in on field conditions, we report on a GPR measurement campaign carried out at Grotte dell’Angelo, Pertosa, (SA), Southern Italy, one of the most famous sites of the Cilento and Vallo di Diano geopark.     

Hyperfine interaction of defects in insulators

The study of point defects, other than impurities, in insulators by electron paramagnetic resonance (E.P.R.) and electron nuclear double resonance (E.N.D.O.R.) is discussed. An introduction to the theory of the hyperfine spectra as obtained by E.P.R. and E.N.D.O.R. is presented in order to provide a basis for the discussion of results and their interpretation. Emphasis is placed on how the theory is used to analyse the spectra and how characteristics of the defects such as electronic state and defect geometry are obtained from the analyses. A résumé of the experimental techniques, especially E.N.D.O.R., is presented to give the reader a reasonably current view of the methods for obtaining data. The presentation of experimental results and recent theoretical developments is divided into a discussion of electron excess defects in cubic lattices and one on trapped hole centres in the same lattice type. This is followed by a discussion of these types of defects in non-cubic lattices. Finally the results of measurements on the hyperfine structure of excited states of electron and hole type in cubic crystals are discussed. In all cases an attempt was made to present the details of an illustrative example in order to indicate the salient features of the topic rather than give an exhaustive discussion. Needs for further measurement and future improvements of the theory are indicated.     

The Eighth Liquid Matter Conference

The Eighth Liquid Matter Conference (LMC8) was held at the Universit?t Wien from 6-10 September 2011. Initiated in 1990, the conferences of this series cover a broad range of highly interdisciplinary topics, ranging from simple liquids to soft matter and biophysical systems. The vast spectrum of scientific subjects presented and discussed at the LMC8 is reflected in the themes of the ten symposia: Ionic and quantum liquids, liquid metals Water, solutions and reaction dynamics Liquid crystals Polymers, polyelectrolytes, biopolymers Colloids Films, foams, surfactants, emulsions, aerosols Confined fluids, interfacial phenomena Supercooled liquids, glasses, gels Non-equilibrium systems, rheology, nanofluids Biofluids, active matter This special issue contains scientific papers, authored by participants of the LMC8, which provide a cross-section of the scientific activities in current liquid matter science, as discussed at the conference, and demonstrate the scientific as well as methodological progress made in this field over the past couple of years. The Eighth Liquid Matter Conference contents The Eighth Liquid Matter ConferenceChristoph Dellago, Gerhard Kahl and Christos N Likos Comparing light-induced colloidal quasicrystals with different rotational symmetriesMichael Schmiedeberg and Holger Stark Hydrogen bond network relaxation in aqueous polyelectrolyte solutions: the effect of temperatureS Sarti, D Truzzolillo and F Bordi Equilibrium concentration profiles and sedimentation kinetics of colloidal gels under gravitational stressS Buzzaccaro, E Secchi, G Brambilla, R Piazza and L Cipelletti The capillary interaction between two vertical cylindersHimantha Cooray, Pietro Cicuta and Dominic Vella Hydrodynamic and viscoelastic effects in polymer diffusionJ Farago, H Meyer, J Baschnagel and A N Semenov A density-functional theory study of microphase formation in binary Gaussian mixturesM Carta, D Pini, A Parola and L Reatto Microcanonical determination of the interface tension of flat and curved interfaces from Monte Carlo simulationsA Tr?ster and K Binder Phase diagrams of particles with dissimilar patches: X-junctions and Y-junctionsJ M Tavares and P I C Teixeira The unbearable heaviness of colloids: facts, surprises, and puzzles in sedimentationRoberto Piazza, Stefano Buzzaccaro and Eleonora Secchi Exploring water and other liquids at negative pressureFrédéric Caupin, Arnaud Arvengas, Kristina Davitt, Mouna El Mekki Azouzi, Kirill I Shmulovich, Claire Ramboz, David A Sessoms and Abraham D Stroock The configurational space of colloidal patchy polymers with heterogeneous sequencesIvan Coluzza and Christoph Dellago Repeated sorption of water in SBA-15 investigated by means of in situ small-angle x-ray scatteringM Erko, D Wallacher, G H Findenegg and O Paris Transition of the hydration state of a surfactant accompanying structural transitions of self-assembled aggregatesM Hishida and K Tanaka The effects of topology on the structural, dynamic and mechanical properties of network-forming materialsMark Wilson Surface tension of an electrolyte-air interface: a Monte Carlo studyAlexandre Diehl, Alexandre P dos Santos and Yan Levin Water and other tetrahedral liquids: order, anomalies and solvationB Shadrack Jabes, Divya Nayar, Debdas Dhabal, Valeria Molinero and Charusita Chakravarty Diffusion coefficient and shear viscosity of rigid water modelsSami Tazi, Alexandru Bo?an, Mathieu Salanne, Virginie Marry, Pierre Turq and Benjamin Rotenberg Phase behaviour of colloidal assemblies on 2D corrugated substratesSamir El Shawish, Emmanuel Trizac and Jure Dobnikar Structural properties of dendrimer-colloid mixturesDominic A Lenz, Ronald Blaak and Christos N Likos Fluid-fluid demixing of off-critical colloid-polymer systems confined between parallel platesE A G Jamie, R P A Dullens and D G A L Aarts Simulations of nematic homopolymer melts using particle-based models with interactions expressed through collective variablesKostas Ch Daoulas, Victor Rühle and Kurt Kremer Smectic shellsTeresa Lopez-Leon, Alberto Fernandez-Nieves, Maurizio Nobili and Christophe Blanc Intrinsic profiles and the structure of liquid surfacesP Tarazona, E Chacón and F Bresme Competing ordered structures formed by particles with a regular tetrahedral patch decorationGünther Doppelbauer, Eva G Noya, Emanuela Bianchi and Gerhard Kahl Heterogeneous crystallization in colloids and complex plasmas: the role of binary mobilitiesH L?wen, E Allahyarov, A Ivlev and G E Morfill Isotope effects in water as investigated by neutron diffraction and path integral molecular dynamicsAnita Zeidler, Philip S Salmon, Henry E Fischer, J?rg C Neuefeind, J Mike Simonson and Thomas E Markland Confined cubic blue phases under shearO Henrich, K Stratford, D Marenduzzo, P V Coveney and M E Cates Depletion-induced biaxial nematic states of boardlike particlesS Belli, M Dijkstra and R van Roij Active Brownian motion tunable by lightIvo Buttinoni, Giovanni Volpe, Felix Kümmel, Giorgio Volpe and Clemens Bechinger Structure and stability of charged clustersMark A Miller, David A Bonhommeau, Christopher J Heard, Yuyoung Shin, Riccardo Spezia and Marie-Pierre Gaigeot Non-equilibrium relaxation and tumbling times of polymers in semidilute solutionChien-Cheng Huang, Gerhard Gompper and Roland G Winkler Thermophoresis of colloids by mesoscale simulationsDaniel Lüsebrink, Mingcheng Yang and Marisol Ripoll Computing the local pressure in molecular dynamics simulationsThomas W Lion and Rosalind J Allen Gradient-driven fluctuations in microgravityA Vailati, R Cerbino, S Mazzoni, M Giglio, C J Takacs and D S Cannell.     

O(1D)+N2O→NO+NO反应的理论研究

应用三体模型及扩展的LEPS势能面(PES),对初始条件为(Ecol=55 kJ/mol,v=0,j=0)的O(1D)+N2O→NO+NO反应体系进行了准经典轨线(QCT)计算.根据计算结果对体系的势能面及反应机理进行详细的分析和讨论,较全面地研究了此反应体系的动力学特征.     

Inverse statistics of smooth signals: the case of two dimensional turbulence

The problem of inverse statistics (statistics of distances for which the signal fluctuations are larger than a certain threshold) in differentiable signals with power law spectrum, E(k) approximately k(-alpha), 3< or =alpha<5, is discussed. We show that for these signals, with random phases, exit-distance moments follow a bifractal distribution. We also investigate two dimensional turbulent flows in the direct cascade regime, which display a more complex behavior. We give numerical evidences that the inverse statistics of 2D turbulent flows is described by a multifractal probability distribution; i.e., the statistics of laminar events is not simply captured by the exponent alpha characterizing the spectrum.     

Parameter scaling in the decoherent quantum-classical transition for chaotic systems

The quantum to classical transition for a system depends on many parameters, including a scale length for its action, variant Planck's over 2 pi, a measure of its coupling to the environment, D, and, for chaotic systems, the classical Lyapunov exponent, lambda. We propose measuring the proximity of quantum and classical evolutions as a multivariate function of (Planck's over 2 pi,lambda,D) and searching for transformations that collapse this hypersurface into a function of a composite parameter zeta= Planck's over 2 pi alpha)lambda beta D gamma. We report results for the quantum Cat Map and Duffing oscillator, showing accurate scaling behavior over a wide parameter range, indicating that this may be used to construct universality classes for this transition.     

Hyperfine interactions in chromium

Conclusion Our experiments suggest that Tantalum in Chromium is non magnetic and does not clamp the S.D.W. The various transitions of chromium have been observed. On the other hand rhodium seems to couple with the S.D.W. but the interaction, if any, between rhodium and the S.D.W. does not strongly modify the shape of the S.D.W. TDPAC experiments are therefore a useful tool to investigate the magnetic properties of impurities dissolved in a chromium matrix. This is extremely interesting because the formation of a localized moment in an itinerant antiferromagnet as well as their interaction with the S.D.W.'s have never been studied. It is also interesting to note that S.D.W.'s could couple magnetic impurities at very large distances and then could give rise to a special type of spin glass.     

Profile Reconstruction in Dynamic Infrared Thermography: Adaptability and Stability

An iterative algorithm for solving geometric inverse problem in dynamic infrared thermography has been proposed in this paper. According to 1D dynamic heat-conduction model and the corrective criterion based on Tailor's series expansion, both 2D and 3D geometrical inverse problems are solved by following an iterative procedure. Computer simulation results show that the reconstructed profile has satisfactory accuracy. The adaptability and stability of the method are also discussed.     

Measurement of spin-lattice relaxation times and concentrations in systems with chemical exchange using the one-pulse sequence: breakdown of the Ernst model for partial saturation in nuclear magnetic resonance spectroscopy

A fundamental problem in Fourier transform NMR spectroscopy is the calculation of observed resonance amplitudes for a repetitively pulsed sample, as first analyzed by Ernst and Anderson in 1966. Applications include determination of spin-lattice relaxation times (T(1)'s) by progressive saturation and correction for partial saturation in order to determine the concentrations of the chemical constituents of a spectrum. Accordingly, the Ernst and Anderson formalism has been used in innumerable studies of chemical and, more recently, physiological systems. However, that formalism implicitly assumes that no chemical exchange occurs. Here, we present an analysis of N sites in an arbitrary chemical exchange network, explicitly focusing on the intermediate exchange rate regime in which the spin-lattice relaxation rates and the chemical exchange rates are comparable in magnitude. As a special case of particular importance, detailed results are provided for a system with three sites undergoing mutual exchange. Specific properties of the N-site network are then detailed. We find that (i) the Ernst and Anderson analysis describing the response of a system to repetitive pulsing is inapplicable to systems with chemical exchange and can result in large errors in T(1) and concentration measurements; (ii) T(1)'s for systems with arbitrary exchange networks may still be correctly determined from a one-pulse experiment using the Ernst formula, provided that a short interpulse delay time and a large flip angle are used; (iii) chemical concentrations for exchanging systems may be correctly determined from a one-pulse experiment either by using a short interpulse delay time with a large flip angle, as for measuring T(1)'s, and correcting for partial saturation by use of the Ernst formula, or directly by using a long interpulse delay time to avoid saturation; (iv) there is a significant signal-to-noise penalty for performing one-pulse experiments under conditions which permit accurate measurements of T(1)'s and chemical concentrations. The present results are analogous to but are much more general than those that we have previously derived for systems with two exchanging sites. These considerations have implications for the design and interpretation of one-pulse experiments for all systems exhibiting chemical exchange in the intermediate exchange regime, including virtually all physiologic samples.     

Phase coexistence and interface structure of a two-component Lennard-Jones fluid in porous media: application of Born—Green—Yvon equation

The Born-Green-Yvon equation with smoothed density approximation is used to calculate the liquid-liquid density profiles of a symmetric Lennard-Jones fluid in a hard sphere disordered matrix. The phase diagrams are evaluated for model systems characterized by different matrix densities and compared with the results of theoretical predictions and the Monte Carlo simulations of Gordon, P. A., and Glandt, E. D., 1996, J. chem. Phys., 105, 4257. It was found that increasing the matrix packing fraction reduces the magnitude of the miscibility gap and smooths the density profiles between two coexisting phases.     

Numerical evaluation of linearized image reconstruction based on finite element method for biomedical photoacoustic imaging

An image reconstruction algorithm for biomedical photoacoustic imaging is discussed. The algorithm solves the inverse problem of the photoacoustic phenomenon in biological media and images the distribution of large optical absorption coefficients, which can indicate diseased tissues such as cancers with angiogenesis and the tissues labeled by exogenous photon absorbers. The linearized forward problem, which relates the absorption coefficients to the detected photoacoustic signals, is formulated by using photon diffusion and photoacoustic wave equations. Both partial differential equations are solved by a finite element method. The inverse problem is solved by truncated singular value decomposition, which reduces the effects of the measurement noise and the errors between forward modeling and actual measurement systems. The spatial resolution and the robustness to various factors affecting the image reconstruction are evaluated by numerical experiments with 2D geometry.     

Order-N sparse minimum-variance open-loop reconstructor for extreme adaptive optics

A scalable sparse minimum-variance open-loop wave-front reconstructor for extreme adaptive optics (ExAO) systems is presented. The reconstructor is based on Ellerbroek's sparse approximation of the wave-front inverse covariance matrix [J. Opt. Soc. Am. A 19, 1803 (2002)]. The baseline of the numerical approach is an iterative conjugate gradient (CG) algorithm for reconstructing a spatially sampled wave front at N grid points on a computational domain of size equal to the telescope's primary mirror's diameter D that uses a multigrid (MG) accelerator to speed up convergence efficiently and enhance its robustness. The combined MGCG scheme is order N and requires only two CG iterations to converge to the asymptotic average Strehl ratio (SR) and root-mean-square reconstruction error. The SR and reconstruction squared error are within standard deviation with figures obtained from a previously proposed MGCG fast-Fourier-transform based minimum-variance reconstructor that incorporates the exact wave-front inverse covariance matrix on a computational domain of size equal to 2D. A cost comparison between the present sparse MGCG algorithm and a Cholesky factorization based algorithm that uses a reordering scheme to preserve sparsity indicates that the latter method is still competitive for real-time ExAO wave-front reconstruction for systems with up to N approximately equal to 10(4) degrees of freedom because the update rate of the Cholesky factor is typically several orders of magnitude lower than the temporal sampling rate.