Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userValid Physical Processes from Numerical Discontinuities in Computational Fluid Dynamics
165
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Due to the limited cell resolution in the representation of flow variables, a piecewise continuous initial reconstruction with discontinuous jump at a cell interface is usually used in modern computational fluid dynamics methods. Starting from the discontinuity, a Riemann problem in the Godunov method is solved for the flux evaluation across the cell interface in a finite volume scheme. With the increasing of Mach number in the CFD simulations, the adaptation of the Riemann solver seems introduce intrinsically a mechanism to develop instabilities in strong shock regions. Theoretically, the Riemann solution of the Euler equations are based on the equilibrium assumption, which may not be valid in the non-equilibrium shock layer. In order to clarify the flow physics from a discontinuity, the unsteady flow behavior of one-dimensional contact and shock wave is studied on a time scale of (0~10000) times of the particle collision time. In the study of the non-equilibrium flow behavior from a discontinuity, the collision-less Boltzmann equation is first used for the time scale within one particle collision time, then the direct simulation Monte Carlo (DSMC) method will be adapted to get the further evolution solution. The transition from the free particle transport to the dissipative Navier-Stokes (NS) solutions are obtained as an increasing of time. The exact Riemann solution becomes a limiting solution with infinite number of particle collisions. For the high Mach number flow simulations, the points in the shock transition region, even though the region is enlarged numerically to the mesh size, should be considered as the points inside a highly non-equilibrium shock layer.
rate research
Read More
In this article, we address both recent advances and open questions in some mathematical and computational issues in geophysical fluid dynamics (GFD) and climate dynamics. The main focus is on 1) the primitive equations (PEs) models and their related mathematical and computational issues, 2) climate variability, predictability and successive bifurcation, and 3) a new dynamical systems theory and its applications to GFD and climate dynamics.
Exact solutions of the linear water-wave problem describing oblique waves over a submerged horizontal cylinder of small (but otherwise fairly arbitrary) cross-section in a two-layer fluid are constructed in the form of convergent series in powers of the small parameter characterizing the thinness of the cylinder. The terms of these series are expressed through the solution of the exterior Neumann problem for the Laplace equation describing the flow of unbounded fluid past the cylinder. The solutions obtained describe trapped modes corresponding to discrete eigenvalues of the problem (lying close to the cut-off frequency of the continuous spectrum) and resonances lying close to the embedded cut-off. We present certain conditions for the submergence of the cylinder in the upper layer when these resonances convert into previously unobserved embedded trapped modes.
Computational Fluid Dynamics (CFD) is a major sub-field of engineering. Corresponding flow simulations are typically characterized by heavy computational resource requirements. Often, very fine and complex meshes are required to resolve physical effects in an appropriate manner. Since all CFD algorithms scale at least linearly with the size of the underlying mesh discretization, finding an optimal mesh is key for computational efficiency. One methodology used to find optimal meshes is goal-oriented adaptive mesh refinement. However, this is typically computationally demanding and only available in a limited number of tools. Within this contribution, we adopt a machine learning approach to identify optimal mesh densities. We generate optimized meshes using classical methodologies and propose to train a convolutional network predicting optimal mesh densities given arbitrary geometries. The proposed concept is validated along 2d wind tunnel simulations with more than 60,000 simulations. Using a training set of 20,000 simulations we achieve accuracies of more than 98.7%. Corresponding predictions of optimal meshes can be used as input for any mesh generation and CFD tool. Thus without complex computations, any CFD engineer can start his predictions from a high quality mesh.
In a particle physics dynamics, we assume a uniform distribution as the physical measure and a measure-theoretic definition of entropy on the velocity configuration space. This distribution is labeled as the physical solution in the remainder of the article. The dynamics is governed by an assumption of a Lagrangian formulation, with the velocity time derivatives as the momenta conjugate to the velocity configurations. From these definitions and assumptions, we show mathematically that a maximum entropy production principle selects the physical measure from among alternate solutions of the Navier-Stokes and Euler equations, but its transformation to an Eulerian frame is not established here, a topic that will be considered separately.
The nonlinear Helmholtz equation (NLH) models the propagation of electromagnetic waves in Kerr media, and describes a range of important phenomena in nonlinear optics and in other areas. In our previous work, we developed a fourth order method for its numerical solution that involved an iterative solver based on freezing the nonlinearity. The method enabled a direct simulation of nonlinear self-focusing in the nonparaxial regime, and a quantitative prediction of backscattering. However, our simulations showed that there is a threshold value for the magnitude of the nonlinearity, above which the iterations diverge. In this study, we numerically solve the one-dimensional NLH using a Newton-type nonlinear solver. Because the Kerr nonlinearity contains absolute values of the field, the NLH has to be recast as a system of two real equations in order to apply Newtons method. Our numerical simulations show that Newtons method converges rapidly and, in contradistinction with the iterations based on freezing the nonlinearity, enables computations for very high levels of nonlinearity. In addition, we introduce a novel compact finite-volume fourth order discretization for the NLH with material discontinuities.The one-dimensional results of the current paper create a foundation for the analysis of multi-dimensional problems in the future.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
