Computation of flow over a rotating body on unstructured Chimera mesh

Flow around moving boundary is ubiquitous in engineering applications. To increse the efficienly of the algorithm to handle moving boundaries is still a major challenge in Computational Fluid Dynamics (CFD). The Chimera grid method is one type of method to handle moving boundaries. A concept of domain de-composition has been proposed in this paper. In this method, sub-domains are meshed independently and governing equations are also solved separately on them. The Chimera grid method was originally used only on structured (curvilinear) meshes. However, in a problem which involves both moving boundary and complex geometry, the number of sub-domains required in a traditional (structured) Chimera method becomes fairly large. Thus the time required in the interior boundary locating, link-building and data exchanging also increases. The use of unstructured Chimera grid can reduce the time consumption significantly by the reduction of domain(block) number. Generally speaking, unstructured Chimera grid method has not been developed. In this paper, a well-known pressure correction scheme - SIMPLEC is modified and implemented on unstructured Chimera mesh. A new interpolation scheme regarding the pressure correction is proposed to prevent the possible decoupling of pressure. A moving-mesh finite volume approach is implemented in an inertial reference frame. This approach is then used to compute incompressible flow around a rotating circular and elliptic cylinder. These numerical examples demonstrate the capability of the proposed scheme in handling moving boundaries. The numerical results are in good agreement with other experimental and computational data in literature. The method proposed in this paper can be efficiently applied to more challenge cases such as free-falling objects or heavy particles in fluid.

Space-time correlations of fluctuating velocities in turbulent shear flows

Space-time correlations or Eulerian two-point two-time correlations of fluctuating velocities are analytically and numerically investigated in turbulent shear flows. An elliptic model for the space-time correlations in the inertial range is developed from the similarity assumptions on the isocorrelation contours: they share a uniform preference direction and a constant aspect ratio. The similarity assumptions are justified using the Kolmogorov similarity hypotheses and verified using the direct numerical simulation DNS of turbulent channel flows. The model relates the space-time correlations to the space correlations via the convection and sweeping characteristic velocities. The analytical expressions for the convection and sweeping velocities are derived from the Navier-Stokes equations for homogeneous turbulent shear flows, where the convection velocity is represented by the mean velocity and the sweeping velocity is the sum of the random sweeping velocity and the shearinduced velocity. This suggests that unlike Taylor’s model where the convection velocity is dominating and Kraichnan and Tennekes’ model where the random sweeping velocity is dominating, the decorrelation time scales of the space-time correlations in turbulent shear flows are determined by the convection velocity, the random sweeping velocity, and the shear-induced velocity. This model predicts a universal form of the spacetime correlations with the two characteristic velocities. The DNS of turbulent channel flows supports the prediction: the correlation functions exhibit a fair good collapse, when plotted against the normalized space and time separations defined by the elliptic model.

Dimensionless Analysis of the Simulation of Bucket Foundations under Dynamic Load

Firstly, the main factors are obtained by use of dimensionless analysis. Secondly, the time scaling factors in centrifuge modeling of bucket foundations under dynamic load are analyzed based on dimensionless analysis and control- ling equation. A simplified method for dealing with the conflict of scaling factors of the inertial and the percolation in sand foundation is presented. The presented method is that the material for experiments is not changed while the effects are modified by perturbation method. Thirdly, the characteristic time of liquefaction state and the characteristic scale of affected zone are analyzed.

Angular distribution and conversion of multi-keV L-shell X-ray sources produced from nanosecond laser irradiated thick-foil targets

An experimental study on the angular distribution and conversion of multi-keV X-ray sources produced from 2 ns-duration 527nm laser irradiated thick-foil targets on Shenguang II laser facility (SG-II) is reported. The angular distributions measured in front of the targets can be fitted with the function of f(theta) = alpha+ (1- alpha)cos(beta) theta (theta is the viewing angle relative to the target normal), where alpha = 0.41 +/- 0.014, beta = 0.77 +/- 0.04 for Ti K-shell X-ray Sources (similar to 4.75 keV for Ti K-shell), and alpha = 0.085 +/- 0.06, beta = 0.59 +/- 0.07 for Ag/Pd/Mo L-shell X-ray Sources (2-2.8 keV for Mo L-shell, 2.8-3.5 keV for Pd L-shell, and 3-3.8 keV for Ag L-shell). The isotropy of the angular-distribution of L-shell emission is worse than that of the K-shell emission at larger viewing angle (>70 degrees), due to its larger optical depth (stronger self-absorption) in the cold plasma side lobe Surrounding the central emission region, and in the central hot plasma region (emission region). There is no observable difference in the angular distributions of the L-shell X-ray emission among Ag, Pd, and Mo. The conversion efficiency of Ag/Pd/Mo L-shell X-ray sources is higher than that of the Ti K-shell X-ray sources, but the gain relative to the K-shell emission is not as high as that by using short pulse lasers. The conversion efficiency of the L-shell X-ray sources decrease, with increasing atomic numbers (or X-ray photon energy), similar to the behavior of the K-shell X-ray Source.

Direct acceleration of solid-density plasma bunch by ultraintense laser

The interaction of a petawatt laser with a small solid-density plasma bunch is studied by particle-in-cell simulation. It is shown that when irradiated by a laser of intensity >10(21) W/cm(2), a dense plasma bunch of micrometer size can be efficiently accelerated. The kinetic energy of the ions in the high-density region of the plasma bunch can exceed ten MeV at a density in the 10(23)-cm(-3) level. Having a flux density orders of magnitude higher than that of the traditional charged-particle pulses, the laser-accelerated plasma bunch can have a wide range of applications. In particular, such a dense energetic plasma bunch impinging on the compressed fuel in inertial fusion can significantly enhance the nuclear-reaction cross section and is thus a promising alternative for fast ignition.

Intense local plasma heating by stopping of ultrashort ultraintense laser pulse in dense plasma

Self-trapping, stopping, and absorption of an ultrashort ultraintense linearly polarized laser pulse in a finite plasma slab of near-critical density is investigated by particle-in-cell simulation. As in the underdense plasma, an electron cavity is created by the pressure of the transmitted part of the light pulse and it traps the latter. Since the background plasma is at near-critical density, no wake plasma oscillation is created. The propagating self-trapped light rapidly comes to a stop inside the slab. Subsequent ion Coulomb explosion of the stopped cavity leads to explosive expulsion of its ions and formation of an extended channel having extremely low plasma density. The energetic Coulomb-exploded ions form shock layers of high density and temperature at the channel boundary. In contrast to a propagating pulse in a lower density plasma, here the energy of the trapped light is deposited onto a stationary and highly localized region of the plasma. This highly localized energy-deposition process can be relevant to the fast ignition scheme of inertial fusion.

谱色散均匀化的计算研究

利用光传输理论对ICF驱动器中使用的光谱色散平滑(SSD)技术作了理论分析，并结合典型的随机相位板(RPP)技术，用计算机模拟了使用光谱色散平滑技术前后激光靶面辐照不均匀性的变化及其技术中不同带宽和调制频率下激光靶面辐照不均匀性的变化。

图像处理在光路自动准直系统中的应用

光路自动准直系统应用于惯性约束聚变的高功率激光装置中的光束自动调整。图像处理是光路自动准直的关键技术之一。针对神光Ⅲ原型装置，结合阈值化、重心法、中值滤波和圆拟合等多种不同的图像处理方法设计了一套合理的准直方案，并且在模拟实验平台上进行了实验验证。实验结果表明，光路自动准直系统能够在15min之内顺利完成光路的自动调整，光束近场调整精度优于近场光斑的±0．5％，光束远场调整精度≤±0．3″，满足了原型装置的总体要求。

利用消衍射透镜列阵及光谱色散平滑实现焦斑均匀辐照

提出将空间域的透镜列阵法和时间域的光谱色散平滑法结合起来实现靶面的均匀辐照．消衍射型透镜列阵能获得边缘陡峭且顶部较平坦的准近场焦斑，光谱色散平滑则能有效地抹平焦斑内部由多光束干涉引起的细密条纹．数值结果显示，通过该方案能获得均匀性较好的焦斑．进一步分析了光谱色散平滑单元中位相调制和光栅的参数对辐照均匀性的影响，发现参数的选取要在焦斑均匀性和能量利用率之间取得合理平衡，以在整体上获得最佳的均匀辐照效果．

谱色散平滑与透镜列阵联用实现均匀照明

为满足激光惯性约束聚变中靶面激光辐照不均匀性低于5%的要求, 在目前使用透镜列阵基础上, 提出了谱色散平滑与透镜列阵联用方案, 对其进行数值计算并分析其平滑效果和应用可行性。结果表明:焦斑的不均匀性从单独使用透镜列阵时的14%降低到与谱色散平滑结合后的3%;对焦斑点功率谱的分析表明谱色散平滑通过抑制焦斑中高频的频谱强度达到平滑效果。该方案可以进一步提高焦斑平滑效果, 计算结果对实际应用有着重要的参考意义。

高功率激光装置前端系统激光脉冲的时间整形及修正技术

针对激光惯性约束聚变实验研究对高功率激光驱动器前端系统复杂时间形状种子激光脉冲的需求, 应用孔径耦合带状线集成波导整形系统设计了满足需要的前端整形激光脉冲。用一种新方法精确计算了孔径耦合带状线电脉冲整形器的耦合系数和孔径宽度的数值关系, 并针对高衬比度整形激光脉冲的需求, 提出了高衬比度双极型集成波导整形系统方案。由该系统可以得到100 ps脉冲前沿、1～3 ns脉冲宽度可调、高衬比度(大于100∶1)、光滑无纹波调制、可精确满足神光II八路及第九路装置需求的前端整形激光脉冲。

ICF中可调法布里-珀罗滤波器间距稳定度研究

在神光Ⅱ第9路ICF高功率激光装置中，采用可调法布里-珀罗（F-P）滤波器对幅度调制效应进行补偿，根据补偿装置的技术要求，提出-种应用nm量级精度的电容式位移传感器对可调F-P滤波器间距稳定度进行监控的系统，详细论述了监控系统的结构与工作原理。给出了电容式位移传感器的驱动电路及数据处理与控制软件的设计方案，并对电容式位移传感器的精度进行了标定。实验结果表明，该位移监控系统能够使可调F—P滤波器的间距稳定度保持在15nm／h以内，使幅度调制效应的调制深度优于4％。

大口径八边形钕玻璃片支撑系统的光机集成分析

控制由机械装夹方式所引入的波前畸变以提高惯性约束聚变(ICF)输出光束的质量，是在大口径钕玻璃片主放大器结构设计中必须考虑的。提出了一种新的有限元变形结果与光学元件面形畸变之间的数据处理方式，并与传统方式进行了对比。基于新的数据处理接口，利用光机集成分析方法对大口径八边形钕玻璃片的支撑系统结构设计参数进行优化。优化的结果保证了由支撑系统引起的透过波前畸变小于十分之一波长，同时波前畸变与设计参数变动的相关性最小。

二元振幅型面板用于光束空间整形

为了提高高功率激光系统的整体效率和充分利用光能,需要对前端注入的高斯光束进行空间整形,实现驱动器终端激光的均匀化输出。采用振幅型二元面板对激光光束进行空间强度整形,利用误差扩散法进行了理论设计,数值摸拟了整形效果,同时讨论了面板加工误差以及空间滤波器的小孔大小等因素带来的影响。根据理论设计,分别加工了反高斯透射率分布和抛物线透射率分布的二元面板,并进行了整形实验,实现了各自的整形功能,并做了误差分析。实验证明二元面板能对激光光束的空间强度分布实现了精确的整形。

一种ICF用新型可阵列反射镜架结构

运用并联机构原理设计了一种新型惯性约束聚变（Inetial Confinement Fusion, ICF）用可阵列反射镜架结构,并用坐标变换方法分析了机构的调整正交性与旋转轴正交性之间的关系.装置的实际测试结果表明调整正交性与分析结果一致,并且反射镜架的稳定性以及调节准确度均能符合ICF装置的使用要求.