Hierarchically penalized additive hazards model with diverging number of parameters

In many applications,covariates can be naturally grouped.For example,for gene expression data analysis,genes belonging to the same pathway might be viewed as a group.This paper studies variable selection problem for censored survival data in the additive hazards model when covariates are grouped.A hierarchical regularization method is proposed to simultaneously estimate parameters and select important variables at both the group level and the within-group level.For the situations in which the number of parameters tends to∞as the sample size increases,we establish an oracle property and asymptotic normality property of the proposed estimators.Numerical results indicate that the hierarchically penalized method performs better than some existing methods such as lasso,smoothly clipped absolute deviation(SCAD)and adaptive lasso.     

On the case when steady converging/diverging flow of a non-Newtonian fluid in a round cone permits an exact solution

This communication considers the steady converging/diverging flow of a non-Newtonian viscous power-law fluid in a round cone. The motion is driven by a sink/source of mass at the origin. It is shown that the problem permits exact similarity solution for a particular value (n=4/3) of the fluid index. In this case a complete set of governing equations can be reduced to an ordinary differential equation, which is solved numerically for different values of the main non-dimensional parameters (the cone angle and the dimensionless sink/source intensity).     

应用阵面追踪法对散心爆轰波传播的数值模拟

考虑到爆轰波阵面曲率及化学反应区宽度的作用 ,由修正Hugoniot关系式解析求解了波后状态 ,应用阵面追踪法 (FTM)数值模拟了曲面散心爆轰波的传播。计算给出的散心爆轰波走时及波阵面上物理量都达到了波阵面曲率的一阶精度。     

Experimental and theoretical study of short fiber orientation in diverging flows

Two types of experiments have been carried out to study the fiber orientation in flow through a divergent channel. First, a reinforced polyamid mold sprue containing two types of orientation was investigated: near the center, the fibers are mostly oriented perpendicular to the flow lines, whereas on the periphery, they are oriented parallel to them. Second, direct observation of copper fibers moving in a corn syrup was performed in a transparent diverging device: the fibers rapidly become oriented transverse to the flow lines. The solution of Stokes equations for the undisturbed fluid motion gives the shear rate and elongation rate, which are then substituted in Jeffery's orientation equations. The resolution shows two types of behavior: in a large area in the center, the fiber tends to a stable equilibrium position which depends strongly on the flow line on which it moves. On the periphery, the fiber follows a shear-like behavior. The strong influence of the elongational component relative to the shear component is demonstrated and the time necessary for orientation is calculated. The theoretical results are found to be in agreement with the observations.     

Axial acoustic radiation force of progressive cylindrical diverging waves on a rigid and a soft cylinder immersed in an ideal compressible fluid

Background and motivation

Previous works investigating the radiation force of diverging spherical progressive waves incident upon spherical particles have demonstrated the direction of reversal of the force when the particle is subjected to a curved wave-front. In this communication, the analysis is extended to the case of diverging cylindrical progressive waves incident upon a rigid or a soft cylinder in a non-viscous fluid with explicit calculations for the radiation force function (which is the radiation force per unit energy density and unit cross-sectional surface) not shown in [F.G. Mitri, Ultrasonics 50 (2010) 620-627].

Method

A closed-form solution presented previously in [F.G. Mitri, Ultrasonics 50 (2010) 620-627] is used to plot the radiation force function with particular emphasis on the difference from the results of incident plane progressive waves versus the size parameter ka (k is the wave number and a is the cylinder’s radius) and the distance of the cylinder from the acoustic source r₀.

Results

Radiation force function calculations for the rigid cylinder unexpectedly reveal that under specific conditions determined by the frequency of the acoustic field, the radius of the cylinder, as well as the distance to the acoustic source, the force becomes attractive (negative force). In addition, the numerical results show that the radiation force on a rigid cylinder does not generally obey the inverse-distance law with respect to the distance from the source.

Conclusion and potential applications

These results suggest that it may be possible, under specific conditions, to pull a cylindrical structure back toward the acoustic source using progressive cylindrical diverging waves. They may also provide a means to predict the radiation force required to manipulate non-destructively a single cylindrical structure. Potential applications include the design of a new generation of acoustic tweezers operating using a single beam of progressive waves (in contrast to the traditional version of acoustical tweezers in which an acoustic standing wave field is produced using two counter-propagating acoustic fields) for investigations in the field of flow cytometry, particle manipulation and entrapment.     

宏观均相多孔电极电化学阻抗谱基础

电化学阻抗谱可用于诊断多孔电极内电荷转移反应,即界面电荷集聚和电荷传导,以及反应物质输运。本文采用复相量方法,在同态假设条件下,重新推演多孔电极阻抗谱模型,厘清传统多孔电极阻抗谱模型中的模糊性表述。(1) 定义多孔电极表征输入参数,包括电极基体电子电导率σ₁ 、电解质离子电导率σ₂、界面电荷传递电导率g_ct、单位面积界面电容C、固相扩散系数D、速度常数k、电极厚度d、特征孔深L_p 和单位体积表面积S_c;(2) 解析阻抗谱特征输出参数,包括场扩散常数K,特征频率ω₀、ω₁、ω₂、ω₃和 ω_max,它们分别相关于界面传导反应、有限场扩散、氧化还原反应、孔内扩散和最小特征孔尺寸,以及分别对应于从传导到扩散和从扩散到饱和的转折频率f_k1 和f_k2;(3) 当参数X和Z同时变化时(X = σ₁和Z = d,S_c,L_p,C,g_ct,D,k),通过阻抗谱特征参数的演变规律,分析了电荷转移反应中X和Ζ参数耦合竞争;(4)为深入分析电荷转移反应中参数X和Z的耦合竞争,引入了分叉频率f_XZ和f_ZX 。f_XZ和f_ZX所处位置可以用于表征参数X和Z影响电荷转移反应的深度和广度。当分叉频率f_XZ和f_ZX不存在时,表明电荷转移反应中参数X和Z在全频率范围内存在耦合竞争。总之,借助于特征频率和分叉频率,本文一方面研究了动力学参数和微观结构参数对多孔电极中电荷转移反应的影响,另一方面分析谱图的变化及其背后的阻抗谱特征演化规律。本文研究结果可为阻抗谱的系统仿真和辨识提供理论基础,可为多孔电极内电荷转移反应的竞争分析提供技术支撑,还可为电化学储能系统的优化设计提供诊断工具。     

In and out forests on combinatorial landscapes

Fitness landscape theory is a mathematical framework for numerical analysis of search algorithms on combinatorial optimization problems. We study a representation of fitness landscape as a weighted directed graph. We consider out forest and in forest structures in this graph and establish important relationships among the forest structures of a directed graph, the spectral properties of the Laplacian matrices, and the numbers of local optima of the landscape. These relationships provide a new approach for computing the numbers of local optima for various problem instances and neighborhood structures.     

SICA for Cox’s proportional hazards model with a diverging number of parameters

The smooth integration of counting and absolute deviation （SICA） penalized variable selection procedure for high-dimensional linear regression models is proposed by Lv and Fan （2009）. In this article, we extend their idea to Cox＇s proportional hazards （PH） model by using a penalized log partial likelihood with the SICA penalty. The number of the regression coefficients is allowed to grow with the sample size. Based on an approximation to the inverse of the Hessian matrix, the proposed method can be easily carried out with the smoothing quasi-Newton （SQN） algorithm. Under appropriate sparsity conditions, we show that the resulting estimator of the regression coefficients possesses the oracle property. We perform an extensive simulation study to compare our approach with other methods and illustrate it on a well known PBC data for predicting survival from risk factors.     

Numerical simulation of tangentially injected turbulent swirling flow in a divergent tube

This paper studies the properties of turbulent swirling decaying flow induced by tangential inlets in a divergent pipe using the realizable k–ε turbulence model and discusses the effects of the injector pressure and injection position. The results of transient solutions show that both the recirculation zone near the wall in upstream of the injectors and the vortex breakdown in downstream of the injectors increase in size during the whole period. A nearly axisymmetric conical breakdown is formed and its internal structure consists of two asymmetric spiral‐like vortices rotating in opposite directions. The stagnation point shifts slowly toward the pipe outlet over time. The maxima of the three velocity components, which are located near the wall, decrease gradually with streamwise direction. It can also be inferred that Mach number approaches 1.0 near the injector outlets. The velocities increase with the increasing injector pressure. However, its increasing trend is not significant. With the increase of the injection position, vortex breakdown moves in downstream direction and the pitch along the streamwise direction increases. Copyright © 2008 John Wiley & Sons, Ltd.     

测量发散波面曲率半径的全息光栅－透镜成像法

介绍一种用全息光栅－透镜成像法精确测量发散球面波、柱面波曲率半径的新方法，分析测量误差，指出提高测量精度的途径。