医用CT仪上基于衰减的曝光量和峰值电压自动控制基本等式

为了给患者以最少的辐射剂量而不影响诊断,现行医用CT仪都有基于噪声指标的自动曝光量控制。峰值电压对于图像质量和辐射剂量影响都很大,因此自动峰值电压控制很重要。本文从图像质量和辐射剂量优化的第一原理出发,给出基于衰减的曝光量和峰值电压同时自动控制的一般等式。此等式可用于单独曝光量自动控制或者单独峰值电压自动控制。对于单独暴光量自动控制,此等式要求对患者体厚度毎增加一厘米,曝光量要增加3.8%。对于单独峰值电压自动控制,此等式要求对患者体厚度每增加一厘米,峰值电压要平均增加1.53%。如果患者体厚度是常数,此等式要求平均百分比曝光量增加是百分比峰值电压减少的2.49倍。这些等式应该可以在现行医用CT仪上用于曝光量和峰值电压的自动控制。      

等效水力隙宽和水力梯度对岩体裂隙网络非线性渗流特性的影响

等效水力隙宽和水力梯度是影响岩体裂隙网络渗流特性的重要因素。制作裂隙网络试验模型,建立高精度渗流试验系统;求解纳维-斯托克斯方程,模拟流体在裂隙网络内的流动状态,研究等效水力隙宽和水力梯度对非线性渗流特性的影响。结果表明,当水力梯度较小时,等效渗透系数保持恒定的常数,流体流动属于达西流动区域,流量与压力具有线性关系,可采用立方定律计算流体流动;当水力梯度较大时,等效渗透系数随着水力梯度的增加而急剧减少,流体流动进入强惯性效应流动区域,流量与压力具有强烈的非线性关系,可采用Forchheimer方程计算流体流动。随着等效水力隙宽的增加,区别线性和非线性流动区域的临界水力梯度呈幂函数关系递减。当水力梯度小于临界水力梯度时,控制方程可选立方定律;当水力梯度大于临界水力梯度时,控制方程可选Forchheimer方程,其参数A和B可根据经验公式计算得到。其研究结果可为临界水力梯度的确定及流体流动控制方程的选取提供借鉴意义。     

TTI介质qP波逆时偏移中伪横波噪声压制方法

在对地下复杂构造介质,特别是盐丘侧翼及岩下区域进行成像时,相对于传统的各向同性逆时偏移和VTI逆时偏移,具有倾斜对称轴的TTI逆时偏移成像效果最优.不仅反射同相轴更加的连续,而且能量得到了更好的聚焦.传统的各向异性介质全弹性波RTM的计算量大且计算效率低.由于目前仍以纵波勘探为主,因此TTI逆时偏移qP波波动方程的选取显得尤为重要.为了提高计算效率,采用将沿着对称轴方向的横波速度设为零的方法,简化得到qP波波动方程.然而,这样会引入一种严重影响成像效果的低速度、低振幅的qSV波人为干扰.本文建立了qP波方程的完全匹配层控制方程,而后借助于辅助波场采用一种高效的压制伪横波噪声传播的方法,通过模型测试验证了该方法的有效性.     

Dynamic model of mesoscale eddies. Eddy parameterization for coarse resolution ocean circulation models

In the framework of the eddy dynamic model developed in two previous papers (Dubovikov, M.S., Dynamical model of mesoscale eddies, Geophys. Astophys. Fluid Dyn., 2003, 97, 311–358; Canuto, V.M. and Dubovikov, M.S., Modeling mesoscale eddies, Ocean Modelling, 2004, 8, 1–30 referred as I–II), we compute the contribution of unresolved mesoscale eddies to the large-scale dynamic equations of the ocean. In isopycnal coordinates, in addition to the bolus velocity discussed in I–II, the mesoscale contribution to the large scale momentum equation is derived. Its form is quite different from the traditional down-gradient parameterization. The model solutions in isopycnal coordinates are transformed to level coordinates to parameterize the eddy contributions to the corresponding large scale density and momentum equations. In the former, the contributions due to the eddy induced velocity and to the residual density flux across mean isopycnals (so called Σ-term) are derived, both contributions being shown to be of the same order. As for the large scale momentum equation, as well as in isopycnal coordinates, the eddy contribution has a form which is quite different from the down-gradient expression.     

Nonlinear pulsations of horizontal jets

A shallow-water model with horizontally nonuniform density is used to study the dynamics of jet flows that arise under the influence of buoyancy and the Coriolis force. Within this approach, the jet is described by a self-similar compactly-localized solution and interpreted as a band of shear flow having a temperature contrast with the ambient fluid. In addition to stationary states, the dynamics of such jets admit cyclonic rotation with a constant angular velocity and transverse nonlinear pulsations. The phase portrait corresponding to this model shows that regimes with pulsating jets develop along closed trajectories bounded by the separatrix loop. The theory predicts that the period for warm jet pulsations is longer than the inertial oscillation period caused by the Earth’s rotation, while for cold jet pulsations, it is shorter. Thus, only warm jets can have a noticeable effect on the atmospheric dynamics in the synoptic range. In particular, they may well be responsible for additional spectral peaks that appear in this range of wind speed fluctuations.     

VTI介质中准P波方程叠前逆时深度偏移

根据具有垂直对称轴的横向各向同性（VTI）介质中的一阶准P波方程,导出了该方程在交错网格中逆时延拓的高阶有限差分格式,给出了其稳定性条件,采用完全匹配层吸收边界条件解决边界反射问题,分别应用下行波最大能量法和归一化互相关成像条件, 实现了VTI介质中准P波方程的叠前逆时深度偏移．各向异性Marmousi模型的试算结果表明,VTI介质准P波方程叠前逆时深度偏移算法不受地下构造倾角和介质横向速度变化的限制,对复杂模型具有良好的成像能力;应用归一化互相关成像条件能得到更好的成像效果．对比该模型的各向异性和各向同性逆时偏移剖面表明,在各向异性地区采集的纵波数据用各向异性偏移算法理论上能得到更好的成像结果．      

A wetting and drying algorithm with a combined pressure/free-surface formulation for non-hydrostatic models

A wetting and drying method for free-surface problems for the three-dimensional, non-hydrostatic Navier–Stokes equations is proposed. The key idea is to use a horizontally fixed mesh and to apply different boundary conditions on the free-surface in wet and dry zones. In wet areas a combined pressure/free-surface kinematic boundary condition is applied, while in dry areas a positive water level and a no-normal flow boundary condition are enforced. In addition, vertical mesh movement is performed to accurately represent the free-surface motion. Non-physical flow in the remaining thin layer in dry areas is naturally prevented if a Manning–Strickler bottom drag is used. The treatment of the wetting and drying processes applied through the boundary condition yields great flexibility to the discretisation used. Specifically, a fully unstructured mesh with any finite element choice and implicit time discretisation method can be applied. The resulting method is mass conservative, stable and accurate. It is implemented within Fluidity-ICOM [1] and verified against several idealized test cases and a laboratory experiment of the Okushiri tsunami.     

Divergence of solutions to perturbation-based advection-dispersion moment equations

Moment equation methods are popular and powerful tools for modeling transport processes in randomly heterogeneous porous media, but the application of these methods to advection-dispersion equations often leads to erroneous oscillations. Perturbative methods, required to close systems of moment equations, become inaccurate for large perturbations; however, little quantitative theory exists for determining when this occurs for advection-dispersion equations. We consider three different methods (asymptotic approximation, Eulerian truncation, and iterative solution) for closing and solving advection-dispersion moment equations describing transport in stratified porous media with random permeability. We obtain approximate analytical expressions for time above which the asymptotic approximation to the mean diverges, in particular quantifying the impact that dispersion has on delaying—but not eliminating—divergence. We demonstrate that Eulerian truncation and iterative solution methods do not eliminate divergent behavior either. Our divergence criteria provide a priori estimates that signal a warning to the practitioner of stochastic advection-dispersion equations to carefully consider whether to apply perturbative approaches.     

Discontinuous Galerkin methods for modeling Hurricane storm surge

Storm surge due to hurricanes and tropical storms can result in significant loss of life, property damage, and long-term damage to coastal ecosystems and landscapes. Computer modeling of storm surge can be used for two primary purposes: forecasting of surge as storms approach land for emergency planning and evacuation of coastal populations, and hindcasting of storms for determining risk, development of mitigation strategies, coastal restoration and sustainability.Storm surge is modeled using the shallow water equations, coupled with wind forcing and in some events, models of wave energy. In this paper, we will describe a depth-averaged (2D) model of circulation in spherical coordinates. Tides, riverine forcing, atmospheric pressure, bottom friction, the Coriolis effect and wind stress are all important for characterizing the inundation due to surge. The problem is inherently multi-scale, both in space and time. To model these problems accurately requires significant investments in acquiring high-fidelity input (bathymetry, bottom friction characteristics, land cover data, river flow rates, levees, raised roads and railways, etc.), accurate discretization of the computational domain using unstructured finite element meshes, and numerical methods capable of capturing highly advective flows, wetting and drying, and multi-scale features of the solution.The discontinuous Galerkin (DG) method appears to allow for many of the features necessary to accurately capture storm surge physics. The DG method was developed for modeling shocks and advection-dominated flows on unstructured finite element meshes. It easily allows for adaptivity in both mesh (h) and polynomial order (p) for capturing multi-scale spatial events. Mass conservative wetting and drying algorithms can be formulated within the DG method.In this paper, we will describe the application of the DG method to hurricane storm surge. We discuss the general formulation, and new features which have been added to the model to better capture surge in complex coastal environments. These features include modifications to the method to handle spherical coordinates and maintain still flows, improvements in the stability post-processing (i.e. slope-limiting), and the modeling of internal barriers for capturing overtopping of levees and other structures. We will focus on applications of the model to recent events in the Gulf of Mexico, including Hurricane Ike.     

Propagation of water waves over permeable rippled beds

Unsteady two-dimensional Navier-Stokes equations and Navier-Stokes type model equations for porous flow were solved numerically to simulate the propagation of water waves over a permeable rippled bed. A boundary-fitted coordinate system was adopted to make the computational meshes consistent with the rippled bed. The accuracy of the numerical scheme was confirmed by comparing the numerical results concerning the spatial distribution of wave amplitudes over impermeable and permeable rippled beds with the analytical solutions. For periodic incident waves, the flow field over the wavy wall is discussed in terms of the steady Eulerian streaming velocity. The trajectories of the fluid particles that are initially located close to the ripples were also determined. One of the main results herein is that under the action of periodic water waves, fluid particles on an impermeable rippled bed initially moved back and forth around the ripple crest, with increasing vertical distance from the rippled wall. After one or two wave periods, they are then lifted towards the next ripple crest. All of the marked particles on a permeable rippled bed were shifted onshore with a much larger displacement than those on an impermeable bed. Finally, the flow fields and the particle motions close to impermeable and permeable beds induced by a solitary wave are elucidated.