scholarly journals Computational accuracy of the upwelling radiation in simulating the inhomogeneous atmosphere by homogeneous sub-layers.

1986 ◽  
Vol 37 (2) ◽  
pp. 83-90
Author(s):  
Tsutomu Takashima
Author(s):  
Takaaki OKUMURA ◽  
Atsushi KUROKAWA ◽  
Hiroo MASUDA ◽  
Toshiki KANAMOTO ◽  
Masanori HASHIMOTO ◽  
...  

2019 ◽  
Vol 11 (01) ◽  
pp. 1950006 ◽  
Author(s):  
Fengbin Liu ◽  
Qiang Wu ◽  
Yumin Cheng

In this study, based on a nonsingular weight function, the improved element-free Galerkin (IEFG) method is presented for solving elastoplastic large deformation problems. By using the improved interpolating moving least-squares (IMLS) method to form the approximation function, and using Galerkin weak form based on total Lagrange formulation of elastoplastic large deformation problems to form the discretilized equations, which is solved with the Newton–Raphson iteration method, we obtain the formulae of the IEFG method for elastoplastic large deformation problems. In numerical examples, the influences of the penalty factor, scale parameter of influence domain and weight functions on the computational accuracy are analyzed, and the numerical solutions show that the IEFG method for elastoplastic large deformation problems has higher computational efficiency and accuracy.


1978 ◽  
Vol 45 (4) ◽  
pp. 711-716 ◽  
Author(s):  
Stephen S.-H. Chang

This paper presents a method for computing the motion and decay of a large dusty, buoyant thermal (cloud) carried by a vortex ring generated from a strong near ground explosion and ascending in an inhomogeneous atmosphere. A system of equations is derived describing the motion of the vortex ring, the thermal, and the pollutants which consist of numerous solid spherical particles. The interior properties and the trajectories of the thermal and the pollutants are obtained. The numerical solution for the thermal trajectory is in excellent agreement with experiment.


Algorithms ◽  
2021 ◽  
Vol 14 (5) ◽  
pp. 129
Author(s):  
Yuan Li ◽  
Ni Zhang ◽  
Yuejiao Gong ◽  
Wentao Mao ◽  
Shiguang Zhang

Compared with continuous elements, discontinuous elements advance in processing the discontinuity of physical variables at corner points and discretized models with complex boundaries. However, the computational accuracy of discontinuous elements is sensitive to the positions of element nodes. To reduce the side effect of the node position on the results, this paper proposes employing partially discontinuous elements to compute the time-domain boundary integral equation of 3D elastodynamics. Using the partially discontinuous element, the nodes located at the corner points will be shrunk into the element, whereas the nodes at the non-corner points remain unchanged. As such, a discrete model that is continuous on surfaces and discontinuous between adjacent surfaces can be generated. First, we present a numerical integration scheme of the partially discontinuous element. For the singular integral, an improved element subdivision method is proposed to reduce the side effect of the time step on the integral accuracy. Then, the effectiveness of the proposed method is verified by two numerical examples. Meanwhile, we study the influence of the positions of the nodes on the stability and accuracy of the computation results by cases. Finally, the recommended value range of the inward shrink ratio of the element nodes is provided.


Author(s):  
Zhihua Niu ◽  
Sun Jin ◽  
Rongrong Wang ◽  
Yansong Zhang

Dynamic analysis is an essential task in the geometry design of suspension systems. Whereas the dynamic simulation based on numerical software like Adams is quite slowly and the existing analytical models of the nonlinear suspension geometry are mostly based on small displacement hypothesis, this paper aims to propose a whole-range dynamic model with high computational efficiency for planar double wishbone suspensions and further achieve the fast optimal design of suspension geometry. Selection of the new generalized coordinate and explicit solutions of the basic four-bar mechanism dramatically reduce the complexity of suspension geometry representation and provide analytical solutions for all of the time varying dimensions. By this means, the running speed and computational accuracy of the new model are guaranteed simultaneously. Furthermore, an original Matlab/Simulink implementation is given to maintain the geometric nonlinearity in the solving process of dynamic differential equations. After verifying its accuracy with an ADAMS prototype, the presented whole-range model is used in the vast-parameter optimization of suspension geometry. Since both kinematic and dynamic performances are evaluated in the objective function, the optimization is qualified to give a comprehensive suggestion to the design of suspension geometry.


2020 ◽  
Vol 2 (1) ◽  
Author(s):  
Lingfa Kong ◽  
Yidao Dong ◽  
Wei Liu ◽  
Huaibao Zhang

AbstractAccuracy of unstructured finite volume discretization is greatly influenced by the gradient reconstruction. For the commonly used k-exact reconstruction method, the cell centroid is always chosen as the reference point to formulate the reconstructed function. But in some practical problems, such as the boundary layer, cells in this area are always set with high aspect ratio to improve the local field resolution, and if geometric centroid is still utilized for the spatial discretization, the severe grid skewness cannot be avoided, which is adverse to the numerical performance of unstructured finite volume solver. In previous work [Kong, et al. Chin Phys B 29(10):100203, 2020], we explored a novel global-direction stencil and combined it with the face-area-weighted centroid on unstructured finite volume methods from differential form to realize the skewness reduction and a better reflection of flow anisotropy. Greatly inspired by the differential form, in this research, we demonstrate that it is also feasible to extend this novel method to the unstructured finite volume discretization from integral form on both second and third-order finite volume solver. Numerical examples governed by linear convective, Euler and Laplacian equations are utilized to examine the correctness as well as effectiveness of this extension. Compared with traditional vertex-neighbor and face-neighbor stencils based on the geometric centroid, the grid skewness is almost eliminated and computational accuracy as well as convergence rate is greatly improved by the global-direction stencil with face-area-weighted centroid. As a result, on unstructured finite volume discretization from integral form, the method also has superiorities on both computational accuracy and convergence rate.


2019 ◽  
Vol 2019 ◽  
pp. 1-9
Author(s):  
Haitao Wang ◽  
Xiangyang Zeng ◽  
Ye Lei

Wave-based methods for acoustic simulations within enclosures suffer the numerical dispersion and then usually have evident dispersion error for problems with high wave numbers. To improve the upper limit of calculating frequency for 3D problems, a hybrid smoothed finite element method (hybrid SFEM) is proposed in this paper. This method employs the smoothing technique to realize the reduction of the numerical dispersion. By constructing a type of mixed smoothing domain, the traditional node-based and face-based smoothing techniques are mixed in the hybrid SFEM to give a more accurate stiffness matrix, which is widely believed to be the ultimate cause for the numerical dispersion error. The numerical examples demonstrate that the hybrid SFEM has better accuracy than the standard FEM and traditional smoothed FEMs under the condition of the same basic elements. Moreover, the hybrid SFEM also has good performance on the computational efficiency. A convergence experiment shows that it costs less time than other comparison methods to achieve the same computational accuracy.


Sign in / Sign up

Export Citation Format

Share Document