A primitive-variable Riemann method for solution of the shallow water equations with wetting and drying

A Riemann flux that uses primitive variables rather than conserved variables is developed for the shallow water equations with nonuniform bathymetry. This primitive-variable flux is both conservative and well behaved at zero depth. The unstructured finite-volume discretization used is suitable for highly nonuniform grids that provide resolution of complex geometries and localized flow structures. A source-term discretization is derived for nonuniform bottom that balances the discrete flux integral both for still water and in dry regions. This primitive-variable formulation is uniformly valid in wet and dry regions with embedded wetting and drying fronts. A fully nonlinear implicit scheme and both nonlinear and time-linearized explicit schemes are developed for the time integration. The implicit scheme is solved by a parallel Newton-iterative algorithm with numerically computed flux Jacobians. A concise treatment of characteristic-variable boundary conditions with source terms is also given. Computed results obtained for the one-dimensional dam break on wet and dry beds and for normal-mode oscillations in a circular parabolic basin are in very close agreement with the analytical solutions. Other results for a forced breaking wave with friction interacting with a sloped bottom demonstrate a complex wave motion with wetting, drying and multiple interacting wave fronts. Finally, a highly nonuniform, coastline-conforming unstructured grid is used to demonstrate an unsteady simulation that models an artificial coastal flooding due to a forced wave entering the Gulf of Mexico.     

HL-2M偏滤器CFC/CuCrZr靶板的非晶钎焊工艺研究

在采用商用的STEMET-1101非晶钎料、温度为710~750°C的真空钎焊下,对CFC/OFC(氯氟烃/无氧铜)复合块与CuCrZr(铬锆铜)的钎焊进行了研究。首先,通过X-射线衍射(XRD)和差式扫描量热分析对钎料的结构和熔化行为进行了表征;然后,通过光学显微镜、电子探针微分析和拉伸试验等方法对焊缝的组织形貌、元素成分分布、相结构和力学性能进行了分析;最后,通过高热负荷装置对CFC/OFC/CuCrZr钎焊模块的热疲劳性能进行测试。结果表明,在710~750°C钎焊温度内焊缝由Cu固溶体、(Cu, Ni)3P和Ni(Cu ,Cr)2P金属间化合物组成,焊缝平整无裂纹;特别是在750°C/15min情况下,抑制了焊缝金属间化合物的连续分布,OFC/CuCrZr的焊接强度大于OFC的抗拉强度,CuCrZr/CuCrZr的结合强度为210MPa,并呈现部分韧性断裂。在750°C/15min情况下制备的CFC/Cu/CuCrZr模块可以承受1000次7MW·m^-2的循环热负荷。     

A high-order discontinuous Galerkin method for wave propagation through coupled elastic–acoustic media

We introduce a high-order discontinuous Galerkin (dG) scheme for the numerical solution of three-dimensional (3D) wave propagation problems in coupled elastic–acoustic media. A velocity–strain formulation is used, which allows for the solution of the acoustic and elastic wave equations within the same unified framework. Careful attention is directed at the derivation of a numerical flux that preserves high-order accuracy in the presence of material discontinuities, including elastic–acoustic interfaces. Explicit expressions for the 3D upwind numerical flux, derived as an exact solution for the relevant Riemann problem, are provided. The method supports h-non-conforming meshes, which are particularly effective at allowing local adaptation of the mesh size to resolve strong contrasts in the local wavelength, as well as dynamic adaptivity to track solution features. The use of high-order elements controls numerical dispersion, enabling propagation over many wave periods. We prove consistency and stability of the proposed dG scheme. To study the numerical accuracy and convergence of the proposed method, we compare against analytical solutions for wave propagation problems with interfaces, including Rayleigh, Lamb, Scholte, and Stoneley waves as well as plane waves impinging on an elastic–acoustic interface. Spectral rates of convergence are demonstrated for these problems, which include a non-conforming mesh case. Finally, we present scalability results for a parallel implementation of the proposed high-order dG scheme for large-scale seismic wave propagation in a simplified earth model, demonstrating high parallel efficiency for strong scaling to the full size of the Jaguar Cray XT5 supercomputer.     

Enhanced pinning of superconducting vortices at circular magnetic dots in the magnetic-vortex state

Low-temperature magnetic force microscopy (LT-MFM) was used to study the distribution of superconducting vortices in Nb above a square array of 1 μm-sized circular ferromagnetic dots in a magnetic-vortex state. The force that the MFM tip exerts on the individual vortex in the depinning process was used to estimate the spatial modulation of the pinning potential. It was found, that the superconducting vortices which are preferably located on top of the Py dots experience a pinning force, about 15 times stronger as compared to the pinning force in the pure Nb film. This strong pinning exceeds the repulsive interaction between the superconducting vortices and allows vortex clusters to be located above the dots.     

Probing the discrete motion of vortices with rf excitations

In this work we present experimental results on the rectification of vortices in a superconductor/ferromagnet system under a high frequency drive. The two-dimensional pinning landscape, induced by the stray fields of the ferromagnetic template, has no intrinsic asymmetry. Nevertheless, an asymmetric potential is artificially induced by an applied dc bias. The experimental results unambiguously show a biased, discrete motion of the vortices in the periodic potential at frequencies above 10 MHz. This synchronized motion is very sensitive to the external applied field. Increasing temperature leads to a reduction of the pinning potential, which in turn results in a lower ac power needed to drive the vortex lattice.     

WALL CONDUCTION EFFECT ON THE ONSET OF THERMAL INSTABILITY IN HORIZONTAL RECTANGULAR CHANNELS

The effect of peripheral thermal conduction in the bottom wall on the onset of thermal instability in the thermal entrance region of horizontal, rectangular, bottom-uniform-heated channels was studied by experiment and theoretical simulation. The predicted results of the onset of thermal instability for the cases of uniform peripheral wall heat flux are in agreement with the published experimental data, but the predicted results for the cases of uniform peripheral wall temperature exceed the previous data by about an order of magnitude. The present experimental data reconfirm the present predictions. The effect of increasing the peripheral wall conduction is shown to stabilize the flow.     

Combustion mechanism of ultralean rotating counterflow twin premixed flame

In our previous numerical studies [Nishioka Makihito, Zhenyu Shen, and Akane Uemichi. “Ultra-lean combustion through the backflow of burned gas in rotating counterflow twin premixed flames.” Combustion and Flame 158.11 (2011): 2188–2198. Uemichi Akane, and Makihito Nishioka. “Numerical study on ultra-lean rotating counterflow twin premixed flame of hydrogen–air.” Proceedings of the Combustion Institute 34.1 (2013): 1135–1142]. we found that methane– and hydrogen–air rotating counterflow twin flames (RCTF) can achieve ultralean combustion when backward flow of burned gas occurs due to the centrifugal force created by rotation. In this study, we investigated the mechanisms of ultralean combustion in these flames by the detailed numerical analyses of the convective and diffusive transport of the main species. We found that, under ultralean conditions, the diffusive transport of fuel exceeds its backward convective transport in the flame zone, which is located on the burned-gas side of the stagnation point. In contrast, the relative magnitudes of diffusive and convective transport for oxygen are reversed compared to those for the fuel. The resulting flows for fuel and oxygen lead to what we call a ‘net flux imbalance’. This net flux imbalance increases the flame temperature and concentrations of active radicals. For hydrogen–air RCTF, a very large diffusivity of hydrogen enhances the net flux imbalance, significantly increasing the flame temperature. This behaviour is intrinsic to a very lean premixed flame in which the reaction zone is located in the backflow of its own burned gas.     

Numerical solution for curved crack problem in elastic half-plane using hypersingular integral equation

A hypersingular integral equation for the curved crack problems of an elastic half-plane is introduced. Formulation of the equation is based on the usage of a modified complex potential. The potential is generally expressed in the form of a Cauchy-type integral. The modified complex potential is composed of the principal part and the complementary part. The principal part of the complex potential is actually equivalent to the original complex potential for the curved crack in an infinite plate. The role of the complementary part is to eliminate the boundary traction along the boundary of the half-plane caused by the principal part. From the assumed boundary traction condition, a hypersingular integral equation is obtained for the curved crack problems of an elastic half-plane. The curve length coordinate method is used to obtain a final solution. Several numerical examples are presented that prove the efficiency of the suggested method.     

Microstructure evolution during severe plastic deformation

Radiotracer diffusion studies of severely deformed, ultra-fine grained materials have revealed the presence of ultra-fast transport paths, which include “non-equilibrium” grain boundaries and free volume. Under some experimental conditions, percolating porosity is produced even in pure copper. Micro-cracks may form in metals, if the local maximum shear stress exceeds the shear yield stress. However, their growth and propagation is postponed till late in the deformation process owing to the ductility of metals, the hydrostatic component of the stress system and/or dynamic recovery/recrystallization. In other words, crack growth and propagation is present only when the scope for further deformation is highly restricted. Using this approach, the load required for equal channel angular pressing, the change in the slope of the Hall–Petch plot with decreasing grain size and the theoretical limit for the smallest grain size attainable in a metal in a severe plastic deformation process are predicted and validated by experimental results. Experimentally successful prevention of percolated crack formation by the superposition of a hydrostatic pressure is also accounted for using this model.