scholarly journals Magnetohydrodynamic simulation code CANS+: Assessments and applications

Author(s):  
Yosuke Matsumoto ◽  
Yuta Asahina ◽  
Yuki Kudoh ◽  
Tomohisa Kawashima ◽  
Jin Matsumoto ◽  
...  

Abstract We present a new magnetohydrodynamic (MHD) simulation package with the aim of providing accurate numerical solutions to astrophysical phenomena where discontinuities, shock waves, and turbulence are inherently important. The code implements the Harten–Lax–van Leer–discontinuitues (HLLD) approximate Riemann solver, the fifth-order-monotonicity-preserving interpolation (MP5) scheme, and the hyperbolic divergence cleaning method for a magnetic field. This choice of schemes has significantly improved numerical accuracy and stability, and saved computational costs in multidimensional problems. Numerical tests of one- and two-dimensional problems show the advantages of using the high-order scheme by comparing with results from a standard second-order total variation diminishing monotonic upwind scheme for conservation laws (MUSCL) scheme. The present code enables us to explore the long-term evolution of a three-dimensional accretion disk around a black hole, in which compressible MHD turbulence causes continuous mass accretion via nonlinear growth of the magneto-rotational instability (MRI). Numerical tests with various computational cell sizes exhibits a convergent picture of the early nonlinear growth of the MRI in a global model, and indicates that the MP5 scheme has more than twice the resolution of the MUSCL scheme in practical applications.

Author(s):  
Tianjiao Wang ◽  
Jun Zhao ◽  
Chuanxin Weng ◽  
Tong Wang ◽  
Yayun Liu ◽  
...  

Shape memory polymers (SMPs) that change shapes as designed by external stimuli have become one of the most promising materials as actuators, sensors, and deployable devices. However, their practical applications...


Author(s):  
Sheng Yu-ming ◽  
Li Chao ◽  
Xia Ming-yao ◽  
Zou Jin-feng

Abstract In this study, elastoplastic model for the surrounding rock of axisymmetric circular tunnel is investigated under three-dimensional (3D) principal stress states. Novel numerical solutions for strain-softening surrounding rock were first proposed based on the modified 3D Hoek–Brown criterion and the associated flow rule. Under a 3D axisymmetric coordinate system, the distributions for stresses and displacement can be effectively determined on the basis of the redeveloped stress increment approach. The modified 3D Hoek–Brown strength criterion is also embedded into finite element software to characterize the yielding state of surrounding rock based on the modified yield surface and stress renewal algorithm. The Euler implicit constitutive integral algorithm and the consistent tangent stiffness matrix are reconstructed in terms of the 3D Hoek–Brown strength criterion. Therefore, the numerical solutions and finite element method (FEM) models for the deep buried tunnel under 3D principal stress condition are presented, so that the stability analysis of surrounding rock can be conducted in a direct and convenient way. The reliability of the proposed solutions was verified by comparison of the principal stresses obtained by the developed numerical approach and FEM model. From a practical point of view, the proposed approach can also be applied for the determination of ground response curve of the tunnel, which shows a satisfying accuracy compared with the measuring data.


Author(s):  
Athanasios Donas ◽  
Ioannis Famelis ◽  
Peter C Chu ◽  
George Galanis

The aim of this paper is to present an application of high-order numerical analysis methods to a simulation system that models the movement of a cylindrical-shaped object (mine, projectile, etc.) in a marine environment and in general in fluids with important applications in Naval operations. More specifically, an alternative methodology is proposed for the dynamics of the Navy’s three-dimensional mine impact burial prediction model, Impact35/vortex, based on the Dormand–Prince Runge–Kutta fifth-order and the singly diagonally implicit Runge–Kutta fifth-order methods. The main aim is to improve the time efficiency of the system, while keeping the deviation levels of the final results, derived from the standard and the proposed methodology, low.


2019 ◽  
Vol 85 (4) ◽  
Author(s):  
Benjamin D. G. Chandran ◽  
Jean C. Perez

We present three-dimensional direct numerical simulations and an analytic model of reflection-driven magnetohydrodynamic (MHD) turbulence in the solar wind. Our simulations describe transverse, non-compressive MHD fluctuations within a narrow magnetic flux tube that extends from the photosphere, through the chromosphere and corona and out to a heliocentric distance  $r$ of 21 solar radii  $(R_{\odot })$ . We launch outward-propagating ‘ $\boldsymbol{z}^{+}$ fluctuations’ into the simulation domain by imposing a randomly evolving photospheric velocity field. As these fluctuations propagate away from the Sun, they undergo partial reflection, producing inward-propagating ‘ $\boldsymbol{z}^{-}$ fluctuations’. Counter-propagating fluctuations subsequently interact, causing fluctuation energy to cascade to small scales and dissipate. Our analytic model incorporates dynamic alignment, allows for strongly or weakly turbulent nonlinear interactions and divides the $\boldsymbol{z}^{+}$ fluctuations into two populations with different characteristic radial correlation lengths. The inertial-range power spectra of $\boldsymbol{z}^{+}$ and $\boldsymbol{z}^{-}$ fluctuations in our simulations evolve toward a $k_{\bot }^{-3/2}$ scaling at $r>10R_{\odot }$ , where $k_{\bot }$ is the wave-vector component perpendicular to the background magnetic field. In two of our simulations, the $\boldsymbol{z}^{+}$ power spectra are much flatter between the coronal base and $r\simeq 4R_{\odot }$ . We argue that these spectral scalings are caused by: (i) high-pass filtering in the upper chromosphere; (ii) the anomalous coherence of inertial-range $\boldsymbol{z}^{-}$ fluctuations in a reference frame propagating outwards with the $\boldsymbol{z}^{+}$ fluctuations; and (iii) the change in the sign of the radial derivative of the Alfvén speed at $r=r_{\text{m}}\simeq 1.7R_{\odot }$ , which disrupts this anomalous coherence between $r=r_{\text{m}}$ and $r\simeq 2r_{\text{m}}$ . At $r>1.3R_{\odot }$ , the turbulent heating rate in our simulations is comparable to the turbulent heating rate in a previously developed solar-wind model that agreed with a number of observational constraints, consistent with the hypothesis that MHD turbulence accounts for much of the heating of the fast solar wind.


1998 ◽  
Vol 4 (S2) ◽  
pp. 378-379
Author(s):  
Z. W. Chen ◽  
D. B. Wittry

A monochromatic x-ray microprobe based on a laboratory source has recently been developed in our laboratory and used for fluorescence excitation. This technique provides high sensitivity (ppm to ppb), nondestructive, quantitative microanalysis with minimum sample preparation and does not require a high vacuum specimen chamber. It is expected that this technique (MMXRF) will have important applications in materials science, geological sciences and biological science.Three-dimensional focusing of x-rays can be obtained by using diffraction from doubly curved crystals. In our MMXRF setup, a small x-ray source was produced by the bombardment of a selected target with a focused electron beam and a toroidal mica diffractor with Johann pointfocusing geometry was used to focus characteristic x-rays from the source. In the previous work ∼ 108 photons/s were obtained in a Cu Kα probe of 75 μm × 43 μm in the specimen plane using the fifth order reflection of the (002) planes of mica.


Author(s):  
Khaled E. Zaazaa ◽  
Brian Whitten ◽  
Brian Marquis ◽  
Erik Curtis ◽  
Magdy El-Sibaie ◽  
...  

Accurate prediction of railroad vehicle performance requires detailed formulations of wheel-rail contact models. In the past, most dynamic simulation tools used an offline wheel-rail contact element based on look-up tables that are used by the main simulation solver. Nowadays, the use of an online nonlinear three-dimensional wheel-rail contact element is necessary in order to accurately predict the dynamic performance of high speed trains. Recently, the Federal Railroad Administration, Office of Research and Development has sponsored a project to develop a general multibody simulation code that uses an online nonlinear three-dimensional wheel-rail contact element to predict the contact forces between wheel and rail. In this paper, several nonlinear wheel-rail contact formulations are presented, each using the online three-dimensional approach. The methods presented are divided into two contact approaches. In the first Constraint Approach, the wheel is assumed to remain in contact with the rail. In this approach, the normal contact forces are determined by using the technique of Lagrange multipliers. In the second Elastic Approach, wheel/rail separation and penetration are allowed, and the normal contact forces are determined by using Hertz’s Theory. The advantages and disadvantages of each method are presented in this paper. In addition, this paper discusses future developments and improvements for the multibody system code. Some of these improvements are currently being implemented by the University of Illinois at Chicago (UIC). In the accompanying “Part 2” and “Part 3” to this paper, numerical examples are presented in order to demonstrate the results obtained from this research.


2017 ◽  
Vol 2017 ◽  
pp. 1-8 ◽  
Author(s):  
Wanzeng Kong ◽  
Jinshuai Yu ◽  
Ying Cheng ◽  
Weihua Cong ◽  
Huanhuan Xue

With 3D imaging of the multisonar beam and serious interference of image noise, detecting objects based only on manual operation is inefficient and also not conducive to data storage and maintenance. In this paper, a set of sonar image automatic detection technologies based on 3D imaging is developed to satisfy the actual requirements in sonar image detection. Firstly, preprocessing was conducted to alleviate the noise and then the approximate position of object was obtained by calculating the signal-to-noise ratio of each target. Secondly, the separation of water bodies and strata is realized by maximum variance between clusters (OTSU) since there exist obvious differences between these two areas. Thus image segmentation can be easily implemented on both. Finally, the feature extraction is carried out, and the multidimensional Bayesian classification model is established to do classification. Experimental results show that the sonar-image-detection technology can effectively detect the target and meet the requirements of practical applications.


1980 ◽  
Vol 7 (1) ◽  
pp. 93-104 ◽  
Author(s):  
A.W. Peterson ◽  
T. Blench

This paper, for river engineers and their environmental counterparts, presents and explains the origin and potential of four-dimensional charts that smooth most of the world's numerical data obtained from the equilibrium dimensions of sand rivers, gravel rivers, and laboratory flumes. These charts aim to provide a practical service comparable with that provided by factual plots on the comprehensive classic three-dimensional Stanton friction-factor diagram for circular pipes and clean Newtonian fluid. In the river problems, especially, the existence of different phases (whose transitions are not susceptible to formulation), the inadequacies of textbook theories even for simple phases, and the unavoidable imperfections of both field and laboratory measurements combine to prevent responsible design. The remedy is a graphing of total information backed by references from which its reliability and practicability can be assessed.The references have been chosen to contain principal information in the forms of: (i) usable photos, graphs, and tables; (ii) explanations free from specialized mathematics and speculative arguments; and (iii) papers with discussions, authors' replies, and further useful references (since a major reference list would be too long for this paper). Because condensation has had to be extreme the authors will be glad to attempt answers to discussions and questions on the subject matter, its practical applications, and its implications in teaching and research.


Sign in / Sign up

Export Citation Format

Share Document