Numerical solutions of the Euler equations for complex three-dimensional aerodynamic configurations

1984 ◽  
Author(s):  
J. DEESE ◽  
R. AGARWAL ◽  
R. UNDERWOOD
1989 ◽  
Author(s):  
N. KROLL ◽  
C. ROSSOW ◽  
S. SCHERR ◽  
J. SCHOENE ◽  
G. WICHMANN

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.


2013 ◽  
Vol 80 (2) ◽  
Author(s):  
D. Zhou ◽  
S. H. Lo

The three-dimensional (3D) free vibration of twisted cylinders with sectorial cross section or a radial crack through the height of the cylinder is studied by means of the Chebyshev–Ritz method. The analysis is based on the three-dimensional small strain linear elasticity theory. A simple coordinate transformation is applied to map the twisted cylindrical domain into a normal cylindrical domain. The product of a triplicate Chebyshev polynomial series along with properly defined boundary functions is selected as the admissible functions. An eigenvalue matrix equation can be conveniently derived through a minimization process by the Rayleigh–Ritz method. The boundary functions are devised in such a way that the geometric boundary conditions of the cylinder are automatically satisfied. The excellent property of Chebyshev polynomial series ensures robustness and rapid convergence of the numerical computations. The present study provides a full vibration spectrum for thick twisted cylinders with sectorial cross section, which could not be determined by 1D or 2D models. Highly accurate results presented for the first time are systematically produced, which can serve as a benchmark to calibrate other numerical solutions for twisted cylinders with sectorial cross section. The effects of height-to-radius ratio and twist angle on frequency parameters of cylinders with different subtended angles in the sectorial cross section are discussed in detail.


2021 ◽  
Author(s):  
Ramtin Sabeti ◽  
Mohammad Heidarzadeh

<p>Landslide-generated waves have been major threats to coastal areas and have led to destruction and casualties. Their importance is undisputed, most recently demonstrated by the 2018 Anak Krakatau tsunami, causing several hundred fatalities. The accurate prediction of the maximum initial amplitude of landslide waves (<em>η<sub>max</sub></em>) around the source region is a vital hazard indicator for coastal impact assessment. Laboratory experiments, analytical solutions and numerical modelling are three major methods to investigate the (<em>η<sub>max</sub></em>). However, the numerical modelling approach provides a more flexible and cost- and time-efficient tool. This research presents a numerical simulation of tsunamis due to rigid landslides with consideration of submerged conditions. In particular, this simulation focuses on studying the effect of landslide parameters on <em>η<sub>max</sub>.</em> Results of simulations are compared with our conducted physical experiments at the Brunel University London (UK) to validate the numerical model.</p><p>We employ the fully three-dimensional computational fluid dynamics package, FLOW-3D Hydro for modelling the landslide-generated waves. This software benefit from the Volume of Fluid Method (VOF) as the numerical technique for tracking and locating the free surface. The geometry of the simulation is set up according to the wave tank of physical experiments (i.e. 0.26 m wide, 0.50 m deep and 4.0 m). In order to calibrate the simulation model based on the laboratory measurements, the friction coefficient between solid block and incline is changed to 0.41; likewise, the terminal velocity of the landslide is set to 0.87 m/s. Good agreement between the numerical solutions and the experimental results is found. Sensitivity analyses of landslide parameters (e.g. slide volume, water depth, etc.) on <em>η<sub>max </sub></em>are performed. Dimensionless parameters are employed to study the sensitivity of the initial landslide waves to various landslide parameters.</p>


2017 ◽  
Vol 140 (1) ◽  
Author(s):  
Sofia Sarraf ◽  
Ezequiel López ◽  
Laura Battaglia ◽  
Gustavo Ríos Rodríguez ◽  
Jorge D'Elía

In the boundary element method (BEM), the Galerkin weighting technique allows to obtain numerical solutions of a boundary integral equation (BIE), giving the Galerkin boundary element method (GBEM). In three-dimensional (3D) spatial domains, the nested double surface integration of GBEM leads to a significantly larger computational time for assembling the linear system than with the standard collocation method. In practice, the computational time is roughly an order of magnitude larger, thus limiting the use of GBEM in 3D engineering problems. The standard approach for reducing the computational time of the linear system assembling is to skip integrations whenever possible. In this work, a modified assembling algorithm for the element matrices in GBEM is proposed for solving integral kernels that depend on the exterior unit normal. This algorithm is based on kernels symmetries at the element level and not on the flow nor in the mesh. It is applied to a BIE that models external creeping flows around 3D closed bodies using second-order kernels, and it is implemented using OpenMP. For these BIEs, the modified algorithm is on average 32% faster than the original one.


Energies ◽  
2020 ◽  
Vol 13 (23) ◽  
pp. 6424
Author(s):  
Cheng-Hung Huang ◽  
Chih-Yang Kuo

A non-linear three-dimensional inverse shape design problem was investigated for a pipe type heat exchanger to estimate the design variables of continuous lateral ribs on internal Z-shape lateral fins for maximum thermal performance factor η. The design variables were considered as the positions, heights, and number of ribs while the physical properties of air were considered as a polynomial function of temperature; this makes the problem non-linear. The direct problem was solved using software package CFD-ACE+, and the Levenberg–Marquardt method (LMM) was utilized as the optimization tool because it has been proven to be a powerful algorithm for solving inverse problems. Z-shape lateral fins were found to be the best thermal performance among Z-shape, S-shape, and V-shape lateral fins. The objective of this study was to include continuous lateral ribs to Z-shape lateral fins to further improve η. Firstly, the numerical solutions of direct problem were solved using both polynomial and constant air properties and then compared with the corrected solutions to verify the necessity for using polynomial air properties. Then, four design cases, A, B, C and D, based on various design variables were conducted numerically, and the resultant η values were computed and compared. The results revealed that considering continuous lateral ribs on the surface of Z-shape lateral fins can indeed improve η value at the design working condition Re = 5000. η values of designs A, B and C were approximately 13% higher than that for Z-shape lateral fins, however, when the rib numbers were increased, i.e., design D, the value of η became only 11.5 % higher. This implies that more ribs will not guarantee higher η value.


Nonlinearity ◽  
2008 ◽  
Vol 21 (8) ◽  
pp. T123-T129 ◽  
Author(s):  
J D Gibbon ◽  
M Bustamante ◽  
R M Kerr

Sign in / Sign up

Export Citation Format

Share Document