No Arabic abstract
One of the most important aspects in tsunami studies is the wave behavior when it approaches the coast. Information on physical parameters that characterize waves is often limited because of the diffilculties in achieving accurate measurements at the time of the event. The impact of a tsunami on the coast is governed by nonlinear physics such as turbulence with spatial and temporal variability. The use of the Smoothed Particle Hydrodynamic method (SPH) presents advantages over models based on two-dimensional Shallow Waters Equations (SWE), because the assumed vertical velocity simplifies hydrodynamics in two dimensions. The study presented here reports numerical SPH simulations of the tsunami event occurred in Coquimbo (Chile) on September 16 of 2015. On the basis of the reconstruction of the physical parameters that characterized this event (flow velocities, direction and water elevations), calibrated by a reference rodel, force values on buildings located on the study coast were numerically calculated, and compared with an estimate of the Chilean Structural Design Standard. Finally, discussion and conclusions of the comparison of both methodologies are presented, including an influence analysis of the topographical detail of the model in the estimation of hydrodynamic forces.
We present a three--pronged numerical approach to the dynamo problem at low magnetic Prandtl numbers $P_M$. The difficulty of resolving a large range of scales is circumvented by combining Direct Numerical Simulations, a Lagrangian-averaged model, and Large-Eddy Simulations (LES). The flow is generated by the Taylor-Green forcing; it combines a well defined structure at large scales and turbulent fluctuations at small scales. Our main findings are: (i) dynamos are observed from $P_M=1$ down to $P_M=10^{-2}$; (ii) the critical magnetic Reynolds number increases sharply with $P_M^{-1}$ as turbulence sets in and then saturates; (iii) in the linear growth phase, the most unstable magnetic modes move to small scales as $P_M$ is decreased and a Kazantsev $k^{3/2}$ spectrum develops; then the dynamo grows at large scales and modifies the turbulent velocity fluctuations.
This paper presents an extension of the hybrid scheme proposed by Wang et al. (J. Comput. Phys. 229 (2010) 169-180) for numerical simulation of compressible isotropic turbulence to flows with higher turbulent Mach numbers. The scheme still utilizes an 8th-order compact scheme with built-in hyperviscosity for smooth regions and a 7th-order WENO scheme for highly compressive regions, but now both in their conservation formulations and for the latter with the Roe type characteristic-wise reconstruction. To enhance the robustness of the WENO scheme without compromising its high-resolution and accuracy, the recursive-order-reduction procedure is adopted, where a new type of reconstruction-failure-detection criterion is constructed. To capture the upwind direction properly in extreme conditions, the global Lax-Friedrichs numerical flux is used. In addition, a new form of cooling function is proposed, which is proved to be positivity-preserving. With these techniques, the new scheme not only inherits the good properties of the original one but also extends largely the computable range of turbulent Mach number, which has been further confirmed by numerical results.
Corridors of size-selected crescent-shaped dunes, known as barchans, are commonly found in water, air, and other planetary environments. The growth of barchans results from the interplay between a fluid flow and a granular bed, but their size regulation involves intricate exchanges between different barchans within a field. One size-regulating mechanism is the binary interaction between nearby dunes, when two dunes exchange mass via the near flow field or by direct contact (collision). In a recent Letter (Assis and Franklin, Geophys. Res. Lett., 2020), we identified five different patterns arising from binary interactions of subaqueous barchans, and proposed classification maps. In this paper, we further inquire into binary exchanges by investigating the motion of individual grains while barchans interact with each other. The experiments were conducted in a water channel where the evolution of pairs of barchans in both aligned and off-centered configurations was recorded by conventional and high-speed cameras. Based on image processing, we obtained the morphology of dunes and motion of grains for all interaction patterns. We present the trajectories of individual grains, from which we show the origin and destination of moving grains, and their typical lengths and velocities. We also show that grains from the impacting dune spread with a diffusion-like component over the target barchan, and we propose a diffusion length. Our results provide new insights into the size-regulating mechanisms of barchans and barchanoid forms found on Earth and other planets.
We propose a super-resolution (SR) simulation system that consists of a physics-based meteorological simulation and an SR method based on a deep convolutional neural network (CNN). The CNN is trained using pairs of high-resolution (HR) and low-resolution (LR) images created from meteorological simulation results for different resolutions so that it can map LR simulation images to HR ones. The proposed SR simulation system, which performs LR simulations, can provide HR prediction results in much shorter operating cycles than those required for corresponding HR simulation prediction system. We apply the SR simulation system to urban micrometeorology, which is strongly affected by buildings and human activity. Urban micrometeorology simulations that need to resolve urban buildings are computationally costly and thus cannot be used for operational real-time predictions even when run on supercomputers. We performed HR micrometeorology simulations on a supercomputer to obtain datasets for training the CNN in the SR method. It is shown that the proposed SR method can be used with a spatial scaling factor of 4 and that it outperforms conventional interpolation methods by a large margin. It is also shown that the proposed SR simulation system has the potential to be used for operational urban micrometeorology predictions.
The intent of this paper is to discuss the history and origins of Lagrangian hydrodynamic methods for simulating shock driven flows. The majority of the pioneering research occurred within the Manhattan Project. A range of Lagrangian hydrodynamic schemes were created between 1943 and 1948 by John von Neumann, Rudolf Peierls, Tony Skyrme, and Robert Richtmyer. These schemes varied significantly from each other; however, they all used a staggered-grid and finite difference approximations of the derivatives in the governing equations, where the first scheme was by von Neumann. These ground-breaking schemes were principally published in Los Alamos laboratory reports that were eventually declassified many decades after authorship, which motivates us to document the work and describe the accompanying history in a paper that is accessible to the broader scientific community. Furthermore, we seek to correct historical omissions on the pivotal contributions made by Peierls and Skyrme to creating robust Lagrangian hydrodynamic methods for simulating shock driven flows. Understanding the history of Lagrangian hydrodynamic methods can help explain the origins of many modern schemes and may inspire the pursuit of new schemes.