Non linear waves in two dimensional keplerian flows

Abstract

Accretion discs in astrophysics are fundamental for converting gravitational binding energy into observed electromagnetic radiation. We study the behavior of waves in a two dimensional supersonic Keplerian flow inside a given gravitational potential. We present the effects of shearing and rotation on short waves, and the numerical study of the dynamical stability of such flows with respect to various perturbations. We show that a large class of dynamical effects, due to pressure and associated to short time scales, may be excited.     

The Relegation Algorithm

The relegation algorithm extends the method of normalization by Lie transformations. Given a Hamiltonian that is a power series = ₀+ ₁+ ... of a small parameter , normalization constructs a map which converts the principal part ₀into an integral of the transformed system — relegation does the same for an arbitrary function [G]. If the Lie derivative induced by [G] is semi-simple, a double recursion produces the generator of the relegating transformation. The relegation algorithm is illustrated with an elementary example borrowed from galactic dynamics; the exercise serves as a standard against which to test software implementations. Relegation is also applied to the more substantial example of a Keplerian system perturbed by radiation pressure emanating from a rotating source.This revised version was published online in October 2005 with corrections to the Cover Date.     

Normalization and the detection of integrability: The generalized Van Der Waals potential

Deprit and Miller have conjectured that normalization of integrable Hamiltonians may produce normal forms exhibiting degenerate equilibria to very high order. Several examples in the class of coupled elliptic oscillators are known. In order to test the utility of normalization as a detector of integrability we normalize, to high order, a perturbed Keplerian system known to have several integrable limits; the generalized van der Waals Hamiltonian for a hydrogen atom. While the separable limits give rise to high order degeneracy we find a non-separable, integrable limit for which the normal form does not exhibit degeneracy. We conclude that normalization may, in certain cases, indicate integrability but is not guaranteed to uncover all integrable limits.     

A sequential sparse polynomial chaos expansion using Bayesian regression for geotechnical reliability estimations

Polynomial chaos expansions (PCEs) have been widely employed to estimate failure probabilities in geotechnical engineering. However, PCEs suffer from two deficiencies: (a) PCE coefficients are solved by the least-square minimization method which easily causes overfitting issues; (b) building a high order PCE is often computationally expensive. In order to overcome the aforementioned drawbacks, the Bayesian regression technique is employed to evaluate PCE coefficients, which not only provides a sparse solution but also avoids overfitting. With the aid of the predictive means and variances given by Bayesian analysis, a learning function is proposed to sequentially select the most informative samples that are critical to build a PCE. This sequential learning scheme can highly enhance the computational efficiency of PCEs. Besides, importance sampling (IS) is incorporated into the sequential learning (SL)-PCEs to deal with geotechnical problems with small failure probabilities. The proposed method of SL-PCE-IS is applied to three illustrative examples, which shows that the improved PCE method is more effective and efficient than the common PCEs method, leading to accurate estimations of small failure probabilities using fewer training samples.     

没水倾斜板式防波堤消波性能分析

提出了一种将波浪中倾斜板问题等效化简为波浪中水平板单元组的方法,该方法建立在使用分离变量法求解水工结构边值问题的基础上,并使用伽辽金法精确求解连续边界条件,确定考虑衰减波态的速度势函数,从而求解没水倾斜板结构的消波性能。等效化简法计算精度于边界元法相当,且计算单元数量少、开销低。基于二维线性势波理论,对没水倾斜板式防波堤消波性能分析显示,没水板的倾斜角度、没水深度与板长是结构消波性能的控制因素:没水倾斜板防波堤的消波性能优于没水水平板防波堤的消波性能,随着没水板结构的倾斜角度增大,没水倾斜板结构的波浪透射系数显著减小,且长板优于短板,浅板优于深板;与前人的水槽实验对比显示,相对没水深度与波陡影响结构的消波性能,且波陡造成的波浪破碎贡献了显著的波能消耗。该结论对板式防波堤的结构配置、优化设计有重要意义。     

Buoyancy and The Sensible Heat Flux Budget Within Dense Canopies

In contrast to atmospheric surface-layer (ASL) turbulence, a linear relationship between turbulent heat fluxes (F_T) and vertical gradients of mean air temperature within canopies is frustrated by numerous factors, including local variation in heat sources and sinks and large-scale eddy motion whose signature is often linked with the ejection-sweep cycle. Furthermore, how atmospheric stability modifies such a relationship remains poorly understood, especially in stable canopy flows. To date, no explicit model exists for relating F_T to the mean air temperature gradient, buoyancy, and the statistical properties of the ejection-sweep cycle within the canopy volume. Using third-order cumulant expansion methods (CEM) and the heat flux budget equation, a “diagnostic” analytical relationship that links ejections and sweeps and the sensible heat flux for a wide range of atmospheric stability classes is derived. Closure model assumptions that relate scalar dissipation rates with sensible heat flux, and the validity of CEM in linking ejections and sweeps with the triple scalar-velocity correlations, were tested for a mixed hardwood forest in Lavarone, Italy. We showed that when the heat sources (S_T) and F_T have the same sign (i.e. the canopy is heating and sensible heat flux is positive), sweeps dominate the sensible heat flux. Conversely, if S_T and F_T are opposite in sign, standard gradient-diffusion closure model predict that ejections must dominate the sensible heat flux.     

地下水模型的Neumann展开Monte-Carlo随机有限元法

讨论了求解地下水模型(水流模型和水质模型)的Neumann展开Monte-Carlo随机有限元法。从基本的随机变量入手,避免了过程中随机变量的增多问题,给出结点水头(浓度)的均值、方差和水头在某区间的概率计算方法;改进了矩阵求逆的效率,对输入随机变量较多、随机变量变异较大的非稳定地下水问题特别有效。同时选取二维承压地下水水流问题(有解析解)作为例子,进行了随机数值模拟实验。      

Wave characteristics past a flexible fishnet

The scattering of surface waves by a flexible fishnet is studied analytically. The fishnet is modelled as a porous flexible barrier displaced solely by hydrodynamic force like a catenary. The objective is to investigate how a flexible permeable barrier affects the passing waves in the way they are transmitted and reflected, as observed by the fact that the water inside a fishfarm surrounded by fishnets is significantly calmer than that outside. The boundary value problems are solved by defining the reflection coefficient in terms of velocity potential and then the full solutions are obtained by suitable application of the eigenfunction expansion method and the least squares approximation method. The variations of the reflection coefficient, hydrodynamic pressure, barrier deformation and surface wave elevation are determined with respect to the barrier length, porosity and stiffness. It is observed that as the fishnet gets more flexible, its deformation increases and the reflection coefficient decreases, whereas as the fishnet gets more porous, more water can pass through it and thus the reflection coefficient, barrier deformation and the hydrodynamic force are reduced. The flexibility of the barrier behaves like its porosity by allowing more wave energy to act on it through its deformation and hence reduce the reflection and hydrodynamic force of the incident waves acting on the barrier.